A novel species of torque teno mini virus (TTMV) in gingival tissue from chronic periodontitis patients.

A new species of torque teno mini virus, named TTMV-222, was detected in gingival tissue from periodontitis patients using a viral metagenomics method. The 2803-nucleotide genome of TTMV-222 is closely related to TTMV1-CBD279, with 62.6% overall nucleotide similarity. Genetic analyses of the new virus genome revealed a classic genomic organization but a weak identity with known sequences. The prevalence of TTMV-222 in the periodontitis group (n = 150) was significantly higher than that in the healthy group (n = 150) (p = 0.032), suggesting that the new virus may be associated with inflammation in chronic periodontitis patients. However, this finding requires further investigation.

Periodontitis is a common, multifactorial, chronic oral disease that can lead to the early loss of teeth. While the incidence of periodontal disease is high, the etiology remains unclear. Multiple infectious agents are believed to play a role in the pathogenesis of this disease. Bacteria, such as Porphyromonas gingivalis, Tannerella forsythia, Aggregatibacter actinomycetemcomitans, and Treponema denticola, are thought to be major causes of periodontal disease12345, but the bacterial theory cannot fully explain all of the clinical symptoms associated with periodontal disease678. Despite the presence of abundant periodontopathic bacteria in the oral cavity, periodontitis is a site-specific disease. The disease course includes remission periods that can be interrupted by episodes of clinical relapses9. The latent biological basis of this phenomenon is not clearly understood. The possible role of viruses in periodontal disease episodes has begun to attract more attention10. Viruses have been implicated in the pathogenesis of periodontal diseases since the mid-1990s1112131415. Previous data have shown that individual periodontal lesions harbor millions of genomic copies of herpes viruses, papilloma viruses, human T-lymphotropic virus-111, hepatitis B virus12, hepatitis C virus13, Human Immunodeficiency Virus (HIV) and torque teno virus (TTV)14. Chikungunya virus was also reported to be associated with Indian gingivitis patients15. Recent studies have discussed the possible role of bacteriophages in periodontal disease patients, suggesting that oral phages may be more highly expressed in periodontally healthy subjects and that periodontitis may be beneficial for the expression of some lytic phages116. All of these studies have provided new insight into the possible cause of this disease.

Viral metagenomics is able to characterize viral communities and to detect previously uncharacterized viruses based on sequence similarities to previously sequenced viral genomes17181920. This strategy has been used to discover a wide diversity of viruses in both marine and freshwater environments and in human and animal samples212223242526272829. In the present study, we characterized a novel human anellovirus from chronic periodontitis patients in China using viral metagenomics technology. Epidemiological strategies were adopted to explore the correlation between the newly identified virus and chronic periodontitis.

Results

Discovery of a novel TTMV

A total of 48 chronic periodontitis subjects (24 men and 24 women) with a mean age of 42.2 years of age were included in this study. All of them received a periodontal examination at the Ninth People’s Hospital, School of Medicine, Shanghai Jiao Tong University in China. Each diagnosis was based on clinical and radiographic findings according to the inclusion standards described in the Methods. The mean pocket depth (PD), clinical attachment loss (CAL), gingival index (GI), and plaque index (PLI) of all of the participants were 8.55 mm, 4.52 mm, 2.27, and 2.58, respectively. Samples were taken from each of the patients. The viral particles and their nucleic acids in these samples were enriched by filtration and nuclease treatment before the nucleic acids were extracted for use in random RT-PCR. The samples were then sequenced using the MiSeq platform. Based on the best BLAST score (E < 0.001), a novel viral sequence (222 bp), identified as an anellovirus, was further analyzed. This sequence had a 39% amino acid homology with a known human anellovirus.

Full-length nucleotide sequences of the newly identified virus

Based on the novel viral 222-bp sequencent, primers (AL1 5′-A*C*TCGTGTGTCTCCTCCTCG-3′, AR1 5′-C*G*GAGACGGCATCACATCAG-3′, AR2 5′-TGATACCCGGCTGATCTTGG-3′) were designed for amplification by inverse PCR and Sanger sequencing. The complete genome of the new isolate was 2803 bp in length and contained three open reading frames (ORFs) (Fig. 1). ORF1 is 1941 bp (nt316–2256) in length and encodes a 646-amino acid protein. ORF2 is 267 bp (nt170–436) in length and encodes an 88-amino acid protein. This ORF partially overlaps with ORF1 (nt316–436). ORF3 is 273 bp (nt1979–2251) in length, encodes a 90-amino acid protein, and also overlaps with ORF1 (nt1979–2251). These characteristics are consistent with the composition of the human TTMV strains that have been reported3031. The viral genome of the new isolate was deposited into GenBank (accession number: KU041847) and was named TTMV-222.

Figure 1

The genome organization of TTMV-222.

Annotations and illustritions were made using Vector NTI 10.

Sequence analysis of TTMV-222

A phylogenetic tree (Fig. 2) was constructed with ORF1 from TTMV-222 and TTMV reference species31. TTMV-222 is related to TTMV1-CBD279, TTMV2-NLC023, and TTMV3-NLC026, all of which belong to the family torque teno mini virus (TTMV), genus Betatorquevirus (Fig. 2). Homology analysis revealed that TTMV-222 has 62.6% and 55.8% nucleotide homologies with the TTMV1-CBD279 (GenBank accession No. AB026931) and TTMV2-NLC023 strains (GenBank accession No. AB038629), respectively (Table 1). In ORF1, TTMV-222 shares the highest nucleotide homology (60.4%) with the TTMV1-CBD279 strain; TTMV-222 ORF1 also has 47.3% amino acid homology with this strain (Table 1). The nucleotide and amino acid homologies of TTMV-222 with reference strains of TTMV for ORF1, ORF2, ORF3, and the complete genome are listed in Table 1. The nucleotide homologies between different TTMV species for ORF1 are shown in Table S1 in the supplemental material.

Figure 2

ORF1 nucleotide phylogenetic tree constructed using the neighbor-joining method with nucleotide p distances and 1,000 bootstrap replicates in the Molecular Evolutionary Genetics Analysis program (MEGA, version 4.0, USA).

Bootstrap values are indicated at each branching point. The solid triangles indicate the noval TTMV species.

Table 1

Sequence similarity of TTMV-222 with other anellovirus strains.

TTMV strain (Genbank accession NO.)	Complete genome, %	Tree Shrew TTMV-222			(KU041847)
		Nucleotides, %			Amino acids, %
		ORF1	ORF2	ORF3	ORF1	ORF2	ORF3
TTMV1-CBD279 (AB026931)	62.6	60.4	60.3	56.2	47.3	43.7	44.4
TTMV2-NLC023 (AB038629)	55.8	57.7	62.7	57.1	41.3	46.5	35.6
TTMV3-NLC026 (AB038630)	56.8	57.5	56.2	61.6	38.3	48.8	40.4
TTMV-LIL-y1 (EF538880)	60.2	59.1	52.9	46.4	34.8	49.4	28.9
TTMV_LY1 (JX134044)	54.3	51.1	49.8	52.8	22.8	28.4	32.2
TTMV_LY2 (JX134045)	54.9	51.5	50.4	45.3	31.9	33.3	29.2
TLMV9-NLC030 (AB038631)	52.6	49.1	49.2	59.6	31.2	32.9	38.9
TLMV6-CBD203 (AB026929)	55.8	51.7	51.6	50.6	32.2	36.6	31.1

Investigation of TTMV-222 prevalence

A case control study was conducted to investigate the association between the prevalence of the newly identified virus and chronic periodontitis. We enrolled 150 chronic periodontitis patients (mean age 45) and 150 periodontal healthy participants (mean age 40) in a case control study. Each diagnosis was based on clinical and radiographic findings according to the inclusion standards described in the Methods. Of the 300 total subjects, no differences in the age or gender distribution were noted between the periodontitis group and the control group (p = 0.647 and 0.729, respectively). Mean PD values, mean CAL, mean PLI, mean GI and the percentage of subjects with bleeding on probing (BOP) were significantly higher in participants with periodontitis than in the healthy group. These data are shown in Tables S2 and S3 in the supplemental material.

Based on the amplified full-length sequence, primers (DL1 5′-T*G*AGTGAAACCACCGAAGTC-3′, DR1 5′-C*G*TTACTTGTTGTCCACCAG-3′, DL2 5′-ACCACGGATTATTCTGCGGC-3′, and DR2 5′-AAAAGACCATGCTCCCCCTC-3′) were designed to screen the samples from the chronic periodontitis patients and periodontally healthy subjects for the novel TTMV. The TTMV-222 virus was detected in 6.00% (9/150) and 1.33% (2/150) of chronic periodontitis patients and periodontally healthy subjects, respectively. The difference in the TTMV-222 prevalence between the two groups was statistically significant (p = 0.032), suggesting that detection of TTMV-222 may be associated with chronic periodontitis, although this finding requires further investigation. Relationships between other factors and detection of TTMV-222 were also detected (Tables 2 and 3). No correlations with age or gender distribution (p = 0.496 and 0.290, respectively) were noted.

Table 2

Relationships between age and TTMV-222.

Group	TTMV-222(+) (mean ± SD)	TTMV-222(−) (mean ± SD)	p
Age	45.00 ± 15.80	42.52 ± 11.69	0.496

Obtained by Student’s t-test (two groups). ap < 0.05.

SD: standard deviation.

Table 3

Relationships between gender, priodontitis and TTMV-222.

Group	n	% with TTMV-222(+)	p
Gender			0.290
Male	7	4.9
Female	4	2.6
Priodontitis			0.032
Positive	9	6.0
Negative	2	1.3

Obtained by Chi-squared test. ap < 0.05.

Discussion

Chronic periodontitis, which can damage the periodontal tissue and result in the early loss of human teeth, is among the most prevalent infectious diseases worldwide32. The causative agents of periodontal disease are generally thought to be periopathogenic bacteria, although the pure bacterial pathogenicity theory cannot cover all of the manifestations of periodontal diseases. Increasing evidence also support a role for pathogens other than bacteria in the pathogenesis of periodontal disease7910. Recently, research into viruses as periodontal pathogenic factors has garnered much attention116333435, especially for herpesviruses. Published evidence has strongly suggested that Epstein-Barr virus (EBV) and human cytomegalovirus (HCMV) are implicated in the pathogenesis of periodontal disease257910. Furthermore, bacteriophages may also play a role in periodontitis, and periodontitis may be beneficial for the expression of some lytic phages16. Nevertheless, the composition of the viral community in the periodontal environment and its association with periodontal diseases is still poorly understood. To investigate the makeup of this community, we applied an unbiased viral metagenomics method to characterize the periodontal virome171920.

Human anelloviruses were first discovered in 1997 and have been classified into three genera in the family Anelloviridae: TTV, genus Alphatorquevirus; TTMV, genus Betatorquevirus; and torque teno midi virus (TTMDV), genus Gammatorquevirus3637. Their genomes consist of small, single-stranded circular DNA, and their capsids have no lipid envelopes. The genomes of anelloviruses range in size from 2.1 to 3.9 kb in length and contain three or four overlapping ORFs and a short, GC-rich un-translated region38. Recent studies have indicated the existence of multiple anellovirus co-infections and interactions between TTVs and the host immune system3940.

In the present study, a new human anellovirus, TTMV-222, was detected by viral metagenomics in gingival tissue from periodontitis patients. Phylogenetic analysis showed that this newly characterized genome belonged to the TTMV genus. Our first phylogenetic tree analysis, which was based on the known TTMV genomic sequences, showed a large cluster of strains that diverged considerably from each other in ORF1 by over 42% at the nucleotide level and over 67% at the amino acid level4142. Three visible clusters were evident in the phylogenetic tree: TTV, TTMDV, and TTMV. TTMV-222 was genetically closest to the TTMV1-CBD279 strain (GenBank accession No. AB026931) and the TTMV2-NLC023 strain (GenBank accession No. AB038629). At the nucleotide level, the full-length sequence of TTMV-222 shared 37.4% and 44.2% sequence heterogeneity while the TTMV-222 ORF1 shared 39.6% and 42.3% sequence heterogeneity with the TTMV1-CBD279 and TTMV2-NLC023 strains, respectively. At the amino acid level, TTMV-222 ORF1 shared 52.7% and 58.7% heterogeneity with the TTMV1-CBD279 and TTMV2-NLC023 strains, respectively (Table 1). The A% nucleotide similarities between TTMV species (TTMV1-9,TTMV-222) for ORF1 are also shown in Table S1. We propose that TTMV-222 may be a new species of torque teno mini virus.

In addition, we investigated the prevalence of TTMV-222 in chronic periodontitis patients and periodontally healthy participants. The TTMV-222 prevalence in the chronic periodontitis group was significantly higher than in the healthy group (p = 0.032), suggesting that TTMV-222 may be more likely to be found in chronic periodontitis patients than in healthy individuals. However, the detection rate was too low to provide a robust connection between chronic periodontitis and TTMV. Several reasons might account for this finding. First and foremost, all samples from chronic periodontitis patients were collected after initial therapy, including supergingival and subgingival scaling. To the best of our knowledge, these treatments would certainly relieve the severity of periodontitis and would reduce the number of viruses4344. Secondly, the reported prevalence rates of TTMV are significantly different among geographical regions. Data from Brazil45, Russia46, Japan47 and Pakistan48 have described prevalence rates varying from 5% to 90%, suggesting that the TTMV prevalence might be related to ethnicity. The epidemic situation of other races requires further investigation. The potential robust association between TTMV and chronic periodontitis requires further investigation.

In summary, this report describes a novel TTMV species from periodontal pockets, which we named TTMV-222. The present study demonstrated the effectiveness of viral metagenomics for the discovery of emerging viruses in the periodontal environment and provided new evidence to implicate viruses in periodontal diseases.

Methods

Sample collection

Patients attending the dental clinics at the Ninth People’s Hospital, School of Medicine, Shanghai Jiao Tong University, China between June 2014 and June 2015 were invited to participate in this study. All participants underwent a comprehensive clinical periodontal examination including X-rays. All participants were systemically healthy, with a minimum of 20 teeth present (excluding third molars). A total of 48 infected gingival epithelium samples of the periodontal pocket from 48 severe periodontitis patients (24 men and 24 women), aged 18–65 years, were collected for high-throughput sequencing. Inclusion criteria for this study were: mean CAL ≥3 mm, mean PD ≥6 mm, BOP, and alveolar bone resorption detected using pantomography49. Exclusion criteria were: diabetes, heart disease, human immunodeficiency virus infection, pregnancy, heavy cigarette smoking (>15 cigarettes/day), previous periodontal treatment (within the previous 6 months), and previous antibiotic intake (within the previous 6 months). Biopsy specimens were periodontal pockets or the gingival epithelium and connective tissue facing the sulcus, which were obtained during the surgical phase aimed at pocket elimination. The samples were rinsed five times with phosphate-buffered saline (PBS) to wash away blood, saliva, and plaque and were then immediately placed in an Eppendorf tube containing 500 μl PBS and transported to the laboratory and stored at −80 °C for analysis. All participants signed a written informed consent form before their inclusion in the study. The Ethics Committee of the Ninth People’s Hospital, School of Medicine, Shanghai Jiao Tong University approved this study. The study was performed in accordance with the principles of the Declaration of Helsinki.

Subjects aged 18–65 years with severe periodontitis (150 total) and 150 periodontal healthy volunteers also aged 18–65 years were recruited for the epidemiological investigation. In the periodontitis group, the inclusion and exclusion standards were the same as those described above, and biopsy specimens were obtained in the same manner. In the control group, participants showed no evident clinical signs of gingival inflammation and had no detectable PD ≤4 mm, no CAL, and no BOP and were considered periodontally healthy. Tissue samples were the gingival epithelium and connective tissue facing the sulcus from periodontally healthy sites. Biopsies were obtained from the sulcular region during a tooth extraction procedure. Some of the subjects in this group exhibited isolated sites with BOP, but for the study purposes, only sites with no BOP were sampled. A full-mouth clinical examination was conducted in each subject. PD (in mm), CAL (in mm), BOP, the GI, and the PLI of each subject were recorded.

Identification of a novel anellovirus by viral metagenomics

The samples were vortexed with small magnetic beads for 5 min and freeze-thawed three times. After centrifugation (10 min at 13,000 × g), 300 μl of supernatant was collected and filtered through a 0.45-μm filter to remove eukaryotic and bacterial cell-sized particles. To digest unprotected nucleic acids (not in viral capsids), a mixture of DNases (Turbo DNase from Ambion, MA, USA, Baseline-ZERO from Epicentre, IL, USA, and benzonase from Novagen, Darmstadt, Germany) and RNase (Thermo Fisher Scientific, MA, USA) was added to the filtrates, which were enriched in viral particles, and the samples were incubated at 37 °C for 90 min50. Viral nucleic acids were then extracted using a QIAamp viral RNA extraction kit (QIAGEN, Dusseldorf, Germany) according to the manufacturer’s instructions, protected from degradation by the addition of an RNase inhibitor (Thermo Fisher Scientific, MA, USA), and stored at −80 °C for future processing. Viral nucleic acid libraries containing both DNA and RNA viral sequences were constructed by random RT-PCR amplification based on the Nextera XT DNA Sample Preparation Kit (Illumina, CA, USA). The libraries were sequenced using the MiSeq platform. The sequencing reads were assigned to 48 data bin-based barcodes. Trimmed sequences from each group were assembled into contigs using the method described by Eric51 in 2011, with a criterion at least 95% identity over 35-bp to merge two fragments. The assembled contigs and singlet sequences were then compared to GenBank using BLASTx. Sequences with E values of ≤10−5 in a BLASTx search were classified as likely originating from a eukaryotic virus, bacterium, phage, eukaryote, other, or unknown based on the taxonomic origin of the sequence with the best E value.

Amplifications of the full-length newly discovered human anellovirus

The full-length anellovirus was amplified by inverted nested PCR based on the sequence obtained from the MiSeq analysis. An overlapping PCR fragment including the remainder of the circular genome (containing the GC-rich region) was obtained and then sequenced. Amplification was performed using Takara LA Taq polymerase and GC buffer I (LA PCR Kit Ver.2.1, TaKaRa, Dalian, China) and was performed as follows: 94 °C for 3 min, followed by five cycles of 1 min at 94 °C, 1 min at 60 °C, and 3.5 min at 72 °C, and then 30 cycles of denaturation at 94 °C for 30 s, annealing at 55 °C for 30 s, and extension at 72 °C for 3.5 min, with an additional 1 s in the extension stage every cycle and a final extension step of 10 min at 72 °C. The second-round cycling conditions were the same as above. The PCR products were excised from 1% agarose gels containing ethidium bromide (0.5 g/ml) and were purified using an AxyPrep DNA Gel Extraction Kit (Axygene, Silicon Valley, USA) according to the manufacturer’s instructions.

Phylogenetic analysis and sequence similarity analysis

The genome sequence of the newly discovered human anellovirus isolated in this study was aligned using ClustalW with the default settings, and the aligned sequences were trimmed to match the human anellovirus sequences obtained from the human anelloviruses detected using ClustalW. Sequence analyses were performed with MegAlign software (DNAStar Inc., Madison, WI, USA). A phylogenetic tree was constructed using the neighbor-joining method with nucleotide p distances and 1,000 bootstrap replicates in the Molecular Evolutionary Genetics Analysis program (MEGA, version 4.0, USA) and by inputting the alignment of the ORF1 genome sequences31. Bootstrap values are indicated at each branching point. A % similarity analysis was also conducted using the MEGA program.

Prevalence investigation

DNA from 300 clinical samples was extracted using a QIAamp DNA Mini kit (QIAGEN, Dusseldorf, Germany) according to the manufacturer’s instructions. Nested PCR assays were performed with a final volume of 50 μl of reaction mixture consisting of 2 μl of extracted DNA from the clinical sample, 25 μl of PrimerSTAR Max Premix, and 10 pmol of primers. The nested PCR parameters were 3 min at 94 °C, followed by 30 s at 94 °C, 30 s at 55 °C, and 50 s at 72 °C for 35 cycles, with a final extension of 5 min at 72 °C. The second-round cycling conditions were the same as described above. Amplification products were resolved by electrophoresis on 1% agarose gels and then sequenced. The detection frequency (%) of the new anellovirus-positive subjects was calculated.

Statistical analyses

All statistical analyses were conducted using SPSS software (ver. 20.0; IBM Corporation, Armonk, NY, USA) with the significance level set to p < 0.05. Descriptive statistics (means, standard deviations, and percentages) were calculated for subjects’ socio-demographic characteristics (gender and age). Chi-squared and Fisher exact tests were performed to examine differences in the prevalence of the newly discovered virus and in one of the clinical indexes (BOP) between the case group and the control group. Student’s t-tests were used to examine differences in other clinical indexes (PD, GI, PLI, and CAL).

Additional Information

How to cite this article: Zhang, Y. et al. A novel species of torque teno mini virus (TTMV) in gingival tissue from chronic periodontitis patients. Sci. Rep. 6, 26739; doi: 10.1038/srep26739 (2016).