Correction to: Nuclear, chloroplast, and mitochondrial data of a US cannabis DNA database.

The original version of this article contained a mistake. In page 10 of the original article, the significant level (p > 0.01) is incorrect. The correct significant level is (p < 0.01). The original article has been corrected.

Nuclear, chloroplast, and mitochondrial data of a US cannabis DNA database.

As Cannabis sativa (marijuana) is a controlled substance in many parts of the world, the ability to track biogeographical origin of cannabis could provide law enforcement with investigative leads regarding its trade and distribution. Population substructure and inbreeding may cause cannabis plants to become more genetically related. This genetic relatedness can be helpful for intelligence purposes. Analysis of autosomal, chloroplast, and mitochondrial DNA allows for not only prediction of biogeographical origin of a plant but also discrimination between individual plants. A previously validated, 13-autosomal STR multiplex was used to genotype 510 samples. Samples were analyzed from four different sites: 21 seizures at the US-Mexico border, Northeastern Brazil, hemp seeds purchased in the US, and the Araucania area of Chile. In addition, a previously reported multi-loci system was modified and optimized to genotype five chloroplast and two mitochondrial markers. For this purpose, two methods were designed: a homopolymeric STR pentaplex and a SNP triplex with one chloroplast (Cscp001) marker shared by both methods for quality control. For successful mitochondrial and chloroplast typing, a novel real-time PCR quantitation method was developed and validated to accurately estimate the quantity of the chloroplast DNA (cpDNA) using a synthetic DNA standard. Moreover, a sequenced allelic ladder was also designed for accurate genotyping of the homopolymeric STR pentaplex. For autosomal typing, 356 unique profiles were generated from the 425 samples that yielded full STR profiles and 25 identical genotypes within seizures were observed. Phylogenetic analysis and case-to-case pairwise comparisons of 21 seizures at the US-Mexico border, using the Fixation Index (F ST ) as genetic distance, revealed the genetic association of nine seizures that formed a reference population. For mitochondrial and chloroplast typing, subsampling was performed, and 134 samples were genotyped. Complete haplotypes (STRs and SNPs) were observed for 127 samples. As expected, extensive haplotype sharing was observed; five distinguishable haplotypes were detected. In the reference population, the same haplotype was observed 39 times and two unique haplotypes were also detected. Haplotype sharing was observed between the US border seizures, Brazil, and Chile, while the hemp samples generated a distinct haplotype. Phylogenetic analysis of the four populations was performed, and results revealed that both autosomal and lineage markers could discern population substructure.

Evaluation of a 13-loci STR multiplex system for Cannabis sativa genetic identification.

Marijuana (Cannabis sativa) is the most commonly used illicit substance in the USA. The development of a validated method using Cannabis short tandem repeats (STRs) could aid in the individualization of samples as well as serve as an intelligence tool to link multiple cases. For this purpose, a modified 13-loci STR multiplex method was optimized and evaluated according to ISFG and SWGDAM guidelines. A real-time PCR quantification method for C. sativa was developed and validated, and a sequenced allelic ladder was also designed to accurately genotype 199 C. sativa samples from 11 U.S. Customs and Border Protection seizures. Distinguishable DNA profiles were generated from 127 samples that yielded full STR profiles. Four duplicate genotypes within seizures were found. The combined power of discrimination of this multilocus system is 1 in 70 million. The sensitivity of the multiplex STR system is 0.25 ng of template DNA. None of the 13 STR markers cross-reacted with any of the studied species, except for Humulus lupulus (hops) which generated unspecific peaks. Phylogenetic analysis and case-to-case pairwise comparison of 11 cases using F st as genetic distance revealed the genetic association of four groups of cases. Moreover, due to their genetic similarity, a subset of samples (N = 97) was found to form a homogeneous population in Hardy-Weinberg and linkage equilibrium. The results of this research demonstrate the applicability of this 13-loci STR system in associating Cannabis cases for intelligence purposes.

The impact of storm events on a riverbed system and its hydraulic conductivity at a site of induced infiltration.

The spatial and temporal variability of riverbed vertical hydraulic conductivity (K(v)) was investigated at a site of induced infiltration, associated with a municipal well field, to assess the impact of high-stage events on scour and subsequently the riverbed K(v). Such impacts are important when considering the potential loss of riverbank filtration capacity due to storm events. The study site, in and along the Great Miami River in southwest Ohio, overlaid a highly productive glacial-outwash aquifer. A three-layer model for this system was conceptualized: a top layer of transient sediment, a second layer comprising large sediment resistant to scour, but clogged with finer sediment (the armor/colmation layer), and a third layer that was transitional to the underlying higher-K(v) aquifer. One location was studied in detail to confirm and quantify the conceptual model. Methods included seepage meters, heat-flow modeling, grain-size analyses, laboratory permeameter tests, slug tests and the use of scour chains and pressure-load cells to directly measure the amount of sediment scour and re-deposition. Seepage meter measured riverbed K(v) ranged from 0.017 to 1.7 m/d with a geometric mean of 0.19 m/d. Heat-transport model-calibrated estimates were even lower, ranging from 0.0061 to 0.046 m/d with a mean of 0.017 m/d. The relatively low K(v) was indicative of the clogged armor layer. In contrast, slug tests in the underlying riverbed sediment yielded K(v) values an order of magnitude greater. There was a linear relationship between scour chain measured scour and event intensity with a maximum scour of only 0.098 m. Load-cell pressure sensor data over a 7-month period indicated a total sediment-height fluctuation of 0.42 m and a maximum storm-event scour of 0.28 m. Scour data indicated that the assumed armor/colmation layer almost always remained intact. Based on measured layer conductivities and thicknesses, the overall K(v) of this conceptualized system was 1.6 m/d. Sensitivity analyses indicated that even complete scour of the armor/colmation layer would likely increase the overall K(v) only by a factor of 1.5. Most scour events observed removed only the transient sediment, having very little effect on the entire system indicating low risk of losing filtration capacity during storms. The research, however, focused on the point bar, depositional side of the river. More research of the entire river profile is necessary.

Developmental and internal validation of a novel 13 loci STR multiplex method for Cannabis sativa DNA profiling.

Marijuana (Cannabis sativa L.) is a plant cultivated and trafficked worldwide as a source of fiber (hemp), medicine, and intoxicant. The development of a validated method using molecular techniques such as short tandem repeats (STRs) could serve as an intelligence tool to link multiple cases by means of genetic individualization or association of cannabis samples. For this purpose, a 13 loci STR multiplex method was developed, optimized, and validated according to relevant ISFG and SWGDAM guidelines. The STR multiplex consists of 13 previously described C. sativa STR loci: ANUCS501, 9269, 4910, 5159, ANUCS305, 9043, B05, 1528, 3735, CS1, D02, C11, and H06. A sequenced allelic ladder consisting of 56 alleles was designed to accurately genotype 101 C. sativa samples from three seizures provided by a U.S. Customs and Border Protection crime lab. Using an optimal range of DNA (0.5-1.0ng), validation studies revealed well-balanced electropherograms (inter-locus balance range: 0.500-1.296), relatively balanced heterozygous peaks (mean peak height ratio of 0.83 across all loci) with minimal artifacts and stutter ratio (mean stutter of 0.021 across all loci). This multi-locus system is relatively sensitive (0.13ng of template DNA) with a combined power of discrimination of 1 in 55 million. The 13 STR panel was found to be species specific for C. sativa; however, non-specific peaks were produced with Humulus lupulus. The results of this research demonstrate the robustness and applicability of this 13 loci STR system for forensic DNA profiling of marijuana samples.

Development and validation of a multi-locus DNA metabarcoding method to identify endangered species in complex samples.

DNA metabarcoding provides great potential for species identification in complex samples such as food supplements and traditional medicines. Such a method would aid Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES) enforcement officers to combat wildlife crime by preventing illegal trade of endangered plant and animal species. The objective of this research was to develop a multi-locus DNA metabarcoding method for forensic wildlife species identification and to evaluate the applicability and reproducibility of this approach across different laboratories. A DNA metabarcoding method was developed that makes use of 12 DNA barcode markers that have demonstrated universal applicability across a wide range of plant and animal taxa and that facilitate the identification of species in samples containing degraded DNA. The DNA metabarcoding method was developed based on Illumina MiSeq amplicon sequencing of well-defined experimental mixtures, for which a bioinformatics pipeline with user-friendly web-interface was developed. The performance of the DNA metabarcoding method was assessed in an international validation trial by 16 laboratories, in which the method was found to be highly reproducible and sensitive enough to identify species present in a mixture at 1% dry weight content. The advanced multi-locus DNA metabarcoding method assessed in this study provides reliable and detailed data on the composition of complex food products, including information on the presence of CITES-listed species. The method can provide improved resolution for species identification, while verifying species with multiple DNA barcodes contributes to an enhanced quality assurance.

Neighboring group participation in the transition state of human purine nucleoside phosphorylase.

The X-ray crystal structures of human purine nucleoside phosphorylase (PNP) with bound inosine or transition-state analogues show His257 within hydrogen bonding distance of the 5'-hydroxyl. The mutants His257Phe, His257Gly, and His257Asp exhibited greatly decreased affinity for Immucillin-H (ImmH), binding this mimic of an early transition state as much as 370-fold (Km/Ki) less tightly than native PNP. In contrast, these mutants bound DADMe-ImmH, a mimic of a late transition state, nearly as well as the native enzyme. These results indicate that His257 serves an important role in the early stages of transition-state formation. Whereas mutation of His257 resulted in little variation in the PNP x DADMe-ImmH x SO4 structures, His257Phe x ImmH x PO4 showed distortion at the 5'-hydroxyl, indicating the importance of H-bonding in positioning this group during progression to the transition state. Binding isotope effect (BIE) and kinetic isotope effect (KIE) studies of the remote 5'-(3)H for the arsenolysis of inosine with native PNP revealed a BIE of 1.5% and an unexpectedly large intrinsic KIE of 4.6%. This result is interpreted as a moderate electronic distortion toward the transition state in the Michaelis complex with continued development of a similar distortion at the transition state. The mutants His257Phe, His257Gly, and His257Asp altered the 5'-(3)H intrinsic KIE to -3, -14, and 7%, respectively, while the BIEs contributed 2, 2, and -2%, respectively. These surprising results establish that forces in the Michaelis complex, reported by the BIEs, can be reversed or enhanced at the transition state.

Targeting a novel Plasmodium falciparum purine recycling pathway with specific immucillins.

Plasmodium falciparum is unable to synthesize purine bases and relies upon purine salvage and purine recycling to meet its purine needs. We report that purines formed as products of polyamine synthesis are recycled in a novel pathway in which 5'-methylthioinosine is generated by adenosine deaminase. The action of P. falciparum purine nucleoside phosphorylase is a convergent step of purine salvage, converting both 5'-methylthioinosine and inosine to hypoxanthine. We used accelerator mass spectrometry to verify that 5'-methylthioinosine is an active nucleic acid precursor in P. falciparum. Prior studies have shown that inhibitors of purine salvage enzymes kill malaria, but potent malaria-specific inhibitors of these enzymes have not been described previously. 5'-Methylthio-immucillin-H, a transition state analogue inhibitor that is selective for malarial relative to human purine nucleoside phosphorylase, kills P. falciparum in culture. Immucillins are currently in clinical trials for other indications and may also have application as anti-malarials.

Nucleophilic participation in the transition state for human thymidine phosphorylase.

Recombinant human thymidine phosphorylase catalyzes the reaction of arsenate with thymidine to form thymine and 2-deoxyribose 1-arsenate, which rapidly decomposes to 2-deoxyribose and inorganic arsenate. The transition-state structure of this reaction was determined using kinetic isotope effect analysis followed by computer modeling. Experimental kinetic isotope effects were determined at physiological pH and 37 degrees C. The extent of forward commitment to catalysis was determined by pulse-chase experiments to be 0.70%. The intrinsic kinetic isotope effects for [1'-(3)H]-, [2'R-(3)H]-, [2'S-(3)H]-, [4'-(3)H]-, [5'-(3)H]-, [1'-(14)C]-, and [1-(15)N]-thymidines were determined to be 0.989 +/- 0.002, 0.974 +/- 0.002, 1.036 +/- 0.002, 1.020 +/- 0.003, 1.061 +/- 0.003, 1.139 +/- 0.005, and 1.022 +/- 0.005, respectively. A computer-generated model, based on density functional electronic structure calculations, was fit to the experimental isotope effect. The structure of the transition state confirms that human thymidine phosphorylase proceeds through an S(N)2-like transition state with bond orders of 0.50 to the thymine leaving group and 0.33 to the attacking oxygen nucleophile. The reaction differs from the dissociative transition states previously reported for N-ribosyl transferases and is the first demonstration of a nucleophilic transition state for an N-ribosyl transferase. The large primary (14)C isotope effect of 1.139 can occur only in nucleophilic displacements and is the largest (14)C primary isotope effect reported for an enzymatic reaction. A transition state structure with substantial bond order to the attacking nucleophile and leaving group is confirmed by the slightly inverse 1'-(3)H isotope effect, demonstrating that the transition state is compressed by the impinging steric bulk of the nucleophile and leaving group.

Binding causes the remote [5'-3H]thymidine kinetic isotope effect in human thymidine phosphorylase.

The remote 5'-3H V/K kinetic isotope effect (KIE) observed in human thymidine phosphorylase (6.1%) is significantly larger than can be explained by the reaction chemistry. One hypothesis connects the 5'-3H KIE in purine nucleoside phosphorylase to that enzyme's SN1 transition state. The transition state of thymidine phosphorylase, however, is an SN2 nucleophilic displacement. Here we report equilibrium binding isotope effects sufficiently large to explain the presence of this substantial KIE in thymidine phosphorylase.

Structure-based design of novel inhibitors of 3-deoxy-D-manno-octulosonate 8-phosphate synthase.

3-Deoxy-D-manno-octulosonate 8-phosphate (KDO8P) is the phosphorylated precursor of KDO, an essential sugar of the lipopolysaccharide of Gram negative bacteria. KDO8P is produced by a specific synthase (KDO8PS) by condensing arabinose 5-phosphate (A5P) and phosphoenolpyruvate (PEP), with release of inorganic phosphate. As KDO8PS is present in bacteria and plants, but not in mammalian cells, and mutations that inactivate KDO8PS also block cell replication, KDO8PS is a promising target for the design of new antimicrobials that act by blocking lipopolysaccharide biosynthesis. Previous studies have shown that a compound mimicking an intermediate of the condensation reaction is a good ligand and a powerful inhibitor. Here we report on the crystallographic investigation of the binding to KDO8PS of new derivatives of this original inhibitor. The structures of the enzyme in complex with these compounds, and also with the PEP analogs, 2-phosphoglyceric acid (2-PGA) and Z-methyl-PEP, point to future strategies for the design of novel inhibitors of KDO8PS.

KpsF is the arabinose-5-phosphate isomerase required for 3-deoxy-D-manno-octulosonic acid biosynthesis and for both lipooligosaccharide assembly and capsular polysaccharide expression in Neisseria meningitidis.

We have identified and defined the function of kpsF of Neisseria meningitidis and the homologues of kpsF in encapsulated K1 and K5 Escherichia coli. KpsF was shown to be the arabinose-5-phosphate isomerase, an enzyme not previously identified in prokaryotes, that mediates the interconversion of ribulose 5-phosphate and arabinose 5-phosphate. KpsF is required for 3-deoxy-d-manno-octulosonic acid (Kdo) biosynthesis in N. meningitidis. Mutation of kpsF or the gene encoding the CMP-Kdo synthetase (kpsU/kdsB) in N. meningitidis resulted in expression of a lipooligosaccharide (LOS) structure that contained only lipid A and reduced capsule expression in the five invasive disease-associated meningococcal serogroups (A, B, C, Y, and W-135). The step linking meningococcal capsule and LOS biosynthesis was shown to be Kdo production as the expression of capsule was wild type in a Kdo transferase (kdtA) mutant. Thus, in addition to lipooligosaccharide assembly, Kdo is required for meningococcal capsular polysaccharide expression. Furthermore, N. meningitidis, unlike enteric Gram-negative bacteria, can survive and synthesize only unglycosylated lipid A.