Open Access

A genome-wide screen identifies a single β-defensin gene cluster in the chicken: implications for the origin and evolution of mammalian defensins

  • Yanjing Xiao1,
  • Austin L Hughes2,
  • Junko Ando3,
  • Yoichi Matsuda3,
  • Jan-Fang Cheng4,
  • Donald Skinner-Noble1 and
  • Guolong Zhang1Email author
BMC Genomics20045:56

DOI: 10.1186/1471-2164-5-56

Received: 13 May 2004

Accepted: 13 August 2004

Published: 13 August 2004

Abstract

Background

Defensins comprise a large family of cationic antimicrobial peptides that are characterized by the presence of a conserved cysteine-rich defensin motif. Based on the spacing pattern of cysteines, these defensins are broadly divided into five groups, namely plant, invertebrate, α-, β-, and θ-defensins, with the last three groups being mostly found in mammalian species. However, the evolutionary relationships among these five groups of defensins remain controversial.

Results

Following a comprehensive screen, here we report that the chicken genome encodes a total of 13 different β-defensins but with no other groups of defensins being discovered. These chicken β-defensin genes, designated as Gallinacin 1–13, are clustered densely within a 86-Kb distance on the chromosome 3q3.5-q3.7. The deduced peptides vary from 63 to 104 amino acid residues in length sharing the characteristic defensin motif. Based on the tissue expression pattern, 13 β-defensin genes can be divided into two subgroups with Gallinacin 1–7 being predominantly expressed in bone marrow and the respiratory tract and the remaining genes being restricted to liver and the urogenital tract. Comparative analysis of the defensin clusters among chicken, mouse, and human suggested that vertebrate defensins have evolved from a single β-defensin-like gene, which has undergone rapid duplication, diversification, and translocation in various vertebrate lineages during evolution.

Conclusions

We conclude that the chicken genome encodes only β-defensin sequences and that all mammalian defensins are evolved from a common β-defensin-like ancestor. The α-defensins arose from β-defensins by gene duplication, which may have occurred after the divergence of mammals from other vertebrates, and θ-defensins have arisen from α-defensins specific to the primate lineage. Further analysis of these defensins in different vertebrate lineages will shed light on the mechanisms of host defense and evolution of innate immunity.

Background

Defensins constitute a large family of small, cysteine-rich, cationic peptides that are capable of killing a broad spectrum of pathogens, including various bacteria, fungi, and certain enveloped viruses [15]. These peptides play a critical role in host defense and disease resistance by protecting the hosts against infections. Transgenic mice expressing human enteric defensin HD5 are fully protected against the doses of Salmonella typhimurium that are otherwise lethal to the wide-type mice [6]. Conversely, mice deficient in the matrilysin gene, which is responsible for activating enteric defensins, become more susceptible to oral infection with S. typhimurium [7].

Defensins have been identified in species ranging from plants, insects to animals and humans [15]. Characterized by the presence of 6–8 cysteine residues in relatively defined positions, all defensins are structurally related in that they form 3–4 intramolecular disulfide bonds and 2–3 antiparallel β-sheets with or without an α-helix. Based on the spacing pattern of cysteines, these peptides are broadly divided into five groups; namely plant, invertebrate, α-, β-, and θ-defensins [15]. Alignment of all known defensin sequences revealed the consensus defensin motif of each group as follows: plant defensin: C-X8–11-C-X3–5-C-X3-C-X9–12-C-X4–11-C-X1-C-X3-C; invertebrate defensin: C-X5–16-C-X3-C-X9–10-C-X4–7-C-X1-C; α-defensin: C-X1-C-X3–4-C-X9-C-X6–10-C-C; and β-defensin: C-X4–8-C-X3–5-C-X9–13-C-X4–7-C-C. The α- and β-defensins are unique to vertebrate animals with α-defensins only being found in rodents and primates, while β-defensins are present in all mammalian species investigated [13]. On the other hand, θ-defensins have only been found in certain primates as a result of posttranslational ligation of two α-defensin-like sequences [810]. A pseudogene for θ-defensin is also present in humans [11].

Analysis of human and mouse genomes indicated that β-defensins form 4–5 distinct clusters on different chromosomes with each cluster consisting of multiple defensin genes [12]. Interestingly, the single mammalian α-defensin locus is located on a β-defensin cluster with θ-defensins residing in the center of α-defensins [12]. Studies with mammalian defensins suggested a rapid duplication followed by positive selection and diversification within each group [1318]. However, the evolutionary relationships among three groups of mammalian defensins and among plant, invertebrate, and mammalian defensins remain controversial. Similarity in spatial structure and biological functions favors the notion that all mammalian defensins are evolutionarily related [19], although a phylogenetic analysis suggested a closer relationship between β- and insect defensins than between α- and β-defensins [16].

Existence of a large number of expressed sequence tag (EST) sequences and recent completion of chicken genome sequencing at a 6.6× coverage [20] provided a timely opportunity to discover a complete repertoire of defensin-related sequences in birds for studying the evolutionary relationship between invertebrate and mammalian defensins. Here we report identification of a single β-defensin cluster that is composed of 13 genes located on the chicken chromosome 3q3.5-q3.7. Evolutionary and comparative analyses of these chicken β-defensins with mammalian homologues strongly suggested that all mammalian defensins have evolved from a common β-defensin-like ancestor, which has undergone rapid duplication, positive diversifying selection, and chromosomal translocations, thereby giving rising to multiple gene clusters on different chromosomal regions.

Results and Discussion

Discovery of novel chicken defensins

To identify novel defensin genes in the chicken, all five groups of known defensin-like peptide sequences from plants, invertebrates, and vertebrates were first queried individually against the translated chicken nonredundant (NR), EST, high throughput genomic sequence (HTGS), and whole-genome shortgun sequence (WGS) databases in the GenBank by using the TBLASTN program[21]. All potential hits were then examined manually for the presence of the characteristic cysteine motifs. For every novel defensin identified, additional iterative BLAST searches were performed until no more novel sequences could be found. In addition to three known chicken β-defensins (Gal 1–3) [22, 23], nine novel putative sequences, namely Gal 4–12, have been found in the EST database with at least two hits for each, and such sequences have also been confirmed in genomic sequences (Table 1). Because of the fact that mammalian defensins tend to form clusters [12, 14, 15, 18], all chicken HTGS and WGS sequences containing defensin sequences were also retrieved from GenBank, translated into six open reading frames, and manually curated. As a result, an additional putative β-defensin, Gal13, was identified in several genomic clones (Table 1). The open reading frame of Gal13 was predicted by GENSCAN [24] and confirmed by directly sequencing of RT-PCR product amplified from chicken kidney.
Table 1

Identification of chicken β-defensins

Gene

EST1,2

HTGS

WGS

Gene Size (bp)3

    

E 1

I 1

E 2

I 2

E 3

I 3

E 4

Gal1

BX260462

AC110874

AADN01058097

70

972

88

482

127

496

217

Gal2

BX540940

AC110874

AADN01058097

66

1113

143

183

121

674

204

Gal3

 

AC110874

AADN01058097

AADN01058096

  

53

980

109

1180

215

Gal4

BU451960

 

AADN01058096

  

136

461

127

117

141

Gal5

BU389548

 

AADN01058096

  

290

445

127

355

187

Gal6

CF251501

AC110874

AADN01058097

AADN01058098

52

704

86

705

130

249

234

Gal7

CF251115

AC110874

AADN01058098

50

656

86

201

130

234

248

Gal8

BU242665

AC110874

AADN01058098

71

915

91

259

134

706

494

Gal9

BX270804

AC110874

AADN01058098

  

220

1592

130

781

343

Gal10

AW198592

AC110874

AADN01058099

  

118

268

133

1719

381

Gal11

BM440069

AC110874

AADN01058101

  

63

966

129

1001

460

Gal12

BX257296

AC110874

AADN01058102

  

84

396

420

  

Gal13

 

AC110874

AADN01058102

  

61

1016

118

3322

91

1 Abbreviations: EST, expressed sequence tag; HTGS, high throughput genomic sequence; WGS, whole-genome shortgun sequence; E, exon;I, intron. 2 One EST sequence entry is given only for the exemplary purpose. In each case, more than two independent EST sequences have been found, except for Gal3 and Gal13, both of which have no EST sequences. Gal3 was found through homology cloning [23], and Gal13 was predicted by us from the genomic sequence. 3 All Gal genes are predicted to consist of four exons separated by three introns, except for Gal12, whose last two exons are fused together. The absence of additional sequence information at the 5'-untranslated regions of the cDNA sequences prevented prediction of the sizes of first exon and intron for Gal 3–5 and Gal 9–13 genes.

No other sequence containing β-defensin-like six-cysteine motif has been found in NR, EST or genomic databases, suggesting that 13 Gal genes constitute the entire repertoire of the β-defensin family encoded in the chicken genome. Although it is highly unlikely, we could not rule out the possibility that additional defensin-related genes with distant homology might be uncovered in the chicken by different computational search methods such as the use of Hidden Markov models [12, 15]. It is noted that none of other groups of defensins have been discovered in the chicken, indicating that plant, invertebrate, α-, and θ-defensins are absent in the chicken lineage.

Similar to Gal 1–3, 10 novel β-defensins, deduced from either EST or genomic sequences, vary from 63 to 104 amino acid residues in length. Alignment of these peptides revealed a conservation of the signal sequence at the N-terminus and the characteristic six-cysteine defensin motif at the C-terminus (Figure 1). Consistent with the fact that all β-defensins are a group of secreted molecules in response to infections, the signal sequences of all chicken defensins are hydrophobic and rich in leucines. In addition, the mature C-terminal sequences are all positively charged due to the presence of excess arginines and lysines. Interestingly, Gal11 contains two tandem, but highly divergent, copies of the six-cysteine motif at the C-terminus, and is the only defensin having such sequences. Functional significance for existence of such two defensin motifs remains to be studied.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-56/MediaObjects/12864_2004_Article_158_Fig1_HTML.jpg
Figure 1

Multiple sequence alignment of chicken β-defensins. The intervening region between signal and mature peptide sequence is the short propiece. The conserved residues are shaded. Also shown is the length of each peptide. Notice the six-cysteine defensin motif is highly conserved. The six cysteines in the second tandem copy of the defensin motif in Gal11 are boxed.

Evolutionary analysis of vertebrate β-defensins

Phylogenetic analysis of vertebrate β-defensins showed that chicken defensins clustered with various different groups of mammalian β-defensins (Figure 2). However, the bootstrap support for these patterns was very weak (less than 50% in all cases). The clustering of certain chicken β-defensins with mammalian homologues suggests that major subfamilies of β-defensins arose before the last common ancestor of birds and mammals, estimated to have occurred about 310 million years ago [25]. This in turn implies that some duplication of β-defensin genes must have taken place before the divergence of birds and mammals. The apparent lack of α-defensins in the chicken and other non-mammalian species (G. Zhang, unpublished data) suggests that α-defensins may have evolved after mammals diverged from other vertebrates.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-56/MediaObjects/12864_2004_Article_158_Fig2_HTML.jpg
Figure 2

Phylogenetic relationship of vertebrate β-defensins. The tree was constructed by the neighbor-joining method and the reliability of each branch was assessed by using 1000 bootstrap replications. Numbers on the branches indicate the percentage of 1000 bootstrap samples supporting the branch. Only branches supported by a bootstrap value of at least 50% are indicated. Chicken β-defensins are highlighted in yellow. Abbreviations: BNBD, bovine neutrophil β-defensin; LAP, lingual antimicrobial peptide; EBD, enteric β-defensin; TAP, tracheal antimicrobial peptide; PBD, porcine β-defensin; DEFB/Defb, β-defensin; Gal, Gallinacin; GAPDH, glyceraldehyde-3-phosphate dehydrogenase.

Comparison of the numbers of synonymous and nonsynonymous nucleotide substitutions provides a powerful test of the hypothesis that positive Darwinian selection has acted to favor changes at the amino acid level [26]. This approach has previously been applied to both α- and β-defensins of mammals and has revealed positive selection acting on the mature defensin but not on other regions of the gene [16, 17]. In the comparison of the chicken β-defensin sequences, synonymous sites were saturated with changes or nearly so, making it impossible to test the hypothesis of positive selection in every case. In pairwise comparisons among all sequences, mean pS in the propeptide region was 0.551 ± 0.036 (S.E.), while mean pN was 0.369 ± 0.040. In the mature defensin region, mean pS was 0.673 ± 0.027, while mean pN was 0.534 ± 0.051. Mean pN in the mature defensin was significantly greater than that in the propeptide (z-test; P < 0.05), indicating lesser functional constraint on the amino acid sequence of the former. The high mean pS shows that chicken β-defensin genes have not duplicated recently, unlike β-defensin genes of the bovine [16]. In the comparison between the most closely related pair of sequences (Gal6 and Gal7), mean pS in the mature defensin was 0.221 ± 0.082, while mean pN was 0.331 ± 0.076. While these values are not significantly different at the 5% level, the fact that pN was higher than pS suggested that positive selection may have acted to diversify the mature defensin region between these two genes.

Genomic organization and chromosomal localization of the chicken β-defensin gene cluster

Searching through HTGS database led to identification of two overlapping bacterial artificial chromosome (BAC) sequences, TAM31-54I5 (accession no. AC110874) and CH261-162O9 (accession no. AC146292), both of which were sequenced and deposited earlier by one of us (J.F. Chen). Alignment of these two sequences allowed to re-order three DNA fragments in AC110874 and to construct a continuous, gap-free genomic contig that includes 11 Gal genes except for Gal4 and Gal5. Later search of chicken WGS sequences released on February 29, 2004 confirmed the order of the genomic contig that we assembled and also revealed the locations of two remaining genes, Gal4 and Gal5, both of which reside on a WGS (accession no. AADN01058096) that overlaps with AC110874 (Figure 3). The position and orientation of each Gal gene were obtained by comparing its cDNA with the assembled DNA sequence. As shown in Figure 3, all 13 Gal genes were clustered densely within a distance of 86.0 Kb on the genome. It was also confirmed by aligning such a contig with the chicken genome assembly, in which 13 Gal genes are located on six WGS contigs (Table 1) of chromosome 3 that are only ~3.3 Mb from the distal end. Consistent with this, the Gal gene cluster was physically mapped to the tip of chicken chromosome 3 at the region of q3.5-q3.7 by fluorescence in situ hybridization (FISH) using the TAM31-54I5 BAC DNA as probe (Figure 4).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-56/MediaObjects/12864_2004_Article_158_Fig3_HTML.jpg
Figure 3

Genomic organization of the chicken β-defensin gene cluster. The horizontal lines at the bottom represent the three overlapping genomic clones that were used to assemble the continuous, gap-free contig. The position of each gene is represented by a solid vertical bar and the width of each bar is proportional to the size of each gene. The direction of transcription is indicated by the triangle above each gene. The genes with solid triangles are transcribed in the direction opposite to the ones with open triangles. Slanted lines refer to the sequences omitted. Note that the three fragments of AC110874 sequence have been re-ordered and the gaps have been filled following alignment with AC146292.

https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-56/MediaObjects/12864_2004_Article_158_Fig4_HTML.jpg
Figure 4

Chromosomal localization of the chicken β-defensin gene cluster by fluorescence in situ hybridization. The BAC clone TAM31-54I4, which harbors 11 Gal genes, was mapped to chicken chromosome 3q3.5-q3.7. Arrows indicate the hybridization signals.

Comparing the cDNA with genomic sequences also revealed the structure of each Gal gene. Unlike most mammalian β-defensin genes, which primarily consist of two exons and one intron, the Gal genes were found to be composed of four short exons separated by three introns with variable lengths ranging from 117 bp to 3,322 bp (Table 1). Gal12 is an exception, in which the last two exons have been fused together. While the first exon of the Gal genes encodes 5'-untranslated region (UTR) and the majority of the last exon encodes 3'-UTR as well as a few C-terminal amino acids, two internal exons resemble mammalian β-defensin genes in that one exon encodes the signal and pro-sequence and the other encodes the mature sequence with six-cysteine motif [19, 2729]. Apparently, the first two and the last two exons of the Gal genes have joined together during the evolution as a result of exon shuffling, which occurred in many other evolutionarily conserved gene families [30], including invertebrate defensins [5]. The fusion of defensin exons in mammals is presumably adaptive because it allows a faster mobilization of such host defense molecules to better cope with invading microbes.

Tissue expression patterns of chicken β-defensins

It has been shown that Gal1 and Gal2 are expressed in bone marrow and lung, while Gal3 is more preferentially expressed in bone marrow, tongue, trachea, and bursa of Fabricius [23]. To study the tissue expression patterns of novel Gal genes that we identified, RT-PCR was performed with a panel of 32 different chicken tissues. Similar to Gal 1–3, Gal 4–7 are highly restricted to bone marrow cells with Gal5 also expressed in tongue, trachea, lung, and brain at lower levels (Figure 5). By contrast, the six remaining genes, Gal 8–13, were not found in bone marrow, but instead in liver, kidney, testicle, ovary, and male and female reproductive tracts (Figure 5). These results clearly suggested that all chicken β-defensin genes can be divided into two subgroups. Seven genes (Gal 1–7) are predominantly expressed in bone marrow and the respiratory tract, whereas the other six genes (Gal 8–13) are more restricted to liver and the urogenital tract. However, the functional significance and transcriptional regulatory mechanisms of these genes during inflammation and infection remain to be investigated.
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-56/MediaObjects/12864_2004_Article_158_Fig5_HTML.jpg
Figure 5

Tissue expression patterns of 10 novel chicken β-defensins by RT-PCR. See Materials and Methods for details. The number of PCR cycles was optimized for each gene, and the specificity of each PCR product was confirmed by sequencing. The house-keeping gene, GAPDH, was used for normalization of the template input.

Comparative analysis of chicken and mammalian β-defensin gene clusters

To study the origin and evolution of mammalian defensins, a comparative analysis of β-defensin gene clusters in the chicken, mouse, and human was performed by employing additional, more phylogenetically conserved gene markers surrounding the defensin clusters. As shown in Figure 6, two genes, CTSB (Cathepsin B, accession no. NP_680093) and a human EST sequence (accession no. BE072524) immediately located centromeric to chicken defensins, were also found to be conserved in the defensin gene clusters on human chromosome 8p22 and mouse chromosome 14C3. Similarly, another gene, HARL2754 (accession no. XP_372011) that is 6-Kb telemetric to Gal4 is also conserved in another defensin cluster in human (8p23) or mouse (8A1.3) (Figure 6).
https://static-content.springer.com/image/art%3A10.1186%2F1471-2164-5-56/MediaObjects/12864_2004_Article_158_Fig6_HTML.jpg
Figure 6

Comparative analysis of defensin clusters among the chicken, mouse, and human. The gene clusters were drawn proportionally according to their sizes. Each vertical line/bar represents the position of a gene, and the width of each line/bar is proportional to the size of each gene. Three highly conserved genes (CTSB, BE072524, and HARL2754) surrounding the defensin clusters in the chicken, mouse, and human were connected by solid lines. The position of the α-defensin locus (DEFA) was indicated as an open square. Note that the human θ-defensin pseudogene resides in the DEFA locus. The positions and orders of defensin genes in human and mouse were drawn based on the genome assemblies released in July 2003 and October 2003, respectively. Abbreviations: GGA, chicken chromosome; MMA, mouse chromosome; HSA, human chromosome; Tel, telomere; Cen, centromere.

These results strongly suggested that all vertebrate β-defensins are evolved from a single gene. This conclusion is further supported by the fact that there are three highly similar β-defensin-like sequences present in the largely finished zebrafish genome (G. Zhang, unpublished data). In addition, a group of homologous β-defensin-like sequences, namely crotamine and myotoxins, have been found in several Crotalus snakes [31], which are presumably derived from a single ancestral gene. The appearance of multiple β-defensin gene clusters on different chromosomal regions in mammalian species [12] is apparently a result of rapid gene duplication, positive diversifying selection, and chromosomal translocation following divergence of mammals from other vertebrate lineages.

In addition to the structural conservation between β-defensin-like sequences in the rattlesnake and mammals [32], a growing body of evidence suggests that their functions appear to be largely conserved in that both are capable of interacting negatively-charged lipid membranes followed by formation of ion channels or pores [3234]. It is noteworthy that the conservation of Cathepsin B (CTSB) adjacent to β-defensins is perhaps not surprising, given the recent finding that cathepsins are involved in the cleavage and inactivation of β-defensins [35].

Conclusions

We have showed that chicken genome encodes a total of 13 different β-defensin genes clustered densely within a 86-Kb distance on the chromosome 3q3.5-q3.7, but with no α-defensin genes. These peptides exhibit homology to different subgroups of mammalian β-defensins-, consistent with the hypothesis that α-defensins and β-defensins arose by gene duplication after the divergence of birds and mammals. The θ-defensins are specific to primates; and thus appear to have arisen from α-defensins by gene duplication specific to the primate lineage. Apparently, the evolution of defensins is rapid and driven by duplication and positive diversifying selection. Collectively, this study represents the first large-scale detailed investigation of defensins in non-mammalian vertebrates. There is no doubt that further analysis of these defensin genes will lead to a better understanding of host defense mechanisms and evolution of innate immunity.

Methods

Computational search for novel chicken defensins

To identify novel defensins in the chicken, all known cysteine-containing defensin-like peptide sequences discovered in plants, invertebrates, birds, and mammals were individually queried against the translated chicken NR, EST, HTGS, and WGS databases in the GenBank by using the TBLASTN program [21] with default settings on the NCBI web site [36]. All potential hits were then examined for the presence of the characteristic defensin motif. For every novel defensin identified, additional iterative BLAST searches were performed until no more novel sequences could be revealed. Because mammalian defensins tend to form clusters [12, 14, 15, 18], all chicken genomic sequences containing defensin sequences were also retrieved from the GenBank and translated into six open reading frames and curated manually for the presence of the defensin motif in order to discover potential sequences with distant homology.

Alignment and phylogenetic analysis of chicken β-defensins

Multiple sequence alignment was constructed by using the ClustalW program (version 1.82) [37]. A phylogenetic tree of amino acid sequences of mature β-defensins was constructed by the neighbor-joining method [38]. So that a comparable data set would be used for all pairwise comparisons, any site at which the alignment postulated a gap in any sequence was excluded from the analysis. To maximize the number of sites available for analysis, certain sequences with large deletions were excluded from the analysis. Because the sequences were very short (25 aligned sites), no correction for multiple hits was applied. The reliability of clustering patterns within the tree was assessed by bootstrapping; 1000 bootstrap pseudo-samples were used. The proportion of synonymous nucleotide differences per synonymous site (pS) and the proportion of nonsynonymous nucleotide differences per nonsynonymous site (pN) were estimated by the method of Nei and Gojobori [26]. Again, no correction for multiple hits was applied because a small number of sites were examined.

Assembly of the chicken β-defensin gene cluster

To generate a continuous defensin gene cluster, the HTGS and WGS sequences containing the putative defensin genes were retrieved from the GenBank, aligned to generate a longer contig, which was confirmed later by searching through the assembled chicken genome released on February 29, 2004, by using the BLAT program [39] under the UCSC Genome Browser web site [40]. The relative positions, orientations, and structural organizations of individual genes were determined by comparing its cDNA sequence to the continuous genomic contig that we assembled.

Chromosome localization of the chicken β-defensin gene cluster

Fluorescence in situ hybridization (FISH) was used for chromosomal assignment of the chicken β-defensin gene cluster by using the BAC clone TAM31-54I4 as probe, which harbors 11 Gal genes. Metaphase chromosome speads were prepared from mitogen-stimulated chicken splenocyte culture as we described [41, 42]. The BAC clone was labeled by nick translation with biotin 16-dUTP (Roche Diagnostics), hybridized to metaphase chromosome DNA, followed by detection with FITC-labeled avidin (Roche Diagnostics) and staining with propidium iodide to simultaneously induce the R-banding.

RT-PCR analysis of the tissue expression patterns of chicken β-defensins

Total RNA was extracted with Trizol (Invitrogen) from a total of 32 different tissues from healthy, 2-month-old chickens (see Figure 5). A total of 4 μg RNA from each tissue were reverse transcribed with random hexamers and Superscript II reverse transcriptase by using a first-strand cDNA synthesis kit (Invitrogen) according to the instructions. The subsequent PCR was carried out with 1/40 of the first-strand cDNA and gene-specific primers for each β-defensin and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as described [28, 43]. Every pair of primers were designed to locate on different exons to aid in distinguishing PCR products amplified from cDNA vs. genomic DNA (Table 2). The PCR program used was: 94°C denaturation for 2 min, followed by different cycles of 94°C denaturation for 20 sec, 55°C annealing for 20 sec, and 72°C extension for 40 sec, followed by a final extension at 72°C for 5 min. The number of PCR cycle was optimized for each gene to ensure linear amplification (Table 2). A half of the PCR products were analyzed by electrophoresis on 1.2% agarose gels containing 0.5 μg/ml ethidium bromide. The specificity of each PCR product was confirmed by cloning of the PCR product into T/A cloning vector, followed by sequencing of the recombinant plasmid.
Table 2

Primer sequences used for RT-PCR analysis of novel chicken β-defensins

Gene

Primer Sequence

Product Size (bp)

Cycles Used

 

Sense

Antisense

cDNA

Genomic

 

Gal4

CATCTCAGTGTCGTTTCTCTGC

ACAATGGTTCCCCAAATCCAAC

321

899

36

Gal5

CTGCCAGCAAGAAAGGAACCTG

TGAACGTGAAGGGACATCAGAG

300

1100

36

Gal6

AGGATTTCACATCCCAGCCGTG

CAGGAGAAGCCAGTGAGTCATC

249

1203

36

Gal7

CTGCTGTCTGTCCTCTTTGTGG

CATTTGGTAGATGCAGGAAGGA

230

665

35

Gal8

ACAGTGTGAGCAGGCAGGAGGGA

CTCTTCTGTTCAGCCTTTGGTG

261

967

35

Gal9

GCAAAGGCTATTCCACAGCAG

AGCATTTCAGCTTCCCACCAC

211

1802

33

Gal10

TGGGGCACGCAGTCCACAAC

ATCAGCTCCTCAAGGCAGTG

298

2285

33

Gal11

ACTGCATCCGTTCCAAAGTCTG

TCGGGCAGCTTCTCTACAAC

301

1299

33

Gal12

CCCAGCAGGACCAAAGCAATG

GTGAATCCACAGCCAATGAGAG

335

731

36

Gal13

CATCGTTGTCATTCTCCTCCTC

ACTTGCAGCGTGTGGGAGTTG

175

4514

50

GAPDH

GCACGCCATCACTATCTTCC

CATCCACCGTCTTCTGTGTG

356

876

30

Note added in proof

Following submission of this manuscript, Lynn et al. reported independently discovery of seven novel chicken β-defensins in the chicken EST database by using homology search strategies [44]. Consistent with our conclusion, they also revealed occurrence of positive selection particularly in the mature region of chicken β-defensins following evolutionary analysis. Moreover, albeit the use of a different nomenclature, they confirmed that the expressions of Gal 4–7 are primarily in bone marrow, while other genes are more restricted to liver and the genitourinary tract.

List of abbreviations

Abbreviations: 

Gal, Gallinacin

NR: 

nonredundant

EST: 

expressed sequence tag

HTGS: 

high throughput genomic sequence

WGS: 

whole-genome shortgun sequence

BAC: 

bacterial artificial chromosome

FISH: 

fluorescence in situ hybridization

UTR: 

untranslated region

GAPDH: 

glyceraldehyde-3-phosphate dehydrogenase.

Declarations

Acknowledgements

This work was supported in part by Oklahoma Center for Advancement of Science and Technology Grant HR-136 (to GZ) and the Oklahoma Agricultural Experiment Station.

Authors’ Affiliations

(1)
Department of Animal Science, Oklahoma State University
(2)
Department of Biological Sciences, University of South Carolina
(3)
Center for Advanced Science and Technology, Hokkaido University
(4)
Department of Genome Sciences, Lawrence Berkeley National Laboratory

References

  1. Ganz T: Defensins: antimicrobial peptides of innate immunity. Nat Rev Immunol. 2003, 3: 710-720. 10.1038/nri1180.View ArticlePubMedGoogle Scholar
  2. Lehrer RI, Ganz T: Defensins of vertebrate animals. Curr Opin Immunol. 2002, 14: 96-102. 10.1016/S0952-7915(01)00303-X.View ArticlePubMedGoogle Scholar
  3. Schutte BC, McCray P.B.,Jr.: Beta-defensins in lung host defense. Annu Rev Physiol. 2002, 64: 709-748. 10.1146/annurev.physiol.64.081501.134340.View ArticlePubMedGoogle Scholar
  4. Thomma BP, Cammue BP, Thevissen K: Plant defensins. Planta. 2002, 216: 193-202. 10.1007/s00425-002-0902-6.View ArticlePubMedGoogle Scholar
  5. Froy O, Gurevitz M: Arthropod and mollusk defensins--evolution by exon-shuffling. Trends Genet. 2003, 19: 684-687. 10.1016/j.tig.2003.10.010.View ArticlePubMedGoogle Scholar
  6. Salzman NH, Ghosh D, Huttner KM, Paterson Y, Bevins CL: Protection against enteric salmonellosis in transgenic mice expressing a human intestinal defensin. Nature. 2003, 422: 522-526. 10.1038/nature01520.View ArticlePubMedGoogle Scholar
  7. Wilson CL, Ouellette AJ, Satchell DP, Ayabe T, Lopez-Boado YS, Stratman JL, Hultgren SJ, Matrisian LM, Parks WC: Regulation of intestinal alpha-defensin activation by the metalloproteinase matrilysin in innate host defense. Science. 1999, 286: 113-117. 10.1126/science.286.5437.113.View ArticlePubMedGoogle Scholar
  8. Tang YQ, Yuan J, Miller CJ, Selsted ME: Isolation, characterization, cDNA cloning, and antimicrobial properties of two distinct subfamilies of alpha-defensins from rhesus macaque leukocytes. Infect Immun. 1999, 67: 6139-6144.PubMed CentralPubMedGoogle Scholar
  9. Tang YQ, Yuan J, Osapay G, Osapay K, Tran D, Miller CJ, Ouellette AJ, Selsted ME: A cyclic antimicrobial peptide produced in primate leukocytes by the ligation of two truncated alpha-defensins. Science. 1999, 286: 498-502. 10.1126/science.286.5439.498.View ArticlePubMedGoogle Scholar
  10. Nguyen TX, Cole AM, Lehrer RI: Evolution of primate theta-defensins: a serpentine path to a sweet tooth. Peptides. 2003, 24: 1647-1654. 10.1016/j.peptides.2003.07.023.View ArticlePubMedGoogle Scholar
  11. Cole AM, Hong T, Boo LM, Nguyen T, Zhao C, Bristol G, Zack JA, Waring AJ, Yang OO, Lehrer RI: Retrocyclin: a primate peptide that protects cells from infection by T- and M-tropic strains of HIV-1. Proc Natl Acad Sci USA. 2002, 99: 1813-1818. 10.1073/pnas.052706399.PubMed CentralView ArticlePubMedGoogle Scholar
  12. Schutte BC, Mitros JP, Bartlett JA, Walters JD, Jia HP, Welsh MJ, Casavant TL, McCray P.B.,Jr.: Discovery of five conserved beta -defensin gene clusters using a computational search strategy. Proc Natl Acad Sci USA. 2002, 99: 2129-2133. 10.1073/pnas.042692699.PubMed CentralView ArticlePubMedGoogle Scholar
  13. Antcheva N, Boniotto M, Zelezetsky I, Pacor S, Falzacappa MV, Crovella S, Tossi A: Effects of positively selected sequence variations in human and Macaca fascicularis beta-defensins 2 on antimicrobial activity. Antimicrob Agents Chemother. 2004, 48: 685-688. 10.1128/AAC.48.2.685-688.2004.PubMed CentralView ArticlePubMedGoogle Scholar
  14. Maxwell AI, Morrison GM, Dorin JR: Rapid sequence divergence in mammalian beta-defensins by adaptive evolution. Mol Immunol. 2003, 40: 413-421. 10.1016/S0161-5890(03)00160-3.View ArticlePubMedGoogle Scholar
  15. Semple CA, Rolfe M, Dorin JR: Duplication and selection in the evolution of primate beta-defensin genes. Genome Biol. 2003, 4: R31-10.1186/gb-2003-4-5-r31.PubMed CentralView ArticlePubMedGoogle Scholar
  16. Hughes AL: Evolutionary diversification of the mammalian defensins. Cell Mol Life Sci. 1999, 56: 94-103. 10.1007/s000180050010.View ArticlePubMedGoogle Scholar
  17. Hughes AL, Yeager M: Coordinated amino acid changes in the evolution of mammalian defensins. J Mol Evol. 1997, 44: 675-682.View ArticlePubMedGoogle Scholar
  18. Morrison GM, Semple CA, Kilanowski FM, Hill RE, Dorin JR: Signal sequence conservation and mature peptide divergence within subgroups of the murine beta-defensin gene family. Mol Biol Evol. 2003, 20: 460-470. 10.1093/molbev/msg060.View ArticlePubMedGoogle Scholar
  19. Liu L, Zhao C, Heng HH, Ganz T: The human beta-defensin-1 and alpha-defensins are encoded by adjacent genes: two peptide families with differing disulfide topology share a common ancestry. Genomics. 1997, 43: 316-320. 10.1006/geno.1997.4801.View ArticlePubMedGoogle Scholar
  20. Chicken Genome Sequencing Project at Washington University. [http://www.genome.wustl.edu/projects/chicken/]
  21. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. 10.1006/jmbi.1990.9999.View ArticlePubMedGoogle Scholar
  22. Harwig SS, Swiderek KM, Kokryakov VN, Tan L, Lee TD, Panyutich EA, Aleshina GM, Shamova OV, Lehrer RI: Gallinacins: cysteine-rich antimicrobial peptides of chicken leukocytes. FEBS Lett. 1994, 342: 281-285. 10.1016/0014-5793(94)80517-2.View ArticlePubMedGoogle Scholar
  23. Zhao C, Nguyen T, Liu L, Sacco RE, Brogden KA, Lehrer RI: Gallinacin-3, an inducible epithelial beta-defensin in the chicken. Infect Immun. 2001, 69: 2684-2691. 10.1128/IAI.69.4.2684-2691.2001.PubMed CentralView ArticlePubMedGoogle Scholar
  24. Burge C, Karlin S: Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997, 268: 78-94. 10.1006/jmbi.1997.0951.View ArticlePubMedGoogle Scholar
  25. Hughes AL, Nei M: Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection. Nature. 1988, 335: 167-170. 10.1038/335167a0.View ArticlePubMedGoogle Scholar
  26. Nei M, Gojobori T: Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol. 1986, 3: 418-426.PubMedGoogle Scholar
  27. Mallow EB, Harris A, Salzman N, Russell JP, DeBerardinis RJ, Ruchelli E, Bevins CL: Human enteric defensins. Gene structure and developmental expression. J Biol Chem. 1996, 271: 4038-4045. 10.1074/jbc.271.8.4038.View ArticlePubMedGoogle Scholar
  28. Zhang G, Hiraiwa H, Yasue H, Wu H, Ross CR, Troyer D, Blecha F: Cloning and characterization of the gene for a new epithelial beta- defensin. Genomic structure, chromosomal localization, and evidence for its constitutive expression. J Biol Chem. 1999, 274: 24031-24037. 10.1074/jbc.274.34.24031.View ArticlePubMedGoogle Scholar
  29. Ouellette AJ, Darmoul D, Tran D, Huttner KM, Yuan J, Selsted ME: Peptide localization and gene structure of cryptdin 4, a differentially expressed mouse paneth cell alpha-defensin. Infect Immun. 1999, 67: 6643-6651.PubMed CentralPubMedGoogle Scholar
  30. Patthy L: Genome evolution and the evolution of exon-shuffling--a review. Gene. 1999, 238: 103-114. 10.1016/S0378-1119(99)00228-0.View ArticlePubMedGoogle Scholar
  31. Bieber AL, Nedelkov D: Structural, biological and biochemical studies of myotoxin and homologous myotoxins. J Toxicol -Toxin Rev. 1997, 16: 33-52.View ArticleGoogle Scholar
  32. Nicastro G, Franzoni L, de Chiara C, Mancin AC, Giglio JR, Spisni A: Solution structure of crotamine, a Na+ channel affecting toxin from Crotalus durissus terrificus venom. Eur J Biochem. 2003, 270: 1969-1979.View ArticlePubMedGoogle Scholar
  33. Lencer WI, Cheung G, Strohmeier GR, Currie MG, Ouellette AJ, Selsted ME, Madara JL: Induction of epithelial chloride secretion by channel-forming cryptdins 2 and 3. Proc Natl Acad Sci U S A. 1997, 94: 8585-8589. 10.1073/pnas.94.16.8585.PubMed CentralView ArticlePubMedGoogle Scholar
  34. White SH, Wimley WC, Selsted ME: Structure, function, and membrane integration of defensins. Curr Opin Struct Biol. 1995, 5: 521-527. 10.1016/0959-440X(95)80038-7.View ArticlePubMedGoogle Scholar
  35. Taggart CC, Greene CM, Smith SG, Levine RL, McCray P. B., Jr., O'Neill S, McElvaney NG: Inactivation of human beta-defensins 2 and 3 by elastolytic cathepsins. J Immunol. 2003, 171: 931-937.View ArticlePubMedGoogle Scholar
  36. NCBI BLAST Programs. [http://www.ncbi.nlm.nih.gov/BLAST/]
  37. Higgins DG, Thompson JD, Gibson TJ: Using CLUSTAL for multiple sequence alignments. Methods Enzymol. 1996, 266: 383-402. 10.1016/S0076-6879(96)66024-8.View ArticlePubMedGoogle Scholar
  38. Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
  39. Kent WJ: BLAT--the BLAST-like alignment tool. Genome Res. 2002, 12: 656-664. 10.1101/gr.229202. Article published online before March 2002.PubMed CentralView ArticlePubMedGoogle Scholar
  40. UCSC Genome Browser. [http://genome.ucsc.edu/]
  41. Matsuda Y, Chapman VM: Application of fluorescence in situ hybridization in genome analysis of the mouse. Electrophoresis. 1995, 16: 261-272.View ArticlePubMedGoogle Scholar
  42. Suzuki T, Kurosaki T, Shimada K, Kansaku N, Kuhnlein U, Zadworny D, Agata K, Hashimoto A, Koide M, Koike M, Takata M, Kuroiwa A, Minai S, Namikawa T, Matsuda Y: Cytogenetic mapping of 31 functional genes on chicken chromosomes by direct R-banding FISH. Cytogenet Cell Genet. 1999, 87: 32-40. 10.1159/000015388.View ArticlePubMedGoogle Scholar
  43. Zhang G, Wu H, Shi J, Ganz T, Ross CR, Blecha F: Molecular cloning and tissue expression of porcine beta-defensin-1. FEBS Lett. 1998, 424: 37-40. 10.1016/S0014-5793(98)00134-3.View ArticlePubMedGoogle Scholar
  44. Lynn DJ, Higgs R, Gaines S, Tierney J, James T, Lloyd AT, Fares MA, Mulcahy G, O'Farrelly C: Bioinformatic discovery and initial characterisation of nine novel antimicrobial peptide genes in the chicken. Immunogenetics. 2004, 56: 170-177. 10.1007/s00251-004-0675-0.View ArticlePubMedGoogle Scholar

Copyright

© Xiao et al; licensee BioMed Central Ltd. 2004

This article is published under license to BioMed Central Ltd. This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.