RFPlasmid: predicting plasmid sequences from short-read assembly data using machine learning

Antimicrobial-resistance (AMR) genes in bacteria are often carried on plasmids and these plasmids can transfer AMR genes between bacteria. For molecular epidemiology purposes and risk assessment, it is important to know whether the genes are located on highly transferable plasmids or in the more stable chromosomes. However, draft whole-genome sequences are fragmented, making it difficult to discriminate plasmid and chromosomal contigs. Current methods that predict plasmid sequences from draft genome sequences rely on single features, like k-mer composition, circularity of the DNA molecule, copy number or sequence identity to plasmid replication genes, all of which have their drawbacks, especially when faced with large single-copy plasmids, which often carry resistance genes. With our newly developed prediction tool RFPlasmid, we use a combination of multiple features, including k-mer composition and databases with plasmid and chromosomal marker proteins, to predict whether the likely source of a contig is plasmid or chromosomal. The tool RFPlasmid supports models for 17 different bacterial taxa, including Campylobacter , Escherichia coli and Salmonella , and has a taxon agnostic model for metagenomic assemblies or unsupported organisms. RFPlasmid is available both as a standalone tool and via a web interface.

Screening of possibly anthrax-contaminated burial sites in eastern and southern Ukraine

Aim. The aim of this study was to screen soil samples of 17 anthrax burial sites in Eastern and Southern Ukraine for the presence of B. anthracis. Methods. Soil samples were collected from anthrax grave sites located in Kharkiv, Sumy and Mykolaiv regions (diseased animals dated from 1946 to 2003). Isolation of B. anthracis from collected soil samples was performed with the GABRI method. From single colonies without hemolysis, that were inactivated with peracetic acid- containing 2 % Terralin PAA solution, DNA was extracted and analyzed by qPCR for the presence of chromosomal marker dhp61, as well as the markers pagA and capC located on virulence plasmids pXO1 and pXO2, respectively. Results. Eleven fi eld trips were conducted from July, 2016 to October, 2018 in which 369 soil samples from 17 burial sites in Kharkiv, Sumy and Mykolaiv oblasts were collected from different depths of presumed anthrax carcass sites. In most cases (12 out of 17 cases), the current status of these burial sites was deteriorated and not prop- erly accounted for. It was possible to obrain viable B. anthracis isolate was obtained from 50 cm depth at the grave site near Koviagy village, Valky district, Kharkiv region (49.92373°N, 35.48951°E). This isolate was named KhR/ VD/Kov2-2-05-3 and deposited in the Collection of Animal Infectious Pathogens of the National Scientifi c Center “Institute of Experimental and Clinical Veterinary Medicine”, Kharkiv, Ukraine. The contamination level of soil at the isolation site reached about 10 4 CFU per g as determined by plate counting. qPCR analysis of this isolate identi- fi ed both the dhp61 B. anthracis chromosomal and the pagA virulence plasmid marker. However, the plasmid pXO2 marker, required for capsule-formation could not be detected. Conclusions. The anthrax burial sites were created between the 1920s and 1960s, however, only approximate locations could be found and demarcated. In most cases the status of the sites was unsuitable for sampling. Nevertheless, isolation of B. anthracis in one case in the Valky district shows that old anthrax burial sites (13.500 exist in Ukraine) still pose a risk as potential source of the infection and therefore require more attention and surveillance, for which a surveillance plan will be developed.

Male recombination map of the autosomal genome in German Holstein

Background Recombination is a process by which chromosomes are broken and recombine to generate new combinations of alleles, therefore playing a major role in shaping genome variation. Recombination frequencies ($$\theta$$ θ ) between markers are used to construct genetic maps, which have important implications in genomic studies. Here, we report a recombination map for 44,696 autosomal single nucleotide polymorphisms (SNPs) according to the coordinates of the most recent bovine reference assembly. The recombination frequencies were estimated across 876 half-sib families with a minimum number of 39 and maximum number of 4236 progeny, comprising over 367 K genotyped German Holstein animals. Results Genome-wide, over 8.9 million paternal recombination events were identified by investigating adjacent markers. The recombination map spans 24.43 Morgan (M) for a chromosomal length of 2486 Mbp and an average of ~ 0.98 cM/Mbp, which concords with the available pedigree-based linkage maps. Furthermore, we identified 971 putative recombination hotspot intervals (defined as $$\theta$$ θ > 2.5 standard deviations greater than the mean). The hotspot regions were non-uniformly distributed as sharp and narrow peaks, corresponding to ~ 5.8% of the recombination that has taken place in only ~ 2.4% of the genome. We verified genetic map length by applying a likelihood-based approach for the estimation of recombination rate between all intra-chromosomal marker pairs. This resulted in a longer autosomal genetic length for male cattle (25.35 cM) and in the localization of 51 putatively misplaced SNPs in the genome assembly. Conclusions Given the fact that this map is built on the coordinates of the ARS-UCD1.2 assembly, our results provide the most updated genetic map yet available for the cattle genome.

RFPlasmid: Predicting plasmid sequences from short read assembly data using machine learning

Antimicrobial resistance (AMR) genes in bacteria are often carried on plasmids and these plasmids can transfer AMR genes between bacteria. For molecular epidemiology purposes and risk assessment, it is important to know if the genes are located on highly transferable plasmids or in the more stable chromosomes. However, draft whole genome sequences are fragmented, making it difficult to discriminate plasmid and chromosomal contigs. Current methods that predict plasmid sequences from draft genome sequences rely on single features, like k-mer composition, circularity of the DNA molecule, copy number or sequence identity to plasmid replication genes, all of which have their drawbacks, especially when faced with large single copy plasmids, which often carry resistance genes. With our newly developed prediction tool RFPlasmid, we use a combination of multiple features, including k-mer composition and databases with plasmid and chromosomal marker proteins, to predict if the likely source of a contig is plasmid or chromosomal. The tool RFPlasmid supports models for 17 different bacterial species, including Campylobacter, E. coli, and Salmonella, and has a species agnostic model for metagenomic assemblies or unsupported organisms. RFPlasmid is available both as standalone tool and via web interface.

Cytogenetic characterization of Ameivula ocellifera (Spix, 1825) (Squamata, Teiidae) from the brazilian northeast

Ameivula is as a new genus of Teiidae family that emerged after extensive revision of species that comprised the former complex of species called Cnemidophorus group. Its species has a wide distribution from the northeast of Brazil to northern Argentina. Cytogenetic studies in the Teiidae family have shown that karyotypical data are important tools in phylogenetic and systematic studies within this group allowing to determine the position of species in the family. Thus, this study aimed to describe the karyotype of Ameivula ocellifera (Spix, 1825) from Picos, Piauí state in the Brazilian Northeast. Specimens were collected from August 2014 to October 2015 using interception traps and pitfalls, mounted randomly along the Caatinga area. The animals were collected and transported to Federal Institute of Piauí, campus Picos, where was carried out all laboratory procedures. Individuals analyzed showed a diploid number of 2n = 50 for both sexes, with karyotype composed by 30 macrochromosomes and 20 microchromosomes of telocentric and subtelocentric types. There were no heteromorphic sex chromosomes in the studied specimens. C-band technique evidenced the heterochromatic blocks in pericentromeric and telomeric regions of chromosomes. The nucleolar organizing regions appeared as a simple unit located at the terminal portion of the long arm of chromosomal pair number 5. The chromosomal characteristics of A. ocellifera analyzed do not show divergences regarding individuals from other regions. However, the nucleolar organizing regions seems to be a good chromosomal marker that permits to distinguish the species already studied.

Bovine abortion by a vaccine strain of Bacillus anthracis

ABSTRACT: This paper reports the abortion of a male Aberdeen Angus bovine by a vaccine strain of Bacillus anthracis, describing the pathological and microbiological findings and the genome sequence. Necropsy findings included multifocal areas of hemorrhage in different organs. Histologically, various organs showed hemorrhage, fibrin exudation, necrosis associated with countless bacillary bacterial clumps and severe neutrophilic inflammatory infiltrate. In the microbiological examination, numerous rough, nonhemolytic, gray and dry colonies with irregular edges were isolated from liver, lung and abomasum content samples. Gram staining revealed square-ended Gram-positive rods arranged in chains. B. anthracis identification was confirmed by detection of the molecular chromosomal marker Ba813. The genomes from the isolated B. anthracis (named SPV842_15) and from the isolated vaccinal strain (Brazilian vaccinal strain), which was recovered from a commercial vaccine used in the pregnant cow, were sequenced. Genomic comparisons displayed a high level of nucleotide identity in the comparisons between B. anthracis SPV842_15 and the B. anthracis Brazilian vaccinal strain (98,2%). Furthermore, in both strains, only the plasmid pX01 sequence was detected. Although, vaccination against anthrax is characterized by an elevated protective profile and very low residual virulence, immunization with Sterne strains can cause abortion in cattle, presumably by the plasmid pX01 toxins in rare or special situations.

Introduction of Noninvasive Prenatal Testing for Blood Group and Platelet Antigens from Cell-Free Plasma DNA Using Digital PCR

Background: Noninvasive prenatal testing (NIPT) for fetal antigens is a common standard for targeted immune prophylaxis in RhD-mediated hemolytic disease of the fetus and newborn, and is most frequently done by quantitative PCR (qPCR). A similar approach is considered for other blood group and human platelet alloantigens (HPA). Because of a higher sensitivity compared to qPCR for rare molecule detection, we established and validated digital PCR (dPCR) assays for the detection of RHD exons 3, 5 and 7, KEL1, HPA-1a, and HPA-5b from cell-free DNA (cfDNA) in plasma. The dPCR assays for the Y-chromosomal marker amelogenin and autosomal SNPs were implemented as controls for the proof of fetal DNA. Methods: Validation was performed on dilution series of mixed plasma samples from volunteer donors with known genotypes. After preamplification of the target loci, two-color (FAM and VIC) TaqManTM probe chemistry and chip-based dPCR were applied. The assays for RHD included GAPDH as an internal control. For the diallelic markers KEL1/2, HPA-1a/b, HPA-5a/b, and AMEL-X/Y and 3 autosomal SNPs, the probes enabled allelic discrimination in the two fluorescence channels. The dPCR protocol for NIPT was applied to plasma samples from pregnant women. Results: The RHD exon 5 assay allowed the detection of a 0.05% RHD target in an RhD-negative background, whereas the exon 7 assay required at least a 0.25% target. The exon 3 assay showed the highest background and required at least a 2.5% RHD target for reliable detection. The dPCR assays for the diallelic markers revealed similar sensitivity and enabled the detection of at least a 0.5% target allele. The HPA-1a assay was the most sensitive and allowed target detection in plasma mixtures containing only 0.05% HPA-1a. The plasma samples from 13 pregnant women at different gestational ages showed unambiguous positive and negative results for the analyzed targets. Conclusion: Analysis of cfDNA from maternal plasma using dPCR is suitable for the detection of fetal alleles. Because of the high sensitivity of the assays, the NIPT protocol for RhD, KEL1, and HPA can also be applied to earlier stages of pregnancy.

Sex determination using circulating cell-free fetal DNA in small volume of maternal plasma in elephants

The genetic sexing of animals having long gestation periods offers significant benefits in regard to breeding management among their populations living in captivity. In our study, a new increased-sensitivity PCR method for fetal sexing was developed and tested successfully on elephants, from only a small volume of maternal plasma. Suitable sensitivity was obtained by using short, reduced amplicon lengths with fluorescent labelling for capillary electrophoresis detection. The fundamental principle for this technique was based on the detection of two Y-specific markers (AmelY and SRY), the presence of which indicates the mother is carrying a male fetus and the absence of these markers designates a female fetus. As a reaction control, the X-chromosomal marker (PlpX) was used. To the best of our knowledge, this is the first report on this topic, confirming the presence of fetal cell-free DNA from the plasma of a pregnant captive elephant, and demonstrating a new opportunity for non-invasive assessment in fetal sex determination.

The recombination-cold region as an epidemiological marker of recombinogenic opportunistic pathogen Mycobacterium avium

Background The rapid identification of lineage remains a challenge in the genotyping of clinical isolates of recombinogenic pathogens. The chromosome of Mycobacterium avium subsp. hominissuis (MAH), an agent of Mycobacterium avium complex (MAC) lung disease, is often mosaic and is composed of chromosomal segments originating from different lineages. This makes it difficult to infer the MAH lineage in a simple experimental set-up. To overcome this difficulty, we sought to identify chromosomal marker genes containing lineage-specific alleles by genome data mining. Results We conducted genetic population structure analysis, phylogenetic analysis, and a survey of historical recombination using data from 125 global MAH isolates. Six MAH lineages (EA1, EA2, SC1, SC2, SC3, and SC4) were identified in the current dataset. One P-450 gene (locus_tag MAH_0788/MAV_0940) in the recombination-cold region was found to have multiple alleles that could discriminate five lineages. By combining the information about allele type from one additional gene, the six MAH lineages as well as other M. avium subspecies were distinguishable. A recombination-cold region of 116 kb contains an insertion hotspot and is flanked by a mammalian cell-entry protein operon where allelic variants have previously been reported to occur. Hence, we speculate that the acquisition of lineage- or strain-specific insertions has introduced homology breaks in the chromosome, thereby reducing the chance of interlineage recombination. Conclusions The allele types of the newly identified marker genes can be used to predict major lineages of M. avium. The single nucleotide polymorphism typing approach targeting multiallelic loci in recombination-cold regions will facilitate the epidemiological study of MAC, and may also be useful for equivalent studies of other nontuberculous mycobacteria potentially carrying mosaic genomes.