Systematic interpretation of molecular beacon polymerase chain reaction for identifying rpoB mutations in Mycobacterium tuberculosis isolates with mixed resistant and susceptible bacteria

Detection of multidrug-resistant tuberculosis (MDR-TB), a frequent cause of treatment failure, takes 2 or more weeks to identify by culture. RIF-resistance is a hallmark of MDR-TB, and detection of mutations in the rpoB gene of Mycobacterium tuberculosis using molecular beacon probes with real-time quantitative polymerase chain reaction (qPCR) is a novel approach that takes ≤2 days. However, qPCR identification of resistant isolates, particularly for isolates with mixed RIF-susceptible and RIF-resistant bacteria, is reader dependent and limits its clinical use. The aim of this study was to develop an objective, reader-independent method to define rpoB mutants using beacon qPCR. This would facilitate the transition from a research protocol to the clinical setting, where high-throughput methods with objective interpretation are required. For this, DNAs from 107 M. tuberculosis clinical isolates with known susceptibility to RIF by culture-based methods were obtained from 2 regions where isolates have not previously been subjected to evaluation using molecular beacon qPCR: the Texas–Mexico border and Colombia. Using coded DNA specimens, mutations within an 81-bp hot spot region of rpoB were established by qPCR with 5 beacons spanning this region. Visual and mathematical approaches were used to establish whether the qPCR cycle threshold of the experimental isolate was significantly higher (mutant) compared to a reference wild-type isolate. Visual classification of the beacon qPCR required reader training for strains with a mixture of RIF-susceptible and RIF-resistant bacteria. Only then had the visual interpretation by an experienced reader had 100% sensitivity and 94.6% specificity versus RIF-resistance by culture phenotype and 98.1% sensitivity and 100% specificity versus mutations based on DNA sequence. The mathematical approach was 98% sensitive and 94.5% specific versus culture and 96.2% sensitive and 100% specific versus DNA sequence. Our findings indicate the mathematical approach has advantages over the visual reading, in that it uses a Microsoft Excel template to eliminate reader bias or inexperience, and allows objective interpretation from high-throughput analyses even in the presence of a mixture of RIF-resistant and RIF-susceptible isolates without the need for reader training.^

Exploration of genetic susceptibility to pancreatic cancer: A GWAS approach

Pancreatic cancer is the 4th most common cause for cancer death in the United States, accompanied by less than 5% five-year survival rate based on current treatments, particularly because it is usually detected at a late stage. Identifying a high-risk population to launch an effective preventive strategy and intervention to control this highly lethal disease is desperately needed. The genetic etiology of pancreatic cancer has not been well profiled. We hypothesized that unidentified genetic variants by previous genome-wide association study (GWAS) for pancreatic cancer, due to stringent statistical threshold or missing interaction analysis, may be unveiled using alternative approaches. To achieve this aim, we explored genetic susceptibility to pancreatic cancer in terms of marginal associations of pathway and genes, as well as their interactions with risk factors. We conducted pathway- and gene-based analysis using GWAS data from 3141 pancreatic cancer patients and 3367 controls with European ancestry. Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Using the logistic kernel machine (LKM) test, we analyzed 17906 genes defined by University of California Santa Cruz (UCSC) database. Using the likelihood ratio test (LRT) in a logistic regression model, we analyzed 177 pathways and 17906 genes for interactions with risk factors in 2028 pancreatic cancer patients and 2109 controls with European ancestry. After adjusting for multiple comparisons, six pathways were marginally associated with risk of pancreatic cancer ( P < 0.00025): Fc epsilon RI signaling, maturity onset diabetes of the young, neuroactive ligand-receptor interaction, long-term depression (Ps < 0.0002), and the olfactory transduction and vascular smooth muscle contraction pathways (P = 0.0002; Nine genes were marginally associated with pancreatic cancer risk (P < 2.62 × 10−5), including five reported genes (ABO, HNF1A, CLPTM1L, SHH and MYC), as well as four novel genes (OR13C4, OR 13C3, KCNA6 and HNF4 G); three pathways significantly interacted with risk factors on modifying the risk of pancreatic cancer (P < 2.82 × 10−4): chemokine signaling pathway with obesity ( P < 1.43 × 10−4), calcium signaling pathway (P < 2.27 × 10−4) and MAPK signaling pathway with diabetes (P < 2.77 × 10−4). However, none of the 17906 genes tested for interactions survived the multiple comparisons corrections. In summary, our current GWAS study unveiled unidentified genetic susceptibility to pancreatic cancer using alternative methods. These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer, once confirmed, will shed promising light on the prevention and treatment of this disease. ^

Pathway Semantics: An Algebraic Data Driven Algorithm to Generate Hypotheses about Molecular Patterns Underlying Disease Progression

The overarching goal of the Pathway Semantics Algorithm (PSA) is to improve the in silico identification of clinically useful hypotheses about molecular patterns in disease progression. By framing biomedical questions within a variety of matrix representations, PSA has the flexibility to analyze combined quantitative and qualitative data over a wide range of stratifications. The resulting hypothetical answers can then move to in vitro and in vivo verification, research assay optimization, clinical validation, and commercialization. Herein PSA is shown to generate novel hypotheses about the significant biological pathways in two disease domains: shock / trauma and hemophilia A, and validated experimentally in the latter. The PSA matrix algebra approach identified differential molecular patterns in biological networks over time and outcome that would not be easily found through direct assays, literature or database searches. In this dissertation, Chapter 1 provides a broad overview of the background and motivation for the study, followed by Chapter 2 with a literature review of relevant computational methods. Chapters 3 and 4 describe PSA for node and edge analysis respectively, and apply the method to disease progression in shock / trauma. Chapter 5 demonstrates the application of PSA to hemophilia A and the validation with experimental results. The work is summarized in Chapter 6, followed by extensive references and an Appendix with additional material.

A Genomic Approach to Identify The Notch Pathway as a Putative Tumor Suppressor in Endometrial Cancer

Endometrial cancer is the most common gynecological malignancy and the fourth most frequently diagnosed cancer among women. The molecular changes that distinguish normal endometrium from endometrial carcinoma are not thoroughly understood. Identification of these changes could potentially aid in identifying at-risk women who are especially prone to develop endometrial cancer, such as obese women and women with Lynch Syndrome. A microarray analysis was performed using normal endometrium from thin and obese women and cancerous endometrium from obese women. We validated the differential expression of ten genes whose expression was significantly up-regulated or down-regulated using qRT-PCR. All of the genes had distinct expression levels depending on the endometrial carcinoma histotype. As a result, they could serve as molecular markers to distinguish between normal endometrium and endometrial cancer, as well as between low grade endometrial carcinomas and high grade endometrial carcinomas. Two of the ten genes validated, HEYL and HES1, are down-stream targets of the Notch signaling pathway. HEYL and HES1 were identified by microarray and qRT-PCR to have a significant decrease in expression in endometrial carcinomas compared to normal endometrium. We further analyzed the differential expression of other components of the Notch signaling pathway, Notch4 and Jagged1. They were also identified by qRT-PCR to be significantly down-regulated in endometrial carcinomas compared to normal endometrium. Therefore, we believe the Notch signaling pathway to act as a tumor suppressor in endometrial carcinomas.

Molecular studies of human high grade B cell non -Hodgkin's lymphomas

The non-Hodgkin's B cell lymphomas are a diverse group of neoplastic diseases. The incidence rate of the malignant tumors has been rising rapidly over the past twenty years in the United States and worldwide. The lack of insight to pathogenesis of the disease poses a significant problem in the early detection and effective treatment of the human malignancies. These studies attempted to investigate the molecular basis of pathogenesis of the human high grade B cell non-Hodgkin's lymphomas with a reverse genetic approach. The specific objective was to clone gene(s) which may play roles in development and progression of human high grade B cell non-Hodgkin's lymphomas.^ The messenger RNAs from two high grade B cell lymphoma lines, CJ and RR, were used for construction of cDNA libraries. Differential screening of the derived cDNA libraries yielded a 1.4 kb cDNA clone. The gene, designated as NHL-B1.4, was shown to be highly amplified and over-expressed in the high grade B cell lymphoma lines. It was not expressed in the peripheral blood lymphoid cells from normal donors. However, it was inducible in peripheral blood T lymphocytes by a T cell mitogen, PHA, but could not be activated in normal B cells by B cell mitogen PMA. Further molecular characterization revealed that the gene may have been rearranged in the RR and some other B cell lymphoma lines. The coding capacity of the cDNA has been confirmed by a rabbit reticulocyte lysate and wheat germ protein synthesis system. A recombinant protein with a molecular weight of approximate 30 kDa was visualized in autoradiogram. Polyclonal antisera have been generated by immunization of two rabbits with the NHL-B1.4 recombinant protein produced in the E. coli JM109. The derived antibody can recognize a natural protein with molecular weight of 49 kDa in cell lysate of activated peripheral T lymphocytes of normal donors and both the cell lysate and supernatant of RR B cell lymphoma lines. The possible biologic functions of the molecule has been tested preliminarily in a B lymphocyte proliferation assay. It was found that the Q-sepharose chromatograph purified supernatant of COS cell transfection could increase tritiated thymidine uptake by B lymphocytes but not by T lymphocytes. The B cell stimulatory activity of the supernatant of COS cell tranfection could be neutralized by the polyclonal antisera, indicating that the NHL-B1.4 gene product may be a molecule with BCGF-like activity.^ The expression profiles of NHL-B1.4 in normal and neoplastic lymphoid cells were consistent with the current B lymphocyte activation model and autocrine hypothesis of high grade B cell lymphomagenesis. These results suggested that the NHL-B1.4 cDNA may be a disease-related gene of human high grade B cell lymphomas, which may codes for a postulated B cell autocrine growth factor. ^

Phosphorylation levels of JAK2 after addition of increasing Epo-concentrations: model approach (Figure 5a)

The present dataset contain source data for Figure 5a from Schilling et al., 2009. Cell fate decisions are regulated by the coordinated activation of signalling pathways such as the extracellular signal-regulated kinase (ERK) cascade, but contributions of individual kinase isoforms are mostly unknown. The authors combined quantitative data from erythropoietin-induced pathway activation in primary erythroid progenitor (colony-forming unit erythroid stage, CFU-E) cells with mathematical modelling, in order to predict and experimentally confirmed a distributive ERK phosphorylation mechanism in CFU-E cells. The authors found evidences that double-phosphorylated ERK1 attenuates proliferation beyond a certain activation level, whereas activated ERK2 enhances proliferation with saturation kinetics. Phosphorylation levels of JAK2 at 7 min after stimulation for Epo concentrations ranging from 0.1 to 1000 U/ml were simulated.

Liver proteome profiling of juvenile Chinese sturgeon (Acipenser sinensis) using GeLC-MS/MS approach

Chinese sturgeon (Acipenser sinensis), mainly distributed in the Yangtze River, has been listed as a grade I protected animal in China because of a dramatic decline in population owing to loss of natural habitat for reproduction and interference by human activities. Understanding the proteome profile of Chinese sturgeon liver would provide an invaluable resource for protecting and increasing the stocks of this species. In this study, we have analyzed proteome profiles of juvenile Chinese sturgeon liver using a one-dimensional gel electrophoresis coupled to LC-MS/MS approach. A total of 1059 proteins and 2084 peptides were identified. The liver proteome was found to be associated with diverse biological processes, cellular components and molecular functions. The proteome profile identified a variety of significant pathways including carbohydrate metabolism, fatty acid metabolism and amino acid metabolism pathways. It also established a network for protein biosynthesis, folding and catabolic processes. The proteome profile established in this study can be used for understanding the development of Chinese sturgeon and studying the molecular mechanisms of action under environmental or chemical stress, providing very useful omics information that can be applied to preserve this species.

Response of three krill species to hypoxia and warming: An experimental approach to oxygen minimum zones expansion in coastal ecosystems

To understand the adaptation of euphausiid (krill) species to oxygen minimum zones (OMZ), respiratory response and stress experiments combining hypoxia/reoxygenation exposure with warming were conducted. Experimental krill species were obtained from the Antarctic (South Georgia area), the Humboldt Current system (HCS, Chilean coast), and the Northern California Current system (NCCS, Oregon). Euphausia mucronata from the HCS shows oxyconforming or oxygen partial pressure (pO2)-dependent respiration below 80% air saturation (18 kPa). Normoxic subsurface oxygenation in winter posed a "high oxygen stress" for this species. The NCCS krill, Euphausia pacifica, and the Antarctic krill, Euphausia superba maintain respiration rates constant down to low critical pO2 values of 6 kPa (30% air saturation) and 11 kPa (55% air saturation), respectively. Antarctic krill had the lowest antioxidant enzyme activities, but the highest concentrations of the molecular antioxidant glutathione (GSH) and was not affected by 6 h exposure to moderate hypoxia. Temperate krill species had higher SOD (superoxide dismutase) values in winter than in summer, which relate to higher winter metabolic rate (E. pacifica). In all species, antioxidant enzyme activities remained constant during hypoxic exposure at habitat temperature. Warming by 7°C above habitat temperature in summer increased SOD activities and GSH levels in E. mucronata (HCS), but no oxidative damage occurred. In winter, when the NCCS is well mixed and the OMZ is deeper, +4°C of warming combined with hypoxia represents a lethal condition for E. pacifica. In summer, when the OMZ expands upwards (100 m subsurface), antioxidant defences counteracted hypoxia and reoxygenation effects in E. pacifica, but only at mildly elevated temperature (+2°C). In this season, experimental warming by +4°C reduced antioxidant activities and the hypoxia combination again caused mortality of exposed specimens. We conclude that a climate change scenario combining warming and hypoxia represents a serious threat to E. pacifica and, as a consequence, NCCS food webs.

Study on the effects of near-future ocean acidification on marine yeasts: a microcosm approach

Marine yeasts play an important role in biodegradation and nutrient cycling and are often associated with marine flora and fauna. They show maximum growth at pH levels lower than present-day seawater pH. Thus, contrary to many other marine organisms, they may actually profit from ocean acidification. Hence, we conducted a microcosm study, incubating natural seawater from the North Sea at present-day pH (8.10) and two near-future pH levels (7.81 and 7.67). Yeasts were isolated from the initial seawater sample and after 2 and 4 weeks of incubation. Isolates were classified by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) and representative isolates were identified by partial sequencing of the large subunit rRNA gene. From the initial seawater sample, we predominantly isolated a yeast-like filamentous fungus related to Aureobasidium pullulans, Cryptococcus sp., Candida sake, and various cold-adapted yeasts. After incubation, we found more different yeast species at near-future pH levels than at present-day pH. Yeasts reacting to low pH were related to Leucosporidium scottii, Rhodotorula mucilaginosa, Cryptococcus sp., and Debaryomyces hansenii. Our results suggest that these yeasts will benefit from seawater pH reductions and give a first indication that the importance of yeasts will increase in a more acidic ocean.

Stopping power dependence of nitrogen sputtering yields in copper nitride films under swift-ion irradiation: Exciton model approach

Nitrogen sputtering yields as high as 104 atoms/ion, are obtained by irradiating N-rich-Cu3N films (N concentration: 33 ± 2 at.%) with Cu ions at energies in the range 10?42 MeV. The kinetics of N sputtering as a function of ion fluence is determined at several energies (stopping powers) for films deposited on both, glass and silicon substrates. The kinetic curves show that the amount of nitrogen release strongly increases with rising irradiation fluence up to reaching a saturation level at a low remaining nitrogen fraction (5?10%), in which no further nitrogen reduction is observed. The sputtering rate for nitrogen depletion is found to be independent of the substrate and to linearly increase with electronic stopping power (Se). A stopping power (Sth) threshold of ?3.5 keV/nm for nitrogen depletion has been estimated from extrapolation of the data. Experimental kinetic data have been analyzed within a bulk molecular recombination model. The microscopic mechanisms of the nitrogen depletion process are discussed in terms of a non-radiative exciton decay model. In particular, the estimated threshold is related to a minimum exciton density which is required to achieve efficient sputtering rates.

Quantum molecular dynamics and ab initio studies of the crystal structure of hydrogen and deuterium

Hydrogen isotopes play a critical role both in inertial and magnetic confinemen Nuclear Fusion. Since the preferent fuel needed for this technology is a mixture of deuterium and tritium. The study of these isotopes particularly at very low temperatures carries a technological interest in other applications. The present line promotes a deep study on the structural configuration that hydrogen and deuterium adopt at cryogenic temperatures and at high pressures. Typical conditions occurring in present Inertial Fusion target designs. Our approach is aims to determine the crystal structure characteristics, phase transitions and other parameters strongly correlated to variations of temperature and pressure.

Ab initio Methodology for Calculating Elastic Constant and Sound Velocity for Hydrogen and Deuterium in Molecular Solid Phase

Hydrogen isotopes play a critical role both in inertial and magnetic confinement Nuclear Fusion. Since the preferent fuel needed for this technology is a mixture of deuterium and tritium. The study of these isotopes particularly at very low temperatures carries a technological interest in other applications. The present line promotes a deep study on the structural configuration that hydrogen and deuterium adopt at cryogenic temperatures and at high pressures. Typical conditions occurring in present Inertial Fusion target designs. Our approach is aims to determine the crystal structure characteristics, phase transitions and other parameters strongly correlated to variations of temperature and pressure. With this results is possible calculated the elastic constant and sound velocity for hydrogen and deuterium in molecular solid phase.

Sensitization profiles and cross-reactivity among plant allergens. Molecular mechanisms of food allergy development and the role of enhancers

La prevalencia de las alergias está aumentando desde mediados del siglo XX, y se estima que actualmente afectan a alrededor del 2-8 % de la población, pero las causas de este aumento aún no están claras. Encontrar el origen del mecanismo por el cual una proteína inofensiva se convierte en capaz de inducir una respuesta alérgica es de vital importancia para prevenir y tratar estas enfermedades. Aunque la caracterización de alérgenos relevantes ha ayudado a mejorar el manejo clínico y a aclarar los mecanismos básicos de las reacciones alérgicas, todavía queda un largo camino para establecer el origen de la alergenicidad y reactividad cruzada. El objetivo de esta tesis ha sido caracterizar las bases moleculares de la alergenicidad tomando como modelo dos familias de panalergenos (proteínas de transferencia de lípidos –LTPs- y taumatinas –TLPs-) y estudiando los mecanismos que median la sensibilización y la reactividad cruzada para mejorar tanto el diagnóstico como el tratamiento de la alergia. Para ello, se llevaron a cabo dos estrategias: estudiar la reactividad cruzada de miembros de familias de panalérgenos; y estudiar moléculas-co-adyuvantes que pudieran favorecer la capacidad alergénica de dichas proteínas. Para estudiar la reactividad cruzada entre miembros de la misma familia de proteínas, se seleccionaron LTPs y TLPs, descritas como alergenos, tomando como modelo la alergia a frutas. Por otra parte, se estudiaron los perfiles de sensibilización a alérgenos de trigo relacionados con el asma del panadero, la enfermedad ocupacional más relevante de origen alérgico. Estos estudios se llevaron a cabo estandarizando ensayos tipo microarrays con alérgenos y analizando los resultados por la teoría de grafos. En relación al estudiar moléculas-co-adyuvantes que pudieran favorecer la capacidad alergénica de dichas proteínas, se llevaron a cabo estudios sobre la interacción de los alérgenos alimentarios con células del sistema inmune humano y murino y el epitelio de las mucosas, analizando la importancia de moléculas co-transportadas con los alérgenos en el desarrollo de una respuesta Th2. Para ello, Pru p 3(LTP y alérgeno principal del melocotón) se selección como modelo para llevarlo a cabo. Por otra parte, se analizó el papel de moléculas activadoras del sistema inmune producidas por patógenos en la inducción de alergias alimentarias seleccionando el modelo kiwi-alternaria, y el papel de Alt a 1, alérgeno mayor de dicho hongo, en la sensibilización a Act d 2, alérgeno mayor de kiwi. En resumen, el presente trabajo presenta una investigación innovadora aportando resultados de gran utilidad tanto para la mejora del diagnóstico como para nuevas investigaciones sobre la alergia y el esclarecimiento final de los mecanismos que caracterizan esta enfermedad. ABSTRACT Allergies are increasing their prevalence from mid twentieth century, and they are currently estimated to affect around 2-8% of the population but the underlying causes of this increase remain still elusive. The understanding of the mechanism by which a harmless protein becomes capable of inducing an allergic response provides us the basis to prevent and treat these diseases. Although the characterization of relevant allergens has led to improved clinical management and has helped to clarify the basic mechanisms of allergic reactions, it seems justified in aspiring to molecularly dissecting these allergens to establish the structural basis of their allergenicity and cross-reactivity. The aim of this thesis was to characterize the molecular basis of the allergenicity of model proteins belonging to different families (Lipid Transfer Proteins –LTPs-, and Thaumatin-like Proteins –TLPs-) in order to identify mechanisms that mediate sensitization and cross reactivity for developing new strategies in the management of allergy, both diagnosis and treatment, in the near future. With this purpose, two strategies have been conducted: studies of cross-reactivity among panallergen families and molecular studies of the contribution of cofactors in the induction of the allergic response by these panallergens. Following the first strategy, we studied the cross-reactivity among members of two plant panallergens (LTPs , Lipid Transfer Proteins , and TLPs , Thaumatin-like Proteins) using the peach allergy as a model. Similarly, we characterized the sensitization profiles to wheat allergens in baker's asthma development, the most relevant occupational disease. These studies were performed using allergen microarrays and the graph theory for analyzing the results. Regarding the second approach, we analyzed the interaction of plant allergens with immune and epithelial cells. To perform these studies , we examined the importance of ligands and co-transported molecules of plant allergens in the development of Th2 responses. To this end, Pru p 3, nsLTP (non-specific Lipid Transfer Protein) and peach major allergen, was selected as a model to investigate its interaction with cells of the human and murine immune systems as well as with the intestinal epithelium and the contribution of its ligand in inducing an allergic response was studied. Moreover, we analyzed the role of pathogen associated molecules in the induction of food allergy. For that, we selected the kiwi- alternaria system as a model and the role of Alt a 1 , major allergen of the fungus, in the development of Act d 2-sensitization was studied. In summary, this work presents an innovative research providing useful results for improving diagnosis and leading to further research on allergy and the final clarification of the mechanisms that characterize this disease.

Quantitative modeling of stochastic systems in molecular biology by using stochastic Petri nets

An integrated understanding of molecular and developmental biology must consider the large number of molecular species involved and the low concentrations of many species in vivo. Quantitative stochastic models of molecular interaction networks can be expressed as stochastic Petri nets (SPNs), a mathematical formalism developed in computer science. Existing software can be used to define molecular interaction networks as SPNs and solve such models for the probability distributions of molecular species. This approach allows biologists to focus on the content of models and their interpretation, rather than their implementation. The standardized format of SPNs also facilitates the replication, extension, and transfer of models between researchers. A simple chemical system is presented to demonstrate the link between stochastic models of molecular interactions and SPNs. The approach is illustrated with examples of models of genetic and biochemical phenomena where the UltraSAN package is used to present results from numerical analysis and the outcome of simulations.

Detection of 400-year-old Yersinia pestis DNA in human dental pulp: An approach to the diagnosis of ancient septicemia

Ancient septicemic plague epidemics were reported to have killed millions of people for 2 millenniums. However, confident diagnosis of ancient septicemia solely on the basis of historical clinical observations is not possible. The lack of suitable infected material has prevented direct demonstration of ancient septicemia; thus, the history of most infections such as plague remains hypothetical. The durability of dental pulp, together with its natural sterility, makes it a suitable material on which to base such research. We hypothesized that it would be a lasting refuge for Yersinia pestis, the plague agent. DNA extracts were made from the dental pulp of 12 unerupted teeth extracted from skeletons excavated from 16th and 18th century French graves of persons thought to have died of plague (“plague teeth”) and from 7 ancient negative control teeth. PCRs incorporating ancient DNA extracts and primers specific for the human β-globin gene demonstrated the absence of inhibitors in these preparations. The incorporation of primers specific for Y. pestis rpoB (the RNA polymerase β-subunit-encoding gene) and the recognized virulence-associated pla (the plasminogen activator-encoding gene) repeatedly yielded products that had a nucleotide sequence indistinguishable from that of modern day isolates of the bacterium. The specific pla sequence was obtained from 6 of 12 plague skeleton teeth but 0 of 7 negative controls (P < 0.034, Fisher exact test). A nucleic acid-based confirmation of ancient plague was achieved for historically identified victims, and we have confirmed the presence of the disease at the end of 16th century in France. Dental pulp is an attractive target in the quest to determine the etiology of septicemic illnesses detected in ancient corpses. Molecular techniques could be applied to this material to resolve historical outbreaks.