Detecting signatures of selection through haplotype differentiation among hierarchically structured populations

Exome-wide comparative analyses revealed differentiating genomic regions for performance traits in Indian native buffaloes

In India, 20 breeds of buffalo have been identified and registered, yet limited studies have been conducted to explore the performance potential of these breeds, especially in the Indian native breeds. This study is a maiden attempt to delineate the important variants and unique genes through exome sequencing for milk yield, milk composition, fertility, and adaptation traits in Indian local breeds of buffalo. In the present study, whole exome sequencing was performed on Chhattisgarhi (n = 3), Chilika (n = 4), Gojri (n = 3), and Murrah (n = 4) buffalo breeds and after stringent quality control, 4333, 6829, 4130, and 4854 InDels were revealed, respectively. Exome-wide FST along 100-kb sliding windows detected 27, 98, 38, and 35 outlier windows in Chhattisgarhi, Chilika, Gojri, and Murrah, respectively. The comparative exome analysis of InDels and subsequent gene ontology revealed unique breed specific genes for milk yield (CAMSAP3), milk composition (CLCN1, NUDT3), fertility (PTGER3) and adaptation (KCNA3, TH) traits. Study provides insight into mechanism of how these breeds have evolved under natural selection, the impact of these events on their respective genomes, and their importance in maintaining purity of these breeds for the traits under study. Additionally, this result will underwrite to the genetic acquaintance of these breeds for breeding application, and in understanding of evolution of these Indian local breeds.

Selection signatures for high altitude adaptation in livestock: A review

High altitude adapted livestock species (cattle, yak, goat, sheep, and horse) has critical role in the human socioeconomic sphere and acts as good source of animal source products including milk, meat, and leather, among other things. These species sustain production and reproduction even in harsh environments on account of adaptation resulting from continued evolution of beneficial traits. Selection pressure leads to various adaptive strategies in livestock whose footprints are evident at the different genomic sites as the "Selection Signature". Scrutiny of these signatures provides us crucial insight into the evolutionary process and domestication of livestock adapted to diverse climatic conditions. These signatures have the potential to change the sphere of animal breeding and further usher the selection programmes in right direction. Technological revolution and recent strides made in genomic studies has opened the routes for the identification of selection signatures. Numerous statistical approaches and bioinformatics tools have been developed to detect the selection signature. Consequently, studies across years have identified candidate genes under selection region found associated with numerous traits which have a say in adaptation to high-altitude environment. This makes it pertinent to have a better understanding about the selection signature, the ways to identify and how to utilize them for betterment of livestock populations as well as farmers. This review takes a closer look into the general concept, various methodologies, and bioinformatics tools commonly employed in selection signature studies and summarize the results of recent selection signature studies related to high-altitude adaptation in various livestock species. This review will serve as an informative and useful insight for researchers and students in the field of animal breeding and evolutionary biology.

Dissecting Selective Signatures and Candidate Genes in Grandparent Lines Subject to High Selection Pressure for Broiler Production and in a Local Russian Chicken Breed of Ushanka

Breeding improvements and quantitative trait genetics are essential to the advancement of broiler production. The impact of artificial selection on genomic architecture and the genetic markers sought remains a key area of research. Here, we used whole-genome resequencing data to analyze the genomic architecture, diversity, and selective sweeps in Cornish White (CRW) and Plymouth Rock White (PRW) transboundary breeds selected for meat production and, comparatively, in an aboriginal Russian breed of Ushanka (USH). Reads were aligned to the reference genome bGalGal1.mat.broiler.GRCg7b and filtered to remove PCR duplicates and low-quality reads using BWA-MEM2 and bcftools software; 12,563,892 SNPs were produced for subsequent analyses. Compared to CRW and PRW, USH had a lower diversity and a higher genetic distinctiveness. Selective sweep regions and corresponding candidate genes were examined based on ZFST, hapFLK, and ROH assessment procedures. Twenty-seven prioritized chicken genes and the functional projection from human homologs suggest their importance for selection signals in the studied breeds. These genes have a functional relationship with such trait categories as body weight, muscles, fat metabolism and deposition, reproduction, etc., mainly aligned with the QTLs in the sweep regions. This information is pivotal for further executing genomic selection to enhance phenotypic traits.

Common signatures of selection reveal target loci for breeding across soybean populations

Understanding the underlying genetic bases of yield-related selection and distinguishing these changes from genetic drift are critical for both improved understanding and future success of plant breeding. Soybean [Glycine max (L.) Merr.] is a key species for world food security, yet knowledge of the mechanism of selective breeding in soybean, such as the century-long program of artificial selection in U.S. soybean germplasm, is currently limited to certain genes and loci. Here, we identify genome-wide signatures of selection in separate populations of soybean subjected to artificial selection for increased yield by multiple breeding programs in the United States. We compared the alternative soybean breeding population (AGP) created by USDA-ARS to the conventional public soybean lines (CGP) developed at three different stages of breeding (ancestral, intermediate, and elite) to identify shared signatures of selection and differentiate these from drift. The results showed a strong selection for specific haplotypes identified by single site frequency and haplotype homozygosity methods. A set of common selection signatures was identified in both AGP and CGP that supports the hypothesis that separate breeding programs within similar environments coalesce on the fixation of the same key haplotypes. Signatures unique to each breeding program were observed. These results raise the possibility that selection analysis can allow the identification of favorable alleles to enhance directed breeding approaches.

Seascapes Shaped the Local Adaptation and Population Structure of South China Coast Yellowfin Seabream (Acanthopagrus latus)

Understanding the genetic composition and regional adaptation of marine species under environmental heterogeneity and fishing pressure is crucial for responsible management. In order to understand the genetic diversity and adaptability of yellowfin seabream (Acanthopagrus latus) along southern China coast, this study was conducted a seascape genome analysis on yellowfin seabream from the ecologically diverse coast, spanning over 1600 km. A total of 92 yellowfin seabream individuals from 15 sites were performed whole-genome resequencing, and 4,383,564 high-quality single nucleotide polymorphisms (SNPs) were called. By conducting a genotype-environment association analysis, 29,951 adaptive and 4,328,299 neutral SNPs were identified. The yellowfin seabream exhibited two distinct population structures, despite high gene flow between sites. The seascape genome analysis revealed that genetic structure was influenced by a variety of factors including salinity gradients, habitat distance, and ocean currents. The frequency of allelic variation at the candidate loci changed with the salinity gradient. Annotation of these loci revealed that most of the genes are associated with osmoregulation, such as kcnab2a, kcnk5a, and slc47a1. These genes are significantly enriched in pathways associated with ion transport including G protein-coupled receptor activity, transmembrane signaling receptor activity, and transporter activity. Overall, our findings provide insights into how seascape heterogeneity affects adaptive evolution, while providing important information for regional management in yellowfin seabream populations.

Genomics and conservation: Guidance from training to analyses and applications

Environmental change is intensifying the biodiversity crisis and threatening species across the tree of life. Conservation genomics can help inform conservation actions and slow biodiversity loss. However, more training, appropriate use of novel genomic methods and communication with managers are needed. Here, we review practical guidance to improve applied conservation genomics. We share insights aimed at ensuring effectiveness of conservation actions around three themes: (1) improving pedagogy and training in conservation genomics including for online global audiences, (2) conducting rigorous population genomic analyses properly considering theory, marker types and data interpretation and (3) facilitating communication and collaboration between managers and researchers. We aim to update students and professionals and expand their conservation toolkit with genomic principles and recent approaches for conserving and managing biodiversity. The biodiversity crisis is a global problem and, as such, requires international involvement, training, collaboration and frequent reviews of the literature and workshops as we do here.

Deciphering climate resilience in Indian cattle breeds by selection signature analyses

The signature of selection is a crucial concept in evolutionary biology that refers to the pattern of genetic variation which arises in a population due to natural selection. In the context of climate adaptation, the signature of selection can reveal the genetic basis of adaptive traits that enable organisms to survive and thrive in changing environmental conditions. Breeds living in diverse agroecological zones exhibit genetic "footprints" within their genomes that mirror the influence of climate-induced selective pressures, subsequently impacting phenotypic variance. It is assumed that the genomes of animals residing in these regions have been altered through selection for various climatic adaptations. These regions are known as signatures of selection and can be identified using various summary statistics. We examined genotypic data from eight different cattle breeds (Gir, Hariana, Kankrej, Nelore, Ongole, Red Sindhi, Sahiwal, and Tharparkar) that are adapted to diverse regional climates. To identify selection signature regions in this investigation, we used four intra-population statistics: Tajima's D, CLR, iHS, and ROH. In this study, we utilized Bovine 50 K chip data and four genome scan techniques to assess the genetic regions of positive selection for high-temperature adaptation. We have also performed a genome-wide investigation of genetic diversity, inbreeding, and effective population size in our target dataset. We identified potential regions for selection that are likely to be caused by adverse climatic conditions. We observed many adaptation genes in several potential selection signature areas. These include genes like HSPB2, HSPB3, HSP20, HSP90AB1, HSF4, HSPA1B, CLPB, GAP43, MITF, and MCHR1 which have been reported in the cattle populations that live in varied climatic regions. The findings demonstrated that genes involved in disease resistance and thermotolerance were subjected to intense selection. The findings have implications for marker-assisted breeding, understanding the genetic landscape of climate-induced adaptation, putting breeding and conservation programs into action.

Identification of Signatures of Positive Selection That Have Shaped the Genomic Landscape of South African Pig Populations

South Africa boasts a diverse range of pig populations, encompassing intensively raised commercial breeds, as well as indigenous and village pigs reared under low-input production systems. The aim of this study was to investigate how natural and artificial selection have shaped the genomic landscape of South African pig populations sampled from different genetic backgrounds and production systems. For this purpose, the integrated haplotype score (iHS), as well as cross population extended haplotype homozygosity (XP-EHH) and Lewontin and Krakauer's extension of the Fst statistic based on haplotype information (HapFLK) were utilised. Our results revealed several population-specific signatures of selection associated with the different production systems. The importance of natural selection in village populations was highlighted, as the majority of genomic regions under selection were identified in these populations. Regions under natural and artificial selection causing the distinct genetic footprints of these populations also allow for the identification of genes and pathways that may influence production and adaptation. In the context of intensively raised commercial pig breeds (Large White, Kolbroek, and Windsnyer), the identified regions included quantitative loci (QTLs) associated with economically important traits. For example, meat and carcass QTLs were prevalent in all the populations, showing the potential of village and indigenous populations' ability to be managed and improved for such traits. Results of this study therefore increase our understanding of the intricate interplay between selection pressures, genomic adaptations, and desirable traits within South African pig populations.

Whole-genome selective sweep analyses identifies the region and candidate gene associated with white earlobe color in Mediterranean chickens

We compared the genomes of multiple domestic chicken breeds with red and white earlobes to identify the differentiated regions between groups of breeds differing in earlobe color. This was done using a selective sweep mapping approach based on whole-genome sequence data. The most significant selective sweep was identified on chromosome 11, where the white earlobe chicken breeds originated from Mediterranean share a common haplotype, and where multiple candidate genes are located. The most plausible functional candidate gene is the Melanocortin 1 Receptor (MC1R), a receptor known to regulate pigmentation in the skin and hair, and it is also the gene with the strongest positional support from the haplotype-based analyses. It, however, still needs to be explored experimentally to identify effects also on chicken earlobe color variation. Our study is the first exploration of the genetic basis of white earlobe color in Mediterranean chickens using a selective sweep mapping method based on whole-genome sequencing data and shows its value for identifying likely functional genes mediating the pigmentation in earlobe. It also indicates a potential novel role of MC1R in birds and exemplifies how selection on fancy traits has influenced the genome during formation of the modern chicken breeds.

A meta-analysis of genetic and phenotypic diversity of European local pig breeds reveals genomic regions associated with breed differentiation for production traits

BACKGROUND

Intense selection of modern pig breeds has resulted in genetic improvement of production traits while the performance of local pig breeds has remained lower. As local pig breeds have been bred in extensive systems, they have adapted to specific environmental conditions, resulting in a rich genotypic and phenotypic diversity. This study is based on European local pig breeds that have been genetically characterized using DNA-pool sequencing data and phenotypically characterized using breed level phenotypes related to stature, fatness, growth, and reproductive performance traits. These data were analyzed using a dedicated approach to detect signatures of selection linked to phenotypic traits in order to uncover potential candidate genes that may underlie adaptation to specific environments.

RESULTS

Analysis of the genetic data of European pig breeds revealed four main axes of genetic variation represented by the Iberian and three modern breeds (i.e. Large White, Landrace, and Duroc). In addition, breeds clustered according to their geographical origin, for example French Gascon and Basque breeds, Italian Apulo Calabrese and Casertana breeds, Spanish Iberian, and Portuguese Alentejano breeds. Principal component analysis of the phenotypic data distinguished the larger and leaner breeds with better growth potential and reproductive performance from the smaller and fatter breeds with low growth and reproductive efficiency. Linking the signatures of selection with phenotype identified 16 significant genomic regions associated with stature, 24 with fatness, 2 with growth, and 192 with reproduction. Among them, several regions contained candidate genes with possible biological effects on stature, fatness, growth, and reproductive performance traits. For example, strong associations were found for stature in two regions containing, respectively, the ANXA4 and ANTXR1 genes, for fatness in a region containing the DNMT3A and POMC genes and for reproductive performance in a region containing the HSD17B7 gene.

CONCLUSIONS

In this study on European local pig breeds, we used a dedicated approach for detecting signatures of selection that were supported by phenotypic data at the breed level to identify potential candidate genes that may have adapted to different living environments and production systems.

Mapping and Functional Dissection of the Rumpless Trait in Piao Chicken Identifies a Causal Loss of Function Mutation in the Novel Gene Rum

Rumpless chickens exhibit an abnormality in their tail development. The genetics and biology of this trait has been studied for decades to illustrate a broad variation in both the types of inheritance and the severity in the developmental defects of the tail. In this study, we created a backcross pedigree by intercrossing Piao (rumpless) with Xianju (normal) to investigate the genetic mechanisms and molecular basis of the rumpless trait in Piao chicken. Through genome-wide association and linkage analyses, the candidate region was fine-mapped to 798.5 kb (chromosome 2: 86.9 to 87.7 Mb). Whole-genome sequencing analyses identified a single variant, a 4.2 kb deletion, which was completely associated with the rumpless phenotype. Explorations of the expression data identified a novel causative gene, Rum, that produced a long, intronless transcript across the deletion. The expression of Rum is embryo-specific, and it regulates the expression of MSGN1, a key factor in regulating T-box transcription factors required for mesoderm formation and differentiation. These results provide genetic and molecular experimental evidence for a novel mechanism regulating tail development in chicken and report the likely causal mutation for the tail abnormity in the Piao chicken. The novel regulatory gene, Rum, will, due to its role in fundamental embryo development, be of interest for further explorations of a potential role in tail and skeletal development also in other vertebrates.

Signatures of selection in riverine buffalo populations revealed by genome-wide SNP data

The detection of selection signatures assists in understanding domestication, evolution, and the identification of genomic regions related to adaptation and production traits in buffaloes. The emergence of high-throughput technologies like Next Generation Sequencing and SNP genotyping had expanded our ability to detect these signatures of selection. In this study, we sought to identify signatures of selection in five buffalo populations (Brazilian Murrah, Bulgarian Murrah, Indian Murrah, Nili-Ravi, and Kundi) using Axiom Buffalo 90 K Genotyping Array data. Using seven different methodologies (Tajima's D, CLR, ROH, iHS, FST, FLK and hapFLK), we identified selection signatures in 374 genomic regions, spanning a total of 381 genes and 350 quantitative trait loci (QTLs). Among these, several candidate genes were associated with QTLs for milk production, reproduction, growth and carcass traits. The genes and QTLs reported in this study provide insight into selection signals shaping the genome of buffalo breeds. Our findings can aid in further genomic association studies, genomic prediction, and the implementation of breeding programmes in Indian buffaloes.

A genomic analysis of the University of Nebraska Replicated Recurrent Selection program

The inbred-hybrid system of maize breeding closely resembles a reciprocal full-sib (RFS) selection program. Studying changes in genetic variation as a result of RFS selection can help illuminate long-standing questions regarding the relative roles of selection and genetic drift and help understand the nature of adaptation occurring in selection programs. The University of Nebraska-Lincoln Replicated Recurrent Selection (UNL-RpRS) program underwent eight cycles of replicated RFS and S1-progeny selection, making it a powerful system to study genomic changes accompanying selection for inter-population performance. The objectives of this study were to identify regions of the genome under selection after eight cycles of selection and evaluate the effect eight cycles of selection for inter-population full-sib performance had in expanding genome-wide and localized population structure. We address these questions with a large set of individuals sampled from the UNL-RpRS program with dense genotyping-by-sequence data. We found evidence of parallel selection signatures in the UNL-RpRS program, with a region on chromosome 7 being implicated in three of the four selection systems studied. Regions that displayed selection signatures across independently run selection programs represent regions likely to be capitalizing on standing genetic variation and support a soft sweep model of adaptation. We did not find selection to be a strong force in diverging populations undergoing RFS. This could be due to the nature of adaptation occurring in these populations, underlying gene action, or a result of unstable genetic topographies.

Candidate genes for domestication and resistance to cold climate according to whole genome sequencing data of Russian cattle and sheep breeds

It is known that different species of animals, when living in the same environmental conditions, can form similar phenotypes. The study of the convergent evolution of several species under the influence of the same environmental factor makes it possible to identify common mechanisms of genetic adaptation. Local cattle and sheep breeds have been formed over thousands of years under the influence of domestication, as well as selection aimed at adaptation to the local environment and meeting human needs. Previously, we identified a number of candidate genes in genome regions potentially selected during domestication and adaptation to the climatic conditions of Russia, in local breeds of cattle and sheep using whole genome genotyping data. However, these data are of low resolution and do not reveal most nucleotide substitutions. The aim of the work was to create, using the whole genome sequencing data, a list of genes associated with domestication, selection and adaptation in Russian cattle and sheep breeds, as well as to identify candidate genes and metabolic pathways for selection for cold adaptation. We used our original data on the search for signatures of selection in the genomes of Russian cattle (Yakut, Kholmogory, Buryat, Wagyu) and sheep (Baikal, Tuva) breeds. We used the HapFLK, DCMS, FST and PBS methods to identify DNA regions with signatures of selection. The number of candidate genes in potentially selective regions was 946 in cattle and 151 in sheep. We showed that the studied Russian cattle and sheep breeds have at least 10 genes in common, apparently involved in the processes of adaptation/selection, including adaptation to a cold climate, including the ASTN2, PM20D1, TMEM176A, and GLIS1 genes. Based on the intersection with the list of selected genes in at least two Arctic/Antarctic mammal species, 20 and 8 genes, have been identified in cattle and sheep, respectively, that are potentially involved in cold adaptation. Among them, the most promising for further research are the ASPH, NCKAP5L, SERPINF1, and SND1 genes. Gene ontology analysis indicated the existence of possible common biochemical pathways for adaptation to cold in domestic and wild mammals associated with cytoskeleton disassembly and apoptosis.

Tropical forest loss and geographic location drive the functional genomic diversity of an endangered palm tree

Human activity has diminished forests in different terrestrial ecosystems. This is well illustrated in the Brazilian Atlantic Forest, which still hosts high levels of species richness and endemism, even with only 28% of its original extent remaining. The consequences of such forest loss in remaining populations can be investigated with several approaches, including the genomic perspective, which allows a broader understanding of how human disturbance influences the genetic variability in natural populations. In this context, our study investigated the genomic responses of Euterpe edulis Martius, an endangered palm tree, in forest remnants located in landscapes presenting different forest cover amount and composed by distinct bird assemblage that disperse its seeds. We sampled 22 areas of the Brazilian Atlantic Forest in four regions using SNP markers inserted into transcribed regions of the genome of E. edulis, distinguishing neutral loci from those putatively under natural selection (outlier). We demonstrate that populations show patterns of structure and genetic variability that differ between regions, as a possible reflection of deforestation and biogeographic histories. Deforested landscapes still maintain high neutral genetic diversity due to gene flow over short distances. Overall, we not only support previous evidence with microsatellite markers, but also show that deforestation can influence the genetic variability outlier, in the scenario of selective pressures imposed by these stressful environments. Based on our findings, we suggest that, to protect genetic diversity in the long term, it is necessary to reforest and enrich deforested areas, using seeds from populations in the same management target region.

Chloroquine resistance evolution in Plasmodium falciparum is mediated by the putative amino acid transporter AAT1

Malaria parasites break down host haemoglobin into peptides and amino acids in the digestive vacuole for export to the parasite cytoplasm for growth: interrupting this process is central to the mode of action of several antimalarial drugs. Mutations in the chloroquine (CQ) resistance transporter, pfcrt, located in the digestive vacuole membrane, confer CQ resistance in Plasmodium falciparum, and typically also affect parasite fitness. However, the role of other parasite loci in the evolution of CQ resistance is unclear. Here we use a combination of population genomics, genetic crosses and gene editing to demonstrate that a second vacuolar transporter plays a key role in both resistance and compensatory evolution. Longitudinal genomic analyses of the Gambian parasites revealed temporal signatures of selection on a putative amino acid transporter (pfaat1) variant S258L, which increased from 0% to 97% in frequency between 1984 and 2014 in parallel with the pfcrt1 K76T variant. Parasite genetic crosses then identified a chromosome 6 quantitative trait locus containing pfaat1 that is selected by CQ treatment. Gene editing demonstrated that pfaat1 S258L potentiates CQ resistance but at a cost of reduced fitness, while pfaat1 F313S, a common southeast Asian polymorphism, reduces CQ resistance while restoring fitness. Our analyses reveal hidden complexity in CQ resistance evolution, suggesting that pfaat1 may underlie regional differences in the dynamics of resistance evolution, and modulate parasite resistance or fitness by manipulating the balance between both amino acid and drug transport.

Genetic diversity and signatures of selection in four indigenous horse breeds of Iran

Indigenous Iranian horse breeds were evolutionarily affected by natural and artificial selection in distinct phylogeographic clades, which shaped their genomes in several unique ways. The aims of this study were to evaluate the genetic diversity and genomewide selection signatures in four indigenous Iranian horse breeds. We evaluated 169 horses from Caspian (n = 21), Turkmen (n = 29), Kurdish (n = 67), and Persian Arabian (n = 52) populations, using genomewide genotyping data. The contemporary effective population sizes were 59, 98, 102, and 113 for Turkmen, Caspian, Persian Arabian, and Kurdish breeds, respectively. By analysis of the population genetic structure, we classified the north breeds (Caspian and Turkmen) and west/southwest breeds (Persian Arabian and Kurdish) into two phylogeographic clades reflecting their geographic origin. Using the de-correlated composite of multiple selection signal statistics based on pairwise comparisons, we detected a different number of significant SNPs under putative selection from 13 to 28 for the six pairwise comparisons (FDR < 0.05). The identified SNPs under putative selection coincided with genes previously associated with known QTLs for morphological, adaptation, and fitness traits. Our results showed HMGA2 and LLPH as strong candidate genes for height variation between Caspian horses with a small size and the other studied breeds with a medium size. Using the results of studies on human height retrieved from the GWAS catalog, we suggested 38 new putative candidate genes under selection. These results provide a genomewide map of selection signatures in the studied breeds, which represent valuable information for formulating genetic conservation and improved breeding strategies for the breeds.

Comparative genomic analysis uncovers candidate genes related with milk production and adaptive traits in goat breeds

During the process of animal domestication, both natural and artificial selection cause variation in allele frequencies among populations. Identifying genomic areas of selection in domestic animals may aid in the detection of genomic areas linked to ecological and economic traits. We studied genomic variation in 140 worldwide goat individuals, including 75 Asian, 30 African and 35 European goats. We further carried out comparative population genomics to detect genomic regions under selection for adaptability to harsh conditions in local Asian ecotypes and also milk production traits in European commercial breeds. In addition, we estimated the genetic distances among 140 goat individuals. The results showed that among all studied goat groups, local breeds from West and South Asia emerged as an independent group. Our search for selection signatures in local goats from West and South Asia revealed candidate genes related to adaptation to hot climate (HSPB6, HSF4, VPS13A and NBEA genes) and immune response (IL7, IL5, IL23A and LRFN5) traits. Furthermore, selection signatures in European commercial goats involved several milk production related genes, such as VPS13C, NCAM2, TMPRSS15, CSN3 and ABCG2. The identified candidate genes could be the fundamental genetic resource for enhancement of goat production and environmental-adaptive traits, and as such they should be used in goat breeding programs to select more efficient breeds.

Selection-driven trait loss in independently evolved cavefish populations

Laboratory studies have demonstrated that a single phenotype can be produced by many different genotypes; however, in natural systems, it is frequently found that phenotypic convergence is due to parallel genetic changes. This suggests a substantial role for constraint and determinism in evolution and indicates that certain mutations are more likely to contribute to phenotypic evolution. Here we use whole genome resequencing in the Mexican tetra, Astyanax mexicanus, to investigate how selection has shaped the repeated evolution of both trait loss and enhancement across independent cavefish lineages. We show that selection on standing genetic variation and de novo mutations both contribute substantially to repeated adaptation. Our findings provide empirical support for the hypothesis that genes with larger mutational targets are more likely to be the substrate of repeated evolution and indicate that features of the cave environment may impact the rate at which mutations occur.

Selection and adaptive introgression guided the complex evolutionary history of the European common bean

Domesticated crops have been disseminated by humans over vast geographic areas. Common bean (Phaseolus vulgaris L.) was introduced in Europe after 1492. Here, by combining whole-genome profiling, metabolic fingerprinting and phenotypic characterisation, we show that the first common bean cultigens successfully introduced into Europe were of Andean origin, after Francisco Pizarro's expedition to northern Peru in 1529. We reveal that hybridisation, selection and recombination have shaped the genomic diversity of the European common bean in parallel with political constraints. There is clear evidence of adaptive introgression into the Mesoamerican-derived European genotypes, with 44 Andean introgressed genomic segments shared by more than 90% of European accessions and distributed across all chromosomes except PvChr11. Genomic scans for signatures of selection highlight the role of genes relevant to flowering and environmental adaptation, suggesting that introgression has been crucial for the dissemination of this tropical crop to the temperate regions of Europe.

Genomic analysis of two phlebotomine sand fly vectors of Leishmania from the New and Old World

Phlebotomine sand flies are of global significance as important vectors of human disease, transmitting bacterial, viral, and protozoan pathogens, including the kinetoplastid parasites of the genus Leishmania, the causative agents of devastating diseases collectively termed leishmaniasis. More than 40 pathogenic Leishmania species are transmitted to humans by approximately 35 sand fly species in 98 countries with hundreds of millions of people at risk around the world. No approved efficacious vaccine exists for leishmaniasis and available therapeutic drugs are either toxic and/or expensive, or the parasites are becoming resistant to the more recently developed drugs. Therefore, sand fly and/or reservoir control are currently the most effective strategies to break transmission. To better understand the biology of sand flies, including the mechanisms involved in their vectorial capacity, insecticide resistance, and population structures we sequenced the genomes of two geographically widespread and important sand fly vector species: Phlebotomus papatasi, a vector of Leishmania parasites that cause cutaneous leishmaniasis, (distributed in Europe, the Middle East and North Africa) and Lutzomyia longipalpis, a vector of Leishmania parasites that cause visceral leishmaniasis (distributed across Central and South America). We categorized and curated genes involved in processes important to their roles as disease vectors, including chemosensation, blood feeding, circadian rhythm, immunity, and detoxification, as well as mobile genetic elements. We also defined gene orthology and observed micro-synteny among the genomes. Finally, we present the genetic diversity and population structure of these species in their respective geographical areas. These genomes will be a foundation on which to base future efforts to prevent vector-borne transmission of Leishmania parasites.

Whole-genome sequencing of cryopreserved resources from French Large White pigs at two distinct sampling times reveals strong signatures of convergent and divergent selection between the dam and sire lines

BACKGROUND

Numerous genomic scans for positive selection have been performed in livestock species within the last decade, but often a detailed characterization of the detected regions (gene or trait under selection, timing of selection events) is lacking. Cryopreserved resources stored in reproductive or DNA gene banks offer a great opportunity to improve this characterization by providing direct access to recent allele frequency dynamics, thereby differentiating between signatures from recent breeding objectives and those related to more ancient selection constraints. Improved characterization can also be achieved by using next-generation sequencing data, which helps narrowing the size of the detected regions while reducing the number of associated candidate genes.

METHODS

We estimated genetic diversity and detected signatures of recent selection in French Large White pigs by sequencing the genomes of 36 animals from three distinct cryopreserved samples: two recent samples from dam (LWD) and sire (LWS) lines, which had diverged from 1995 and were selected under partly different objectives, and an older sample from 1977 prior to the divergence.

RESULTS

French LWD and LWS lines have lost approximately 5% of the SNPs that segregated in the 1977 ancestral population. Thirty-eight genomic regions under recent selection were detected in these lines and the corresponding selection events were further classified as convergent between lines (18 regions), divergent between lines (10 regions), specific to the dam line (6 regions) or specific to the sire line (4 regions). Several biological functions were found to be significantly enriched among the genes included in these regions: body size, body weight and growth regardless of the category, early life survival and calcium metabolism more specifically in the signatures in the dam line and lipid and glycogen metabolism more specifically in the signatures in the sire line. Recent selection on IGF2 was confirmed and several other regions were linked to a single candidate gene (ARHGAP10, BMPR1B, GNA14, KATNA1, LPIN1, PKP1, PTH, SEMA3E or ZC3HAV1, among others).

CONCLUSIONS

These results illustrate that sequencing the genome of animals at several recent time points generates considerable insight into the traits, genes and variants under recent selection in a population. This approach could be applied to other livestock populations, e.g. by exploiting the rich biological resources stored in cryobanks.

Selective footprints and genes relevant to cold adaptation and other phenotypic traits are unscrambled in the genomes of divergently selected chicken breeds

BACKGROUND

The genomes of worldwide poultry breeds divergently selected for performance and other phenotypic traits may also be affected by, and formed due to, past and current admixture events. Adaptation to diverse environments, including acclimation to harsh climatic conditions, has also left selection footprints in breed genomes.

RESULTS

Using the Chicken 50K_CobbCons SNP chip, we genotyped four divergently selected breeds: two aboriginal, cold tolerant Ushanka and Orloff Mille Fleur, one egg-type Russian White subjected to artificial selection for cold tolerance, and one meat-type White Cornish. Signals of selective sweeps were determined in the studied breeds using three methods: (1) assessment of runs of homozygosity islands, (2) F_ST based population differential analysis, and (3) haplotype differentiation analysis. Genomic regions of true selection signatures were identified by two or more methods or in two or more breeds. In these regions, we detected 540 prioritized candidate genes supplemented them with those that occurred in one breed using one statistic and were suggested in other studies. Amongst them, SOX5, ME3, ZNF536, WWP1, RIPK2, OSGIN2, DECR1, TPO, PPARGC1A, BDNF, MSTN, and beta-keratin genes can be especially mentioned as candidates for cold adaptation. Epigenetic factors may be involved in regulating some of these important genes (e.g., TPO and BDNF).

CONCLUSION

Based on a genome-wide scan, our findings can help dissect the genetic architecture underlying various phenotypic traits in chicken breeds. These include genes representing the sine qua non for adaptation to harsh environments. Cold tolerance in acclimated chicken breeds may be developed following one of few specific gene expression mechanisms or more than one overlapping response known in cold-exposed individuals, and this warrants further investigation.

Discovering novel clues of natural selection on four worldwide goat breeds

In goat breeds, the domestication followed by artificial selection for economically important traits have shaped genetic variation within populations, leading to the fixation of specific alleles for specific traits. This led to the formation and evolution of many different breeds specialised and raised for a particular purpose. However, and despite the intensity of artificial selection, natural selection continues acting, possibly leaving a more diluted contribution over time, whose traces may be more difficult to capture. In order to explore selection footprints as response of environmental adaptation, we analysed a total of 993 goats from four transboundary goats breeds (Angora, Boer, Nubian and Saanen) genotyped with the SNP chip 50 K using outlier detection, runs of homozygosity and haplotype-based detection methods. Our results showed that all methods identified footprints on chromosome 6 (from 30 to 49 Mb) for two specific populations of Nubian goats sampled in Egypt. In Angora and Saanen breeds, we detected two selective sweeps using HapFLK, on chromosome 21 (from 52 to 55 Mb) and chromosome 25 (from 1 to 5 Mb) respectively. The analysis of runs of homozygosity showed some hotspots in all breeds. The overall investigation of the selected regions detected combining the different approaches and the gene ontology exploration revealed both novel and well-known loci related to adaptation, especially for heat stress. Our findings can help to better understand the balance between the two selective pressures in commercial goat breeds providing new insights on the molecular mechanisms of adaptation.

Selection signatures for fiber production in commercial species: A review

Natural fibers derived from diverse animal species have gained increased attention in recent years due to their favorable environmental effects, long-term sustainability benefits, and remarkable physical and mechanical properties that make them valuable raw materials used for textile and non-textile production. Domestication and selective breeding for the economically significant fiber traits play an imperative role in shaping the genomes and, thus, positively impact the overall productivity of the various fiber-producing species. These selection pressures leave unique footprints on the genome due to alteration in the allelic frequencies at specific loci, characterizing selective sweeps. Recent advances in genomics have enabled the discovery of selection signatures across the genome using a variety of methods. The increased demand for 'green products' manufactured from natural fibers necessitates a detailed investigation of the genomes of the various fiber-producing plant and animal species to identify the candidate genes associated with important fiber attributes such as fiber diameter/fineness, color, length, and strength, among others. The objective of this review is to present a comprehensive overview of the concept of selection signature and selective sweeps, discuss the main methods used for its detection, and address the selection signature studies conducted so far in the diverse fiber-producing animal species.

Genomic insight into the influence of selection, crossbreeding, and geography on population structure in poultry

BACKGROUND

In poultry, the population structure of local breeds is usually complex mainly due to unrecorded breeding. Local chicken breeds offer an interesting proxy to understand the complexity of population structure in the context of human-mediated development of diverse morphologies and varieties. We studied 37 traditional Dutch chicken breeds to investigate population structure and the corresponding genomic impact using whole-genome sequence data.

RESULTS

Looking at the genetic differences between breeds, the Dutch chicken breeds demonstrated a complex and admixed subdivided structure. The dissection of this complexity highlighted the influence of selection adhering to management purposes, as well as the role of geographic distance within subdivided breed clusters. Identification of signatures of genetic differentiation revealed genomic regions that are associated with diversifying phenotypic selection between breeds, including dwarf size (bantam) and feather color. In addition, with a case study of a recently developed bantam breed developed by crossbreeding, we provide a genomic perspective on the effect of crossbreeding.

CONCLUSIONS

This study demonstrates the complex population structure of local traditional Dutch chicken, and provides insight into the genomic basis and the factors involved in the formation of this complexity.

Fast and accurate joint inference of coancestry parameters for populations and/or individuals

We introduce a fast, new algorithm for inferring from allele count data the FST parameters describing genetic distances among a set of populations and/or unrelated diploid individuals, and a tree with branch lengths corresponding to FST values. The tree can reflect historical processes of splitting and divergence, but seeks to represent the actual genetic variance as accurately as possible with a tree structure. We generalise two major approaches to defining FST, via correlations and mismatch probabilities of sampled allele pairs, which measure shared and non-shared components of genetic variance. A diploid individual can be treated as a population of two gametes, which allows inference of coancestry coefficients for individuals as well as for populations, or a combination of the two. A simulation study illustrates that our fast method-of-moments estimation of FST values, simultaneously for multiple populations/individuals, gains statistical efficiency over pairwise approaches when the population structure is close to tree-like. We apply our approach to genome-wide genotypes from the 26 worldwide human populations of the 1000 Genomes Project. We first analyse at the population level, then a subset of individuals and in a final analysis we pool individuals from the more homogeneous populations. This flexible analysis approach gives advantages over traditional approaches to population structure/coancestry, including visual and quantitative assessments of long-standing questions about the relative magnitudes of within- and between-population genetic differences.

Contrasting Phylogeographic Patterns of Mitochondrial and Genome-Wide Variation in the Groundwater Amphipod Crangonyx islandicus That Survived the Ice Age in Iceland

The analysis of phylogeographic patterns has often been based on mitochondrial DNA variation, but recent analyses dealing with nuclear DNA have in some instances revealed mito-nuclear discordances and complex evolutionary histories. These enigmatic scenarios, which may involve stochastic lineage sorting, ancestral hybridization, past dispersal and secondary contacts, are increasingly scrutinized with a new generation of genomic tools such as RADseq, which also poses additional analytical challenges. Here, we revisited the previously inconclusive phylogeographic history, showing the mito-nuclear discordance of an endemic groundwater amphipod from Iceland, Crangonyx islandicus, which is the only metazoan known to have survived the Pleistocene beneath the glaciers. Previous studies based on three DNA markers documented a mitochondrial scenario with the main divergence occurring between populations in northern Iceland and an ITS scenario with the main divergence between the south and north. We used double digest restriction-site-associated DNA sequencing (ddRADseq) to clarify this mito-nuclear discordance by applying several statistical methods while estimating the sensitivity to different analytical approaches (data-type, differentiation indices and base call uncertainty). A majority of nuclear markers and methods support the ITS divergence. Nevertheless, a more complex scenario emerges, possibly involving introgression led by male-biased dispersal among northern locations or mitochondrial capture, which may have been further strengthened by natural selection.

Signatures of selection in indigenous Chinese cattle genomes reveal adaptive genes and genetic variations to cold climate

Cold climate shapes the genome of animals and drives them to carry sufficient genetic variations to adapt to changes in temperature. However, limited information is available about the genome-wide pattern of adaptations to cold environments in cattle. In the present study, we used 777K SNP genotyping (15 Chinese cattle breeds, 198 individuals) and whole genome resequencing data (54 cattle breeds of the world, 432 individuals) to disentangle divergent selection signatures, especially between the cold-adapted (annual average temperature of habitat, 6.24 °C to 10.3 °C) and heat-adapted (20.2 °C to 24.73 °C) Chinese native cattle breeds. Genomic analyses revealed a set of candidate genes (e.g., UQCR11, DNAJC18, EGR1, and STING1) were functionally associated with thermogenesis and energy metabolism. We also characterized the adaptive loci of cattle exposed to cold temperatures. Our study finds new candidate genes and provides new insights into adaptations to cold climates in cattle.

Integrating genome-wide association studies and population genomics analysis reveals the genetic architecture of growth and backfat traits in pigs

Growth and fat deposition are complex traits, which can affect economical income in the pig industry. Due to the intensive artificial selection, a significant genetic improvement has been observed for growth and fat deposition in pigs. Here, we first investigated genomic-wide association studies (GWAS) and population genomics (e.g., selection signature) to explore the genetic basis of such complex traits in two Large White pig lines (n = 3,727) with the GeneSeek GGP Porcine HD array (n = 50,915 SNPs). Ten genetic variants were identified to be associated with growth and fatness traits in two Large White pig lines from different genetic backgrounds by performing both within-population GWAS and cross-population GWAS analyses. These ten significant loci represented eight candidate genes, i.e., NRG4, BATF3, IRS2, ANO1, ANO9, RNF152, KCNQ5, and EYA2. One of them, ANO1 gene was simultaneously identified for both two lines in BF100 trait. Compared to single-population GWAS, cross-population GWAS was less effective for identifying SNPs with population-specific effect, but more powerful for detecting SNPs with population-shared effects. We further detected genomic regions specifically selected in each of two populations, but did not observe a significant enrichment for the heritability of growth and backfat traits in such regions. In summary, the candidate genes will provide an insight into the understanding of the genetic architecture of growth-related traits and backfat thickness, and may have a potential use in the genomic breeding programs in pigs.

Candidate pigmentation genes related to feather color variation in an indigenous chicken breed revealed by whole genome data

Chicken plumage color is an inheritable phenotype that was naturally and artificially selected for during domestication. The Baicheng You chicken is an indigenous Chinese chicken breed presenting three main feather colors, lavender, black, and yellow plumages. To explore the genetic mechanisms underlying the pigmentation in Baicheng You chickens, we re-sequenced the whole genome of Baicheng You chicken with the three plumage colors. By analyzing the divergent regions of the genome among the chickens with different feather colors, we identified some candidate genomic regions associated with the feather colors in Baicheng You chickens. We found that EGR1, MLPH, RAB17, SOX5, and GRM5 genes were the potential genes for black, lavender, and yellow feathers. MLPH, GRM5, and SOX5 genes have been found to be related to plumage colors in birds. Our results showed that EGR1 is a most plausible candidate gene for black plumage, RAB17, MLPH, and SOX5 for lavender plumage, and GRM5 for yellow plumage in Baicheng You chicken.

Insights into breeding history, hotspot regions of selection, and untapped allelic diversity for bread wheat breeding

Breeding has increasingly altered the genetics of crop plants since the domestication of their wild progenitors. It is postulated that the genetic diversity of elite wheat breeding pools is too narrow to cope with future challenges. In contrast, plant genetic resources (PGRs) of wheat stored in genebanks are valuable sources of unexploited genetic diversity. Therefore, to ensure breeding progress in the future, it is of prime importance to identify the useful allelic diversity available in PGRs and to transfer it into elite breeding pools. Here, a diverse collection consisting of modern winter wheat cultivars and genebank accessions was investigated based on reduced-representation genomic sequencing and an iSelect single nucleotide polymorphism (SNP) chip array. Analyses of these datasets provided detailed insights into population structure, levels of genetic diversity, sources of new allelic diversity, and genomic regions affected by breeding activities. We identified 57 regions representing genomic signatures of selection and 827 regions representing private alleles associated exclusively with genebank accessions. The presence of known functional wheat genes, quantitative trait loci, and large chromosomal modifications, i.e., introgressions from wheat wild relatives, provided initial evidence for putative traits associated within these identified regions. These findings were supported by the results of ontology enrichment analyses. The results reported here will stimulate further research and promote breeding in the future by allowing for the targeted introduction of novel allelic diversity into elite wheat breeding pools.

Arctic introgression and chromatin regulation facilitated rapid Qinghai-Tibet Plateau colonization by an avian predator

The Qinghai-Tibet Plateau (QTP), possesses a climate as cold as that of the Arctic, and also presents uniquely low oxygen concentrations and intense ultraviolet (UV) radiation. QTP animals have adapted to these extreme conditions, but whether they obtained genetic variations from the Arctic during cold adaptation, and how genomic mutations in non-coding regions regulate gene expression under hypoxia and intense UV environment, remain largely unknown. Here, we assemble a high-quality saker falcon genome and resequence populations across Eurasia. We identify female-biased hybridization with Arctic gyrfalcons in the last glacial maximum, that endowed eastern sakers with alleles conveying larger body size and changes in fat metabolism, predisposing their QTP cold adaptation. We discover that QTP hypoxia and UV adaptations mainly involve independent changes in non-coding genomic variants. Our study highlights key roles of gene flow from Arctic relatives during QTP hypothermia adaptation, and cis-regulatory elements during hypoxic response and UV protection.

Assessing Genetic Diversity and Searching for Selection Signatures by Comparison between the Indigenous Livni and Duroc Breeds in Local Livestock of the Central Region of Russia

Indigenous pig breeds are mainly associated with the adaptive capacity that is necessary to respond adequately to climate change, food security, and livelihood needs, and natural resources conservation. Livni pigs are an indigenous fat-type breed farmed in a single farm in the Orel region and located in the Central European part of the Russian Federation. To determine the genomic regions and genes that are affected by artificial selection, we conducted the comparative study of two pig breeds with different breeding histories and breeding objectives, i.e., the native fat-type Livni and meat-type Duroc breeds using the Porcine GGP HD BeadChip, which contains ~80,000 SNPs. To check the Livni pigs for possible admixture, the Landrace and the Large White breeds were included into the study of genetic diversity as these breeds participated in the formation of the Livni pigs. We observed the highest level of genetic diversity in Livni pigs compared to commercial breeds (UHE = 0.409 vs. 0.319–0.359, p < 0.001; AR = 1.995 vs. 1.894–1.964, p < 0.001). A slight excess of heterozygotes was found in all of the breeds. We identified 291 candidate genes, which were localized within the regions under putative selection, including 22 and 228 genes, which were specific for Livni and Duroc breeds, respectively, and 41 genes common for both breeds. A detailed analysis of the molecular functions identified the genes, which were related to the formation of meat and fat traits, and adaptation to environmental stress, including extreme temperatures, which were different between breeds. Our research results are useful for conservation and sustainable breeding of Livni breed, which shows a high level of genetic diversity. This makes Livni one of the valuable national pig genetic resources.

Deciphering signatures of natural selection via deep learning

Identifying genomic regions influenced by natural selection provides fundamental insights into the genetic basis of local adaptation. However, it remains challenging to detect loci under complex spatially varying selection. We propose a deep learning-based framework, DeepGenomeScan, which can detect signatures of spatially varying selection. We demonstrate that DeepGenomeScan outperformed principal component analysis- and redundancy analysis-based genome scans in identifying loci underlying quantitative traits subject to complex spatial patterns of selection. Noticeably, DeepGenomeScan increases statistical power by up to 47.25% under nonlinear environmental selection patterns. We applied DeepGenomeScan to a European human genetic dataset and identified some well-known genes under selection and a substantial number of clinically important genes that were not identified by SPA, iHS, Fst and Bayenv when applied to the same dataset.

Natural adaptation and human selection of northeast African sheep genomes

African sheep manifest diverse but distinct physio-anatomical traits, which are the outcomes of natural- and human-driven selection. Here, we generated 34.8 million variants from 150 indigenous northeast African sheep genomes sequenced at an average depth of ∼54× for 130 samples (Ethiopia, Libya) and ∼20× for 20 samples (Sudan). These represented sheep from diverse environments, tail morphology and post-Neolithic introductions to Africa. Phylogenetic and model-based admixture analysis provided evidence of four genetic groups corresponding to altitudinal geographic origins, tail morphotypes and possible historical introduction and dispersal of the species into and across the continent. Running admixture at higher levels of K (6 ≤ K ≤ 25), revealed cryptic levels of genome intermixing as well as distinct genetic backgrounds in some populations. Comparative genomic analysis identified targets of selection that spanned conserved haplotype structures overlapping clusters of genes and gene families. These were related to hypoxia responses, ear morphology, caudal vertebrae and tail skeleton length, and tail fat-depot structures. Our findings provide novel insights underpinning morphological variation and response to human-driven selection and environmental adaptation in African indigenous sheep.

Using comparative genomics to detect mutations regulating plumage variations in graylag (A. anser) and swan geese (A. cygnoides)

Although graylag geese (A. anser) showed similar plumages of white, grey, and white with grey patches compared to those in swan geese (A. cygnoides), it was believed the substantial molecular mechanism for plumage variations were different. To date, studies on genes responsible for diverse plumages among graylag geese were limited and causal mutations remain unknown. In this study, genomes from 57 individuals belonging to six breeds showing different plumages were sequenced at ∼10X depth. Firstly, the allele frequency differences (AFD) of variants on the scaffold394 (NW_013185915.1) between grey and white goose breeds (A. anser) was calculated and a genomic region between 768,290-779,889 bp was detected to carry candidate variants associated with plumages, including one SNP (g. 775,151G > T, ∼18.6 kb upstream of EDNRB2) found to be fixed in white geese. This region was overlapped with the one detected by the haplotype-based sweep analysis, in which significant signals defined a candidate region of 736,610-820,622 bp on the same scaffold. Results from the transcriptomic data showed that expression levels of EDNRB2 and many other melanogenesis-related genes were significantly decreased among white geese compared to that in grey geese, especially at late embryonic stages (>E15). Modifications at transcriptional levels might result in abnormal melanocyte developments and thus the white plumages when they grow up. In addition, a frameshift mutation (C > -) in exon4 of MLANA gene on scaffold176 (NW_013185876.1) was suggested as the causal mutation for sex-linked dilution phenotype in graylag geese although this requires more demonstration experiments. Together with observed white plumages caused by EDNRB2 mutations in coding regions among swan geese and chicken, our study provided new examples to study the parallel evolution.

Selection signature analyses and genome-wide association reveal genomic hotspot regions that reflect differences between breeds of horse with contrasting risk of degenerative suspensory ligament desmitis

Degenerative suspensory ligament desmitis is a progressive idiopathic condition that leads to scarring and rupture of suspensory ligament fibers in multiple limbs in horses. The prevalence of degenerative suspensory ligament desmitis is breed related. Risk is high in the Peruvian Horse, whereas pony and draft breeds have low breed risk. Degenerative suspensory ligament desmitis occurs in families of Peruvian Horses, but its genetic architecture has not been definitively determined. We investigated contrasts between breeds with differing risk of degenerative suspensory ligament desmitis and identified associated risk variants and candidate genes. We analyzed 670k single nucleotide polymorphisms from 10 breeds, each of which was assigned one of the four breed degenerative suspensory ligament desmitis risk categories: control (Belgian, Icelandic Horse, Shetland Pony, and Welsh Pony), low risk (Lusitano, Arabian), medium risk (Standardbred, Thoroughbred, Quarter Horse), and high risk (Peruvian Horse). Single nucleotide polymorphisms were used for genome-wide association and selection signature analysis using breed-assigned risk levels. We found that the Peruvian Horse is a population with low effective population size and our breed contrasts suggest that degenerative suspensory ligament desmitis is a polygenic disease. Variant frequency exhibited signatures of positive selection across degenerative suspensory ligament desmitis breed risk groups on chromosomes 7, 18, and 23. Our results suggest degenerative suspensory ligament desmitis breed risk is associated with disturbances to suspensory ligament homeostasis where matrix responses to mechanical loading are perturbed through disturbances to aging in tendon (PIN1), mechanotransduction (KANK1, KANK2, JUNB, SEMA7A), collagen synthesis (COL4A1, COL5A2, COL5A3, COL6A5), matrix responses to hypoxia (PRDX2), lipid metabolism (LDLR, VLDLR), and BMP signaling (GREM2). Our results do not suggest that suspensory ligament proteoglycan turnover is a primary factor in disease pathogenesis.

Genetic and Genomic Characterization of a New Beef Cattle Composite Breed (Purunã) Developed for Production in Pasture-Based Systems

Purunã is a composite beef cattle breed, developed in Southern Brazil by crossing the Angus, Charolais, Canchim, and Caracu breeds. The goal of this study was to perform the first genetic characterization of the Purunã breed, based on both pedigree and genomic information. For this, 100 randomly selected animals were genotyped, and 11,205 animals born from 1997 to 2019 had pedigree information. The genetic analyses performed were principal component analysis, admixture, phylogenic tree, pedigree and genomic inbreeding, linkage disequilibrium (LD), effective population size (Ne), consistency of the gametic phase, runs of homozygosity (ROH), heterozygosity-enriched regions (HERs), and functional analyses of the ROH and HER regions identified. Our findings indicate that Purunã is more genetically related to the Charolais, Canchim, and Angus breeds than Caracu or Nellore. The levels of inbreeding were shown to be small based on all the metrics evaluated and ranged from −0.009 to 0.029. A low (−0.12–0.31) correlation of the pedigree-based inbreeding compared to all the genomic inbreeding coefficients evaluated was observed. The LD average was 0.031 (±0.0517), and the consistency of the gametic phase was shown to be low for all the breed pairs, ranging from 0.42 to 0.27 to the distance of 20 Mb. The Ne values based on pedigree and genomic information were 158 and 115, respectively. A total of 1,839 ROHs were found, and the majority of them are of small length (&lt;4 Mb). An important homozygous region was identified on BTA5 with pathways related to behavioral traits (sensory perception, detection of stimulus, and others), as well as candidate genes related to heat tolerance (MY O 1A), feed conversion rate (RDH5), and reproduction (AMDHD1). A total of 1,799 HERs were identified in the Purunã breed with 92.3% of them classified within the 0.5–1 Mb length group, and 19 HER islands were identified in the autosomal genome. These HER islands harbor genes involved in growth pathways, carcass weight (SDCBP), meat and carcass quality (MT2A), and marbling deposition (CISH). Despite the genetic relationship between Purunã and the founder breeds, a multi-breed genomic evaluation is likely not feasible due to their population structure and low consistency of the gametic phase among them.

Genetic architecture of a composite beef cattle population

Composite breeds are widely used in the beef industry. Composites allow producers to combine desirable traits from the progenitor breeds and simplify herd management, without repeated crossbreeding and maintenance of purebreds. In this study, genomic information was used to evaluate the genetic composition and characteristics of a three-breed beef cattle composite. This composite population referred to as Composite Gene Combination (CGC) consisted of 50% Red Angus, 25% Charolais, 25% Tarentaise. A total of 248 animals were used in this study CGC (n=79), Red Angus (n=61), Charolais (n=79) and Tarentaise (n=29). All animals were genotyped with 777k HD panel. Principal component and ADMIXTURE analyses were carried out to evaluate the genetic structure of CGC animals. The ADMIXTURE revealed the proportion of Tarentaise increased to approximately 57% while Charolais decreased to approximately 5%, and Red Angus decreased to 38% across generations. To evaluate these changes in the genomic composition across different breeds and in CGC across generations runs of homozygosity (ROH) were conducted. This analysis showed Red Angus to have the highest total length of ROH segments per animal with a mean of 349.92 Mb and lowest in CGC with a mean of 141.10 Mb. Furthermore, it showed the formation of new haplotypes in CGC around the sixth generation. Selection signatures were evaluated through Fst and HapFlk analyses. Several selection sweeps in CGC were identified especially in chromosomes 5 and 14 which have previously been reported to be associated with coat color and growth traits. The study supports our previous findings that progenitor combinations are not stable over generations and that either direct or natural selection plays a role in modifying the progenitor proportions. Furthermore, the results showed that Tarentaise contributed useful attributes to the composite in a cool semi-arid environment and suggests a re-exploration of this breed's role may be warranted.

Resequencing the Yaroslavl cattle genomes reveals signatures of selection and a rare haplotype on BTA28 likely to be related to breed phenotypes

The genomes of local livestock could shed light on their genetic history, mechanisms of adaptations to environments and unique genetics. Herein we look into the genetics and adaptations of the Russian native dairy Yaroslavl cattle breed using 22 resequenced individuals and comparing them with two related breeds (Russian Kholmogory and Holstein), and to the taurine set of the 1000 Bull Genomes Project (Run 9). HapFLK analysis with Kholmogory and Holstein breeds (using Yakut cattle as outgroup) resulted in 22 regions under selection (q‐value < 0.01) on 11 chromosomes assigned to Yaroslavl cattle, including a strong signature of selection in the region of the KIT gene on BTA6. The F_ST (fixation index) with the 1000 Bull Genomes Dataset showed 48 non‐overlapping top (0.1%) F_ST regions of which three overlapped HapFLK regions. We identified 1982 highly differentiated (F_ST > 0.40) missense mutations in the Yaroslavl genomes. These genes were enriched in the epidermal growth factor and calcium‐binding functional categories. The top F_ST intervals contained eight genes with allele frequencies quite different between the Yaroslavl and Kholmogory breeds and the rest of the 1000 Bull Genomes Dataset, including KAT6B, which had a nearly Yaroslavl breed‐specific deleterious missense mutation with the highest F_ST in our dataset (0.99). This gene is a part of a long haplotype containing other genes from F_ST and hapFLK analyses and with a negative association with weight and carcass traits according to the genotyping of 30 phenotyped Yaroslavl cattle individuals. Our work provides the industry with candidate genetic variants to be focused on in breed improvement efforts.

Early life microbial exposures shape the Crassostrea gigas immune system for lifelong and intergenerational disease protection

BACKGROUND

The interaction of organisms with their surrounding microbial communities influences many biological processes, a notable example of which is the shaping of the immune system in early life. In the Pacific oyster, Crassostrea gigas, the role of the environmental microbial community on immune system maturation - and, importantly, protection from infectious disease - is still an open question.

RESULTS

Here, we demonstrate that early life microbial exposure durably improves oyster survival when challenged with the pathogen causing Pacific oyster mortality syndrome (POMS), both in the exposed generation and in the subsequent one. Combining microbiota, transcriptomic, genetic, and epigenetic analyses, we show that the microbial exposure induced changes in epigenetic marks and a reprogramming of immune gene expression leading to long-term and intergenerational immune protection against POMS.

CONCLUSIONS

We anticipate that this protection likely extends to additional pathogens and may prove to be an important new strategy for safeguarding oyster aquaculture efforts from infectious disease. tag the videobyte/videoabstract in this section Video Abstract.

Detection of common copy number of variation underlying selection pressure in Middle Eastern horse breeds using whole-genome sequence data

Dareshouri, Arabian, and Akhal-Teke are three Middle Eastern horse breeds that have been selected for endurance and adaptation to harsh climates. Deciphering the genetic characteristics of these horses by tracing selection footprints and copy number of variations will be helpful in improving our understanding of equine breeds' development and adaptation. For this purpose, we sequenced the whole-genome of four Dareshouri horses using Illumina Hiseq panels and compared them with publicly available whole-genome sequences of Arabian (n=3) and Akhal-Teke (n=3) horses . Three tests of FLK, hapFLK, and pooled heterozygosity were applied using a sliding window (window size=100kb, step size=50kb) approach to detect putative selection signals. Copy number variation analysis was applied to investigate copy number of variants (CNVs), and the results were used to suggest selection signatures involving CNVs. Whole-genome sequencing demonstrated 8,837,950 single nucleotide polymorphisms (SNPs) in autosomal chromosomes. We suggested 58 genes and three quantitative trait loci (QTLs), including some related to horse gait, insect bite hypersensitivity, and withers height, based on selective signals detected by adjusted p-value of Mahalanobis distance based on the rank-based P-values (Md-rank-P) method. We proposed 12 genomic regions under selection pressure involving CNVs which were previously reported to be associated with metabolism energy (SLC5A8), champagne dilution in horses (SLC36A1), and synthesis of polyunsaturated fatty acids (FAT2). Only 10 Middle Eastern horses were tested in this study; therefore, the conclusions are speculative. Our findings are useful to better understanding the evolution and adaptation of Middle Eastern horse breeds.

Selection Signature Analyses Revealed Genes Associated With Adaptation, Production, and Reproduction in Selected Goat Breeds in Kenya

Artificial and natural selection in livestock is expected to leave unique footprints on their genomes. Goat breeds in Kenya have evolved for survival, breeding, and production in various harsh ecological areas, and their genomes are likely to have acquired unique alleles for adaptation to such diverse production environments and other traits of economic importance. To investigate signals of selection for some selected goat breeds in Kenya, Alpine (n = 29), Galla (n = 12), Saanen (n = 24), and Toggenburg (n = 31) were considered. A total of 53,347 single-nucleotide polymorphisms (SNPs) generated using the Illumina GoatSNP50 BeadChip were analyzed. After quality control, 47,663 autosomal single-nucleotide polymorphisms remained for downstream analyses. Several complementary approaches were applied for the following analyses: integrated Haplotype Score (iHS), cross-population-extended haplotype homozygosity (XP-EHH), hapFLK, and FLK. A total of 404 top genomic regions were identified across all the four breeds, based on the four complementary analyses. Out of the 16 identified putative selection signature regions by the intersection of multiple-selective signal analyses, most of the putative regions were found to overlap significantly with the iHS and XP-EHH analyses on chromosomes 3, 4, 10, 15, 22, and 26. These regions were enriched with some genes involved in pathways associated directly or indirectly with environmental adaptation regulating immune responses (e.g., HYAL1 and HYAL3), milk production (e.g., LEPR and PDE4B), and adaptability (e.g., MST1 and PCK). The results revealed few intersect between breeds in genomic selection signature regions. In general, this did not present the typical classic selection signatures as predicted due to the complex nature of the traits. The results support that some various selection pressures (e.g., environmental challenges, artificial selection, and genome admixture challenges) have molded the genome of goat breeds in Kenya. Therefore, the research provides new knowledge on the conservation and utilization of these goat genetic resources in Kenya. In-depth research is needed to detect precise genes connected with adaptation and production in goat breeds in Kenya.

Resequencing and Signatures of Selective Scans Point to Candidate Genetic Variants for Hair Length Traits in Long-Haired and Normal-Haired Tianzhu White Yak

Tianzhu white yak is a rare local yak breed with a pure white coat in China. In recent years, breeders have discovered long-haired individuals characterized by long hair on the forehead in the Tianzhu white yak, and the length and density of the hair on these two parts of the body are higher than that of the normal Tianzhu white yak. To elucidate the genetic mechanism of hair length in Tianzhu white yak, we re-sequence the whole genome of long-haired Tianzhu White yak (LTWY) (n = 10) and normal Tianzhu White yak (NTWY) (n = 10). Then, fixation index (F ST), θπ ratio, cross-population composite likelihood ratio (XP-CLR), integrated haplotype score (iHS), cross-population extended haplotype homozygosity (XP-EHH), and one composite method, the de-correlated composite of multiple signals (DCMS) were performed to discover the loci and genes related to long-haired traits. Based on five single methods, we found two hotspots of 0.2 and 1.1 MB in length on chromosome 6, annotating two (FGF5, CFAP299) and four genes (ATP8A1, SLC30A9, SHISA3, TMEM33), respectively. Function enrichment analysis of genes in two hotspots revealed Ras signaling pathway, MAPK signaling pathway, PI3K-Akt signaling pathway, and Rap1 signaling pathway were involved in the process of hair length differences. Besides, the DCMS method further found that four genes (ACOXL, PDPK1, MAGEL2, CDH1) were associated with hair follicle development. Henceforth, our work provides novel genetic insights into the mechanisms of hair growth in the LTWY.

Selection signatures for heat tolerance in Brazilian horse breeds

Since domestication, horse breeds have adapted to their environments and differentiated from one another. This paper uses two methods to detect selection signatures in 23 horse breeds, eight of which are Brazilian (610 animals), both cold-blooded and warm-blooded, from temperate and tropical regions. These animals were genotyped using the GGP Equine BeadChip and we analysed the data by Principal Component Analysis (PCA). The samples were separated into groups based on their geographical area of origin and PCA results studied. The genomic regions under selection were detected by hapFLK and PCAdapt methodologies, identifying six regions under selection with at least one Brazilian horse breed. These regions contain genes associated with heat tolerance, skin colour, body size, energy production/metabolism, genes involved in protein degradation/turnover/DNA repair, genes reducing the impact of oxidative stress/cellular repair, and transcriptional regulation. This work confirmed LCORL and NCAPG gene regions in previous studies associated with body size on Equine Chromosome Autosome 3 (ECA3). On the same ECA3, a region implicating genes linked to coat colour was identified, also previously related to heat stress. Regions with genes coding heat shock proteins were found on ECA1 and 2, and many candidate genes for oxidation-reduction which are a natural response to heat stress. However, a larger sample size and whole-genome SNPs are needed to understand better and identify new candidate regions as well as their functional relation with heat tolerance.

ddRAD-seq reveals the genetic structure and detects signals of selection in Italian brown trout

Background

Brown trout is one of the most widespread fresh-water fish species in Europe. The evolutionary history of and phylogenetic relationships between brown trout populations are complex, and this is especially true for Italian populations, which are heavily influenced in different ways by stocking practices. The characterization of the genetic structure of Italian brown trout populations may give information on the risk of losing endemic Italian populations due to lack of genetic diversity or to admixture with stocking populations. The identification of signatures of selection, and the information deriving from dense genotyping data will help genotype-informed breeding programs. We used a ddRAD-seq approach to obtain more than 100,000 single nucleotide polymorphisms (SNPs), and to characterize the population structure and signatures of selection in 90 brown trout samples.

Results

Italian brown trout populations are genetically differentiated, although the stocking practices have introduced strong admixture in endemic Italian trout, especially with the Atlantic lineage. Most of the analysed populations showed high levels of kinship and inbreeding. We detected putative signatures of selection using different approaches, and investigated if the regions were enriched for functional categories. Several regions putatively under selection and characterized by a reduction in heterozygosity across all the studied populations are enriched for genes involved in the response to viral infections.

Conclusions

Our results, which show evidence of admixture with the Atlantic lineage (commonly used for stocking), confirm the need for controlling stocking practices, in order to avoid the erosion of the endemic gene pool; given the apparently high levels of kinship and inbreeding in local populations, our results also show the need to take action for increasing gene diversity. In addition, we used the genetically-distinct lineages to detect signatures of selection and we identified putative signatures of selection in several regions associated with resistance to infectious diseases. These constitute candidate regions for the study of resistance to infections in wild and farmed trout.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00698-7.

Selective signatures and high genome-wide diversity in traditional Brazilian manioc (Manihot esculenta Crantz) varieties

Knowledge about genetic diversity is essential to promote effective use and conservation of crops, because it enables farmers to adapt their crops to specific needs and is the raw material for breeding. Manioc (Manihot esculenta ssp. esculenta) is one of the world's major food crops and has the potential to help achieve food security in the context of on-going climate changes. We evaluated single nucleotide polymorphisms in traditional Brazilian manioc varieties conserved in the gene bank of the Luiz de Queiroz College of Agriculture, University of São Paulo. We assessed genome-wide diversity and identified selective signatures contrasting varieties from different biomes with samples of manioc's wild ancestor M. esculenta ssp. flabellifolia. We identified signatures of selection putatively associated with resistance genes, plant development and response to abiotic stresses that might have been important for the crop's domestication and diversification resulting from cultivation in different environments. Additionally, high neutral genetic diversity within groups of varieties from different biomes and low genetic divergence among biomes reflect the complexity of manioc's evolutionary dynamics under traditional cultivation. Our results exemplify how smallholder practices contribute to conserve manioc's genetic resources, maintaining variation of potential adaptive significance and high levels of neutral genetic diversity.

Genome resequencing reveals rapid, repeated evolution in the Colorado potato beetle

Insecticide resistance and rapid pest evolution threatens food security and the development of sustainable agricultural practices, yet the evolutionary mechanisms that allow pests to rapidly adapt to control tactics remains unclear. Here we examine how a global super-pest, the Colorado potato beetle (CPB), Leptinotarsa decemlineata, rapidly evolves resistance to insecticides. Using whole genome resequencing and transcriptomic data focused on its ancestral and pest range in North America, we assess evidence for three, non-mutually exclusive models of rapid evolution: pervasive selection on novel mutations, rapid regulatory evolution, and repeated selection on standing genetic variation. Population genomic analysis demonstrates that CPB is geographically structured, even among recently established pest populations. Pest populations exhibit similar levels of nucleotide diversity, relative to non-pest populations, and show evidence of recent expansion. Genome scans provide clear signatures of repeated adaptation across CPB populations, with especially strong evidence of selection on insecticide resistance genes in different populations. Analyses of gene expression show that constitutive upregulation of candidate insecticide resistance genes drives distinctive population patterns. CPB evolves insecticide resistance repeatedly across agricultural regions, leveraging similar genetic pathways but different genes, demonstrating a polygenic trait architecture for insecticide resistance that can evolve from standing genetic variation. Despite expectations, we do not find support for strong selection on novel mutations, or rapid evolution from selection on regulatory genes. These results suggest that integrated pest management practices must mitigate the evolution of polygenic resistance phenotypes among local pest populations, in order to maintain the efficacy and sustainability of novel control techniques.