Genome-Wide Runs of Homozygosity Revealed Selection Signatures in Bos indicus

Genome-wide runs of homozygosity (ROH) are suitable for understanding population history, calculating genomic inbreeding, deciphering genetic architecture of complex traits and diseases as well as identifying genes linked with agro-economic traits. Autozygosity and ROH islands, genomic regions with elevated ROH frequencies, were characterized in 112 animals of seven Indian native cattle breeds (B. indicus) using BovineHD BeadChip. In total, 4138 ROH were detected. The average number of ROH per animal was maximum in draft breed, Kangayam (63.62 ± 22.71) and minimum in dairy breed, Sahiwal (24.62 ± 11.03). The mean ROH length was maximum in Vechur (6.97 Mb) and minimum in Hariana (4.04 Mb). Kangayam revealed the highest ROH based inbreeding (F_{ROH > 1Mb} = 0.113 ± 0.059), whereas Hariana (F_{ROH > 1Mb} = 0.042 ± 0.031) and Sahiwal (F_{ROH > 1Mb} = 0.043 ± 0.048) showed the lowest. The high standard deviation observed in each breed highlights a considerable variability in autozygosity. Out of the total autozygous segments observed in each breed except Vechur, > 80% were of short length (< 8 Mb) and contributed almost 50% of the genome proportion under ROH. However, in Vechur cattle, long ROH contributed 75% of the genome proportion under ROH. ROH patterns revealed Hariana and Sahiwal breeds as less consanguineous, while recent inbreeding was apparent in Vechur. Maximum autozygosity observed in Kangayam is attributable to both recent and ancient inbreeding. The ROH islands were harbouring higher proportion of QTLs for production traits (20.68% vs. 14.64%; P≤ 0.05) but lower for reproductive traits (11.49% vs. 15.76%; P≤ 0.05) in dairy breeds compared to draft breed. In draft cattle, genes associated with resistant to diseases/higher immunity (LYZL1, SVIL, and GPX4) and stress tolerant (CCT4) were identified in ROH islands; while in dairy breeds, for milk production (PTGFR, CSN1S1, CSN2, CSN1S2, and CSN3). Significant difference in ROH islands among large and short statured breeds was observed at chromosome 3 and 5 involving genes like PTGFR and HMGA2 responsible for milk production and stature, respectively. PCA analysis on consensus ROH regions revealed distinct clustering of dairy, draft and short stature cattle breeds.

The Global Durum Wheat Panel (GDP): An International Platform to Identify and Exchange Beneficial Alleles

Representative, broad and diverse collections are a primary resource to dissect genetic diversity and meet pre-breeding and breeding goals through the identification of beneficial alleles for target traits. From 2,500 tetraploid wheat accessions obtained through an international collaborative effort, a Global Durum wheat Panel (GDP) of 1,011 genotypes was assembled that captured 94-97% of the original diversity. The GDP consists of a wide representation of Triticum turgidum ssp. durum modern germplasm and landraces, along with a selection of emmer and primitive tetraploid wheats to maximize diversity. GDP accessions were genotyped using the wheat iSelect 90K SNP array. Among modern durum accessions, breeding programs from Italy, France and Central Asia provided the highest level of genetic diversity, with only a moderate decrease in genetic diversity observed across nearly 50 years of breeding (1970-2018). Further, the breeding programs from Europe had the largest sets of unique alleles. LD was lower in the landraces (0.4 Mbp) than in modern germplasm (1.8 Mbp) at r ² = 0.5. ADMIXTURE analysis of modern germplasm defined a minimum of 13 distinct genetic clusters (k), which could be traced to the breeding program of origin. Chromosome regions putatively subjected to strong selection pressure were identified from fixation index (F _st ) and diversity reduction index (DRI) metrics in pairwise comparisons among decades of release and breeding programs. Clusters of putative selection sweeps (PSW) were identified as co-localized with major loci controlling phenology (Ppd and Vrn), plant height (Rht) and quality (gliadins and glutenins), underlining the role of the corresponding genes as driving elements in modern breeding. Public seed availability and deep genetic characterization of the GDP make this collection a unique and ideal resource to identify and map useful genetic diversity at loci of interest to any breeding program.

New Insights into the Genetic Basis of Monge's Disease and Adaptation to High-Altitude

Human high-altitude (HA) adaptation or mal-adaptation is explored to understand the physiology, pathophysiology, and molecular mechanisms that underlie long-term exposure to hypoxia. Here, we report the results of an analysis of the largest whole-genome-sequencing of Chronic Mountain Sickness (CMS) and nonCMS individuals, identified candidate genes and functionally validated these candidates in a genetic model system (Drosophila). We used PreCIOSS algorithm that uses Haplotype Allele Frequency score to separate haplotypes carrying the favored allele from the noncarriers and accordingly, prioritize genes associated with the CMS or nonCMS phenotype. Haplotypes in eleven candidate regions, with SNPs mostly in nonexonic regions, were significantly different between CMS and nonCMS subjects. Closer examination of individual genes in these regions revealed the involvement of previously identified candidates (e.g., SENP1) and also unreported ones SGK3, COPS5, PRDM1, and IFT122 in CMS. Remarkably, in addition to genes like SENP1, SGK3, and COPS5 which are HIF-dependent, our study reveals for the first time HIF-independent gene PRDM1, indicating an involvement of wider, nonHIF pathways in HA adaptation. Finally, we observed that down-regulating orthologs of these genes in Drosophila significantly enhanced their hypoxia tolerance. Taken together, the PreCIOSS algorithm, applied on a large number of genomes, identifies the involvement of both new and previously reported genes in selection sweeps, highlighting the involvement of multiple hypoxia response systems. Since the overwhelming majority of SNPs are in nonexonic (and possibly regulatory) regions, we speculate that adaptation to HA necessitates greater genetic flexibility allowing for transcript variability in response to graded levels of hypoxia.

Twenty years of artificial directional selection have shaped the genome of the Italian Large White pig breed

In this study, we investigated at the genome-wide level if 20 years of artificial directional selection based on boar genetic evaluation obtained with a classical BLUP animal model shaped the genome of the Italian Large White pig breed. The most influential boars of this breed (n = 192), born from 1992 (the beginning of the selection program of this breed) to 2012, with an estimated breeding value reliability of >0.85, were genotyped with the Illumina Porcine SNP60 BeadChip. After grouping the boars in eight classes according to their year of birth, filtered single nucleotide polymorphisms (SNPs) were used to evaluate the effects of time on genotype frequency changes using multinomial logistic regression models. Of these markers, 493 had a PBonferroni  < 0.10. However, there was an increasing number of SNPs with a decreasing level of allele frequency changes over time, representing a continuous profile across the genome. The largest proportion of the 493 SNPs was on porcine chromosome (SSC) 7, SSC2, SSC8 and SSC18 for a total of 204 haploblocks. Functional annotations of genomic regions, including the 493 shifted SNPs, reported a few Gene Ontology terms that might underly the biological processes that contributed to increase performances of the pigs over the 20 years of the selection program. The obtained results indicated that the genome of the Italian Large White pigs was shaped by a directional selection program derived by the application of methodologies assuming the infinitesimal model that captured a continuous trend of allele frequency changes in the boar population.

Comparative selection signature analyses identify genomic footprints in Reggiana cattle, the traditional breed of the Parmigiano-Reggiano cheese production system

Reggiana is an autochthonous cattle breed reared mainly in the province of Reggio Emilia, located in the North of Italy. Reggiana cattle (originally a triple-purpose population largely diffused in the North of Italy) are characterised by a typical solid red coat colour. About 2500 cows of this breed are currently registered to its herd book. Reggiana is now considered a dual-purpose breed even if it is almost completely dedicated to the production of a mono-breed branded Protected Designation of Origin Parmigiano-Reggiano cheese, which is the main driver of the sustainable conservation of this local genetic resource. In this study, we provided the first overview of genomic footprints that characterise Reggiana and define the diversity of this local cattle breed. A total of 168 Reggiana sires (all bulls born over 35 years for which semen was available) and other 3321 sires from 3 cosmopolitan breeds (Brown, Holstein and Simmental) were genotyped with the Illumina BovineSNP50 panel. ADMIXTURE analysis suggested that Reggiana breed might have been influenced, at least in part, by the other three breeds included in this study. Selection signatures in the Reggiana genome were identified using three statistical approaches based on allele frequency differences among populations or on properties of haplotypes segregating in the populations (fixation index (FST); integrated haplotype score; cross-population extended haplotype homozygosity). We identified several regions under peculiar selection in the Reggiana breed, particularly on bovine chromosome (BTA) 6 in the KIT gene region, that is known to be involved in coat colour pattern distribution, and within the region of the LAP3, NCAPG and LCORL genes, that are associated with stature, conformation and carcass traits. Another already known region that includes the PLAG1 gene (BTA14), associated with conformation traits, showed a selection signature in the Reggiana cattle. On BTA18, a signal of selection included the MC1R gene that causes the red coat colour in cattle. Other selection sweeps were in regions, with high density of quantitative trait loci for milk production traits (on BTA20) and in several other large regions that might have contributed to shape and define the Reggiana genome (on BTA17 and BTA29). All these results, overall, indicate that the Reggiana genome might still contain several signs of its multipurpose and non-specialised utilisation, as already described for other local cattle populations, in addition to footprints derived by its ancestral origin and by its adaptation to the specialised Parmigiano-Reggiano cheese production system.

Genomic Diversity, Population Structure, and Signature of Selection in Five Chinese Native Sheep Breeds Adapted to Extreme Environments

Through long term natural and artificial selection, domestic sheep (Ovis aries) have become adapted to a diverse range of agro-ecological environments and display multiple phenotypic traits. Characterization of diversity and selection signature is essential for genetic improvement, understanding of environmental adaptation, as well as utilization and conservation of sheep genetic resources. Here, we aimed to assess genomic diversity, population structure, and genomic selection among five Chinese native sheep breeds using 600K high density SNP genotypes. A total of 96 animals of the five breeds were selected from different geographical locations with extremely dry or humid conditions. We found a high proportion of informative SNPs, ranging from 93.3% in Yabuyi to 95.5% in Wadi, Hu, and Hetian sheep. The average pairwise population differentiation (F_ST) between the breeds was 0.048%, ranging from 0.022% to 0.054%, indicating their low to moderate differentiation. PCA, ADMIXTURE, and phylogenetic tree analyses revealed a clustering pattern of the five Chinese sheep breeds according to their geographical distribution, tail type, coat color, body size, and breeding history. The genomic regions under putative selection identified by F_ST and XP-EHH approaches frequently overlapped across the breeds, and spanned genes associated with adaptation to extremely dry or humid environments, innate and adaptive immune responses, and growth, wool, milk, and reproduction traits. The present study offers novel insight into genomic adaptation to dry and humid climates in sheep among other domestic animals and provides a valuable resource for further investigation. Moreover, it contributes useful information to sustainable utilization and conservation of sheep genetic resources.

Changed Patterns of Genomic Variation Following Recent Domestication: Selection Sweeps in Farmed Atlantic Salmon

The introduction of wild Atlantic salmon into captivity, and their subsequent artificial selection for production traits, has caused phenotypic differences between domesticated fish and their wild counterparts. Identification of regions of the genome underling these changes offers the promise of characterizing the early biological consequences of domestication. In the current study, we sequenced a population of farmed European Atlantic salmon and compared the observed patterns of SNP variation to those found in conspecific wild populations. This identified 139 genomic regions that contained significantly elevated SNP homozygosity in farmed fish when compared to their wild counterparts. The most extreme was adjacent to versican, a gene involved in control of neural crest cell migration. To control for false positive signals, a second and independent dataset of farmed and wild European Atlantic salmon was assessed using the same methodology. A total of 81 outlier regions detected in the first dataset showed significantly reduced homozygosity within the second one, strongly suggesting the genomic regions identified are enriched for true selection sweeps. Examination of the associated genes identified a number previously characterized as targets of selection in other domestic species and that have roles in development, behavior and olfactory system. These include arcvf, sema6, errb4, id2-like, and 6n1-like genes. Finally, we searched for evidence of parallel sweeps using a farmed population of North American origin. This failed to detect a convincing overlap to the putative sweeps present in European populations, suggesting the factors that drive patterns of variation under domestication and early artificial selection were largely independent. This is the first analysis on domestication of aquaculture species exploiting whole-genome sequence data and resulted in the identification of sweeps common to multiple independent populations of farmed European Atlantic salmon.

Genomic analyses reveal selection footprints in rice landraces grown under on-farm conservation conditions during a short-term period of domestication

Traditional rice landraces grown under on-farm conservation conditions by indigenous farmers are extremely important for future crop improvement. However, little is known about how the natural selection and agriculture practices of indigenous farmers interact to shape and change the population genetics of rice landraces grown under on-farm conservation conditions during the domestication. In this study, we sequenced DNA from 108 core on-farm conserved rice landraces collected from the ethnic minority regions of Yunnan, China, including 56 accessions collected in 1980 and 52 accessions collected in 2007 and obtained 2,771,245 of credible SNPs. Our findings show that most genetic diversity was retained during the 27 years of domestication by on-farm conservation. However, SNPs with marked allele frequency differences were found in some genome regions, particularly enriched in genic regions, indicating changes in genic regions may have played a much more prominent role in the short-term domestication of 27 years. We identified 186 and 183 potential selective-sweep regions in the indica and japonica genomes, respectively. We propose that on-farm conserved rice landraces during the short-term domestication had a highly polygenic basis with many loci responding to selection rather than a few loci with critical changes in response to selection. Moreover, loci affecting important agronomic traits and biotic or abiotic stress responses have been particularly targeted in selection. A genome-wide association study identified 90 significant signals for six traits, 13 of which were in regions of selective sweeps. Moreover, we observed a number of significant and interesting associations between loci and environmental factors, which implies adaptation to local environment. Our results provide insights into short-term evolutionary processes and shed light on the underlying mechanisms of on-farm conservation.

Genome-wide analysis of polymorphisms identified domestication-associated long low-diversity region carrying important rice grain size/weight quantitative trait loci

Rice grain size and weight are major determinants of grain quality and yield and so have been under rigorous selection since domestication. However, the genetic basis for contrasting grain size/weight trait among Indian germplasms and their association with domestication-driven evolution is not well understood. In this study, two long (LGG) and two short grain (SGG) genotypes were resequenced. LGG (LGR and PB 1121) differentiated from SGG (Sonasal and Bindli) by 504 439 single nucleotide polymorphisms (SNPs) and 78 166 insertion-and-deletion polymorphisms. The LRK gene cluster was different and a truncation mutation in the LRK8 kinase domain was associated with LGG. Phylogeny with 3000 diverse rice accessions revealed that the four sequenced genotypes belonged to the japonica group and were at the edge of the clades indicating them to be the potential source of genetic diversity available in Indian rice germplasm. Six SNPs were significantly associated with grain size/weight and the top four of these could be validated in mapping a population, suggesting this study as a valuable resource for high-throughput genotyping. A contiguous long low-diversity region (LDR) of approximately 6 Mb carrying a major grain weight quantitative trait loci (harbouring OsTOR gene) was identified on Chromosome 5. This LDR was identified as an evolutionary important site with significant positive selection and multiple selection sweeps, and showed association with many domestication-related traits, including grain size/weight. The aus population retained more allelic variations in the LDR than the japonica and indica populations, suggesting it to be one of the divergence loci. All the data and analyses can be accessed from the RiceSzWtBase database.

Population genomic structures and signatures of selection define the genetic uniqueness of several fancy and meat rabbit breeds

Following the recent domestication process of the European rabbit (Oryctolagus cuniculus), many different breeds and lines, distinguished primarily by exterior traits such as coat colour, fur structure and body size and shape, have been constituted. In this study, we genotyped, with a high-density single-nucleotide polymorphism panel, a total of 645 rabbits from 10 fancy breeds (Belgian Hare, Champagne d'Argent, Checkered Giant, Coloured Dwarf, Dwarf Lop, Ermine, Giant Grey, Giant White, Rex and Rhinelander) and three meat breeds (Italian White, Italian Spotted and Italian Silver). ADMIXTURE analysis indicated that breeds with similar phenotypic traits (e.g. coat colour and body size) shared common ancestries. Signatures of selection using two haplotype-based approaches (iHS and XP-EHH), combined with the results obtained with other methods previously reported that we applied to the same breeds, we identified a total of 5079 independent genomic regions with some signatures of selection, covering about 1777 Mb of the rabbit genome. These regions consistently encompassed many genes involved in pigmentation processes (ASIP, EDNRA, EDNRB, KIT, KITLG, MITF, OCA2, TYR and TYRP1), coat structure (LIPH) and body size, including two major genes (LCORL and HMGA2) among many others. This study revealed novel genomic regions under signatures of selection and further demonstrated that population structures and signatures of selection, left into the genome of these rabbit breeds, may contribute to understanding the genetic events that led to their constitution and the complex genetic mechanisms determining the broad phenotypic variability present in these untapped rabbit genetic resources.

Genetic Diversity Estimation and Genome-Wide Selective Sweep Analysis of the Bazhou Yak

The Bazhou yak, a major native meat yak breed in Xinjiang, China, is renowned for its fast growth rate, strong adaptability, and particularly high intramuscular fat (IMF) content. However, limited knowledge regarding its phylogenetic history and genomic composition has hindered its long-term conservation and utilization. This study evaluated the genetic diversity, population phylogenetics, and genome-wide selective sweep analysis (GWSA) of 100 newly obtained Bazhou yaks through genome resequencing, as well as 340 public yak genomes from nine other populations on the Qinghai-Tibet Plateau. The results revealed moderate diversity, lower genomic inbreeding levels, and rapid linkage disequilibrium (LD) decay in Bazhou yaks. Principal component analysis (PCA) and phylogenetic analysis showed a clear separation of Bazhou yaks from other yak populations, indicating the Bazhou yak as an independent genetic population. Furthermore, less genetic differentiation was found between the Bazhou yak and the Huanhu yak, while ADMIXTURE analysis revealed a common ancestral lineage between Bazhou yaks and Huanhu yaks, indicating an important genetic contribution of the Qinghai yak population to Bazhou yaks. The GWSA identified a total of 833 selected genes in Bazhou yaks using the top 5% interaction windows of both parameters (FST, Pi ratio, and XP-EHH). A significant number of these genes are related to fat synthesis and deposition, such as MTOR, APOA1, and GPAT4. In summary, this study sheds light on the phylogenetic status and distinctive genomic features of Bazhou yaks, which facilitates our understanding of the genetic basis of the IMF phenotype in Bazhou yaks.

High-integrity Pueraria montana var. lobata genome and population analysis revealed the genetic diversity of Pueraria genus

Pueraria montana var. lobata (P. lobata) is a traditional medicinal plant belonging to the Pueraria genus of Fabaceae family. Pueraria montana var. thomsonii (P. thomsonii) and Pueraria montana var. montana (P. montana) are its related species. However, evolutionary history of the Pueraria genus is still largely unknown. Here, a high-integrity, chromosome-level genome of P. lobata and an improved genome of P. thomsonii were reported. It found evidence for an ancient whole-genome triplication and a recent whole-genome duplication shared with Fabaceae in three Pueraria species. Population genomics of 121 Pueraria accessions demonstrated that P. lobata populations had substantially higher genetic diversity, and P. thomsonii was probably derived from P. lobata by domestication as a subspecies. Selection sweep analysis identified candidate genes in P. thomsonii populations associated with the synthesis of auxin and gibberellin, which potentially play a role in the expansion and starch accumulation of tubers in P. thomsonii. Overall, the findings provide new insights into the evolutionary and domestication history of the Pueraria genome and offer a valuable genomic resource for the genetic improvement of these species.

Genome wide association and haplotype analyses for the crease depth trait in bread wheat (Triticum aestivum L.)

Wheat grain has a complex structure that includes a crease on one side, and tissues within the crease region play an important role in nutrient transportation during wheat grain development. However, the genetic architecture of the crease region is still unclear. In this study, 413 global wheat accessions were resequenced and a method was developed for evaluating the phenotypic data of crease depth (CD). The CD values exhibited continuous and considerable large variation in the population, and the broad-sense heritability was 84.09%. CD was found to be positively correlated with grain-related traits and negatively with quality-related traits. Analysis of differentiation of traits between landraces and cultivars revealed that grain-related traits and CD were simultaneously improved during breeding improvement. Moreover, 2,150.8-Mb genetic segments were identified to fall within the selective sweeps between the landraces and cultivars; they contained some known functional genes for quality- and grain-related traits. Genome-wide association study (GWAS) was performed using around 10 million SNPs generated by genome resequencing and 551 significant SNPs and 18 QTLs were detected significantly associated with CD. Combined with cluster analysis of gene expression, haplotype analysis, and annotated information of candidate genes, two promising genes TraesCS3D02G197700 and TraesCS5A02G292900 were identified to potentially regulate CD. To the best of our knowledge, this is the first study to provide the genetic basis of CD, and the genetic loci identified in this study may ultimately assist in wheat breeding programs.

Genomic selection analysis of morphological and adaptation traits in Chinese indigenous dog breeds

The significant morphological differences and abundant germplasm resources of Chinese indigenous dog breeds can be attributed to the diverse geographical environment, including plateaus, mountains, and a long history of raising dogs. The combination of both natural and artificial selection during the past several thousand years has led to hundreds of dog breeds with distinct morphological traits and environmental adaptations. China is one of the earliest countries to domesticate dogs and there are more than 50 ancient indigenous dog breeds. In this study, the run of homozygosity (ROH) and proportion of the autosomal genome covered by ROHs (FROH) were calculated for 10 dog breeds that are the most representative Chinese indigenous dogs based on 170K SNP microarray. The results of FROH showed that the Chuandong hound dogs (HCSSC) have the highest level of inbreeding among the tested breeds. The inbreeding in HCSSC occurred more recently than the Liangshan dogs (SCLSQ) dogs because of more numbers of long ROHs in HCSSC dogs, and the former also have higher inbreeding degree. In addition, there are significant differences in the inbreeding degree among different subpopulations of the same breed, such as the Thin dogs from Shaanxi and Shandong province. To explore genome-wide selection signatures among different breeds, including coat color, ear shape, and altitude adaptability, we performed genome selection analyses of FST and cross population extended haplotype homozygosity (XP-EHH). For the coat color, the FST analysis between Xiasi dogs (XSGZ) and HCSSC dogs was performed and identified multiple genes involved in coat color, hair follicle, and bone development, including MC1R, KITLG, SOX5, RSPO2, and TBX15. For the plateau adaptability, we performed FST and XP-EHH analyses between dogs from Tibet (Tibetan Mastiffs and Nyingchi dogs) and plain regions (Guangxi Biwei dogs GXBWQ and Guandong Sharpei dogs). The results showed the EPAS1 gene in dogs from Tibet undergo strong selection. Multiple genes identified for selection signals based on different usage of dogs. Furthermore, the results of ear shape analyses showed that MSRB3 was likely to be the main gene causing the drop ear of domestic dogs. Our study provides new insights into further understanding of Chinese indigenous dogs.