High-throughput DNA sequence analysis elucidates novel insight into the genetic basis of adaptation in local sheep.

Understanding how evolutionary factors related to climate adaptation and human selection have influenced the genetic architecture of domesticated animals is of great interest in biology. In the current study, by using 304 whole genomes from different geographical regions (including Europe, north Africa, Southwest Asia, east Asia, west Africa, south Asia, east Africa, Australia and Turkey), We evaluate global sheep population dynamics in terms of genetic variation and population structure. We further conducted comparative population analysis to study the genetic underpinnings of climate adaption to local environments and also morphological traits. In order to identify genomic signals under selection, we applied fixation index (FST) and also nucleotide diversity (θπ) statistical measurements. Our results revealed several candidate genes on different chromosomes under selection for local climate adaptation (e.g. HOXC12, HOXC13, IRF1, FGD2 and GNAQ), body size (PDGFA, HMGA2, PDE3A) and also morphological related traits (RXFP2). The discovered candidate genes may offer newel insights into genetic underpinning of regional adaptation and commercially significant features in local sheep.

The transcriptomic insight into the differential susceptibility of African Swine Fever in inbred pigs.

African swine fever (ASF) is a global threat to animal health and food security. ASF is typically controlled by strict biosecurity, rapid diagnosis, and culling of affected herds. Much progress has been made in developing modified live virus vaccines against ASF. There is host variation in response to ASF infection in the field and under controlled conditions. To better understand the dynamics underlying this host differential morbidity, whole transcriptome profiling was carried out in twelve immunized and five sham immunized pigs. Seventeen MHC homozygous inbred Large white Babraham pigs were sampled at three time points before and after the challenge. The changes in the transcriptome profiles of infected animals were surveyed over time. In addition, the immunization effect on the host response was studied as well among the contrasts of all protection subgroups. The results showed two promising candidate genes to distinguish between recovered and non-recovered pigs after infection with a virulent African swine fever virus (ASFV) pre-infection: HTRA3 and GFPT2 (padj < 0.05). Variant calling on the transcriptome assemblies showed a two-base pair insertion into the ACOX3 gene closely located to HTRA3 that may regulate its expression as a putative genomic variant for ASF. Several significant DGEs, enriched gene ontology (GO) terms, and KEGG pathways at 1 day and 7 days post-infection, compared to the pre-infection, indicate a significant inflammation response immediately after ASF infection. The presence of the virus was confirmed by the mapping of RNA-Seq reads on two whole viral genome sequences. This was concordant with a higher virus load in the non-recovered animals 7 days post-infection. There was no transcriptome signature on the immunization at pre-infection and 1 day post-infection. More samples and data from additional clinical trials may support these findings.

Strong selection signatures for Aleutian disease tolerance acting on novel candidate genes linked to immune and cellular responses in American mink (Neogale vison).

Aleutian disease (AD) is a multi-systemic infectious disease in American mink (Neogale vison) caused by Aleutian mink disease virus (AMDV). This study aimed to identify candidate regions and genes underlying selection for response against AMDV using whole-genome sequence (WGS) data. Three case-control selection signatures studies were conducted between animals (N = 85) producing high versus low antibody levels against AMDV, grouped by counter immunoelectrophoresis (CIEP) test and two enzyme-linked immunosorbent assays (ELISA). Within each study, selection signals were detected using fixation index (FST) and nucleotide diversity (θπ ratios), and validated by cross-population extended haplotype homozygosity (XP-EHH) test. Within- and between-studies overlapping results were then evaluated. Within-studies overlapping results indicated novel candidate genes related to immune and cellular responses (e.g., TAP2, RAB32), respiratory system function (e.g., SPEF2, R3HCC1L), and reproduction system function (e.g., HSF2, CFAP206) in other species. Between-studies overlapping results identified three large segments under strong selection pressure, including two on chromosome 1 (chr1:88,770-98,281 kb and chr1:114,133-120,473) and one on chromosome 6 (chr6:37,953-44,279 kb). Within regions with strong signals, we found novel candidate genes involved in immune and cellular responses (e.g., homologous MHC class II genes, ITPR3, VPS52) in other species. Our study brings new insights into candidate regions and genes controlling AD response.

Population genetic analysis and scans for adaptation and contemporary selection footprints provide genomic insight into aus, indica and japonica rice cultivars diversification.

Following domestication, rice cultivars have been spread worldwide to different climates and have experienced selection pressures to improve desirable traits. This has resulted in diverse cultivars that display variations in phenotypic traits, such as stress tolerance, grain size, and yield. To better understand the genomic composition arising from cultivar's development and local adaptation, high-density genotypes (containing 286,183 single-nucleotide polymorphisms after the quality control) of 1284 rice cultivars of aus, indica, and temperate and tropical japonica were scanned for diversifying signatures by applying a pairwise comparison of fixation index (Fst) test. Each cultivar's population was investigated for contemporary selection using the integrated haplotype score test. Signatures of diversifying selection among the pairwise comparisons were found in genomic regions mainly involved in response to stress (pathogens, drought, heat, cold) and development and morphology of various structures, such as root, pollen, spikelet, and grain. The most significant diversification signal between indica and japonica cultivars was detected at the location of ROX2 gene. Aus with indica comparison detected the most divergent signal at important candidate genes of OsEXPA8 and OsEXPA9, whereas temperate with tropical japonica comparison resulted in two well-known candidate genes OsHCT4 and OsGpx4. Recent selection analysis detected different patterns of contemporary selection in genomic regions related to rice breeding standard criteria such as stress tolerance, seed germination, starch content, and flowering time. Our findings highlight the underlying molecular basis of adaptive divergence and propose that modern rice breeding may provide additional diversification among rice cultivars.

Genome-wide association study on abdomen depth, head width, hip width, and withers height in native cattle of Guilan (Bos indicus).

Native breeds in any country are a national capital, and their preservation is of great importance. Native Cattle of Guilan (NCG) is one of the few pure native breeds in Iran and the West Asia region. During the last decade, NCG population has decreased by more than 40%. This study aimed to identify significant single nucleotide polymorphisms (SNPs) and candidate genes associated with meat production traits in NCG using a genome-wide association study (GWAS). The blood and hair samples were collected from 72 NCG individuals and genotyped using the Illumina Bovine SNP50 chip. The results of the genomic scan showed that several SNPs were associated with abdominal depth, head width, hip width, and withers height in NCG. Several candidate genes were identified, including multiple epidermal growth factor-like domains 11 (MEGF11), Methionine Sulfoxide Reductase A (MSRA), chondroitin sulfate synthase 3 (CHSY3), Cyclin-Dependent Kinase 7 (CDK7), and Parkin (PRKN) genes, which are involved in muscle growth, meat tenderness, differentiation of fat cells, fat metabolism, and adipogenesis. These genes can contribute to meat quantity and quality in NCG. This study provided valuable insights into the genetics of NCG and the identification of effective genes associated with meat production traits. The results of this study could be used for the preservation and sustainable use of this breed of native cattle, as an important genetic resource in Iran.

Epidemiology, pathogenesis, and diagnosis of Aleutian disease caused by Aleutian mink disease virus: A literature review with a perspective of genomic breeding for disease control in American mink (Neogale vison).

Aleutian disease (AD) is a multi-systemic infectious disease in American mink (Neogale vison) caused by the Aleutian mink disease virus (AMDV). Commonly referred to as mink plasmacytosis, AD is an economically significant disease in mink-breeding countries. Aleutian disease mainly induces weight loss, lower fertility, and dropped pelt quality in adults and can result in acute interstitial pneumonia with high mortality rates in kits. In this review, we employed the scientific literature on AD over the last 70 years to discuss the historical and contemporary status of AD outbreaks and seroprevalence in mink farming countries. We also explained different forms of AD and the differences between the pathogenicity of the virus in kits and adults. The application of the available AD serological tests in AD control strategies was argued. We explained how selection programs could help AD control and proposed different approaches to selecting animals for building AD-tolerant herds. The advantages of genomic selection for AD tolerance over traditional breeding strategies were discussed in detail. We also explained how genomic selection could help AD control by selecting tolerant animals for the next generation based on genome-wide single nucleotide polymorphisms (SNP) data and the challenges of implementing genomic selection for AD tolerance in the mink industry. This review collected the information required for designing successful breeding programs for AD tolerance. Examples of the application of information are presented, and data gaps are highlighted. We showed that AD tolerance is necessary to be among the traits that animals are selected for in the mink industry.

Expanding the application of haplotype-based genomic predictions to the wild: A case of antibody response against Teladorsagia circumcincta in Soay sheep.

BACKGROUND: Genomic prediction of breeding values (GP) has been adopted in evolutionary genomic studies to uncover microevolutionary processes of wild populations or improve captive breeding strategies. While recent evolutionary studies applied GP with individual single nucleotide polymorphism (SNP), haplotype-based GP could outperform individual SNP predictions through better capturing the linkage disequilibrium (LD) between the SNP and quantitative trait loci (QTL). This study aimed to evaluate the accuracy and bias of haplotype-based GP of immunoglobulin (Ig) A (IgA), IgE, and IgG against Teladorsagia circumcincta in lambs of an unmanaged sheep population (Soay breed) based on Genomic Best Linear Unbiased Prediction (GBLUP) and five Bayesian [BayesA, BayesB, BayesCπ, Bayesian Lasso (BayesL), and BayesR] methods. RESULTS: The accuracy and bias of GPs using SNP, haplotypic pseudo-SNP from blocks with different LD thresholds (0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.00), or the combinations of pseudo-SNPs and non-LD clustered SNPs were obtained. Across methods and marker sets, higher ranges of genomic estimated breeding values (GEBV) accuracies were observed for IgA (0.20 to 0.49), followed by IgE (0.08 to 0.20) and IgG (0.05 to 0.14). Considering the methods evaluated, up to 8% gains in GP accuracy of IgG were achieved using pseudo-SNPs compared to SNPs. Up to 3% gain in GP accuracy for IgA was also obtained using the combinations of the pseudo-SNPs with non-clustered SNPs in comparison to fitting individual SNP. No improvement in GP accuracy of IgE was observed using haplotypic pseudo-SNPs or their combination with non-clustered SNPs compared to individual SNP. Bayesian methods outperformed GBLUP for all traits. Most scenarios yielded lower accuracies for all traits with an increased LD threshold. GP models using haplotypic pseudo-SNPs predicted less-biased GEBVs mainly for IgG. For this trait, lower bias was observed with higher LD thresholds, whereas no distinct trend was observed for other traits with changes in LD. CONCLUSIONS: Haplotype information improves GP performance of anti-helminthic antibody traits of IgA and IgG compared to fitting individual SNP. The observed gains in the predictive performances indicate that haplotype-based methods could benefit GP of some traits in wild animal populations.

Comparison of allele-specific expression in Sistani cattle and its crossbreed with Holstein, Simmental, and Montbeliarde breeds.

The present study has analyzed the allelic-specific expression in Purebred Sistani (Bos Indicus) and their crossbreed with Holstein, Simmental, and Montbeliarde breeds (Bos Taurus). The blood samples were taken from the caudal vein of purebred Sistani cows and crossbreed Sistani's with Holstein, Simental, and Montbeliarde (4 treatments). We discovered 152,496 (Purebred Sistani), 134,285 (Sistani × Simmental), 163,362 (Sistani × Montbeliarde), and 177,042 (Sistani × Holstein) SNPs on the assembled transcriptomes. In the Purebred Sistani, 8295 (5%), Sistani × Holstein crossbreed 11,900 (7%), Sistani × Simmental crossbreed 13,187 (10%), and Sistani × Montbeliarde crossbreed 16,666 (10%) number of SNPs were identified as ASE-SNPs. In the present study, 12 SNPs types identify, of which four were transition and eight were transversion. The most common SNPs were transition types. These SNPs were present in purebred Sistani 71.84%, Sistani × Holstein crossbreed 72.65%, Sistani × Simmental crossbreed 72.60%, and Sistani × Montbeliarde crossbreed 71.94%. Ontology analysis of the expressed genes in these cows revealed the involvement of these genes in different Biological classifications. Conducting such studies in parts of the world, such as the Sistan region, where it is not possible to record accurate records of cows, is a suitable and economical method for identifying genes with different expressions.

Multi-Omics Integration and Network Analysis Reveal Potential Hub Genes and Genetic Mechanisms Regulating Bovine Mastitis.

Mastitis, inflammation of the mammary gland, is the most prevalent disease in dairy cattle that has a potential impact on profitability and animal welfare. Specifically designed multi-omics studies can be used to prioritize candidate genes and identify biomarkers and the molecular mechanisms underlying mastitis in dairy cattle. Hence, the present study aimed to explore the genetic basis of bovine mastitis by integrating microarray and RNA-Seq data containing healthy and mastitic samples in comparative transcriptome analysis with the results of published genome-wide association studies (GWAS) using a literature mining approach. The integration of different information sources resulted in the identification of 33 common and relevant genes associated with bovine mastitis. Among these, seven genes-CXCR1, HCK, IL1RN, MMP9, S100A9, GRO1, and SOCS3-were identified as the hub genes (highly connected genes) for mastitis susceptibility and resistance, and were subjected to protein-protein interaction (PPI) network and gene regulatory network construction. Gene ontology annotation and enrichment analysis revealed 23, 7, and 4 GO terms related to mastitis in the biological process, molecular function, and cellular component categories, respectively. Moreover, the main metabolic-signalling pathways responsible for the regulation of immune or inflammatory responses were significantly enriched in cytokine-cytokine-receptor interaction, the IL-17 signaling pathway, viral protein interaction with cytokines and cytokine receptors, and the chemokine signaling pathway. Consequently, the identification of these genes, pathways, and their respective functions could contribute to a better understanding of the genetics and mechanisms regulating mastitis and can be considered a starting point for future studies on bovine mastitis.

Genome-Wide Selection Signatures and Human-Mediated Introgression Events in Bos taurus indicus-influenced Composite Beef Cattle.

Genetic introgression from interbreeding hybridization of European Bos taurus taurus (EBT) and Indian Bos taurus indicus (IBI) cattle breeds have been widely used to combine the climatic resilience of the IBI cattle and the higher productivity of EBT when forming new composite beef cattle (CB) populations. The subsequent breeding strategies have shifted their initial genomic compositions. To uncover population structure, signatures of selection, and potential introgression events in CB populations, high-density genotypes [containing 492,954 single nucleotide polymorphisms (SNPs) after the quality control] of 486 individuals from 15 cattle breeds, including EBT, IBI, and CB populations, along with two Bos grunniens genotypes as outgroup were used in this study. Then, in-depth population genetics analyses were performed for three CB breeds of Beefmaster, Brangus, and Santa Gertrudis. Neighbor-joining, principal components, and admixture analyses confirmed the historical introgression of EBT and IBI haplotypes into CB breeds. The fdM statistics revealed that only 12.9% of CB populations' genetic components are of IBI origin. The results of signatures of selection analysis indicated different patterns of selection signals in the three CB breeds with primary pressure on pathways involved in protein processing and stress response in Beefmaster, cell proliferation regulation and immune response in Brangus, and amino acids and glucose metabolisms in Santa Gertrudis. An average of >90% of genomic regions underlying selection signatures were of EBT origin in the studied CB populations. Investigating the CB breeds' genome allows the estimation of EBT and IBI ancestral proportions and the locations within the genome where either taurine or indicine origin alleles are under selective pressure. Such findings highlight various opportunities to control the selection process more efficiently and explore complementarity at the genomic level in CB populations.

Identification of Candidate Genes for Pigmentation in Camels Using Genotyping-by-Sequencing.

The coat color of dromedary is usually uniform and varies from black to white, although dark- to light-brown colors are the most common phenotypes. This project was designed to gain knowledge on novel color-related variants using genotyping-by-sequencing (GBS). The association between the SNPs and coat color was tested using MLM (mixed linear models) with kinship matrix. Three GWAS models including white color vs. non-white color, black vs. non-black color, and light-brown vs. dark-brown color were performed. There were no distinct genetic clusters detected based on the color phenotypes. However, admixture occurred among all individuals of the four different coat color groups. We identified nine significant SNPs associated with white color after Bonferroni correction, located close to ANKRD26, GNB1, TSPYL4, TEKT5, DEXI, CIITA, TVP23B, CLEC16A, TMPRSS13, FXYD6, MPZL3, ANKRD26, HFM1, CDC7, TGFBR3, and HACE1 genes in neighboring flanking regions. The 13 significant SNPs associated with black color and the candidate genes were: CAPN7, CHRM4, CIITA, CLEC16A, COL4A4, COL6A6, CREB3L1, DEXI, DGKZ, DGKZ, EAF1, HDLBP, INPP5F, MCMBP, MDK, SEC23IP, SNAI1, TBX15, TEKT5, TMEM189, trpS, TSPYL4, TVP23B, and UBE2V1. The SNAI1 gene interacted with MCIR, ASIP and KIT genes. These genes play a key role in the melanin biosynthetic and pigmentation biological process and melanogenesis biological pathway. Further research using a larger sample size and pedigree data will allow confirmation of associated SNPs and the identified candidate genes.

The Candidate Chromosomal Regions Responsible for Milk Yield of Cow: A GWAS Meta-Analysis.

Milk yield (MY) is highly heritable and an economically important trait in dairy livestock species. To increase power to detect candidate genomic regions for this trait, we carried out a meta-analysis of genome-wide association studies (GWAS). In the present study, we identified 19 studies in PubMed for the meta-analysis. After review of the studies, 16 studies passed the filters for meta-analysis, and the number of chromosomes, detected markers and their positions, number of animals, and p-values were extracted from these studies and recorded. The final data set based on 16 GWAS studies had 353,698 cows and 3950 markers and was analyzed using METAL software. Our findings revealed 1712 significant (p-value < 2.5 × 10−6) genomic loci related to MY, with markers associated with MY found on all autosomes and sex chromosomes and the majority of them found on chromosome 14. Furthermore, gene ontology (GO) annotation was used to explore biological functions of the genes associated with MY; therefore, different regions of this chromosome may be suitable as genomic regions for further research into gene expression.

Gene-Set Enrichment Analysis for Identifying Genes and Biological Activities Associated with Growth Traits in Dromedaries.

Growth is an important heritable economic trait for dromedaries and necessary for planning a successful breeding program. Until now, genome-wide association studies (GWAS) and QTL-mapping have identified significant single nucleotide polymorphisms (SNPs) associated with growth in domestic animals, but in dromedaries, the number of studies is very low. This project aimed to find biological themes affecting growth in dromedaries. In the first step, 99 candidate SNPs were chosen from a previously established set of SNPs associated with body weight, gain, and birth weight in Iranian dromedaries. Next, 0.5 kb upstream and downstream of each candidate SNP were selected from NCBI (assembly accession: GCA_000803125.3). The annotation of fragments with candidate SNPs regarding the reference genome was retrieved using the Blast2GO tool. Candidate SNPs associated with growth were mapped to 22 genes, and 25 significant biological themes were identified to be related to growth in dromedaries. The main biological functions included calcium ion binding, protein binding, DNA-binding transcription factor activity, protein kinase activity, tropomyosin binding, myosin complex, actin-binding, ATP binding, receptor signaling pathway via JAK-STAT, and cytokine activity. EFCAB5, MTIF2, MYO3A, TBX15, IFNL3, PREX1, and TMOD3 genes are candidates for improving growth in camel breeding programs.

Weighted Single-Step GWAS for Body Mass Index and Scans for Recent Signatures of Selection in Yorkshire Pigs.

Controlling extra fat deposition is economically favorable in modern swine industry. Understanding the genetic architecture of fat deposition traits such as body mass index (BMI) can help in improving genomic selection for such traits. We utilized a weighted single-step genome-wide association study (WssGWAS) to detect genetic regions and candidate genes associated with BMI in a Yorkshire pig population. Three extended haplotype homozygosity (EHH)-related statistics were also incorporated within a de-correlated composite of multiple signals (DCMS) framework to detect recent selection signatures signals. Overall, the full pedigree consisted of 7016 pigs, of which 5561 had BMI records and 598 pigs were genotyped with an 80 K single nucleotide polymorphism (SNP) array. Results showed that the most significant windows (top 15) explained 9.35% of BMI genetic variance. Several genes were detected in regions previously associated with pig fat deposition traits and treated as potential candidate genes for BMI in Yorkshire pigs: FTMT, SRFBP1, KHDRBS3, FOXG1, SOD3, LRRC32, TSKU, ACER3, B3GNT6, CCDC201, ADCY1, RAMP3, TBRG4, CCM2. Signature of selection analysis revealed multiple candidate genes previously associated with various economic traits. However, BMI genetic variance explained by regions under selection pressure was minimal (1.31%). In conclusion, candidate genes associated with Yorkshire pigs' BMI trait were identified using WssGWAS. Gene enrichment analysis indicated that the identified candidate genes were enriched in the insulin secretion pathway. We anticipate that these results further advance our understanding of the genetic architecture of BMI in Yorkshire pigs and provide information for genomic selection for fat deposition in this breed.

Allele Specific Expression (ASE) analysis between Bos Taurus and Bos Indicus cows using RNA-Seq data at SNP level and gene level.

In the current study, allele specific expression analysis was performed in two subspecies cows (Bos taurus and Bos indicus) at SNP and gene levels. RNA-Seq data of 21,078,477 and 20940063 paired end reads from pooling of whole blood samples (Leukocyte) from 40 US Holstein (Bos Taurus) and 45 Cholistani cows (Bos indicus) obtained from SRA database in NCBI. Quality control and trimming of row RNA-Seq data were processed by FASTQC and Trimmomatic softwares. The transcriptome was assembled by TopHat2 software in two cow's population by aligning and mapping the RNA-Seq reads on bovine reference genome. The SNPs were discovered by Samtools software and ASE analysis was performed by Chi-square test. Results showed that 50183 and 137954 SNPs were discovered on the assembled transcriptome of Holstein and Cholistani cow samples, respectively, and 15308 SNPs were common in both breeds. 10158 SNPs from 50183 (20%) in Holstein and 31523 SNPs from 137954 (23%) in Cholistani cows were identified as ASE-SNPs. Reference allele and alternative allele count in Holstein and Cholistani cows were 3041 and 7155, respectively. Among 131 discovered SNPs in 41 genes with different expression in Holstein and Cholistani cows, 31 ASE-SNPs (5 in Holstein; 26 in Cholistani cows) were discovered.

Genomic prediction for growth using a low-density SNP panel in dromedary camels.

For thousands of years, camels have produced meat, milk, and fiber in harsh desert conditions. For a sustainable development to provide protein resources from desert areas, it is necessary to pay attention to genetic improvement in camel breeding. By using genotyping-by-sequencing (GBS) method we produced over 14,500 genome wide markers to conduct a genome- wide association study (GWAS) for investigating the birth weight, daily gain, and body weight of 96 dromedaries in the Iranian central desert. A total of 99 SNPs were associated with birth weight, daily gain, and body weight (p-value < 0.002). Genomic breeding values (GEBVs) were estimated with the BGLR package using (i) all 14,522 SNPs and (ii) the 99 SNPs by GWAS. Twenty-eight SNPs were associated with birth weight, daily gain, and body weight (p-value < 0.001). Annotation of the genomic region (s) within ± 100 kb of the associated SNPs facilitated prediction of 36 candidate genes. The accuracy of GEBVs was more than 0.65 based on all 14,522 SNPs, but the regression coefficients for birth weight, daily gain, and body weight were 0.39, 0.20, and 0.23, respectively. Because of low sample size, the GEBVs were predicted using the associated SNPs from GWAS. The accuracy of GEBVs based on the 99 associated SNPs was 0.62, 0.82, and 0.57 for birth weight, daily gain, and body weight. This report is the first GWAS using GBS on dromedary camels and identifies markers associated with growth traits that could help to plan breeding program to genetic improvement. Further researches using larger sample size and collaboration of the camel farmers and more profound understanding will permit verification of the associated SNPs identified in this project. The preliminary results of study show that genomic selection could be the appropriate way to genetic improvement of body weight in dromedary camels, which is challenging due to a long generation interval, seasonal reproduction, and lack of records and pedigrees.

Genome-wide diversity and global migration patterns in dromedaries follow ancient caravan routes.

Dromedaries have been essential for the prosperity of civilizations in arid environments and the dispersal of humans, goods and cultures along ancient, cross-continental trading routes. With increasing desertification their importance as livestock species is rising rapidly, but little is known about their genome-wide diversity and demographic history. As previous studies using few nuclear markers found weak phylogeographic structure, here we detected fine-scale population differentiation in dromedaries across Asia and Africa by adopting a genome-wide approach. Global patterns of effective migration rates revealed pathways of dispersal after domestication, following historic caravan routes like the Silk and Incense Roads. Our results show that a Pleistocene bottleneck and Medieval expansions during the rise of the Ottoman empire have shaped genome-wide diversity in modern dromedaries. By understanding subtle population structure we recognize the value of small, locally adapted populations and appeal for securing genomic diversity for a sustainable utilization of this key desert species.

Genome-Wide Diversity, Population Structure and Demographic History of Dromedaries in the Central Desert of Iran.

The development of camel husbandry for good production in a desert climate is very important, thus we need to understand the genetic basis of camels and give attention to genomic analysis. We assessed genome-wide diversity, linkage disequilibrium (LD), effective population size (Ne) and relatedness in 96 dromedaries originating from five different regions of the central desert of Iran using genotyping-by-sequencing (GBS). A total of 14,522 Single Nucleotide Polymorphisms (SNPs) with an average minor allele frequency (MAF) of 0.19 passed quality control and filtering steps. The average observed heterozygosity in the population was estimated at 0.25 ± 0.03. The mean of LD at distances shorter than 40 kb was low (r2 = 0.089 ± 0.234). The camels sampled from the central desert of Iran exhibited higher relatedness than Sudanese and lower than Arabian Peninsula dromedaries. Recent Ne of Iran's camels was estimated to be 89. Predicted Tajima's D (1.28) suggested a bottleneck or balancing selection in dromedary camels in the central desert of Iran. A general decrease in effective and census population size poses a threat for Iran's dromedaries. This report is the first SNP calling report on nearly the chromosome level and a first step towards understanding genomic diversity, population structure and demography in Iranian dromedaries.

Signatures of selection analysis using whole-genome sequence data reveals novel candidate genes for pony and light horse types.

Natural selection and domestication have shaped modern horse populations, resulting in a vast range of phenotypically diverse breeds. Horse breeds are classified into three types (pony, light, and draft) generally based on their body type. Understanding the genetic basis of horse type variation and selective pressures related to the evolutionary trend can be particularly important for current selection strategies. Whole-genome sequences were generated for 14 pony and 32 light horses to investigate the genetic signatures of selection of the horse type in pony and light horses. In the overlapping extremes of the fixation index and nucleotide diversity results, we found novel genomic signatures of selective sweeps near key genes previously implicated in body measurements including C4ORF33, CRB1, CPN1, FAM13A, and FGF12 that may influence variation in pony and light horse types. This study contributes to a better understanding of the genetic background of differences between pony and light horse types.