Molecular characterization of Mycobacterium tuberculosis isolates from Tehran, Iran by restriction fragment length polymorphism analysis and spoligotyping

Abstract: INTRODUCTION Characterization of Mycobacterium tuberculosis (MTB) isolates by DNA fingerprinting has contributed to tuberculosis (TB) control. The aim of this study was to determine the genetic diversity of MTB isolates from Tehran province in Iran. METHODS MTB isolates from 60 Iranian and 10 Afghan TB patients were fingerprinted by standard IS6110-restriction fragment length polymorphism (RFLP) analysis and spoligotyping. RESULTS The copy number of IS6110 ranged from 10-24 per isolate. The isolates were classified into 22 clusters showing ≥ 80% similarity by RFLP analysis. Fourteen multidrug-resistant (MDR) isolates were grouped into 4 IS6110-RFLP clusters, with 10 isolates [71% (95% CI: 45-89%)] in 1 cluster, suggesting a possible epidemiological linkage. Eighteen Iranian isolates showed ≥ 80% similarity with Afghan isolates. There were no strains with identical fingerprints. Spoligotyping of 70 isolates produced 23 distinct patterns. Sixty (85.7%) isolates were grouped into 13 clusters, while the remaining 10 isolates (14.2%) were not clustered. Ural (formerly Haarlem4) (n = 22, 31.4%) was the most common family followed by Central Asian strain (CAS) (n = 18, 25.7%) and T (n = 9, 12.8%) families. Only 1strain was characterized as having the Beijing genotype. Among 60 Iranian and 10 Afghan MTB isolates, 25% (95% CI: 16-37) and 70% (95% CI: 39-89) were categorized as Ural lineage, respectively. CONCLUSIONS A higher prevalence of Ural family MTB isolates among Afghan patients than among Iranian patients suggests the possible transmission of this lineage following the immigration of Afghans to Iran.

Results of An International Study of Quantitative BCR/ABL Assays Using Defined RNA Samples

Molecular monitoring of BCR/ABL transcripts by real time quantitative reverse transcription PCR (qRT-PCR) is an essential technique for clinical management of patients with BCR/ABL-positive CML and ALL. Though quantitative BCR/ABL assays are performed in hundreds of laboratories worldwide, results among these laboratories cannot be reliably compared due to heterogeneity in test methods, data analysis, reporting, and lack of quantitative standards. Recent efforts towards standardization have been limited in scope. Aliquots of RNA were sent to clinical test centers worldwide in order to evaluate methods and reporting for e1a2, b2a2, and b3a2 transcript levels using their own qRT-PCR assays. Total RNA was isolated from tissue culture cells that expressed each of the different BCR/ABL transcripts. Serial log dilutions were prepared, ranging from 100 to 10-5, in RNA isolated from HL60 cells. Laboratories performed 5 independent qRT-PCR reactions for each sample type at each dilution. In addition, 15 qRT-PCR reactions of the 10-3 b3a2 RNA dilution were run to assess reproducibility within and between laboratories. Participants were asked to run the samples following their standard protocols and to report cycle threshold (Ct), quantitative values for BCR/ABL and housekeeping genes, and ratios of BCR/ABL to housekeeping genes for each sample RNA. Thirty-seven (n=37) participants have submitted qRT-PCR results for analysis (36, 37, and 34 labs generated data for b2a2, b3a2, and e1a2, respectively). The limit of detection for this study was defined as the lowest dilution that a Ct value could be detected for all 5 replicates. For b2a2, 15, 16, 4, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. For b3a2, 20, 13, and 4 labs showed a limit of detection at the 10-5, 10-4, and 10-3 dilutions, respectively. For e1a2, 10, 21, 2, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. Log %BCR/ABL ratio values provided a method for comparing results between the different laboratories for each BCR/ABL dilution series. Linear regression analysis revealed concordance among the majority of participant data over the 10-1 to 10-4 dilutions. The overall slope values showed comparable results among the majority of b2a2 (mean=0.939; median=0.9627; range (0.399 - 1.1872)), b3a2 (mean=0.925; median=0.922; range (0.625 - 1.140)), and e1a2 (mean=0.897; median=0.909; range (0.5174 - 1.138)) laboratory results (Fig. 1-3)). Thirty-four (n=34) out of the 37 laboratories reported Ct values for all 15 replicates and only those with a complete data set were included in the inter-lab calculations. Eleven laboratories either did not report their copy number data or used other reporting units such as nanograms or cell numbers; therefore, only 26 laboratories were included in the overall analysis of copy numbers. The median copy number was 348.4, with a range from 15.6 to 547,000 copies (approximately a 4.5 log difference); the median intra-lab %CV was 19.2% with a range from 4.2% to 82.6%. While our international performance evaluation using serially diluted RNA samples has reinforced the fact that heterogeneity exists among clinical laboratories, it has also demonstrated that performance within a laboratory is overall very consistent. Accordingly, the availability of defined BCR/ABL RNAs may facilitate the validation of all phases of quantitative BCR/ABL analysis and may be extremely useful as a tool for monitoring assay performance. Ongoing analyses of these materials, along with the development of additional control materials, may solidify consensus around their application in routine laboratory testing and possible integration in worldwide efforts to standardize quantitative BCR/ABL testing.

Side effects of genome structural changes.

The first extensive catalog of structural human variation was recently released. It showed that large stretches of genomic DNA that vary considerably in copy number were extremely abundant. Thus it is conceivable that they play a major role in functional variation. Consistently, genomic insertions and deletions were shown to contribute to phenotypic differences by modifying not only the expression levels of genes within the aneuploid segments but also of normal copy-number neighboring genes. In this report, we review the possible mechanisms behind this latter effect.

Expression of Staphylococcus aureus clumping factor A in Lactococcus lactis subsp. cremoris using a new shuttle vector.

Staphylococcus aureus harbors redundant adhesins mediating tissue colonization and infection. To evaluate their intrinsic role outside of the staphylococcal background, a system was designed to express them in Lactococcus lactis subsp. cremoris 1363. This bacterium is devoid of virulence factors and has a known genetic background. A new Escherichia coli-L. lactis shuttle and expression vector was constructed for this purpose. First, the high-copy-number lactococcal plasmid pIL253 was equipped with the oriColE1 origin, generating pOri253 that could replicate in E. coli. Second, the lactococcal promoters P23 or P59 were inserted at one end of the pOri253 multicloning site. Gene expression was assessed by a luciferase reporter system. The plasmid carrying P23 (named pOri23) expressed luciferase constitutively at a level 10,000 times greater than did the P59-containing plasmid. Transcription was absent in E. coli. The staphylococcal clumping factor A (clfA) gene was cloned into pOri23 and used as a model system. Lactococci carrying pOri23-clfA produced an unaltered and functional 130-kDa ClfA protein attached to their cell walls. This was indicated both by the presence of the protein in Western blots of solubilized cell walls and by the ability of ClfA-positive lactococci to clump in the presence of plasma. ClfA-positive lactococci had clumping titers (titer of 4,112) similar to those of S. aureus Newman in soluble fibrinogen and bound equally well to solid-phase fibrinogen. These experiments provide a new way to study individual staphylococcal pathogenic factors and might complement both classical knockout mutagenesis and modern in vivo expression technology and signature tag mutagenesis.

Rapid evolution of cancer/testis genes on the X chromosome.

BACKGROUND: Cancer/testis (CT) genes are normally expressed only in germ cells, but can be activated in the cancer state. This unusual property, together with the finding that many CT proteins elicit an antigenic response in cancer patients, has established a role for this class of genes as targets in immunotherapy regimes. Many families of CT genes have been identified in the human genome, but their biological function for the most part remains unclear. While it has been shown that some CT genes are under diversifying selection, this question has not been addressed before for the class as a whole. RESULTS: To shed more light on this interesting group of genes, we exploited the generation of a draft chimpanzee (Pan troglodytes) genomic sequence to examine CT genes in an organism that is closely related to human, and generated a high-quality, manually curated set of human:chimpanzee CT gene alignments. We find that the chimpanzee genome contains homologues to most of the human CT families, and that the genes are located on the same chromosome and at a similar copy number to those in human. Comparison of putative human:chimpanzee orthologues indicates that CT genes located on chromosome X are diverging faster and are undergoing stronger diversifying selection than those on the autosomes or than a set of control genes on either chromosome X or autosomes. CONCLUSION: Given their high level of diversifying selection, we suggest that CT genes are primarily responsible for the observed rapid evolution of protein-coding genes on the X chromosome.

Soft tissue sarcomas with complex genomic profiles.

Soft tissue sarcomas (STS) with complex genomic profiles (50% of all STS) are predominantly composed of spindle cell/pleomorphic sarcomas, including leiomyosarcoma, myxofibrosarcoma, pleomorphic liposarcoma, pleomorphic rhabdomyosarcoma, malignant peripheral nerve sheath tumor, angiosarcoma, extraskeletal osteosarcoma, and spindle cell/pleomorphic unclassified sarcoma (previously called spindle cell/pleomorphic malignant fibrous histiocytoma). These neoplasms show, characteristically, gains and losses of numerous chromosomes or chromosome regions, as well as amplifications. Many of them share recurrent aberrations (e.g., gain of 5p13-p15) that seem to play a significant role in tumor progression and/or metastatic dissemination. In this paper, we review the cytogenetic, molecular genetic, and clinicopathologic characteristics of the most common STS displaying complex genomic profiles. Features of diagnostic or prognostic relevance will be discussed when needed.

Genomic, Epigenomic and Transcriptomic Molecular Characterization of Splenic Marginal Zone Lymphoma Reveals Recurrent Abnormalities in miR-182

Splenic marginal zone lymphoma (SMZL) is a low grade B-cell non-Hodgkin's lymphoma. The molecular pathology of this entity remains poorly understood. To characterise this lymphoma at the molecular level, we performed an integrated analysis of 1) genome wide genetic copy number alterations 2) gene expression profiles and 3) epigenetic DNA methylation profiles.We have previously shown that SMZL is characterised by recurrent alterations of chromosomes 7q, 6q, 3q, 9q and 18; however, gene resolution oligonucleotide array comparative genomic hybridisation did not reveal evidence of cryptic amplification or deletion in these regions. The most frequently lost 7q32 region contains a cluster of miRNAs. qRT-PCR revealed that three of these (miR-182/96/183) show underexpression in SMZL, and miR-182 is somatically mutated in >20% of cases of SMZL, as well as in >20% of cases of follicular lymphoma, and between 5-15% of cases of chronic lymphocytic leukaemia, MALT-lymphoma and hairy cell leukaemia. We conclude that miR-182 is a strong candidate novel tumour suppressor miRNA in lymphoma.The overall gene expression signature of SMZL was found to be strongly distinct fromthose of other lymphomas. Functional analysis of gene expression data revealed SMZL to be characterised by abnormalities in B-cell receptor signalling (especially through the CD19/21-PI3K/AKT pathway) and apoptotic pathways. In addition, genes involved in the response to viral infection appeared upregulated. SMZL shows a unique epigenetic profile, but analysis of differentially methylated genes showed few with methylation related transcriptional deregulation, suggesting that DNA methylation abnormalities are not a critical component of the SMZL malignant phenotype.

Recurrent loss of heterozygosity in 1p36 associated with TNFRSF14 mutations in IRF4 translocation negative pediatric follicular lymphomas.

Pediatric follicular lymphoma is a rare disease that differs genetically and clinically from its adult counterpart. With the exception of pediatric follicular lymphoma with IRF4-translocation, the genetic events associated with these lymphomas have not yet been defined. We applied array-comparative genomic hybridization and molecular inversion probe assay analyses to formalin-fixed paraffin-embedded tissues from 18 patients aged 18 years and under with IRF4 translocation negative follicular lymphoma. All evaluable cases lacked t(14;18). Only 6 of 16 evaluable cases displayed chromosomal imbalances with gains or amplifications of 6pter-p24.3 (including IRF4) and deletion and copy number neutral-loss of heterozygosity in 1p36 (including TNFRSF14) being most frequent. Sequencing of TNFRSF14 located in the minimal region of loss in 1p36.32 showed nine mutations in 7 cases from our series. Two subsets of pediatric follicular lymphoma were delineated according to the presence of molecular alterations, one with genomic aberrations associated with higher grade and/or diffuse large B-cell lymphoma component and more widespread disease, and another one lacking genetic alterations associated with more limited disease.

Deletion, insertion and translocation of DNA sequences contribute to chromosome size polimorphism in Plasmodium berghei

Extensive chromosome size polymorphism in Plasmodium berghei in vivo mitotic multiplication. Size differences between homologous chromosomes mainly involve rearrangements in the subtelomeric regions while internal chromosomal regions are more conserved. Size differences are almost exclusively due to differences in the copy number of a 2.3 kb subtelomeric repeat unit. Not only deletion of 2.3 kb repeats occurs, but addition of new copies of this repeat sometimes results in the formation of enlarged chromosomes. Even chromosomes which originally lack 2.3 kb repeats, can acquire these during mitotic multiplication. In one karyotype mutant, 2.3 kb repeats were inserted within one of the original telomeres of chromosome 4, creating an internal stretch oftelomeric repeats. Chromosome translocation can contribute to chromosome size polymorphism as well We found a karyotype mutant in which chromosome 7 with a size of about 1.4 Mb is translocated to chromosome 13/14 with a size of about 3 Mb, resulting in a rearranged chromosome, which was shown to contain a junction between internal DNA sequences of chromosome 13/14 and subtelomeric 2.3 kb repeats of chromosome 7. In this mutant a new chromosome of 1.4 Mb is present which consists of part of chromosome 13/14.

Subtelomeric structure of Plasmodium falciparum chromosomes

Previous studies of subtelomeric regions in Plasmodium berghei led to the identification of subtelomeric repeats (2.3kb long) present in a variable number at many chromosomal ends. Both loss and increase in 2.3kb-repeat copy number are involved in chromosome-size polymorphisms. Subtelomeric losses leading to chromosome-size polymorphisms have been described by several authors in P.falciparum where the structure of subtelomeric regions is not known in detail. We therefore undertook their characterisation, by means of chromosome walking and jumping techniques, starting from the telomere-flanking sequence present in pPftel.1, the P.falciparum telomeric clone described by Vernick and McCutchan (1988). The results indicate that at least 20 (out of 28) chromosomal ends in P.falciparum 3D7 chromosomes share a subtelomeric region, about 40kb long, covering (but not limited to) the Rep20 region. Non repetitive, AT-rich portions flanking the Rep20 region on both sides are also conserved at most chromosomal ends.

Low-grade and high-grade mammary carcinomas in WAP-T transgenic mice are independent entities distinguished by Met expression.

Mammary carcinomas developing in SV40 transgenic WAP-T mice arise in two distinct histological phenotypes: as differentiated low-grade and undifferentiated high-grade tumors. We integrated different types of information such as histological grading, analysis of aCGH-based gene copy number and gene expression profiling to provide a comprehensive molecular description of mammary tumors in WAP-T mice. Applying a novel procedure for the correlation of gene copy number with gene expression on a global scale, we observed in tumor samples a global coherence between genotype and transcription. This coherence can be interpreted as a matched transcriptional regulation inherited from the cells of tumor origin and determined by the activity of cancer driver genes. Despite common recurrent genomic aberrations, e.g. gain of chr. 15 in most WAP-T tumors, loss of chr. 19 frequently occurs only in low-grade tumors. These tumors show features of "basal-like" epithelial differentiation, particularly expression of keratin 14. The high-grade tumors are clearly separated from the low-grade tumors by strong expression of the Met gene and by coexpression of epithelial (e.g. keratin 18) and mesenchymal (e.g. vimentin) markers. In high-grade tumors, the expression of the nonmutated Met protein is associated with Met-locus amplification and Met activity. The role of Met as a cancer driver gene is supported by the contribution of active Met signaling to motility and growth of mammary tumor-derived cells. Finally, we discuss the independent origin of low- and high-grade tumors from distinct cells of tumor origin, possibly luminal progenitors, distinguished by Met gene expression and Met signaling.

Variation between geographical populations of Lutzomyia (Nyssomyia) whitmani (Antunes & Coutinho, 1939) sensu lato (Diptera:Psychodidae:Phlebotominae) in Brazil

Phylogenetic analysis of morphometric and biological characters indicated that there are two distinct forms of Lutzomyia whitmani in Brazil: one is present both north and south of the River Amazonas in the State of Pará while the other occurs in northeast Brazil, in the State of Ceará, and further south, including the type locality in State of Bahia. The Amazonian form is reportedly neither strongly anthropophilic nor synanthropic, and it is the vector of Leishmania shawi; whereas the southern form is often collected peridomestically, while biting man, and has been found infected with Le.(V.) braziliensis. The ratio of the length of the genital filaments to that the genital pump was found to be consistently smaller in males of the Amazonian populations. A middle repetitive DNA element was isolated by differentially screening a genomic library made using Amazonian material, and the sequence was diagnostic for this form of Lu. whitmani (being absent or occurring in low copy number in the southern form). The total evidence suggests there are at least two, geographically-isolated forms of Lu. whitmani, which may represent different cryptic species.

Thirty new loci for age at menarche identified by a meta-analysis of genome-wide association studies.

To identify loci for age at menarche, we performed a meta-analysis of 32 genome-wide association studies in 87,802 women of European descent, with replication in up to 14,731 women. In addition to the known loci at LIN28B (P = 5.4 × 10⁻⁶⁰) and 9q31.2 (P = 2.2 × 10⁻³³), we identified 30 new menarche loci (all P < 5 × 10⁻⁸) and found suggestive evidence for a further 10 loci (P < 1.9 × 10⁻⁶). The new loci included four previously associated with body mass index (in or near FTO, SEC16B, TRA2B and TMEM18), three in or near other genes implicated in energy homeostasis (BSX, CRTC1 and MCHR2) and three in or near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and gene-set enrichment pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing.

arrayMap 2014: an updated cancer genome resource.

Somatic copy number aberrations (CNA) represent a mutation type encountered in the majority of cancer genomes. Here, we present the 2014 edition of arrayMap (http://www.arraymap.org), a publicly accessible collection of pre-processed oncogenomic array data sets and CNA profiles, representing a vast range of human malignancies. Since the initial release, we have enhanced this resource both in content and especially with regard to data mining support. The 2014 release of arrayMap contains more than 64,000 genomic array data sets, representing about 250 tumor diagnoses. Data sets included in arrayMap have been assembled from public repositories as well as additional resources, and integrated by applying custom processing pipelines. Online tools have been upgraded for a more flexible array data visualization, including options for processing user provided, non-public data sets. Data integration has been improved by mapping to multiple editions of the human reference genome, with the majority of the data now being available for the UCSC hg18 as well as GRCh37 versions. The large amount of tumor CNA data in arrayMap can be freely downloaded by users to promote data mining projects, and to explore special events such as chromothripsis-like genome patterns.

Sample preparation for two-dimensional gel electrophoresis.

The choice of sample preparation protocol is a critical influential factor for isoelectric focusing which in turn affects the two-dimensional gel result in terms of quality and protein species distribution. The optimal protocol varies depending on the nature of the sample for analysis and the properties of the constituent protein species (hydrophobicity, tendency to form aggregates, copy number) intended for resolution. This review explains the standard sample buffer constituents and illustrates a series of protocols for processing diverse samples for two-dimensional gel electrophoresis, including hydrophobic membrane proteins. Current methods for concentrating lower abundance proteins, by removal of high abundance proteins, are also outlined. Finally, since protein staining is becoming increasingly incorporated into the sample preparation procedure, we describe the principles and applications of current (and future) pre-electrophoretic labelling methods.