Conceived and designed the experiments: WBX. Performed the experiments: YZ STX HLW ZZ YXJ CYL XJZ LWS JHZ PSL YH DXF ZYZ CYW XQF HYZ LL XDS WT YW YL HG HT JTM SYG SW YF FB JFL YS SJZ YYM SWW SDZ FCL ZRD ZHY. Analyzed the data: YZ STX HLW. Contributed reagents/materials/analysis tools: PAR DF YJ WJB WBX. Wrote the paper: YZ PAR DF YJ WJB WBX.
The incidence of measles in China from 1991 to 2008 was reviewed, and the nucleotide sequences from 1507 measles viruses (MeV) isolated during 1993 to 2008 were phylogenetically analyzed. The results showed that measles epidemics peaked approximately every 3 to 5 years with the range of measles cases detected between 56,850 and 140,048 per year. The Chinese MeV strains represented three genotypes; 1501 H1, 1 H2 and 5 A. Genotype H1 was the predominant genotype throughout China continuously circulating for at least 16 years. Genotype H1 sequences could be divided into two distinct clusters, H1a and H1b. A 4.2% average nucleotide divergence was found between the H1a and H1b clusters, and the nucleotide sequence and predicted amino acid homologies of H1a viruses were 92.3%–100% and 84.7%–100%, H1b were 97.1%–100% and 95.3%–100%, respectively. Viruses from both clusters were distributed throughout China with no apparent geographic restriction and multiple co-circulating lineages were present in many provinces. Cluster H1a and H1b viruses were co-circulating during 1993 to 2005, while no H1b viruses were detected after 2005 and the transmission of that cluster has presumably been interrupted. Analysis of the nucleotide and predicted amino acid changes in the N proteins of H1a and H1b viruses showed no evidence of selective pressure. This study investigated the genotype and cluster distribution of MeV in China over a 16-year period to establish a genetic baseline before MeV elimination in Western Pacific Region (WPR). Continuous and extensive MeV surveillance and the ability to quickly identify imported cases of measles will become more critical as measles elimination goals are achieved in China in the near future. This is the first report that a single endemic genotype of measles virus has been found to be continuously circulating in one country for at least 16 years.
Measles continues to be a leading cause of childhood morbidity and mortality in developing countries and an outbreak threat in the majority of countries, despite the availability of an effective vaccine for over 40 years
Measles virus (MeV) is a negative-sense, single-stranded RNA virus in the genus
Source: National Notifiable Diseases Reporting System (NNDRS). Blue bars indicate the reported measles cases and red solid diamonds indicate the incidence (/100,000 ) of each year. X-axis denotes year, left-hand y-axis denotes reported number of cases and right-hand y-axis denotes the incidence per 100,000.
Molecular epidemiological studies of MeV have made significant contributions to measles control efforts, by providing a means to confirm the sources of virus or suggest a source for unknown-source cases as well as to establish links, trace the transmission pathway, allow discrimination of imported from indigenous measles cases, or classify suspected cases as vaccine reactions
Measles is a vaccine-preventable disease, but is still a major killer of children worldwide. Measles caused many outbreaks in China before the 1965 introduction of liquid measles vaccine (MV) administered as 1 dose to infants aged >8 months. The National Expanded Program on Immunization (EPI) was established in 1978 to allow routine immunization with MV to cover the entire country
The WHO Regional Committee of the Western Pacific Region (WPR) formally declared a measles elimination goal in 2005 and established a target date of 2012 for regional measles elimination
This study did not involve human participants or human experimentation; the only human materials used were throat swab or urine samples collected from patients with acute, febrile maculopapular rash at the instigation of the Ministry of Health P. R. of China for public health purposes, and written informed consent for the use of their clinical samples was obtained from all patients involved in this study. This study was approved by the second session of the Ethics Review Committee of the Chinese Center for Disease Control and Prevention.
Reported numbers and incidence of measles cases and deaths in this report were obtained from the National Notifiable Diseases Reporting System of China CDC (NNDRS). Population denominators for calculation of incidence and mortality rates were determined on the basis of data reported by the National Bureau of Statistics. Epidemiologic data were analyzed by using Microsoft Excel.
Throat swab or urine samples were collected from patients with acute, febrile maculopapular rash. All clinical samples were collected within five days of rash onset and transported to the laboratory in accordance with standard protocols
RNA was extracted from 250 μl of infected cell lysate using a Trizol reagent, following the manufacturer's instructions. For all virus isolates, RT-PCR amplification was performed using previously described primers to amplify a 600 bp fragment in the N gene which included the 450 bp fragment recommended for genotyping
Sequences of the PCR products were derived by automated sequencing and the BigDye terminator v3.0 chemistry according to the manufacturer's protocol in both sense and antisense strands by an automated ABI PRISM™ 3100 DNA Sequencer (Perkin Elmer). Sequence proof reading and editing was conducted with Sequencer™ (Gene Codes Corporation). Sequence data were analyzed using version 7.0 of Bioedit and phylogenetic analyses were performed with the evolution model of Maximum Composite Likelihood using Mega 4. The robustness of the groupings was assessed using bootstrap re-sampling of 500 replicates and the trees were visualized with Mega programs.
Synonymous and non-synonymous substitutions: to determine whether the detected mutations in wild type MeV were the consequence of selection pressure, the frequencies of synonymous (dS) or non-synonymous (dN) was calculated, which is termed as ω, the (average) dN/dS rate ratio. An ω<1 suggests purifying (negative) selection, whereas an ω>1 is indicative of positive selection because that the rate of fixation is higher than the background rate of mutation, which defies neutrality (ω = 1)
The nucleotide sequences of 135 viruses, representative of the 333 Chinese MeV of 2008 that were isolated in this study, were deposited in the GenBank database under accession numbers JN380206-JN380341. An additional 692 viruses used in this study were previously deposited in GenBank (accession numbers HAF045191-HAF045218, DQ356683- DQ356873, EU557194- EU 557238, FJ602549- FJ 602674, GU237175- GU 237481).
The system of reporting infectious diseases, including measles, to the National Notifiable Disease Reporting System (NNDRS) had been established in the 1950's in China. The cases of disease were reported by posting a card monthly from the hospitals to the national level (sequentially, from county Centers for Disease Control (CDC), to prefecture CDC and provincial CDC, to the national level). Since 2004, the NNDRS has been updated to allow hospitals or local CDCs to directly report all suspected measles cases immediately to the national level through a web-based, real-time reporting system.
With the attainment of childhood immunization goals of 2 doses of measles vaccine since 1986, the measles morbidity and mortality have decreased markedly. Subsequently, measles epidemics occurred about every 3–5 years, with the incidence fluctuating around 8/100,000. The average annual reported incidence in 1994–2004 was below 7.5/100,000 (average, 5.84/100,000; range, 4.5[1998]–7.4[1994] /100,000). But measles cases resurged in 2005–2008 and the average annual reported incidence (average, 8.95/100,000; range, 7.7[2006]–10.1[2008]/100,000) represented a 53.3% increase compared with 1994–2004 (
A total of 1507 measles viral isolates were identified from throat swab or urine specimens obtained from 30 of 31 provinces in China (except for Tibet), between 1993 and 2008 (
These trees are based on the WHO standard sequence window within the N gene. Panel A. Phylogenetic tree of 87 representative measles isolates from China during 1993–2008 compared to the WHO reference sequences for each genotype. Sequences from Chinese viruses of H1a cluster from 1993–2008 are indicated by blue and sequences from Chinese viruses of H1b cluster from 1994–2005 are indicated by fluorescence green, and WHO reference strains are indicated by red. All isolates from China, including 3 WHO reference strains(Hunan.China93-7/H1, Beijing.China94-1/H2, MVi/Menglian.Yunnan.CHN/47.09/d11), are indicated by solid rounded dots. Panel B. Phylogenetic tree of 5 measles vaccine viruses from China compared to the A genotype strains of Edmonston wild type and other vaccine strains used worldwide. Sequences from viruses isolated in China are indicated by blue, and A genotype wild type Edmonston strain is indicated by red. Two Chinese measles vaccines are indicated by green solid triangles.
Another phylogenetic tree was constructed with all H1 sequences between 1993 and 2008, which showed that all the H1 sequences could be divided into 2 groups that we propose as H1a (with 1445 strains) and H1b (with 56 strains) clusters. Based on the WHO standard 450 nucleotide sequencing window in the N gene, a 4.2% nucleotide divergence was found between the H1a and H1b clusters, and the nucleotide sequence and predicted amino acid homologies of H1a viruses were 92.3%–100% and 84.7%–100%, H1b were 97.1%–100% and 95.3%–100%, respectively. The greatest nucleotide variation among all H1 strains was 8.6% between a H1a strain, Mvi/Guangdong.PRC/19.06/5 and a H1b strain, MVi/Hebei.PRC/13.05/2.
Two clusters of measles viruses were found in China during three time periods (1993 to 1994, 1995 to 2005, and 2006 to 2008). Nucleotide sequences from 14 viruses were available between 1993 and 1994, those from 569 viruses were available between 1995 and 2005, and those from 920 viruses were available between 2006 and 2008. MeV isolates were obtained from most provinces since the year 2000 (
The geographical distribution of Chinese MeV genotypes and clusters are shown in
The provinces where the measles viruses of the indicated genotypes (clusters) were found are shown. The location within each province is not indicated. Genotype H2 viruses in Beijing may be classified as imports, and 5 genotype A viruses in four provinces are vaccine associated (see the text).
The number on the top of each column represents the number of the isolates in each time period.
Five genotype A viruses were found from 1993 to 2007: Shandong in 1993 (MVi/Shandong,CHN/9353); Hunan in 1996(MVi/Hunan.PRC/15.96/10); Henan in 1996(MVi/Henan.PRC/7.99/1); Xinjiang in 2003 (MVi/Xinjiang.PRC/24.03/1) and Shandong in 2007(MVi/Shandong.PRC/14.07/1). There was one nucleotide difference between MVi/Shandong,CHN/9353 and Edmonston-wt.USA/54/A, and one nucleotide difference between MVi/Shandong,CHN/9353 and vaccine strain AIK-C and Zagreb as well. The other four genotype A strains were genetically close to the Chinese measles vaccine strain, S191, with identical or one nucleotide difference comparing with S191.
The frequencies of synonymous (dS) and non-synonymous (dN) were calculated using the Nei-Gojobori method with MEGA Version 4.0 based on the N gene 3′ region. Eighty one H1a and nine H1b sequences shown in phylogenetic tree (
This study summarizes the measles surveillance data and documents epidemic measles cycles occurring approximately every 3 to 5 years (1992, 1997, 2001, 2005, 2008). China experienced lower incidence for 11 years during 1994–2004 after two dose measles vaccine schedule were initiated in 1986. However, this was followed by measles resurgence during 2005–2008. Thus, the 3–5 years measles epidemic cycles were possibly associated with the fluctuation of immunity coming from the combination of natural infection and vaccine-induced immunity.
China has committed to eliminate endemic MeV transmission by 2012
Among the wild-type viruses of two genotypes detected in China between 1993 and 2008 (H2 and H1), genotypes H2 had been not found since 1993–1994. H2 genotype virus, Beijing.China94-1/H2 (later designated as the WHO reference strain) was isolated from an unvaccinated 22-year-old male from Beijing, but no epidemiologic data were available to indicate the source of infection. However, it was found to be closely genetically similar to contemporaneously circulating Vietnamese measles viruses
The data presented here show that genotype H1 was the predominant circulating genotype throughout China continuously for at least 16 years. Genotype H1 has the quite great degree of intratypic variation among all genotypes, with sequences falling into two major clusters, previously designated as cluster 1 and cluster 2
A 4.2% average nucleotide divergence was found between the two clusters, cluster H1a and H1b. These two clusters of viruses were co-circulating during 1993 to 2005, while cluster H1a viruses was continuously circulating after 2006. In order to rapidly increase the immunity level of population, block the spread of the virus and reduce the incidence, the province-specific SIAs had been conducted since 2004. The transmission of H1b viruses was possibly interrupted by vaccination in 2005 and that only the H1a viruses continued to circulate. With the implementation of province-specific SIAs in 2004–2009 and following synchronized nationwide SIAs throughout of China, the reported measles incidence decreased to a historically low level of 39.5–28.9 cases per million in 2009–2010. The imported D4, D9 and d11 genotype measles viruses were found, however H1 still is predominant genotype circulating in China mainland during 2009–2010(unpublished data). More and more imported non-H1 genotype measles viruses could be detected at the measles elimination stage in China in the future.
The molecular epidemiological pattern of measles viruses in China appears to be similar to mumps virus
In this study, analysis of the predicted amino acid substitution in the 450 nt sequencing window of the MeV N gene showed no evidence of selective pressure in the past 16 years. Of course, the major target of neutralizing antibodies to measles virus in the H protein and it will be important to monitor genetic changes that may affect antigenic sites on this protein.
Although measles virological surveillance has been successfully implemented in China, and now included 30 of 31 provinces, systematic collection from each province and each prefecture should be encouraged to establish a baseline of complete chronological and geographical distribution of Chinese MeV for the entire country. Quickly identifying imported cases of measles will become more critical as measles elimination goals are achieved in China in the near future.
The list of measles virus isolates in China in 1993–2008.
(DOC)
Click here for additional data file.
We thank all the provincial and prefecture measles laboratory staffs and epidemiologists in mainland China for providing clinical specimens, isolates and epidemiologic data; We thank WHO headquarters, WPRO, US CDC for the technical and financial support. We also acknowledge anonymous reviewers for comments that improved the manuscript.