# A Novel Method to Identify Routes of Hepatitis C Virus Transmission

Cyrille Féray, Julie Bouscaillou, Bruno Falissard, Mostafa Mohamed, Naglaa Arafa, Iman Bakr, Mostafa El-Hoseiny, Mai El Daly, Sherif El-Kafrawy, Sabine Plancoulaine, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Cyrille Féray, Julie Bouscaillou, Bruno Falissard, Mostafa Mohamed, Naglaa Arafa, et al. A Novel Method to Identify Routes of Hepatitis C Virus Transmission. PLoS ONE, 2014, 9 (1), pp.e86098. 10.1371/journal.pone.0086098. http://dx.doi.org/10.1371/journal.pone.0086098

HAL Id: hal-05046878

https://hal.science/hal-05046878v1

Submitted on 25 Apr 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A Novel Method to Identify Routes of Hepatitis C Virus Transmission

Cyrille Féray¹*, Julie Bouscaillou², Bruno Falissard³, Mostafa K. Mohamed⁴, Naglaa Arafa⁴, Iman Bakr⁴, Mostafa El-Hoseiny⁴, Mai El Daly⁵,⁶, Sherif El-Kafrawy⁵,⁶, Sabine Plancoulaine⁷, Mohamed Abdel-Hamid⁶,⁸, Valérie Thiers⁹, Arnaud Fontanet²,¹⁰

¹Inserm 955, Hopital Henri Mondor, Créteil, France, ²Institut Pasteur, Unité d'Epidémiologie des Maladies Emergentes, Paris, France, ³Unité INSERM 669, La maison de Solenn, Paris, France, ⁴Department of Community, Environmental and Occupational Medicine, Faculty of Medicine, Ain Shams University, Cairo, Egypt, ⁵National Liver Institute, Menoufia University, Shebeen El-Kom, Menoufia, Egypt, ⁶Viral Hepatitis Research Laboratory, National Hepatology and Tropical Medicine Research Institute, Cairo, Egypt, ⁷Unité INSERM 550, Paris, France, ⁸Department of Microbiology, Faculty of Medicine, Minia University, Minya, Egypt, ⁹Unité INSERM 785, Centre Hépatobiliaire, Villejuif, France, ¹⁰Conservatoire National des Arts et Métiers, Chaire Santé et Développement, Paris, France

**Abstract**

**Background:** We propose a new approach based on genetic distances among viral strains to infer about risk exposures and location of transmission at population level.

**Methods:** We re-analysed 133 viral sequences obtained during a cross-sectional survey of 4020 subjects living in a hepatitis C virus (HCV) endemic area in 2002. A permutation test was used to analyze the correlation between matrices of genetic distances in the NS5b region of all pairwise combinations of the 133 viral strains and exposure status (jointly exposed or not) to several potential HCV risk factors.

**Results:** Compared to subjects who did not share the same characteristics or iatrogenic exposures, the median Kimura genetic distances of viral strains were significantly smaller between brothers and sisters (0.031 versus 0.102, P<0.001), mother and child (0.044 versus 0.102, P<0.001), father and child (0.045 versus 0.102, P<0.001), or subjects exposed to periodontal treatment (0.084 versus 0.102, P=0.02). Conversely, viral strains were more divergent between subjects exposed to blood transfusions (0.216 versus 0.102, P=0.04) or tooth filling or extraction (0.108, versus 0.097, P=0.05), suggesting acquisition of the virus outside of the village.

**Conclusion:** This method provided insights on where infection took place (household, village) for several socio-demographic characteristics or iatrogenic procedures, information of great relevance for targeting prevention interventions. This method may have interesting applications for virologists and epidemiologists studying transmission networks in health-care facilities or among intravenous drug users.

**Citation:** Féray C, Bouscaillou J, Falissard B, Mohamed MK, Arafa N, et al. (2014) A Novel Method to Identify Routes of Hepatitis C Virus Transmission. PLoS ONE 9(1): e86098. doi:10.1371/journal.pone.0086098

**Editor:** Naglaa H. Shoukry, University of Montreal, Canada

**Received** July 23, 2013; **Accepted** December 4, 2013; **Published** January 23, 2014

**Copyright:** © 2014 Féray et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was supported by the French Agency for the research on HIV and hepatitis (ANRS 1211 and ANRS 12107). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: cyrille.feray@gmail.com

**Introduction**

Hepatitis C is a major health concern in Egypt [1–2] with more than 10% of the population being infected. Many studies have demonstrated that initial spread of the infection in this country was related to the massive parenteral anti-schistosomiasis treatment campaigns during the 1960s–1970s [3–5]. This was confirmed through coalescent analysis [6] which modeled the rapid spread of the epidemic from the beginning of the twentieth century. In the cohort of Zawiat Razin [4], the hepatitis C virus (HCV) prevalence in the village follows a pattern typical of Nile Delta villages in Egypt with the prevalence of HCV antibodies increasing with age, from 2.7% in those <20 years, to reach more than 30% for those above 45 years. Using oral questionnaire, HCV infections could be attributed to iatrogenic factors in subjects

older than 20 years, but not in younger subjects. We recently published an analysis of the same data taking into account the clustering of infections by families [7]. Based on HCV antibody status, there was highly significant evidence for intra-familial clustering of HCV infections, especially among father-offspring, mother-offspring, and sib-sib pairs.

However, how much of this clustering is related to shared at risk behaviours, predisposing genetic factors [8] or intra-household transmission could not be studied based on serological data alone.

To overcome this limitation, we developed a new approach for analyzing viral sequencing data at the household and community level. The basic hypothesis is that highly similar strains reflect either a direct transmission event or exposure to a shared source, whereas divergent strains suggest separate acquisition of HCV in time and place. To explore this hypothesis, we computed
correlation coefficients relating genetic distances and exposure to a given risk factor for all pairwise combinations of subjects with viral sequences available. We used a permutation method for significance testing, the Mantel's test [9] as it does not require independence among observations.

# Materials and Methods

## Population

A large epidemiological study was conducted between May and November 2002 in Zawiat Razin, a village of around 20,000 inhabitants in Menofia Governorate (Nile Delta). Practical details have been described elsewhere [4]. After informed consent was obtained (from the head of household for children less than 18 years old), all residents older than three years of age and living in one sector of the village (25% of the total village population) were invited to participate. A questionnaire was administered on socio-demographic characteristics (age, sex, duration of residence in the village, marital status, educational level) and iatrogenic exposures (hospitalisation, injections, dentistry, surgical and obstetrical procedures, blood transfusions) possibly associated with the transmission of HCV in the village. A second-round follow-up of the cohort was conducted between June 2003 and March 2004 [7]. First-degree relationships (sib-sib, mother-child, father-child), and spouse relationships living in the same household were recorded. Other persons living in the same household were considered as not family-related for the purpose of this study. The starting population was 4020 individuals. As previously reported [4], 475 (11.8%) were positive for anti-HCV antibodies, 456 could be tested for HCV RNA and 273/456 (60%) were serum HCV RNA positive. Because this study was interested in household transmission, we selected for sequencing viremic patients only when someone else in the household was viremic leading to a subgroup of 133 subjects.

All participants (or head of the household for children aged less than 18 years old) signed an informed consent. The study protocol was reviewed and approved by the committee for biomedical research at Institut Pasteur, Paris, and by the Institutional Review Board of the National Hepatology & Tropical Medicine Research Institute, Cairo.

## Methods

**Amplification and sequencing.** Partial amplification of the NS5B (nucleotide position 7915–8303 330 bp) was performed as described elsewhere [10]. Molecular characterization was first carried out by sequence analysis and phylogenetic studies of HCV NS5B region. A partial amplification of E1 region (nucleotide position 482–1320, 830 bp was also performed as described previously [11] in a subset of positive samples. Positive samples were purified, and bidirectional sequence analysis was performed. Consensus nucleotide sequences were used for phylogenetic studies.

**Phylogenetic analysis.** The NS5B or E1 sequences were added to a set of representative HCV genotypes and subtypes for the same gene region available in the databanks (n>100). The DNA alignments were generated with Clustal W and followed by further phylogenetic analysis using the Neighbor Joining program using 1,000 bootstrap replication with in the Mega package [12]. Positions for which an ambiguity was observed were considered as missing for the calculation of nucleotidic distances.

**Correlation matrix analysis.** Traditional statistics cannot be applied since genetic distances are pairwise defined and these distances are thus not independent from each other. Mantel's test [9] is a generalized regression permutation procedure classically

used to compare two distance matrices. In our approach we correlated the genetic distance matrix consisting of pairwise Kimura two-parameter distances (d) with phenotypic matrix. They were calculated between each pair of subjects for a set of given risk factors and arranged in dissimilarity matrices where the distance was 0 if both subjects were exposed and otherwise 1 (Figure S1 in supplemental data). The studied variables were those listed in Table 1.

The correlation between each distance matrix and each phenotypic matrix was evaluated using the Spearman correlation coefficient (R0), which ranges from -1.0 for a perfect negative to 1.0 for a perfect positive correlation between two matrices. As variables arranged into matrices are not independent (e.g. the distance between case 1 and 3 is not independent of the distance between case 1 and 2, because case 1 is involved in both), the significance of the correlation is determined by a permutation test. The rows and columns of one matrix were randomly permuted 5000 times and the Spearman correlation was calculated for each permutation. The measure of significance is given by the ratio N/ 5000, where N is the number of times that R0 is exceeded by correlation coefficients calculated with permuted matrices. If the original matrices are correlated, the disruption caused by the permutations should produce correlation coefficients below R0.

Analyses were conducted using STATA 11.0 software (Stata Corporation, College Station, Texas, USA).

## Results

### Study Population

The mean age of the subgroup with sequencing data was 36 years (±15; range 9–71 years) with 77 (58%) subjects less than 41 years. Seventy-eight (59%) were male. 122 (92%) were born in the village. These 133 subjects lived in 54 different houses; 38 of them had no familial relationship (spouse, child, parent, brother/sister) with any of the other included subjects. Among the others, we could analyze viral distances in 14 married couples, in 22 mother-child pairs, in 15 father-child pairs and 28 sib-sib pairs (Note: The same individual can belong to more than one pair, e.g., a child

<table><caption>Table 1. Number of participants exposed to each risk factor.</caption><thead><tr><td>Possible risk factor</td><td>n</td></tr></thead><tbody><tr><td>Born in the village</td><td>122</td></tr><tr><td>Infected brother/sister</td><td>41</td></tr><tr><td>Infected mother</td><td>22</td></tr><tr><td>Infected father</td><td>15</td></tr><tr><td>Infected spouse</td><td>28</td></tr><tr><td>Infected non familial household member</td><td>62</td></tr><tr><td>Surgery</td><td>51</td></tr><tr><td>Stitches</td><td>64</td></tr><tr><td>Abcess drainage</td><td>33</td></tr><tr><td>Intravenous injections</td><td>18</td></tr><tr><td>Anti-schistomiasis treatment</td><td>21</td></tr><tr><td>Periodontal treatment</td><td>21</td></tr><tr><td>Tooth filling or extraction</td><td>68</td></tr><tr><td>Urinary catheter</td><td>11</td></tr><tr><td>Blood transfusion</td><td>8</td></tr><tr><td>Instrumental delivery</td><td>7</td></tr></tbody></table>

doi:10.1371/journal.pone.0086098.t001
<figure><img src="image_3.png" /><figcaption>**Figure 1. Phylogenetic tree from the HCV sequences of 105 patients infected by the genotype 4a.** Same shape filled with same color means same house and same family. Same shape empty means same house, but not the same family. Absence of symbol means no connection to other participant with sequence.</figcaption></figure>

doi:10.1371/journal.pone.0086098.g001

may be in a mother-child, father-child, and sib-sib pairs simultaneously).

Other risk factors for this population are described in Table 1. We only considered risk factors for which the percentage of exposed ranged between 5% and 95%, so that genetic distances estimates for each group (jointly exposed or not) could be based on sufficient numbers. Of note, many patients were exposed to more than one risk factor.

## Phylogenetic Analysis

Comparison with published reference sequences from Genbank showed that 130 (98%) of the Egyptian sequences clustered with genotype 4, and that 3 (2%) clustered with genotype 1. Most cases (105 out of 133; 77%) clustered with the subtype 4a Egyptian

reference sequence ED43. The remaining sequences formed 5 distinct groups with strong bootstrap support. Two of these groups clustered with previously described Egyptian sequences identified as subtypes 4m and 4o (n = 2 and 3, respectively), and one with an undescribed yet subtype (provisionally called subtype 4u; n = 20). The last group (n = 3) clustered as a monophyletic group (bootstrap value 99%) corresponding to genotype 1g (n = 2) and 1b (n = 1). The NS5B phylogenetic classification of the Egyptian isolates was confirmed by partial amplification of the E1 region (position 482–1320) in 66 subjects. Altogether these data indicate the presence of a new subtype of type 4. Sequence access numbers are EF694393–EF694525 for the NS5B region and EF694321–EF694386 for the E1 region.
<table><caption>**Table 2.** Pairwise NS5b HCV genetic distances among 133 subjects according to common exposure patterns.</caption><thead><tr><th rowspan="2"></th><th>Yes</th><th>No</th><th rowspan="2">rho</th><th rowspan="2">p</th></tr><tr><th>1<sup>st</sup> quart. - Median - 3<sup>rd</sup> quart.</th><th>1<sup>st</sup> quart. - Median - 3<sup>rd</sup> quart.</th></tr></thead><tbody><tr><td colspan="5"><strong>Household/Family ties</strong></td></tr><tr><td>Same house, no family tie</td><td>0.049 - 0.078 - 0.112</td><td>0.081 - 0.102 - 0.214</td><td>0.04</td><td>&lt;0.001</td></tr><tr><td>Sib-sib</td><td>0.019 - 0.031 - 0.052</td><td>0.081 - 0.102 - 0.214</td><td>0.07</td><td>&lt;0.001</td></tr><tr><td>Mother-child</td><td>0.026 - 0.044 - 0.096</td><td>0.081 - 0.102 - 0.214</td><td>0.05</td><td>&lt;0.001</td></tr><tr><td>Father-child</td><td>0.032 - 0.045 - 0.052</td><td>0.081 - 0.102 - 0.214</td><td>0.05</td><td>&lt;0.001</td></tr><tr><td>Spouses</td><td>0.032 - 0.061 - 0.209</td><td>0.081 - 0.102 - 0.214</td><td>0.02</td><td>0.03</td></tr><tr><td>Entire life in the village</td><td>0.081 - 0.101 - 0.212</td><td>0.082 - 0.154 - 0.221</td><td>0.05</td><td>0.21</td></tr><tr><td colspan="5"><strong>Iatrogenic exposures</strong></td></tr><tr><td>Surgery</td><td>0.083 - 0.105 - 0.220</td><td>0.081 - 0.102 - 0.213</td><td>-0.04</td><td>0.21</td></tr><tr><td>Stitches</td><td>0.083 - 0.110 - 0.229</td><td>0.080 - 0.101 - 0.206</td><td>-0.09</td><td>0.05</td></tr><tr><td>Blood transfusion</td><td>0.114 - 0.216 - 0.467</td><td>0.081 - 0.102 - 0.214</td><td>-0.04</td><td>0.04</td></tr><tr><td>Bladder catheterization</td><td>0.094 - 0.106 - 0.122</td><td>0.081 - 0.102 - 0.214</td><td>-0</td><td>0.46</td></tr><tr><td>IV cannula</td><td>0.083 - 0.101 - 0.205</td><td>0.081 - 0.102 - 0.217</td><td>0.01</td><td>0.38</td></tr><tr><td>10+ injections</td><td>0.080 - 0.099 - 0.207</td><td>0.083 - 0.106 - 0.222</td><td>0.06</td><td>0.16</td></tr><tr><td>Treatment for abcess</td><td>0.077 - 0.095 - 0.211</td><td>0.083 - 0.102 - 0.214</td><td>0.03</td><td>0.24</td></tr><tr><td>Periodontal treatment</td><td>0.073 - 0.084 - 0.102</td><td>0.083 - 0.102 - 0.215</td><td>0.07</td><td>0.02</td></tr><tr><td>Tooth filling or extraction</td><td>0.084 - 0.108 - 0.219</td><td>0.079 - 0.097 - 0.208</td><td>-0.11</td><td>0.05</td></tr><tr><td>Anti-schistomiasis treatment</td><td>0.080 - 0.113 - 0.195</td><td>0.081 - 0.102 - 0.215</td><td>0.01</td><td>0.36</td></tr><tr><td colspan="5"><strong>Obstetrical history*</strong></td></tr><tr><td>5+ delivery</td><td>0.087 - 0.106 - 0.233</td><td>0.083 - 0.102 - 0.218</td><td>-0.05</td><td>0.25</td></tr><tr><td>Traumatic delivery</td><td>0.076 - 0.109 - 0.113</td><td>0.083 - 0.102 - 0.218</td><td>0.03</td><td>0.29</td></tr><tr><td>Abortion</td><td>0.094 - 0.207 - 0.237</td><td>0.083 - 0.101 - 0.216</td><td>-0.11</td><td>0.07</td></tr></tbody></table>

doi:10.1371/journal.pone.0086098.t002

Figure 1 shows the phylogenetic tree corresponding to the 105 subjects infected by HCV genotype 4a. The proximity of strains according to families and houses is suggested although bootstrap value for each familial or household cluster was always below 50%.

## Matrix Correlation

From the 133 subjects, we constructed a genetic matrix of 8778 combinations of genetic distances defined in the NS5B region and 18 corresponding phenotypic matrices of the same size.

The genetic (Kimura) distances between the strains infecting two subjects having a sib-sib, mother-child, father-child or spouse relationship, living in the same house without familial relationship, or submitted to periodontal treatment were lower than the distances between the strains infecting two subjects who did not share the same relationships (Table 2).

Conversely, the genetic distances between strains infecting two subjects who have been exposed to stitching, blood transfusion, or tooth filling or extraction were significantly higher than in those not sharing these risk factors (the difference was marginally significant for abortion).

We did not evidence a relationship between genetic distance and a common exposure to antischistomiasis therapy, injections, instrumental deliveries or the other studied risk factors.

The distribution of viral genetic distances according to the risk factors is shown in the figures 2 to 4. Among non-exposed (continuous line), the three peaks correspond to genetic distance among individuals belonging to the same subtype (left), to the same genotype, but not the same subtype (middle), and to different genotypes (right). The same applies to the peaks observed among

exposed (dashed lines), but low numbers may affect the visibility of the peaks. Genetic relatedness of strains translated into shifting to the left of the distribution due to lower distances in the first peak (within subtype) among exposed compared to non-exposed (e.g., siblings). Genetic remoteness of the strains translated into shifting to the right of the distribution due to a higher right peak following the introduction of a new genotype, genotype 1, in the village (e.g., blood transfusion).

The same analysis was repeated using genetic distances derived from the E1 region of 66 subjects. There was a high correlation between the pairwise genetic distances calculated from NS5b and E1 regions (Fig. S1). The risk factor analysis using the E1 region pairwise distances found the same results as with the NS5b region, except for non-related household members and spouses, which no longer showed the same level of closeness, and for stitches which were no longer more distant (Table 3). The same analysis was performed in the group of 105 subjects with infection by genotype 4a and the same relationships were evidenced in this group (data not shown).

## Discussion

The proposed approach is an original application of the matrix correlation test to the issue of viral transmission among a limited population of infected subjects exposed to multiple risk factors. It allows identification of risk factors that would, in a given community, lead to either homogeneous strains (when transmission occurs among members of that community or through exposure to a common source of infection), or divergent strains
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Distributions of viral genetic distances depending on household/family ties (n = 133). Among non-exposed (continuous line), the three peaks correspond to genetic distance among individuals belonging to the same subtype (left), to the same genotype, but not the same subtype (middle), and to different genotypes (right).
doi:10.1371/journal.pone.0086098.g002</figcaption></figure>

(when infections are contracted outside of the community or at different times). We demonstrated that two relatives living in the same house were infected by less divergent strains than two subjects living in different houses or belonging to different families. Also, two subjects having been exposed to periodontal treatment were infected by more similar strains compared to others. Conversely, two subjects exposed to blood transfusions, or tooth filling or extraction had more divergent strains than others.

The approach was particularly powerful to document intrafamilial transmission of HCV (i.e., transmission of HCV among two members of the same family). Such demonstration has proven to be difficult in the past with classical epidemiological studies (reviewed in [13]). Many studies relied on clustering of HCV infections based on positive serology only. In such circumstances, one cannot exclude common at risk behaviours shared among members of the same family. More convincing are studies using sequences homologies to document intrahousehold transmission. Such studies were few, and with limited statistical power. In that same village, it took us four years of follow-up to identify 25 seroconverters among 3184 subjects. Of these, 17 were viremic, and 2 only had another family member infected with the same strain indicating intrafamilial transmission [14]. With the new approach described in this paper, we obtained similar evidence for intrahousehold transmission with a much smaller number of individuals and a cross-sectional survey only. Our approach, based on sampling of households with multiple cases, is quick to realize (cross-sectional) and pose no methodological challenges (no need for complex strategies in choosing the controls of case-control studies).

HCV is classically considered as poorly contagious within families or between spouses [15–19]. In western countries, clear parenteral routes (intravenous drug usage and transfusions before 1991), are the main routes of transmission. However, subjects with no identifiable route of infection represent 10% to 20% of all cases in most Western countries series [20–21] and were frequent in the studied Egyptian subjects, particularly among children [4]. Our

results indicate that an important source of contamination comes from the household and is particularly strong among siblings and between parents and children pointing these relationships as targets for screening and prevention. Further behavioural studies are needed to clearly identify the route of infection which may involve either classical modes of transmission (e.g., shared medical injections) or minute amounts of blood, or other body fluids such as saliva. Transmission between spouses was likely in this study and others performed in other rural communities of Egypt [22], but remains an issue of debate [17–19].

Other risk factors identified through this analysis were periodontal treatment, tooth filling or extraction, and blood transfusion. All, except for tooth filling or extraction, were risk factors for HCV infection in previous reports from the same community [4,7], with the additional information provided here that periodontal treatment was most likely performed within the village since subjects reporting this practice had more similar strains, while tooth filling or extraction, and blood transfusions were performed outside of the village since subjects reporting these practices had more divergent strains. This information is particularly useful for public health practitioners when trying to identify where infection takes place in order to propose corrective measures. In this context, it would be useful to examine what type of periodontal treatment is performed within the village, and how the material is cleaned between patients. Local practices include cleaning of gums by non medical staff with pieces of metal for aesthetic purposes, and may be associated with HCV transmission if not properly sterilised between individuals. Of note, periodontal treatment has also been reported associated with HCV infection in a case-control study analyzing incident infection in Cairo [23]. Bleeding during periodontal treatment is frequent. A study reported the detection of HCV RNA on instruments after dental treatment in HCV infected patients [24].

This approach should not be used as a replacement of classical epidemiological studies, such as case-control or cohort studies, for the identification of risk factors. Indeed, it will not identify
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Distributions of viral genetic distances depending on obstetric risk factors (n = 133). Same interpretation as for Figure 2.
doi:10.1371/journal.pone.0086098.g003</figcaption></figure>

procedures taking place both within and outside of the commu-
nity, as shown here by the lack of significant relationships for risk
factors such as treatment for schistosomiasis, intravenous injec-
tions, or obstetrical risks. The fact that history of treatment for
schistosomiasis did not appear as a risk factor in our analysis may
also be related to the very old timing of this exposure which took
place 30 to 40 years prior to this survey. Although the strains that
spread during the mass treatment campaigns might have been
closely related initially, they had time to diverge since, and cannot
be recognized during this analysis. Conversely, the approach may

be biased towards emphasizing intrafamilial transmission if strains
infecting two siblings experience lower divergence over time due to
similar HLA-driven immune response targeting epitopic sites.

This same Mantel test has been successfully used in different
fields with the purpose of comparing matrix variables [25–26]. It
has also been applied in the field of virology for the study of
cellular tropism of HIV [27–28] or HCV quasispecies [29–30] or
for the study of relations between variability of genotype 1b and
the severity of fibrosis after liver transplantation [31]. It has been
used to explain spatial patterns in genetic distances among

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Distributions of viral genetic distances depending on obstetric risk factors (n = 55). Same interpretation as for Figure 2.
doi:10.1371/journal.pone.0086098.g004</figcaption></figure>
<table><caption>**Table 3.** Pairwise E1 HCV genetic distances among 66 subjects according to common exposure patterns.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Yes</th><th colspan="2">No</th></tr><tr><th>1<sup>st</sup> quart. - Median - 3<sup>rd</sup> quart.</th><th>1<sup>st</sup> quart. - Median - 3<sup>rd</sup> quart.</th><th>rho</th><th>p</th></tr></thead><tbody><tr><td colspan="5"><strong>Household/Family ties</strong></td></tr><tr><td>Same house, no family tie</td><td>0,080 - 0,113 - 0,195</td><td>0,086 - 0,110 - 0,227</td><td>0,03</td><td>0,01</td></tr><tr><td>Sib-sib</td><td>0,024 - 0,049 - 0,079</td><td>0,086 - 0,110 - 0,227</td><td>0,06</td><td>0,001</td></tr><tr><td>Mother-child</td><td>0,026 - 0,044 - 0,210</td><td>0,086 - 0,110 - 0,227</td><td>0,02</td><td>0,05</td></tr><tr><td>Father-child</td><td>0,032 - 0,061 - 0,219</td><td>0,086 - 0,110 - 0,227</td><td>0,05</td><td>0,01</td></tr><tr><td>Spouses</td><td>0,083 - 0,105 - 0,220</td><td>0,086 - 0,110 - 0,227</td><td>0,02</td><td>0,10</td></tr><tr><td>Entire life in the village</td><td>0,091 - 0,170 - 0,232</td><td>0,112 - 0,190 - 0,242</td><td>0,05</td><td>0,31</td></tr><tr><td colspan="5"><strong>Iatrogenic exposures</strong></td></tr><tr><td>Surgery</td><td>0,083 - 0,105 - 0,260</td><td>0,083 - 0,120 - 0,224</td><td>-0,04</td><td>0,25</td></tr><tr><td>Stitches</td><td>0,073 - 0,110 - 0,227</td><td>0,084 - 0,120 - 0,222</td><td>-0,1</td><td>0,1</td></tr><tr><td>Blood transfusion</td><td>0,114 - 0,216 - 0,467</td><td>0,079 - 0,123 - 0,214</td><td>-0,04</td><td>0,05</td></tr><tr><td>Bladder catheterization</td><td>0,083 - 0,105 - 0,220</td><td>0,086 - 0,110 - 0,227</td><td>0</td><td>0,46</td></tr><tr><td>IV cannula</td><td>0,087 - 0,110 - 0,233</td><td>0,086 - 0,110 - 0,227</td><td>0,009</td><td>0,58</td></tr><tr><td>10+ injections</td><td>0,080 - 0,106 - 0,259</td><td>0,086 - 0,110 - 0,227</td><td>0</td><td>0,21</td></tr><tr><td>Treatment for abcess</td><td>0,081 - 0,106 - 0,253</td><td>0,086 - 0,110 - 0,227</td><td>0,02</td><td>0,34</td></tr><tr><td>Periodontal treatment</td><td>0,080 - 0,102 - 0,228</td><td>0,086 - 0,110 - 0,227</td><td>0,005</td><td>0,09</td></tr><tr><td>Tooth filling or extraction</td><td>0,084 - 0,108 - 0,219</td><td>0,071 - 0,090 - 0,217</td><td>-0,15</td><td>0,07</td></tr><tr><td>Anti-schistomiasis treatment</td><td>0,087 - 0,103 - 0,233</td><td>0,081 - 0,102 - 0,215</td><td>0,02</td><td>0,36</td></tr><tr><td colspan="5"><strong>Obstetrical history*</strong></td></tr><tr><td>5+ delivery</td><td>0,083 - 0,156 - 0,267</td><td>0,087 - 0,106 - 0,243</td><td>-0,06</td><td>0,34</td></tr><tr><td>Traumatic delivery</td><td>0,061 - 0,102 - 0,209</td><td>0,087 - 0,106 - 0,239</td><td>0,02</td><td>0,37</td></tr><tr><td>Abortion</td><td>0,091 - 0,212 - 0,243</td><td>0,087 - 0,105 - 0,233</td><td>-0,1</td><td>0,17</td></tr></tbody></table>

doi:10.1371/journal.pone.0086098.t003

different populations of human, animals [32] or bacteria. However it has never been applied to the issue of viral transmission.

In conclusion, this study has shown an interesting application of a new statistical approach to the analysis of viral genetic distances to better understand transmission routes of infection. It has clearly documented the distinct participation of intrafamilial and intrahousehold components in the transmission of HCV in a rural area of Egypt. It has also pointed at periodontal treatment as a potential source of HCV transmission within the village. Other applications of this technique may be envisaged for the investigation of transmission networks in health care facilities or among intravenous drug users.

## References

1. Kamel MA, Ghaffar YA, Wasef MA, Wright M, Clark LC, et al. (1992) High HCV prevalence in Egyptian blood donors. Lancet. 340(8816): 427.

2. Abdel-Aziz F, Habib M, Mohamed MK, Abdel-Hamid M, Gamil F, et al. (2000) Hepatitis C virus (HCV) infection in a community in the Nile Delta: population description and HCV prevalence. Hepatology (Baltimore, Md). 32(1): 111–5.

3. Frank C, Mohamed MK, Strickland GT, Lavanchy D, Arthur RR, et al. (2000) The role of parenteral antischistosomal therapy in the spread of hepatitis C virus in Egypt. Lancet. 355(9207): 887–91.

4. Arafa N, El Hoseiny M, Rekacewicz C, Bakr I, El-Kafrawy S, et al. (2005) Changing pattern of hepatitis C virus spread in rural areas of Egypt. Journal of hepatology. 43(3): 418–24.

5. Mohamed MK, Abdel-Hamid M, Mikhail NN, Abdel-Aziz F, Medhat A, et al. (2005) Intrafamilial transmission of hepatitis C in Egypt. Hepatology (Baltimore, Md. 42(3): 683–7.

## Supporting Information

**Figure S1** Scatterplot of correlation between the pairwise genetic distances derived from NS5b and E1 regions in a subsample of 66 participants. (TIF)

## Acknowledgments

We thank Pierre Legendre for helpful discussions.

## Author Contributions

Conceived and designed the experiments: CF MKM AF. Performed the experiments: NA IB MED SEK MAH VT SP. Analyzed the data: CF JB BF MKM SP MAH VT AF MEH. Contributed reagents/materials/ analysis tools: MED SEK MAH VT. Wrote the paper: CF JB AF.

6. Pybus OG, Drummond AJ, Nakano T, Robertson BH, Rambaut A (2003) The epidemiology and iatrogenic transmission of hepatitis C virus in Egypt: a Bayesian coalescent approach. Molecular biology and evolution. 20(3): 381–7.

7. Plancoulaine S, Mohamed MK, Arafa N, Bakr I, Rekacewicz C, et al. (2008) Dissection of familial correlations in hepatitis C virus (HCV) seroprevalence suggests intrafamilial viral transmission and genetic predisposition to infection. Gut. 57(9): 1268–74.

8. Laouenan C, Plancoulaine S, Mohamed MK, Arafa N, Bakr I, et al. (2009) Evidence for a dominant major gene conferring predisposition to hepatitis C virus infection in endemic conditions. Hum Genet. 126(5): 697–705.

9. Mantel N (1967) The detection of disease clustering and a generalized regression approach. Cancer research. 27(2): 209–20.

10. Delarocque-Astagneau E, Baffoy N, Thiers V, Simon N, de Valk H, et al. (2002) Outbreak of hepatitis C virus infection in a hemodialysis unit: potential
transmission by the hemodialysis machine? Infect Control Hosp Epidemiol. 23(6): 328-34.

11. McKeating JA, Zhang LQ, Logvinoff C, Flint M, Zhang J, et al. (2004) Diverse hepatitis C virus glycoproteins mediate viral infection in a CD81-dependent manner. Journal of virology. 78(16): 8496-505.

12. Kumar S, Tamura K, Nei M (2004) MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 5(2): 150-63.

13. Ackerman Z, Ackerman E, Paltiel O (2000) Intrafamilial transmission of hepatitis C virus: a systematic review. Journal of viral hepatitis. 7(2): 93-103.

14. Mostafa A, Taylor SM, El-Daly M, El Hoseiny M, Bakr I, et al. (2010) Is the hepatitis C virus epidemic over in Egypt? Incidence and risk factors of new hepatitis C virus infections. Liver Int. 30: 560-6.

15. Diago M, Zapater R, Tuset C, Carbonell P, Gonzalez C, et al. (1996) Intrafamily transmission of hepatitis C virus: sexual and non-sexual contacts. Journal of hepatology. 25(2): 125-8.

16. Meisel H, Reip A, Faltus B, Lu M, Porst H, et al. (1995) Transmission of hepatitis C virus to children and husbands by women infected with contaminated anti-D immunoglobulin. Lancet. 345(8959): 1209-11.

17. Dienstag JL (1997) Sexual and perinatal transmission of hepatitis C. Hepatology (Baltimore, Md). 26(3 Suppl 1): 66S-70S.

18. Stroffolini T, Lorenzoni U, Menniti-Ippolito F, Infantolino D, Chiaramonte M (2001) Hepatitis C virus infection in spouses: sexual transmission or common exposure to the same risk factors? The American journal of gastroenterology. 96(11): 3138-41.

19. Tahan V, Karaca C, Yildirim B, Bozbas A, Ozaras R, et al. (2005) Sexual transmission of HCV between spouses. The American journal of gastroenterology. 100(4): 821-4.

20. Zeuzem S, Teuber G, Lee JH, Ruster B, Roth WK (1996) Risk factors for the transmission of hepatitis C. Journal of hepatology. 24(2 Suppl): 3-10.

21. Murphy EL, Bryzman SM, Glynn SA, Ameti DI, Thomson RA, et al. (2000) Risk factors for hepatitis C virus infection in United States blood donors. NHLBI Retrovirus Epidemiology Donor Study (REDS). Hepatology (Baltimore, Md). 31(3): 756-62.

22. Magder LS, Fix AD, Mikhail NN, Mohamed MK, Abdel-Hamid M, et al. (2005) Estimation of the risk of transmission of hepatitis C between spouses in Egypt based on seroprevalence data. International journal of epidemiology. 34(1): 160-5.

23. Paez Jimenez A, Sharaf Eldin N, Rimlinger F, El-Daly M, El-Hariri H, et al. (2010) HCV iatrogenic and intrafamilial transmission in Greater Cairo, Egypt. Gut. 59(11): 1554-60.

24. Piazza M, Borgia G, Picciotto L, Nappa S, Cicciarello S, et al. (1995) Detection of hepatitis C virus-RNA by polymerase chain reaction in dental surgeries. Journal of medical virology. 45(1): 40-2.

25. Sokal RR, Oden NL, Thomson BA (1992) Origins of the Indo-Europeans: genetic evidence. Proceedings of the National Academy of Sciences of the United States of America. 89(16): 7669-73.

26. Waddle DM (1994) Matrix correlation tests support a single origin for modern humans. Nature. 368(6470): 452-4.

27. Poss M, Rodrigo AG, Gosink JJ, Learn GH, de Vange Panteleeff D, et al. (1998) Evolution of envelope sequences from the genital tract and peripheral blood of women infected with clade A human immunodeficiency virus type 1. Journal of virology. 72(10): 8240-51.

28. Collins KR, Quinones-Mateu ME, Wu M, Luzze H, Johnson JL, et al. (2002) Human immunodeficiency virus type 1 (HIV-1) quasispecies at the sites of Mycobacterium tuberculosis infection contribute to systemic HIV-1 heterogeneity. Journal of virology. 76(4): 1697-706.

29. Roque Afonso AM, Jiang J, Penin F, Tareau C, Samuel D, et al. (1999) Nonrandom distribution of hepatitis C virus quasispecies in plasma and peripheral blood mononuclear cell subsets. Journal of virology. 73(11): 9213-21.

30. Ducoulombier D, Roque-Afonso AM, Di Liberto G, Penin F, Kara R, et al. (2004) Frequent compartmentalization of hepatitis C virus variants in circulating B cells and monocytes. Hepatology (Baltimore, Md). 39(3): 817-25.

31. Gigou M, Roque-Afonso AM, Falissard B, Penin F, Dussaix E, et al. (2001) Genetic clustering of hepatitis C virus strains and severity of recurrent hepatitis after liver transplantation. Journal of virology. 75(23): 11292-7.

32. Legendre P, Lapointe FJ (1995) Matching behavioral evolution to brain morphology. Brain, behavior and evolution. 45(2): 110-21.