Genome-wide association analyses reveal significant loci and strong candidate genes for growth and fatness traits in two pig populations

Background

Recently, genome-wide association studies (GWAS) have been reported on various pig traits. We performed a GWAS to analyze 22 traits related to growth and fatness on two pig populations: a White Duroc × Erhualian F₂ intercross population and a Chinese Sutai half-sib population.

Results

We identified 14 and 39 loci that displayed significant associations with growth and fatness traits at the genome-wide level and chromosome-wide level, respectively. The strongest association was between a 750 kb region on SSC7 (SSC for Sus scrofa) and backfat thickness at the first rib. This region had pleiotropic effects on both fatness and growth traits in F₂ animals and contained a promising candidate gene HMGA1 (high mobility group AT-hook 1). Unexpectedly, population genetic analysis revealed that the allele at this locus that reduces fatness and increases growth is derived from Chinese indigenous pigs and segregates in multiple Chinese breeds. The second strongest association was between the region around 82.85 Mb on SSC4 and average backfat thickness. PLAG1 (pleiomorphic adenoma gene 1), a gene under strong selection in European domestic pigs, is proximal to the top SNP and stands out as a strong candidate gene. On SSC2, a locus that significantly affects fatness traits mapped to the region around the IGF2 (insulin-like growth factor 2) gene but its non-imprinting inheritance excluded IGF2 as a candidate gene. A significant locus was also detected within a recombination cold spot that spans more than 30 Mb on SSCX, which hampered the identification of plausible candidate genes. Notably, no genome-wide significant locus was shared by the two experimental populations; different loci were observed that had both constant and time-specific effects on growth traits at different stages, which illustrates the complex genetic architecture of these traits.

Conclusions

We confirm several previously reported QTL and provide a list of novel loci for porcine growth and fatness traits in two experimental populations with Chinese Taihu and Western pigs as common founders. We showed that distinct loci exist for these traits in the two populations and identified HMGA1 and PLAG1 as strong candidate genes on SSC7 and SSC4, respectively.

Electronic supplementary material

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

Mitochondrial and nuclear DNA analyses reveal population differentiation in Brazilian Creole sheep

Using ND5 sequences from mtDNA and 10 nuclear markers, we investigated the genetic differentiation of two South American Creole sheep phenotypes that historically have been bred in different biomes in southern Brazil. In total, 18 unique mtDNA haplotypes were detected, none of which was shared between the two phenotypes. Bayesian analysis also indicated two different groups (k = 2). Thus, these varieties are supported as being genotypically distinct. This situation could have resulted either from geographical isolation, associated with differences in the cultural habits of sheep farmers and in the way that flocks were managed, or more likely, from the introduction of different stocks four centuries ago.     

Genome-wide association study to identify genetic loci associated with gastrointestinal nematode resistance in Katahdin sheep

Resistance to gastrointestinal nematodes has previously been shown to be a moderately heritable trait in some breeds of sheep, but the mechanisms of resistance are not well understood. Selection for resistance currently relies upon faecal egg counts (FEC), blood packed cell volumes and FAMACHA visual indicator scores of anaemia. Identifying genomic markers associated with disease resistance would potentially improve the selection process and provide a more reliable means of classifying and understanding the biology behind resistant and susceptible sheep. A GWAS was conducted to identify possible genetic loci associated with resistance to Haemonchus contortus in Katahdin sheep. Forty animals were selected from the top and bottom 10% of estimated breeding values for FEC from a total pool of 641 sires and ram lambs. Samples were genotyped using Applied Biosystems™ Axiom™ Ovine Genotyping Array (50K) consisting of 51 572 SNPs. Following quality control, 46 268 SNPs were included in subsequent analyses. Analyses were conducted using a linear regression model in plink v1.90 and a single-locus mixed model in snp and variation suite . Genome-wide significance was determined by a Bonferroni correction for multiple testing. Using linear regression, loci on chromosomes 2, 3, 16, 23 and 24 were significantly associated at the genome level with FEC estimated breeding values, and we identified a region on chromosome 2 that was significant using both statistical analyses. We suggest a potential role for the gene DIS3L2 for gastrointestinal nematode resistance in Katahdin sheep, although further research is needed to validate these findings.     

Selection signature analysis reveals genes associated with tail type in Chinese indigenous sheep

Fat‐tailed sheep have commercial value because consumers prefer high‐protein and low‐fat food and producers care about feed conversion rate. However, fat‐tailed sheep still have some scientific significance, as the fat tail is commonly regarded as a characteristic of environmental adaptability. Finding the candidate genes associated with fat tail formation is essential for breeding and conservation. To identify these candidate genes, we applied FST and hapFLK approaches in fat‐ and thin‐tailed sheep with available 50K SNP genotype data. These two methods found 6.24 Mb of overlapped regions and 43 genes that may associated with fat tail development. Gene annotation showed that HOXA11, BMP2, PPP1CC, SP3, SP9, WDR92, PROKR1 and ETAA1 may play important roles in fat tail formation. These findings provide insight into tail fat development and a guide for molecular breeding and conservation.     

Genome-wide analysis reveals population structure and selection in Chinese indigenous sheep breeds

Background

Traditionally, Chinese indigenous sheep were classified geographically and morphologically into three groups: Mongolian, Kazakh and Tibetan. Herein, we aimed to evaluate the population structure and genome selection among 140 individuals from ten representative Chinese indigenous sheep breeds: Ujimqin, Hu, Tong, Large-Tailed Han and Lop breed (Mongolian group); Duolang and Kazakh (Kazakh group); and Diqing, Plateau-type Tibetan, and Valley-type Tibetan breed (Tibetan group).

Results

We analyzed the population using principal component analysis (PCA), STRUCTURE and a Neighbor-Joining (NJ)-tree. In PCA plot, the Tibetan and Mongolian groups were clustered as expected; however, Duolang and Kazakh (Kazakh group) were segregated. STRUCTURE analyses suggested two subpopulations: one from North China (Kazakh and Mongolian groups) and the other from the Southwest (Tibetan group). In the NJ-tree, the Tibetan group formed an independent branch and the Kazakh and Mongolian groups were mixed. We then used the d_i statistic approach to reveal selection in Chinese indigenous sheep breeds. Among the 599 genome sequence windows analyzed, sixteen (2.7%) exhibited signatures of selection in four or more breeds. We detected three strong selection windows involving three functional genes: RXFP2, PPP1CC and PDGFD. PDGFD, one of the four subfamilies of PDGF, which promotes proliferation and inhibits differentiation of preadipocytes, was significantly selected in fat type breeds by the Rsb (across pairs of populations) approach. Two consecutive selection regions in Duolang sheep were obviously different to other breeds. One region was in OAR2 including three genes (NPR2, SPAG8 and HINT2) the influence growth traits. The other region was in OAR 6 including four genes (PKD2, SPP1, MEPE, and IBSP) associated with a milk production quantitative trait locus. We also identified known candidate genes such as BMPR1B, MSRB3, and three genes (KIT, MC1R, and FRY) that influence lambing percentage, ear size and coat phenotypes, respectively.

Conclusions

Based on the results presented here, we propose that Chinese native sheep can be divided into two genetic groups: the thin type (Tibetan group), and the fat type (Mongolian and Kazakh group). We also identified important genes that drive valuable phenotypes in Chinese indigenous sheep, especially PDGFD, which may influence fat deposition in fat type sheep.

Electronic supplementary material

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

Genome-wide analyses of light-regulated genes in Aspergillus nidulans reveal a complex interplay between different photoreceptors and novel photoreceptor functions

Fungi sense light of different wavelengths using blue-, green-, and red-light photoreceptors. Blue light sensing requires the “white-collar” proteins with flavin as chromophore, and red light is sensed through phytochrome. Here we analyzed genome-wide gene expression changes caused by short-term, low-light intensity illumination with blue-, red- or far-red light in Aspergillus nidulans and found that more than 1100 genes were differentially regulated. The largest number of up- and downregulated genes depended on the phytochrome FphA and the attached HOG pathway. FphA and the white-collar orthologue LreA fulfill activating but also repressing functions under all light conditions and both appear to have roles in the dark. Additionally, we found about 100 genes, which are red-light induced in the absence of phytochrome, suggesting alternative red-light sensing systems. We also found blue-light induced genes in the absence of the blue-light receptor LreA. We present evidence that cryptochrome may be part of this regulatory cue, but that phytochrome is essential for the response. In addition to in vivo data showing that FphA is involved in blue-light sensing, we performed spectroscopy of purified phytochrome and show that it responds indeed to blue light.     

Combining gene expression and genetic analyses to identify candidate genes involved in cold responses in pea

Cold stress affects plant growth and development. In order to better understand the responses to cold (chilling or freezing tolerance), we used two contrasted pea lines. Following a chilling period, the Champagne line becomes tolerant to frost whereas the Terese line remains sensitive.     

Natural variation of salinity response,population structure and candidate genes associated with salinity tolerance in perennial ryegrass accessions

Natural variation in salinity response, effects of population structure on growth and physiological traits and gene–trait association were examined in 56 global collections of diverse perennial ryegrass (Lolium perenne L.) accessions. Three population structure groups were identified with 66 simple sequence repeat markers, which on average accounted for 9 and 11% of phenotypic variation for the control and salinity treatment at 300 mm NaCl. Group 1 (10 accessions) had greater plant height, leaf dry weight and water content, chlorophyll index, K⁺ concentration and K⁺/Na⁺ than group 2 (39 accessions) and group 3 (7 accessions) under salinity stress, while group 3 had higher Na⁺ than groups 1 and 2. Eighty‐seven single nucleotide polymorphisms were detected from four partial candidate genes encoding aquaporin and Na⁺/H⁺ antiporter in both plasma and tonoplast membranes. Overall, rapid decay of linkage disequilibrium was observed within 500 bp. Significant associations were found between the putative LpTIP1 and Na⁺ for the control and between the putative LpNHX1 and K⁺/Na⁺ under the control and salinity treatments after controlling population structure. These results indicate that population structure influenced phenotypic traits, and allelic variation in LpNHX1 may affect salinity tolerance of perennial ryegrass.     

Genome-wide local ancestry approach identifies genes and variants associated with chemotherapeutic susceptibility in African Americans

Chemotherapeutic agents are used in the treatment of many cancers, yet variable resistance and toxicities among individuals limit successful outcomes. Several studies have indicated outcome differences associated with ancestry among patients with various cancer types. Using both traditional SNP-based and newly developed gene-based genome-wide approaches, we investigated the genetics of chemotherapeutic susceptibility in lymphoblastoid cell lines derived from 83 African Americans, a population for which there is a disparity in the number of genome-wide studies performed. To account for population structure in this admixed population, we incorporated local ancestry information into our association model. We tested over 2 million SNPs and identified 325, 176, 240, and 190 SNPs that were suggestively associated with cytarabine-, 5'-deoxyfluorouridine (5'-DFUR)-, carboplatin-, and cisplatin-induced cytotoxicity, respectively (p≤10(-4)). Importantly, some of these variants are found only in populations of African descent. We also show that cisplatin-susceptibility SNPs are enriched for carboplatin-susceptibility SNPs. Using a gene-based genome-wide association approach, we identified 26, 11, 20, and 41 suggestive candidate genes for association with cytarabine-, 5'-DFUR-, carboplatin-, and cisplatin-induced cytotoxicity, respectively (p≤10(-3)). Fourteen of these genes showed evidence of association with their respective chemotherapeutic phenotypes in the Yoruba from Ibadan, Nigeria (p<0.05), including TP53I11, COPS5 and GAS8, which are known to be involved in tumorigenesis. Although our results require further study, we have identified variants and genes associated with chemotherapeutic susceptibility in African Americans by using an approach that incorporates local ancestry information.     

Whole-genome resequencing identified candidate genes associated with the number of ribs in Hu sheep

The number of ribs is an important economic trait in the sheep industry when the sheep are raised for mutton. However, in sheep, the genetic mechanisms regulating rib number are poorly understood. In the present study, we aimed to identify important candidate genes that affect the increase in rib number in sheep. Whole-genome resequencing of 36 Hu sheep with an increased number of ribs (R14) and 36 sheep with normal (R13) rib numbers was carried out. Analysis using three methods (fixation index (F_ST), Fisher's exact test, and Chi-squared test) showed that 219 single nucleotide polymorphism sites overlapped among the results of the three methods, which represented 206 genes. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses indicated that the genes were mainly associated with regulation of developmental process, inorganic anion transport, cellular biosynthetic process, tight junction, the oxytocin signaling pathway, and arrhythmogenic right ventricular cardiomyopathy. Four mutations were selected according to the significantly selected genomic regions and important pathways for genotyping and association analysis. The result demonstrated that three synonymous mutations correlated significantly with the rib number. Importantly, we revealed that the CPOX (encoding coproporphyrinogen oxidase), KCNH1 (encoding potassium voltage-gated channel subfamily H member 1), and CPQ (encoding carboxypeptidase Q) genes have a combined effect on rib number in Hu sheep. Our results identified candidate molecular markers for rib number in sheep breeding.     

Landscape genetic analyses reveal cryptic population structure and putative selection gradients in a large-scale estuarine environment

Disentangling the relative contributions of selective and neutral processes underlying phenotypic and genetic variation under natural, environmental conditions remains a central challenge in evolutionary ecology. However, much of the variation that could be informative in this area of research is likely to be cryptic in nature; thus, the identification of wild populations suitable for study may be problematic. We use a landscape genetics approach to identify such populations of three-spined stickleback inhabiting the Saint Lawrence River estuary. We sampled 1865 adult fish over multiple years. Individuals were genotyped for nine microsatellite loci, and georeferenced multilocus data were used to infer population groupings, as well as locations of genetic discontinuities, under a Bayesian model framework ( geneland ). We modelled environmental data using nonparametric multiple regression to explain genetic differentiation as a function of spatio-ecological effects. Additionally, we used genotype data to estimate dispersal and gene flow to parameterize a simple model predicting adaptive vs. plastic divergence between demes. We demonstrate a bipartite division of the genetic landscape into freshwater and maritime zones, independent of geographical distance. Moreover, we show that the greatest proportion of genetic variation (31.5%) is explained by environmental differences. However, the potential for either adaptive or plastic divergence between demes is highly dependent upon the strength of migration and selection. Consequently, we highlight the utility of landscape genetics as a tool for hypothesis generation and experimental design, to identify focal populations and putative selection gradients, in order to distinguish between phenotypic plasticity and local adaptation.     

Evolutionary analyses of the major variant surface antigen-encoding genes reveal population structure of Plasmodium falciparum within and between continents

Malaria remains a major public health problem in many countries. Unlike influenza and HIV, where diversity in immunodominant surface antigens is understood geographically to inform disease surveillance, relatively little is known about the global population structure of PfEMP1, the major variant surface antigen of the malaria parasite Plasmodium falciparum. The complexity of the var multigene family that encodes PfEMP1 and that diversifies by recombination, has so far precluded its use in malaria surveillance. Recent studies have demonstrated that cost-effective deep sequencing of the region of var genes encoding the PfEMP1 DBLα domain and subsequent classification of within host sequences at 96% identity to define unique DBLα types, can reveal structure and strain dynamics within countries. However, to date there has not been a comprehensive comparison of these DBLα types between countries. By leveraging a bioinformatic approach (jumping hidden Markov model) designed specifically for the analysis of recombination within var genes and applying it to a dataset of DBLα types from 10 countries, we are able to describe population structure of DBLα types at the global scale. The sensitivity of the approach allows for the comparison of the global dataset to ape samples of Plasmodium Laverania species. Our analyses show that the evolution of the parasite population emerging out of Africa underlies current patterns of DBLα type diversity. Most importantly, we can distinguish geographic population structure within Africa between Gabon and Ghana in West Africa and Uganda in East Africa. Our evolutionary findings have translational implications in the context of globalization. Firstly, DBLα type diversity can provide a simple diagnostic framework for geographic surveillance of the rapidly evolving transmission dynamics of P. falciparum. It can also inform efforts to understand the presence or absence of global, regional and local population immunity to major surface antigen variants. Additionally, we identify a number of highly conserved DBLα types that are present globally that may be of biological significance and warrant further characterization.     

Genome-wide association study revealed a promising region and candidate genes for eggshell quality in an F2 resource population

Background

Eggshell is subject to quality loss with aging process of laying hens, and damaged eggshells result in economic losses of eggs. However, the genetic architecture underlying the dynamic eggshell quality remains elusive. Here, we measured eggshell quality traits, including eggshell weight (ESW), eggshell thickness (EST) and eggshell strength (ESS) at 11 time points from onset of laying to 72 weeks of age and conducted comprehensive genome-wide association studies (GWAS) in 1534 F₂ hens derived from reciprocal crosses between White Leghorn (WL) and Dongxiang chickens (DX).

Results

ESWs at all ages exhibited moderate SNP-based heritability estimates (0.30 ~ 0.46), while the estimates for EST (0.21 ~ 0.31) and ESS (0.20 ~ 0.27) were relatively low. Eleven independent univariate genome-wide screens for each trait totally identified 1059, 1026 and 1356 significant associations with ESW, EST and ESS, respectively. Most significant loci were in a region spanning from 57.3 to 71.4 Mb of chromosome 1 (GGA1), which together account for 8.4 ~ 16.5 % of the phenotypic variance for ESW from 32 to 72 weeks of age, 4.1 ~ 6.9 % and 2.95 ~ 16.1 % for EST and ESS from 40 to 72 weeks of age. According to linkage disequilibrium (LD) and conditional analysis, the significant SNPs in this region were in extremely strong linkage disequilibrium status. Ultimately, two missense SNPs in GGA1 and one in GGA4 were considered as promising loci on three independent genes including ITPR2, PIK3C2G, and NCAPG. The homozygotes of advantageously effective alleles on PIK3C2G and ITPR2 possessed the best eggshell quality and could partly counteract the negative effect of aging process. NCAPG had certain effect on eggshell quality for young hens.

Conclusions

Identification of the promising region as well as potential candidate genes will greatly advance our understanding of the genetic basis underlying dynamic eggshell quality and has the practical significance in breeding program for the improvement of eggshell quality, especially at the later part of laying cycle.

Electronic supplementary material

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

Differential expression of six genes and correlation with fatness traits in a unique broiler population

Previous results from genome wide association studies (GWASs) in chickens divergently selected for abdominal fat content of Northeast Agricultural University (NEAUHLF) showed that many single nucleotide polymorphism (SNP) variants were associated with abdominal fat content. Of them, six top significant SNPs at the genome level were located within SRD5A3, SGCZ, DLC1, GBE1, GALNT9 and DNAJB6 genes. Here, expression levels of these six candidate genes were investigated in abdominal fat and liver tissue between fat and lean broilers from the 14th generation population of NEAUHLF. The results showed that expression levels of SRD5A3, SGCZ and DNAJB6 in the abdominal fat and SRD5A3, DLC1, GALNT9, DNAJB6 and GBE1 in the liver tissue differed significantly between the fat and lean birds, and were correlated with abdominal fat traits. The findings will provide important references for further function investigation of the six candidate genes involved in abdominal fat deposition in chickens.     

Genome-wide scans identify known and novel regions associated with prolificacy and reproduction traits in a sub-Saharan African indigenous sheep (Ovis aries)

Maximizing the number of offspring born per female is a key functionality trait in commercial- and/or subsistence-oriented livestock enterprises. Although the number of offspring born is closely associated with female fertility and reproductive success, the genetic control of these traits remains poorly understood in sub-Saharan Africa livestock. Using selection signature analysis performed on Ovine HD BeadChip data from the prolific Bonga sheep in Ethiopia, 41 candidate regions under selection were identified. The analysis revealed one strong selection signature on a candidate region on chromosome X spanning BMP15, suggesting this to be the primary candidate prolificacy gene in the breed. The analysis also identified several candidate regions spanning genes not reported before in prolific sheep but underlying fertility and reproduction in other species. The genes associated with female reproduction traits included SPOCK1 (age at first oestrus), GPR173 (mediator of ovarian cyclicity), HB-EGF (signalling early pregnancy success) and SMARCAL1 and HMGN3a (regulate gene expression during embryogenesis). The genes involved in male reproduction were FOXJ1 (sperm function and successful fertilization) and NME5 (spermatogenesis). We also observed genes such as PKD2L2, MAGED1 and KDM3B, which have been associated with diverse fertility traits in both sexes of other species. The results confirm the complexity of the genetic mechanisms underlying reproduction while suggesting that prolificacy in the Bonga sheep, and possibly African indigenous sheep is partly under the control of BMP15 while other genes that enhance male and female fertility are essential for reproductive fitness.

     

Comparative population genomics identified genomic regions and candidate genes associated with fruit domestication traits in peach

Crop evolution is a long‐term process involving selection by natural evolutionary forces and anthropogenic influences; however, the genetic mechanisms underlying the domestication and improvement of fruit crops have not been well studied to date. Here, we performed a population structure analysis in peach (Prunus persica) based on the genome‐wide resequencing of 418 accessions and confirmed the presence of an obvious domestication event during evolution. We identified 132 and 106 selective sweeps associated with domestication and improvement, respectively. Analysis of their tissue‐specific expression patterns indicated that the up‐regulation of selection genes during domestication occurred mostly in fruit and seeds as opposed to other organs. However, during the improvement stage, more up‐regulated selection genes were identified in leaves and seeds than in the other organs. Genome‐wide association studies (GWAS) using 4.24 million single nucleotide polymorphisms (SNPs) revealed 171 loci associated with 26 fruit domestication traits. Among these loci, three candidate genes were highly associated with fruit weight and the sorbitol and catechin content in fruit. We demonstrated that as the allele frequency of the SNPs associated with high polyphenol composition decreased during peach evolution, alleles associated with high sugar content increased significantly. This indicates that there is genetic potential for the breeding of more nutritious fruit with enhanced bioactive polyphenols without disturbing a harmonious sugar and acid balance by crossing with wild species. This study also describes the development of the genomic resources necessary for evolutionary research in peach and provides the large‐scale characterization of key agronomic traits in this crop species.