Menu
September 22, 2019  |  

Noise-Cancelling Repeat Finder: Uncovering tandem repeats in error-prone long-read sequencing data

Tandem DNA repeats can be sequenced with long-read technologies, but cannot be accurately deciphered due to the lack of computational tools taking high error rates of these technologies into account. Here we introduce Noise-Cancelling Repeat Finder (NCRF) to uncover putative tandem repeats of specified motifs in noisy long reads produced by Pacific Biosciences and Oxford Nanopore sequencers. Using simulations, we validated the use of NCRF to locate tandem repeats with motifs of various lengths and demonstrated its superior performance as compared to two alternative tools. Using real human whole-genome sequencing data, NCRF identified long arrays of the (AATGG)n repeat involved in heat shock stress response.


September 22, 2019  |  

Reconstitution of eukaryotic chromosomes and manipulation of DNA N6-methyladenine alters chromatin and gene expression

DNA N6-adenine methylation (6mA) has recently been reported in diverse eukaryotes, spanning unicellular organisms to metazoans. Yet the functional significance of 6mA remains elusive due to its low abundance, difficulty of manipulation within native DNA, and lack of understanding of eukaryotic 6mA writers. Here, we report a novel DNA 6mA methyltransferase in ciliates, termed MTA1. The enzyme contains an MT-A70 domain but is phylogenetically distinct from all known RNA and DNA methyltransferases. Disruption of MTA1 in vivo leads to the genome-wide loss of 6mA in asexually growing cells and abolishment of the consensus ApT dimethylated motif. Genes exhibit subtle changes in chromatin organization or RNA expression upon loss of 6mA, depending on their starting methylation level. Mutants fail to complete the sexual cycle, which normally coincides with a peak of MTA1 expression. Thus, MTA1 functions in a developmental stage-specific manner. We determine the impact of 6mA on chromatin organization in vitro by reconstructing complete, full-length ciliate chromosomes harboring 6mA in native or ectopic positions. Using these synthetic chromosomes, we show that 6mA directly disfavors nucleosomes in vitro in a local, quantitative manner, independent of DNA sequence. Furthermore, the chromatin remodeler ACF can overcome this effect. Our study identifies a novel MT-A70 protein necessary for eukaryotic 6mA methylation and defines the impact of 6mA on chromatin organization using epigenetically defined synthetic chromosomes.


September 22, 2019  |  

Genomic Tandem Quadruplication is Associated with Ketoconazole Resistance in Malassezia pachydermatis.

Malassezia pachydermatis is a commensal yeast found on the skin of dogs. However, M. pachydermatis is also considered an opportunistic pathogen and is associated with various canine skin diseases including otitis externa and atopic dermatitis, which usually require treatment using an azole antifungal drug, such as ketoconazole. In this study, we isolated a ketoconazole-resistant strain of M. pachydermatis, designated “KCTC 27587,” from the external ear canal of a dog with otitis externa and analyzed its resistance mechanism. To understand the mechanism underlying ketoconazole resistance of the clinical isolate M. pachydermatis KCTC 27587, the whole genome of the yeast was sequenced using the PacBio platform and was compared with M. pachydermatis type strain CBS 1879. We found that a ~84-kb region in chromosome 4 of M. pachydermatis KCTC 27587 was tandemly quadruplicated. The quadruplicated region contains 52 protein coding genes, including the homologs of ERG4 and ERG11, whose overexpression is known to be associated with azole resistance. Our data suggest that the quadruplication of the ~84-kb region may be the cause of the ketoconazole resistance in M. pachydermatis KCTC 27587.


September 22, 2019  |  

Complete genome sequence of Leuconostoc citreum EFEL2700, a host strain for transformation of pCB vectors.

Leuconostoc citreum is an important lactic acid bacterium used as a starter culture for producing kimchi, the traditional Korean fermented vegetables. An efficient host strain for plasmid transformation, L. citreum EFEL2700, was isolated from kimchi, and it has been frequently used for genetic engineering of L. citreum. In this study, we report the whole genome sequence of the strain and its genetic characteristics. Genome assembly yielded 5 contigs (1 chromosome and 4 plasmids), and the complete genome contained 1,923,830 base pairs (bp) with a G?+?C content of 39.0%. Average nucleotide identity analysis showed high homology (= 99%) to the reference strain L. citreum KM 20. The smallest plasmid (4.3 kbp) was used as an Escherichia coli shuttle vector (pCB) for heterologous gene expression, and L. citreum EFEL2700 showed the highest transformation efficiency, 6.7?×?104 CFU µg-1 DNA. Genetic analysis of the genome enabled the construction of primary metabolic pathway showing a typical hetero-type lactic acid fermentation. Notably, no core genes for primary metabolism were observed in plasmid 4 and it could be eliminated to create an efficient host for gene transformation. This report will facilitate the understanding and application of L. citreum EFEL2700 as a food-grade microbial cell factory.Copyright © 2018. Published by Elsevier B.V.


September 22, 2019  |  

Implications of stx loss for clinical diagnostics of Shiga toxin-producing Escherichia coli.

The dynamics related to the loss of stx genes from Shiga toxin-producing Escherichia coli remain unclear. Current diagnostic procedures have shortcomings in the detection and identification of STEC. This is partly owing to the fact that stx genes may be lost during an infection or in the laboratory. The aim of the present study was to provide new insight into in vivo and in vitro stx loss in order to improve diagnostic procedures. Results from the study support the theory that loss of stx is a strain-related phenomenon and not induced by patient factors. It was observed that one strain could lose stx both in vivo and in vitro. Whole genome comparison of stx-positive and stx-negative isolates from the same patient revealed that different genomic rearrangements, such as complete or partial loss of the parent prophage, may be factors in the loss of stx. Of diagnostic interest, it was shown that patients can be co-infected with different E. coli pathotypes. Therefore, identification of eae-positive, but stx-negative isolates should not be interpreted as “Shiga toxin-lost” E. coli without further testing. Growth and recovery of STEC were supported by different selective agar media for different strains, arguing for inclusion of several media in STEC diagnostics.


September 22, 2019  |  

Improved nucleic acid extraction protocols for Ganoderma boninense, G. miniatocinctum and G. tornatum.

The first and most crucial step of all molecular techniques is to isolate high quality and intact nucleic acids. However, DNA and RNA isolation from fungal samples are usually difficult due to the cell walls that are relatively unsusceptible to lysis and often resistant to traditional extraction procedures. Although there are many extraction protocols for Ganoderma species, different extraction protocols have been applied to different species to obtain high yields of good quality nucleic acids, especially for genome and transcriptome sequencing. Ganoderma species, mainly G. boninense causes the basal stem rot disease, a devastating disease that plagues the oil palm industry. Here, we describe modified DNA extraction protocols for G. boninense, G. miniatocinctum and G. tornatum, and an RNA extraction protocol for G. boninense. The modified salting out DNA extraction protocol is suitable for G. boninense and G. miniatocinctum while the modified high salt and low pH protocol is suitable for G. tornatum. The modified DNA and RNA extraction protocols were able to produce high quality genomic DNA and total RNA of?~?140 to 160 µg/g and?~?80 µg/g of mycelia respectively, for Single Molecule Real Time (PacBio Sequel® System) and Illumina sequencing. These protocols will benefit those studying the oil palm pathogens at nucleotide level.


September 22, 2019  |  

The genome of the tegu lizard Salvator merianae: combining Illumina, PacBio, and optical mapping data to generate a highly contiguous assembly.

Reptiles are a species-rich group with great phenotypic and life history diversity but are highly underrepresented among the vertebrate species with sequenced genomes.Here, we report a high-quality genome assembly of the tegu lizard, Salvator merianae, the first lacertoid with a sequenced genome. We combined 74X Illumina short-read, 29.8X Pacific Biosciences long-read, and optical mapping data to generate a high-quality assembly with a scaffold N50 value of 55.4 Mb. The contig N50 value of this assembly is 521 Kb, making it the most contiguous reptile assembly so far. We show that the tegu assembly has the highest completeness of coding genes and conserved non-exonic elements (CNEs) compared to other reptiles. Furthermore, the tegu assembly has the highest number of evolutionarily conserved CNE pairs, corroborating a high assembly contiguity in intergenic regions. As in other reptiles, long interspersed nuclear elements comprise the most abundant transposon class. We used transcriptomic data, homology- and de novo gene predictions to annotate 22,413 coding genes, of which 16,995 (76%) likely have human orthologs as inferred by CESAR-derived gene mappings. Finally, we generated a multiple genome alignment comprising 10 squamates and 7 other amniote species and identified conserved regions that are under evolutionary constraint. CNEs cover 38 Mb (1.8%) of the tegu genome, with 3.3 Mb in these elements being squamate specific. In contrast to placental mammal-specific CNEs, very few of these squamate-specific CNEs (<20 Kb) overlap transposons, highlighting a difference in how lineage-specific CNEs originated in these two clades.The tegu lizard genome together with the multiple genome alignment and comprehensive conserved element datasets provide a valuable resource for comparative genomic studies of reptiles and other amniotes.


September 22, 2019  |  

Regulation of yeast-to-hyphae transition in Yarrowia lipolytica.

The yeast Yarrowia lipolytica undergoes a morphological transition from yeast-to-hyphal growth in response to environmental conditions. A forward genetic screen was used to identify mutants that reliably remain in the yeast phase, which were then assessed by whole-genome sequencing. All the smooth mutants identified, so named because of their colony morphology, exhibit independent loss of DNA at a repetitive locus made up of interspersed ribosomal DNA and short 10- to 40-mer telomere-like repeats. The loss of repetitive DNA is associated with downregulation of genes with stress response elements (5′-CCCCT-3′) and upregulation of genes with cell cycle box (5′-ACGCG-3′) motifs in their promoter region. The stress response element is bound by the transcription factor Msn2p in Saccharomyces cerevisiae We confirmed that the Y. lipolyticamsn2 (Ylmsn2) ortholog is required for hyphal growth and found that overexpression of Ylmsn2 enables hyphal growth in smooth strains. The cell cycle box is bound by the Mbp1p/Swi6p complex in S. cerevisiae to regulate G1-to-S phase progression. We found that overexpression of either the Ylmbp1 or Ylswi6 homologs decreased hyphal growth and that deletion of either Ylmbp1 or Ylswi6 promotes hyphal growth in smooth strains. A second forward genetic screen for reversion to hyphal growth was performed with the smooth-33 mutant to identify additional genetic factors regulating hyphal growth in Y. lipolytica Thirteen of the mutants sequenced from this screen had coding mutations in five kinases, including the histidine kinases Ylchk1 and Ylnik1 and kinases of the high-osmolarity glycerol response (HOG) mitogen-activated protein (MAP) kinase cascade Ylssk2, Ylpbs2, and Ylhog1 Together, these results demonstrate that Y. lipolytica transitions to hyphal growth in response to stress through multiple signaling pathways.IMPORTANCE Many yeasts undergo a morphological transition from yeast-to-hyphal growth in response to environmental conditions. We used forward and reverse genetic techniques to identify genes regulating this transition in Yarrowia lipolytica We confirmed that the transcription factor Ylmsn2 is required for the transition to hyphal growth and found that signaling by the histidine kinases Ylchk1 and Ylnik1 as well as the MAP kinases of the HOG pathway (Ylssk2, Ylpbs2, and Ylhog1) regulates the transition to hyphal growth. These results suggest that Y. lipolytica transitions to hyphal growth in response to stress through multiple kinase pathways. Intriguingly, we found that a repetitive portion of the genome containing telomere-like and rDNA repeats may be involved in the transition to hyphal growth, suggesting a link between this region and the general stress response. Copyright © 2018 Pomraning et al.


September 22, 2019  |  

Genomic and transcriptomic comparisons of closely related malaria parasites differing in virulence and sequestration pattern.

Background: Malaria parasite species differ greatly in the harm they do to humans. While P. falciparum kills hundreds of thousands per year, P. vivax kills much less often and P. malariae is relatively benign. Strains of the rodent malaria parasite Plasmodium chabaudi show phenotypic variation in virulence during infections of laboratory mice. This make it an excellent species to study genes which may be responsible for this trait. By understanding the mechanisms which underlie differences in virulence we can learn how parasites adapt to their hosts and how we might prevent disease. Methods: Here we present a complete reference genome sequence for a more virulent P. chabaudi strain, PcCB, and perform a detailed comparison with the genome of the less virulent PcAS strain. Results: We found the greatest variation in the subtelomeric regions, in particular amongst the sequences of the pir gene family, which has been associated with virulence and establishment of chronic infection. Despite substantial variation at the sequence level, the repertoire of these genes has been largely maintained, highlighting the requirement for functional conservation as well as diversification in host-parasite interactions. However, a subset of pir genes, previously associated with increased virulence, were more highly expressed in PcCB, suggesting a role for this gene family in virulence differences between strains. We found that core genes involved in red blood cell invasion have been under positive selection and that the more virulent strain has a greater preference for reticulocytes, which has elsewhere been associated with increased virulence. Conclusions: These results provide the basis for a mechanistic understanding of the phenotypic differences between Plasmodium chabaudi strains, which might ultimately be translated into a better understanding of malaria parasites affecting humans.


September 22, 2019  |  

Detection and visualization of complex structural variants from long reads.

With applications in cancer, drug metabolism, and disease etiology, understanding structural variation in the human genome is critical in advancing the thrusts of individualized medicine. However, structural variants (SVs) remain challenging to detect with high sensitivity using short read sequencing technologies. This problem is exacerbated when considering complex SVs comprised of multiple overlapping or nested rearrangements. Longer reads, such as those from Pacific Biosciences platforms, often span multiple breakpoints of such events, and thus provide a way to unravel small-scale complexities in SVs with higher confidence.We present CORGi (COmplex Rearrangement detection with Graph-search), a method for the detection and visualization of complex local genomic rearrangements. This method leverages the ability of long reads to span multiple breakpoints to untangle SVs that appear very complicated with respect to a reference genome. We validated our approach against both simulated long reads, and real data from two long read sequencing technologies. We demonstrate the ability of our method to identify breakpoints inserted in synthetic data with high accuracy, and the ability to detect and plot SVs from NA12878 germline, achieving 88.4% concordance between the two sets of sequence data. The patterns of complexity we find in many NA12878 SVs match known mechanisms associated with DNA replication and structural variant formation, and highlight the ability of our method to automatically label complex SVs with an intuitive combination of adjacent or overlapping reference transformations.CORGi is a method for interrogating genomic regions suspected to contain local rearrangements using long reads. Using pairwise alignments and graph search CORGi produces labels and visualizations for local SVs of arbitrary complexity.


September 22, 2019  |  

Emerging multidrug-resistant hybrid pathotype shiga toxin-producing Escherichia coli O80 and related strains of clonal complex 165, Europe.

Enterohemorrhagic Escherichia coli serogroup O80, involved in hemolytic uremic syndrome associated with extraintestinal infections, has emerged in France. We obtained circularized sequences of the O80 strain RDEx444, responsible for hemolytic uremic syndrome with bacteremia, and noncircularized sequences of 35 O80 E. coli isolated from humans and animals in Europe with or without Shiga toxin genes. RDEx444 harbored a mosaic plasmid, pR444_A, combining extraintestinal virulence determinants and a multidrug resistance-encoding island. All strains belonged to clonal complex 165, which is distantly related to other major enterohemorrhagic E. coli lineages. All stx-positive strains contained eae-?, ehxA, and genes characteristic of pR444_A. Among stx-negative strains, 1 produced extended-spectrum ß-lactamase, 1 harbored the colistin-resistance gene mcr1, and 2 possessed genes characteristic of enteropathogenic and pyelonephritis E. coli. Because O80-clonal complex 165 strains can integrate intestinal and extraintestinal virulence factors in combination with diverse drug-resistance genes, they constitute dangerous and versatile multidrug-resistant pathogens.


September 22, 2019  |  

Genome of the small hive beetle (Aethina tumida, Coleoptera: Nitidulidae), a worldwide parasite of social bee colonies, provides insights into detoxification and herbivory.

The small hive beetle (Aethina tumida; ATUMI) is an invasive parasite of bee colonies. ATUMI feeds on both fruits and bee nest products, facilitating its spread and increasing its impact on honey bees and other pollinators. We have sequenced and annotated the ATUMI genome, providing the first genomic resources for this species and for the Nitidulidae, a beetle family that is closely related to the extraordinarily species-rich clade of beetles known as the Phytophaga. ATUMI thus provides a contrasting view as a neighbor for one of the most successful known animal groups.We present a robust genome assembly and a gene set possessing 97.5% of the core proteins known from the holometabolous insects. The ATUMI genome encodes fewer enzymes for plant digestion than the genomes of wood-feeding beetles but nonetheless shows signs of broad metabolic plasticity. Gustatory receptors are few in number compared to other beetles, especially receptors with known sensitivity (in other beetles) to bitter substances. In contrast, several gene families implicated in detoxification of insecticides and adaptation to diverse dietary resources show increased copy numbers. The presence and diversity of homologs involved in detoxification differ substantially from the bee hosts of ATUMI.Our results provide new insights into the genomic basis for local adaption and invasiveness in ATUMI and a blueprint for control strategies that target this pest without harming their honey bee hosts. A minimal set of gustatory receptors is consistent with the observation that, once a host colony is invaded, food resources are predictable. Unique detoxification pathways and pathway members can help identify which treatments might control this species even in the presence of honey bees, which are notoriously sensitive to pesticides.


September 22, 2019  |  

Approaches for surveying cosmic radiation damage in large populations of Arabidopsis thaliana seeds-Antarctic balloons and particle beams.

The Cosmic Ray Exposure Sequencing Science (CRESS) payload system is a proof of concept experiment to assess the genomic impact of space radiation on seeds. CRESS was designed as a secondary payload for the December 2016 high-altitude, high-latitude, and long-duration balloon flight carrying the Boron And Carbon Cosmic Rays in the Upper Stratosphere (BACCUS) experimental hardware. Investigation of the biological effects of Galactic Cosmic Radiation (GCR), particularly those of ions with High-Z and Energy (HZE), is of interest due to the genomic damage this type of radiation inflicts. The biological effects of upper-stratospheric mixed radiation above Antarctica (ANT) were sampled using Arabidopsis thaliana seeds and were compared to those resulting from a controlled simulation of GCR at Brookhaven National Laboratory (BNL) and to laboratory control seed. The payload developed for Antarctica exposure was broadly designed to 1U CubeSat specifications (10cmx10cmx10cm, =1.33kg), maintained 1 atm internal pressure, and carried an internal cargo of four seed trays (about 580,000 seeds) and twelve CR-39 Solid-State Nuclear Track Detectors (SSNTDs). The irradiated seeds were recovered, sterilized and grown on Petri plates for phenotypic screening. BNL and ANT M0 seeds showed significantly reduced germination rates and elevated somatic mutation rates when compared to non-irradiated controls, with the BNL mutation rate also being significantly higher than that of ANT. Genomic DNA from mutants of interest was evaluated with whole-genome sequencing using PacBio SMRT technology. Sequence data revealed the presence of an array of genome structural variants in the genomes of M0 and M1 mutant plants.


September 22, 2019  |  

First draft genome for red sea bream of family Sparidae.

Reference genomes for all organisms on earth are now attainable owing to advances in genome sequencing technologies (Goodwin et al., 2016). Generally, species that contribute considerably to the economy or human welfare are sequenced and are considered more important than others. Furthermore, coastal indigenous people mainly depend on marine species for their food sources, which has resulted in the extinction of several marine species (Cisneros-Montemayor et al., 2016). Of these, an extinction risk assessment of marine fishes, mainly for sea breams (Family: Sparidae), has recently been conducted by way of a global extinction risk assessment from the dataset of the International Union for Conservation of Nature’s Red List Process, which mentions that around 25 species are threatened/near-threatened according to their body weight (Comeros-Raynal et al., 2016). Another report clearly showed the benefit of worldwide aquaculture production, which contributed to 47% of total seafood production, and also highlighted the over-fishing of sea breams (FAO, 2018). The Republic of Korea is the fourth largest seafood producer in the world, producing 3.3 million tons in 2015 and exporting seafood worth $1.6 billion in 2016; therefore, aquaculture- associated research is fundamental for Korea. In the present study, the red sea bream (Pagrus major), which belongs to the family Sparidae, which comprises 35 genera, 132 species, and 10 subspecies (de la Herran et al., 2001; NCBI, 2018), was assessed.


September 22, 2019  |  

Chemical Synergy between Ionophore PBT2 and Zinc Reverses Antibiotic Resistance.

The World Health Organization reports that antibiotic-resistant pathogens represent an imminent global health disaster for the 21st century. Gram-positive superbugs threaten to breach last-line antibiotic treatment, and the pharmaceutical industry antibiotic development pipeline is waning. Here we report the synergy between ionophore-induced physiological stress in Gram-positive bacteria and antibiotic treatment. PBT2 is a safe-for-human-use zinc ionophore that has progressed to phase 2 clinical trials for Alzheimer’s and Huntington’s disease treatment. In combination with zinc, PBT2 exhibits antibacterial activity and disrupts cellular homeostasis in erythromycin-resistant group A Streptococcus (GAS), methicillin-resistant Staphylococcus aureus (MRSA), and vancomycin-resistant Enterococcus (VRE). We were unable to select for mutants resistant to PBT2-zinc treatment. While ineffective alone against resistant bacteria, several clinically relevant antibiotics act synergistically with PBT2-zinc to enhance killing of these Gram-positive pathogens. These data represent a new paradigm whereby disruption of bacterial metal homeostasis reverses antibiotic-resistant phenotypes in a number of priority human bacterial pathogens.IMPORTANCE The rise of bacterial antibiotic resistance coupled with a reduction in new antibiotic development has placed significant burdens on global health care. Resistant bacterial pathogens such as methicillin-resistant Staphylococcus aureus and vancomycin-resistant Enterococcus are leading causes of community- and hospital-acquired infection and present a significant clinical challenge. These pathogens have acquired resistance to broad classes of antimicrobials. Furthermore, Streptococcus pyogenes, a significant disease agent among Indigenous Australians, has now acquired resistance to several antibiotic classes. With a rise in antibiotic resistance and reduction in new antibiotic discovery, it is imperative to investigate alternative therapeutic regimens that complement the use of current antibiotic treatment strategies. As stated by the WHO Director-General, “On current trends, common diseases may become untreatable. Doctors facing patients will have to say, Sorry, there is nothing I can do for you.” Copyright © 2018 Bohlmann et al.


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.