A ∼ 4.1 kb deletion in IRX1 gene upstream is completely associated with rumplessness in Piao chicken

Piao chicken, a Chinese indigenous rumpless chicken breed, lacks pygostyle, caudal vertebra, uropygial gland and tail feathers. The rumplessness in Piao chicken presents an autosomal dominant inheritance pattern. However, the molecular genetic mechanisms underlying the rumplessness in Piao chicken remains unclear. In this study, whole-genome resequencing was performed for 146 individuals from 10 chicken breeds, including 9 tailed chicken breeds and Piao rumpless breed. Tailbone CT scan for Piao chickens and WL chickens, revealed that some Piao chicken tails were normal in number, and for a few Piao chickens tail length and tail bone numbers were between the rumpless and the normal tailed chickens. The results showed that the rumpless phenotype has not been completely fixed in Piao chicken breed. Using selection signature analysis and structural variation detection, we found a 4174 bp deletion located in the upstream region of IRX1 gene on chromosome 2 related to rumpless phenotype. Structural variation genotyping showed that the deletion was present in all 32 rumpless Piao chickens (del/del, wild/del) and absent from all 112 tailed chickens included in the dataset for the other 9 breeds and 2 tailed Piao chickens (wild/wild). In summary, all rumpless Piao chickens tested here carry this deletion mutation, to show a complete linkage association with rumplessness trait. We suggested that the 4174 bp deletion could be causative for rumpless phenotype in Piao chicken since this is the only mutation to show the complete linkage disequilibrium with rumplessness on whole genome level across all of 146 chickens from the 10 breeds. This study could facilitate a better understanding of the genetic characteristics of Piao chicken.     

Genome-wide analyses reveal population structure and identify candidate genes associated with tail fatness in local sheep from a semi-arid area

Under a climate change perspective, the genetic make-up of local livestock breeds showing adaptive traits should be explored and preserved as a priority. We used genotype data from the ovine 50 k Illumina BeadChip for assessing breed autozygosity based on runs of homozygosity (ROH) and fine-scale genetic structure and for detecting genomic regions under selection in 63 Tunisian sheep samples. The average genomic inbreeding coefficients based on ROH were estimated at 0.017, 0.021, and 0.024 for Barbarine (BAR, n = 26), Noire de Thibar (NDT, n = 23), and Queue fine de l'Ouest (QFO, n = 14) breeds, respectively. The genomic relationships among individuals based on identity by state (IBS) distance matrix highlighted a recent introgression of QFO into the BAR and a genetic differentiation of NDT samples, possibly explained by past introgression of European gene pools. Genome-wide scan for ROH across breeds and within the BAR sample set identified an outstanding signal on chromosome 13 (46.58–49.61 Mbp). These results were confirmed using F_ST index, differentiating fat vs. thin-tailed individuals. Candidate genes under selection pressure (CDS2, PROKR1, and BMP2) were associated to lipid storage and probably preferentially selected in fat-tailed BAR animals. Our findings suggest paying more attention to preserve the genetic integrity and adaptive alleles of local sheep breeds.     

A Chromosome-level Genome Assembly of Wild Castor Provides New Insights into its Adaptive Evolution in Tropical Desert

Wild castor grows in the high-altitude tropical desert of the African Plateau,a region known for high ultraviolet radiation,strong light,and extremely dry condition.To investigate the potential genetic basis of adaptation to both highland and tropical deserts,we generated a chromosome-level genome sequence assembly of the wild castor accession WT05,with a genome size of 316 Mb,a scaffold N50 of 31.93 Mb,and a contig N50 of 8.96 Mb,respectively.Compared with cultivated castor and other Euphorbiac...     

Genome-Wide Detection of Selective Signatures in Chicken through High Density SNPs

Chicken is recognized as an excellent model for studies of genetic mechanism of phenotypic and genomic evolution, with large effective population size and strong human-driven selection. In the present study, we performed Extended Haplotype Homozygosity (EHH) tests to identify significant core regions employing 600K SNP Chicken chip in an F2 population of 1,534 hens, which was derived from reciprocal crosses between White Leghorn and Dongxiang chicken. Results indicated that a total of 49,151 core regions with an average length of 9.79 Kb were identified, which occupied approximately 52.15% of genome across all autosomes, and 806 significant core regions attracted us mostly. Genes in candidate regions may experience positive selection and were considered to have possible influence on beneficial economic traits. A panel of genes including AASDHPPT, GDPD5, PAR3, SOX6, GPC1 and a signal pathway of AKT1 were detected with the most extreme P-values. Further enrichment analyses indicated that these genes were associated with immune function, sensory organ development and neurogenesis, and may have experienced positive selection in chicken. Moreover, some of core regions exactly overlapped with genes excavated in our previous GWAS, suggesting that these genes have undergone positive selection may affect egg production. Findings in our study could draw a comparatively integrate genome-wide map of selection signature in the chicken genome, and would be worthy for explicating the genetic mechanisms of phenotypic diversity in poultry breeding.     

Comparative genetic mapping and genomic region collinearity analysis of the powdery mildew resistance gene Pm41

Key message

By applying comparative genomics analyses, a high-density genetic linkage map narrowed the powdery mildew resistance gene Pm41 originating from wild emmer in a sub-centimorgan genetic interval.

Abstract

Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici, results in large yield losses worldwide. A high-density genetic linkage map of the powdery mildew resistance gene Pm41, originating from wild emmer (Triticum turgidum var. dicoccoides) and previously mapped to the distal region of chromosome 3BL bin 0.63–1.00, was constructed using an F_5:6 recombinant inbred line population derived from a cross of durum wheat cultivar Langdon and wild emmer accession IW2. By applying comparative genomics analyses, 19 polymorphic sequence-tagged site markers were developed and integrated into the Pm41 genetic linkage map. Ultimately, Pm41 was mapped in a 0.6 cM genetic interval flanked by markers XWGGC1505 and XWGGC1507, which correspond to 11.7, 19.2, and 24.9 kb orthologous genomic regions in Brachypodium, rice, and sorghum, respectively. The XWGGC1506 marker co-segregated with Pm41 and could be served as a starting point for chromosome landing and map-based cloning as well as marker-assisted selection of Pm41. Detailed comparative genomics analysis of the markers flanking the Pm41 locus in wheat and the putative orthologous genes in Brachypodium, rice, and sorghum suggests that the gene order is highly conserved between rice and sorghum. However, intra-chromosome inversions and re-arrangements are evident in the wheat and Brachypodium genomic regions, and gene duplications are also present in the orthologous genomic regions of Pm41 in wheat, indicating that the Brachypodium gene model can provide more useful information for wheat marker development.     

A genome-wide perspective on the diversity and selection signatures in indigenous goats using 53 K single nucleotide polymorphism array

Tibetan goats, Taihang goats, Jining grey goats, and Meigu goats are the representative indigenous goats in China, found in Qinghai-Tibet Plateau, Western pastoral area, Northern and Southern agricultural regions. Very few studies have conducted a comprehensive analysis of the genomic diversity and selection of these breeds. We genotyped 96 unrelated individuals, using goat 53 K Illumina BeadChip array, of the following goat breeds: Tibetan (TG), Taihang (THG), Jining grey (JGG), and Meigu (MGG). A total of 45 951 single nucleotide polymorphisms were filtered to estimate the genetic diversity and selection signatures. All breeds had a high proportion (over 95%) of polymorphic loci. The observed and excepted heterozygosity ranged from 0.338 (MGG) to 0.402 (JGG) and 0.339 (MGG) to 0.395 (JGG), respectively. Clustering analysis displayed a genetically distinct lineage for each breed, and their Fst were greater than 0.25, indicating that they had a higher genetic differentiation between groups. Furthermore, effective population size reduced in all four populations, indicating a loss of genetic diversity. In addition, runs of homozygosity were mainly distributed in 5–10 Mb. Lastly, we identified signature genes, which were closely related to high-altitude adaptation (ADIRF) and prolificity (CNTROB, SMC3, and PTEN). This study provides a valuable resource for future studies on genome-wide perspectives on the diversity and selection signatures of Chinese indigenous goats.     

A combination of genome-wide association study and selection signature analysis dissects the genetic architecture underlying bone traits in chickens

The bones of chicken play an important role in supporting and protecting the body. The growth and development of bones have a substantial influence on the health and production performance in chickens. However, genetic architecture underlying chicken bone traits is not well understood. The objectives of this study are to dissect the genetic basis of bone traits in chickens and to identify valuable genes and genetic markers for chicken breeding. We performed a combination of genome-wide association study (GWAS) and selection signature analysis (fixation index values and nucleotide diversity ratios) in an F₂ crossbred experimental population with different genetic backgrounds (broiler × layer) to identify candidate genes and significant variants related to femur, shank, keel length, chest width, metatarsal claw weight, metatarsal length, and metatarsal circumference. A total of 545 individuals were genotyped based on the whole genome re-sequencing method (26 F₀ individuals were re-sequenced at 10 × coverage; 519 F₂ individuals were re-sequenced at 3 × coverage). A total of 2 028 112 single-nucleotide polymorphisms (SNPs) remained to carry out analysis after quality control and imputation. The integration of GWAS and selection signature analysis indicated that all significant SNPs responsible for bone traits were mainly localized on chicken chromosomes 1, 4, and 27. Finally, we identified 21 positional candidate genes that might regulate chicken bone growth and development, including LRCH1, RB1, FNDC3A, MLNR, CAB39L, FOXO1, LHFP, TRPC4, POSTN, SMAD9, RBPJ, PPARGC1A, SLIT2, NCAPG, NKX3-2, CPZ, SPOP, NGFR, SOST, ZNF652, and HOXB3. Additionally, an array of uncharacterized genes was identified. The findings provide an in-depth understanding of the genetic architecture of chicken bone traits and offer a molecular basis for applying genomics in practical chicken breeding.     

Genome-wide association study of economically important traits in Charolais and Limousin beef cows

Genomic selection has proven effective for advancing genetic gain for key profit traits in dairy cattle production systems. However, its impact to-date on genetic improvement programs for beef cattle has been less effective. Despite this, the technology is thought to be particularly useful for low heritability traits such as those associated with reproductive efficiency. The objective of this study was to identify genetic variants associated with key determinants of reproductive and overall productive efficiency in beef cows. The analysis employed a large dataset derived from the national genetic evaluation program in Ireland for two of the most predominant beef breeds, viz. Charolais (n = 5 244 cows) and Limousin (n = 7 304 cows). Single nucleotide polymorphisms (SNPs) were identified as being statistically significantly associated (adj. P < 0.05) with both reproductive and productive traits for both breed types. However, there was little across breed commonality, with only two SNPs (rs110240246 and rs110344317; adj. P < 0.05) located within the genomic regions of the LCORL and MSTN genes respectively, identified in both Charolais and Limousin populations, associated with traits including carcass weight, cull-cow weight and live-weight. Significant SNPs within the MSTN gene were also associated with both reproduction and production related traits within each breed. Finally, traits including calving difficulty, calf mortality and calving interval were associated with SNPs within genomic regions comprising genes involved in cellular growth and lipid metabolism. Genetic variants identified as associated with both important reproductive efficiency and production related traits from this study warrant further analyses for their potential incorporation into breeding programmes to support the sustainability of beef cattle production.     

Identification and functional characterization of copy number variations in diverse chicken breeds

Background

The detection and functional characterization of genomic structural variations are important for understanding the landscape of genetic variation in the chicken. A recently recognized aspect of genomic structural variation, called copy number variation (CNV), is gaining interest in chicken genomic studies. The aim of the present study was to investigate the pattern and functional characterization of CNVs in five characteristic chicken breeds, which will be important for future studies associating phenotype with chicken genome architecture.

Results

Using a commercial 385 K array-based comparative genomic hybridization (aCGH) genome array, we performed CNV discovery using 10 chicken samples from four local Chinese breeds and the French breed Houdan chicken. The female Anka broiler was used as a reference. A total of 281 copy number variation regions (CNVR) were identified, covering 12.8 Mb of polymorphic sequences or 1.07% of the entire chicken genome. The functional annotation of CNVRs indicated that these regions completely or partially overlapped with 231 genes and 1032 quantitative traits loci, suggesting these CNVs have important functions and might be promising resources for exploring differences among various breeds. In addition, we employed quantitative PCR (qPCR) to further validate several copy number variable genes, such as prolactin receptor, endothelin 3 (EDN3), suppressor of cytokine signaling 2, CD8a molecule, with important functions, and the results suggested that EDN3 might be a molecular marker for the selection of dark skin color in poultry production. Moreover, we also identified a new CNVR (chr24: 3484617–3512275), encoding the sortilin-related receptor gene, with copy number changes in only black-bone chicken.

Conclusions

Here, we report a genome-wide analysis of the CNVs in five chicken breeds using aCGH. The association between EDN3 and melanoblast proliferation was further confirmed using qPCR. These results provide additional information for understanding genomic variation and related phenotypic characteristics.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-934) contains supplementary material, which is available to authorized users.     

Genomic prediction for rust resistance in diverse wheat landraces

Key message

We have demonstrated that genomic selection in diverse wheat landraces for resistance to leaf, stem and strip rust is possible, as genomic breeding values were moderately accurate. Markers with large effects in the Bayesian analysis confirmed many known genes, while also discovering many previously uncharacterised genome regions associated with rust scores.

Abstract

Genomic selection, where selection decisions are based on genomic estimated breeding values (GEBVs) derived from genome-wide DNA markers, could accelerate genetic progress in plant breeding. In this study, we assessed the accuracy of GEBVs for rust resistance in 206 hexaploid wheat (Triticum aestivum) landraces from the Watkins collection of phenotypically diverse wheat genotypes from 32 countries. The landraces were genotyped for 5,568 SNPs using an Illumina iSelect 9 K bead chip assay and phenotyped for field-based leaf rust (Lr), stem rust (Sr) and stripe rust (Yr) responses across multiple years. Genomic Best Linear Unbiased Prediction (GBLUP) and a Bayesian Regression method (BayesR) were used to predict GEBVs. Based on fivefold cross-validation, the accuracy of genomic prediction averaged across years was 0.35, 0.27 and 0.44 for Lr, Sr and Yr using GBLUP and 0.33, 0.38 and 0.30 for Lr, Sr and Yr using BayesR, respectively. Inclusion of PCR-predicted genotypes for known rust resistance genes increased accuracy more substantially when the marker was diagnostic (Lr34/Sr57/Yr18) for the presence-absence of the gene rather than just linked (Sr2). Investigation of the impact of genetic relatedness between validation and reference lines on accuracy of genomic prediction showed that accuracy will be higher when each validation line had at least one close relationship to the reference lines. Overall, the prediction accuracies achieved in this study are encouraging, and confirm the feasibility of genomic selection in wheat. In several instances, estimated marker effects were confirmed by published literature and results of mapping experiments using Watkins accessions.     

Generation of selectable marker-free transgenic eggplant resistant to Alternaria solani using the R/RS site-specific recombination system

Key message

Marker-free transgenic eggplants, exhibiting enhanced resistance to Alternaria solani , can be generated on plant growth regulators (PGRs)- and antibiotic-free MS medium employing the multi-auto-transformation (MAT) vector, pMAT21 - wasabi defensin , wherein isopentenyl transferase ( ipt ) gene is used as a positive selection marker.

Abstract

Use of the selection marker genes conferring antibiotic or herbicide resistance in transgenic plants has been considered a serious problem for environment and the public. Multi-auto-transformation (MAT) vector system has been one of the tools to excise the selection marker gene and produce marker-free transgenic plants. Ipt gene was used as a selection marker gene. Wasabi defensin gene, isolated from Wasabia japonica (a Japanese horseradish which has been a potential source of antimicrobial proteins), was used as a gene of interest. Wasabi defensin gene was cloned from the binary vector, pEKH-WD, to an ipt-type MAT vector, pMAT21, by gateway cloning technology and transferred to Agrobacterium tumefaciens strain EHA105. Infected cotyledon explants of eggplant were cultured on PGRs- and antibiotic-free MS medium. Extreme shooty phenotype/ipt shoots were produced by the explants infected with the pMAT21-wasabi defensin (WD). The same PGRs- and antibiotic-free MS medium was used in subcultures of the ipt shoots. Subsequently, morphologically normal shoots emerged from the Ipt shoots. Molecular analyses of genomic DNA from transgenic plants confirmed the integration of the WD gene and excision of the selection marker (ipt gene). Expression of the WD gene was confirmed by RT-PCR and Northern blot analyses. In vitro whole plant and detached leaf assay of the marker-free transgenic plants exhibited enhanced resistance against Alternaria solani.     

BGVD:An Integrated Database for Bovine Sequencing Variations and Selective Signatures

Next-generation sequencing has yielded a vast amount of cattle genomic data for global characterization of population genetic diversity and identification of genomic regions under natural and artificial selection. However, efficient storage, querying, and visualization of such large datasets remain challenging. Here, we developed a comprehensive database, the Bovine Genome Variation Database (BGVD). It provides six main functionalities: gene search, variation search, genomic signature search, Genome Browser, alignment search tools, and the genome coordinate conversion tool. BGVD contains information on genomic variations comprising ~60.44 M SNPs, ~6.86 M indels, 76,634 CNV regions, and signatures of selective sweeps in 432 samples from modern cattle worldwide. Users can quickly retrieve distribution patterns of these variations for 54 cattle breeds through an interactive source of breed origin map, using a given gene symbol or genomic region for any of the three versions of the bovine reference genomes (ARS-UCD1.2, UMD3.1.1, and Btau 5.0.1). Signals of selection sweep are displayed as Manhattan plots and Genome Browser tracks. To further investigate and visualize the relationships between variants and signatures of selection, the Genome Browser integrates all variations, selection data, and resources, from NCBI, the UCSC Genome Browser, and Animal QTLdb. Collectively, all these features make the BGVD a useful archive for in-depth data mining and analyses of cattle biology and cattle breeding on a global scale. BGVD is publicly available at http://animal.nwsuaf.edu.cn/BosVar.     

Genetic mapping of a putative Thinopyrum intermedium-derived stripe rust resistance gene on wheat chromosome 1B

Key message

Stripe rust resistance transferred from Thinopyrum intermedium into common wheat was controlled by a single dominant gene, which mapped to chromosome 1B near Yr26 and was designated YrL693.

Abstract

Stripe rust caused by Puccinia striiformis f. sp. tritici (Pst) is a highly destructive disease of wheat (Triticum aestivum). Stripe rust resistance was transferred from Thinopyrum intermedium to common wheat, and the resulting introgression line (L693) exhibited all-stage resistance to the widely virulent and predominant Chinese pathotypes CYR32 and CYR33 and to the new virulent pathotype V26. There was no cytological evidence that L693 had alien chromosomal segments from Th. intermedium. Genetic analysis of stripe rust resistance was performed by crossing L693 with the susceptible line L661. F₁, F₂, and F_2:3 populations from reciprocal crosses showed that resistance was controlled by a single dominant gene. A total 479 F_2:3 lines and 781 pairs of genomic simple sequence repeat (SSR) primers were employed to determine the chromosomal location of the resistance gene. The gene was linked to six publicly available and three recently developed wheat genomic SSR markers. The linked markers were localized to wheat chromosome 1B using Chinese Spring nulli-tetrasomic lines, and the resistance gene was localized to chromosome 1B based on SSR and wheat genomic information. A high-density genetic map was also produced. The pedigree, molecular marker data, and resistance response indicated that the stripe rust resistance gene in L693 is a novel gene, which was temporarily designated YrL693. The SSR markers that co-segregate with this gene (Xbarc187-1B, Xbarc187-1B-1, Xgwm18-1B, and Xgwm11-1B) have potential application in marker-assisted breeding of wheat, and YrL693 will be useful for broadening the genetic basis of stripe rust resistance in wheat.     

A machine learning approach for the identification of population-informative markers from high-throughput genotyping data: application to several pig breeds

Single nucleotide polymorphisms (SNPs) able to describe population differences can be used for important applications in livestock, including breed assignment of individual animals, authentication of mono-breed products and parentage verification among several other applications. To identify the most discriminating SNPs among thousands of markers in the available commercial SNP chip tools, several methods have been used. Random forest (RF) is a machine learning technique that has been proposed for this purpose. In this study, we used RF to analyse PorcineSNP60 BeadChip array genotyping data obtained from a total of 2737 pigs of 7 Italian pig breeds (3 cosmopolitan-derived breeds: Italian Large White, Italian Duroc and Italian Landrace, and 4 autochthonous breeds: Apulo-Calabrese, Casertana, Cinta Senese and Nero Siciliano) to identify breed informative and reduced SNP panels using the mean decrease in the Gini Index and the Mean Decrease in Accuracy parameters with stability evaluation. Other reduced informative SNP panels were obtained using Delta, Fixation index and principal component analysis statistics, and their performances were compared with those obtained using the RF-defined panels using the RF classification method and its derived Out Of Bag rates and correct prediction proportions. Therefore, the performances of a total of six reduced panels were evaluated. The correct assignment of the animals to its breed was close to 100% for all tested approaches. Porcine chromosome 8 harboured the largest number of selected SNPs across all panels. Many SNPs were included in genomic regions in which previous studies identified signatures of selection or genes (e.g. ESR1, KITL and LCORL) that could contribute to explain, at least in part, phenotypically or economically relevant traits that might differentiate cosmopolitan and autochthonous pig breeds. Random forest used as preselection statistics highlighted informative SNPs that were not the same as those identified by other methods. This might be due to specific features of this machine learning methodology. It will be interesting to explore if the adaptation of RF methods for the identification of selection signature regions could be able to describe population-specific features that are not captured by other approaches.     

Development of wheat–Aegilops speltoides recombinants and simple PCR-based markers for Sr32 and a new stem rust resistance gene on the 2S#1 chromosome

Key message

Wheat– Aegilops speltoides recombinants carrying stem rust resistance genes Sr32 and SrAes1t effective against Ug99 and PCR markers for marker-assisted selection.

Abstract

Wild relatives of wheat are important resources for new rust resistance genes but underutilized because the valuable resistances are often linked to negative traits that prevent deployment of these genes in commercial wheats. Here, we report ph1b-induced recombinants with reduced alien chromatin derived from E.R. Sears’ wheat–Aegilops speltoides 2D-2S#1 translocation line C82.2, which carries the widely effective stem rust resistance gene Sr32. Infection type assessments of the recombinants showed that the original translocation in fact carries two stem rust resistance genes, Sr32 on the short arm and a previously undescribed gene SrAes1t on the long arm of chromosome 2S#1. Recombinants with substantially shortened alien chromatin were produced for both genes, which confer resistance to stem rust races in the TTKSK (Ug99) lineage and representative races of all Australian stem rust lineages. Selected recombinants were back crossed into adapted Australian cultivars and PCR markers were developed to facilitate the incorporation of these genes into future wheat varieties. Our recombinants and those from several other labs now show that Sr32, Sr39, and SrAes7t on the short arm and Sr47 and SrAes1t on the long arm of 2S#1 form two linkage groups and at present no rust races are described that can distinguish these resistance specificities.     

The Strong Selective Sweep Candidate Gene ADRA2C Does Not Explain Domestication Related Changes In The Stress Response Of Chickens

Analysis of selective sweeps to pinpoint causative genomic regions involved in chicken domestication has revealed a strong selective sweep on chromosome 4 in layer chickens. The autoregulatory α-adrenergic receptor 2C (ADRA2C) gene is the closest to the selective sweep and was proposed as an important gene in the domestication of layer chickens. The ADRA2C promoter region was also hypermethylated in comparison to the non-selected ancestor of all domesticated chicken breeds, the Red Junglefowl, further supporting its relevance. In mice the receptor is involved in the fight-or-flight response as it modulates epinephrine release from the adrenals. To investigate the involvement of ADRA2C in chicken domestication, we measured gene expression in the adrenals and radiolabeled receptor ligand in three brain regions comparing the domestic White Leghorn strain with the wild ancestor Red Junglefowl. In adrenals ADRA2C was twofold greater expressed than the related receptor gene ADRA2A, indicating that ADRA2C is the predominant modulator of epinephrine release but no strain differences were measured. In hypothalamus and amygdala, regions associated with the stress response, and in striatum, receptor binding pIC50 values ranged between 8.1–8.4, and the level was not influenced by the genotyped allele. Because chicken strains differ in morphology, physiology and behavior, differences attributed to a single gene may be lost in the noise caused by the heterogeneous genetic background. Therefore an F₁₀ advanced intercross strain between White Leghorn and Red Junglefowl was used to investigate effects of ADRA2C alleles on fear related behaviors and fecundity. We did not find compelling genotype effects in open field, tonic immobility, aerial predator, associative learning or fecundity. Therefore we conclude that ADRA2C is probably not involved in the domestication of the stress response in chicken, and the strong selective sweep is probably caused by selection of some unknown genetic element in the vicinity of the gene.     

Detection of selection signatures in dairy and beef cattle using high-density genomic information

Background

Artificial selection for economically important traits in cattle is expected to have left distinctive selection signatures on the genome. Access to high-density genotypes facilitates the accurate identification of genomic regions that have undergone positive selection. These findings help to better elucidate the mechanisms of selection and to identify candidate genes of interest to breeding programs.

Results

Information on 705 243 autosomal single nucleotide polymorphisms (SNPs) in 3122 dairy and beef male animals from seven cattle breeds (Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental) were used to detect selection signatures by applying two complementary methods, integrated haplotype score (iHS) and global fixation index (F_ST). To control for false positive results, we used false discovery rate (FDR) adjustment to calculate adjusted iHS within each breed and the genome-wide significance level was about 0.003. Using the iHS method, 83, 92, 91, 101, 85, 101 and 86 significant genomic regions were detected for Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental cattle, respectively. None of these regions was common to all seven breeds. Using the F_ST approach, 704 individual SNPs were detected across breeds. Annotation of the regions of the genome that showed selection signatures revealed several interesting candidate genes i.e. DGAT1, ABCG2, MSTN, CAPN3, FABP3, CHCHD7, PLAG1, JAZF1, PRKG2, ACTC1, TBC1D1, GHR, BMP2, TSG1, LYN, KIT and MC1R that play a role in milk production, reproduction, body size, muscle formation or coat color. Fifty-seven common candidate genes were found by both the iHS and global F_ST methods across the seven breeds. Moreover, many novel genomic regions and genes were detected within the regions that showed selection signatures; for some candidate genes, signatures of positive selection exist in the human genome. Multilevel bioinformatic analyses of the detected candidate genes suggested that the PPAR pathway may have been subjected to positive selection.

Conclusions

This study provides a high-resolution bovine genomic map of positive selection signatures that are either specific to one breed or common to a subset of the seven breeds analyzed. Our results will contribute to the detection of functional candidate genes that have undergone positive selection in future studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0127-3) contains supplementary material, which is available to authorized users.