A multiplex bead-based suspension array assay for interrogation of phylogenetically informative single nucleotide polymorphisms for Bacillus anthracis

Single nucleotide polymorphisms (SNPs) are abundant in genomes of all species and represent informative DNA markers extensively used to analyze phylogenetic relationships between strains. Medium to high throughput, open methodologies able to test many SNPs in a minimum time are therefore in great need. By using the versatile Luminex® xTAG technology, we developed an efficient multiplexed SNP genotyping assay to score 13 phylogenetically informative SNPs within the genome of Bacillus anthracis. The Multiplex Oligonucleotide Ligation-PCR procedure (MOL-PCR) described by Deshpande et al., 2010 has been modified and adapted for simultaneous interrogation of 13 biallelic canonical SNPs in a 13-plex assay. Changes made to the originally published method include the design of allele-specific dual-priming-oligonucleotides (DPOs) as competing detection probes (MOLigo probes) and use of asymmetric PCR reaction for signal amplification and labeling of ligation products carrying SNP targets. These innovations significantly reduce cross-reactivity observed when initial MOLigo probes were used and enhance hybridization efficiency onto the microsphere array, respectively. When evaluated on 73 representative samples, the 13-plex assay yielded unambiguous SNP calls and lineage affiliation. Assay limit of detection was determined to be 2 ng of genomic DNA. The reproducibility, robustness and easy-of-use of the present method were validated by a small-scale proficiency testing performed between four European laboratories. While cost-effective compared to other singleplex methods, the present MOL-PCR method offers a high degree of flexibility and scalability. It can easily accommodate newly identified SNPs to increase resolving power to the canSNP typing of B. anthracis.     

A Real-Time PCR Array for Hierarchical Identification of Francisella Isolates

A robust, rapid and flexible real-time PCR assay for hierarchical genetic typing of clinical and environmental isolates of Francisella is presented. Typing markers were found by multiple genome and gene comparisons, from which 23 canonical single nucleotide polymorphisms (canSNPs) and 11 canonical insertion-deletion mutations (canINDELs) were selected to provide phylogenetic guidelines for classification from genus to isolate level. The specificity of the developed assay, which uses 68 wells of a 96-well real-time PCR format with a detection limit of 100 pg DNA, was assessed using 62 Francisella isolates of diverse genetic and geographical origins. It was then successfully used for typing 14 F. tularensis subsp. holarctica isolates obtained from tularemia patients in Sweden in 2008 and five more genetically diverse Francisella isolates of global origins. When applied to human ulcer specimens for direct pathogen detection the results were incomplete due to scarcity of DNA, but sufficient markers were identified to detect fine-resolution differences among F. tularensis subsp. holarctica isolates causing infection in the patients. In contrast to other real-time PCR assays for Francisella, which are typically designed for specific detection of a species, subspecies, or strain, this type of assay can be easily tailored to provide appropriate phylogenetic and/or geographical resolution to meet the objectives of the analysis.     

Microevolution of the Chromosomal Region of Acute Disease Antigen A (adaA) in the Query (Q) Fever Agent Coxiella burnetii

The acute disease antigen A (adaA) gene is believed to be associated with Coxiella burnetii strains causing acute Q fever. The detailed analysis of the adaA genomic region of 23 human- and 86 animal-derived C. burnetii isolates presented in this study reveals a much more polymorphic appearance and distribution of the adaA gene, resulting in a classification of C. burnetii strains of better differentiation than previously anticipated. Three different genomic variants of the adaA gene were identified which could be detected in isolates from acute and chronic patients, rendering the association of adaA positive strains with acute Q fever disease disputable. In addition, all adaA positive strains in humans and animals showed the occurrence of the QpH1 plasmid. All adaA positive isolates of acute human patients except one showed a distinct SNP variation at position 431, also predominant in sheep strains, which correlates well with the observation that sheep are a major source of human infection. Furthermore, the phylogenetic analysis of the adaA gene revealed three deletion events and supported the hypothesis that strain Dugway 5J108-111 might be the ancestor of all known C. burnetii strains. Based on our findings, we could confirm the QpDV group and we were able to define a new genotypic cluster. The adaA gene polymorphisms shown here improve molecular typing of Q fever, and give new insights into microevolutionary adaption processes in C. burnetii.     

Single-Nucleotide Polymorphism Phylotyping of Escherichia coli

We describe a rapid and easily automated phylogenetic grouping technique based on analysis of bacterial genome single-nucleotide polymorphisms (SNPs). We selected 13 SNPs derived from a complete sequence analysis of 11 essential genes previously used for multilocus sequence typing (MLST) of 30 Escherichia coli strains representing the genetic diversity of the species. The 13 SNPs were localized in five genes, trpA, trpB, putP, icdA, and polB, and were selected to allow recovery of the main phylogenetic groups (groups A, B1, E, D, and B2) and subgroups of the species. In the first step, we validated the SNP approach in silico by extracting SNP data from the complete sequences of the five genes for a panel of 65 pathogenic strains belonging to different E. coli pathovars, which were previously analyzed by MLST. In the second step, we determined these SNPs by dideoxy single-base extension of unlabeled oligonucleotide primers for a collection of 183 commensal and extraintestinal clinical E. coli isolates and compared the SNP phylotyping method to previous well-established typing methods. This SNP phylotyping method proved to be consistent with the other methods for assigning phylogenetic groups to the different E. coli strains. In contrast to the other typing methods, such as multilocus enzyme electrophoresis, ribotyping, or PCR phylotyping using the presence/absence of three genomic DNA fragments, the SNP typing method described here is derived from a solid phylogenetic analysis, and the results obtained by this method are more meaningful. Our results indicate that similar approaches may be used for a wide variety of bacterial species.     

In Vitro and In Vivo Infectious Potential of Coxiella burnetii: A Study on Belgian Livestock Isolates

Q-fever is a zoonosis caused by the gram-negative obligate intracellular pathogen Coxiella burnetii. Since its discovery, and particularly following the recent outbreaks in the Netherlands, C. burnetii appeared as a clear public health concern. In the present study, the infectious potential displayed by goat and cattle isolates of C. burnetii was compared to a reference strain (Nine Mile) using both in vitro (human HeLa and bovine macrophage cells) and in vivo (BALB/c mice) models. The isolates had distant genomic profiles with one - the goat isolate - being identical to the predominant strain circulating in the Netherlands during the 2007–2010 outbreaks. Infective doses were established with ethidium monoazide-PCR for the first time here applied to C. burnetii. This method allowed for the preparation of reproducible and characterized inocula thanks to its capacity to discriminate between live and dead cells. Globally, the proliferative capacity of the Nine Mile strain in cell lines and mice was higher compared to the newly isolated field strains. In vitro, the bovine C. burnetii isolate multiplied faster in a bovine macrophage cell line, an observation tentatively explained by the preferential specificity of this strain for allogeneic host cells. In the BALB/c mouse model, however, the goat and bovine isolates multiplied at about the same rate indicating no peculiar hypervirulent behavior in this animal model.     

Insights to Genetic Characterization Tools for Epidemiological Tracking of Francisella tularensis in Sweden

Tularaemia, caused by the bacterium Francisella tularensis, is endemic in Sweden and is poorly understood. The aim of this study was to evaluate the effectiveness of three different genetic typing systems to link a genetic type to the source and place of tularemia infection in Sweden. Canonical single nucleotide polymorphisms (canSNPs), MLVA including five variable number of tandem repeat loci and PmeI-PFGE were tested on 127 F. tularensis positive specimens collected from Swedish case-patients. All three typing methods identified two major genetic groups with near-perfect agreement. Higher genetic resolution was obtained with canSNP and MLVA compared to PFGE; F. tularensis samples were first assigned into ten phylogroups based on canSNPs followed by 33 unique MLVA types. Phylogroups were geographically analysed to reveal complex phylogeographic patterns in Sweden. The extensive phylogenetic diversity found within individual counties posed a challenge to linking specific genetic types with specific geographic locations. Despite this, a single phylogroup (B.22), defined by a SNP marker specific to a lone Swedish sequenced strain, did link genetic type with a likely geographic place. This result suggests that SNP markers, highly specific to a particular reference genome, may be found most frequently among samples recovered from the same location where the reference genome originated. This insight compels us to consider whole-genome sequencing (WGS) as the appropriate tool for effectively linking specific genetic type to geography. Comparing the WGS of an unknown sample to WGS databases of archived Swedish strains maximizes the likelihood of revealing those rare geographically informative SNPs.     

Effect of a glucose impulse on the CcpA regulon in Staphylococcus aureus

Background

A low genetic diversity in Francisella tularensis has been documented. Current DNA based genotyping methods for typing F. tularensis offer a limited and varying degree of subspecies, clade and strain level discrimination power. Whole genome sequencing is the most accurate and reliable method to identify, type and determine phylogenetic relationships among strains of a species. However, lower cost typing schemes are necessary in order to enable typing of hundreds or even thousands of isolates.

Results

We have generated a high-resolution phylogenetic tree from 40 Francisella isolates, including 13 F. tularensis subspecies holarctica (type B) strains, 26 F. tularensis subsp. tularensis (type A) strains and a single F. novicida strain. The tree was generated from global multi-strain single nucleotide polymorphism (SNP) data collected using a set of six Affymetrix GeneChip^® resequencing arrays with the non-repetitive portion of LVS (type B) as the reference sequence complemented with unique sequences of SCHU S4 (type A). Global SNP based phylogenetic clustering was able to resolve all non-related strains. The phylogenetic tree was used to guide the selection of informative SNPs specific to major nodes in the tree for development of a genotyping assay for identification of F. tularensis subspecies and clades. We designed and validated an assay that uses these SNPs to accurately genotype 39 additional F. tularensis strains as type A (A1, A2, A1a or A1b) or type B (B1 or B2).

Conclusion

Whole-genome SNP based clustering was shown to accurately identify SNPs for differentiation of F. tularensis subspecies and clades, emphasizing the potential power and utility of this methodology for selecting SNPs for typing of F. tularensis to the strain level. Additionally, whole genome sequence based SNP information gained from a representative population of strains may be used to perform evolutionary or phylogenetic comparisons of strains, or selection of unique strains for whole-genome sequencing projects.     

Single Nucleotide Polymorphism Typing of Mycobacterium ulcerans Reveals Focal Transmission of Buruli Ulcer in a Highly Endemic Region of Ghana

Buruli ulcer (BU) is an emerging necrotizing disease of the skin and subcutaneous tissue caused by Mycobacterium ulcerans. While proximity to stagnant or slow flowing water bodies is a risk factor for acquiring BU, the epidemiology and mode of M. ulcerans transmission is poorly understood. Here we have used high-throughput DNA sequencing and comparisons of the genomes of seven M. ulcerans isolates that appeared monomorphic by existing typing methods. We identified a limited number of single nucleotide polymorphisms (SNPs) and developed a real-time PCR SNP typing method based on these differences. We then investigated clinical isolates of M. ulcerans on which we had detailed information concerning patient location and time of diagnosis. Within the Densu river basin of Ghana we observed dominance of one clonal complex and local clustering of some of the variants belonging to this complex. These results reveal focal transmission and demonstrate, that micro-epidemiological analyses by SNP typing has great potential to help us understand how M. ulcerans is transmitted.     

Genomic DNA Enrichment Using Sequence Capture Microarrays: a Novel Approach to Discover Sequence Nucleotide Polymorphisms (SNP) in Brassica napus L

Targeted genomic selection methodologies, or sequence capture, allow for DNA enrichment and large-scale resequencing and characterization of natural genetic variation in species with complex genomes, such as rapeseed canola (Brassica napus L., AACC, 2n=38). The main goal of this project was to combine sequence capture with next generation sequencing (NGS) to discover single nucleotide polymorphisms (SNPs) in specific areas of the B. napus genome historically associated (via quantitative trait loci –QTL– analysis) to traits of agronomical and nutritional importance. A 2.1 million feature sequence capture platform was designed to interrogate DNA sequence variation across 47 specific genomic regions, representing 51.2 Mb of the Brassica A and C genomes, in ten diverse rapeseed genotypes. All ten genotypes were sequenced using the 454 Life Sciences chemistry and to assess the effect of increased sequence depth, two genotypes were also sequenced using Illumina HiSeq chemistry. As a result, 589,367 potentially useful SNPs were identified. Analysis of sequence coverage indicated a four-fold increased representation of target regions, with 57% of the filtered SNPs falling within these regions. Sixty percent of discovered SNPs corresponded to transitions while 40% were transversions. Interestingly, fifty eight percent of the SNPs were found in genic regions while 42% were found in intergenic regions. Further, a high percentage of genic SNPs was found in exons (65% and 64% for the A and C genomes, respectively). Two different genotyping assays were used to validate the discovered SNPs. Validation rates ranged from 61.5% to 84% of tested SNPs, underpinning the effectiveness of this SNP discovery approach. Most importantly, the discovered SNPs were associated with agronomically important regions of the B. napus genome generating a novel data resource for research and breeding this crop species.     

Genetic variation in TLR10 is not associated with chronic Q fever,despite the inhibitory effect of TLR10 on Coxiella burnetii-induced cytokines in vitro

Coxiella burnetii, the causative agent of Q fever, is recognized by TLR2. TLR10 can act as an inhibitory receptor on TLR2-derived immune responses. Therefore, we investigated the role of TLR10 on C. burnetii-induced cytokine production and assessed whether genetic polymorphisms in TLR10 influences the development of chronic Q fever. HEK293 cells, transfected with TLR2, TLR10 or TLR2/TLR10, and human peripheral blood mononuclear cells (PBMCs) in the presence of anti-TLR10, were stimulated with C. burnetii. In both assays, the absence of TLR10 resulted in increased cytokine responses after C. burnetii stimulation. In addition, the effect of single nucleotide polymorphisms (SNPs) in TLR10 was examined in healthy volunteers whose PBMCs were stimulated with C. burnetii Nine Mile or the Dutch outbreak isolate C. burnetii 3262. Individuals bearing SNPs in TLR10 displayed increased cytokine production upon C. burnetii 3262 stimulation. Furthermore, 139 chronic Q fever patients and 220 controls were genotyped for TLR10 N241H, I775V and I369L. None of these polymorphisms were associated with increased susceptibility to chronic Q fever. In conclusion, TLR10 has an inhibitory effect on in vitro cytokine production by C. burnetii, but the presence of TLR10 polymorphisms does not lead to an increased risk of developing chronic Q fever.     

Human and mouse switch-like genes share common transcriptional regulatory mechanisms for bimodality

Background

The increasing number of sequenced prokaryotic genomes contains a wealth of genomic data that needs to be effectively analysed. A set of statistical tools exists for such analysis, but their strengths and weaknesses have not been fully explored. The statistical methods we are concerned with here are mainly used to examine similarities between archaeal and bacterial DNA from different genomes. These methods compare observed genomic frequencies of fixed-sized oligonucleotides with expected values, which can be determined by genomic nucleotide content, smaller oligonucleotide frequencies, or be based on specific statistical distributions. Advantages with these statistical methods include measurements of phylogenetic relationship with relatively small pieces of DNA sampled from almost anywhere within genomes, detection of foreign/conserved DNA, and homology searches. Our aim was to explore the reliability and best suited applications for some popular methods, which include relative oligonucleotide frequencies (ROF), di- to hexanucleotide zero'th order Markov methods (ZOM) and 2.order Markov chain Method (MCM). Tests were performed on distant homology searches with large DNA sequences, detection of foreign/conserved DNA, and plasmid-host similarity comparisons. Additionally, the reliability of the methods was tested by comparing both real and random genomic DNA.

Results

Our findings show that the optimal method is context dependent. ROFs were best suited for distant homology searches, whilst the hexanucleotide ZOM and MCM measures were more reliable measures in terms of phylogeny. The dinucleotide ZOM method produced high correlation values when used to compare real genomes to an artificially constructed random genome with similar %GC, and should therefore be used with care. The tetranucleotide ZOM measure was a good measure to detect horizontally transferred regions, and when used to compare the phylogenetic relationships between plasmids and hosts, significant correlation (R ² = 0.4) was found with genomic GC content and intra-chromosomal homogeneity.

Conclusion

The statistical methods examined are fast, easy to implement, and powerful for a number of different applications involving genomic sequence comparisons. However, none of the measures examined were superior in all tests, and therefore the choice of the statistical method should depend on the task at hand.     

Presence and Persistence of Coxiella burnetii in the Environments of Goat Farms Associated with a Q Fever Outbreak

Q fever is a zoonotic disease caused by inhalation of the bacterium Coxiella burnetii. Ruminant livestock are common reservoirs for C. burnetii, and bacteria present in aerosols derived from the waste of infected animals can infect humans. The significance of infection from material deposited in the environment versus transmission directly from infected animals is not known. In 2011, an outbreak of Q fever cases on farms in Washington and Montana was associated with infected goats. A study was undertaken to investigate the quantity and spatial distribution of C. burnetii in the environment of these goat farms. Soil, vacuum, and sponge samples collected on seven farms epidemiologically linked to the outbreak were tested for the presence of C. burnetii DNA by quantitative PCR. Overall, 70.1% of the samples were positive for C. burnetii. All farms had positive samples, but the quantity of C. burnetii varied widely between samples and between farms. High quantities of C. burnetii DNA were in goat housing/birthing areas, and only small quantities were found in samples collected more than 50 m from these areas. Follow-up sampling at one of the farms 1 year after the outbreak found small quantities of C. burnetii DNA in air samples and large quantities of C. burnetii persisting in soil and vacuum samples. The results suggest that the highest concentrations of environmental C. burnetii are found in goat birthing areas and that contamination of other areas is mostly associated with human movement.     

Detection of Coxiella burnetii in Ambient Air after a Large Q Fever Outbreak

One of the largest Q fever outbreaks ever occurred in the Netherlands from 2007–2010, with 25 fatalities among 4,026 notified cases. Airborne dispersion of Coxiella burnetii was suspected but not studied extensively at the time. We investigated temporal and spatial variation of Coxiella burnetii in ambient air at residential locations in the most affected area in the Netherlands (the South-East), in the year immediately following the outbreak. One-week average ambient particulate matter < 10 μm samples were collected at eight locations from March till September 2011. Presence of Coxiella burnetii DNA was determined by quantitative polymerase chain reaction. Associations with various spatial and temporal characteristics were analyzed by mixed logistic regression. Coxiella burnetii DNA was detected in 56 out of 202 samples (28%). Airborne Coxiella burnetii presence showed a clear seasonal pattern coinciding with goat kidding. The spatial variation was significantly associated with number of goats on the nearest goat farm weighted by the distance to the farm (OR per IQR: 1.89, CI: 1.31–2.76). We conclude that in the year after a large Q fever outbreak, temporal variation of airborne Coxiella burnetii is suggestive to be associated with goat kidding, and spatial variation with distance to and size of goat farms. Aerosol measurements show to have potential for source identification and attribution of an airborne pathogen, which may also be applicable in early stages of an outbreak.     

Comparative Genomic Analysis of Mycobacterium tuberculosis Drug Resistant Strains from Russia

Tuberculosis caused by multidrug-resistant (MDR) and extensively drug-resistant (XDR) Mycobacterium tuberculosis (MTB) strains is a growing problem in many countries. The availability of the complete nucleotide sequences of several MTB genomes allows to use the comparative genomics as a tool to study the relationships of strains and differences in their evolutionary history including acquisition of drug-resistance. In our work, we sequenced three genomes of Russian MTB strains of different phenotypes – drug susceptible, MDR and XDR. Of them, MDR and XDR strains were collected in Tomsk (Siberia, Russia) during the local TB outbreak in 1998–1999 and belonged to rare KQ and KY families in accordance with IS6110 typing, which are considered endemic for Russia. Based on phylogenetic analysis, our isolates belonged to different genetic families, Beijing, Ural and LAM, which made the direct comparison of their genomes impossible. For this reason we performed their comparison in the broader context of all M. tuberculosis genomes available in GenBank. The list of unique individual non-synonymous SNPs for each sequenced isolate was formed by comparison with all SNPs detected within the same phylogenetic group. For further functional analysis, all proteins with unique SNPs were ascribed to 20 different functional classes based on Clusters of Orthologous Groups (COG). We have confirmed drug resistant status of our isolates that harbored almost all known drug-resistance associated mutations. Unique SNPs of an XDR isolate CTRI-4^XDR, belonging to a Beijing family were compared in more detail with SNPs of additional 14 Russian XDR strains of the same family. Only type specific mutations in genes of repair, replication and recombination system (COG category L) were found common within this group. Probably the other unique SNPs discovered in CTRI-4^XDR may have an important role in adaptation of this microorganism to its surrounding and in escape from antituberculosis drugs treatment.     

Identification of Novel Coxiella burnetii Genotypes from Ethiopian Ticks

Background

Coxiella burnetii, the etiologic agent of Q fever, is a highly infectious zoonotic bacterium. Genetic information about the strains of this worldwide distributed agent circulating on the African continent is limited. The aim of the present study was the genetic characterization of C. burnetii DNA samples detected in ticks collected from Ethiopian cattle and their comparison with other genotypes found previously in other parts of the world.

Methodology/Principal Findings

A total of 296 tick samples were screened by real-time PCR targeting the IS1111 region of C. burnetii genome and from the 32 positive samples, 8 cases with sufficient C. burnetii DNA load (Amblyomma cohaerens, n = 6; A. variegatum, n = 2) were characterized by multispacer sequence typing (MST) and multiple-locus variable-number tandem repeat analysis (MLVA). One novel sequence type (ST), the proposed ST52, was identified by MST. The MLVA-6 discriminated the proposed ST52 into two newly identified MLVA genotypes: type 24 or AH was detected in both Amblyomma species while type 26 or AI was found only in A. cohaerens.

Conclusions/Significance

Both the MST and MLVA genotypes of the present work are closely related to previously described genotypes found primarily in cattle samples from different parts of the globe. This finding is congruent with the source hosts of the analyzed Ethiopian ticks, as these were also collected from cattle. The present study provides genotype information of C. burnetii from this seldom studied East-African region as well as further evidence for the presumed host-specific adaptation of this agent.     

Detection of <Emphasis Type="Italic">Coxiella burnetii</Emphasis> DNA and anti-<Emphasis Type="Italic">Coxiella burnetii</Emphasis> IgG antibodies in precolostral blood samples of stillborn calves in an endemically infected Holstein dairy herd

Coxiella burnetii (C. burnetii), an intracellular zoonotic bacterium causing Q fever, occurs widely in cattle herds. After invasion of the pregnant uterus and initial localization in the placenta, active C. burnetii infections may spread to the fetus hematogenously or by the amniotic-oral route and thus may cause abortion, premature delivery, stillbirth, and weak offspring (APSW) complex. In a case-control study, we investigated precolostral blood samples of 56 stillborn calves and 30 live births from a dairy herd endemically infected with C. burnetii “C-cluster” strains and an increased stillbirth rate in primiparous cows. Within the group of the stillborn calves, four precolostral blood samples (7.1%) were tested positive for C. burnetii DNA by PCR and one serum sample (1.8%) positive for anti-C. burnetii IgG antibodies by a commercial ELISA test, respectively. Neither C. burnetii DNA nor anti-C. burnetii IgG antibodies were detected in the samples of calves being born alive. In conclusion, we demonstrated that coxiellaemia and precolostral seroconversion occurred sporadically in stillborn calves from this endemically infected herd. Due to the low detection rates, C. burnetii could not be confirmed to be the cause of the increased stillbirth rate.     

Complete genome sequences of probiotic strains Bifidobacterium animalis subsp. lactis B420 and Bi-07

We present the complete genomes of Bifidobacterium animalis subsp. lactis B420 and Bi-07. Comparative genomic analysis with the type strain DSMZ10140 revealed 40 to 55 single nucleotide polymorphisms (SNPs) and an indel in a clustered regularly interspaced short palindromic repeat (CRISPR) locus. These genetic differences provide a molecular basis for strain typing within the two main phylogenetic groups of this monomorphic species.     

PCR Detection of Coxiella burnetii from Different Clinical Specimens, Especially Bovine Milk, on the Basis of DNA Preparation with a Silica Matrix

For PCR detection of Coxiella burnetii in various clinical specimens we developed a sample preparation method in which silica binding of DNA was used. This method was found to be fast, easily performed with large numbers of samples, and equally sensitive for all of the specimens tested (livers, spleens, placentas, heart valves, milk, blood). The DNA preparation method described here can also be used as an initial step in any PCR-based examination of specimens. The procedure was tested with more than 600 milk samples, which were taken from 21 cows that were seropositive for C. burnetii and reportedly had fertility problems (and therefore were suspected of shedding the agent through milk intermittently or continuously). Of the 21 cows tested, 6 were shedding C. burnetii through milk. Altogether, C. burnetii DNA was detected in 6% of the samples. There was no correlation between the shedding pattern and the serological results.     

First historical genome of a crop bacterial pathogen from herbarium specimen: Insights into citrus canker emergence

Over the past decade, ancient genomics has been used in the study of various pathogens. In this context, herbarium specimens provide a precious source of dated and preserved DNA material, enabling a better understanding of plant disease emergences and pathogen evolutionary history. We report here the first historical genome of a crop bacterial pathogen, Xanthomonas citri pv. citri (Xci), obtained from an infected herbarium specimen dating back to 1937. Comparing the 1937 genome within a large set of modern genomes, we reconstructed their phylogenetic relationships and estimated evolutionary parameters using Bayesian tip-calibration inferences. The arrival of Xci in the South West Indian Ocean islands was dated to the 19^th century, probably linked to human migrations following slavery abolishment. We also assessed the metagenomic community of the herbarium specimen, showed its authenticity using DNA damage patterns, and investigated its genomic features including functional SNPs and gene content, with a focus on virulence factors.     

Genome Sequences of Populus tremula Chloroplast and Mitochondrion: Implications for Holistic Poplar Breeding

Complete Populus genome sequences are available for the nucleus (P. trichocarpa; section Tacamahaca) and for chloroplasts (seven species), but not for mitochondria. Here, we provide the complete genome sequences of the chloroplast and the mitochondrion for the clones P. tremula W52 and P. tremula x P. alba 717-1B4 (section Populus). The organization of the chloroplast genomes of both Populus clones is described. A phylogenetic tree constructed from all available complete chloroplast DNA sequences of Populus was not congruent with the assignment of the related species to different Populus sections. In total, 3,024 variable nucleotide positions were identified among all compared Populus chloroplast DNA sequences. The 5-prime part of the LSC from trnH to atpA showed the highest frequency of variations. The variable positions included 163 positions with SNPs allowing for differentiating the two clones with P. tremula chloroplast genomes (W52, 717-1B4) from the other seven Populus individuals. These potential P. tremula-specific SNPs were displayed as a whole-plastome barcode on the P. tremula W52 chloroplast DNA sequence. Three of these SNPs and one InDel in the trnH-psbA linker were successfully validated by Sanger sequencing in an extended set of Populus individuals. The complete mitochondrial genome sequence of P. tremula is the first in the family of Salicaceae. The mitochondrial genomes of the two clones are 783,442 bp (W52) and 783,513 bp (717-1B4) in size, structurally very similar and organized as single circles. DNA sequence regions with high similarity to the W52 chloroplast sequence account for about 2% of the W52 mitochondrial genome. The mean SNP frequency was found to be nearly six fold higher in the chloroplast than in the mitochondrial genome when comparing 717-1B4 with W52. The availability of the genomic information of all three DNA-containing cell organelles will allow a holistic approach in poplar molecular breeding in the future.