Search

Background DNA sequence composition affects meiotic recombination. However, the correlation between tandem repeat composition and meiotic recombination in common wheat (Triticum aestivum L.) is unclear. Results Non-denaturing fluorescent in situ hybridization (ND-FISH) with oligonucleotide (oligo) probes derived from tandem repeats and single-copy FISH were used to investigate recombination in three kinds of the long arm of wheat 5A chromosome (5AL). 5AL535–18/275 arm carries the tandem repeats pTa-535, Oligo-18, and pTa-275, 5AL119.2–18/275 arm carries the tandem repeats pSc119.2, Oligo-18 and pTa-275, and 5AL119.2 arm carries the tandem repeats pSc119.2. In the progeny of 5AL535–18/275 × 5AL119.2, double recombination occurred between pSc119.2 and pTa-535 clusters (119–535 interval), and between pTa-535 and Oligo-18/pTa-275 clusters (535–18 interval). The recombination rate in the 119–535 interval in the progeny of 5AL535–18/275 × 5AL119.2–18/275 was higher than that in the progeny of 5AL535–18/275 × 5AL119.2. Recombination in the 119–535 interval produced 5AL119 + 535 segments with pTa-535 and pSc119.2 tandem repeats and 5ALNo segments without these repeats. The 5AL119 + 535 and 5ALNo segments were localized between the signal sites of the single-copy probes SC5A-479 and SC5A-527. The segment between SC5A-479 and SC5A-527 in the metaphase 5ALNo was significantly longer than that in the metaphase 5AL119 + 535. Conclusion The structural variations caused by tandem repeats might be one of the factors affecting meiotic recombination in wheat. Meiotic recombination aggregated two kinds of tandemly repeated clusters into the same chromosome, making the metaphase chromosome more condensed. To conclude, our study provides a robust tool to measure meiotic recombination and select parents for wheat breeding programs.

Background Telomeric DNA is typically comprised of G-rich tandem repeat motifs and maintained by telomerase (Greider CW, Blackburn EH; Cell 51:887–898; 1987). In eukaryotes lacking telomerase, a variety of DNA repair and DNA recombination based pathways for telomere maintenance have evolved in organisms normally dependent upon telomerase for telomere elongation (Webb CJ, Wu Y, Zakian VA; Cold Spring Harb Perspect Biol 5:a012666; 2013); collectively called Alternative Lengthening of Telomeres (ALT) pathways. By measuring (TTAGGG) n tract lengths from the same large DNA molecules that were optically mapped, we simultaneously analyzed telomere length dynamics and subtelomere-linked structural changes at a large number of specific subtelomeric loci in the ALT-positive cell lines U2OS, SK-MEL-2 and Saos-2. Results Our results revealed loci-specific ALT telomere features. For example, while each subtelomere included examples of single molecules with terminal (TTAGGG) n tracts as well as examples of recombinant telomeric single molecules, the ratio of these molecules was subtelomere-specific, ranging from 33:1 (19p) to 1:25 (19q) in U2OS. The Saos-2 cell line shows a similar percentage of recombinant telomeres. The frequency of recombinant subtelomeres of SK-MEL-2 (11%) is about half that of U2OS and Saos-2 (24 and 19% respectively). Terminal (TTAGGG) n tract lengths and heterogeneity levels, the frequencies of telomere signal-free ends, and the frequency and size of retained internal telomere-like sequences (ITSs) at recombinant telomere fusion junctions all varied according to the specific subtelomere involved in a particular cell line. Very large linear extrachromosomal telomere repeat (ECTR) DNA molecules were found in all three cell lines; these are in principle capable of templating synthesis of new long telomere tracts via break-induced repair (BIR) long-tract DNA synthesis mechanisms and contributing to the very long telomere tract length and heterogeneity characteristic of ALT cells. Many of longest telomere tracts (both end-telomeres and linear ECTRs) displayed punctate CRISPR/Cas9-dependent (TTAGGG) n labeling patterns indicative of interspersion of stretches of non-canonical telomere repeats. Conclusion Identifying individual subtelomeres and characterizing linked telomere (TTAGGG) n tract lengths and structural changes using our new single-molecule methodologies reveals the structural consequences of telomere damage, repair and recombination mechanisms in human ALT cells in unprecedented molecular detail and significant differences in different ALT-positive cell lines.

Background Chloroplast (cp) genome information would facilitate the development and utilization of Taxodium resources. However, cp genome characteristics of Taxodium were poorly understood. Results We determined the complete cp genome sequences of T. distichum, T. mucronatum, and T. ascendens. The cp genomes are 131,947 bp to 132,613 bp in length, encode 120 genes with the same order, and lack typical inverted repeat (IR) regions. The longest small IR, a 282 bp trnQ-containing IR, were involved in the formation of isomers. Comparative analysis of the 3 cp genomes showed that 91.57% of the indels resulted in the periodic variation of tandem repeat (TR) motifs and 72.46% single nucleotide polymorphisms (SNPs) located closely to TRs, suggesting a relationship between TRs and mutational dynamics. Eleven hypervariable regions were identified as candidates for DNA barcode development. Hypothetical cp open reading frame 1(Ycf1) was the only one gene that has an indel in coding DNA sequence, and the indel is composed of a long TR. When extended to cupressophytes, ycf1 genes have undergone a universal insertion of TRs accompanied by extreme length expansion. Meanwhile, ycf1 also located in rearrangement endpoints of cupressophyte cp genomes. All these characteristics highlight the important role of repeats in the evolution of cp genomes. Conclusions This study added new evidence for the role of repeats in the dynamics mechanism of cp genome mutation and rearrangement. Moreover, the information of TRs and hypervariable regions would provide reliable molecular resources for future research focusing on the infrageneric taxa identification, phylogenetic resolution, population structure and biodiversity for the genus Taxodium and Cupressophytes.

Background Simple sequence repeats (SSRs) are tandem repeats of DNA that have been used to develop robust genetic markers. These molecular markers are powerful tools for basic and applied studies such as molecular breeding. In the model plants in Nicotiana genus e.g. N. benthamiana, a comprehensive assessment of SSR content has become possible now because several Nicotiana genomes have been sequenced. We conducted a genome-wide SSR characterization and marker development across seven Nicotiana genomes. Results Here, we initially characterized 2,483,032 SSRs (repeat units of 1–10 bp) from seven genomic sequences of Nicotiana and developed SSR markers using the GMATA® software package. Of investigated repeat units, mono-, di- and tri-nucleotide SSRs account for 98% of all SSRs in Nicotiana. More complex SSR motifs, although rare, are highly variable between Nicotiana genomes. A total of 1,224,048 non-redundant Nicotiana (NIX) markers were developed, of which 99.98% are novel. An efficient and uniform genotyping protocol for NIX markers was developed and validated. We created a web-based database of NIX marker information including amplicon sizes of alleles in each genome for downloading and online analysis. Conclusions The present work constitutes the first deep characterization of SSRs in seven genomes of Nicotiana, and the development of NIX markers for these SSRs. Our online marker database and an efficient genotyping protocol facilitate the application of these markers. The NIX markers greatly expand Nicotiana marker resources, thus providing a useful tool for future research and breeding. We demonstrate a novel protocol for SSR marker development and utilization at the whole genome scale that can be applied to any lineage of organisms. The Tobacco Markers & Primers Database (TMPD) is available at http://biodb.sdau.edu.cn/tmpd/index.html Electronic supplementary material The online version of this article (10.1186/s12864-018-4878-4) contains supplementary material, which is available to authorized users.

Background The central tenet of ‘genome content’ has been that the ‘non-coding’ parts are highly enriched with ‘microsatellites’ or ‘Simple Sequence Repeats’ (SSRs). We presume that the presence and change in number of repeat unit (n) of SSRs in different genomic locations may or may not become beneficial, depending on the position of SSRs in a gene. Very few studies have looked into the existence of SSRs in the hair-pin precursors of miRNAs (pre-miRNAs). The interplay between SSRs and miRNAs is not yet clearly understood. Results Considering the potential significance of SSRs in pre-miRNAs, we analysed the miRNA hair-pin precursors of 171 organisms, which revealed a noticeable (29.8%) existence of SSRs in their pre-miRNAs. The maintenance of SSRs in pre-miRNAs even in the complex, highly evolved phyla like Chordata and Magnoliophyta shed light upon its diverse functions. Putative effects of SSRs in either regulating the biogenesis or function of miRNAs were more underlined based on computational and experimental analysis. A preliminary computational analysis to explore the relevance of such SSRs maintained in pre-miRNA sequences led to the detection of splicing regulatory elements (SREs) either in or near to the SSRs. The absence of SSRs correspondingly decreased the detection of SREs. Conclusion The present study is the first implication for the possible involvement of SSRs in shaping the SREs to undergo Alternative Splicing events to produce miRNA isoforms in accordance with different stress environments. This part of work well demonstrates the importance of studying such consistently maintained SSRs residing in pre-miRNAs and can enhance more and more research towards deciphering the exact function of SSRs in the near future. Electronic supplementary material The online version of this article (10.1186/s12863-018-0615-x) contains supplementary material, which is available to authorized users.

Background Chromocenters are defined as a punctate condensed blocks of chromatin in the interphase cell nuclei of certain cell types with unknown biological significance. In recent years a progress in revealing of chromocenters protein content has been made although the details of DNA content within constitutive heterochromatin still remain unclear. It is known that these regions are enriched in tandem repeats (TR) and transposable elements. Quick improvement of genome sequencing does not help to assemble the heterochromatic regions due to lack of appropriate bioinformatics techniques. Results Chromocenters DNA have been isolated by a biochemical approach from mouse liver cells nuclei and sequenced on the Illumina MiSeq resulting in ChrmC dataset. Analysis of ChrmC dataset by the bioinformatics tools available revealed that the major component of chromocenter DNA are TRs: ~ 66% MaSat and ~ 4% MiSat. Other previously classified TR families constitute ~ 1% of ChrmC dataset. About 6% of chromocenters DNA are mostly unannotated sequences. In the contigs assembled with IDBA_UD there are many fragments of heterochromatic Y-chromosome, rDNA and other pseudo-genes and non-coding DNA. A protein coding sfi1 homolog gene fragment was also found in contigs. The Sfi1 homolog gene is located on the chromosome 11 in the reference genome very close to the Golden Pass Gap (a ~ 3 Mb empty region reserved to the pericentromeric region) and proves the purity of chromocenters isolation. The second major fraction are non-LTR retroposons (SINE and LINE) with overwhelming majority of LINE - ~ 11% of ChrmC. Most of the LINE fragments are from the ~ 2 kb region at the end of the 2nd ORF and its’ flanking region. The precise LINEs’ segment of ~ 2 kb is the necessary mouse constitutive heterohromatin component together with TR. The third most abundant fraction are ERVs. The ERV distribution in chromocenters differs from the whole genome: IAP (ERV2 class) is the most numerous in ChrmC while MaLR (ERV3 class) prevails in the reference genome. IAP and its LTR also prevail in TR containing contigs extracted from the WGS dataset. In silico prediction of IAP and LINE fragments in chromocenters was confirmed by direct fluorescent in situ hybridization (FISH). Conclusion Our data of chromocenters’ DNA (ChrmC) sequencing demonstrate that IAP with LTR and a precise ~ 2 kb fragment of LINE represent a substantial fraction of mouse chromocenters (constitutive heteroсhromatin) along with TRs. Electronic supplementary material The online version of this article (10.1186/s12864-018-4534-z) contains supplementary material, which is available to authorized users.

Fil: Milani, Diogo. Universidad Estatal Paulista Julio de Mesquita Filho. Instituto de Biociencias. Departamento de Biología; Brasil. Fil: Ramos, Érica. Universidad Estatal Paulista Julio de Mesquita Filho. Instituto de Biociencias. Departamento de Morfología; Brasil. Fil: Loreto, Vilma. Universidad Federal de Pernambuco. Centro de Biociencias. Departamento de Genética; Brasil. Fil: Marti, Dardo Andrea. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico (Nordeste). Instituto de Biología Subtropical; Argentina. Fil: Marti, Dardo Andrea. Universidad Nacional de Misiones. Facultad de Ciencias Exactas, Química y Naturales. Instituto de Biología Subtropical; Argentina. Fil: Lima Cardoso, Adauto. Universidad Estatal Paulista Julio de Mesquita Filho. Instituto de Biociencias. Departamento de Morfología; Brasil. Fil: Martinez de Moraes, Karen Cristiane. Universidad Estatal Paulista Julio de Mesquita Filho. Instituto de Biociencias. Departamento de Biología; Brasil. Fil: Martins, César. Universidad Estatal Paulista Julio de Mesquita Filho. Instituto de Biociencias. Departamento de Morfología; Brasil. Fil: Cabral de Mello, Diogo Cavalcanti. Universidad Estatal Paulista Julio de Mesquita Filho. Instituto de Biociencias. Departamento de Biología; Brasil. Background: Satellite DNAs (satDNAs) are organized in repetitions directly contiguous to one another, forming long arrays and composing a large portion of eukaryote genomes. These sequences evolve according to the concerted evolution model, and homogenization of repeats is observed at the intragenomic level. Satellite DNAs are the primary component of heterochromatin, located primarily in centromeres and telomeres. Moreover, satDNA enrichment in specific chromosomes has been observed, such as in B chromosomes, that can provide clues about composition, origin and evolution of this chromosome. In this study, we isolated and characterized a satDNA in A and B chromosomes of Abracris flavolineata by integrating cytogenetic, molecular and genomics approaches at intra- and inter-population levels, with the aim to understand the evolution of satDNA and composition of B chromosomes. Results: AflaSAT-1 satDNA was shared with other species and in A. flavolineata, was associated with another satDNA, AflaSAT-2. Chromosomal mapping revealed centromeric blocks variable in size in almost all chromosomes (except pair 11) of A complement for both satDNAs, whereas for B chromosome, only a small centromeric signal occurred. In distinct populations, variable number of AflaSAT-1 chromosomal sites correlated with variability in copy number. Instead of such variability, low sequence diversity was observed in A complement, but monomers from B chromosome were more variable, presenting also exclusive mutations. AflaSAT-1 was transcribed in five tissues of adults in distinct life cycle phases. Conclusions: The sharing of AflaSAT-1 with other species is consistent with the library hypothesis and indicates common origin in a common ancestor; however, AflaSAT-1 was highly amplified in the genome of A. flavolineata. At the population level, homogenization of repeats in distinct populations was documented, but dynamic expansion or elimination of repeats was also observed. Concerning the B chromosome, our data provided new information on the composition in A. flavolineata. Together with previous results, the sequences of heterochromatic nature were not likely highly amplified in the entire B chromosome. Finally, the constitutive transcriptional activity suggests a possible unknown functional role, which should be further investigated.

Background Identification of the types of strains of Mycobacterium tuberculosis (M. tuberculosis) complex causing tuberculosis (TB) could contribute to TB control program of specific geographic region as well as it could add knowledge onto the existing literature on TB worldwide. The objective of the present study was to identify the species and strains of M. tuberculosis complex causing pulmonary tuberculosis in central Ethiopia. Methods A health institution- based cross-sectional study was conducted on 338 smear positive TB cases visiting three hospitals between October 2012 and September 2013. Morning and spot sputum samples were collected before the starting of treatment regimens. Thus, a total of 338 pooled sputum samples collected from these cases. Samples were cultured on Löwenstein Jensen media and the isolates were identified by the region of difference (RD) 9 based polymerase chain reaction (PCR) and spoligotyping. Result Of the total isolates 98.6% of the isolates were identified to be M. tuberculosis while the remaining 1.4% were identified as M. africanum. Further, typing of M. tuberculosis using spoligotyping lead to the identification of 90 different strains of M. tuberculosis. Of these strains, 32 were clustered consisting of more than one isolate while the remaining 58 strains were unique consisting of single isolate. Thus, 79.3% (223/281) of the isolates were found in the clustered while only 20.6% (58/281) of the strains were unique. Forty-five of the spolgotyping patterns were registeredin the SITVIT2 or SpolDB4 database in while the remaining 45 were notfound in the database and hence were orphan strains. The dominant strains were SIT53, SIT149, and SIT54, consisting of 43, 37 and 34 isolates, respectively. Classification of the spoligotype patterns using TB-insight RUN TB-Lineage showed that 86.8, 6.4, 5, 1.4% ofthe isolatesbelonged to the Euro-American lineage, East-African-Indian, Indo-oceanic and M. africanum, respectively. Conclusion The identification of clustered and new strains using spolygotyping in present study does not give conclusive finding as spoligotyping has low discriminatory power. Thus, further identification of these isolates using mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VENTR) and or whole genome sequencing (WGS) recommended.

Background A prominent and distinctive feature of the rye (Secale cereale) chromosomes is the presence of massive blocks of subtelomeric heterochromatin, the size of which is correlated with the copy number of tandem arrays. The rapidity with which these regions have formed over the period of speciation remains unexplained. Results Using a BAC library created from the short arm telosome of rye chromosome 1R we uncovered numerous arrays of the pSc200 and pSc250 tandem repeat families which are concentrated in subtelomeric heterochromatin and identified the adjacent DNA sequences. The arrays show significant heterogeneity in monomer organization. 454 reads were used to gain a representation of the expansion of these tandem repeats across the whole rye genome. The presence of multiple, relatively short monomer arrays, coupled with the mainly star-like topology of the monomer phylogenetic trees, was taken as indicative of a rapid expansion of the pSc200 and pSc250 arrays. The evolution of subtelomeric heterochromatin appears to have included a significant contribution of illegitimate recombination. The composition of transposable elements (TEs) within the regions flanking the pSc200 and pSc250 arrays differed markedly from that in the genome a whole. Solo-LTRs were strongly enriched, suggestive of a history of active ectopic exchange. Several DNA motifs were over-represented within the LTR sequences. Conclusion The large blocks of subtelomeric heterochromatin have arisen from the combined activity of TEs and the expansion of the tandem repeats. The expansion was likely based on a highly complex network of recombination mechanisms. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2667-5) contains supplementary material, which is available to authorized users.

Background Listeria (L.) monocytogenes causes fatal infections in many species including ruminants and humans. In ruminants, rhombencephalitis is the most prevalent form of listeriosis. Using multilocus variable number tandem repeat analysis (MLVA) we recently showed that L. monocytogenes isolates from ruminant rhombencephalitis cases are distributed over three genetic complexes (designated A, B and C). However, the majority of rhombencephalitis strains and virtually all those isolated from cattle cluster in MLVA complex A, indicating that strains of this complex may have increased neurotropism and neurovirulence. The aim of this study was to investigate whether ruminant rhombencephalitis strains have an increased ability to propagate in the bovine hippocampal brain-slice model and can be discriminated from strains of other sources. For this study, forty-seven strains were selected and assayed on brain-slice cultures, a bovine macrophage cell line (BoMac) and a human colorectal adenocarcinoma cell line (Caco-2). They were isolated from ruminant rhombencephalitis cases (n = 21) and other sources including the environment, food, human neurolisteriosis cases and ruminant/human non-encephalitic infection cases (n = 26). Results All but one L. monocytogenes strain replicated in brain slices, irrespectively of the source of the isolate or MLVA complex. The replication of strains from MLVA complex A was increased in hippocampal brain-slice cultures compared to complex C. Immunofluorescence revealed that microglia are the main target cells for L. monocytogenes and that strains from MLVA complex A caused larger infection foci than strains from MLVA complex C. Additionally, they caused larger plaques in BoMac cells, but not CaCo-2 cells. Conclusions Our brain slice model data shows that all L. monocytogenes strains should be considered potentially neurovirulent. Secondly, encephalitis strains cannot be conclusively discriminated from non-encephalitis strains with the bovine organotypic brain slice model. The data indicates that MLVA complex A strains are particularly adept at establishing encephalitis possibly by virtue of their higher resistance to antibacterial defense mechanisms in microglia cells, the main target of L. monocytogenes. Electronic supplementary material The online version of this article (doi:10.1186/s12866-015-0454-0) contains supplementary material, which is available to authorized users.

Background Subtelomeres are located close to the ends of chromosomes and organized by tandemly repetitive sequences, duplicated copies of genes, pseudogenes and retrotransposons. Transcriptional activity of tandemly organized DNA at terminal chromosomal regions and the distribution of subtelomere-derived non-coding RNAs are poorly investigated. Here we aimed to analyze transcriptional activity of subtelomeric tandem repeat in somatic tissues and cultured cells of birds. We focused on tissue-specific differences of subtelomeric repeats transcription, structure of the resulting transcripts and the behavior of subtelomere-derived RNA during mitosis. Results Transcriptional activity of short subtelomeric PO41 (“pattern of 41 bp”) tandem repeat in the somatic and cultured cells of chicken (Gallus gallus domesticus) and Japanese quail (Coturnix coturnix japonica) was examined using RNA fluorescence in situ hybridization approach. We discovered transcripts from both strands of the PO41 repeat in chicken MDCC-MSB1 cells as well as in chicken and Japanese quail somatic tissues, such as tissues of cerebellum, telencephalon, muscles, oviduct, small and large intestine. Normal somatic and transformed cells demonstrate similar distribution of PO41 repeat transcripts in interphase nuclei. We revealed one or two major foci of PO41 repeat transcripts associated with RNA polymerase II, representing nascent RNA, and dispersed PO41 repeat transcripts localized in euchromatin or interchromatin space, representing released RNA. During mitosis PO41 non-coding RNA distribute between condensed chromosomes till anaphase, when they concentrate at the cleavage plane. At telophase, clusters of PO41 RNA surround terminal regions of chromosomes. Treatments with RNases of different substrate specificity indicate that PO41 repeat transcripts are single-stranded RNAs with short double-stranded regions due to appearance of inverted repeats. Conclusion Transcription of a subtelomeric tandem repeat in avian somatic cells is reported here for the first time. PO41 repeat transcription is conserved among Galliformes and has similar pattern in somatic tissues. We demonstrated redistribution of non-coding PO41 RNA occurring during the cell cycle. Potential regulatory role of the PO41 repeat transcripts in RNA-dependent process of subtelomere heterochromatin maintenance is discussed. Electronic supplementary material The online version of this article (doi:10.1186/s13039-014-0102-7) contains supplementary material, which is available to authorized users.

Background Perennial ryegrass (Lolium perenne L.) is one of the most important forage and turf grass species of temperate regions worldwide. Its mitochondrial genome is inherited maternally and contains genes that can influence traits of agricultural importance. Moreover, the DNA sequence of mitochondrial genomes has been established and compared for a large number of species in order to characterize evolutionary relationships. Therefore, it is crucial to understand the organization of the mitochondrial genome and how it varies between and within species. Here, we report the first de novo assembly and annotation of the complete mitochondrial genome from perennial ryegrass. Results Intact mitochondria from perennial ryegrass leaves were isolated and used for mtDNA extraction. The mitochondrial genome was sequenced to a 167-fold coverage using the Roche 454 GS-FLX Titanium platform, and assembled into a circular master molecule of 678,580 bp. A total of 34 proteins, 14 tRNAs and 3 rRNAs are encoded by the mitochondrial genome, giving a total gene space of 48,723 bp (7.2%). Moreover, we identified 149 open reading frames larger than 300 bp and covering 67,410 bp (9.93%), 250 SSRs, 29 tandem repeats, 5 pairs of large repeats, and 96 pairs of short inverted repeats. The genes encoding subunits of the respiratory complexes – nad1 to nad9, cob, cox1 to cox3 and atp1 to atp9 – all showed high expression levels both in absolute numbers and after normalization. Conclusions The circular master molecule of the mitochondrial genome from perennial ryegrass presented here constitutes an important tool for future attempts to compare mitochondrial genomes within and between grass species. Our results also demonstrate that mitochondria of perennial ryegrass contain genes crucial for energy production that are well conserved in the mitochondrial genome of monocotyledonous species. The expression analysis gave us first insights into the transcriptome of these mitochondrial genes in perennial ryegrass.

8 pages, 3 figures, 1 table, 1 additional file., Background: RNA-binding proteins of the PUF family share a conserved domain consisting of tandemly repeated 36-40 amino acid motifs (typically eight) known as Puf repeats. Proteins containing tandem repeats are often dominant targets of humoral responses during infectious diseases. Thus, we considered of interest to analyze whether Leishmania PUF proteins result antigenic during visceral leishmaniasis (VL)., Findings: Here, employing whole-genome databases, we report the composition, and structural features, of the PUF family in Leishmania infantum. Additionally, the 10 genes of the L. infantum PUF family were cloned and used to express the Leishmania PUFs in bacteria as recombinant proteins. Finally, the antigenicity of these PUF proteins was evaluated by determining levels of specific antibodies in sera from experimentally infected hamsters. The Leishmania PUFs were all recognized by the sera, even though with different degree of reactivity and/or frequency of recognition. The reactivity of hamster sera against recombinant LiPUF1 and LiPUF2 was particularly prominent, and these proteins were subsequently assayed against sera from human patients. High antibody responses against rLiPUF1 and rLiPUF2 were found in sera from VL patients, but these proteins resulted also recognized by sera from Chagas’ disease patients., Conclusion: Our results suggest that Leishmania PUFs are targets of the humoral response during L. infantum infection and may represent candidates for serodiagnosis and/or vaccine reagents; however, it should be kept in mind the cross-reactivity of LiPUFs with antibodies induced against other trypanosomatids such as Trypanosoma cruzi., This work was funded by grants from the Ministerio de Ciencia y Tecnología (BFU2006-08346), and the Spanish Ministry of Science and Innovation and the Instituto de Salud Carlos III within the Network of Tropical Diseases Research (RICET RD06/0021/0008 - FEDER). Also, an institutional grant from Fundación Ramón Areces is acknowledged.

Background Expressed sequence tag (EST) databases provide a valuable source of genetic data in organisms whose genome sequence information is not yet compiled. We used a published EST database for the waterflea Daphnia magna (Crustacea:Cladocera) to isolate variable number of tandem repeat (VNTR) markers for linkage mapping, Quantitative Trait Loci (QTL), and functional studies. Findings Seventy-four polymorphic markers were isolated and characterised. Analyses of repeat structure, putative gene function and polymorphism indicated that intragenic tandem repeats are not distributed randomly in the mRNA sequences; instead, dinucleotides are more frequent in non-coding regions, whereas trinucleotides (and longer motifs involving multiple-of-three nucleotide repeats) are preferentially situated in coding regions. We also observed differential distribution of repeat motifs across putative genetic functions. This indicates differential selective constraints and possible functional significance of VNTR polymorphism in at least some genes. Conclusion Databases of VNTR markers situated in genes whose putative function can be inferred from homology searches will be a valuable resource for the genetic study of functional variation and selection.

Background Halibuts are commercially important flatfish species confined to the North Pacific and North Atlantic Oceans. We have determined the complete mitochondrial genome sequences of four specimens each of Atlantic halibut (Hippoglossus hippoglossus), Pacific halibut (Hippoglossus stenolepis) and Greenland halibut (Reinhardtius hippoglossoides), and assessed the nucleotide variability within and between species. Results About 100 variable positions were identified within the four specimens in each halibut species, with the control regions as the most variable parts of the genomes (10 times that of the mitochondrial ribosomal DNA). Due to tandem repeat arrays, the control regions have unusually large sizes compared to most vertebrate mtDNAs. The arrays are highly heteroplasmic in size and consist mainly of different variants of a 61-bp motif. Halibut mitochondrial genomes lacking arrays were also detected. Conclusion The complexity, distribution, and biological role of the heteroplasmic tandem repeat arrays in halibut mitochondrial control regions are discussed. We conclude that the most plausible explanation for array maintenance includes both the slipped-strand mispairing and DNA recombination mechanisms.

Background The evolution of individual microsatellite loci is often complex and homoplasy is common but often goes undetected. Sequencing alleles at a microsatellite locus can provide a more complete picture of the common evolutionary mechanisms occurring at that locus and can reveal cases of homoplasy. Within species homoplasy can lead to an underestimate of differentiation among populations and among species homoplasy can produce a misleading interpretation regarding shared alleles and hybridization. This is especially problematic with cryptic species. Results By sequencing alleles from three cryptic species of the sweetpotato whitefly (Bemisia tabaci), designated MEAM1, MED, and NW, the evolution of the putatively dinucleotide Bem6 (CA8)imp microsatellite locus is inferred as one of primarily stepwise mutation occurring at four distinct heptaucleotide tandem repeats. In two of the species this pattern yields a compound tandem repeat. Homoplasy was detected both among species and within species. Conclusions In the absence of sequencing, size homoplasious alleles at the Bem6 locus lead to an overestimate of alleles shared and hybridization among cryptic species of Bemisia tabaci. Furthermore, the compound heptanucleotide motif structure of a putative dinucleotide microsatellite has implications for the nomenclature of heptanucleotide tandem repeats with step-wise evolution.

Background H. pylori (Helicobacter pylori) is the major causative agent of chronic active gastritis. The population of H. pylori shows a high genomic variability among isolates. And the polymorphism of repeat-units of genomics had participated the important process of evolution. Its long term colonization of the stomach caused different clinical outcomes, which may relate to the high degree of genetic variation of H. pylori. A variety of molecular typing tools have been developed to access genetic relatedness in H. pylori isolates. However, there is still no standard genotyping system of this bacterium. The MLVA (Multi-locus of variable number of tandem repeat analysis) method is useful for performing phylogenetic analysis and is widely used in bacteria genotyping; however, there's little application in H. pylori analysis. This article is the first application of the MLVA method to investigate H. pylori from different districts and ethnic groups of China. Results MLVA of 12 VNTR loci with high discrimination power based on 30 candidates were performed on a collection of 202 strains of H. pylori which originated from five regions of China and Japan. Phylogenetic tree was constructed using MLVA profiles. 12 VNTR loci presented with high various polymorphisms, and the results demonstrated very close relationships between genotypes and ethnic groups. Conclusions This study used MLVA methodology providing a new perspective on the ethnic groups and distribution characteristics of H. pylori.

Background Multiple-locus variable-number tandem repeat (VNTR) analysis (MLVA) is a prominent subtyping method to resolve closely related microbial isolates to provide information for establishing genetic patterns among isolates and to investigate disease outbreaks. The usefulness of MLVA was recently demonstrated for the avian major pathogen Chlamydophila psittaci. In the present study, we developed a similar method for another pathogen of birds: the filamentous fungus Aspergillus fumigatus. Results We selected 10 VNTR markers located on 4 different chromosomes (1, 5, 6 and 8) of A. fumigatus. These markers were tested with 57 unrelated isolates from different hosts or their environment (53 isolates from avian species in France, China or Morocco, 3 isolates from humans collected at CHU Henri Mondor hospital in France and the reference strain CBS 144.89). The Simpson index for individual markers ranged from 0.5771 to 0.8530. A combined loci index calculated with all the markers yielded an index of 0.9994. In a second step, the panel of 10 markers was used in different epidemiological situations and tested on 277 isolates, including 62 isolates from birds in Guangxi province in China, 95 isolates collected in two duck farms in France and 120 environmental isolates from a turkey hatchery in France. A database was created with the results of the present study http://minisatellites.u-psud.fr/MLVAnet/. Three major clusters of isolates were defined by using the graphing algorithm termed Minimum Spanning Tree (MST). The first cluster comprised most of the avian isolates collected in the two duck farms in France, the second cluster comprised most of the avian isolates collected in poultry farms in China and the third one comprised most of the isolates collected in the turkey hatchery in France. Conclusions MLVA displayed excellent discriminatory power. The method showed a good reproducibility. MST analysis revealed an interesting clustering with a clear separation between isolates according to their geographic origin rather than their respective hosts.

Background About forty human diseases are caused by repeat instability mutations. A distinct subset of these diseases is the result of extreme expansions of polymorphic trinucleotide repeats; typically CAG repeats encoding poly-glutamine (poly-Q) tracts in proteins. Polymorphic repeat length variation is also apparent in human poly-Q encoding genes from normal individuals. As these coding sequence repeats are subject to selection in mammals, it has been suggested that normal variations in some of these typically highly conserved genes are implicated in morphological differences between species and phenotypic variations within species. At present, poly-Q encoding genes in non-human mammalian species are poorly documented, as are their functions and propensities for polymorphic variation. Results The current investigation identified 178 bovine poly-Q encoding genes (Q ≥ 5) and within this group, 26 genes with orthologs in both human and mouse that did not contain poly-Q repeats. The bovine poly-Q encoding genes typically had ubiquitous expression patterns although there was bias towards expression in epithelia, brain and testes. They were also characterised by unusually large sizes. Analysis of gene ontology terms revealed that the encoded proteins were strongly enriched for functions associated with transcriptional regulation and many contributed to physical interaction networks in the nucleus where they presumably act cooperatively in transcriptional regulatory complexes. In addition, the coding sequence CAG repeats in some bovine genes impacted mRNA splicing thereby generating unusual transcriptional diversity, which in at least one instance was tissue-specific. The poly-Q encoding genes were prioritised using multiple criteria for their likelihood of being polymorphic and then the highest ranking group was experimentally tested for polymorphic variation within a cattle diversity panel. Extensive and meiotically stable variation was identified. Conclusions Transcriptional diversity can potentially be generated in poly-Q encoding genes by the impact of CAG repeat tracts on mRNA alternative splicing. This effect, combined with the physical interactions of the encoded proteins in large transcriptional regulatory complexes suggests that polymorphic variations of proteins in these complexes have strong potential to affect phenotype.

Background In most protein-coding genes, greater sequence variation is observed in noncoding regions (introns and untranslated regions) than in coding regions due to selective constraints. During characterization of genes and transcripts encoding small secreted salivary gland proteins (SSSGPs) from the Hessian fly, we found exactly the opposite pattern of conservation in several families of genes: the non-coding regions were highly conserved, but the coding regions were highly variable. Results Seven genes from the SSSGP-1 family are clustered as one inverted and six tandem repeats within a 15 kb region of the genome. Except for SSSGP-1A2, a gene that encodes a protein identical to that encoded by SSSGP-1A1, the other six genes consist of a highly diversified, mature protein-coding region as well as highly conserved regions including the promoter, 5'- and 3'-UTRs, a signal peptide coding region, and an intron. This unusual pattern of highly diversified coding regions coupled with highly conserved regions in the rest of the gene was also observed in several other groups of SSSGP-encoding genes or cDNAs. The unusual conservation pattern was also found in some of the SSSGP cDNAs from the Asian rice gall midge, but not from the orange wheat blossom midge. Strong positive selection was one of the forces driving for diversification whereas concerted homogenization was likely a mechanism for sequence conservation. Conclusion Rapid diversification in mature SSSGPs suggests that the genes are under selection pressure for functional adaptation. The conservation in the noncoding regions of these genes including introns also suggested potential mechanisms for sequence homogenization that are not yet fully understood. This report should be useful for future studies on genetic mechanisms involved in evolution and functional adaptation of parasite genes.

Background Shigella flexneri is one of the causative agents of shigellosis, a major cause of childhood mortality in developing countries. Multilocus variable-number tandem repeat (VNTR) analysis (MLVA) is a prominent subtyping method to resolve closely related bacterial isolates for investigation of disease outbreaks and provide information for establishing phylogenetic patterns among isolates. The present study aimed to develop an MLVA method for S. flexneri and the VNTR loci identified were tested on 242 S. flexneri isolates to evaluate their variability in various serotypes. The isolates were also analyzed by pulsed-field gel electrophoresis (PFGE) to compare the discriminatory power and to evaluate the usefulness of MLVA as a tool for phylogenetic analysis of S. flexneri. Results Thirty-six VNTR loci were identified by exploring the repeat sequence loci in genomic sequences of Shigella species and by testing the loci on nine isolates of different subserotypes. The VNTR loci in different serotype groups differed greatly in their variability. The discriminatory power of an MLVA assay based on four most variable VNTR loci was higher, though not significantly, than PFGE for the total isolates, a panel of 2a isolates, which were relatively diverse, and a panel of 4a/Y isolates, which were closely-related. Phylogenetic groupings based on PFGE patterns and MLVA profiles were considerably concordant. The genetic relationships among the isolates were correlated with serotypes. The phylogenetic trees constructed using PFGE patterns and MLVA profiles presented two distinct clusters for the isolates of serotype 3 and one distinct cluster for each of the serotype groups, 1a/1b/NT, 2a/2b/X/NT, 4a/Y, and 6. Isolates that had different serotypes but had closer genetic relatedness than those with the same serotype were observed between serotype Y and subserotype 4a, serotype X and subserotype 2b, subserotype 1a and 1b, and subserotype 3a and 3b. Conclusions The 36 VNTR loci identified exhibited considerably different degrees of variability among S. flexneri serotype groups. VNTR locus could be highly variable in a serotype but invariable in others. MLVA assay based on four highly variable loci could display a comparable resolving power to PFGE in discriminating isolates. MLVA is also a prominent molecular tool for phylogenetic analysis of S. flexneri; the resulting data are beneficial to establish clear clonal patterns among different serotype groups and to discern clonal groups among isolates within the same serotype. As highly variable VNTR loci could be serotype-specific, a common MLVA protocol that consists of only a small set of loci, for example four to eight loci, and that provides high resolving power to all S. flexneri serotypes may not be obtainable.

Background Alpha-isopropylmalate synthase (α-IPMS) is the key enzyme that catalyzes the first committed step in the leucine biosynthetic pathway. The gene encoding α-IPMS in Mycobacterium tuberculosis, leuA, is polymorphic due to the insertion of 57-bp repeat units referred to as Variable Number of Tandem Repeats (VNTR). The role of the VNTR found within the M. tuberculosis genome is unclear. To investigate the role of the VNTR in leuA, we compared two α-IPMS proteins with different numbers of amino acid repeats, one with two copies and the other with 14 copies. We have cloned leuA with 14 copies of the repeat units into the pET15b expression vector with a His6-tag at the N-terminus, as was previously done for the leuA gene with two copies of the repeat units. Results The recombinant His6-α-IPMS proteins with two and 14 copies (α-IPMS-2CR and α-IPMS-14CR, respectively) of the repeat units were purified by immobilized metal ion affinity chromatography and gel filtration. Both enzymes were found to be dimers by gel filtration. Both enzymes work well at pH values of 7–8.5 and temperatures of 37–42°C. However, α-IPMS-14CR tolerates pH values and temperatures outside of this range better than α-IPMS-2CR does. α-IPMS-14CR has higher affinity than α-IPMS-2CR for the two substrates, α-ketoisovalerate and acetyl CoA. Furthermore, α-IPMS-2CR was feedback inhibited by the end product l-leucine, whereas α-IPMS-14CR was not. Conclusion The differences in the kinetic properties and the l-leucine feedback inhibition between the two M. tuberculosis α-IPMS proteins containing low and high numbers of VNTR indicate that a large VNTR insertion affects protein structure and function. Demonstration of l-leucine binding to α-IPMS-14CR would confirm whether or not α-IPMS-14CR responds to end-product feedback inhibition.

A genomic analysis of the annual fish Nothobranchius furzeri, a vertebrate with the shortest known life span in captivity and which may provide a new model organism for aging research., Background The annual fish Nothobranchius furzeri is the vertebrate with the shortest known life span in captivity. Fish of the GRZ strain live only three to four months under optimal laboratory conditions, show explosive growth, early sexual maturation and age-dependent physiological and behavioral decline, and express aging related biomarkers. Treatment with resveratrol and low temperature significantly extends the maximum life span. These features make N. furzeri a promising new vertebrate model for age research. Results To contribute to establishing N. furzeri as a new model organism, we provide a first insight into its genome and a comparison to medaka, stickleback, tetraodon and zebrafish. The N. furzeri genome contains 19 chromosomes (2n = 38). Its genome of between 1.6 and 1.9 Gb is the largest among the analyzed fish species and has, at 45%, the highest repeat content. Remarkably, tandem repeats comprise 21%, which is 4-12 times more than in the other four fish species. In addition, G+C-rich tandem repeats preferentially localize to centromeric regions. Phylogenetic analysis based on coding sequences identifies medaka as the closest relative. Genotyping of an initial set of 27 markers and multi-locus fingerprinting of one microsatellite provides the first molecular evidence that the GRZ strain is highly inbred. Conclusions Our work presents a first basis for systematic genomic and genetic analyses aimed at understanding the mechanisms of life span determination in N. furzeri.

Background The complete proteome of the starlet sea anemone, Nematostella vectensis, provides insights into gene invention dating back to the Cnidarian-Bilaterian ancestor. With the addition of the complete proteomes of Hydra magnipapillata and Monosiga brevicollis, the investigation of proteins having unique features in early metazoan life has become practical. We focused on the properties and the evolutionary trends of tandem repeat (TR) sequences in Cnidaria proteomes. Results We found that 11-16% of N. vectensis proteins contain tandem repeats. Most TRs cover 150 amino acid segments that are comprised of basic units of 5-20 amino acids. In total, the N. Vectensis proteome has about 3300 unique TR-units, but only a small fraction of them are shared with H. magnipapillata, M. brevicollis, or mammalian proteomes. The overall abundance of these TRs stands out relative to that of 14 proteomes representing the diversity among eukaryotes and within the metazoan world. TR-units are characterized by a unique composition of amino acids, with cysteine and histidine being over-represented. Structurally, most TR-segments are associated with coiled and disordered regions. Interestingly, 80% of the TR-segments can be read in more than one open reading frame. For over 100 of them, translation of the alternative frames would result in long proteins. Most domain families that are characterized as repeats in eukaryotes are found in the TR-proteomes from Nematostella and Hydra. Conclusions While most TR-proteins have originated from prediction tools and are still awaiting experimental validations, supportive evidence exists for hundreds of TR-units in Nematostella. The existence of TR-proteins in early metazoan life may have served as a robust mode for novel genes with previously overlooked structural and functional characteristics.

Background Identification of approximate tandem repeats is an important task of broad significance and still remains a challenging problem of computational genomics. Often there is no single best approach to periodicity detection and a combination of different methods may improve the prediction accuracy. Discrete Fourier transform (DFT) has been extensively used to study primary periodicities in DNA sequences. Here we investigate the application of DFT method to identify and study alphoid higher order repeats. Results We used method based on DFT with mapping of symbolic into numerical sequence to identify and study alphoid higher order repeats (HOR). For HORs the power spectrum shows equidistant frequency pattern, with characteristic two-level hierarchical organization as signature of HOR. Our case study was the 16 mer HOR tandem in AC017075.8 from human chromosome 7. Very long array of equidistant peaks at multiple frequencies (more than a thousand higher harmonics) is based on fundamental frequency of 16 mer HOR. Pronounced subset of equidistant peaks is based on multiples of the fundamental HOR frequency (multiplication factor n for nmer) and higher harmonics. In general, nmer HOR-pattern contains equidistant secondary periodicity peaks, having a pronounced subset of equidistant primary periodicity peaks. This hierarchical pattern as signature for HOR detection is robust with respect to monomer insertions and deletions, random sequence insertions etc. For a monomeric alphoid sequence only primary periodicity peaks are present. The 1/fβ – noise and periodicity three pattern are missing from power spectra in alphoid regions, in accordance with expectations. Conclusion DFT provides a robust detection method for higher order periodicity. Easily recognizable HOR power spectrum is characterized by hierarchical two-level equidistant pattern: higher harmonics of the fundamental HOR-frequency (secondary periodicity) and a subset of pronounced peaks corresponding to constituent monomers (primary periodicity). The number of lower frequency peaks (secondary periodicity) below the frequency of the first primary periodicity peak reveals the size of nmer HOR, i.e., the number n of monomers contained in consensus HOR.

Comparative genomic analysis of the Impact locus, which is imprinted in Glires but not in other mammals, reveals features required for genomic imprinting., Background Imprinted genes are exclusively expressed from one of the two parental alleles in a parent-of-origin-specific manner. In mammals, nearly 100 genes are documented to be imprinted. To understand the mechanism behind this gene regulation and to identify novel imprinted genes, common features of DNA sequences have been analyzed; however, the general features required for genomic imprinting have not yet been identified, possibly due to variability in underlying molecular mechanisms from locus to locus. Results We performed a thorough comparative genomic analysis of a single locus, Impact, which is imprinted only in Glires (rodents and lagomorphs). The fact that Glires and primates diverged from each other as recent as 70 million years ago makes comparisons between imprinted and non-imprinted orthologues relatively reliable. In species from the Glires clade, Impact bears a differentially methylated region, whereby the maternal allele is hypermethylated. Analysis of this region demonstrated that imprinting was not associated with the presence of direct tandem repeats nor with CpG dinucleotide density. In contrast, a CpG periodicity of 8 bp was observed in this region in species of the Glires clade compared to those of carnivores, artiodactyls, and primates. Conclusions We show that tandem repeats are dispensable, establishment of the differentially methylated region does not rely on G+C content and CpG density, and the CpG periodicity of 8 bp is meaningful to the imprinting. This interval has recently been reported to be optimal for de novo methylation by the Dnmt3a-Dnmt3L complex, suggesting its importance in the establishment of imprinting in Impact and other genes.

Background The homologous 4q and 10q subtelomeric regions include two distinctive polymorphic arrays of 3.3 kb repeats, named D4Z4. An additional BlnI restriction site on the 10q-type sequence allows to distinguish the chromosomal origin of the repeats. Reduction in the number of D4Z4 repeats below a threshold of 10 at the 4q locus is tightly linked to Facioscapulohumeral Muscular Dystrophy (FSHD), while similar contractions at 10q locus, are not pathogenic. Sequence variations due to the presence of BlnI-sensitive repeats (10q-type) on chromosome 4 or viceversa of BlnI-resistant repeats (4q-type) on chromosome 10 are observed in both alleles. Results We analysed DNA samples from 116 healthy subiects and 114 FSHD patients and determined the size distributions of polymorphic 4q and 10q alleles, the frequency and the D4Z4 repeat assortment of variant alleles, and finally the telomeric sequences both in standard and variant alleles. We observed the same frequency and types of variant alleles in FSHD patients and controls, but we found marked differences between the repeat arrays of the 4q and 10q chromosomes. In particular we detected 10q alleles completely replaced by the 4q subtelomeric region, consisting in the whole set of 4q-type repeats and the distal telomeric markers. However the reciprocal event, 10q-type subtelomeric region on chromosome 4, was never observed. At 4q locus we always identified hybrid alleles containing a mixture of 4q and 10q-type repeats. Conclusion The different size distribution and different structure of 10q variant alleles as compared with 4q suggests that these loci evolved in a different manner, since the 4q locus is linked to FSHD, while no inheritable disease is associated with mutations in 10qter genomic region. Hybrid alleles on chromosome 4 always retain a minimum number of 4q type repeats, as they are probably essential for maintaining the structural and functional properties of this subtelomeric region. In addition we found: i) several instances of variant alleles that could be misinterpreted and interfere with a correct diagnosis of FSHD; ii) the presence of borderline alleles in the range of 30–40 kb that carried a qA type telomere and were not associated with the disease.

The identification and analysis of repetitive patterns are active areas of biological and computational research. Tandem repeats in telomeres play a role in cancer and hypervariable trinucleotide tandem repeats are linked to over a dozen major neurodegenerative genetic disorders. In this paper, we present an algorithm to identify the exact and inexact repeat patterns in DNA sequences based on orthogonal exactly periodic subspace decomposition technique. Using the new measure our algorithm resolves the problems like whether the repeat pattern is of period P or its multiple (i.e., 2P, 3P, etc.), and several other problems that were present in previous signal-processing-based algorithms. We present an efficient algorithm of O(NL(w) log L(w)), where N is the length of DNA sequence and L(w) is the window length, for identifying repeats. The algorithm operates in two stages. In the first stage, each nucleotide is analyzed separately for periodicity, and in the second stage, the periodic information of each nucleotide is combined together to identify the tandem repeats. Datasets having exact and inexact repeats were taken up for the experimental purpose. The experimental result shows the effectiveness of the approach.

Background Tandem repeat arrays showing variation between sequences within a population, between strains or across species may have functional effects. The increasing availability of genomic sequence data makes routine description of observed variation possible, creating a need for tools to describe such variability. Results We present a set of programs that facilitate the identification of tandem repeats showing variation across multiple sequences or genomes, and the prediction of potentially polymorphic tandem repeats. The VNTRfinder (Variable Number of Tandem Repeats finder) program enables the detection of sequence length variation between arrays of inter-specific or intra-specific tandem repeats. In the absence of comparable sequences to explore observed variation, predictions are provided describing which tandem repeats are more likely to be variable, to help guide and focus further experimental evaluation. Conclusion These tools represent a resource for researchers interested in tandem repeats in nucleotide sequences that are most likely to be of clinical and evolutionary interest. The tools are available at http://bioinformatics.rcsi.ie/vntrfinder/. Downloadable versions for UNIX/LINUX and WINDOWS which permit the consideration of longer and more numerous sequences are also available.

Background The multilocus variable-number tandem repeat (VNTR) analysis (MLVA) technique has been developed for fine typing of many bacterial species. The genomic sequences of Neisseria meningitidis strains Z2491, MC58 and FAM18 have been available for searching potential VNTR loci by computer software. In this study, we developed and evaluated a MLVA method for molecular subtyping and phylogenetic analysis of N. meningitidis strains. Results A total of 12 VNTR loci were identified for subtyping and phylogenetic analysis of 100 N. meningitidis isolates, which had previously been characterized by pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing. The number of alleles ranges from 3 to 40 for the 12 VNTR loci; theoretically, the numbers of alleles can generate more than 5 × 1011 MLVA types. In total, 93 MLVA types were identified in the 100 isolates, indicating that MLVA is powerful in discriminating N. meningitidis strains. In phylogenetic analysis with the minimal spanning tree method, clonal relationships, established with MLVA types, agreed well with those built with ST types. Conclusion Our study indicates that the MLVA method has a higher degree of resolution than PFGE in discriminating N. meningitidis isolates and may be a useful tool for phylogenetic studies of strains evolving over different time scales.

Background The genome of Bacillus anthracis, the etiological agent of anthrax, is highly monomorphic which makes differentiation between strains difficult. A Multiple Locus Variable-number tandem repeats (VNTR) Analysis (MLVA) assay based on 20 markers was previously described. It has considerable discrimination power, reproducibility, and low cost, especially since the markers proposed can be typed by agarose-gel electrophoresis. However in an emergency situation, faster genotyping and access to representative databases is necessary. Results Genotyping of B. anthracis reference strains and isolates from France and Italy was done using a 25 loci MLVA assay combining 21 previously described loci and 4 new ones. DNA was amplified in 4 multiplex PCR reactions and the length of the resulting 25 amplicons was estimated by automated capillary electrophoresis. The results were reproducible and the data were consistent with other gel based methods once differences in mobility patterns were taken into account. Some alleles previously unresolved by agarose gel electrophoresis could be resolved by capillary electrophoresis, thus further increasing the assay resolution. One particular locus, Bams30, is the result of a recombination between a 27 bp tandem repeat and a 9 bp tandem repeat. The analysis of the array illustrates the evolution process of tandem repeats. Conclusion In a crisis situation of suspected bioterrorism, standardization, speed and accuracy, together with the availability of reference typing data are important issues, as illustrated by the 2001 anthrax letters event. In this report we describe an upgrade of the previously published MLVA method for genotyping of B. anthracis and apply the method to the typing of French and Italian B. anthracis strain collections. The increased number of markers studied compared to reports using only 8 loci greatly improves the discrimination power of the technique. An Italian strain belonging to the B branch was described, and two new branches, D and E, are proposed. Owing to the upgrading achieved here, precise genotyping can now be produced either by automated capillary electrophoresis, or by the more accessible but slower and for some markers slightly less accurate agarose gel methodology.

Summary The genomes of many species are dominated by short sequences repeated consecutively. It is estimated that over 10% of the human genome consists of tandemly repeated sequences. Finding repeated regions in long sequences is important in sequence analysis. We develop a software, LocRepeat, that finds regions of pseudo-periodic repeats in a long sequence. We use the definition of Li et al. 1 for the pseudo-periodic partition of a region and extend the algorithm that can select the repeated region from a given long sequence and give the pseudo-periodic partition of the region. Availability LocRepeat is available at http://www.cs.cityu.edu.hk/~lwang/software/LocRepeat