Menu
April 21, 2020

The Anaplasma ovis genome reveals a high proportion of pseudogenes.

The genus Anaplasma is made up of organisms characterized by small genomes that are undergoing reductive evolution. Anaplasma ovis, one of the seven recognized species in this genus, is an understudied pathogen of sheep and other ruminants. This tick-borne agent is thought to induce only mild clinical disease; however, small deficits may add to larger economic impacts due to the wide geographic distribution of this pathogen.In this report we present the first complete genome sequence for A. ovis and compare the genome features with other closely related species. The 1,214,674?bp A. ovis genome encodes 933 protein coding sequences, the split operon arrangement for ribosomal RNA genes, and more pseudogenes than previously recognized for other Anaplasma species. The metabolic potential is similar to other Anaplasma species. Anaplasma ovis has a small repertoire of surface proteins and transporters. Several novel genes are identified.Analyses of these important features and significant gene families/genes with potential to be vaccine candidates are presented in a comparative context. The availability of this genome will significantly facilitate research for this pathogen.


April 21, 2020

Structural variation of centromeric endogenous retroviruses in human populations and their impact on cutaneous T-cell lymphoma, Sézary syndrome, and HIV infection.

Human Endogenous Retroviruses type K HML-2 (HK2) are integrated into 117 or more areas of human chromosomal arms while two newly discovered HK2 proviruses, K111 and K222, spread extensively in pericentromeric regions, are the first retroviruses discovered in these areas of our genome.We use PCR and sequencing analysis to characterize pericentromeric K111 proviruses in DNA from individuals of diverse ethnicities and patients with different diseases.We found that the 5′ LTR-gag region of K111 proviruses is missing in certain individuals, creating pericentromeric instability. K111 deletion (-/- K111) is seen in about 15% of Caucasian, Asian, and Middle Eastern populations; it is missing in 2.36% of African individuals, suggesting that the -/- K111 genotype originated out of Africa. As we identified the -/-K111 genotype in Cutaneous T-cell lymphoma (CTCL) cell lines, we studied whether the -/-K111 genotype is associated with CTCL. We found a significant increase in the frequency of detection of the -/-K111 genotype in Caucasian patients with severe CTCL and/or Sézary syndrome (n?=?35, 37.14%), compared to healthy controls (n?=?160, 15.6%) [p?=?0.011]. The -/-K111 genotype was also found to vary in HIV-1 infection. Although Caucasian healthy individuals have a similar frequency of detection of the -/- K111 genotype, Caucasian HIV Long-Term Non-Progressors (LTNPs) and/or elite controllers, have significantly higher detection of the -/-K111 genotype (30.55%; n?=?36) than patients who rapidly progress to AIDS (8.5%; n?=?47) [p?=?0.0097].Our data indicate that pericentromeric instability is associated with more severe CTCL and/or Sézary syndrome in Caucasians, and appears to allow T-cells to survive lysis by HIV infection. These findings also provide new understanding of human evolution, as the -/-K111 genotype appears to have arisen out of Africa and is distributed unevenly throughout the world, possibly affecting the severity of HIV in different geographic areas.


April 21, 2020

Sequencing a Juglans regia?×?J. microcarpa hybrid yields high-quality genome assemblies of parental species.

Members of the genus Juglans are monecious wind-pollinated trees in the family Juglandaceae with highly heterozygous genomes, which greatly complicates genome sequence assembly. The genomes of interspecific hybrids are usually comprised of haploid genomes of parental species. We exploited this attribute of interspecific hybrids to avoid heterozygosity and sequenced an interspecific hybrid Juglans microcarpa?×?J. regia using a novel combination of single-molecule sequencing and optical genome mapping technologies. The resulting assemblies of both genomes were remarkably complete including chromosome termini and centromere regions. Chromosome termini consisted of arrays of telomeric repeats about 8?kb long and heterochromatic subtelomeric regions about 10?kb long. The centromeres consisted of arrays of a centromere-specific Gypsy retrotransposon and most contained genes, many of them transcribed. Juglans genomes evolved by a whole-genome-duplication dating back to the Cretaceous-Paleogene boundary and consist of two subgenomes, which were fractionated by numerous short gene deletions evenly distributed along the length of the chromosomes. Fractionation was shown to be asymmetric with one subgenome exhibiting greater gene loss than the other. The asymmetry of the process is ongoing and mirrors an asymmetry in gene expression between the subgenomes. Given the importance of J. microcarpa?×?J. regia hybrids as potential walnut rootstocks, we catalogued disease resistance genes in the parental genomes and studied their chromosomal distribution. We also estimated the molecular clock rates for woody perennials and deployed them in estimating divergence times of Juglans genomes and those of other woody perennials.


April 21, 2020

Identification of Candidate Genes for the Plateau Adaptation of a Tibetan Amphipod, Gammarus lacustris, Through Integration of Genome and Transcriptome Sequencing.

The amphipod Gammarus lacustris has been distributing in the Tibetan region with well-known uplifts of the Tibetan plateau. It is hence considered as a good model for investigating stress adaptations of the plateau. Here, we sequenced the whole-genome and full-length transcriptome of G. lacustris, and compared the transcriptome results with its counterpart Gammarus pisinnus from a nearby plain. Our main goal was to provide a genomic resource for investigation of genetic mechanisms, by which G. lacustris adapted to living on the plateau. The final draft genome assembly of G. lacustris was 5.07 gigabases (Gb), and it contained 443,304 scaffolds (>2 kb) with an N50 of 2,578 bp. A total of 8,858 unigenes were predicted in the full-length transcriptome of G. lacustris, with an average gene length of 1,811 bp. Compared with the G. pisinnus transcriptome, 2,672 differentially expressed genes (DEGs) were up-regulated and 2,881 DEGs were down-regulated in the G. lacustris transcriptome. Along with these critical DEGs, several enriched metabolic pathways, such as oxidative phosphorylation, ribosome, cell energy homeostasis, glycolysis and gluconeogenesis, were predicted to play essential roles in the plateau adaptation. In summary, the present study provides a genomic basis for understanding the plateau adaption of G. lacustris, which lays a fundamental basis for further biological and ecological studies on other resident aquatic species in the Tibetan plateau.


April 21, 2020

IMOS: improved Meta-aligner and Minimap2 On Spark.

Long reads provide valuable information regarding the sequence composition of genomes. Long reads are usually very noisy which renders their alignments on the reference genome a daunting task. It may take days to process datasets enough to sequence a human genome on a single node. Hence, it is of primary importance to have an aligner which can operate on distributed clusters of computers with high performance in accuracy and speed.In this paper, we presented IMOS, an aligner for mapping noisy long reads to the reference genome. It can be used on a single node as well as on distributed nodes. In its single-node mode, IMOS is an Improved version of Meta-aligner (IM) enhancing both its accuracy and speed. IM is up to 6x faster than the original Meta-aligner. It is also implemented to run IM and Minimap2 on Apache Spark for deploying on a cluster of nodes. Moreover, multi-node IMOS is faster than SparkBWA while executing both IM (1.5x) and Minimap2 (25x).In this paper, we purposed an architecture for mapping long reads to a reference. Due to its implementation, IMOS speed can increase almost linearly with respect to the number of nodes in a cluster. Also, it is a multi-platform application able to operate on Linux, Windows, and macOS.


April 21, 2020

Complete genome sequence of the halophilic PHA-producing bacterium Halomonas sp. SF2003: insights into its biotechnological potential.

A halophilic Gram-negative eubacterium was isolated from the Iroise Sea and identified as an efficient producer of polyhydroxyalkanoates (PHA). The strain, designated SF2003, was found to belong to the Halomonas genus on the basis of 16S rRNA gene sequence similarity. Previous biochemical tests indicated that the Halomonas sp. strain SF2003 is capable of supporting various culture conditions which sometimes can be constraining for marine strains. This versatility could be of great interest for biotechnological applications. Therefore, a complete bacterial genome sequencing and de novo assembly were performed using a PacBio RSII sequencer and Hierarchical Genome Assembly Process software in order to predict Halomonas sp. SF2003 metabolisms, and to identify genes involved in PHA production and stress tolerance. This study demonstrates the complete genome sequence of Halomonas sp. SF2003 which contains a circular 4,36 Mbp chromosome, and replaces the strain in a phylogenetic tree. Genes related to PHA metabolism, carbohydrate metabolism, fatty acid metabolism and stress tolerance were identified and a comparison was made with metabolisms of relative species. Genes annotation highlighted the presence of typical genes involved in PHA biosynthesis such as phaA, phaB and phaC and enabled a preliminary analysis of their organization and characteristics. Several genes of carbohydrates and fatty acid metabolisms were also identified which provided helpful insights into both a better knowledge of the intricacies of PHA biosynthetic pathways and of production purposes. Results show the strong versatility of Halomonas sp. SF2003 to adapt to various temperatures and salinity which can subsequently be exploited for industrial applications such as PHA production.


April 21, 2020

Replicon-Based Typing of IncI-Complex Plasmids, and Comparative Genomics Analysis of IncI?/K1 Plasmids.

IncI-complex plasmids can be divided into seven subgroups IncI1, IncI2, IncI?, IncB/O, IncK1, IncK2, and IncZ. In this study, a replicon-based scheme was proposed for typing IncI-complex plasmids into four separately clustering subgroups IncI2, IncI1/B/O, IncI?/K1 and IncK2/Z, the last three of which were combined from IncI1 and IncB/O, IncI? and IncK1, and IncK2 and IncZ, respectively. Four IncI?/K1 plasmids p205880-NR2, p14E509-CTXM, p11011-CTXM and p61806-CTXM were fully sequenced and compared with IncI?/K1 reference pCT, IncI2 reference R721, IncI1/B/O reference R64 and IncK2/Z reference pO26-CRL-125. These plasmids shared conserved gene organization in the replication and conjugal transfer regions, but displaying considerable sequence diversity among different subgroups. Remarkable modular differences were observed in the maintenance and transfer leading regions. p205880-NR2 contained no resistance genes or accessory modules, while the other seven plasmids acquired one or more accessory modules, which harbored mobile elements [including unit transposons, insertion sequence (IS)-based transposition units and individual IS elements] and associated resistance markers (especially including those involved in resistance to ß-lactams, aminoglycosides, tetracyclins, phenicols, streptomycins, trimethoprims, sulphonamides, tunicamycins and erythromycins). Data presented here provided a deeper insight into diversification and evolution of IncI-complex plasmids.


April 21, 2020

The complexity of the Fragaria x ananassa (octoploid) transcriptome by single-molecule long-read sequencing.

Strawberry (Fragaria x ananassa) is an allopolyploid species with diverse and complex transcripts. The regulatory mechanisms of fruit development and maturation have been extensively studied; however, little is known about the signaling mechanisms that direct this process in octoploid strawberry (Fragaria x ananassa). Here, we used long-read sequencing (LRS) technology and RNA-seq analysis to investigate the diversity and complexity of the polyploid transcriptome and differentially expressed transcripts along four successive fruit developmental stages of cultivated strawberry. We obtained a reference transcriptome with 119,897 unique full-length isoforms, including 2017 new isoforms and 2510 long noncoding RNAs. Based on the genome of the plausible progenitor (Fragaria vesca), 20,229 alternative splicing (AS) events were identified. Using this transcriptome, we found 17,485 differentially expressed transcripts during strawberry fruit development, including 527 transcription factors (TFs) belonging to 41 families. The expression profiles of all members of the auxin, ABA pathway, and anthocyanin biosynthesis gene families were also examined, and many of them were highly expressed at the ripe fruit stage, strongly indicating that the role of those genes is in the regulation of fruit ripening. We produce a high-quality reference transcriptome for octoploid strawberry, including much of the full-length transcript diversity, to help understand the regulatory mechanisms of fruit development and maturation of polyploid species, particularly via elucidation of the biochemical pathways involved in auxin, ABA, and anthocyanin biosynthesis.


April 21, 2020

Genome plasticity favours double chromosomal Tn4401b-blaKPC-2 transposon insertion in the Pseudomonas aeruginosa ST235 clone.

Pseudomonas aeruginosa Sequence Type 235 is a clone that possesses an extraordinary ability to acquire mobile genetic elements and has been associated with the spread of resistance genes, including genes that encode for carbapenemases. Here, we aim to characterize the genetic platforms involved in resistance dissemination in blaKPC-2-positive P. aeruginosa ST235 in Colombia.In a prospective surveillance study of infections in adult patients attended in five ICUs in five distant cities in Colombia, 58 isolates of P. aeruginosa were recovered, of which, 27 (46.6%) were resistant to carbapenems. The molecular analysis showed that 6 (22.2%) and 4 (14.8%) isolates harboured the blaVIM and blaKPC-2 genes, respectively. The four blaKPC-2-positive isolates showed a similar PFGE pulsotype and belonged to ST235. Complete genome sequencing of a representative ST235 isolate shows a unique chromosomal contig of 7097.241?bp with eight different resistance genes identified and five transposons: a Tn6162-like with ant(2?)-Ia, two Tn402-like with ant(3?)-Ia and blaOXA-2 and two Tn4401b with blaKPC-2. All transposons were inserted into the genomic islands. Interestingly, the two Tn4401b copies harbouring blaKPC-2 were adjacently inserted into a new genomic island (PAGI-17) with traces of a replicative transposition process. This double insertion was probably driven by several structural changes within the chromosomal region containing PAGI-17 in the ST235 background.This is the first report of a double Tn4401b chromosomal insertion in P. aeruginosa, just within a new genomic island (PAGI-17). This finding indicates once again the great genomic plasticity of this microorganism.


April 21, 2020

Complete genome sequence of Caulobacter flavus RHGG3T, a type species of the genus Caulobacter with plant growth-promoting traits and heavy metal resistance.

Caulobacter flavus RHGG3T, a novel type species in the genus Caulobacter, originally isolated from rhizosphere soil of watermelon (Citrullus lanatus), has the ability to improve the growth of watermelon seedling and tolerate heavy metals. In vitro, C. flavus RHGG3T was able to solubilize phosphate (80.56 mg L-1), produce indole-3-acetic acid (IAA) (11.58 mg L-1) and was resistant to multiple heavy metals (copper, zinc, cadmium, cobalt and lead). Inoculating watermelon with this strain increased shoot and root length by 22.1% and 43.7%, respectively, and the total number of lateral roots by 55.9% compared to non-inoculated watermelon. In this study, we present the complete genome sequence of C. flavus RHGG3T, which was comprised of a single circular chromosome of 5,659,202 bp with a G?+?C content of 69.25%. An annotation analysis revealed that the C. flavus RHGG3T genome contained 5172 coding DNA sequences, 9 rRNA and 55 tRNA genes. Genes related to plant growth promotion (PGP), such as those associated with phosphate solubilization, nitrogen fixation, IAA, phenazine, volatile compounds, spermidine and cobalamin synthesis, were found in the C. flavus RHGG3T genome. Some genes responsible for heavy metal tolerance were also identified. The genome sequence of strain RHGG3T reported here provides new insight into the molecular mechanisms underlying the promotion of plant growth and the resistance to heavy metals in C. flavus. This study will be valuable for further exploration of the biotechnological applications of strain RHGG3T in agriculture.


April 21, 2020

Genomic characterisation of the new Dickeya fangzhongdai species regrouping plant pathogens and environmental isolates.

The Dickeya genus is part of the Pectobacteriaceae family that is included in the newly described enterobacterales order. It comprises a group of aggressive soft rot pathogens with wide geographic distribution and host range. Among them, the new Dickeya fangzhongdai species groups causative agents of maceration-associated diseases that impact a wide variety of crops and ornamentals. It affects mainly monocot plants, but D. fangzhongdai strains have also been isolated from pear trees and water sources. Here, we analysed which genetic novelty exists in this new species, what are the D. fangzhongdai-specific traits and what is the intra-specific diversity.The genomes of eight D. fangzhongdai strains isolated from diverse environments were compared to 31 genomes of strains belonging to other Dickeya species. The D. fangzhongdai core genome regroups approximately 3500 common genes, including most genes that encode virulence factors and regulators characterised in the D. dadantii 3937 model strain. Only 38 genes are present in D. fangzhongdai and absent in all other Dickeyas. One of them encodes a pectate lyase of the PL10 family of polysaccharide lyases that is found only in a few bacteria from the plant environment, soil or human gut. Other D. fangzhongdai-specific genes with a known or predicted function are involved in regulation or metabolism. The intra-species diversity analysis revealed that seven of the studied D. fangzhongdai strains were grouped into two distinct clades. Each clade possesses a pool of 100-150 genes that are shared by the clade members, but absent from the other D. fangzhongdai strains and several of these genes are clustered into genomic regions. At the strain level, diversity resides mainly in the arsenal of T5SS- and T6SS-related toxin-antitoxin systems and in secondary metabolite biogenesis pathways.This study identified the genome-specific traits of the new D. fangzhongdai species and highlighted the intra-species diversity of this species. This diversity encompasses secondary metabolites biosynthetic pathways and toxins or the repertoire of genes of extrachromosomal origin. We however didn’t find any relationship between gene content and phenotypic differences or sharing of environmental habitats.


April 21, 2020

Whole Genome Sequencing and Comparative Genomics Analyses of Pandoraea sp. XY-2, a New Species Capable of Biodegrade Tetracycline.

Few bacteria are resistant to tetracycline and can even biodegrade tetracycline in the environment. In this study, we isolated a bacterium Pandoraea sp. XY-2, which could biodegrade 74% tetracycline at pH 7.0 and 30°C within 6 days. Thereafter, we determined the whole genome sequence of Pandoraea sp. XY-2 genome is a single circular chromosome of 5.06 Mb in size. Genomic annotation showed that two AA6 family members-encoding genes and nine glutathione S-transferase (GSTs)-encoding genes could be relevant to tetracycline biodegradation. In addition, the average nucleotide identities (ANI) analysis between the genomes of Pandoraea sp. XY-2 and other Pandoraea spp. revealed that Pandoraea sp. XY-2 belongs to a new species. Moreover, comparative genome analysis of 36 Pandoraea strains identified the pan and specific genes, numerous single nucleotide polymorphisms (SNPs), insertions, and deletion variations (InDels) and different syntenial relationships in the genome of Pandoraea sp. XY-2. Finally, the evolution and the origin analysis of genes related to tetracycline resistance revealed that the six tetA(48) genes and two specificgenes tetG and tetR in Pandoraea sp. XY-2 were acquired by horizontal gene transfer (HGT) events from sources related to Paraburkholderia, Burkholderia, Caballeronia, Salmonella, Vibrio, Proteobacteria, Pseudomonas, Acinetobacter, Flavimaricola, and some unidentified sources. As a new species, Pandoraea sp. XY-2 will be an excellent resource for the bioremediation of tetracycline-contaminated environment.


April 21, 2020

Direct pathway cloning of the sodorifen biosynthetic gene cluster and recombinant generation of its product in E. coli.

Serratia plymuthica WS3236 was selected for whole genome sequencing based on preliminary genetic and chemical screening indicating the presence of multiple natural product pathways. This led to the identification of a putative sodorifen biosynthetic gene cluster (BGC). The natural product sodorifen is a volatile organic compound (VOC) with an unusual polymethylated hydrocarbon bicyclic structure (C16H26) produced by selected strains of S. plymuthica. The BGC encoding sodorifen consists of four genes, two of which (sodA, sodB) are homologs of genes encoding enzymes of the non-mevalonate pathway and are thought to enhance the amounts of available farnesyl pyrophosphate (FPP), the precursor of sodorifen. Proceeding from FPP, only two enzymes are necessary to produce sodorifen: an S-adenosyl methionine dependent methyltransferase (SodC) with additional cyclisation activity and a terpene-cyclase (SodD). Previous analysis of S. plymuthica found sodorifen production titers are generally low and vary significantly among different producer strains. This precludes studies on the still elusive biological function of this structurally and biosynthetically fascinating bacterial terpene.Sequencing and mining of the S. plymuthica WS3236 genome revealed the presence of 38 BGCs according to antiSMASH analysis, including a putative sodorifen BGC. Further genome mining for sodorifen and sodorifen-like BGCs throughout bacteria was performed using SodC and SodD as queries and identified a total of 28 sod-like gene clusters. Using direct pathway cloning (DiPaC) we intercepted the 4.6 kb candidate sodorifen BGC from S. plymuthica WS3236 (sodA-D) and transformed it into Escherichia coli BL21. Heterologous expression under the control of the tetracycline inducible PtetO promoter firmly linked this BGC to sodorifen production. By utilizing this newly established expression system, we increased the production yields by approximately 26-fold when compared to the native producer. In addition, sodorifen was easily isolated in high purity by simple head-space sampling.Genome mining of all available genomes within the NCBI and JGI IMG databases led to the identification of a wealth of sod-like pathways which may be responsible for producing a range of structurally unknown sodorifen analogs. Introduction of the S. plymuthica WS3236 sodorifen BGC into the fast-growing heterologous expression host E. coli with a very low VOC background led to a significant increase in both sodorifen product yield and purity compared to the native producer. By providing a reliable, high-level production system, this study sets the stage for future investigations of the biological role and function of sodorifen and for functionally unlocking the bioinformatically identified putative sod-like pathways.


April 21, 2020

Full-length transcript sequencing and comparative transcriptomic analysis to evaluate the contribution of osmotic and ionic stress components towards salinity tolerance in the roots of cultivated alfalfa (Medicago sativa L.).

Alfalfa is the most extensively cultivated forage legume. Salinity is a major environmental factor that impacts on alfalfa’s productivity. However, little is known about the molecular mechanisms underlying alfalfa responses to salinity, especially the relative contribution of the two important components of osmotic and ionic stress.In this study, we constructed the first full-length transcriptome database for alfalfa root tips under continuous NaCl and mannitol treatments for 1, 3, 6, 12, and 24?h (three biological replicates for each time points, including the control group) via PacBio Iso-Seq. This resulted in the identification of 52,787 full-length transcripts, with an average length of 2551?bp. Global transcriptional changes in the same 33 stressed samples were then analyzed via BGISEQ-500 RNA-Seq. Totals of 8861 NaCl-regulated and 8016 mannitol-regulated differentially expressed genes (DEGs) were identified. Metabolic analyses revealed that these DEGs overlapped or diverged in the cascades of molecular networks involved in signal perception, signal transduction, transcriptional regulation, and antioxidative defense. Notably, several well characterized signalling pathways, such as CDPK, MAPK, CIPK, and PYL-PP2C-SnRK2, were shown to be involved in osmotic stress, while the SOS core pathway was activated by ionic stress. Moreover, the physiological shifts of catalase and peroxidase activity, glutathione and proline content were in accordance with dynamic transcript profiles of the relevant genes, indicating that antioxidative defense system plays critical roles in response to salinity stress.Overall, our study provides evidence that the response to salinity stress in alfalfa includes both osmotic and ionic components. The key osmotic and ionic stress-related genes are candidates for future studies as potential targets to improve resistance to salinity stress via genetic engineering.


April 21, 2020

Construction of JRG (Japanese reference genome) with single-molecule real-time sequencing

In recent genome analyses, population-specific reference panels have indicated important. However, reference panels based on short-read sequencing data do not sufficiently cover long insertions. Therefore, the nature of long insertions has not been well documented. Here, we assembled a Japanese genome using single-molecule real-time sequencing data and characterized insertions found in the assembled genome. We identified 3691 insertions ranging from 100?bps to ~10,000?bps in the assembled genome relative to the international reference sequence (GRCh38). To validate and characterize these insertions, we mapped short-reads from 1070 Japanese individuals and 728 individuals from eight other populations to insertions integrated into GRCh38. With this result, we constructed JRGv1 (Japanese Reference Genome version 1) by integrating the 903 verified insertions, totaling 1,086,173 bases, shared by at least two Japanese individuals into GRCh38. We also constructed decoyJRGv1 by concatenating 3559 verified insertions, totaling 2,536,870 bases, shared by at least two Japanese individuals or by six other assemblies. This assembly improved the alignment ratio by 0.4% on average. These results demonstrate the importance of refining the reference assembly and creating a population-specific reference genome. JRGv1 and decoyJRGv1 are available at the JRG website.


Talk with an expert

If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.