scalepopgen: Bioinformatic Workflow Resources Implemented in Nextflow for Comprehensive Population Genomic Analyses.

Population genomic analyses such as inference of population structure and identifying signatures of selection usually involve the application of a plethora of tools. The installation of tools and their dependencies, data transformation, or series of data preprocessing in a particular order sometimes makes the analyses challenging. While the usage of container-based technologies has significantly resolved the problems associated with the installation of tools and their dependencies, population genomic analyses requiring multistep pipelines or complex data transformation can greatly be facilitated by the application of workflow management systems such as Nextflow and Snakemake. Here, we present scalepopgen, a collection of fully automated workflows that can carry out widely used population genomic analyses on the biallelic single nucleotide polymorphism data stored in either variant calling format files or the plink-generated binary files. scalepopgen is developed in Nextflow and can be run locally or on high-performance computing systems using either Conda, Singularity, or Docker. The automated workflow includes procedures such as (i) filtering of individuals and genotypes; (ii) principal component analysis, admixture with identifying optimal K-values; (iii) running TreeMix analysis with or without bootstrapping and migration edges, followed by identification of an optimal number of migration edges; (iv) implementing single-population and pair-wise population comparison-based procedures to identify genomic signatures of selection. The pipeline uses various open-source tools; additionally, several Python and R scripts are also provided to collect and visualize the results. The tool is freely available at https://github.com/Popgen48/scalepopgen.

Genetic basis of ear length in sheep breeds sampled across the region from the Middle East to the Alps.

Ear length in sheep (Ovis aries) shows a wide range of natural variation, from the absence of an outer ear structure (anotia), to small outer ears (microtia), to regular ear length. Up until now, the underlying genetics of this phenotype has been studied in four sheep breeds from China, Jordan and Italy. These studies revealed a broad range of genes significantly associated with ear length, potentially indicating genetic heterogeneity across breeds or geographic regions. In the current study, we performed genome-wide SNP genotyping and haplotype-based mapping, in a population of 340 individuals, to identify loci influencing ear length variation in additional sheep breeds from Slovenia, Croatia, Cyprus and Greece. Additionally, two previously described candidate variants were also genotyped in our mapping population. The mapping model without candidate variant genotypes revealed only one genome-wide significant signal, which was located next to HMX1 on OAR6. This region was previously described as being associated with ear length variation in the Altay and Awassi sheep breeds. The mapping model including the candidate duplication genotype near HMX1 as a fixed effect explained the phenotypic variance on OAR6 and revealed an additional genome-wide significant locus on OAR13 associated with ear length. Our results, combined with published evidence, suggest that a duplication in the evolutionarily conserved region near HMX1 is the major regulator of ear length in sheep breeds descended from a larger region from Central Asia, to the Middle East, Cyprus, Greece and to the Alps. This distribution suggests an ancient origin of the derived allele.

Whole-genome SNP genotyping unveils ancestral and recent introgression in wild and domestic goats.

After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats.

Quantitative trait locus for calving traits on Bos taurus autosome 18 in Holstein cattle is embedded in a complex genomic region.

Although the quantitative trait locus (QTL) on chromosome 18 (BTA18) associated with paternal calving ease and stillbirth in Holstein Friesian cattle and its cross has been known for over 20 years, to our knowledge, the exact causal genetic sequence has yet escaped identification. The aim of this study was to re-examine the region of the published QTL on BTA18 and to investigate the possible reasons behind this elusiveness. For this purpose, we carried out a combined linkage disequilibrium and linkage analysis using genotyping data of 2,697 German Holstein Friesian (HF) animals and subsequent whole-genome sequencing (WGS) data analyses and genome assembly of HF samples. We confirmed the known QTL in the 95% confidence interval of 1.089 Mbp between 58.34 and 59.43 Mbp on BTA18. Additionally, these 4 SNPs in the near-perfect linkage disequilibrium with the QTL haplotype were identified: rs381577268 (on 57,816,137 bp, C/T), rs381878735 (on 59,574,329 bp, A/T), rs464221818 (on 59,329,176 bp, C/T), and rs472502785 (on 59,345,689 bp, T/C). Search for the causal mutation using short and long-read sequences, and methylation data of the BTA18 QTL region did not reveal any candidates though. The assembly showed problems in the region, as well as an abundance of segmental duplications within and around the region. Taking the QTL of BTA18 in Holstein cattle as an example, the data presented in this study comprehensively characterize the genomic features that could also be relevant for other such elusive QTL in various other cattle breeds and livestock species as well.

Fine-mapping and identification of candidate causal genes for tail length in the Merinolandschaf breed.

Docking the tails of lambs in long-tailed sheep breeds is a common practice worldwide. But this practice is associated with pain. Breeding for a shorter tail could offer an alternative. Therefore, this study aimed to analyze the natural tail length variation in the Merinolandschaf and to identify causal alleles for the short tail phenotype segregating within long-tailed breeds. We used SNP-based association analysis and haplotype-based mapping in 362 genotyped (Illumina OvineSNP50) and phenotyped Merinolandschaf lambs. Genome-wide significant regions were capture sequenced in 48 lambs and comparatively analyzed in various long and short-tailed sheep breeds and wild sheep subspecies. Here we show a SNP located in the first exon of HOXB13 and a SINE element located in the promotor of HOXB13 as promising candidates. These results enable more precise breeding towards shorter tails, improve animal welfare by amplification of ancestral alleles and contribute to a better understanding of differential embryonic development.

Evaluation of Arabian Vascular Plant Barcodes (rbcL and matK): Precision of Unsupervised and Supervised Learning Methods towards Accurate Identification.

Arabia is the largest peninsula in the world, with >3000 species of vascular plants. Not much effort has been made to generate a multi-locus marker barcode library to identify and discriminate the recorded plant species. This study aimed to determine the reliability of the available Arabian plant barcodes (>1500; rbcL and matK) at the public repository (NCBI GenBank) using the unsupervised and supervised methods. Comparative analysis was carried out with the standard dataset (FINBOL) to assess the methods and markers' reliability. Our analysis suggests that from the unsupervised method, TaxonDNA's All Species Barcode criterion (ASB) exhibits the highest accuracy for rbcL barcodes, followed by the matK barcodes using the aligned dataset (FINBOL). However, for the Arabian plant barcode dataset (GBMA), the supervised method performed better than the unsupervised method, where the Random Forest and K-Nearest Neighbor (gappy kernel) classifiers were robust enough. These classifiers successfully recognized true species from both barcode markers belonging to the aligned and alignment-free datasets, respectively. The multi-class classifier showed high species resolution following the two classifiers, though its performance declined when employed to recognize true species. Similar results were observed for the FINBOL dataset through the supervised learning approach; overall, matK marker showed higher accuracy than rbcL. However, the lower rate of species identification in matK in GBMA data could be due to the higher evolutionary rate or gaps and missing data, as observed for the ASB criterion in the FINBOL dataset. Further, a lower number of sequences and singletons could also affect the rate of species resolution, as observed in the GBMA dataset. The GBMA dataset lacks sufficient species membership. We would encourage the taxonomists from the Arabian Peninsula to join our campaign on the Arabian Barcode of Life at the Barcode of Life Data (BOLD) systems. Our efforts together could help improve the rate of species identification for the Arabian Vascular plants.

Whole genome sequencing reveals a complex introgression history and the basis of adaptation to subarctic climate in wild sheep.

To predict species responses to anthropogenic disturbances and climate change, it is reasonable to use species with high sensitivity to such factors. Snow sheep (Ovis nivicola) could represent a good candidate for this; as the only large herbivore species adapted to the cold and alpine habitats of northeastern Siberia, it plays a crucial role in its ecosystem. Despite having an extensive geographical distribution among all ovine species, it is one of the least studied. In this study, we sequenced and analysed six genomes of snow sheep in combination with all other wild sheep species to infer key aspects of their evolutionary history and unveil the genetic basis of their adaptation to subarctic environments. Despite their large census population size, snow sheep genomes showed remarkably low heterozygosity, which could reflect the effect of isolation and historical bottlenecks that we inferred using the pairwise sequential Markovian coalescent and runs of homozygosity. F4 -statistics indicated instances of introgression involving snow sheep with argali (Ovis ammon) and Dall (Ovis dalli) sheep, suggesting that these species might have been more widespread during the Pleistocene. Furthermore, the introgressed segments, which were identified using mainly minimum relative node depth, covered genes associated with immunity, adipogenesis and morphology-related traits, representing potential targets of adaptive introgression. Genes related to mitochondrial functions and thermogenesis associated with adipose tissue were identified to be under selection. Overall, our data suggest introgression as a mechanism facilitating adaptation in wild sheep species and provide insights into the genetic mechanisms underlying cold adaptation in snow sheep.

Introgression contributes to distribution of structural variations in cattle.

Structural variations (SVs) are an important source of phenotypic diversity in cattle. Here, 72 whole genome sequences representing taurine and zebu cattle were used to identify SVs. Applying multiple approaches, 16,738 SVs were identified. A comparison against the Database of Genomic Variants archives revealed that 1575 SVs were novel in our data. A novel duplication covering the entire GALNT15 gene, was observed only in N'Dama. A duplication, which was previously reported only in zebu and associated with navel length, was also observed in N'Dama. Investigation of a novel deletion located upstream of CAST13 gene and identified only in Italian cattle and zebu, revealed its introgressed origin in the former. Overall, our data highlights how the SVs distribution in cattle is also shaped by forces such as demographical differences and gene flow. The cattle SVs of this study and its meta-data can be visualized on an interactive genome browser at https://tinyurl.com/svCowArs.

A de novo frameshift mutation in ZEB2 causes polledness, abnormal skull shape, small body stature and subfertility in Fleckvieh cattle.

Polledness in cattle is an autosomal dominant trait. Previous studies have revealed allelic heterogeneity at the polled locus and four different variants were identified, all in intergenic regions. In this study, we report a case of polled bull (FV-Polled1) born to horned parents, indicating a de novo origin of this polled condition. Using 50K genotyping and whole genome sequencing data, we identified on chromosome 2 an 11-bp deletion (AC_000159.1:g.52364063_52364073del; Del11) in the second exon of ZEB2 gene as the causal mutation for this de novo polled condition. We predicted that the deletion would shorten the protein product of ZEB2 by almost 91%. Moreover, we showed that all animals carrying Del11 mutation displayed symptoms similar to Mowat-Wilson syndrome (MWS) in humans, which is also associated with genetic variations in ZEB2. The symptoms in cattle include delayed maturity, small body stature and abnormal shape of skull. This is the first report of a de novo dominant mutation affecting only ZEB2 and associated with a genetic absence of horns. Therefore our results demonstrate undoubtedly that ZEB2 plays an important role in the process of horn ontogenesis as well as in the regulation of overall development and growth of animals.

Genomic diversity and population structure of the indigenous Greek and Cypriot cattle populations.

BACKGROUND: The indigenous cattle populations from Greece and Cyprus have decreased to small numbers and are currently at risk of extinction due to socio-economic reasons, geographic isolation and crossbreeding with commercial breeds. This study represents the first comprehensive genome-wide analysis of 10 indigenous cattle populations from continental Greece and the Greek islands, and one from Cyprus, and compares them with 104 international breeds using more than 46,000 single nucleotide polymorphisms (SNPs). RESULTS: We estimated several parameters of genetic diversity (e.g. heterozygosity and allelic diversity) that indicated a severe loss of genetic diversity for the island populations compared to the mainland populations, which is mainly due to the declining size of their population in recent years and subsequent inbreeding. This high inbreeding status also resulted in higher genetic differentiation within the Greek and Cyprus cattle group compared to the remaining geographical breed groups. Supervised and unsupervised cluster analyses revealed that the phylogenetic patterns in the indigenous Greek breeds were consistent with their geographical origin and historical information regarding crosses with breeds of Anatolian or Balkan origin. Cyprus cattle showed a relatively high indicine ancestry. Greek island populations are placed close to the root of the tree as defined by Gir and the outgroup Yak, whereas the mainland breeds share a common historical origin with Busa. Unsupervised clustering and D-statistics analyses provided strong support for Bos indicus introgression in almost all the investigated local cattle breeds along the route from Anatolia up to the southern foothills of the Alps, as well as in most cattle breeds along the Apennine peninsula to the southern foothills of the Alps. CONCLUSIONS: All investigated Cyprus and Greek breeds present complex mosaic genomes as a result of historical and recent admixture events between neighbor and well-separated breeds. While the contribution of some mainland breeds to the genetic diversity pool seems important, some island and fragmented mainland breeds suffer from a severe decline of population size and loss of alleles due to genetic drift. Conservation programs that are a compromise between what is feasible and what is desirable should focus not only on the still highly diverse mainland breeds but also promote and explore the conservation possibilities for island breeds.

The First Draft Genome Assembly of Snow Sheep (Ovis nivicola).

The snow sheep, Ovis nivicola, which is endemic to the mountain ranges of northeastern Siberia, are well adapted to the harsh cold climatic conditions of their habitat. In this study, using long reads of Nanopore sequencing technology, whole-genome sequencing, assembly, and gene annotation of a snow sheep were carried out. Additionally, RNA-seq reads from several tissues were also generated to supplement the gene prediction in snow sheep genome. The assembled genome was ∼2.62 Gb in length and was represented by 7,157 scaffolds with N50 of about 2 Mb. The repetitive sequences comprised of 41% of the total genome. BUSCO analysis revealed that the snow sheep assembly contained full-length or partial fragments of 97% of mammalian universal single-copy orthologs (n = 4,104), illustrating the completeness of the assembly. In addition, a total of 20,045 protein-coding sequences were identified using comprehensive gene prediction pipeline. Of which 19,240 (∼96%) sequences were annotated using protein databases. Moreover, homology-based searches and de novo identification detected 1,484 tRNAs; 243 rRNAs; 1,931 snRNAs; and 782 miRNAs in the snow sheep genome. To conclude, we generated the first de novo genome of the snow sheep using long reads; these data are expected to contribute significantly to our understanding related to evolution and adaptation within the Ovis genus.

Are scurs in heterozygous polled (Pp) cattle a complex quantitative trait?

BACKGROUND: Breeding genetically hornless, i.e. polled, cattle provides an animal welfare-friendly and non-invasive alternative to the dehorning of calves. However, the molecular regulation of the development of horns in cattle is still poorly understood. Studying genetic characters such as polledness and scurs, can provide valuable insights into this process. Scurs are hornlike formations that occur occasionally in a wide variety of sizes and forms as an unexpected phenotype when breeding polled cattle. METHODS: We present a unique dataset of 885 Holstein-Friesian cattle with polled parentage. The horn phenotype was carefully examined, and the phenotypic heterogeneity of the trait is described. Using a direct gene test for polledness, the polled genotype of the animals was determined. Subsequently, the existence of a putative scurs locus was investigated using high-density genotype data of a selected subset of 232 animals and two mapping approaches: mixed linear model-based association analyses and combined linkage disequilibrium and linkage analysis. RESULTS: The results of an exploratory data analysis indicated that the expression of scurs depends on age at phenotyping, sex and polled genotype. Scurs were more prevalent in males than in females. Moreover, homozygous polled animals did not express any pronounced scurs and we found that the Friesian polled allele suppresses the development of scurs more efficiently than the Celtic polled allele. Combined linkage and linkage disequilibrium mapping revealed four genome-wide significant loci that affect the development of scurs, one on BTA5 and three on BTA12. Moreover, suggestive associations were detected on BTA16, 18 and 23. The mixed linear model-based association analysis supports the results of the combined linkage and linkage disequilibrium analysis. None of the mapping approaches provided convincing evidence for a monogenic inheritance of scurs. CONCLUSIONS: Our results contradict the initial and still broadly accepted model for the inheritance of horns and scurs. We hypothesise an oligogenetic model to explain the development of scurs and polledness.

Adaptive introgression from indicine cattle into white cattle breeds from Central Italy.

Cattle domestication occurred at least twice independently and gave rise to the modern taurine and indicine cattle breeds. European cattle diversity is generally dominated by taurine cattle, although elevated levels of indicine ancestry have been recorded in several breeds from southern Europe. Here we use genome-wide high-density SNP genotyping data to investigate the taurine and indicine ancestry in southern European cattle, based on a dataset comprising 508 individuals from 23 cattle breeds of taurine, indicine and mixed ancestry, including three breeds from Central Italy known to exhibit the highest levels of indicine introgression among southern European breeds. Based on local genomic ancestry analyses, we reconstruct taurine and indicine ancestry genome-wide and along chromosomes. We scrutinise local genomic introgression signals and identify genomic regions that have introgressed from indicine into taurine cattle under positive selection, harbouring genes with functions related to body size and feed efficiency. These findings suggest that indicine-derived traits helped enhance Central Italian cattle through adaptive introgression. The identified genes could provide genomic targets for selection for improved cattle performance. Our findings elucidate the key role of adaptive introgression in shaping the phenotypic features of modern cattle, aided by cultural and livestock exchange among historic human societies.

Genomic relatedness and diversity of Swedish native cattle breeds.

BACKGROUND: Native cattle breeds are important genetic resources given their adaptation to the local environment in which they are bred. However, the widespread use of commercial cattle breeds has resulted in a marked reduction in population size of several native cattle breeds worldwide. Therefore, conservation management of native cattle breeds requires urgent attention to avoid their extinction. To this end, we genotyped nine Swedish native cattle breeds with genome-wide 150 K single nucleotide polymorphisms (SNPs) to investigate the level of genetic diversity and relatedness between these breeds. RESULTS: We used various SNP-based approaches on this dataset to connect the demographic history with the genetic diversity and population structure of these Swedish cattle breeds. Our results suggest that the Väne and Ringamåla breeds originating from southern Sweden have experienced population isolation and have a low genetic diversity, whereas the Fjäll breed has a large founder population and a relatively high genetic diversity. Based on the shared ancestry and the constructed phylogenetic trees, we identified two major clusters in Swedish native cattle. In the first cluster, which includes Swedish mountain cattle breeds, there was little differentiation among the Fjäll, Fjällnära, Swedish Polled, and Bohus Polled breeds. The second cluster consists of breeds from southern Sweden: Väne, Ringamåla and Swedish Red. Interestingly, we also identified sub-structuring in the Fjällnära breed, which indicates different breeding practices on the farms that maintain this breed. CONCLUSIONS: This study represents the first comprehensive genome-wide analysis of the genetic relatedness and diversity in Swedish native cattle breeds. Our results show that different demographic patterns such as genetic isolation and cross-breeding have shaped the genomic diversity of Swedish native cattle breeds and that the Swedish mountain breeds have retained their authentic distinct gene pool without significant contribution from any of the other European cattle breeds that were included in this study.

Deciphering the patterns of genetic admixture and diversity in southern European cattle using genome-wide SNPs.

The divergence between indicine cattle (Bos indicus) and taurine cattle (Bos taurus) is estimated to have occurred approximately 250,000 years ago, but a small number of European cattle breeds still display shared ancestry with indicine cattle. Additionally, following the divergence of African and European taurine, the gene flow between African taurine and southern European cattle has also been proposed. However, the extent to which non-European cattle ancestry is diffused across southern European cattle has not been investigated thoroughly. Also, in recent times, many local breeds have suffered severe reductions in effective population size. Therefore, in the present study, we investigated the pattern of genetic diversity in various European cattle based on single nucleotide polymorphisms (SNP) identified from whole-genome sequencing data. Additionally, we also employed unlinked and phased SNP-based approaches on high-density SNP array data to characterize non-European cattle ancestry in several southern European cattle breeds. Using heterozygosity-based parameters, we concluded that, on average, nucleotide diversity is greater in southern European cattle than western European (British and commercial) cattle. However, an abundance of long runs of homozygosity (ROH) and the pattern of Linkage disequilibrium decay suggested recent bottlenecks in Maltese and Romagnola. High nucleotide diversity outside ROH indicated a highly diverse founder population for southern European and African taurine. We also show that Iberian cattle display shared ancestry with African cattle. Furthermore, we show that Podolica is an ancient cross-bred between Indicine zebu and European taurine. Additionally, we also inferred similar ancestry profile of non-European cattle ancestry in different Balkan and Italian cattle breeds which might be an indication of the common origin of indicine ancestry in these breeds. Finally, we discuss several plausible demographic scenarios which might account for the presence of non-European cattle ancestry in these cattle breeds.

Genome-Wide Characterization of Selection Signatures and Runs of Homozygosity in Ugandan Goat Breeds.

Both natural and artificial selection are among the main driving forces shaping genetic variation across the genome of livestock species. Selection typically leaves signatures in the genome, which are often characterized by high genetic differentiation across breeds and/or a strong reduction in genetic diversity in regions associated with traits under intense selection pressure. In this study, we evaluated selection signatures and genomic inbreeding coefficients, FROH, based on runs of homozygosity (ROH), in six Ugandan goat breeds: Boer (n = 13), and the indigenous breeds Karamojong (n = 15), Kigezi (n = 29), Mubende (n = 29), Small East African (n = 29), and Sebei (n = 29). After genotyping quality control, 45,294 autosomal single nucleotide polymorphisms (SNPs) remained for further analyses. A total of 394 and 6 breed-specific putative selection signatures were identified across all breeds, based on marker-specific fixation index (FST-values) and haplotype differentiation (hapFLK), respectively. These regions were enriched with genes involved in signaling pathways associated directly or indirectly with environmental adaptation, such as immune response (e.g., IL10RB and IL23A), growth and fatty acid composition (e.g., FGF9 and IGF1), and thermo-tolerance (e.g., MTOR and MAPK3). The study revealed little overlap between breeds in genomic regions under selection and generally did not display the typical classic selection signatures as expected due to the complex nature of the traits. In the Boer breed, candidate genes associated with production traits, such as body size and growth (e.g., GJB2 and GJA3) were also identified. Furthermore, analysis of ROH in indigenous goat breeds showed very low levels of genomic inbreeding (with the mean FROH per breed ranging from 0.8% to 2.4%), as compared to higher inbreeding in Boer (mean FROH = 13.8%). Short ROH were more frequent than long ROH, except in Karamojong, providing insight in the developmental history of these goat breeds. This study provides insights into the effects of long-term selection in Boer and indigenous Ugandan goat breeds, which are relevant for implementation of breeding programs and conservation of genetic resources, as well as their sustainable use and management.

Distribution and Functionality of Copy Number Variation across European Cattle Populations.

Copy number variation (CNV), which is characterized by large-scale losses or gains of DNA fragments, contributes significantly to genetic and phenotypic variation. Assessing CNV across different European cattle populations might reveal genetic changes responsible for phenotypic differences, which have accumulated throughout the domestication history of cattle as consequences of evolutionary forces that act upon them. To explore pattern of CNVs across European cattle, we genotyped 149 individuals, that represent different European regions, using the Illumina Bovine HD Genotyping array. A total of 9,944 autosomal CNVs were identified in 149 samples using a Hidden Markov Model (HMM) as employed in PennCNV. Animals originating from several breeds of British Isles, and Balkan and Italian regions, on average, displayed higher abundance of CNV counts than Dutch or Alpine animals. A total of 923 CNV regions (CNVRs) were identified by aggregating CNVs overlapping in at least two animals. The hierarchical clustering of CNVRs indicated low differentiation and sharing of high-frequency CNVRs between European cattle populations. Various CNVRs identified in the present study overlapped with olfactory receptor genes and genes related to immune system. In addition, we also detected a CNV overlapping the Kit gene in English longhorn cattle which has previously been associated with color-sidedness. To conclude, we provide a comprehensive overview of CNV distribution in genome of European cattle. Our results indicate an important role of purifying selection and genomic drift in shaping CNV diversity that exists between different European cattle populations.

A genome-wide approach to screen for genetic variants in broilers (Gallus gallus) with divergent feed conversion ratio.

Feed conversion ratio (FCR) is an economically important trait in broilers and feed accounts for a significant proportion of the costs involved in broiler production. To explore the contribution of functional variants to FCR trait, we analyzed coding and non-coding single-nucleotide variants (SNVs) across the genome by exome sequencing in seven pairs of full-sibs broilers with divergent FCR and with a sequence coverage at an average depth of fourfold. We identified 192,119 high-quality SNVs, including 30,380 coding SNVs (cSNVs) in the experimental population. We discovered missense SNVs in PGM2, NOX4, TGFBR3, and TMX4, and synonymous SNVs in TSNAX, ITA, HSP90B1, and COL18A1 associated with FCR. Haplotype analyses of genome-wide significant SNVs in PGM2, PHKG1, DGKZ, and SOD2 were also observed with suggestive evidence of haplotype association with FCR. Single-variant and FCR QTL-related genes-based association analyses of SNVs identified newly associated genes for FCR in the regions subjected to targeted exome sequencing. The top seven SNVs were next evaluated in independent replication data sets where SNV chr. 3: 13,990,160 (c. 961G>C) at TMX4 was replicated (p < 0.05). Collectively, we have detected SNVs associated with FCR in broiler as well as identification of SNVs in known FCR QTL region. These findings should facilitate the discovery of causative variants for FCR and contribute to marker-assisted selection.