key: cord-1004978-qhea8t2c
authors: Xiong, P.; Zeng, X.; Song, M. S.; Jia, S. W.; Zhong, M. H.; Xiao, L. L.; Lan, W.; Cai, C.; Wu, X. W.; Gong, F. L.; Wang, W.
title: Lack of association between HLA‐A, ‐B and ‐DRB1 alleles and the development of SARS: a cohort of 95 SARS‐recovered individuals in a population of Guangdong, southern China
date: 2008-01-03
journal: Int J Immunogenet
DOI: 10.1111/j.1744-313x.2007.00741.x
sha: b3fce64f2cc6dd580e65be72e26c0531c414f313
doc_id: 1004978
cord_uid: qhea8t2c

Severe acute respiratory syndrome (SARS), caused by infection with a novel coronavirus (SARS‐CoV), was the first major novel infectious disease at the beginning of the 21st century, with China especially affected. SARS was characterized by high infectivity, morbidity and mortality, and the confined pattern of the disease spreading among the countries of South‐East and East Asia suggested the existence of susceptible factor(s) in these populations. Studies in the populations of Hong Kong and Taiwan showed an association of human leucocyte antigen (HLA) polymorphisms with the development and/or severity of SARS, respectively. The aim of the present study was to define the genotypic patterns of HLA‐A, ‐B and ‐DRB1 loci in SARS patients and a co‐resident population of Guangdong province, southern China, where the first SARS case was reported. The samples comprised 95 cases of recovered SARS patients and 403 unrelated healthy controls. HLA ‐A, ‐B and ‐DRB1 alleles were genotyped using polymerase chain reaction with sequence‐specific primers. The severity of the disease was assessed according to the history of lung infiltration, usage of assisted ventilation and occurrence of lymphocytopenia. Although the allelic frequencies of A23, A34, B60, DRB1*12 in the SARS group were slightly higher, and A33, ‐B58 and ‐B61 were lower than in the controls, no statistical significance was found when the Pc value was considered. Similarly, no association of HLA alleles with the severity of the disease was detected. Thus, variations in the major histocompatibility complex are unlikely to have contributed significantly to either the susceptibility or the severity of SARS in the population of Guangdong.

Severe acute respiratory syndrome (SARS) originated in southern China in November 2002 , reached Hong Kong in February 2003 , and thereafter spread rapidly to 29 countries/regions on five continents. By the end of the epidemic, the global cumulative patient total was 8098 with 774 deaths. Seven Asian countries/regions were among the top 10 populations affected, China and Hong Kong accounting for 87% of all cases and 84% of deaths from SARS (Lam et al., 2003) . SARS was mostly confined to the populations of South-East and East Asia (Mainland Chinese, Hong Kong residents, Vietnamese, Singaporeans and Taiwanese). The disease is caused by infection with a novel single-stranded RNA coronavirus (Peiris et al., 2003) and starts with an influenza-like illness characterized by non-specific, systemic symptoms, followed by the rapid development of an unspecific bronchopneumonia associated with lower respiratory tract or gastrointestinal symptoms. Most patients recovered after 1 or 2 weeks, but some developed an acute respiratory distress syndrome . The confined geographical and ethnic spread pattern of the disease suggested the possibility of genetic variations that had influenced susceptibility to infection with SARS-CoV and the development of the disease itself.

Recent reports have suggested that a single nucleotide polymorphism (SNP) in the 3′ UTR of the OAS-1 and MxA promoter region is associated with host susceptibility to SARS infection in the Han Chinese population (Hamano et al., 2005; He et al., 2006) . Also, an association study that examined polymorphisms in the cytokine response genes [interferon (IFN)-gamma, tumour necrosis factor (TNF)-alpha and interleukin (IL)-10] showed that the 874 A/T IFN-gamma SNP was a risk factor for SARS infection, although no association was seen with either TNF-alpha or IL-10 polymorphisms (Chong et al., 2006) . Chan et al. (2006) reported that individuals who were homozygous for CLEC4M (a binding receptor for SARS Co-V) were less susceptible to SARS (Chan et al., 2006) . Conversely, association studies failed to demonstrate that ACE2 SNPs were implicated in SARS , and that the ACE I/D polymorphism was neither directly related to increased susceptibility to SARScoronavirus infection nor the clinical outcome of SARS patients (Chan et al., 2005) .

The human leucocyte antigen (HLA) system is widely used to investigate the aetiology of infectious diseases and autoimmune disorders. Studies in Hong Kong (Ng et al., 2004) and Taiwan (Lin et al., 2003) populations showed an association of HLA-B*0703-DRB1*0301 with the development of SARS, and HLA-B46 with the severity of SARS, respectively. However, the data on which these reports were based were limited in number, and no comparable data from Mainland China have been available. In this study, samples from 95 SARS patients and 403 healthy controls were typed for HLA-A, -B and -DRB1 alleles to determine if an association existed between HLA variations and SARS.

The cohort of 95 SARS-recovered individuals comprised patients admitted to the Guangdong Provincial TCM Hospital in Guangzhou (82 cases) and the Peking University Shenzhen Hospital in Shenzhen (13 cases) between March and June 2003. All patients were diagnosed according to the World Health Organization (WHO) case definition of probable SARS case (WHO, 2003) . All of the SARS subjects were ethnic Chinese and none were biological relatives (only the first infected patient in a family was included if several family members had been diagnosed). The age range of the patients was 19-72 years, and the male:female ratio was 0.76. All 95 SARS patients showed radiographic evidence of lung infiltrates and were seropositive for SARS-CoV antibody during the SARS episode. A sex-and age-matched healthy control group of 403 healthy volunteers resident in the same geographical area were also genotyped. The mean age of the controls was 45.6 years, and the male:female ratio was 0.73. Informed written consent was obtained from all tested subjects.

The hospital records of the 95 SARS cases indicated that 63 patients (66.3%) had lymphocytopenia (< 1 × 10 9 L -1 ), with two patients (2.11%) showing diffuse infiltration of both lungs. Auxiliary ventilation was applied by mask in 30 cases (31.6%), with intubation additionally used in two cases (2.1%). Lymphocytopenia, lung involvement and auxiliary ventilation were considered as indicative of disease severity in our study.

Peripheral blood samples were collected in EDTA vacutainers, with genomic DNA prepared by standard phenol/ chloroform/isoamyl alcohol extraction and ethanol precipitation for immediate storage at -20 °C. DNA quantity and quality were assessed by spectroscope (Pharmacia Biotech, Uppsala, Sweden) absorbance at 260 and 280 nm. The identification of HLA-A, -B and -DRB1 alleles in all subjects was performed by DNA typing, using polymerase chain reaction with sequence-specific primers (PCR-SSP) based on the protocol of Bunce et al. (1995) , with slight modification to correspond with serological typing results. The conversion was based on the WHO nomenclature table of correspondence (Marsh et al., 2002) .

Briefly, PCR mixtures consisted of 67 mm Tris base pH 8.8, 16.6 mm ammonium sulphate, 2 mm magnesium chloride, 0.01% (v/v) Tween 20, 200 pmol of each dNTP, 1-4 pmol of each allele-specific primer, 0.1 pmol of DRB1 control primers or 2 pmol of APC control primers, between 0.1 and 0.01 pg DNA and 0. 1875 U Taq polymerase. The primer concentrations in the primer mixes were empirically titrated with positive and negative samples to ensure maximum efficiency. PCR amplications were carried out in an MJ Research PTC-100 (Waltham, MA, USA) machine. The cycling parameters for 13 μl reactions in rapid-cycling PCR machines were: 1 min at 96 °C followed by five cycles of 25 s at 96 °C, 45 s at 70 °C, 45 s at 72 °C, followed by 21 cycles of 25 s at 96 °C, 50 s at 65 °C, 45 s at 72 °C, followed by four cycles of 25 s at 96 °C, 60 s at 55 °C and 20 s at 72 °C.

To assess the association of HLA variations with SARS infection, the frequencies of HLA class I and II alleles were estimated by direct counting in both the patient and the control groups. The frequencies of HLA alleles (HLA antigens) in patients with SARS and healthy controls were each tested for deviation from Hardy-Weinberg equilibrium, and compared using a χ 2 test (continuity correction for the χ 2 test was used where a small number of alleles was observed (E < 5) and a two-tail Fisher's exact test was used to calculate the P-value when E < 1. Corrected Pvalues (Pc) were obtained by multiplying the number of antigens or alleles (Edward, 1974 : reference not listed) at each locus: A (n = 15), B (n = 32) and DRB1 (n = 13). A Pc value < 0.05 was considered statistically significant. Odds ratios (OR) rather than relative risks (RR) were quoted since ORs are most commonly used in HLA association studies. The values were obtained from standard contingency table analysis using Haldane's modification of Woolf's method (Haldane, 1956 ). An OR value > 1.0 implies that subjects carrying the antigen are at an increased risk of a disease (a positive association), while an OR value < 1.0 implies that subjects carrying the antigen are at a lower risk of a disease than those not carrying the antigen (a negative association).

To assess the association of HLA alleles with the severity of our SARS cases, SARS patients were divided into assisted ventilation (mask and intubation) (n = 30) and non-assisted ventilation groups (n = 65). The patients were also grouped according to the occurrence of lymphocytopenia (< 1 × 10 9 L -1 ) (n = 63). Four grades were scored to indicate the area of lung infiltration as indicated by the chest radiographs: grade 1 indicated 1/4 (n = 25), 2 for 2/4 (n = 58), 3 for 3/4 (n = 10) and 4 for 4/4 (n = 2). The association between lymphocytopenia and HLA variation was analysed using a one-way anova using the SNK test for multiple comparisons. Mordridit score analysis was applied when comparing associations between lung infiltrates and HLA. All analyses were performed using spss 11.0. All probabilities were two tailed with a P-value < 0.05 considered statistically significant.

The were lower. However, after correction of the Pvalue by multiplying the number of alleles at each locus, all statistical significance was lost. With regard to the association of HLA polymorphism with severity of SARS, the auxiliary ventilation, areas of lung infiltrate and lymphocytopenia all seemed to occur randomly with HLA-A, -B, and -DRB1 alleles, since none of the associations tested were significant (data not shown).

No deviation from Hardy-Weinberg equilibrium was observed when the frequencies of HLA alleles were tested in the SARS patients or the healthy controls (P > 0.05). 

Classical HLA classes I and II molecules are essential to present antigenic peptide (including viral peptide) for Tcell recognition and thus initiate an adaptive immune response. To assist the survival of a population, the genetic cluster of HLA loci manifests a striking level of polymorphism to protect the host from any pathogen mutation eluding host immune surveillance. Antigenic peptides binding to different alleles of HLA molecules vary in their affinity, with the HLA allotype thus affecting the host immune response to a given antigen. Studies have indicated that HLA variations are associated with susceptibility or resistance to a wide range of infections, including malaria, tuberculosis, leprosy, HIV and hepatitis virus persistence (De Vries et al., 1979; Hill, 1998; Liu et al., 2003) . SARS has been shown to be caused by infection with SARS-CoV, with major lung involvement and consequent high infectivity and mortality. Although many questions remain to be answered, immunological mechanisms are believed to be involved in the development of the disease. As previously noted, in the 2002-03 SARS epidemic, the disease was mostly confined to South-East and East Asian populations, suggestive of greater susceptibility to SARS-CoV.

Guangdong province, where the SARS epidemic started, has a population of 85.2 million. The first case of SARS that fulfilled the WHO case definition was reported in Foshan, a city about 20 km from Guangzhou (formerly Canton), the provincial capital Xu et al., 2004) . Phylogenetic studies have indicated that the populations of Guangdong, Hong Kong and Taiwan populations largely share common origins (Chu et al., 1998) . For example, over 98% of the Hong Kong population are Chinese, and the majority have ancestral origins in neighbouring Guangdong province and are of Han Chinese ethnicity (Chang & Hawkins, 1997) . The present-day population of Taiwan comprises indigenous Taiwanese groups, and people of Minnan and Hakka ethnicity who are the descendants of settlers from the provinces of Fujian and Guangdong during the last few centuries (Lin et al., 2001) . It is of interest that the SARS cases in Taiwan mainly consisted of southern Chinese origin Taiwanese rather than the indigenous peoples (Lin et al., 2003) .

Similarities in the genetic profiles of these three populations, Guangdong, Hong Kong and Taiwan offers the opportunity of identifying possible associations between HLA alleles and SARS infection. Hypothetical SARSassociated HLA alleles among the three populations are listed in Table 4 , and it is obvious that the SARS-associated HLA alleles that have been reported in any single population could not be verified in either of the other two populations. It also seems that the Guangzhou population was more susceptible to the SARS epidemic irrespective of their HLA profiles. The mutation of SARS-CoV could give rise to variations and discrepancies in the results of the HLA association studies in the three populations, and large-scale studies are essential for the demonstration of an association of HLA variation with susceptibility to a disease. On the other hand, mutation of SARS-CoV might result in changes of its T-cell receptor epitope, which would require a specific HLA allele(s) for presentation. Because of the poor fidelity of RNA-dependent RNA polymerase, genetic variation typically gives rise to a heterogeneous virus pool in RNA populations, including coronaviruses. This feature makes viruses highly adaptable and contributes to difficulties in preventing and controlling viral diseases. The virus isolated from the Guangzhou patients is the prototype of the SARS-CoV, with a reported moderate mutation rate, that subsequently was found in Hong Kong and countries around the world (Zhao et al., 2004) . The Guangzhou SARS-CoV has the highest homology (> 99%) with the sequences of the SARS-CoV that were isolated from SARS cases in Hong Kong, Beijing, Singapore, Canada and the USA, suggesting that all of these SARS-CoV are closely related , therefore mutation of SARS-CoV could not be considered as a major contributory factor in the lack of HLA association observed in the present study.

During the SARS episode in Guangdong, six of the 95 SARS patients died of the disease. No samples from these deceased patients were available for our study. However, a history of lymphocytopenia, lung involvement and auxiliary ventilation were considered as surrogate measures of the severity of the disease, and no significant association of HLA variations with these parameters was observed among the patients, who instead displayed a random distribution of HLA-A, -B and -DRB1 alleles.

In conclusion, the confined regional spread of the 2002-03 SARS epidemic appears to have been mostly dependent on factors such as relative population density, the frequency and extent of travel of affected individuals, and the rapid introduction and uptake of protective measures, rather than to HLA variation in the population.

Phototyping: comprehensive DNA typing for HLA-A, B, C, DRB1, DRB3, DRB4, DRB5 & DQB1 by PCR with 144 primer mixes utilizing sequence-specific primers (PCR-SSP)

Absence of association between angiotensin converting enzyme polymorphism and development of adult respiratory distress syndrome in patients with severe acute respiratory syndrome: a case control study

Homozygous L-SIGN (CLEC4M.) plays a protective role in SARS coronavirus infection

HLA class I and class II frequencies of a Hong Kong Chinese population based on bone marrow donor registry data

The interferon gamma gene polymorphism +874 A/T is associated with severe acute respiratory syndrome

Genetic relationship of populations in China

Genetic control of survival in epidemics

HLA and disease, the detection of associations

The estimation and significance of the logarithm of a ratio of frequencies

Polymorphisms of interferon-inducible genes OAS-1 and MxA associated with SARS in the Vietnamese population

Association of SARS susceptibility with single nucleic acid polymorphisms of OAS1 and MxA genes: a case-control study

The immunogenetics of human infectious diseases

Identification of an alternative 5′-untranslated exon and new polymorphisms of angiotensin-converting enzyme 2 gene: lack of association with SARS in the Vietnamese population

The origin of Minnan and Hakka, the so called 'Taiwanese', inferred by HLA study

Association of HLA class I with severe acute respiratory syndrome coronavirus infection

Association of polymorphisms in human leukocyte antigen class I and transporter associated with antigen processing genes with resistance to human immunodeficiency virus type 1 infection

Nomenclature for factors of the HLA system

Association of human-leukocyte-antigen class I (B*0703) and class II (DRB1*0301) genotypes with susceptibility and resistance to the development of severe acute respiratory syndrome

Coronavirus as a possible cause of severe acute respiratory syndrome. Lancet, 361, 1319. World Health Organization (WHO) (2003) Case definitions for surveillance of severe acute respiratory syndrome (SARS)

Epidemiologic clues to SARS origin in China

Moderate mutation rate in the SARS coronavirus genome and its implications

Epidemiology and cause of severe acute respiratory syndrome (SARS) in Guangdong, People's Republic of China

Beijing Natural Science Foundation (KZ200610025014). We thank Professor Alan H Bittles of Edith Cowan University, Australia, for kind comments and language editing of the manuscript.