Mapping quantitative trait loci in Gallus gallus using principal components

Mapeamento de locos de características quantitativas em Gallus gallus com utilização de componentes principais

BREEDING, GENETIC AND REPRODUCTION

Mapping quantitative trait loci in Gallus gallus using principal components

Mapeamento de locos de características quantitativas em Gallus gallus com utilização de componentes principais

Luís Fernando Batista Pinto^I; Irineu Umberto Packer^II,*; Mônica Corrêa Ledur^III; Ana Silvia Alves Meira Tavares Moura^IV; Kátia Nones^II; Luiz Lehmann Coutinho^II

^IDepartamento de Produção Animal, Escola de Medicina Veterinária, UFBA. Av. Ademar de Barros, 500, 40170-110, Salvador, BA, Brazil

^IIDepartamento de Zootecnia - ESALQ/USP. Av. Pádua Dias, 11, 13418-900, Piracicaba, SP, Brazil

^IIIEmbrapa Suínos e Aves. BR 153, Km 110, 89700-000, Concórdia, SC, Brazil

^IVDepartamento de Produção Animal, FMVZ/UNESP, 18618-000, Botucatu, SP, Brazil

ABSTRACT

This study aimed at mapping QTL (quantitative trait loci) using linear combinations of characteristics of economical interest in Gallus gallus. A total of 350 F2 chickens from an initial crossing among males from a broiler line (TT) with females from a layer line (CC) were used. It was conducted a QTL mapping in chromosomes of Gallus gallus (GGA1, GGA3, GGA5, GGA8, GGA11, and GGA13) for 20 performance and carcass traits. For detecting QTL, it was used the likelihood ratio test between a reduced model (including fixed effects of sex, hatch and random effect of infinitesimal genetic value) and a full model (including all the previous effects plus QTL effects). When original characterists were analyzed, that is, before the formation of linear combinations, six significant QTLs were mapped at 1% in the genome, four in the GGA1 (live weight at 35 days of age and at 42 days of age, abdominal fat and heart weight); and two on GGA3 (live weight at 35 and 42 days of age); three significant QTLs at 5% in the genome, one on GGA1 (head weight), one on GGA3 (wings weight), and one on GGA8 (gizzard weight); besides seven suggestive linkages for several traits. When QTLs were mapped for principal components, many mapped QTLs were confirmed for original traits, in addition to finding three QTLs and eight suggestive linkages not mapped for the original traits.

Key Words: broilers, canonical variables, chickens, genome

RESUMO

Este estudo teve por objetivo mapear QTL (quantitative trait loci) utilizando combinações lineares de características de interesse econômico em Gallus gallus. Foram estudadas 350 aves F₂ oriundas de um cruzamento inicial entre machos de uma linhagem de corte (TT) e fêmeas de uma linhagem de postura (CC). Foi conduzido mapeamento de QTL nos cromossomos Gallus gallus (GGA1, GGA3, GGA5, GGA8, GGA11 e GGA13), para 20 características de desempenho e de carcaça. Para detectar QTL foi utilizado o teste da razão de verossimilhanças entre um modelo reduzido (incluindo efeitos fixos de sexo, incubação e o efeito aleatório de valor genético infinitesimal) e um completo (incluindo todos os efeitos anteriores mais os efeitos de QTL). Ao analisar as características originais, ou seja, antes da formação das combinações lineares, foram mapeados seis QTLs significativos a 1% no genoma, quatro no GGA1 (peso vivo aos 35 dias, peso vivo aos 42 dias, gordura abdominal e peso do coração) e dois no GGA3 (peso vivo aos 35 e aos 42 dias de idade); três QTLs significativos a 5% no genoma, sendo um no GGA1 (peso da cabeça), um no GGA3 (peso das asas), e um no GGA8 (peso da moela); além de sete ligações sugestivas para diversas características. Ao mapear QTLs para os componentes principais, foram confirmados vários QTLs mapeados para as características originais, além de encontrar três QTLs e oito ligações sugestivas não mapeadas para as características originais.

Palavras-chave: frangos, galinha, genoma, variáveis canônicas

Introduction

Poultry industry has been able to reach high indices of productivity partly because of advances in animal breeding. Maintaining high levels of genetic gain in poultry is a challenge. In order to do so, molecular genetics may be applied to better understand the complex genetic architecture that controls quantitative traits as described by De Koning (2008).

A Brazilian population of F₂ animals has been repeatedly analyzed in QTL mapping studies in poultry. Some QTLs mapped in this population have been described by Nones et al. (2006), Ambo et al. (2009) and Campos et al. (2009). Procedures seeking to reduce the occurrence of false positives or negatives were not employed in these studies, except of course to adjust the level of type I error to an acceptable level.

The literature has suggested using principal components, as described by Mangin et al. (1998), to reduce false negatives (Type II error). This analysis consists in obtaining new characteristics, known as principal components or canonical variables, which are simply linear combinations of the characteristics of interest. According to Mangin et al. (1998), QTL mapping for principal components may find QTLs not mapped when each characteristic is analyzed independently. Therefore, the objective of the present study was to map QTLs for linear combinations of several characteristics of economic interest in poultry in order to independently identify QTLs that are not mapped using characteristics of economical interest.

Material and Methods

The reference population in the present study was developed at Embrapa Suínos e Aves in Concórdia, SC, Brazil using the line-cross model described by Knott & Haley (1992). A broiler line (TT) was selected for six generations for body weight, feed conversion, feed consumption, yields of carcass and parts, viability, fertility, eclodibility, reduced body fat and reduced occurrence of metabolic diseases. A layer line (CC) was selected for eight generations for egg production, egg weight, feed conversion, viability, sexual maturity, fertility, eclodibility, egg quality and reduced body weight. The F₂ chicken population used in this study was developed by crossing the broiler line (TT) with the layer line (CC). Four F₂ families (350 chickens) were selected for genotyping chromosomes 1, 3, 5, 8, 11 and 13 using microsatellite markers.

The F₂ population was raised as broilers, and received feed and water ad libitum. Animals received a starter ration from 1 to 21 days of age (21% crude protein and 3,150 Kcal/kg metabolizable energy), growth ration from 22 to 35 days of age (20% crude protein and 3,200 Kcal/kg ME), and a finishing ration from 36 to 41 days of age (18.5% crude protein and 3,200 Kcal/kg ME). Chicks were kept in collective boxes until 34 days of age. From the 35^th day until the 41^st day, animals were kept in individual cages to evaluate individual feed intake.

Linkage maps (Table 1) were obtained using CRIMAP software according to Green et al. (1990).

Live weight at 1, 35 and 42 (LW1, LW35, and LW42) days of age, weight gain (WG35/41), feed intake (FI35/41) and feed conversion (FC35/41) from day 35 to 41 were performance traits evaluated. Body weight at 42 days of age was measured after 6 hours of fasting and transportation to slaughter house. Carcass traits evaluated after four hours of refrigeration consisted of weight (g) of lungs, liver, heart, gizzard, breast, legs (weight of thigh and drumstick), carcass (without viscera, feet and head), residual carcass (weight of carcass without breast, wings and legs), wings, head, feet and abdominal fat and intestinal length (cm). Hematocrit values were determined in the laboratory using blood collected at slaughter.

To test sex effects, hatch, family and the covariable LW42, an analysis of variance was used, in which the residuals were initially evaluated assuming normal distribution, homogeneous variance and outliers using SAS Guided Data Analysis (SAS, 2002). Means were adjusted using the following model:

Y_ijk = S_i + H_j + F_k + c_ijkLW42 + e_ijk in which Y_ijk is the phenotypic value of a characteristic of the ijk^th individual, S_i is the fixed effect of the i^th sex, H_j is the fixed effect of the j^th hatch, F_k is the fixed effect of the k^th family of full sibs, c_ijk is the coefficient corresponding to the effect of the covariable LW42 of the ijk^th individual and e_ijk is the error, NID (0, σ). The different factors were evaluated using GLM procedure (SAS, 2002). The covariable LW42 was not included in the model for weights at 1, 35 and 42 days of age. Means, standard deviations, coefficients of variance, as well as maximal and minimal values were obtained using PROC MEANS (SAS, 2002).

Sequences of principle components analyses were carried out to separate the 20 traits into four groups. In this way, a characteristic was removed after each analysis in order to maximize the variance explained by the first principle component. It was established that each group would have a minimum of four traits, since simultaneously mapping multiple traits works reasonably well with up to three traits and the results are easier to interpret (Gilbert and Le Roy, 2003). Group 1 (G1) was formed by breast, thigh, wings, carcass and residual carcass weights as well as LW42 and LW35. Group 2 (G2) included head, feet, liver, heart and gizzard weights. Group three (G3) consisted of lungs weight, abdominal fat as well as weight gain and feed intake from 35 to 41 days of age. Lastly, group 4 (G4) included weight at 1 day of age, intestinal length, hematocrit values and feed conversion from 35 to 41 days of age.

Canonical variables were obtained from the original traits using PROC PRINCOMP (SAS, 2002). A correlation matrix was used to obtain the eigenvectors, since the original traits sometimes present variances of very different magnitudes. Principle components analysis seeks to obtain canonical variables (principle components) so that the first canonical variable corresponds to the greatest eigenvalue, the second canonical variable to the second greatest and so on (Table 2). Each characteristic used in the principal components analysis has a determined contribution to the canonical variable in question (Table 3). These values will be useful in interpreting the QTL map of the canonical variables, as it will be shown.

For QTL mapping, a likelihood ratio test was performed between two hierarchical models. The complete model is as it follows:

Y_ijk = µ + S_i + H_j + a_ijkLW42 + b_ijk I + c_ijk A + d_ijk D + e_ijk in which Y_ijk = the phenotypic value of a characteristic of the ijk^th animal, µ = the general mean of the phenotype considered, S_i = the fixed effect of the i^th sex, H_j = the fixed effect of the j^th hatch, a_ijk = the coefficient associated with the covariable LW42 of the ijk^th animal, b_ijk = the coefficient associated with the random effect of infinitesimal genetic value (I) of the ijk^th animal, c_ijk = the coefficient associated with the additive effect (A) of the QTL in the ijk^th animal, d_ijk = the coefficient associated with the dominant effect (D) of the QTL in the ijk^th animal and e_ijk = the random residual of the model in the ijk^th animal. The reduced model is nothing more than the complete model without the parameter of interest under study. It is important to mention that in order to avoid super-parametrization of the model, it was opted for maintaining only the significant effects in the analysis of variance.

The Qxpak software (Perez-Enciso & Misztal, 2004) was used for QTL mapping. Significance levels of mapped QTLs were estimated according to Lander and Kruglyak (1995), and the confidence interval for QTL location was obtained as described by Mangin et al. (1994). Phenotypic variation explained by the QTL was calculated as described by Sorensen et al. (2003).

Results and Discussion

Traits showed coefficients of variation higher than averages observed by Mohallem et al. (2008) (Table 4), which is due to the elevated effect of sex, family and hatch, but that is mainly due to the fact of dealing with an F2 population originating from crossing phenotypically distinct lines. Differences between minimal and maximal values were expected and desired, since variability in the F2 population is needed to map QTL.

Average live weight at 42 days of age is lower than the average observed in commercial lines by Havenstein et al. (2003), which was expected from a broiler × layer cross. On the other hand, this value is near that observed by Zerehdaran et al. (2004) in an experimental F₉ population whose founder lines varied for live weight.

QTLs of small effect, i.e. those that explain reduced percentages of phenotypic variance, are more difficult to map (Gilbert & Le Roy, 2003). In these cases, larger samples, denser genetic maps or alternative QTL mapping methods are needed for increasing power of statistical test. Mangin et al. (1998) demonstrated that using canonical variable is a good approach for increasing the power of detecting QTLs. To verify the effect of this alternative methodology, QTL of original traits, i.e. before obtaining linear combination, were mapped (Table 5) for consequent comparison to principal component analysis results. Percentages of phenotypic variance explained by the QTLs varied from 1.37 to 7.98%. It may be expected that some QTLs may not have been mapped for original traits in the present study (350 animals and medium density of a marker every 17.7 cM) due to the reduced power of the statistical test.

Adjustment of infinitesimal genetic value is extremely important to the present analysis, since non-isolation of the polygenic effect may obtain biased estimates of the QTL effects. The infinitesimal effect was not adjusted in works by Nones et al. (2006), Ambo et al. (2009) and Campos et al. (2009) which may have led these studies to detect false positive QTLs or to overestimate the effects of QTLs. For example, Campos et al. (2009) identified QTL for abdominal fat at 115cM on chromosome 3, explaining 4.08% phenotypic variance of this trait with significant additive, dominance and imprinting effects. This may be the QTL that was detected in the present study at 115 cM, although with a significant additive effect and explaining approximately 2.33% phenotypic variance. These differences may result from adjustment in the infinitesimal genetic effect used in the present study and not considered in the work by Campos et al. (2009).

It is important to stress that canonical variables were used to identify new QTLs that had not been mapped in analysis of the original traits, i.e. before obtaining linear combinations. Consequently, interest of the present study is restricted to increasing the power of the test as described in Weller et al. (1996), Magin et al. (1998) and Gilbert & Le Roy (2003), and not in any specific linear combination. Two highly significant linkages, four significant linkages and ten suggestive linkages were mapped for canonical variables (Table 6). Since the criteria for suggestive linkage leads to one false positive throughout the genome (Lander & Kruglyak, 1995), it is possible to conclude that at least nine of these suggestive linkages are real QTLs. Concerning principal components, only QTL positions were described because both dominance and imprinting effects were not significant, while additive effects are associated with linear combinations, but no with original traits.

Four QTLs and one suggestive linkage were mapped for Group 1 canonical variables. The first QTL is located at 79 cM of GGA1 and it is associated with PC1 (Table 6). Apparently, this QTL has an effect on all traits used in forming this component since contributions of seven original traits for PC1 are very similar (Table 3). However, QTL results from original traits (Table 5) show that QTLs in this region of GGA1 only exist for LW35 and LW42. Based on the considerable fall in nominal probability, evidenced by the significance level of the test for LW35 and LW42 versus the test for PC1 in GGA1, it is possible that the QTL mapped here is the same one found in analyses of LW42 and LW35, and that the other traits included in the linear combination reduced the nominal probability.

The second QTL for Group 1 was mapped for 376 cM of GGA1 and it is associated with PC4 (Table 6). It is a novel QTL since this region of GGA1 did not have any suggestive linkages in the analysis of the original traits composing Group 1. This principal component greatly contributes to the weight of wings, thighs and drumsticks and breast (Table 3). However, this region of GGA1 should be refined to identify and to explain the effect of this QTL and the characteristic(s) which it affects.

The third QTL for G1 was mapped at 62 cM of GGA3 and it is associated with PC1. Its effect was highly significant and near the region where other highly significant QTLs for live weight at 35 and 42 days of age have been mapped (Table 5). Since the probability of the nominal test was not altered, it is possible that this QTL has a pleiotropic effect on some traits of this principal component or even on all of them. The fourth QTL of G1 was mapped at 126 cM of GGA3 and it is associated with PC3. The power of the test can be seen now since this region of GGA3 did not have QTL for the traits that compose G1. This principal component received great contribution from the residual carcass weight and LW35 (Table 3). Consequently, inclusion of more markers in this region of GGA3 may help identify the traits influenced by this QTL.

Suggestive linkage for canonical variables of G1 was mapped at 152 cM from GGA3 and it is associated with PC4. Wings, thighs and breast are traits with important contribution for PC4 (Table 3). Since a QTL was mapped on GGA3 for wings at 151 cM and suggestive linkage for thighs at 146 cM (Table 5), this suggestive linkage at 152 cM for PC4 confirmed results found with original trait.

Two QTLs and a suggestive linkage were mapped for G2 (Table 6). The first QTL is located at 83 cM of GGA1 and it is associated with PC1. This component received great contribution from head, feet, liver, heart and gizzard weights (Table 3). Furthermore, a QTL for heart weight was mapped at this region at 74cM (Table 5). Evidence from Navarro et al. (2005) indicates QTL for heart weight near this region. Since there was no reduction in nominal probability, it is possible that this QTL also has an effect on some other trait(s) that contribute to the variation explained by this principal component. The second QTL influencing traits of G2 is located at 111 cM of GGA3 and also has an effect on PC1. It is a new QTL, since no QTL were mapped in GGA3 in the analysis of the original traits composing G2. However, it is possible that this QTL has a pleiotropic effect on some traits in this group. Suggestive linkage with some influence on G2 was mapped at 117cM from GGA1 and it is associated with PC2 (Table 6). This component was mainly influenced by heart and gizzard weights (Table 3) which did not have any QTL in this region of GGA1. It is possible that this QTL has a pleiotropic effect on these two traits, although Navarro et al. (2005) mapped QTL only for gizzard weight in this same region. Consequently, inclusion of more markers in this region of GGA1 or increasing the F2 population may help make the effects of this QTL clear.

Five suggestive linkages for the G3 canonical variable were mapped (Table 6). The first for PC3 at 91cM of GGA1 confirms the QTL mapped at 91 cM for weight of abdominal fat (Table 5). This component is strongly influenced by variation in lungs weight and to a lesser degree by variations in weight of abdominal fat, weight gain and feed intake (Table 3). The contribution of lungs weight to PC2 can explain the reduced nominal probability, which stopped being highly significant in QTL analysis of weight of abdominal fat to become suggestive in analysis of PC3 in G3. The second suggestive linkage is located at 180 cM of GGA1 and has an effect on PC4. Not even suggestive linkages for original traits composing G3 have been found in this region of GGA1. This suggestive linkage may possibly affect weight gain and/or feed intake. Both of which greatly contribute to the amount of variation PC4 explains in G3 (Table 3). Therefore, refining this GGA3 region may help identify and clarify the effects of this suggestive linkage. The third suggestive linkage of G3 is associated with PC2 and located at 112 cM of GGA3. This linkage may result from the effects of suggestive linkage mapped at 115 cM of GGA3 for abdominal fat (Table 5), since variation in abdominal fat has important contribution to PC2 in G3 (Table 3). The fourth suggestive linkage was mapped at 166 cM of GGA3 and it is associated with PC1. This principal component received great contribution from the four traits forming this group (Table 3). When evaluated singularly, no QTL were mapped for any of the four traits in the group. Thus, the linkage mapped may have an effect on any of the traits in G3 or even a pleiotropic effect on some of them. The last suggestive linkage for G3 was mapped at 110 cM of GGA3 and it is associated with PC4, which received strong contribution from weight gain and feed intake. When these traits were analyzed individually, no QTL were mapped. More studies of this region of GGA3 are needed since feed intake and weight gain are traits of economic interest.

Three suggestive linkages were mapped for G4 canonical variables (Table 6). The first is located at 32 cM of GGA1 and has an effect on PC1. This suggestive linkage receives an important contribution from hematocrit and weight at 1 day (Table 3). However, no QTL were mapped for these traits, when analyzed individually (Table 5). The remaining suggestive linkages for G4 were observed for PC2 at 340 cM of GGA1 and 36 cM of GGA11. These principle components received great contribution from variations in intestine length and feed conversion (Table 3). Both linkages may be false positives or novel QTLs found through linear combinations since no QTLs or suggestive linkages for G4 traits were mapped for these regions of GGA1 and GGA11. Consequently, these regions should be studied further.

Interpretation of QTL analysis results for canonical variables is difficult when the original variables are interesting for studying and not the linear combinations. For example, the QTL for PC1 in G1 at 79cM of GGA1 had a negative additive effect (data not shown). Analysis and comparison of QTL analysis results of the original traits leads to the conclusion that this QTL is the same as that mapped with a positive additive effect for LW35 and LW42. However, the effects of QTLs associated with linear combinations do not have an immediate interpretation. The same problem is observed for the percent of variance explained by the QTL. Thus, in cases in which the linear combination is not very important, QTL results for canonical variables may only be used to identify regions with potential QTLs. In the present study, use of canonical variables allowed the identification of three QTLs not mapped with the original traits and eight new suggestive linkages. These regions may now be refined with a greater number of molecular markers in the expectation of obtaining more information about the effect of these QTLs on the original traits. However, the use of linear combinations will not necessarily increase the number of QTLs or the suggestive linkages found and it may even reduce this number. For example, the suggestive linkages for wings and heart weights respectively located at 96 cM of GGA11 and 57 cM of GGA13 were not identified in the principal components analysis. Consequently, the linear combinations should be carefully defined to maximize the capacity of locating QTLs.

Conclusions

Group composition, for principal component analysis, should be carefully planned, to avoid reduction of the power of the test in the mapping analysis.

All significant effects found for canonical variables must be compared to those found when working with the original traits.

Mapping loci of canonical variables helps identify significant regions that are not mapped in isolated analysis of the original traits.

Acknowledgments

The authors would like to thank FAPESP for granting scholarships to Luis Fernando Batista Pinto and Kátia Nones; Dr. Jerry Dodgson and Dr. Hans Cheng (United States Poultry Genome Project) for donating primers; CNPq for granting scholarships to Ana Silva Alves Meira Tavares Moura and Luiz Lehmann Coutinho; and to EMBRAPA/PRODETAB for their financial contribution.

Received July 28, 2009 and accepted March 17, 2010.

Corresponding author: [e-mail]

References

• AMBO, M.; MOURA, A.S.A.M.T.; LEDUR, M.C. et al. Quantitative trait loci for performance traits in a broiler x layer cross. Animal Genetics, v.40, p.200-208, 2009.
• CAMPOS, R.L.R.; NONES, K.; LEDUR, M.C. et al. Quantitative trait loci associated with fatness in a broiler x layer cross. Animal Genetics, v.40, p729-736, 2009.
• DE KONING, D.J. Application of molecular information in sustainable animal breeding. Revista Brasileira de Zootecnia, v.37, p.122-126, 2008.
• GREEN, P.; FALLS, K.; CROOKS, S. Documentation for CRI-MAP, version 2.4 St. Louis: Washington School of Medicine, 1990. Available in: <http://linkage.rockefeller.edu/soft/crimap/> Access at: 7/23/2009.
• GILBERT, H.; LE ROY, P. Comparison of three multitrait methods for QTL detection. Genetic, Selection and Evolution, v.35, p.281-304, 2003.
• HAVENSTEIN, G.B.; FERKET, P.R.; QURESHI, M.A. Carcass composition and yield of 1957 versus 2001 broilers when fed representative 1957 and 2001 broiler diets. Poultry Science, v.82, p.1509-1518, 2003.
• KNOTT, S.A.; HALEY, C.S. Maximum likelihood mapping of quantitative trait loci using full-sib families. Genetics, v.132, p.1211-1222, 1992.
• LANDER, E.S.; KRUGLYAK, L. Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nature Genetics, v.11, p.241-247, 1995.
• MANGIN, B.; GOFFINET, B.; REBAÏ, A. Constructing confidence intervals for QTL location. Genetics, v.138, p.1301-1308, 1994.
• MANGIN, B.; THOQUET, P.; GRIMSLEY, N. Pleiotropic QTL analysis. Biometrics, v.54, p.88-99, 1998.
• MOHALLEM, D.F.; TAVARES, M.; SILVA, P.L. et al. Avaliação do coeficiente de variação como medida da precisão em experimentos com frangos de corte. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, v.60, p.449-453, 2008.
• NAVARRO, P.; VISSCHER, P.M.; KNOTT, S.A. et al. Mapping of quantitative trait loci affecting organ weights and blood variables in a broiler x layer cross. British Poultry Science, v.46, p.430-442, 2005.
• NONES, K.; LEDUR, M.C.; RUY, D.C. et al. Mapping QTL on chicken chromosome 1 for performance and carcass traits in a broiler x layer cross. Animal Genetics, v.37, p.95-100, 2006.
• PÉREZ-ENCISO, M.; MISZTAL, I. Qxpak: a versatile mixed model application for genetical genomics and QTL analyses. Bioinformatics, v.20, p.2792-2798, 2004.
• SORENSEN, P.; LUND, M.S.; GULDBRANDTSEN, B. et al. A comparison of bivariate and univariate QTL mapping in livestock populations. Genetics, Selection and Evolution, v.35, p.605-622, 2003.
• STATISTICAL ANALYSIS SYSTEM - SAS. The SAS system for windows Cary: SAS Institute Inc., 2002. (CD-ROM).
• WELLER, J.I.; WIGGANS, G.R.; VAN RADEN, P.M. Application of a canonical transformation to detection of quantitative trait loci with the aid of genetic markers in a multi-trait experiment. Theoretical and Applied Genetics, v.92, p.998-1002, 1996.
• ZEREHDARAN, S.; VEREIJKEN, A.L.J.; VAN ARENDONK, J.A.M. et al. Estimation of genetic parameters for fat deposition and carcass traits in broiler. Poultry Science, v.83, p.521-525, 2004.

*

In memoriam.

Publication Dates

• Publication in this collection
  25 Nov 2010
• Date of issue
  Nov 2010

History

• Received
  28 July 2009
• Accepted
  17 Mar 2010