Genome-wide association study of birth weight in sheep

Birth weight is the earliest available growth trait with considerable impacts on lamb survivability and growth performance traits. This study was conducted to perform a genome-wide association study of birth weight in a meat-type sheep. A total of 132 Lori-Bakhtiari sheep were selected based on estimated of breeding values (EBVs) for BW analyses. The selected animals were genotyped using Illumina Ovine SNP50 Bead Chip. After quality control, a total of 41 323 single-nucleotide polymorphisms (SNPs) and 130 sheep were used for subsequent analyses. Plink 1.90 beta software was used for the analyses. Seven SNPs on chromosomes 1, 16, 19 and 22 were detected based on genome-wide unadjusted P-values (P <10⁻⁶), which jointly accounted for 1.2% of total genetic variation. However, based on Bonferroni-adjusted P-values, only three SNPs on chromosome 1 had significant associations with EBVs for birth weight (P <0.05), which jointly explained 0.8% of total genetic variation. A total of seven genes were found in 50 kb intervals from the three significant SNPs on chromosome 1, but only three genes, including RAB6B (a member of RAS oncogene family), Tf serotransferrin and GIGYF2 (a GRB10 interacting GYF protein 2), could be considered as candidate genes for birth weight in future studies. The results of this study may facilitate potential use of the genes involving in growth and production traits for genetic improvement of productivity in sheep.     

Comparative Analysis of Disease-Linked Single Nucleotide Polymorphic Markers from Brassica rapa for Their Applicability to Brassica oleracea

Numerous studies using single nucleotide polymorphisms (SNPs) have been conducted in humans, and other animals, and in major crops, including rice, soybean, and Chinese cabbage. However, the number of SNP studies in cabbage is limited. In this present study, we evaluated whether 7,645 SNPs previously identified as molecular markers linked to disease resistance in the Brassica rapa genome could be applied to B. oleracea. In a BLAST analysis using the SNP sequences of B. rapa and B. oleracea genomic sequence data registered in the NCBI database, 256 genes for which SNPs had been identified in B. rapa were found in B. oleracea. These genes were classified into three functional groups: molecular function (64 genes), biological process (96 genes), and cellular component (96 genes). A total of 693 SNP markers, including 145 SNP markers [BRH—developed from the B. rapa genome for high-resolution melt (HRM) analysis], 425 SNP markers (BRP—based on the B. rapa genome that could be applied to B. oleracea), and 123 new SNP markers (BRS—derived from BRP and designed for HRM analysis), were investigated for their ability to amplify sequences from cabbage genomic DNA. In total, 425 of the SNP markers (BRP-based on B. rapa genome), selected from 7,645 SNPs, were successfully applied to B. oleracea. Using PCR, 108 of 145 BRH (74.5%), 415 of 425 BRP (97.6%), and 118 of 123 BRS (95.9%) showed amplification, suggesting that it is possible to apply SNP markers developed based on the B. rapa genome to B. oleracea. These results provide valuable information that can be utilized in cabbage genetics and breeding programs using molecular markers derived from other Brassica species.     

Imprints of Natural Selection Along Environmental Gradients in Phenology-Related Genes of Quercus petraea

We explored single nucleotide polymorphism (SNP) variation in candidate genes for bud burst from Quercus petraea populations sampled along gradients of latitude and altitude in Western Europe. SNP diversity was monitored for 106 candidate genes, in 758 individuals from 32 natural populations. We investigated whether SNP variation reflected the clinal pattern of bud burst observed in common garden experiments. We used different methods to detect imprints of natural selection (F_ST outlier, clinal variation at allelic frequencies, association tests) and compared the results obtained for the two gradients. F_ST outlier SNPs were found in 15 genes, 5 of which were common to both gradients. The type of selection differed between the two gradients (directional or balancing) for 3 of these 5. Clinal variations were observed for six SNPs, and one cline was conserved across both gradients. Association tests between the phenotypic or breeding values of trees and SNP genotypes identified 14 significant associations, involving 12 genes. The results of outlier detection on the basis of population differentiation or clinal variation were not very consistent with the results of association tests. The discrepancies between these approaches may reflect the different hierarchical levels of selection considered (inter- and intrapopulation selection). Finally, we obtained evidence for convergent selection (similar for gradients) and clinal variation for a few genes, suggesting that comparisons between parallel gradients could be used to screen for major candidate genes responding to natural selection in trees.     

A catalogue of validated single nucleotide polymorphisms in bovine orthologs of mammalian imprinted genes and associations with beef production traits

Genetic (or ‘genomic’) imprinting, a feature of approximately 100 mammalian genes, results in monoallelic expression from one of the two parentally inherited chromosomes. To date, most studies have been directed on imprinted genes in murine or human models; however, there is burgeoning interest in the effects of imprinted genes in domestic livestock species. In particular, attention has focused on imprinted genes that influence foetal growth and development and that are associated with several economically important production traits in cattle, sheep and pigs. We have re-sequenced regions in 20 candidate bovine imprinted genes in order to validate single nucleotide polymorphisms (SNPs) that may influence important production traits in cattle. Putative SNPs detected via re-sequencing were subsequently re-formatted for high-throughput SNP genotyping in 185 cattle samples comprising 138 performance-tested European Bos taurus (all Limousin bulls), 29 African B. taurus and 18 Indian B. indicus samples. Analysis of the resulting genotypic data identified 117 validated SNPs. Preliminary genotype–phenotype association analyses using 83 SNPs that were polymorphic in the Limousin samples with minor allele frequencies ⩾0.05 revealed significant associations between two candidate bovine imprinted genes and a range of important beef production traits: average daily gain, average feed intake, live weight, feed conversion ratio, residual feed intake and residual gain. These genes were the Ras protein-specific guanine nucleotide releasing factor gene (RASGRF1) and the zinc finger, imprinted 2 gene (ZIM2). Despite the relatively small sample size used in these analyses, the observed associations with production traits are supported by the purported biological function of the RASGRF1 and ZIM2 gene products. These results support the hypothesis that imprinted genes contribute significantly to important complex production traits in cattle. Furthermore, these SNPs may be usefully incorporated into future marker-assisted and genomic selection breeding schemes.     

Genetic Dissection of Drought and Heat Tolerance in Chickpea through Genome-Wide and Candidate Gene-Based Association Mapping Approaches

To understand the genetic basis of tolerance to drought and heat stresses in chickpea, a comprehensive association mapping approach has been undertaken. Phenotypic data were generated on the reference set (300 accessions, including 211 mini-core collection accessions) for drought tolerance related root traits, heat tolerance, yield and yield component traits from 1–7 seasons and 1–3 locations in India (Patancheru, Kanpur, Bangalore) and three locations in Africa (Nairobi, Egerton in Kenya and Debre Zeit in Ethiopia). Diversity Array Technology (DArT) markers equally distributed across chickpea genome were used to determine population structure and three sub-populations were identified using admixture model in STRUCTURE. The pairwise linkage disequilibrium (LD) estimated using the squared-allele frequency correlations (r²; when r²<0.20) was found to decay rapidly with the genetic distance of 5 cM. For establishing marker-trait associations (MTAs), both genome-wide and candidate gene-sequencing based association mapping approaches were conducted using 1,872 markers (1,072 DArTs, 651 single nucleotide polymorphisms [SNPs], 113 gene-based SNPs and 36 simple sequence repeats [SSRs]) and phenotyping data mentioned above employing mixed linear model (MLM) analysis with optimum compression with P3D method and kinship matrix. As a result, 312 significant MTAs were identified and a maximum number of MTAs (70) was identified for 100-seed weight. A total of 18 SNPs from 5 genes (ERECTA, 11 SNPs; ASR, 4 SNPs; DREB, 1 SNP; CAP2 promoter, 1 SNP and AMDH, 1SNP) were significantly associated with different traits. This study provides significant MTAs for drought and heat tolerance in chickpea that can be used, after validation, in molecular breeding for developing superior varieties with enhanced drought and heat tolerance.     

Identification of loci controlling forage yield and nutritive value in diploid alfalfa using GBS-GWAS

Key message

We attempted to identify genomic regions controlling forage yield and nutritive value in alfalfa. Several candidate genes and associated genetic markers were identified that could potentially be useful for alfalfa breeding to more efficiently develop improved cultivars.

Abstract

Alfalfa is one of the most widely cultivated forage legumes worldwide and improving alfalfa forage yield and nutritive value is a major global breeding goal. Genotyping-by-sequencing (GBS) provides cost-effective molecular marker genotyping for genome-wide association studies (GWAS). Using more than 15,000 genome-wide single nucleotide polymorphisms (SNP) identified from GBS, we conducted a GWAS to investigate forage yield and nutritive value-related traits. We have detected a number of associations for all the traits evaluated and a number of associations detected were located on the Medicago truncatula genome. The SNP in a coding region of a cell wall biosynthesis gene was associated with several cell wall-related traits, and we suggest that it may be the causative polymorphism. Two other SNPs residing in meristematic development and early growth genes were found to associate with the total biomass yield. None of the SNPs associated with regrowth after harvest or with spring regrowth were mapped to the M. truncatula genome, possibly reflecting the fact that M. truncatula is an annual species related to alfalfa that typically has limited ability to regrow. The alleles we identify with the major impact on forage yield and nutritive value can be rapidly incorporated into our breeding program.

     

Association of IGF1 and KDM5A polymorphisms with performance, fatness and carcass traits in chickens

Two functional and positional candidate genes were selected in a region of chicken chromosome 1 (GGA1), based on their biological roles, and also where several quantitative trait loci (QTL) have been mapped and associated with performance, fatness and carcass traits in chickens. The insulin-like growth factor 1 (IGF1) gene has been associated with several physiological functions related to growth. The lysine (K)-specific demethylase 5A (KDM5A) gene participates in the epigenetic regulation of genes involved with the cell cycle. Our objective was to find associations of selected single-nucleotide polymorphisms (SNPs) in these genes with performance, fatness and carcass traits in 165 F₂ chickens from a resource population. In the IGF1 gene, 17 SNPs were detected, and in the KDM5A gene, nine SNPs were detected. IGF1 SNP c.47673G?>?A was associated with body weight and haematocrit percentage, and also with feed intake and percentages of abdominal fat and gizzard genotype × sex interactions. KDM5A SNP c.34208C?>?T genotype × sex interaction affected body weight, feed intake, percentages of abdominal fat (p?=?0.0001), carcass, gizzard and haematocrit. A strong association of the diplotype × sex interaction (p?<?0.0001) with abdominal fat was observed, and also associations with body weight, feed intake, percentages of carcass, drums and thighs, gizzard and haematocrit. Our findings suggest that the KDM5A gene might play an important role in the abdominal fat deposition in chickens. The IGF1 and KDM5A genes are strong candidates to explain the QTL mapped in this region of GGA1.     

Analyzing Large-Scale Samples Confirms the Association Between the ABCA7 rs3764650 Polymorphism and Alzheimer’s Disease Susceptibility

Large-scale genome-wide association studies (GWAS) have revealed that the ABCA7 rs3764650 polymorphism (or its proxies, namely rs115550680, rs3752246, and rs4147929) is associated with Alzheimer’s disease (AD) susceptibility in individuals of Caucasian ancestry. The following studies have investigated this finding in Chinese (N?=?633 and N?=?1,224), Japanese (N?=?1,735), Korean (N?=?844), African American (N?=?5,896), and Canadian (N?=?1,104) populations. However, these studies reported a weak or negligible association. We hypothesized that these negative results may have been caused by either relatively small sample sizes compared with those used for the previous GWAS in individuals of Caucasian ancestry or the genetic heterogeneity of the rs3764650 polymorphism (or its proxies) in different populations. Here, we reevaluated the association between rs3764650 and AD using large-scale samples from 18 previous studies (N?=?79,381—30,590 cases and 48,791 controls) by searching PubMed, AlzGene, and Google Scholar databases. Using allele, dominant, recessive, and additive models, we did not identify significant heterogeneity among the 18 studies. We observed a significant association between rs3764650 and AD using the allele (P?=?1.76E???26, odds ratio (OR)?=?1.21, 95 % confidence interval (CI) 1.17–1.26), dominant (P?=?4.00E???04, OR?=?1.17, 95 % CI 1.07–1.28), recessive (P?=?3.00E???03, OR?=?1.43, 95 % CI 1.13–1.81), and additive models (P?=?3.00E???03, OR?=?1.49, 95 % CI 1.16–1.91). Collectively, our analysis further supports previous findings that the ABCA7 rs3764650 polymorphism is associated with AD susceptibility. We believe that our findings will be very useful for future genetic studies on AD.     

Development of single nucleotide polymorphism markers in the large and complex rubber tree genome using next-generation sequence data

The development of single nucleotide polymorphism (SNP) markers provides the opportunity to improve many areas of plant breeding and population genetics. Unfortunately, for species such as the rubber tree (Hevea brasiliensis), the use of next-generation sequencing for genomic SNP discovery is very difficult because of the large genome size and the abundance of repeated sequences. Access to a set of validated SNP markers is a significant advantage for rubber researchers who wish to apply SNPs in scientific research. Here, we performed genomic sequencing of H. brasiliensis and generated 10,993,648 short reads, which were assembled into 10,071 contigs (N50 = 3078) by a de novo assembly strategy. A total of 2446 contigs presented no hits in the current H. brasiliensis genome assembly and may therefore be considered novel genomic sequences of rubber tree. A total of 143 putative polymorphic positions were selected, gene annotations were available for 58.7 % of the markers, and all of the sequences could be anchored to the released H. brasiliensis genome. These SNPs were validated in eight genotypes of H. brasiliensis and 15 F1 plants from a mapping population, resulting in 30 (20.9 %) positions correctly classified. The analysis revealed key candidate genes responsible for defence mechanisms and provided markers for further genetic improvement of Hevea in breeding programmes.     

Identification of favorable SNP alleles and candidate genes responsible for inflorescence-related traits via GWAS in chrysanthemum

Key message

81 SNPs were identified for three inflorescence-related traits, in which 15 were highly favorable. Two dCAPS markers were developed for future MAS breeding, and six candidate genes were predicted.

Abstract

Chrysanthemum is a leading ornamental species worldwide and demonstrates a wealth of morphological variation. Knowledge about the genetic basis of its phenotypic variation for key horticultural traits can contribute to its effective management and genetic improvement. In this study, we conducted a genome-wide association study (GWAS) based on two years of phenotype data and a set of 92,617 single nucleotide polymorphisms (SNPs) using a panel of 107 diverse cut chrysanthemums to dissect the genetic control of three inflorescence-related traits. A total of 81 SNPs were significantly associated with the three inflorescence-related traits (capitulum diameter, number of ray florets and flowering time) in at least one environment, with an individual allele explaining 22.72–38.67% of the phenotypic variation. Fifteen highly favorable alleles were identified for the three target traits by computing the phenotypic effect values for the stable associations detected in 2 year-long trials at each locus. Dosage pyramiding effects of the highly favorable SNP alleles and significant linear correlations between highly favorable allele numbers and corresponding phenotypic performance were observed. Two highly favorable SNP alleles correlating to flowering time and capitulum diameter were converted to derived cleaved amplified polymorphic sequence (dCAPS) markers to facilitate future breeding. Finally, six putative candidate genes were identified that contribute to flowering time and capitulum diameter. These results serve as a foundation for analyzing the genetic mechanisms underlying important horticultural traits and provide valuable insights into molecular marker-assisted selection (MAS) in chrysanthemum breeding programs.

     

Genetic variation of wood chemical traits and association with underlying genes in Eucalyptus urophylla

We developed a quantitative and association genetic study with Eucalyptus urophylla using a progeny trial. Based on a sample of 831 trees distributed in 84 half-sib families whose wood was phenotyped by near-infrared spectroscopy, the results showed that traits related to lignin, cellulose, and wood extractives presented significant additive genetic variability with moderate to high narrow sense heritability (h ²?=?0.28 to 0.93). Genetic correlations varied with high standard error and showed low to moderate values. Using three cellulose synthase genes (EuCesA1, EuCesA2, and EuCesA3) and three candidate genes involved in the lignin pathway (EuC4H1, EuC4H2, and EuCAD2), an association study was performed for each of the gene action models (co-dominant, recessive, and dominant) using two methods. Firstly, single-marker association tests were done and 539 tests (49 single nucleotide polymorphisms (SNPs)?×?11 traits) were analyzed. After Bonferroni correction with a significance level of P?=?0.00102, only four SNPs presented significant association with syringyl and syringyl-to-guaiacyl ratio with an adjusted coefficient of determination varying between 2.6 and 4.4 %. Secondly, a model selection method, the backward approach, was implemented. Similar SNPs were detected by both the backward selection and the individual marker approaches. However, the latter detected new associations with other traits, genes, and SNPs and improved the quality of the model as shown by the BIC criteria and the higher adjusted determination coefficient (1.5 to 8.3 %). Our results reveal that cellulose genes can be associated with lignin traits (syringyl-to-guaiacyl ratio) and stress the possible pleiotropic effect of some genes.     

Genome-Wide Association of Stem Water Soluble Carbohydrates in Bread Wheat

Water soluble carbohydrates (WSC) in stems play an important role in buffering grain yield in wheat against biotic and abiotic stresses; however, knowledge of genes controlling WSC is very limited. We conducted a genome-wide association study (GWAS) using a high-density 90K SNP array to better understand the genetic basis underlying WSC, and to explore marker-based breeding approaches. WSC was evaluated in an association panel comprising 166 Chinese bread wheat cultivars planted in four environments. Fifty two marker-trait associations (MTAs) distributed across 23 loci were identified for phenotypic best linear unbiased estimates (BLUEs), and 11 MTAs were identified in two or more environments. Liner regression showed a clear dependence of WSC BLUE scores on numbers of favorable (increasing WSC content) and unfavorable alleles (decreasing WSC), indicating that genotypes with higher numbers of favorable or lower numbers of unfavorable alleles had higher WSC content. In silico analysis of flanking sequences of trait-associated SNPs revealed eight candidate genes related to WSC content grouped into two categories based on the type of encoding proteins, namely, defense response proteins and proteins triggered by environmental stresses. The identified SNPs and candidate genes related to WSC provide opportunities for breeding higher WSC wheat cultivars.     

Association analysis for udder health based on SNP-panel and sequence data in Danish Holsteins

Background

The sensitivity of genome-wide association studies for the detection of quantitative trait loci (QTL) depends on the density of markers examined and the statistical models used. This study compares the performance of three marker densities to refine six previously detected QTL regions for mastitis traits: 54 k markers of a medium-density SNP (single nucleotide polymorphism) chip (MD), imputed 777 k markers of a high-density SNP chip (HD), and imputed whole-genome sequencing data (SEQ). Each dataset contained data for 4496 Danish Holstein cattle. Comparisons were performed using a linear mixed model (LM) and a Bayesian variable selection model (BVS).

Results

After quality control, 587, 7825, and 78 856 SNPs in the six targeted regions remained for MD, HD, and SEQ data, respectively. In general, the association patterns between SNPs and traits were similar for the three marker densities when tested using the same statistical model. With the LM model, 120 (MD), 967 (HD), and 7209 (SEQ) SNPs were significantly associated with mastitis, whereas with the BVS model, 43 (MD), 131 (HD), and 1052 (SEQ) significant SNPs (Bayes factor > 3.2) were observed. A total of 26 (MD), 75 (HD), and 465 (SEQ) significant SNPs were identified by both models. In addition, one, 16, and 33 QTL peaks for MD, HD, and SEQ data were detected according to the QTL intensity profile of SNP bins by post-analysis of the BVS model.

Conclusions

The power to detect significant associations increased with increasing marker density. The BVS model resulted in clearer boundaries between linked QTL than the LM model. Using SEQ data, the six targeted regions were refined to 33 candidate QTL regions for udder health. The comparison between these candidate QTL regions and known genes suggested that NPFFR2, SLC4A4, DCK, LIFR, and EDN3 may be considered as candidate genes for mastitis susceptibility.

Electronic supplementary material

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

Genome-Wide Association Studies Identify the Loci for 5 Exterior Traits in a Large White × Minzhu Pig Population

As one of the main breeding selection criteria, external appearance has special economic importance in the hog industry. In this study, an Illumina Porcine SNP60 BeadChip was used to conduct a genome-wide association study (GWAS) in 605 pigs of the F₂ generation derived from a Large White × Minzhu intercross. Traits under study were abdominal circumference (AC), body height (BH), body length (BL), cannon bone circumference (CBC), chest depth (CD), chest width (CW), rump circumference (RC), rump width (RW), scapula width (SW), and waist width (WW). A total of 138 SNPs (the most significant being MARC0033464) on chromosome 7 were found to be associated with BH, BL, CBC, and RC (P-value  = 4.15E-6). One SNP on chromosome 1 was found to be associated with CD at genome-wide significance levels. The percentage phenotypic variance of these significant SNPs ranged from 0.1–25.48%. Moreover, a conditional analysis revealed that the significant SNPs were derived from a single quantitative trait locus (QTL) and indicated additional chromosome-wide significant association for 25 SNPs on SSC4 (BL, CBC) and 9 SNPs on SSC7 (RC). Linkage analysis revealed two complete linkage disequilibrium haplotype blocks that contained seven and four SNPs, respectively. In block 1, the most significant SNP, MARC0033464, was present. Annotations from pig reference genome suggested six genes (GRM4, HMGA1, NUDT3, RPS10, SPDEF and PACSIN1) in block 1 (495 kb), and one gene (SCUBE3) in block 3 (124 kb). Functional analysis indicated that HMGA1 and SCUBE3 genes are the potential genes controlling BH, BL, and RC in pigs, with an application in breeding programs. We screened several candidate intervals and genes based on SNP location and gene function, and predicted their function using bioinformatics analyses.     

SNP-based pool genotyping and haplotype analysis accelerate fine-mapping of the wheat genomic region containing stripe rust resistance gene <Emphasis Type="Italic">Yr26</Emphasis>

Key message

NGS-assisted super pooling emerging as powerful tool to accelerate gene mapping and haplotype association analysis within target region uncovering specific linkage SNPs or alleles for marker-assisted gene pyramiding.

Abstract

Conventional gene mapping methods to identify genes associated with important agronomic traits require significant amounts of financial support and time. Here, a single nucleotide polymorphism (SNP)-based mapping approach, RNA-Seq and SNP array assisted super pooling analysis, was used for rapid mining of a candidate genomic region for stripe rust resistance gene Yr26 that has been widely used in wheat breeding programs in China. Large DNA and RNA super-pools were genotyped by Wheat SNP Array and sequenced by Illumina HiSeq, respectively. Hundreds of thousands of SNPs were identified and then filtered by multiple filtering criteria. Among selected SNPs, over 900 were found within an overlapping interval of less than 30 Mb as the Yr26 candidate genomic region in the centromeric region of chromosome arm 1BL. The 235 chromosome-specific SNPs were converted into KASP assays to validate the Yr26 interval in different genetic populations. Using a high-resolution mapping population (>?30,000 gametes), we confined Yr26 to a 0.003-cM interval. The Yr26 target region was anchored to the common wheat IWGSC RefSeq v1.0 and wild emmer WEWSeq v.1.0 sequences, from which 488 and 454 kb fragments were obtained. Several candidate genes were identified in the target genomic region, but there was no typical resistance gene in either genome region. Haplotype analysis identified specific SNPs linked to Yr26 and developed robust and breeder-friendly KASP markers. This integration strategy can be applied to accelerate generating many markers closely linked to target genes/QTL for a trait of interest in wheat and other polyploid species.

     

Targeted re-sequencing and genome-wide association analysis for wood property traits in breeding population of Eucalyptus tereticornis × E. grandis

Globally, Eucalyptus plantations occupy 22 million ha area and is one of the preferred hardwood species due to their short rotation, rapid growth, adaptability and wood properties. In this study, we present results of GWAS in parents and 100 hybrids of Eucalyptus tereticornis × E. grandis using 762 genes presumably involved in wood formation. Comparative analysis between parents predicted 32,202 polymorphic SNPs with high average read depth of 269-562× per individual per nucleotide. Seventeen wood related traits were phenotyped across three diverse environments and GWAS was conducted using 13,610 SNPs. A total of 45 SNP-trait associations were predicted across two locations. Seven large effect markers were identified which explained more than 80% of phenotypic variation for fibre area. This study has provided an array of candidate genes which may govern fibre morphology in this genus and has predicted potential SNPs which can guide future breeding programs in tropical Eucalyptus.