A Genome-Wide Association Study Identifies Two Novel Promising Candidate Genes Affecting Escherichia coli F4ab/F4ac Susceptibility in Swine

Abstract

Enterotoxigenic Escherichia coli (ETEC) expressing F4 fimbria is the major pathogenic bacteria causing diarrhoea in neonatal and post-weaning piglets. Previous studies have revealed that the susceptibility to ETEC F4ab/F4ac is an autosomal Mendelian dominant trait and the loci controlling the F4ab/F4ac receptor are located on SSC13q41, between markers SW207 and S0283. To pinpoint these loci and further validate previous findings, we performed a genome-wide association study (GWAS) using a two generation family-based population, consisting of 301 piglets with phenotypes of susceptibility to ETEC F4ab/F4ac by the vitro adhesion test. The DNA of all piglets and their parents was genotyped using the Illumina PorcineSNP60 BeadChip, and 50,972 and 50,483 SNPs were available for F4ab and F4ac susceptibility, respectively, in the association analysis after quality control. In summary, 28 and 18 significant SNPs (p<0.05) were detected associated with F4ab and F4ac susceptibility respectively at genome-wide significance level. From these significant findings, two novel candidate genes, HEG1 and ITGB5, were firstly identified as the most promising genes underlying F4ab/F4ac susceptibility in swine according to their functions and positions. Our findings herein provide a novel evidence for unravelling genetic mechanism of diarrhoea risk in piglets.

About the Authors:

Wei-Xuan Fu

Contributed equally to this work with: Wei-Xuan Fu, Yang Liu

Affiliation: Key Laboratory of Animal Genetics Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, China

Yang Liu

Contributed equally to this work with: Wei-Xuan Fu, Yang Liu

Affiliation: Key Laboratory of Animal Genetics Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, China

Xin Lu

Affiliations Key Laboratory of Animal Genetics Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, China, State Key Laboratory for Infectious Disease Prevention and Control, National Institute for Communicable Disease Control and Prevention, Chinese Centre for Disease Control and Prevention, Beijing, China

Xiao-Yan Niu

Affiliation: Key Laboratory of Animal Genetics Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, China

Xiang-Dong Ding

Affiliation: Key Laboratory of Animal Genetics Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, China

Jian-Feng Liu

* E-mail: [email protected] (J-FL); [email protected] (QZ)

Affiliation: Key Laboratory of Animal Genetics Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, China

Qin Zhang

* E-mail: [email protected] (J-FL); [email protected] (QZ)

Affiliation: Key Laboratory of Animal Genetics Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, China

Introduction

Susceptibility to enterotoxigenic Escherichia (E.) coli (ETEC) with F4 (K88) fimbriae is dominantly inherited in neonatal and pre-weaning piglets, potentially causing the diarrhoea and death of piglets. Antigenically, three fimbriae variants have been identified, i.e., F4ab, F4ac and F4ad [1]. Diarrhoea due to ETEC F4 (mainly F4ab and F4ac) infection is very common over the world where pigs are raised in large numbers [2]. Once the bacteria adhere to the brush border of enterocytes and colonize the small intestine, the enterotoxins produced by them induce an increased secretion of electrolytes into lumen, resulting in diarrhoea. However, not all piglets are susceptible to ETEC F4 and the susceptibility is determined by the existence of a specific receptor of ETEC F4ab/F4ac [3].

Identification of causal mutation(s) affecting F4ab/F4ac susceptibility is a feasible way for controlling piglet diarrhoea. In the past few years, several groups have focused on mapping potential genes of F4ab/F4ac receptor (F4abR/F4acR), which has been mapped on SSC13 initially [4]. Python et al. refined the localization of the F4acR gene in the interval of S0068 and SW1030, close to S0075 and SW225 [5], and four functional genes in the region between SW207 and S0283 were considered as candidates of F4acR in their further research, but no mutations were found in the cDNA sequences of these genes associated with the F4acR genotypes [6]. A linkage and comparative mapping study demonstrated that the most likely region of F4abR/F4acR gene(s) was between SW207 and S0075 [7]. Studies [3], [8], [9] suggested that the MUC4 (mucin 4) gene should be considered as one of the most promising candidate genes for F4abR/F4acR based on three aspects: (1) its physical position falls into the region identified in most of linkage analyses, (2) the protein encoded by the MUC4 gene is one of mucin-like sialoglycoprotein which the ETEC F4 maybe adheres to, and (3) association studies showed strong association between a mutation within the gene and the adhesive phenotypes. However, in a recent study, Rampoldi et al. found that the region of F4abR/F4acR gene(s) was distal to the interval between the MUC4 gene and the LMLN gene (leishmanolysin-like gene) through testing recombination events in three–generation pedigree [10]. Findings from these studies suggest that further endeavours are still expected to explore more convincing evidences on the F4abR/F4acR gene(s).

Although QTL mapping has been very successful in domestic animals for many complex quantitative traits and a few prominent findings have been reported [11], [12], [13], [14], identification of quantitative trait mutations (QTMs) is still a challenge due to the commonly existing limitations of linkage analysis [15]. Recently, the first high-density 60 K porcine SNP array has been developed [16], which offers the prerequisite of genome-wide association study (GWAS) in swine, a powerful approach for high-resolution mapping of loci controlling phenotypic traits. Duijvesteijn et al. [17] reported a GWAS for androstenone levels in pigs, which reveals a cluster of candidate genes on chromosome 6. Moreover, Pryce et al. [18] and Jiang et al. [19] performed GWAS to identify loci affecting milk production traits in dairy cattle in two different populations and obtained very similar results. More recently, Fan et al. [20] performed a GWAS for body composition and structural soundness traits of pigs and identified several genes by functional clustering analysis. Up to now, GWAS has acted as a most commonly used strategy for gene identification for complex traits in animals as well as humans.

Motivated by further clarifying genetic basis of ETEC F4ab/F4ac susceptibility and pursuing more confirmatory evidences of F4abR/F4acR gene(s), we performed a GWAS for high-resolution identification of loci controlling F4ab/F4acsusceptibility in swine using a case-control design. Our study identified 28 and 18 significant SNPs for ETEC F4ab/F4ac susceptibility, respectively. These significant findings indicate two genes, HEG1 and ITGB5, can be treated as novel promising candidates underlying F4ab/F4ac susceptibility in swine according to their functions and positions.

Materials and Methods

Animal resource

The animal resource used for this study is a two generation family-based population, which is a subset of the population of our previous study [21], consisting of 301 pure bred piglets of three breeds, among which 67 were Landrace offspring of 4 boars and 13 sows, 161 were Yorkshire offspring of 7 boars and 29 sows, and 73 were Songliao Black (a Chinese native breed) offspring of 3 boars and 13 sows (see Table 1). The pigs were raised under standard indoor conditions at the experimental farm of the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences. At 35 days of age, all piglets were slaughtered, and their jejunum and ear tissue samples were collected. The ear tissue samples of 69 parents recorded in the original pedigree were also collected. In total, there were 370 pigs in our study, including 301 piglets and 69 parental individuals.

[Figure omitted. See PDF.]

Table 1. Family structure in three swine breeds.

https://doi.org/10.1371/journal.pone.0032127.t001

The whole procedure for the collection of jejunum and ear tissue samples and the slaughter of piglets was carried out in strict accordance with the protocol approved by the Animal Welfare Committee of China Agricultural University (Permit Number: DK996).

Measurement of phenotypes

All 301 piglets were phenotyped for ETEC F4ab/F4ac susceptibility using in vitro adhesion test. Two ETEC F4 strains (195 (F4ab, C83901, O8:K87) and 200 (F4ac, C83907, O149:K91)) and a bovine-origined E. coli strain (238 (C83286, O38:K99)) as negative control were provided by the China Institute of Veterinary Drug Control. The procedures of collecting jejunal epithelial cells, preparing bacterial suspension, in vitro adhesion test and classification of adhesion phenotypes (strongly adhesive, adhesive, weakly adhesive and non-adhesive) herein were described in detail in our previous report [21].

For the case-control design in our GWAS, the four phenotype categories were further classified into two classes. The non-adhesive and weakly adhesive phenotypes were classified to as negative (control) and the other two categories as positive (case). The distribution of the two classes in the three breeds is given in Table 2.

[Figure omitted. See PDF.]

Table 2. Distribution of F4ab/F4ac adhesion phenotypes in three swine breeds.

https://doi.org/10.1371/journal.pone.0032127.t002

Genotyping and quality control

Extracted from ear tissue samples of all piglets and their parents, DNA was genotyped using the Illumina PorcineSNP60 BeadChip containing 62,163 SNPs. Features of the chip have been detailed previously [16]. The genotypes were judged using BeadStudio (Version 3.2.2, Illumina, lnc.) and a custom cluster file developed from the 370 samples.

To assess the technical reliability of the genotyping panel, a randomly selected DNA sample was genotyped twice, and over 99% identity of called genotypes was obtained. This demonstrates the technically robust feature of the 60 K SNP BeadChip panel employed herein.

Like most GWAS using the case-control design [22], [23], quality control procedures of the genotype data were performed as follows. First, only samples with a minimum of 90% call rate were included. Second, out of the initial full set of 62,163 SNPs, we discarded: (1) SNPs with a call rate <90% (n = 3,773) in all piglets; (2) those severely deviating from Hardy–Weinberg equilibrium (HWE) (p<10E-6) in the two control groups (n = 2,179 for F4ab and n = 2,693 for F4ac); and (3) those having a minor allele frequency (MAF)<0.03 in all piglets (n = 7,797). After quality control, 50,972 and 50,483 SNPs were available for F4ab and F4ac respectively in the subsequent analyses, and their distributions in the porcine genome are presented in Table 3.

[Figure omitted. See PDF.]

Table 3. Distribution of SNPs on chromosomes after quality control and the average distances between adjacent SNPs.

https://doi.org/10.1371/journal.pone.0032127.t003

Parentage test

Considering the probability of potential parentage mistakes in the original parentage records, we adopted Cervus (Version 3.0) [24] to infer the most possible parent-offspring pairs with maximum likelihood method using 200 randomly chosen autosomal SNPs with more than 99% call rate. Accordingly, among all 301 piglets, a total of 50 (16.6%) individuals had parentage errors in the original records, including 11 with incorrect paternal records, 31 with incorrect maternal records and 8 with both incorrect paternal and maternal records. Since the parentage information was to be used in the association analysis, parentage correction was further conducted. Among the 50 piglets with parentage errors, 9 were reassigned to the correct parents among the 69 known parents, while 41 were unable to be assigned to any of the known fathers or mothers. Assuming correct sibship information in the original parentage records, they were assigned to be offspring of 14 unknown parents including 3 boars and 11 sows. Hence, in our study, the 301 piglets were actually from 83 parents. The corrected parentage information was used in the subsequent analyses. The information of both original and corrected parentage is given in Table 1.

Association Analysis

Compared with traditional population-based case-control design in GWAS, individuals in cases (piglets with adhesive phenotypes) and controls (piglets with non-adhesive phenotypes) in our studies are related within each breed, and heterogeneity also exists among three different breeds, which may potentially induce confounding in the analysis. To overcome this limitation, we used a recently published program ROADTRIPS (Version 1.2) [25] to perform the association analysis. An important advantage of ROADTRIP is that it can simultaneously deal with data with pedigree structure as well as population admixture in association test. In ROADTRIP, an empirical covariance matrix Ψ constructed using genome-wide SNP data is employed to adjust for potential population admixture as well as relatedness among individuals, while maintaining the advantage of utilizing known pedigree information when it is available.

ROADTRIPS provides three association tests named RM test, RW test and Rχ test, respectively. Compared with the RW and Rχ tests, the RM test can use the phenotypic information for individuals with missing genotypes provided they have a sampled relative who is genotyped at the tested marker. The RM and RW tests can improve power by using this information when partial or complete pedigree information is available. Furthermore, the RM test is the most powerful in a general class of linear statistics under the framework of two-allele disease model for outbred populations.

Considering features of the RM test aforementioned as well as the data structure of our study based on the corrected pedigree above, we adopted the test to detect loci associated with susceptibility to ETEC F4ab/F4ac, and the p values for the RM statistic were derived from an asymptotic chi-square distribution with 1 degree of freedom.

For the significant SNPs detected by the RM test, linkage disequilibrium (LD) patterns between them were quantified as r2 using Haploview (Version 4.2) [26] and the LD blocks were defined by the criteria of Gabriel et al. [27].

Statistical Inference

In the study, the permutation method was adopted to adjust for multiple testing for the number of SNPs tested through constructing a genome-wide empirical distribution of the RM statistic under null hypothesis. The phenotypes of ETEC F4ab/F4ac susceptibility were randomly shuffled 10,000 times; and the empirical critical value was determined by choosing the 95th percentile of the highest test statistic over the 10,000 permutation replicates. We declared a significant SNP at a genome-wide 0.05 significance level if its RM statistic value was larger than the empirical critical value.

Results

The profiles of the p values (in terms of −log10 p) of all tested SNPs for susceptibility to ETEC F4ab/F4ac are shown in Figure 1. The genome-wide significant SNPs detected by the RM test for ETEC F4ab/F4ac susceptibility at the permutation-based critical level are presented in Table 4. In total, 28 and 18 genome-wide significant SNPs (p<0.05) were detected for susceptibility to ETEC F4ab and F4ac, respectively, and all of the 18 significant SNPs for F4ac are also significant for F4ab. Furthermore, based on the most recent porcine genome sequence assembly (Sscrofa9.2), 18 SNPs (15 of them are significant for both F4ab and F4ac) among the 28 significant SNPs are located within an interval of about 2.6 Mb on SSC13, while the positions of other 10 SNPs are not available (see Table 4).

[Figure omitted. See PDF.]

Figure 1. Manhattan plots of genome-wide association for ETEC F4ab/F4ac susceptibility.

Negative log10-transformed p values of all tested SNPs for susceptibility to ETEC F4ab (Fig. 1A) and F4ac (Fig. 1B) are plotted against position on each of 19 chromosomes. Different chromosomes are represented by different colours. Chr. 19 stands for the X chromosome of swine.

https://doi.org/10.1371/journal.pone.0032127.g001

[Figure omitted. See PDF.]

Table 4. Genome-wide significant SNPs for ETEC F4ab and F4ac susceptibility.

https://doi.org/10.1371/journal.pone.0032127.t004

The LD patterns among the 18 significant SNPs with known positions are shown in Figure 2. Four LD blocks were defined with the criteria of Gabriel et al. [27]. Outside of these four blocks, there is merely one significant SNP (H3GA0037388) located about 972 kb away from its nearest significant SNP, which could be a long-distance LD marker, i.e., it has strong LD with but a long physical distance from the causal mutation(s) of ETEC F4ab/F4ac susceptibility.

[Figure omitted. See PDF.]

Figure 2. Linkage disequilibrium (LD) pattern for common significant SNPs on SSC13.

LD blocks are marked with triangles. Values in boxes are LD (r2) between SNP pairs and the boxes are coloured according to the standard Haploview colour scheme: LOD>2 and D′ = 1, red; LOD>2 and D′<1, shades of pink/red; LOD<2 and D′ = 1, blue; LOD<2 and D′<1, white (LOD is the log of the likelihood odds ratio, a measure of confidence in the value of D′).

https://doi.org/10.1371/journal.pone.0032127.g002

Discussion

GWAS has been considered as a promising tool for gene identification for complex traits. So far GWAS in domestic animals are largely focused on economically important growth and production traits, such as milk production in dairy cattle, backfat in swine, etc. In this study, we carried out a GWAS to explore potential causal gene(s) for F4abR and F4acR in swine. To our knowledge, this is the first study aiming at unravelling the genetic mechanism of the ETEC F4ab/F4ac susceptibility in piglets using a case-control design based on a high density SNP chip panel. Findings from our study will lay a preliminary foundation for the follow-up functional validation of F4abR/F4acR candidate gene(s) in swine.

In the vitro adhesion test of our study, misclassification may occur for the weakly adhesive animals but seldom for strongly adhesive, adhesive and non-adhesive animals. Billey et al. [28] considered the weak adhesion was an artefact because it was rarely detected. This was the same case in our study, i.e., the number of weakly adhesive piglets was merely 11 for F4ab, and 24 for F4ac. Python et al. suggested that the receptor for weak adhesion might be different from that for normal and strong adhesion [5], [6]. Hence it was reasonable to treat the weak adhesion as non-adhesion in contrast to the adhesion and strong adhesion. To further examine potential false-positive/false-negative errors raised from misclassification, we performed a GWAS using a subset of original samples without the weakly adhesive ones, and the association results (data unpresented) totally kept unchanged compared with those showed herein.

In the study, according to the critical value determined by permutation tests via 10,000 replicates, 18 genome-wide significant SNPs with known positions in porcine genome were identified for susceptibility to ETEC F4ab/F4ac. The positions of the 18 SNPs are consistent with previously reported QTL regions [3], [9]. These results offer a high possibility that F4abR and F4acR share the common causal mutation(s) as suggested by some previous studies [5], [7], [29]. Similar with the findings by Rampoldi et al. [10], no significant SNP was detected in the MUC4-LMLN region.

The LD patterns of the 18 significant SNPs for ETEC F4ab/F4ac showed that almost all of them are in high LD (r2) level even though there are long distances between some of them (see Figure 2). Further scrutinizing their positions, we found that they are located in either intergenic regions (n = 11) or intronic regions (n = 7) (see Table 4). Under the assumption that the causal mutation(s) are in strong LD with these significant SNPs, we are in effort to tag potential functional genes within the LD regions covered by these significant SNPs. This is also the common strategy adopted by a suite of prior similar studies [6], [30], [31]. The promising candidate genes were then determined in terms of their known functional information in other species. To further pinpoint the true causal mutation(s), follow-up functional validations should be performed by focusing on mining all possible mutations in coding or non-coding regions of these candidate genes.

Based on the most recent porcine genome sequence assembly (Sscrofa9.2, http://www.ensembl.org/Sus_scrofa/Info/Index), one SNP, ASGA0058885, is located within an intron of TNK2 (tyrosine kinase non-receptor 2) and another two SNPs, MARC0012378 and M1GA0017682, are very close to TNK2 (6127 bp and 48 bp away, respectively); Three SNPs, MARC0099692, ALGA0072075 (the most significant SNP) and MARC0002946, are located between ZNF148 (zinc finger protein 148) and HEG1 (HEG homolog 1); Two SNPs, ASGA0089965 and ASGA0091537, are located between HEG1 and MUC13 (mucin 13); Three SNPs, MARC0067282, ASGA0058925 and ALGA0106330, are located in introns of ZNF148, HEG1 and MUC13, respectively; Two SNPs, H3GA0037348 and H3GA0037351, are located in introns of ITGB5 (integrin beta-5). Based on these findings, we could further focus on these genes involved, i.e., TNK2, ZNF148, HEG1, MUC13 and ITGB5.

Among these five genes, the association of TNK2 and MUC13 with susceptibility to ETEC F4ab/F4ac were analysed in previous studies [30], [32]. TNK2 encodes a tyrosine kinase that binds to CDC42 (cell division cycle 42 protein) in its GTP-bound form is inactivated by intrinsic hydrolysis of the nucleotide γ-phosphate, which can be stimulated by GTPase-activating proteins (GAPs) [33]. In humans, several alternatively spliced transcript variants have been identified from TNK2, and the full-length of two transcript variants has been determined. The homologue swine TNK2 gene has a genomic length of about 42 kb and is located in the region of 100,885–100,927 kb on SSC13, which is very close to MUC4, so it makes TNK2 to act as a possible candidate gene similar to MUC4 based on its position as well as its functional aspect. However, according to the results of Joller [32], the sequence variants of both MUC4 and TNK2 were not completely linked to the phenotypes, and none of them was causative for susceptibility to ETEC F4ab/F4ac. A similar result was reported by Rampoldi et al. [10], who also suggested that the causal mutation(s) was downstream of the gene LMLN and might be located around the region containing the MUC13 gene. Therefore, TNK2 should not be considered as a candidate gene for F4abR/F4acR in further research.

Similar to MUC4, MUC13 belongs to the family of secreted and cell surface glycoproteins expressed by ductal and glandular epithelial tissues [34] and plays a role in cell signalling. MUC13 should be a highly possible candidate gene having the causal mutation(s) for susceptibility to ETEC F4ab/F4ac since MUC4 was denied with quite reasonable experimental evidences [10]. Zhang et al. [30] showed that the expression pattern of the porcine MUC13 mRNA in tissues was similar to humans, with highest level in jejunum and moderate levels in trachea, stomach and liver, and the SNPs detected in MUC13 were strongly associated with susceptibility to ETEC F4ab/F4ac in a White Duroc×Erhualian resource population in their initial study. Subsequently, MUC13 was assigned as a positional candidate gene for F4abR/F4acR via the study about a pig–human comparative radiation hybrid (RH) map [31]. However, no causal mutation can be identified in MUC13 so far.

The other three genes, ZNF148, HEG1 and ITGB5, have not been treated as candidate genes for F4abR/F4acR so far. ZNF148 has been shown to be involved in regulation of T cell receptors in human [35], but no evidence is available to indicate its functional relationship with bacterial infection or molecular receptor on epithelial cells. We are strongly in favour of HEG1 and ITGB5 as potential candidate genes for F4abR/F4acR based the following reasons.

From the physical positions of these two genes, MUC13 is flanked by ITGB5 and HEG1, each transcribed from the forward strand. Interestingly, HEG1 and MUC13 were proved sharing same molecular features, suggesting they might be evolutionarily related [36]. However, the expression pattern of HEG1 in human is quite different from that of MUC13 based on the information from the BioGPS (http://biogps.org), which showed a low expression level for HEG1 but the highest expression level for MUC13 in small intestine. The functional information of this gene is limited so far, and it is not clear whether it is possible that HEG1 has the causal mutation(s) responsible for the ETEC F4ab/F4ac susceptibility or not. In spite of these facts above, our study revealed three significant SNPs (including the most significant one) close to HEG1 and one within it, which suggests that HEG1 could be considered as a candidate gene and merit follow-up validation in the future.

Like HEG1, based on the same database, the expression level of ITGB5 is lower than that of MUC13 in small intestine of human, but it is not an essential issue because the expression pattern is not stable in different ages or species. And, some functional information of ITGB5 proved in previous studies shows its potentiality of being F4abR/F4acR gene. The protein type of ITGB5 belongs to the integrin beta chain family, and is associated with alpha-V for compounding integrin αVβ5, which plays an important role in the innate defence system against bacterial infection by influencing the rapid turnover and exfoliation of mucosal epithelial cells [37]. In human, αVβ5 is a major endocytic receptor for vitronectin (Vn) which has an Arg-Gly-Asp (RGD) sequence for binding [38], [39], [40]. Vn plays an important role in bacterial serum resistance, adhesion and internalization mediated by host cell signalling. Furthermore, it has distinct binding sites for pathogens and epithelial cells like a cross-link between bacteria and epithelial cells [41]. In addition, Vn bound to E. coli, Staph. aureus and S. pneumoniae provides a more efficient bacterial adhesion to epithelial cells [42]. Therefore, it is extremely possible that a potential mutation in porcine ITGB5 could affect integrin αVβ5 binding to Vn, and as a result its accompanying E. coli could not adhere to jejunal epithelial cells. Moreover, it has been reported that the polymorphisms of ITGB5 is the host factor which might affect adenovirus infection and decrease lung function in human [43]. Additionally, the ITGB5 subunit was found on both the apical and basal surface of epithelial cells and its expression is essentially oestrous cycle-independent in mice [44]. Based on the evidences above and the results in our study, we confirm that ITGB5 is the most possible functional candidate gene for the F4abR/F4acR.

In this study, individuals from three different swine breeds were involved in the association analysis, so the issue of population stratification should be a major concern which needed to be addressed. Since the RM test itself is immune to population stratification [45], [46], it is safe to assume that the SNPs detected have convincing associations with ETEC F4ab/F4ac susceptibility in swine with a reasonable false discovery rate.

Acknowledgments

We thank visiting professor Geoffrey Kirton for his assistance. We are grateful to two anonymous reviewers for their valuable comments which have greatly improved the presentation of the manuscript.

Author Contributions

Conceived and designed the experiments: WXF YL XYN JFL QZ. Performed the experiments: WXF YL XL XYN. Analyzed the data: WXF YL XL XDD JFL. Contributed reagents/materials/analysis tools: WXF YL XDD JFL QZ. Wrote the paper: WXF YL JFL QZ.

Citation: Fu W-X, Liu Y, Lu X, Niu X-Y, Ding X-D, Liu J-F, et al. (2012) A Genome-Wide Association Study Identifies Two Novel Promising Candidate Genes Affecting Escherichia coli F4ab/F4ac Susceptibility in Swine. PLoS ONE7(3): e32127. https://doi.org/10.1371/journal.pone.0032127

References

1. Guinee PA, Jansen WH (1979) Behavior of Escherichia coli K antigens K88ab, K88ac, and K88ad in immunoelectrophoresis, double diffusion, and hemagglutination. Infect Immun 23: 700–705.PA GuineeWH Jansen1979Behavior of Escherichia coli K antigens K88ab, K88ac, and K88ad in immunoelectrophoresis, double diffusion, and hemagglutination.Infect Immun23700705

2. Marquardt RR, Jin LZ, Kim JW, Fang L, Frohlich AA, et al. (1999) Passive protective effect of egg-yolk antibodies against enterotoxigenic Escherichia coli K88+ infection in neonatal and early-weaned piglets. Fems Immunol Med Mic 23: 283–288.RR MarquardtLZ JinJW KimL. FangAA Frohlich1999Passive protective effect of egg-yolk antibodies against enterotoxigenic Escherichia coli K88+ infection in neonatal and early-weaned piglets.Fems Immunol Med Mic23283288

3. Joller D, Jorgensen CB, Bertschinger HU, Python P, Edfors I, et al. (2009) Refined localization of the Escherichia coli F4ab/F4ac receptor locus on pig chromosome 13. Anim Genet 40: 749–752.D. JollerCB JorgensenHU BertschingerP. PythonI. Edfors2009Refined localization of the Escherichia coli F4ab/F4ac receptor locus on pig chromosome 13.Anim Genet40749752

4. Edfors-Lilja I, Gustafsson U, Duval-Iflah Y, Ellergren H, Johansson M, et al. (1995) The porcine intestinal receptor for Escherichia coli K88ab, K88ac: regional localization on chromosome 13 and influence of IgG response to the K88 antigen. Anim Genet 26: 237–242.I. Edfors-LiljaU. GustafssonY. Duval-IflahH. EllergrenM. Johansson1995The porcine intestinal receptor for Escherichia coli K88ab, K88ac: regional localization on chromosome 13 and influence of IgG response to the K88 antigen.Anim Genet26237242

5. Python P, Jorg H, Neuenschwander S, Hagger C, Stricker C, et al. (2002) Fine-mapping of the intestinal receptor locus for enterotoxigenic Escherichia coli F4ac on porcine chromosome 13. Anim Genet 33: 441–447.P. PythonH. JorgS. NeuenschwanderC. HaggerC. Stricker2002Fine-mapping of the intestinal receptor locus for enterotoxigenic Escherichia coli F4ac on porcine chromosome 13.Anim Genet33441447

6. Python P, Jorg H, Neuenschwander S, Asai-Coakwell M, Hagger C, et al. (2005) Inheritance of the F4ab, F4ac and F4ad E. coli receptors in swine and examination of four candidate genes for F4acR. J Anim Breed Genet 122: Suppl 15–14.P. PythonH. JorgS. NeuenschwanderM. Asai-CoakwellC. Hagger2005Inheritance of the F4ab, F4ac and F4ad E. coli receptors in swine and examination of four candidate genes for F4acR.J Anim Breed Genet122Suppl 1514

7. Jorgensen CB, Cirera S, Anderson SI, Archibald AL, Raudsepp T, et al. (2003) Linkage and comparative mapping of the locus controlling susceptibility towards E.coli F4ab/ac diarrhoea in pigs. Cytogenet Genome Res 102: 157–162.CB JorgensenS. CireraSI AndersonAL ArchibaldT. Raudsepp2003Linkage and comparative mapping of the locus controlling susceptibility towards E.coli F4ab/ac diarrhoea in pigs.Cytogenet Genome Res102157162

8. Erickson AK, Baker DR, Bosworth BT, Casey TA, Benfield DA, et al. (1994) Characterization of porcine intestinal receptors for the K88ac fimbrial adhesin of Escherichia coli as mucin-type sialoglycoproteins. Infect Immun 62: 5404–5410.AK EricksonDR BakerBT BosworthTA CaseyDA Benfield1994Characterization of porcine intestinal receptors for the K88ac fimbrial adhesin of Escherichia coli as mucin-type sialoglycoproteins.Infect Immun6254045410

9. Jacobsen M, Kracht SS, Esteso G, Cirera S, Edfors I, et al. (2010) Refined candidate region specified by haplotype sharing for Escherichia coli F4ab/F4ac susceptibility alleles in pigs. Anim Genet 41: 21–25.M. JacobsenSS KrachtG. EstesoS. CireraI. Edfors2010Refined candidate region specified by haplotype sharing for Escherichia coli F4ab/F4ac susceptibility alleles in pigs.Anim Genet412125

10. Rampoldi A, Jacobsen MJ, Bertschinger HU, Joller D, Burgi E, et al. (2011) The receptor locus for Escherichia coli F4ab/F4ac in the pig maps distal to the MUC4-LMLN region. Mamm Genome 22: 122–129.A. RampoldiMJ JacobsenHU BertschingerD. JollerE. Burgi2011The receptor locus for Escherichia coli F4ab/F4ac in the pig maps distal to the MUC4-LMLN region.Mamm Genome22122129

11. Andersson L, Haley CS, Ellegren H, Knott SA, Johansson M, et al. (1994) Genetic mapping of quantitative trait loci for growth and fatness in pigs. Science 263: 1771–1774.L. AnderssonCS HaleyH. EllegrenSA KnottM. Johansson1994Genetic mapping of quantitative trait loci for growth and fatness in pigs.Science26317711774

12. Georges M, Nielsen D, Mackinnon M, Mishra A, Okimoto R, et al. (1995) Mapping quantitative trait loci controlling milk production in dairy cattle by exploiting progeny testing. Genetics 139: 907–920.M. GeorgesD. NielsenM. MackinnonA. MishraR. Okimoto1995Mapping quantitative trait loci controlling milk production in dairy cattle by exploiting progeny testing.Genetics139907920

13. de Koning DJ, Janss LL, Rattink AP, van Oers PA, de Vries BJ, et al. (1999) Detection of quantitative trait loci for backfat thickness and intramuscular fat content in pigs (Sus scrofa). Genetics 152: 1679–1690.DJ de KoningLL JanssAP RattinkPA van OersBJ de Vries1999Detection of quantitative trait loci for backfat thickness and intramuscular fat content in pigs (Sus scrofa).Genetics15216791690

14. Lagarrigue S, Pitel F, Carre W, Abasht B, Le Roy P, et al. (2006) Mapping quantitative trait loci affecting fatness and breast muscle weight in meat-type chicken lines divergently selected on abdominal fatness. Genet Sel Evol 38: 85–97.S. LagarrigueF. PitelW. CarreB. AbashtP. Le Roy2006Mapping quantitative trait loci affecting fatness and breast muscle weight in meat-type chicken lines divergently selected on abdominal fatness.Genet Sel Evol388597

15. Andersson L (2009) Genome-wide association analysis in domestic animals: a powerful approach for genetic dissection of trait loci. Genetica 136: 341–349.L. Andersson2009Genome-wide association analysis in domestic animals: a powerful approach for genetic dissection of trait loci.Genetica136341349

16. Ramos AM, Crooijmans RP, Affara NA, Amaral AJ, Archibald AL, et al. (2009) Design of a high density SNP genotyping assay in the pig using SNPs identified and characterized by next generation sequencing technology. PLoS One 4: e6524.AM RamosRP CrooijmansNA AffaraAJ AmaralAL Archibald2009Design of a high density SNP genotyping assay in the pig using SNPs identified and characterized by next generation sequencing technology.PLoS One4e6524

17. Duijvesteijn N, Knol EF, Merks JW, Crooijmans RP, Groenen MA, et al. (2010) A genome-wide association study on androstenone levels in pigs reveals a cluster of candidate genes on chromosome 6. BMC Genet 11: 42.N. DuijvesteijnEF KnolJW MerksRP CrooijmansMA Groenen2010A genome-wide association study on androstenone levels in pigs reveals a cluster of candidate genes on chromosome 6.BMC Genet1142

18. Pryce J, Bolormaa S, Chamberlain A, Bowman P, Savin K, et al. (2010) A validated genome-wide association study in 2 dairy cattle breeds for milk production and fertility traits using variable length haplotypes. J Dairy Sci 93: 3331–3345.J. PryceS. BolormaaA. ChamberlainP. BowmanK. Savin2010A validated genome-wide association study in 2 dairy cattle breeds for milk production and fertility traits using variable length haplotypes.J Dairy Sci9333313345

19. Jiang L, Liu J, Sun D, Ma P, Ding X, et al. (2010) Genome wide association studies for milk production traits in Chinese Holstein population. PLoS One 5: e13661.L. JiangJ. LiuD. SunP. MaX. Ding2010Genome wide association studies for milk production traits in Chinese Holstein population.PLoS One5e13661

20. Fan B, Onteru SK, Du ZQ, Garrick DJ, Stalder KJ, et al. (2011) Genome-wide association study identifies Loci for body composition and structural soundness traits in pigs. PLoS One 6: e14726.B. FanSK OnteruZQ DuDJ GarrickKJ Stalder2011Genome-wide association study identifies Loci for body composition and structural soundness traits in pigs.PLoS One6e14726

21. Li Y, Qiu X, Li H, Zhang Q (2007) Adhesive patterns of Escherichia coli F4 in piglets of three breeds. J Genet Genomics 34: 591–599.Y. LiX. QiuH. LiQ. Zhang2007Adhesive patterns of Escherichia coli F4 in piglets of three breeds.J Genet Genomics34591599

22. Burgner D, Davila S, Breunis WB, Ng SB, Li Y, et al. (2009) A genome-wide association study identifies novel and functionally related susceptibility loci for Kawasaki disease. PLoS Genet 5: e1000319.D. BurgnerS. DavilaWB BreunisSB NgY. Li2009A genome-wide association study identifies novel and functionally related susceptibility loci for Kawasaki disease.PLoS Genet5e1000319

23. Pongpanich M, Sullivan PF, Tzeng JY (2010) A quality control algorithm for filtering SNPs in genome-wide association studies. Bioinformatics 26: 1731–1737.M. PongpanichPF SullivanJY Tzeng2010A quality control algorithm for filtering SNPs in genome-wide association studies.Bioinformatics2617311737

24. Kalinowski ST, Taper ML, Marshall TC (2007) Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol 16: 1099–1106.ST KalinowskiML TaperTC Marshall2007Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment.Mol Ecol1610991106

25. Thornton T, McPeek MS (2010) ROADTRIPS: case-control association testing with partially or completely unknown population and pedigree structure. Am J Hum Genet 86: 172–184.T. ThorntonMS McPeek2010ROADTRIPS: case-control association testing with partially or completely unknown population and pedigree structure.Am J Hum Genet86172184

26. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–265.JC BarrettB. FryJ. MallerMJ Daly2005Haploview: analysis and visualization of LD and haplotype maps.Bioinformatics21263265

27. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, et al. (2002) The structure of haplotype blocks in the human genome. Science 296: 2225–2229.SB GabrielSF SchaffnerH. NguyenJM MooreJ. Roy2002The structure of haplotype blocks in the human genome.Science29622252229

28. Billey LO, Erickson AK, Francis DH (1998) Multiple receptors on porcine intestinal epithelial cells for the three variants of Escherichia coli K88 fimbrial adhesin. Vet Microbiol 59: 203–212.LO BilleyAK EricksonDH Francis1998Multiple receptors on porcine intestinal epithelial cells for the three variants of Escherichia coli K88 fimbrial adhesin.Vet Microbiol59203212

29. Bijlsma IG, Bouw J (1987) Inheritance of K88-mediated adhesion of Escherichia coli to jejunal brush borders in pigs: a genetic analysis. Vet Res Commun 11: 509–518.IG BijlsmaJ. Bouw1987Inheritance of K88-mediated adhesion of Escherichia coli to jejunal brush borders in pigs: a genetic analysis.Vet Res Commun11509518

30. Zhang B, Ren J, Yan X, Huang X, Ji H, et al. (2008) Investigation of the porcine MUC13 gene: isolation, expression, polymorphisms and strong association with susceptibility to enterotoxigenic Escherichia coli F4ab/ac. Anim Genet 39: 258–266.B. ZhangJ. RenX. YanX. HuangH. Ji2008Investigation of the porcine MUC13 gene: isolation, expression, polymorphisms and strong association with susceptibility to enterotoxigenic Escherichia coli F4ab/ac.Anim Genet39258266

31. Ren J, Tang H, Yan X, Huang X, Zhang B, et al. (2009) A pig-human comparative RH map comprising 20 genes on pig chromosome 13q41 that harbours the ETEC F4ac receptor locus. J Anim Breed Genet 126: 30–36.J. RenH. TangX. YanX. HuangB. Zhang2009A pig-human comparative RH map comprising 20 genes on pig chromosome 13q41 that harbours the ETEC F4ac receptor locus.J Anim Breed Genet1263036

32. Joller D (2009) Comparative molecular approaches to identify host determinants mediating adhesion of E. coli F4 strains in pigs. D. Joller2009Comparative molecular approaches to identify host determinants mediating adhesion of E. coli F4 strains in pigs.PhD dissertation. PhD dissertation.

33. Manser E, Leung T, Salihuddin H, Tan L, Lim L (1993) A non-receptor tyrosine kinase that inhibits the GTPase activity of p21cdc42. Nature 363: 364–367.E. ManserT. LeungH. SalihuddinL. TanL. Lim1993A non-receptor tyrosine kinase that inhibits the GTPase activity of p21cdc42.Nature363364367

34. Williams SJ, Wreschner DH, Tran M, Eyre HJ, Sutherland GR, et al. (2001) MUC13, a novel human cell surface mucin expressed by epithelial and hemopoietic cells. J Biol Chem 276: 18327–18336.SJ WilliamsDH WreschnerM. TranHJ EyreGR Sutherland2001MUC13, a novel human cell surface mucin expressed by epithelial and hemopoietic cells.J Biol Chem2761832718336

35. Wang Y, Kobori JA, Hood L (1993) The ht beta gene encodes a novel CACCC box-binding protein that regulates T-cell receptor gene expression. Mol Cell Biol 13: 5691–5701.Y. WangJA KoboriL. Hood1993The ht beta gene encodes a novel CACCC box-binding protein that regulates T-cell receptor gene expression.Mol Cell Biol1356915701

36. Lang TG, Hansson GC, Samuelsson T (2006) An inventory of mucin genes in the chicken genome shows that the mucin domain of Muc13 is encoded by multiple exons and that ovomucin is part of a locus of related gel-forming mucins. BMC Genomics 7: -.TG LangGC HanssonT. Samuelsson2006An inventory of mucin genes in the chicken genome shows that the mucin domain of Muc13 is encoded by multiple exons and that ovomucin is part of a locus of related gel-forming mucins.BMC Genomics7-

37. Kim M, Ogawa M, Fujita Y, Yoshikawa Y, Nagai T, et al. (2009) Bacteria hijack integrin-linked kinase to stabilize focal adhesions and block cell detachment. Nature 459: 578-U109.M. KimM. OgawaY. FujitaY. YoshikawaT. Nagai2009Bacteria hijack integrin-linked kinase to stabilize focal adhesions and block cell detachment.Nature459578-U109

38. Finnemann SC, Bonilha VL, Marmorstein AD, Rodriguez-Boulan E (1997) Phagocytosis of rod outer segments by retinal pigment epithelial cells requires alpha(v)beta5 integrin for binding but not for internalization. Proc Natl Acad Sci U S A 94: 12932–12937.SC FinnemannVL BonilhaAD MarmorsteinE. Rodriguez-Boulan1997Phagocytosis of rod outer segments by retinal pigment epithelial cells requires alpha(v)beta5 integrin for binding but not for internalization.Proc Natl Acad Sci U S A941293212937

39. Memmo LM, McKeown-Longo P (1998) The alpha v beta 5 integrin functions as an endocytic receptor for vitronectin. J Cell Sci 111: 425–433.LM MemmoP. McKeown-Longo1998The alpha v beta 5 integrin functions as an endocytic receptor for vitronectin.J Cell Sci111425433

40. Mousa SA, Lorelli W, Campochiaro PA (1999) Role of hypoxia and extracellular matrix-integrin binding in the modulation of angiogenic growth factors secretion by retinal pigmented epithelial cells. J Cell Biochem 74: 135–143.SA MousaW. LorelliPA Campochiaro1999Role of hypoxia and extracellular matrix-integrin binding in the modulation of angiogenic growth factors secretion by retinal pigmented epithelial cells.J Cell Biochem74135143

41. Singh B, Su YC, Riesbeck K (2010) Vitronectin in bacterial pathogenesis: a host protein used in complement escape and cellular invasion. Mol Microbiol 78: 545–560.B. SinghYC SuK. Riesbeck2010Vitronectin in bacterial pathogenesis: a host protein used in complement escape and cellular invasion.Mol Microbiol78545560

42. Chhatwal GS, Preissner KT, Muller-Berghaus G, Blobel H (1987) Specific binding of the human S protein (vitronectin) to streptococci, Staphylococcus aureus, and Escherichia coli. Infect Immun 55: 1878–1883.GS ChhatwalKT PreissnerG. Muller-BerghausH. Blobel1987Specific binding of the human S protein (vitronectin) to streptococci, Staphylococcus aureus, and Escherichia coli.Infect Immun5518781883

43. Kasuga I, Hogg JC, Pare PD, Hayashi S, Sedgwick EG, et al. (2009) Role of genetic susceptibility to latent adenoviral infection and decreased lung function. Resp Med 103: 1672–1680.I. KasugaJC HoggPD PareS. HayashiEG Sedgwick2009Role of genetic susceptibility to latent adenoviral infection and decreased lung function.Resp Med10316721680

44. Aplin JD, Spanswick C, Behzad F, Kimber SJ, Vicovac L (1996) Integrins beta 5, beta 3 and alpha v are apically distributed in endometrial epithelium. Mol Hum Reprod 2: 527–534.JD AplinC. SpanswickF. BehzadSJ KimberL. Vicovac1996Integrins beta 5, beta 3 and alpha v are apically distributed in endometrial epithelium.Mol Hum Reprod2527534

45. Price AL, Zaitlen NA, Reich D, Patterson N (2010) New approaches to population stratification in genome-wide association studies. Nat Rev Genet 11: 459–463.AL PriceNA ZaitlenD. ReichN. Patterson2010New approaches to population stratification in genome-wide association studies.Nat Rev Genet11459463

46. Wu C, DeWan A, Hoh J, Wang Z (2011) A comparison of association methods correcting for population stratification in case-control studies. Ann Hum Genet 75: 418–427.C. WuA. DeWanJ. HohZ. Wang2011A comparison of association methods correcting for population stratification in case-control studies.Ann Hum Genet75418427