The complete mitochondrial genome of Eumeta variegata (Lepidoptera: Psychidae)

The complete mitochondrial genome of Eumeta variegate, largest bagworm moth in Japan, has been sequenced using a nanopore sequencer as a single long read. The genome has a total length of 16,601 bp, consisting of 13 protein-coding genes, 20 tRNA, 2 rRNA genes, and an AT-rich control region. The nucleotide composition was extremely AT-rich, with 42.4% A, 40.4% T, 6.67% G, and 10.6% C. This is the second report of a complete mitochondrial genome of Psychidae, and the sequence information together with a phylogenetic analysis would provide a reference data in the future studies of Lepidoptera and Psychidae.

Base-pairing probability in the microRNA stem region affects the binding and editing specificity of human A-to-I editing enzymes ADAR1-p110 and ADAR2

Adenosine deaminases acting on RNA (ADARs) catalyze the deamination of adenosine (A) to inosine (I). A-to-I RNA editing targets double-stranded RNA (dsRNA), and increases the complexity of gene regulation by modulating base pairing-dependent processes such as splicing, translation, and microRNA (miRNA)-mediated gene silencing. This study investigates the genome-wide binding preferences of the nuclear constitutive isoforms ADAR1-p110 and ADAR2 on human miRNA species by RNA immunoprecipitation of ADAR-bound small RNAs (RIP-seq). Our results suggest that secondary structure predicted by base-pairing probability in the mainly double-stranded region of a pre-miRNA or mature miRNA duplex may determine ADAR isoform preference for binding distinct subpopulations of miRNAs. Furthermore, we identify 31 unique editing sites with statistical significance, 19 sites of which are novel editing sites. Editing sites are enriched in the seed region responsible for target recognition by miRNAs, and isoform-specific nucleotide motifs in the immediate vicinity and opposite of editing sites are consistent with previous studies, and further reveal that ADAR2 may edit A/C bulges more frequently than ADAR1-p110 in the context of miRNA.

Ultralow Input Genome Sequencing Library Preparation from a Single Tardigrade Specimen

Tardigrades are microscopic animals that enter an ametabolic state called anhydrobiosis when facing desiccation and can return to their original state when water is supplied. The genomic sequencing of microscopic animals such as tardigrades risks bacterial contamination that sometimes leads to erroneous interpretations, for example, regarding the extent of horizontal gene transfer in these animals. Here, we provide an ultralow input method to sequence the genome of the tardigrade, Hypsibius dujardini, from a single specimen. By employing rigorous washing and contaminant exclusion along with an efficient extraction of the 50 ~ 200 pg genomic DNA from a single individual, we constructed a library sequenced with a DNA sequencing instrument. These libraries were highly reproducible and unbiased, and an informatics analysis of the sequenced reads with other H. dujardini genomes showed a minimal amount of contamination. This method can be applied to unculturable tardigrades that could not be sequenced using previous methods.

Antithrombotic Mechanisms of Urokinase lmmobilized Polyurethane

Urokinase immobilized polymer is highly antithrombotic, which cannot be explained only by fibrinolysis. We immobllized 10 IU/cm2 of urokinase to polyurethane by using maleic anhydride methylvinyl ether copolymer as a carrier. Then we incubated blood in circular tubes made of this material, measured the clotting factors and observed the surface of the tubes after incubation by scanning electronmicroscopy and immunofluorescence microscopy. After 5 min incubation, the relative activities of factors V, VIII, IX, X and XII, fibrinogeil, plasminogen and α2 plasmin inhibitor decreased, but the activity of factor VII increased. No platelet adhesion to the surface of the urokinase immobilized polyurethane was observed and there was no significant adsorption of serum proteins, including fibrinogen, fibronectin and vWF antigen, or the surface. Urokinase-immobilized polyurethane catalyzed the digestion of clotting factors as well as fibrinolysis and also inhibited platelet adhesion on its surface probably by inhibiting protein adsorption and its clinical application including vessel prosthesis should be developed further.

Abstract 1442: Novel linkage of thymidine catabolism and the glycolytic pathway in human cancer cells

Thymidine phosphorylase (TP), a rate-limiting enzyme in the thymidine catabolic pathway, catalyzes the reversible conversion of thymidine to thymine and 2-deoxy-D-ribose 1-phosphate (DR1P). We previously demonstrated that TP is identical to an anigogenic factor, platelet-derived endothelial cell growth factor (PD-ECGF). Although many reports indicate that TP has angiogenic activity and plays a pivotal role in tumor progression, the mechanism of the TP functions is not fully understood. Here, we find that TP-mediated thymidine catabolism can supply the carbon source in the glycolytic pathway and thus contribute to cell survival under conditions of nutrient deprivation. In certain procaryotes, phosphopentosemutase converts thymidine-derived DR1P to 2-deoxy-D-ribose 5-phosphate (DR5P), which is then converted to glyceraldehyde 3-phosphate (GAP) and acetaldehyde by deoxyriboaldorase (DERA). Given that GAP is an intermediate in the glycolytic pathway and the pentose phosphate pathway (PPP), we hypothesized that the activation of thymidine catabolism by TP could supply the carbon source for those pathways. However, the thymidine catabolic pathway in mammalian cells remains unclarified. We used 13C-labeled thymidine (13C5-thymidine) as a tracer and examined whether 13C5-thymidine-derived metabolites enter the glycolytic pathway and the PPP in TP-expressing cancer cells using capillary electrophoresis time-of-flight mass spectrometry. Thymidine was converted to metabolites, including glucose 6-phosphate, lactate, 5-phospho-α-ribose 1-diphosphate, and serine, via the glycolytic pathway both in vitro and in vivo. These thymidine-derived metabolites were required for the survival of cells under low-glucose conditions. Furthermore, enhanced thymidine catabolism was observed in human gastric cancer. These findings suggest that the thymidine catabolism plays an important role in the survival and growth of cancer cells in a microenvironment where glucose is insufficient.

Citation Format: Sho Tabata, Masatatsu Yamamoto, Hisatsugu Goto, Yasuhiko Nishioka, Masaru Tomita, Tomoyoshi Soga, Tatsuhiko Furukawa, Shin-ichi Akiyama. Novel linkage of thymidine catabolism and the glycolytic pathway in human cancer cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 1442.

Efficient delivery of large DNA from Escherichia coli to Synechococcus elongatus PCC7942 by broad-host-range conjugal plasmid pUB307

Synechococcus elongatus PCC7942, a cyanobacterium that uses light and carbon dioxide to grow, has a high ability to incorporate DNA by transformation. To assess the effective delivery of large DNA in plasmid form, we cloned the endogenous plasmid pANL (46.4 kbp) into a BAC vector of Escherichia coli. The plasmid p38ANL (54.3 kbp) replaced the native plasmid. To assess the delivery of larger DNA into PCC7942, p38ANL was fused to the broad-host-range conjugal transfer plasmid pUB307IP (53.5 kbp). The resulting plasmid pUB307IP501 (107.9 kbp) was transmitted from E. coli to PCC7942 by simple mixing of donor and recipient cultures. PCC7942 transcipients possessed only pUB307IP501, replacing the preexisting pANL. In contrast, the pUB307IP501 plasmid was unable to transform PCC7942, indicating that natural transformation of DNA may be restricted by size limitations. The ability to deliver large DNA by conjugation may lead to genetic engineering in PCC7942.

Sex-specific relevance of diabetes to occlusive vascular and other mortality: a collaborative meta-analysis of individual data from 980 793 adults from 68 prospective studies

BACKGROUND

Several studies have shown that diabetes confers a higher relative risk of vascular mortality among women than among men, but whether this increased relative risk in women exists across age groups and within defined levels of other risk factors is uncertain. We aimed to determine whether differences in established risk factors, such as blood pressure, BMI, smoking, and cholesterol, explain the higher relative risks of vascular mortality among women than among men.

METHODS

In our meta-analysis, we obtained individual participant-level data from studies included in the Prospective Studies Collaboration and the Asia Pacific Cohort Studies Collaboration that had obtained baseline information on age, sex, diabetes, total cholesterol, blood pressure, tobacco use, height, and weight. Data on causes of death were obtained from medical death certificates. We used Cox regression models to assess the relevance of diabetes (any type) to occlusive vascular mortality (ischaemic heart disease, ischaemic stroke, or other atherosclerotic deaths) by age, sex, and other major vascular risk factors, and to assess whether the associations of blood pressure, total cholesterol, and body-mass index (BMI) to occlusive vascular mortality are modified by diabetes.

RESULTS

Individual participant-level data were analysed from 980 793 adults. During 9·8 million person-years of follow-up, among participants aged between 35 and 89 years, 19 686 (25·6%) of 76 965 deaths were attributed to occlusive vascular disease. After controlling for major vascular risk factors, diabetes roughly doubled occlusive vascular mortality risk among men (death rate ratio [RR] 2·10, 95% CI 1·97-2·24) and tripled risk among women (3·00, 2·71-3·33; χ2 test for heterogeneity p<0·0001). For both sexes combined, the occlusive vascular death RRs were higher in younger individuals (aged 35-59 years: 2·60, 2·30-2·94) than in older individuals (aged 70-89 years: 2·01, 1·85-2·19; p=0·0001 for trend across age groups), and, across age groups, the death RRs were higher among women than among men. Therefore, women aged 35-59 years had the highest death RR across all age and sex groups (5·55, 4·15-7·44). However, since underlying confounder-adjusted occlusive vascular mortality rates at any age were higher in men than in women, the adjusted absolute excess occlusive vascular mortality associated with diabetes was similar for men and women. At ages 35-59 years, the excess absolute risk was 0·05% (95% CI 0·03-0·07) per year in women compared with 0·08% (0·05-0·10) per year in men; the corresponding excess at ages 70-89 years was 1·08% (0·84-1·32) per year in women and 0·91% (0·77-1·05) per year in men. Total cholesterol, blood pressure, and BMI each showed continuous log-linear associations with occlusive vascular mortality that were similar among individuals with and without diabetes across both sexes.

INTERPRETATION

Independent of other major vascular risk factors, diabetes substantially increased vascular risk in both men and women. Lifestyle changes to reduce smoking and obesity and use of cost-effective drugs that target major vascular risks (eg, statins and antihypertensive drugs) are important in both men and women with diabetes, but might not reduce the relative excess risk of occlusive vascular disease in women with diabetes, which remains unexplained.

FUNDING

UK Medical Research Council, British Heart Foundation, Cancer Research UK, European Union BIOMED programme, and National Institute on Aging (US National Institutes of Health).

MICOP: Maximal information coefficient-based oscillation prediction to detect biological rhythms in proteomics data

BACKGROUND

Circadian rhythms comprise oscillating molecular interactions, the disruption of the homeostasis of which would cause various disorders. To understand this phenomenon systematically, an accurate technique to identify oscillating molecules among omics datasets must be developed; however, this is still impeded by many difficulties, such as experimental noise and attenuated amplitude.

RESULTS

To address these issues, we developed a new algorithm named Maximal Information Coefficient-based Oscillation Prediction (MICOP), a sine curve-matching method. The performance of MICOP in labeling oscillation or non-oscillation was compared with four reported methods using Mathews correlation coefficient (MCC) values. The numerical experiments were performed with time-series data with (1) mimicking of molecular oscillation decay, (2) high noise and low sampling frequency and (3) one-cycle data. The first experiment revealed that MICOP could accurately identify the rhythmicity of decaying molecular oscillation (MCC > 0.7). The second experiment revealed that MICOP was robust against high-level noise (MCC > 0.8) even upon the use of low-sampling-frequency data. The third experiment revealed that MICOP could accurately identify the rhythmicity of noisy one-cycle data (MCC > 0.8). As an application, we utilized MICOP to analyze time-series proteome data of mouse liver. MICOP identified that novel oscillating candidates numbered 14 and 30 for C57BL/6 and C57BL/6 J, respectively.

CONCLUSIONS

In this paper, we presented MICOP, which is an MIC-based algorithm, for predicting periodic patterns in large-scale time-resolved protein expression profiles. The performance test using artificially generated simulation data revealed that the performance of MICOP for decaying data was superior to that of the existing widely used methods. It can reveal novel findings from time-series data and may contribute to biologically significant results. This study suggests that MICOP is an ideal approach for detecting and characterizing oscillations in time-resolved omics data sets.

Tuning Gene Activity by Inducible and Targeted Regulation of Gene Expression in Minimal Bacterial Cells

Functional genomics studies in minimal mycoplasma cells enable unobstructed access to some of the most fundamental processes in biology. Conventional transposon bombardment and gene knockout approaches often fail to reveal functions of genes that are essential for viability, where lethality precludes phenotypic characterization. Conditional inactivation of genes is effective for characterizing functions central to cell growth and division, but tools are limited for this purpose in mycoplasmas. Here we demonstrate systems for inducible repression of gene expression based on clustered regularly interspaced short palindromic repeats-mediated interference (CRISPRi) in Mycoplasma pneumoniae and synthetic Mycoplasma mycoides, two organisms with reduced genomes actively used in systems biology studies. In the synthetic cell, we also demonstrate inducible gene expression for the first time. Time-course data suggest rapid kinetics and reversible engagement of CRISPRi. Targeting of six selected endogenous genes with this system results in lowered transcript levels or reduced growth rates that agree with lack or shortage of data in previous transposon bombardment studies, and now produces actual cells to analyze. The ksgA gene encodes a methylase that modifies 16S rRNA, rendering it vulnerable to inhibition by the antibiotic kasugamycin. Targeting the ksgA gene with CRISPRi removes the lethal effect of kasugamycin and enables cell growth, thereby establishing specific and effective gene modulation with our system. The facile methods for conditional gene activation and inactivation in mycoplasmas open the door to systematic dissection of genetic programs at the core of cellular life.

Far rapid synthesis of giant DNA in the Bacillus subtilis genome by a conjugation transfer system

Bacillus subtilis offers a platform for giant DNA synthesis, which is mediated by the connection of overlapping DNA segments called domino DNA, in the cloning locus of the host. The domino method was successfully used to produce DNA fragments as large as 3500 kbp. However, domino DNA is limited to <100 kbp because of size restrictions regarding the transformation (TF) of B. subtilis competent cells. A novel conjugal transfer (CT) method was designed to eliminate the TF size limit. The CT method enables rapid and efficient domino reactions in addition to the transfer of giant DNA molecules of up to 875 kbp to another B. subtilis genome within 4 hours. The combined use of the TF and CT should enable significantly rapid giant DNA production.

Changes of liver metabolites following hepatectomy with ischemia reperfusion towards liver regeneration

BACKGROUND

Metabolome analysis is one of the omics which investigates the final product of a central dogma. Changes of liver metabolites during liver regeneration following hepatectomy (Hx) continue to remain unclear. The aim of the present study was to investigate the changes of liver metabolites following Hx with ischemia reperfusion (I/R) towards liver regeneration.

METHODS

Twenty-three patients who underwent Hx were enrolled in this study. Non-tumor tissues were sampled immediately before and after Hx and a comparison was made between the liver samples taken before and after Hx using capillary electrophoresis (CE)-time-of-flight mass spectrometry (TOFMS) as metabolome analysis.

RESULTS

The metabolic pathway showed that there was a significant increase in "lactate" following Hx. There was a significant decrease in metabolites only in the first half of the tricarboxylic acid cycle (TCA) cycle, and adenosine triphosphate (ATP) by anaerobiotic glycolysis did not occur in time for energy consumption of the Hx. Principal component analysis revealed remarkably different component profiles between the samples taken before and after Hx. One hundred and three metabolites were selected as critical metabolites for separating components. Valine and tryptophan increased significantly after Hx and they were regulated by resected liver volume, ischemic time and liver function.

CONCLUSION

The liver metabolites changed remarkably between before and after Hx. Especially, liver valine and tryptophan were increased.

Thymidine catabolism promotes NADPH oxidase-derived reactive oxygen species (ROS) signalling in KB and yumoto cells

Thymidine phosphorylase (TP) is a rate-limiting enzyme in the thymidine catabolic pathway. TP is identical to platelet-derived endothelial cell growth factor and contributes to tumour angiogenesis. TP induces the generation of reactive oxygen species (ROS) and enhances the expression of oxidative stress-responsive genes, such as interleukin (IL)-8. However, the mechanism underlying ROS induction by TP remains unclear. In the present study, we demonstrated that TP promotes NADPH oxidase-derived ROS signalling in cancer cells. NADPH oxidase inhibition using apocynin or small interfering RNAs (siRNAs) abrogated the induction of IL-8 and ROS in TP-expressing cancer cells. Meanwhile, thymidine catabolism induced by TP increased the levels of NADPH and intermediates of the pentose phosphate pathway (PPP). Both siRNA knockdown of glucose 6-phosphate dehydrogenase (G6PD), a rate-limiting enzyme in PPP, and a G6PD inhibitor, dihydroepiandrosterone, reduced TP-induced ROS production. siRNA downregulation of 2-deoxy-D-ribose 5-phosphate (DR5P) aldolase, which is needed for DR5P to enter glycolysis, also suppressed the induction of NADPH and IL-8 in TP-expressing cells. These results suggested that TP-mediated thymidine catabolism increases the intracellular NADPH level via the PPP, which enhances the production of ROS by NADPH oxidase and activates its downstream signalling.

RBC Aggregometer head as a Warning Monitor of Flow Disturbance in Extracorporeal Systems

A simple photoelectric apparatus is described, to monitor the flow condition of whole blood in an extracorporeal system. It can raise an immediate alarm on disturbance of flow. The apparatus comprises an RBC aggregometer head, described previously, for measuring the rate of RBC aggregation in whole blood. Its construction was modified for easy attach-ability and detachability without causing tube damage. In practice, the apparatus was applied by hooking it to the tube of the extracorporeal system in animals and in a clinical case of blood dialysis in a patient with renal failure, in which stoppage of blood flow elicited a dramatic change in the baseline record of light transmission of the mobile blood. The signal of the apparatus was fed to an alarm system via a voltage comparator switch. The apparatus is inexpensive, solid and durable, easy to operate by untrained personnel, and has excellent stability and reproducibility.

Development of a sheathless CE-ESI-MS interface

A sheath-flow interface is the most common ionization technique in CE-ESI-MS. However, this interface dilutes the analytes with the sheath liquid and decreases the sensitivity. In this study, we developed a sheathless CE-MS interface to improve sensitivity. The interface was fabricated by making a small crack approximately 2 cm from the end of a capillary column fixed on a plastic plate, and then covering the crack with a dialysis membrane to prevent metabolite loss during separation. A voltage for CE separation was applied between the capillary inlet and the buffer reservoir. Under optimum conditions, 52 cationic metabolite standards were separated and selectively detected using MS. With a pressure injection of 5 kPa for 15 s (ca. 1.4 nL), the detection limits for the tested compounds were between 0.06 and 1.7 μmol/L (S/N = 3). The method was applied to analysis of cationic metabolites extracted from a small number (12 000) of cancer cells, and the number of peaks detected was about 2.5 times higher than when using conventional sheath-flow CE-MS. Because the interface is easy to construct, it is cost-effective and can be adapted to any commercially available capillaries. This method is a powerful new tool for highly sensitive CE-MS-based metabolomic analysis.

Thymidine Catabolism as a Metabolic Strategy for Cancer Survival

Thymidine phosphorylase (TP), a rate-limiting enzyme in thymidine catabolism, plays a pivotal role in tumor progression; however, the mechanisms underlying this role are not fully understood. Here, we found that TP-mediated thymidine catabolism could supply the carbon source in the glycolytic pathway and thus contribute to cell survival under conditions of nutrient deprivation. In TP-expressing cells, thymidine was converted to metabolites, including glucose 6-phosphate, lactate, 5-phospho-α-D-ribose 1-diphosphate, and serine, via the glycolytic pathway both in vitro and in vivo. These thymidine-derived metabolites were required for the survival of cells under low-glucose conditions. Furthermore, activation of thymidine catabolism was observed in human gastric cancer. These findings demonstrate that thymidine can serve as a glycolytic pathway substrate in human cancer cells.

Human AK2 links intracellular bioenergetic redistribution to the fate of hematopoietic progenitors

AK2 is an adenylate phosphotransferase that localizes at the intermembrane spaces of the mitochondria, and its mutations cause a severe combined immunodeficiency with neutrophil maturation arrest named reticular dysgenesis (RD). Although the dysfunction of hematopoietic stem cells (HSCs) has been implicated, earlier developmental events that affect the fate of HSCs and/or hematopoietic progenitors have not been reported. Here, we used RD-patient-derived induced pluripotent stem cells (iPSCs) as a model of AK2-deficient human cells. Hematopoietic differentiation from RD-iPSCs was profoundly impaired. RD-iPSC-derived hemoangiogenic progenitor cells (HAPCs) showed decreased ATP distribution in the nucleus and altered global transcriptional profiles. Thus, AK2 has a stage-specific role in maintaining the ATP supply to the nucleus during hematopoietic differentiation, which affects the transcriptional profiles necessary for controlling the fate of multipotential HAPCs. Our data suggest that maintaining the appropriate energy level of each organelle by the intracellular redistribution of ATP is important for controlling the fate of progenitor cells.

Reliability of plasma polar metabolite concentrations in a large-scale cohort study using capillary electrophoresis-mass spectrometry

Background

Cohort studies with metabolomics data are becoming more widespread, however, large-scale studies involving 10,000s of participants are still limited, especially in Asian populations. Therefore, we started the Tsuruoka Metabolomics Cohort Study enrolling 11,002 community-dwelling adults in Japan, and using capillary electrophoresis-mass spectrometry (CE-MS) and liquid chromatography–mass spectrometry. The CE-MS method is highly amenable to absolute quantification of polar metabolites, however, its reliability for large-scale measurement is unclear. The aim of this study is to examine reproducibility and validity of large-scale CE-MS measurements. In addition, the study presents absolute concentrations of polar metabolites in human plasma, which can be used in future as reference ranges in a Japanese population.

Methods

Metabolomic profiling of 8,413 fasting plasma samples were completed using CE-MS, and 94 polar metabolites were structurally identified and quantified. Quality control (QC) samples were injected every ten samples and assessed throughout the analysis. Inter- and intra-batch coefficients of variation of QC and participant samples, and technical intraclass correlation coefficients were estimated. Passing-Bablok regression of plasma concentrations by CE-MS on serum concentrations by standard clinical chemistry assays was conducted for creatinine and uric acid.

Results and conclusions

In QC samples, coefficient of variation was less than 20% for 64 metabolites, and less than 30% for 80 metabolites out of the 94 metabolites. Inter-batch coefficient of variation was less than 20% for 81 metabolites. Estimated technical intraclass correlation coefficient was above 0.75 for 67 metabolites. The slope of Passing-Bablok regression was estimated as 0.97 (95% confidence interval: 0.95, 0.98) for creatinine and 0.95 (0.92, 0.96) for uric acid. Compared to published data from other large cohort measurement platforms, reproducibility of metabolites common to the platforms was similar to or better than in the other studies. These results show that our CE-MS platform is suitable for conducting large-scale epidemiological studies.

The first high frequency of recombination-like conjugal transfer from an integrated origin of transfer sequence in Bacillus subtilis 168

Bacillus subtilis 168 was developed as a genome vector to manipulate large DNA fragments. The system is based on the inherent natural transformation (TF) activity. However, DNA size transferred by TF is limited up to approximately 100 kb. A conjugal transfer system capable of transferring DNA fragments considerably larger than those transferred by TF was developed. A well-defined oriT¹¹⁰ sequence and a cognate relaxase gene from the pUB110 plasmid were inserted into the xkdE gene of the B. subtilis genome. Transfer of antibiotic resistance markers distant from the oriT¹¹⁰ locus to the recipient B. subtilis occurred only in the presence of pLS20, a helper plasmid that provides a type IV secretion system. Marker transmission was consistent with the orientation of oriT¹¹⁰ and required a recA-proficient recipient. The first conjugal transfer system of genomic DNA should provide a valuable alternative genetic tool for editing the B. subtilis genome.

Systematic Analysis of the Binding Surfaces between tRNAs and Their Respective Aminoacyl tRNA Synthetase Based on Structural and Evolutionary Data

To determine the mechanism underlying the flow of genetic information, it is important to understand the relationship between a tRNA and its binding enzyme, a member of the aminoacyl-tRNA synthetase (aaRS) family. We have developed a novel method to project the interacting regions of tRNA-aaRS complexes, obtained from their three-dimensional structures, onto two-dimensional space. The interacting surface between each tRNA and its aaRS was successfully identified by determining these interactions with an atomic distance threshold of 3.3 Å. We analyzed their interactions, using 60 mainly bacterial and eukaryotic tRNA-aaRS complexes, and showed that the tRNA sequence regions that interacted most strongly with each aaRS are the anticodon loop and the CCA terminal region, followed by the D-stem. A sequence conservation analysis of the canonical tRNAs was conducted in 83 bacterial, 182 archaeal, and 150 eukaryotic species. Our results show that the three tRNA regions that interact with the aaRS and two additional loop regions (D-loop and TΨC-loop) known to be important for formation of the tRNA L-shaped structure are broadly conserved. We also found sequence conservations near the tRNA discriminator in the Bacteria and Archaea, and an enormous number of noncanonical tRNAs in the Eukaryotes. This is the first global view of tRNA evolution based on its structure and an unprecedented number of sequence data.

Comparative genomic analysis of Clostridium difficile ribotype 027 strains including the newly sequenced strain NCKUH-21 isolated from a patient in Taiwan

Background

Clostridium difficile is a Gram-positive anaerobe and the leading cause of antibiotic-associated diarrhea worldwide. The emergence of ribotype 027 (RT027) strains is associated with increased incidence of infection and mortality. To further understand the relationship between C. difficile NCKUH-21, a RT027 strain isolated from a patient in Taiwan, and other RT027 strains, we performed whole-genome shotgun sequencing on NCKUH-21 and comparative genomic analyses.

Results

The genome size, G+C content, and gene number for the NCKUH-21 strain were determined to be similar to those for other C. difficile strains. The core genome phylogeny indicated that the five RT027 strains R20291, CD196, NCKUH-21, BI1, and 2007855 formed a clade. A pathogenicity locus, tcdR-tcdB-tcdE-orf-tcdA-tcdC, was conserved in the genome. A genomic region highly similar to the Clostridium phage \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upvarphi$$\end{document}φCD38-2 was present in the NCKUH-21 strain but absent in the other RT027 strains and designated as the prophage \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upvarphi$$\end{document}φNCKUH-21. The prophage \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upvarphi$$\end{document}φNCKUH-21 genes were significantly higher in G+C content than the other genes in the NCKUH-21 genome, indicating that the prophage does not match the base composition of the host genome.

Conclusions

This is the first whole-genome analysis of a RT027 C. difficile strain isolated from Taiwan. Due to the high identity with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upvarphi$$\end{document}φCD38-2, the prophage identified in the NCKUH-21 genome has the potential to regulate toxin production. These results provide important information for understanding the pathogenicity of RT027 C. difficile in Taiwan.

Electronic supplementary material

The online version of this article (10.1186/s13099-017-0219-4) contains supplementary material, which is available to authorized users.

The GC Skew Index: A Measure of Genomic Compositional Asymmetry and the Degree of Replicational Selection

Circular bacterial chromosomes have highly polarized nucleotide composition in the two replichores, and this genomic strand asymmetry can be visualized using GC skew graphs. Here we propose and discuss the GC skew index (GCSI) for the quantification of genomic compositional skew, which combines a normalized measure of fast Fourier transform to capture the shape of the skew graph and Euclidean distance between the two vertices in a cumulative skew graph to represent the degree of skew. We calculated GCSI for all available bacterial genomes, and GCSI correlated well with the visibility of GC skew. This novel index is useful for estimating confidence levels for the prediction of replication origin and terminus by methods based on GC skew and for measuring the strength of replicational selection in a genome.

Selection Effects on the Positioning of Genes and Gene Structures from the Interplay of Replication and Transcription in Bacterial Genomes

Bacterial chromosomes are partly shaped by the functional requirements for efficient replication, which lead to strand bias as commonly characterized by the excess of guanines over cytosines in the leading strand. Gene structures are also highly organized within bacterial genomes as a result of such functional constraints, displaying characteristic positioning and structuring along the genome. Here we analyze the gene structures in completely sequenced bacterial chromosomes to observe the positional constraints on gene orientation, length, and codon usage with regard to the positions of replication origin and terminus. Selection on these gene features is different in regions surrounding the terminus of replication from the rest of the genome, but the selection could be either positive or negative depending on the species, and these positional effects are partly attributed to the A-T enrichment near the terminus. Characteristic gene structuring relative to the position of replication origin and terminus is commonly observed among most bacterial species with circular chromosomes, and therefore we argue that the highly organized gene positioning as well as the strand bias should be considered for genomics studies of bacteria.

Evolutionary Analysis of HIV-1 Pol Proteins Reveals Representative Residues for Viral Subtype Differentiation

RNA viruses have been used as model systems to understand the patterns and processes of molecular evolution because they have high mutation rates and are genetically diverse. Human immunodeficiency virus 1 (HIV-1), the etiological agent of acquired immune deficiency syndrome, is highly genetically diverse, and is classified into several groups and subtypes. However, it has been difficult to use its diverse sequences to establish the overall phylogenetic relationships of different strains or the trends in sequence conservation with the construction of phylogenetic trees. Our aims were to systematically characterize HIV-1 subtype evolution and to identify the regions responsible for HIV-1 subtype differentiation at the amino acid level in the Pol protein, which is often used to classify the HIV-1 subtypes. In this study, we systematically characterized the mutation sites in 2,052 Pol proteins from HIV-1 group M (144 subtype A; 1,528 subtype B; 380 subtype C), using sequence similarity networks. We also used spectral clustering to group the sequences based on the network graph structures. A stepwise analysis of the cluster hierarchies allowed us to estimate a possible evolutionary pathway for the Pol proteins. The subtype A sequences also clustered according to when and where the viruses were isolated, whereas both the subtype B and C sequences remained as single clusters. Because the Pol protein has several functional domains, we identified the regions that are discriminative by comparing the structures of the domain-based networks. Our results suggest that sequence changes in the RNase H domain and the reverse transcriptase (RT) connection domain are responsible for the subtype classification. By analyzing the different amino acid compositions at each site in both domain sequences, we found that a few specific amino acid residues (i.e., M357 in the RT connection domain and Q480, Y483, and L491 in the RNase H domain) represent the differences among the subtypes. These residues were located on the surface of the RT structure and in the vicinity of the amino acid sites responsible for RT enzymatic activity or function.

eRP arrangement: a strategy for assembled genomic contig rearrangement based on replication profiling in bacteria

Background

The reduced cost of sequencing has made de novo sequencing and the assembly of draft microbial genomes feasible in any ordinary biology lab. However, the process of finishing and completing the genome remains labor-intensive and computationally challenging in some cases, such as in the study of complete genome sequences, genomic rearrangements, long-range syntenic relationships, and structural variations.

Methods

Here, we show a contig reordering strategy based on experimental replication profiling (eRP) to recapitulate the bacterial genome structure within draft genomes. During the exponential growth phase, the majority of bacteria show a global genomic copy number gradient that is enriched near the replication origin and gradually declines toward the terminus. Therefore, if genome sequencing is performed with appropriate timing, the short-read coverage reflects this copy number gradient, providing information about the contig positions relative to the replication origin and terminus.

Results

We therefore investigated the appropriate timing for genomic DNA sampling and developed an algorithm for the reordering of the contigs based on eRP. As a result, this strategy successfully recapitulates the genomic structure of various structural mutants with draft genome sequencing.

Conclusions

Our strategy was successful for contig rearrangement with intracellular DNA replication behavior mechanisms and can be applied to almost all bacteria because the DNA replication system is highly conserved. Therefore, eRP makes it possible to understand genomic structural information and long-range syntenic relationships using a draft genome that is based on short reads.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4162-z) contains supplementary material, which is available to authorized users.