Molecular diversity of Giardia duodenalis in children under 5 years from the Manhiça district, Southern Mozambique enrolled in a matched case-control study on the aetiology of diarrhoea

Abstract

Giardia duodenalis is an enteric parasite commonly detected in children. Exposure to this organism may lead to asymptomatic or symptomatic infection. Additionally, early-life infections by this protozoan have been associated with impaired growth and cognitive function in poor resource settings. The Global Enteric Multicenter Study (GEMS) in Mozambique demonstrated that G. duodenalis was more frequent among controls than in diarrhoeal cases (≥3 loosing stools in the previous 24 hours). However, no molecular investigation was conducted to ascertain the molecular variability of the parasite. Therefore, we describe here the frequency and genetic diversity of G. duodenalis infections in children younger than five years of age with and without diarrhoea from the Manhiça district in southern Mozambique enrolled in the context of GEMS. Genomic DNA from 757 G. duodenalis-positive stool samples by immunoassay collected between 2007-2012, were reanalysed by multiplex PCR targeting the E1-HP and C1-P21 genes for the differentiation of assemblages A and B. Overall, 47% (353) of the samples were successfully amplified in at least one locus. Assemblage B accounted for 90% (319/353) of all positives, followed by assemblage A (8%, 29/353) and mixed A+B infections (1%, 5/353). No association between the presence of a given assemblage and the occurrence of diarrhoea could be demonstrated. A total of 351 samples were further analysed by a multi-locus sequence genotyping (MLSG) approach at the glutamate dehydrogenase (gdh), ß-giardin (bg) and triose phosphate isomerase (tpi) genes. Overall, 63% (222/351) of samples were genotyped and/or sub-genotyped in at least one of the three markers. Sequence analysis revealed the presence of assemblages A (10%; 23/222) and B (90%; 199/222) with high molecular diversity at the nucleotide level within the latter; no mixed infections were identified under the MLSG scheme. Assemblage A sequences were assigned to sub-assemblages AI (0.5%, 1/222), AII (7%, 15/222) or ambiguous AII/AIII (3%, 7/222). Within assemblage B, sequences were assigned to sub-assemblages BIII (13%, 28/222), BIV (14%, 31/222) and ambiguous BIII/BIV (59%, 132/222). BIII/BIV sequences accumulated the majority of the single nucleotide polymorphisms detected, particularly in the form of double peaks at chromatogram inspection. This study demonstrated that the occurrence of gastrointestinal illness (diarrhoea) was not associated to a given genotype of G. duodenalis in Mozambican children younger than five years of age. The assemblage B of the parasite was responsible for nine out of ten infections detected in this paediatric population. The extremely high genetic diversity observed within assemblage B isolates was compatible with an hyperendemic epidemiological scenario where infections and reinfections were common. The obtained molecular data may be indicative of high coinfection rates by different G. duodenalis assemblages/sub-assemblages and/or genetic recombination events, although the exact contribution of both mechanisms to the genetic diversity of the parasite remains unknown.

Fig 1. Study profile of stool samples analysed for molecular typing.

Fig 2. Distribution of single nucleotide polymorphisms segregated by mutations and double peaks, in Giardia duodenalis assemblage B sequences according to the genetic marker considered.

(A) SNPs at the glutamate dehydrogenase (gdh) locus; (B) SNPs at the triose phosphate isomerase (tpi) locus; (C) SNPs at the ß-giardin (bg) locus and overall figures for all assemblage B sequences.

Fig 3. Phylogenetic tree depicting evolutionary relationships among Giardia duodenalis sequences at the gdh locus from Mozambican children.

The analysis was inferred using the Neighbor-Joining method of the nucleotide sequence covering a 416-bp region (positions 76–491 of GenBank: L40508) of the gene. Bootstrap values lower than 50% are not shown. Cyan filled circles represent reference sequences downloaded from the GenBank database; filled dark green triangles represent sequences generated in the present study. Giardia ardeae was used as the outgroup.