The mutation rate of mycobacterial repetitive unit loci in strains of M. tuberculosisfrom cynomolgus macaque infection
- Mark N Ragheb†1,
- Christopher B Ford†1,
- Michael R Chase1,
- Philana Ling Lin2,
- JoAnne L Flynn3 and
- Sarah M Fortune1, 4, 5Email author
© Ragheb et al.; licensee BioMed Central Ltd. 2013
Received: 20 December 2012
Accepted: 26 February 2013
Published: 5 March 2013
Mycobacterial interspersed repetitive units (MIRUs) are minisatellites within the Mycobacterium tuberculosis (Mtb) genome. Copy number variation (CNV) in MIRU loci is used for epidemiological typing, making the rate of variation important for tracking the transmission of Mtb strains. In this study, we developed and assessed a whole-genome sequencing (WGS) approach to detect MIRU CNV in Mtb. We applied this methodology to a panel of Mtb strains isolated from the macaque model of tuberculosis (TB), the animal model that best mimics human disease. From these data, we have estimated the rate of MIRU variation in the host environment, providing a benchmark rate for future epidemiologic work.
We assessed variation at the 24 MIRU loci used for typing in a set of Mtb strains isolated from infected cynomolgus macaques. We previously performed WGS of these strains and here have applied both read depth (RD) and paired-end mapping (PEM) metrics to identify putative copy number variants. To assess the relative power of these approaches, all MIRU loci were resequenced using Sanger sequencing. We detected two insertion/deletion events both of which could be identified as candidates by PEM criteria. With these data, we estimate a MIRU mutation rate of 2.70 × 10-03 (95% CI: 3.30 × 10-04- 9.80 × 10-03) per locus, per year.
Our results represent the first experimental estimate of the MIRU mutation rate in Mtb. This rate is comparable to the highest previous estimates gathered from epidemiologic data and meta-analyses. Our findings allow for a more rigorous interpretation of data gathered from MIRU typing.
KeywordsMycobacterium tuberculosis Mycobacterial interspersed repetitive units MIRU Molecular epidemiology Copy number variation Whole-genome sequencing Read depth Paired-end mapping Mutation rate
The ability to genetically differentiate among microbial strains facilitates tracing the origins and spread of bacterial pathogens, including Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis (TB). Various methods for genetically typing clinical strains of Mtb have been developed [1–4]. This includes typing strains through copy number variation (CNV) in mycobacterial interspersed repetitive units (MIRUs), which are minisatellite loci in the Mtb genome . This approach, termed MIRU-VNTR (variable number of tandem repeats; hereafter, MIRU) typing, distinguishes genetically divergent strains rapidly and with relatively high accuracy . As a result, MIRU analysis has been employed in a wide array of epidemiological studies [7–10], where an identical MIRU profile between isolates is interpreted as a recent transmission event . This is based on the assumption that over short periods of time, a change in MIRU copy number is unlikely. Therefore, accurate estimation of the MIRU mutation rate is essential to infer the relationship between transmitted strains based on typing profile. The limits of this approach are highlighted by a recent study of a TB outbreak in British Columbia, where MIRU typing identified a single clonal outbreak, yet higher resolution WGS established two separate, simultaneous outbreaks .
Approximating the rate of MIRU CNV based on the profiles of clinical strains has proven to be a challenge. The approaches taken by various groups to estimate a rate have relied primarily on modeling of epidemiological data and meta-analyses [13–16]. The discrepancies in methodology, loci analyzed, and underlying assumptions between different studies have resulted in estimates ranging over two orders of magnitude. For example, Grant et al. (2008) analyzed copy number changes between Mtb lineages and used previous research estimating the time of the most recent common ancestor between lineages to estimate a MIRU mutation rate, yielding a per locus, per year rate of 1.05 × 10-5. Reyes and Tanaka (2010) used an infinite alleles model to define a relative rate and then benchmarked this rate against estimates of the IS6110 mutation rate to infer a MIRU mutation rate (7.00 × 10-4- 1.50 × 10-2). More recently, Aandahl et al. (2012) developed a stepwise mutation model for MIRU evolution and utilized Bayesian statistics to estimate the MIRU mutation rate of previously gathered epidemiological data (3.55 × 10-3), supporting the estimates provided by Reyes and Tanaka (2010).
Here, we seek to experimentally determine a mutation rate by assessing CNV in Mtb strains isolated from cynomolgus macaques , an animal model of Mtb infection that closely recapitulates the course of human disease . The genomes from Mtb isolated from infected macaques were previously sequenced using the Illumina platform . While protocols on detecting single nucleotide polymorphisms (SNPs) and small insertion/deletions (indels) using Illumina sequencing data are well established, it is less clear how to best determine the copy number in minisatellite loci. As WGS becomes increasingly common in epidemiological studies, the ability to establish MIRU copy number from sequencing data becomes important for the analysis of new WGS data in the context of previously existing typing data.
The length of Illumina reads fails to span the majority of MIRU repeats in a locus, which range from one to five or more repeats of 40–100 basepairs each in Mtb. Thus, short read sequencing cannot capture unique sequence and define copy number. Recently, researchers have utilized read depth (RD), a measure of the density of sequencing reads at each nucleotide in the genome, in order to identify CNV [19, 20]. This approach relies on the observation that the absolute number of sequencing reads mapped to a reference genome is proportional to the copy number of a particular strain . However, this approach has been primarily utilized to identify large (> 1 kb) variants in human tandem repeats, and RD has not been successfully applied to assess smaller minisatellite variation in microbes. Similarly, attempts to identify structural variation, including CNV, using mate pair distance have been previously employed in human genome studies [21–23]. This approach, termed paired-end mapping (PEM), utilizes the likely distance between paired reads from Illumina sequencing to identify structural variants. When mapping back to a reference genome, if the distance between paired reads is discordant from the expected value, it is suggestive of CNV relative to the reference genome. PEM has been successfully employed to identify large insertions or deletions. However, bacterial minisatellite CNV produces only small discordances relative to the reference genome which fall within the expected distribution of mate pair distance, making such events difficult to detect.
In this study, we sought to define a WGS methodology useful for identifying MIRU CNV by mapping sequencing reads to a single copy-number reference genome. We validated the WGS analyses by resequencing all MIRU loci via Sanger sequencing, which allowed us to assess the accuracy of using RD and PEM approaches to identify minisatellite variations. We then used our WGS and Sanger resequencing data to estimate a MIRU mutation rate during the course of infection. This rate will help guide the analysis of epidemiological data and provide a preliminary understanding of site-specific mutability in Mtb.
Identifying MIRU CNV’s using WGS
We analyzed WGS data from 16 sequenced strains  for CNV at 22 of the 24 MIRU loci currently standardized for strain typing . The remaining two loci showed poor read density and were excluded from this portion of the analysis. We utilized RD and PEM, both of which have been shown to correlate with copy number, to identify MIRU CNV. We hypothesized that any strain containing a MIRU insertion would exhibit an increase in RD and a decrease in mate pair distance relative to the input strain, while a deletion would exhibit decreased RD and increased mate pair distance [19, 20]. However, either method is complicated by the inherent limitations of mapping short reads.
List of CNVs identified via WGS and Sanger resequencing
Read depth candidate strains
Mate pair candidate strains
3; (H-3 = 4)
2; (G-3 = 1)
Comparing mean mate pair distance between H37Rv and single copy reference genomes
Single copy (MIRU)
H37Rv (+/− 100 bp)
Single copy (+/− 100 bp)
In order to assess the validity of the RD and PEM approaches we used Sanger sequencing to quantify the number of MIRU repeats at each of the 24 standard MIRU loci described previously . To improve our estimate of the MIRU mutation rate, we assessed copy number in 17 additional strains isolated from cynomolgus macaques that were experimentally infected with the Mtb strain Erdman. All 33 strains were assessed at the 24-MIRU loci previously described. From the 792 loci sequenced, two of the four putative indels identified by PEM were confirmed, and no new variants were discovered (Table 1). Locus 802 in strain H-3 increased from three to four repeats and was identified by both RD and PEM, and locus 3192 in strain G-3 decreased from two to one copy and was identified solely by PEM.
Estimation of the MIRU mutation rate during the course of disease
Values of published MIRU mutation rates
Mean MIRU mutation rate
Estimated rate, macaque model
2.70 × 10-3
3.30 × 10-4
9.80 × 10-3
Aandahl et al. 2012
3.55 × 10-3
8.51 × 10-4
1.15 × 10-2
Reyes and Tanaka 2010
7.00 × 10-4
1.50 × 10-2
Wirth et al. 2008
1.23 × 10-4
1.96 × 10-5
7.93 × 10-4
Grant et al. 2008*
1.05 × 10-5
Here we have utilized both conventional and next-generation sequencing approaches to define the number of copy variants at MIRU loci that arose during the course of macaque infection. Approaches to detect MIRU CNV by Illumina sequencing are increasingly important given the expanding use of WGS in molecular epidemiology. However, mapping ambiguities due to the repetitive nature of MIRU loci complicate copy number assessment by WGS. In this work, we sought to reduce mapping ambiguity and identify CNV in Mtb by mapping to a single copy-number genome. Sequencing technologies have advanced since the sequencing of these strains, and longer reads with less coverage variability may improve specificity in future work [31, 32]. Additionally, it is thought that variability in both RD and PEM arises from local disparities in GC content [25, 33]. Future work may exploit advances in correcting for GC bias [19, 32, 34–36], in combination with the approaches described here, to more reliably detect MIRU CNV from WGS data.
With the MIRU CNV identified in strains isolated from cynomolgus macaques, we have estimated the mutation rate at MIRU loci. Our in vivo MIRU mutation rate, 2.70 × 10-03 per locus, per year, is most consistent with the highest published estimates. Variability in these estimates may be partially driven by differences in MIRU loci analyzed, the epidemiology of the strains used, and differences in the assumptions of the models used to estimate a rate. The resulting differences in rate estimates have motivated extensive debate in the literature [15, 37, 38]. Our estimate is derived from an alternative, experimental approach, relying on the cynomolgus macaque model of TB infection to assess the mutation rate of the 24 loci standard used in MIRU typing. Though our analysis is somewhat limited by a relatively small signal, strikingly, our estimate closely aligns with previous rates.
What are the biologic consequences of a high mutation rate at MIRU loci? Interestingly, most MIRUs are located in intergenic regions and are hypothesized to be transcribed as part of a polycistronic operon. Several MIRU elements are located within the coding region of well-described two-component regulatory systems as well as genes essential for virulence and host adaptation . It is interesting to hypothesize that rapid genetic variation at these MIRU loci may have effects on the transcription of the regulon, thus generating population diversity. Some evidence exists that CNV variation at MIRU loci may result in transcriptional changes of the downstream gene within a MIRU locus [39, 40], though further characterization is required to establish the biologic relevance of these loci and the role of variation in MIRU elements.
Repeat variation is a well-established means of generating locus specific mutation in other microbial genomes [41, 42]. The rate established here is comparable to tandem repeat variation rates in other organisms. For example, in Bacillus anthracis the estimated mutation rate of tandem repeat loci is reported to be roughly 10-05 to 10-04 per generation , while in pathogenic E. coli O157:H7, the rate is reported to be roughly 6.4 × 10-04 per generation . The MIRU per generation mutation rate is between 5.6 × 10-06 and 7.64 × 10-05, using the lower (18 hours) and upper (240 hours) estimates of generation times previously described . It is important to note the loci analyzed in this study were selected from the larger set of MIRU loci for their relative stability, making them ideal for typing. The rate of variation in the remaining loci may vary from the rate reported here, especially in loci previously identified as hypermutable . While further work is needed to determine the biological consequences of MIRU variation, we have shown that there is potential for detecting MIRU variants by WGS and that the rate uncovered from macaque infection is consistent with the highest previous estimates.
CNV in typing markers is an essential tool to differentiate and classify clinical strains, and quantitation of marker variation allows for enhanced interpretation of epidemiological data. In this study, we have used the macaque model of Mtb infection to estimate the MIRU mutation rate during the course of disease, and we have explored the use of WGS to assess MIRU copy number. Subsequent Sanger resequencing confirmed two of the four MIRU indels identified by PEM to a reduced copy genome, and from this we have estimated a per locus, per year mutation rate of 2.70 × 10-03. This value agrees with the higher published estimations of MIRU mutation rates. Further assessment of RD and PEM as indicators for copy number may streamline minisatellite detection via WGS.
Preparation of isolates and Illumina sequencing
Infection of macaques and isolation of strains was performed as previously described . Briefly, cynomolgus macaques were infected with a low dose (roughly 25 CFU/macaque), virulent strain (Mtb Erdman) by bronchoscopy. Infected macaques were allowed to progress to either latent or active disease. Nine infected macaques were selected (four active, three latent, two reactivated) and 33 bacterial isolates from 17 different lesions were chosen for study. Colonies were expanded for extraction of genomic DNA as described previously. Minimal expansion occurred between strain isolation and genomic extraction. An Illumina Genome Analyzer (Illumina) was used for WGS of isolated strains. A detailed protocol of WGS data analysis has been previously published .
WGS analysis of MIRU regions
75 basepair paired-end read data gathered from sequencing at the Broad Institute of MIT and Harvard were analyzed for MIRU CNV. Only strains containing sufficient sequencing coverage at the MIRU locus were analyzed. Sequences were mapped to a single copy H37Rv reference genome. This genome was created by reducing MIRU elements to a single copy at the 24 loci analyzed in the H37Rv genome [GenBank: AL123456]. If two MIRUs at the same locus contained a greater than 4 SNP difference, they were treated as unique MIRUs (denoted A and B) and were not collapsed to a single unit. Illumina fastq files were mapped to this reference genome with SSAHA2 using the Solexa defaults and allowing for paired end reads up to 700 bp . Proper pairs were extracted using samtools. RD was calculated by including mate pairs traversing each reference coordinate with a perl script . The two-sided Wilcoxon rank sum test was used to compare the difference in mate pair distances when mapping to the H37Rv genome versus the single copy-number genome (Mathworks, Natick MA). This analysis was done both on reads mapping within the MIRU and reads mapping in the surrounding +/−100 bp. For this comparison Illumina fastq files were mapped to the H37Rv genome using the same parameters as used when mapping to the single copy-number genome.
To identify CNV in MIRUs, RD and mate pair distance values were obtained for 16 strains at 24 MIRU loci and 100 basepairs upstream and downstream of the MIRU element. Values at each coordinate were normalized to its mean RD and mate pair distance values for the window assessed. MATLAB was used to generate plots, mean RD values and standard deviations (MathWorks, Natick MA).
Sanger sequencing of isolates and estimation of MIRU mutation rate
Where n indels = total number of MIRU insertion or deletion events, t (a-i) = duration of infection per macaque (in days), g (a-i) = number of strains analyzed per macaque, and l = total number of MIRU loci analyzed. n indels is divided by the sum of the product of t, g (a-i) , and l(a-i) per macaque. Because infection length and isolates acquired is variable between macaques, the formula must sum the product acquired for each individual macaque. A Poisson distribution was used to model the number of indels and estimate the 95% confidence interval. Estimation of the mutation rate and the confidence interval was generated using MATLAB (MathWorks, Natick MA).
Mycobacterial interspersed repetitive unit
Copy number variation
Variable number of tandem repeats
Single nucleotide polymorphism
This work was supported by a New Innovator’s Award, DP2 0D001378 from the Director’s Office of the National Institute of Health to S.M.F., by a subcontract from National Institute for Allergy and Infectious Diseases (NIAID) U19 AI076217 to S.M.F., by the US National Institutes of Health (NIH) RO1 HL075845 to J.L.F. and by the Bill and Melinda Gates Foundation (J.L.F.). The genome sequencing has been funded in part with federal funds from the National Institute of Allergy and Infectious Disease, US NIH, US Department of Health and Human Services, under contract no. HHSN266200400001C. Development of methodology to analyze genome sequences was funded in part by a Clinical Scientist Development Award to S.M.F. from the Doris Duke Charitable Foundation. We thank D. Gurgil of the Enterprise Research IS group at Partners Healthcare for his support and for provision of the HPC facilities. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
- Eisenach KD, Crawford JT, Bates JH: Genetic relatedness among strains of the mycobacterium tuberculosis complex. Analysis of restriction fragment heterogeneity using cloned DNA probes. Am Rev Respir Dis. 1986, 133 (6): 1065-1068.PubMedGoogle Scholar
- Eisenach KD, Crawford JT, Bates JH: Repetitive DNA sequences as probes for mycobacterium tuberculosis. J Clin Microbiol. 1988, 26 (11): 2240-2245.PubMed CentralPubMedGoogle Scholar
- Thierry D, Matsiota-Bernard P, Pitsouni E, Costopoulos C, Guesdon JL: Use of the insertion element IS6110 for DNA fingerprinting of mycobacterium tuberculosis isolates presenting various profiles of drug susceptibility. FEMS Immunol Med Microbiol. 1993, 6 (4): 287-297. 10.1111/j.1574-695X.1993.tb00341.x.View ArticlePubMedGoogle Scholar
- Mazars E, Lesjean S, Banuls AL, Gilbert M, Vincent V, Gicquel B, Tibayrenc M, Locht C, Supply P: High-resolution minisatellite-based typing as a portable approach to global analysis of Mycobacterium tuberculosis molecular epidemiology. Proc Natl Acad Sci U S A. 2001, 98 (4): 1901-1906. 10.1073/pnas.98.4.1901.PubMed CentralView ArticlePubMedGoogle Scholar
- Supply P, Magdalena J, Himpens S, Locht C: Identification of novel intergenic repetitive units in a mycobacterial two-component system operon. Mol Microbiol. 1997, 26 (5): 991-1003. 10.1046/j.1365-2958.1997.6361999.x.View ArticlePubMedGoogle Scholar
- Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rusch-Gerdes S, Willery E, Savine E, de Haas P, van Deutekom H, Roring S: Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of mycobacterium tuberculosis. J Clin Microbiol. 2006, 44 (12): 4498-4510. 10.1128/JCM.01392-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Spurgiesz RS, Quitugua TN, Smith KL, Schupp J, Palmer EG, Cox RA, Keim P: Molecular typing of mycobacterium tuberculosis by using nine novel variable-number tandem repeats across the Beijing family and low-copy-number IS6110 isolates. J Clin Microbiol. 2003, 41 (9): 4224-4230. 10.1128/JCM.41.9.4224-4230.2003.View ArticlePubMedGoogle Scholar
- Sun YJ, Bellamy R, Lee AS, Ng ST, Ravindran S, Wong SY, Locht C, Supply P, Paton NI: Use of mycobacterial interspersed repetitive unit-variable-number tandem repeat typing to examine genetic diversity of mycobacterium tuberculosis in Singapore. J Clin Microbiol. 2004, 42 (5): 1986-1993. 10.1128/JCM.42.5.1986-1993.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Valcheva V, Mokrousov I, Rastogi N, Narvskaya O, Markova N: Molecular characterization of mycobacterium tuberculosis isolates from different regions of Bulgaria. J Clin Microbiol. 2008, 46 (3): 1014-1018. 10.1128/JCM.01841-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Christianson S, Wolfe J, Orr P, Karlowsky J, Levett PN, Horsman GB, Thibert L, Tang P, Sharma MK: Evaluation of 24 locus MIRU-VNTR genotyping of mycobacterium tuberculosis isolates in Canada. Tuberculosis (Edinb). 2010, 90 (1): 31-38. 10.1016/j.tube.2009.12.003.View ArticleGoogle Scholar
- Yeh RW, Ponce De Leon A, Agasino CB, Hahn JA, Daley CL, Hopewell PC, Small PM: Stability of mycobacterium tuberculosis DNA genotypes. J Infect Dis. 1998, 177 (4): 1107-1111. 10.1086/517406.View ArticlePubMedGoogle Scholar
- Gardy JL, Johnston JC, Ho Sui SJ, Cook VJ, Shah L, Brodkin E, Rempel S, Moore R, Zhao Y, Holt R: Whole-genome sequencing and social-network analysis of a tuberculosis outbreak. N Engl J Med. 2011, 364 (8): 730-739. 10.1056/NEJMoa1003176.View ArticlePubMedGoogle Scholar
- Grant A, Arnold C, Thorne N, Gharbia S, Underwood A: Mathematical modelling of mycobacterium tuberculosis VNTR loci estimates a very slow mutation rate for the repeats. J Mol Evol. 2008, 66 (6): 565-574. 10.1007/s00239-008-9104-6.View ArticlePubMedGoogle Scholar
- Wirth T, Hildebrand F, Allix-Beguec C, Wolbeling F, Kubica T, Kremer K, van Soolingen D, Rusch-Gerdes S, Locht C, Brisse S: Origin, spread and demography of the mycobacterium tuberculosis complex. PLoS Pathog. 2008, 4 (9): e1000160-10.1371/journal.ppat.1000160.PubMed CentralView ArticlePubMedGoogle Scholar
- Reyes JF, Tanaka MM: Mutation rates of spoligotypes and variable numbers of tandem repeat loci in mycobacterium tuberculosis. Infect Genet Evol. 2010, 10 (7): 1046-1051. 10.1016/j.meegid.2010.06.016.View ArticlePubMedGoogle Scholar
- Aandahl RZ, Reyes JF, Sisson SA, Tanaka MM: A model-based Bayesian estimation of the rate of evolution of VNTR loci in mycobacterium tuberculosis. PLoS Comput Biol. 2012, 8 (6): e1002573-10.1371/journal.pcbi.1002573.PubMed CentralView ArticlePubMedGoogle Scholar
- Ford CB, Lin PL, Chase MR, Shah RR, Iartchouk O, Galagan J, Mohaideen N, Ioerger TR, Sacchettini JC, Lipsitch M: Use of whole genome sequencing to estimate the mutation rate of mycobacterium tuberculosis during latent infection. Nat Genet. 2011, 43 (5): 482-486. 10.1038/ng.811.PubMed CentralView ArticlePubMedGoogle Scholar
- Capuano SV, 3rd Croix DA, Pawar S, Zinovik A, Myers A, Lin PL, Bissel S, Fuhrman C, Klein E, Flynn JL: Experimental mycobacterium tuberculosis infection of cynomolgus macaques closely resembles the various manifestations of human M. Tuberculosis infection. Infect Immun. 2003, 71 (10): 5831-5844. 10.1128/IAI.71.10.5831-5844.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Yoon S, Xuan Z, Makarov V, Ye K, Sebat J: Sensitive and accurate detection of copy number variants using read depth of coverage. Genome Res. 2009, 19 (9): 1586-1592. 10.1101/gr.092981.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang ZD, Gerstein MB: Detection of copy number variation from array intensity and sequencing read depth using a stepwise Bayesian model. BMC Bioinforma. 2010, 11: 539-10.1186/1471-2105-11-539.View ArticleGoogle Scholar
- Korbel JO, Urban AE, Affourtit JP, Godwin B, Grubert F, Simons JF, Kim PM, Palejev D, Carriero NJ, Du L: Paired-end mapping reveals extensive structural variation in the human genome. Science. 2007, 318 (5849): 420-426. 10.1126/science.1149504.PubMed CentralView ArticlePubMedGoogle Scholar
- Medvedev P, Fiume M, Dzamba M, Smith T, Brudno M: Detecting copy number variation with mated short reads. Genome Res. 2010, 20 (11): 1613-1622. 10.1101/gr.106344.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Tuzun E, Sharp AJ, Bailey JA, Kaul R, Morrison VA, Pertz LM, Haugen E, Hayden H, Albertson D, Pinkel D: Fine-scale structural variation of the human genome. Nat Genet. 2005, 37 (7): 727-732. 10.1038/ng1562.View ArticlePubMedGoogle Scholar
- Korbel JO, Abyzov A, Mu XJ, Carriero N, Cayting P, Zhang Z, Snyder M, Gerstein MB: PEMer: a computational framework with simulation-based error models for inferring genomic structural variants from massive paired-end sequencing data. Genome Biol. 2009, 10 (2): R23-10.1186/gb-2009-10-2-r23.PubMed CentralView ArticlePubMedGoogle Scholar
- Hormozdiari F, Alkan C, Eichler EE, Sahinalp SC: Combinatorial algorithms for structural variation detection in high-throughput sequenced genomes. Genome Res. 2009, 19 (7): 1270-1278. 10.1101/gr.088633.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Alkan C, Kidd JM, Marques-Bonet T, Aksay G, Antonacci F, Hormozdiari F, Kitzman JO, Baker C, Malig M, Mutlu O: Personalized copy number and segmental duplication maps using next-generation sequencing. Nat Genet. 2009, 41 (10): 1061-1067. 10.1038/ng.437.PubMed CentralView ArticlePubMedGoogle Scholar
- He D, Hormozdiari F, Furlotte N, Eskin E: Efficient algorithms for tandem copy number variation reconstruction in repeat-rich regions. Bioinformatics. 2011, 27 (11): 1513-1520. 10.1093/bioinformatics/btr169.PubMed CentralView ArticlePubMedGoogle Scholar
- Abyzov A, Urban AE, Snyder M, Gerstein M: CNVnator: an approach to discover, genotype, and characterize typical and atypical CNVs from family and population genome sequencing. Genome Res. 2011, 21 (6): 974-984. 10.1101/gr.114876.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen K, Wallis JW, McLellan MD, Larson DE, Kalicki JM, Pohl CS, McGrath SD, Wendl MC, Zhang Q, Locke DP: BreakDancer: an algorithm for high-resolution mapping of genomic structural variation. Nat Methods. 2009, 6 (9): 677-681. 10.1038/nmeth.1363.PubMed CentralView ArticlePubMedGoogle Scholar
- Treangen TJ, Salzberg SL: Repetitive DNA and next-generation sequencing: computational challenges and solutions. Nat Rev Genet. 2012, 13 (1): 36-46.Google Scholar
- Eid J, Fehr A, Gray J, Luong K, Lyle J, Otto G, Peluso P, Rank D, Baybayan P, Bettman B: Real-time DNA sequencing from single polymerase molecules. Science. 2009, 323 (5910): 133-138. 10.1126/science.1162986.View ArticlePubMedGoogle Scholar
- Loman NJ, Constantinidou C, Chan JZ, Halachev M, Sergeant M, Penn CW, Robinson ER, Pallen MJ: High-throughput bacterial genome sequencing: an embarrassment of choice, a world of opportunity. Nat Rev Microbiol. 2012, 10 (9): 599-606. 10.1038/nrmicro2850.View ArticlePubMedGoogle Scholar
- Dohm JC, Lottaz C, Borodina T, Himmelbauer H: Substantial biases in ultra-short read data sets from high-throughput DNA sequencing. Nucleic Acids Res. 2008, 36 (16): e105-10.1093/nar/gkn425.PubMed CentralView ArticlePubMedGoogle Scholar
- Boeva V, Zinovyev A, Bleakley K, Vert JP, Janoueix-Lerosey I, Delattre O, Barillot E: Control-free calling of copy number alterations in deep-sequencing data using GC-content normalization. Bioinformatics. 2011, 27 (2): 268-269. 10.1093/bioinformatics/btq635.PubMed CentralView ArticlePubMedGoogle Scholar
- Miller CA, Hampton O, Coarfa C, Milosavljevic A: ReadDepth: a parallel R package for detecting copy number alterations from short sequencing reads. PLoS One. 2011, 6 (1): e16327-10.1371/journal.pone.0016327.PubMed CentralView ArticlePubMedGoogle Scholar
- Benjamini Y, Speed TP: Summarizing and correcting the GC content bias in high- throughput sequencing. Nucleic Acids Res. 2012, 40 (10): e72-10.1093/nar/gks001.PubMed CentralView ArticlePubMedGoogle Scholar
- Supply P, Niemann S, Wirth T: On the mutation rates of spoligotypes and variable numbers of tandem repeat loci of Mycobacterium tuberculosis. Infect Genet Evol. 2011, 11 (2): 251-252. 10.1016/j.meegid.2010.12.009.View ArticlePubMedGoogle Scholar
- Tanaka MM, Reyes JF: Mutation rate of VNTR loci in mycobacterium tuberculosis: response to supply et al. Infect Genet Evol. 2011, 11 (6): 1189-1190. 10.1016/j.meegid.2011.01.014. discussion 1191–1182View ArticlePubMedGoogle Scholar
- Olsen I, Balasingham SV, Davidsen T, Debebe E, Rodland EA, van Soolingen D, Kremer K, Alseth I, Tonjum T: Characterization of the major formamidopyrimidine-DNA glycosylase homolog in mycobacterium tuberculosis and its linkage to variable tandem repeats. FEMS Immunol Med Microbiol. 2009, 56 (2): 151-161. 10.1111/j.1574-695X.2009.00562.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Tantivitayakul P, Panapruksachat S, Billamas P, Palittapongarnpim P: Variable number of tandem repeat sequences act as regulatory elements in mycobacterium tuberculosis. Tuberculosis (Edinb). 2010, 90 (5): 311-318. 10.1016/j.tube.2010.08.003.View ArticleGoogle Scholar
- Linton D, Gilbert M, Hitchen PG, Dell A, Morris HR, Wakarchuk WW, Gregson NA, Wren BW: Phase variation of a beta-1,3 galactosyltransferase involved in generation of the ganglioside GM1-like lipo-oligosaccharide of campylobacter jejuni. Mol Microbiol. 2000, 37 (3): 501-514.View ArticlePubMedGoogle Scholar
- Erwin AL, Bonthuis PJ, Geelhood JL, Nelson KL, McCrea KW, Gilsdorf JR, Smith AL: Heterogeneity in tandem octanucleotides within haemophilus influenzae lipopolysaccharide biosynthetic gene losA affects serum resistance. Infect Immun. 2006, 74 (6): 3408-3414. 10.1128/IAI.01540-05.PubMed CentralView ArticlePubMedGoogle Scholar
- Keim P, Van Ert MN, Pearson T, Vogler AJ, Huynh LY, Wagner DM: Anthrax molecular epidemiology and forensics: using the appropriate marker for different evolutionary scales. Infect Genet Evol. 2004, 4 (3): 205-213. 10.1016/j.meegid.2004.02.005.View ArticlePubMedGoogle Scholar
- Vogler AJ, Keys C, Nemoto Y, Colman RE, Jay Z, Keim P: Effect of repeat copy number on variable-number tandem repeat mutations in Escherichia coli O157:H7. J Bacteriol. 2006, 188 (12): 4253-4263. 10.1128/JB.00001-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Lin PL, Rodgers M, Smith L, Bigbee M, Myers A, Bigbee C, Chiosea I, Capuano SV, Fuhrman C, Klein E: Quantitative comparison of active and latent tuberculosis in the cynomolgus macaque model. Infect Immun. 2009, 77 (10): 4631-4642. 10.1128/IAI.00592-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Ning Z, Cox AJ, Mullikin JC: SSAHA: a fast search method for large DNA databases. Genome Res. 2001, 11 (10): 1725-1729. 10.1101/gr.194201.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The sequence alignment/Map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.