The complete mitochondrial genome of a cold seep gastropod Phymorhynchus buccinoides (Neogastropoda: Conoidea: Raphitomidae)

Phymorhynchus is a genus of deep-sea snails that are most distributed in hydrothermal vent or cold seep environments. In this study, we presented the complete mitochondrial genome of P. buccinoides, a cold seep snail from the South China Sea. It is the first mitochondrial genome of a cold seep member of the superfamily Conoidea. The mitochondrial genome is 15,764 bp in length, and contains 13 protein-coding genes (PCGs), 2 rRNA genes, and 22 tRNA genes. These genes are encoded on the positive strand, except for 8 tRNA genes that are encoded on the negative strand. The start codon ATG and 3 types of stop codons, TAA, TAG and the truncated termination codon T, are used in the 13 PCGs. All 13 PCGs in the 26 species of Conoidea share the same gene order, while several tRNA genes have been translocated. Phylogenetic analysis revealed that P. buccinoides clustered with Typhlosyrinx sp., Eubela sp., and Phymorhynchus sp., forming the Raphitomidae clade, with high support values. Positive selection analysis showed that a residue located in atp6 (18 S) was identified as the positively selected site with high posterior probabilities, suggesting potential adaption to the cold seep environment. Overall, our data will provide a useful resource on the evolutionary adaptation of cold seep snails for future studies.

Inducible modulation of miR-204 levels in a zebrafish melanoma model

Here, we present miniCoopR-I, an inducible upgrade of the constitutive miniCoopR vector. We developed miniCoopR-I-sponge-204 and miniCoopR-I-pre-miR-204 vectors and we successfully tested them for their ability to achieve time- (embryo/juvenile/adult) and space- (melanocytic lineage) restricted inhibition/overexpression of miR-204, a positive modulator of pigmentation previously discovered by us. Furthermore, melanoma-free survival curves performed on induced fish at the adult stage indicate that miR-204 overexpression accelerates the development of BRAFV600E-driven melanoma. miniCoopR-I allows study of the impact that coding and non-coding modulators of pigmentation exert on melanomagenesis in adult zebrafish, uncoupling it from the impact that they exert on melanogenesis during embryonic development.This article has an associated First Person interview with the first author of the paper.

Improved Brassica oleracea JZS assembly reveals significant changing of LTR-RT dynamics in different morphotypes

Brassica oleracea is an important vegetable crop that has provided ancestor genomes of the two most important Brassica oil crops, Brassica napus and Brassica carinata. The current B. oleracea reference genome (JZS, also named 02-12) displays problems of large mis-assemblies, low sequence continuity, and low assembly integrity, thus limiting genomic analysis. We reported an updated assembly of the B. oleracea reference genome (JZS v2) obtained through single-molecule sequencing and chromosome conformation capture technologies. We assembled an additional 83.16 Mb of genomic sequences, and the updated genome features a contig N50 size of 2.37 Mb, representing an ~ 88-fold improvement. We detected a new round of long terminal repeat retrotransposon (LTR-RT) burst in the new assembly. Comparative analysis with the reported genome sequences of two other genomes of B. oleracea (TO1000 and HDEM) identified extensive gene order and gene structural variation. In addition, we found that the genome-specific amplification of Gypsy-like LTR-RTs occurred around 0-1 million years ago (MYA). In particular, the athila, tat, and Del families were extensively amplified in JZS around 0-1 MYA. Moreover, we identified that the syntenic genes were modified due to the insertion of genome-specific LTR-RTs. These results indicated that the genome-specific LTR-RT dynamics were associated with genome diversification in B. oleracea.

Molecular determinants that govern scaRNA processing by Drosha/DGCR8

The Cajal body (CB) is a subnuclear domain that participates in the biogenesis of many different types of ribonucleoproteins (RNPs), including small nuclear RNPs (snRNPs), small Cajal body-specific RNPs (scaRNPs) and telomerase. Most scaRNAs, the RNA component of scaRNPs, accumulate in CBs. However, there are three scaRNAs (scaRNA 2, 9, and 17) that are known to be processed into small, nucleolar-enriched fragments. Evidence suggests that these fragments are packaged into a new class of RNPs, called regulatory RNPs (regRNPs), and may modify small nucleolar RNP (snoRNP) activity, thus playing a role in rRNA modification. However, the mechanism by which these fragments are produced is unknown. Previous work has reported the involvement of Drosha and DGCR8 in the cleavage of primary-scaRNA9. Here, we expand on that knowledge by identifying sequence elements necessary for the efficient production of these RNA fragments and demonstrate that primary scaRNA 2 and 17 are also processed by the Drosha-DGCR8 complex. Collectively, our work establishes new factors in the scaRNP biogenesis pathway and adds to the ever-expanding list of noncanonical functions for the microprocessor complex.

Sequence comparison of the mitochondrial genomes in two species of the genus Nerita (Gastropoda: Neritimorpha: Neritidae): phylogenetic implications and divergence time estimation for Neritimorpha

Many Nerita species live in warm-water environments, and they are some of the few organisms from the intertidal zone that can live in both freshwater and seawater. Previous comparative studies of the mitogenomes of Nerita species suggest that the genome rearrangements are very conservative. Generally, the species possess a set of similar mitochondrial gene arrangements, but nucleotide sequences can be used to elucidate phylogenetic relationships at various levels of divergence. Here, the mitogenomes of Nerita undata and Nerita balteata were sequenced and found to be 15,583 bp and 15,571 bp, respectively. The mitogenomes of both species contain 13 protein-coding genes, 22 tRNA genes, and two rRNA genes. The nucleotides of the two mitogenomes are highly similar, with the same gene composition and genomic organization as those present in other Nerita species. The tRNA secondary structures were different from those of other gastropods: trnS2 is not folded into typical secondary structures, and the dihydrouridine (DHU) arm simply forms a loop. The phylogenetic analysis showed that Neritimorpha is a sister group of Vetigastropoda and Caenogastropoda. Nerita balteata is a sister group of Nerita versicolor and Nerita undata, and all three species belong to Neritimorpha. This study contributes towards the comparative mitogenomic analysis of Neritidae and phylogenetic considerations among Neritimorpha species. The estimation of divergence time revealed that the two Nerita species were differentiated in the late Paleogene of the Cenozoic Era, and their evolution may be related to environmental changes.

Stable integration of an optimized inducible promoter system enables spatiotemporal control of gene expression throughout avian development

Precisely altering gene expression is critical for understanding molecular processes of embryogenesis. Although some tools exist for transgene misexpression in developing chick embryos, we have refined and advanced them by simplifying and optimizing constructs for spatiotemporal control. To maintain expression over the entire course of embryonic development we use an enhanced piggyBac transposon system that efficiently integrates sequences into the host genome. We also incorporate a DNA targeting sequence to direct plasmid translocation into the nucleus and a D4Z4 insulator sequence to prevent epigenetic silencing. We designed these constructs to minimize their size and maximize cellular uptake, and to simplify usage by placing all of the integrating sequences on a single plasmid. Following electroporation of stage HH8.5 embryos, our tetracycline-inducible promoter construct produces robust transgene expression in the presence of doxycycline at any point during embryonic development in ovo or in culture. Moreover, expression levels can be modulated by titrating doxycycline concentrations and spatial control can be achieved using beads or gels. Thus, we have generated a novel, sensitive, tunable, and stable inducible-promoter system for high-resolution gene manipulation in vivo.

A targeted reactivation of latent HIV-1 using an activator vector in patient samples from acute infection

BACKGROUND

During combined anti-retroviral treatment, a latent HIV reservoir persists within resting memory CD4 T cells that initiates viral recrudescence upon treatment interruption. Strategies for HIV-1 cure have largely focused on latency reversing agents (LRAs) capable of reactivating and eliminating this viral reservoir. Previously investigated LRAs have largely failed to achieve a robust latency reversal sufficient for reduction of latent HIV pool or the potential of virus-free remission in the absence of treatment.

METHODS

We utilize a polyvalent virus-like particle (VLP) formulation called Activator Vector (ACT-VEC) to 'shock' provirus into transcriptional activity. Ex vivo co-culture experiments were used to evaluate the efficacy of ACT-VEC in relation to other LRAs in individuals diagnosed and treated during the acute stage of infection. IFN-γ ELISpot, qRT-PCR and Illumina MiSeq were used to evaluate antigenicity, latency reversal, and diversity of induced virus respectively.

FINDINGS

Using samples from HIV+ patients diagnosed and treated at acute/early infection, we demonstrate that ACT-VEC can reverse latency in HIV infected CD4 T cells to a greater extent than other major recall antigens as stimuli or even mitogens such as PMA/Iono. Furthermore, ACT-VEC activates more latent HIV-1 than clinically tested HDAC inhibitors or protein kinase C agonists.

INTERPRETATION

Taken together, these results show that ACT-VEC can induce HIV reactivation from latently infected CD4 T cells collected from participants on first line combined antiretroviral therapy for at least two years after being diagnosed and treated at acute/early stage of infection. These findings could provide guidance to possible targeted cure strategies and treatments.

FUNDING

NIH and CIHR.

Phylogenetic Reconstruction Based on Synteny Block and Gene Adjacencies

Gene order can be used as an informative character to reconstruct phylogenetic relationships between species independently from the local information present in gene/protein sequences. PhyChro is a reconstruction method based on chromosomal rearrangements, applicable to a wide range of eukaryotic genomes with different gene contents and levels of synteny conservation. For each synteny breakpoint issued from pairwise genome comparisons, the algorithm defines two disjoint sets of genomes, named partial splits, respectively, supporting the two block adjacencies defining the breakpoint. Considering all partial splits issued from all pairwise comparisons, a distance between two genomes is computed from the number of partial splits separating them. Tree reconstruction is achieved through a bottom-up approach by iteratively grouping sister genomes minimizing genome distances. PhyChro estimates branch lengths based on the number of synteny breakpoints and provides confidence scores for the branches. PhyChro performance is evaluated on two data sets of 13 vertebrates and 21 yeast genomes by using up to 130,000 and 179,000 breakpoints, respectively, a scale of genomic markers that has been out of reach until now. PhyChro reconstructs very accurate tree topologies even at known problematic branching positions. Its robustness has been benchmarked for different synteny block reconstruction methods. On simulated data PhyChro reconstructs phylogenies perfectly in almost all cases, and shows the highest accuracy compared with other existing tools. PhyChro is very fast, reconstructing the vertebrate and yeast phylogenies in <15 min.

Generation of Marker-Free pbd-2 Knock-in Pigs Using the CRISPR/Cas9 and Cre/loxP Systems

Porcine β-defensin 2 (PBD-2), expressed by different tissues of pigs, is a multifunctional cationic peptide with antimicrobial, immunomodulatory and growth-promoting abilities. As the latest generation of genome-editing tool, CRISPR/Cas9 system makes it possible to enhance the expression of PBD-2 in pigs by site-specific knock-in of pbd-2 gene into the pig genome. In this study, we aimed to generate marker-free pbd-2 knock-in pigs using the CRISPR/Cas9 and Cre/loxP systems. Two copies of pbd-2 gene linked by a T2A sequence were inserted into the porcine Rosa26 locus through CRISPR/Cas9-mediated homology-directed repair. The floxed selectable marker gene neoR, used for G418 screening of positive cell clones, was removed by cell-penetrating Cre recombinase with a recombination efficiency of 48.3%. Cloned piglets were produced via somatic cell nuclear transfer and correct insertion of pbd-2 genes was confirmed by PCR and Southern blot. Immunohistochemistry and immunofluorescence analyses indicated that expression levels of PBD-2 in different tissues of transgenic (TG) piglets were significantly higher than those of their wild-type (WT) littermates. Bactericidal assays demonstrated that there was a significant increase in the antimicrobial properties of the cell culture supernatants of porcine ear fibroblasts from the TG pigs in comparison to those from the WT pigs. Altogether, our study improved the protein expression level of PBD-2 in pigs by site-specific integration of pbd-2 into the pig genome, which not only provided an effective pig model to study the anti-infection mechanisms of PBD-2 but also a promising genetic material for the breeding of disease-resistant pigs.

Insights into Mobile Genetic Elements of the Biocide-Degrading Bacterium Pseudomonas nitroreducens HBP-1

The sewage sludge isolate Pseudomonas nitroreducens HBP-1 was the first bacterium known to completely degrade the fungicide 2-hydroxybiphenyl. PacBio and Illumina whole-genome sequencing revealed three circular DNA replicons: a chromosome and two plasmids. Plasmids were shown to code for putative adaptive functions such as heavy metal resistance, but with unclarified ability for self-transfer. About one-tenth of strain HBP-1's chromosomal genes are likely of recent horizontal influx, being part of genomic islands, prophages and integrative and conjugative elements (ICEs). P. nitroreducens carries two large ICEs with different functional specialization, but with homologous core structures to the well-known ICEclc of Pseudomonas knackmussii B13. The variable regions of ICEPni1 (96 kb) code for, among others, heavy metal resistances and formaldehyde detoxification, whereas those of ICEPni2 (171 kb) encodes complete meta-cleavage pathways for catabolism of 2-hydroxybiphenyl and salicylate, a protocatechuate pathway and peripheral enzymes for 4-hydroxybenzoate, ferulate, vanillin and vanillate transformation. Both ICEs transferred at frequencies of 10^-6-10^-8 per P. nitroreducens HBP-1 donor into Pseudomonas putida, where they integrated site specifically into tRNA^Gly-gene targets, as expected. Our study highlights the underlying determinants and mechanisms driving dissemination of adaptive properties allowing bacterial strains to cope with polluted environments.

The complete mitochondrial genome of the Columbia lance nematode, Hoplolaimus columbus, a major agricultural pathogen in North America

BACKGROUND

The plant-parasitic nematode Hoplolaimus columbus is a pathogen that uses a wide range of hosts and causes substantial yield loss in agricultural fields in North America. This study describes, for the first time, the complete mitochondrial genome of H. columbus from South Carolina, USA.

METHODS

The mitogenome of H. columbus was assembled from Illumina 300 bp pair-end reads. It was annotated and compared to other published mitogenomes of plant-parasitic nematodes in the superfamily Tylenchoidea. The phylogenetic relationships between H. columbus and other 6 genera of plant-parasitic nematodes were examined using protein-coding genes (PCGs).

RESULTS

The mitogenome of H. columbus is a circular AT-rich DNA molecule 25,228 bp in length. The annotation result comprises 12 PCGs, 2 ribosomal RNA genes, and 19 transfer RNA genes. No atp8 gene was found in the mitogenome of H. columbus but long non-coding regions were observed in agreement to that reported for other plant-parasitic nematodes. The mitogenomic phylogeny of plant-parasitic nematodes in the superfamily Tylenchoidea agreed with previous molecular phylogenies. Mitochondrial gene synteny in H. columbus was unique but similar to that reported for other closely related species.

CONCLUSIONS

The mitogenome of H. columbus is unique within the superfamily Tylenchoidea but exhibits similarities in both gene content and synteny to other closely related nematodes. Among others, this new resource will facilitate population genomic studies in lance nematodes from North America and beyond.

A Highly Efficacious PS Gene Editing System Corrects Metabolic and Neurological Complications of Mucopolysaccharidosis Type I

Our previous study delivered zinc finger nucleases to treat mice with mucopolysaccharidosis type I (MPS I), resulting in a phase I/II clinical trial (ClinicalTrials.gov: NCT02702115). However, in the clinical trial, the efficacy needs to be improved due to the low transgene expression level. To this end, we designed a proprietary system (PS) gene editing approach with CRISPR to insert a promoterless α-l-iduronidase (IDUA) cDNA sequence into the albumin locus of hepatocytes. In this study, adeno-associated virus 8 (AAV8) vectors delivering the PS gene editing system were injected into neonatal and adult MPS I mice. IDUA enzyme activity in the brain significantly increased, while storage levels were normalized. Neurobehavioral tests showed that treated mice had better memory and learning ability. Also, histological analysis showed efficacy reflected by the absence of foam cells in the liver and vacuolation in neuronal cells. No vector-associated toxicity or increased tumorigenesis risk was observed. Moreover, no off-target effects were detected through the unbiased genome-wide unbiased identification of double-stranded breaks enabled by sequencing (GUIDE-seq) analysis. In summary, these results showed the safety and efficacy of the PS in treating MPS I and paved the way for clinical studies. Additionally, as a therapeutic platform, the PS has the potential to treat other lysosomal diseases.

The nuclear and mitochondrial genomes of Frieseomelitta varia - a highly eusocial stingless bee (Meliponini) with a permanently sterile worker caste

BACKGROUND

Most of our understanding on the social behavior and genomics of bees and other social insects is centered on the Western honey bee, Apis mellifera. The genus Apis, however, is a highly derived branch comprising less than a dozen species, four of which genomically characterized. In contrast, for the equally highly eusocial, yet taxonomically and biologically more diverse Meliponini, a full genome sequence was so far available for a single Melipona species only. We present here the genome sequence of Frieseomelitta varia, a stingless bee that has, as a peculiarity, a completely sterile worker caste.

RESULTS

The assembly of 243,974,526 high quality Illumina reads resulted in a predicted assembled genome size of 275 Mb composed of 2173 scaffolds. A BUSCO analysis for the 10,526 predicted genes showed that these represent 96.6% of the expected hymenopteran orthologs. We also predicted 169,371 repetitive genomic components, 2083 putative transposable elements, and 1946 genes for non-coding RNAs, largely long non-coding RNAs. The mitochondrial genome comprises 15,144 bp, encoding 13 proteins, 22 tRNAs and 2 rRNAs. We observed considerable rearrangement in the mitochondrial gene order compared to other bees. For an in-depth analysis of genes related to social biology, we manually checked the annotations for 533 automatically predicted gene models, including 127 genes related to reproductive processes, 104 to development, and 174 immunity-related genes. We also performed specific searches for genes containing transcription factor domains and genes related to neurogenesis and chemosensory communication.

CONCLUSIONS

The total genome size for F. varia is similar to the sequenced genomes of other bees. Using specific prediction methods, we identified a large number of repetitive genome components and long non-coding RNAs, which could provide the molecular basis for gene regulatory plasticity, including worker reproduction. The remarkable reshuffling in gene order in the mitochondrial genome suggests that stingless bees may be a hotspot for mtDNA evolution. Hence, while being just the second stingless bee genome sequenced, we expect that subsequent targeting of a selected set of species from this diverse clade of highly eusocial bees will reveal relevant evolutionary signals and trends related to eusociality in these important pollinators.

AAV-CRISPR Gene Editing Is Negated by Pre-existing Immunity to Cas9

Adeno-associated viral (AAV) vectors are a leading candidate for the delivery of CRISPR-Cas9 for therapeutic genome editing in vivo. However, AAV-based delivery involves persistent expression of the Cas9 nuclease, a bacterial protein. Recent studies indicate a high prevalence of neutralizing antibodies and T cells specific to the commonly used Cas9 orthologs from Streptococcus pyogenes (SpCas9) and Staphylococcus aureus (SaCas9) in humans. We tested in a mouse model whether pre-existing immunity to SaCas9 would pose a barrier to liver genome editing with AAV packaging CRISPR-Cas9. Although efficient genome editing occurred in mouse liver with pre-existing SaCas9 immunity, this was accompanied by an increased proportion of CD8+ T cells in the liver. This cytotoxic T cell response was characterized by hepatocyte apoptosis, loss of recombinant AAV genomes, and complete elimination of genome-edited cells, and was followed by compensatory liver regeneration. Our results raise important efficacy and safety concerns for CRISPR-Cas9-based in vivo genome editing in the liver.

Complete Mitochondrial Genomes of Two Toxin-Accumulated Nassariids (Neogastropoda: Nassariidae: Nassarius) and Their Implication for Phylogeny

The Indo-Pacific nassariids (genus Nassarius) possesses the highest diversity within the family Nassariidae. However, the previous shell or radula-based classification of Nassarius is quite confusing due to the homoplasy of certain morphological characteristics. The toxin accumulators Nassarius glans and Nassarius siquijorensis are widely distributed in the subtidal regions of the Indo-Pacific Ocean. In spite of their biological significance, the phylogenetic positions of N. glans and N. siquijorensis are still undetermined. In the present study, the complete mitochondrial genomes of N. glans and N. siquijorensis were sequenced. The present mitochondrial genomes were 15,296 and 15,337 bp in length, respectively, showing negative AT skews and positive GC skews as well as a bias of AT rich on the heavy strand. They contained 13 protein coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and several noncoding regions, and their gene order was identical to most caenogastropods. Based on the nucleotide sequences combining 13 protein coding genes and two rRNA genes, a well-supported phylogeny of Nassarius was reconstructed, and several morphological synapomorphies were observed corresponding to the phylogenetic framework. In addition, the sister group relationship between N. variciferus and the remaining toxin-accumulated nassariids was determined, suggesting that the phylogeny might be related to their diet. The divergence time estimation analysis revealed a correlation between speciation events of nassariids and glacial cycles during the Pliocene-Pleistocene epoch.

Assessment of faithful interleukin-3 production by novel bicistronic interleukin-3 reporter mice

Interleukin-3 (IL-3) is an important hematopoietic growth factor and immunregulatory cytokine. Although activated T helper cells represent a main source of IL-3, other cell types have been reported to express this cytokine. However, precise identification and quantification of the cells that produce IL-3 in vivo have not been performed. Therefore, we used a CRISPR/Cas approach to engineer mice containing a bicistronic mRNA linking a readily identifiable reporter, enhanced green fluorescent protein (ZsGreen1), to IL-3 expression. To characterize these novel reporter mice, we first examined ZsGreen1 expression by CD4 T cells subsets primed and activated in vitro. We found that activated Th1 cells expressed ∼4-fold higher levels of ZsGreen1 as compared to Th0 and Th2 cells. Endogenous IL-3 expression remained intact although reporter Th1 cells secreted ∼33 % less IL-3 than similarly activated wild-type cells. To characterize the ability of reporter mice to accurately mark IL-3-producing cells in vivo, we infected mice with Nippostrongylus brasiliensis. Low but significant numbers of ZsGreen1+ CD4 T cells were detected in the mesenteric lymph nodes and lung following both primary and secondary infection. No difference in basophil and intestinal mast cell numbers were observed between infected reporter and wild-type mice indicating that reporter mice secreted IL-3 levels in vivo that results in IL-3-driven biological activities which are indistinguishable from those observed in corresponding wild-type mice. These IL-3 reporter mice will be a valuable resource to investigate IL-3-dependent immune responses in vivo.

How the COVID-19 Overcomes the Battle? An Approach to Virus Structure

Coronaviruses primarily cause zoonotic infections, however in the past few decades several interspecies transmissions have occurred, the last one by SARS-CoV-2, causing COVID-19 pandemic, posing serious threat to global health. The SARS-CoV-2 spike (S) protein plays an important role in viral attachment, fusion and entry. However, other structural and non-structural SARS-CoV-2 proteins are potential influencers in virus pathogenicity. Among these proteins; Orf3, Orf8, and Orf10 show the least homology to SARSCoV proteins and therefore should be further studied for their abilities to modulate antiviral and inflammatory responses. Here, we discuss how SARS-COV-2 interacts with our immune system.

The complete mitochondrial genome of the endangered Assam Roofed Turtle, Pangshura sylhetensis (Testudines: Geoemydidae): Genomic features and phylogeny

The Assam Roofed Turtle, Pangshura sylhetensis is an endangered and least studied species endemic to India and Bangladesh. The present study decodes the first complete mitochondrial genome of P. sylhetensis (16,568 bp) by using next-generation sequencing. The assembly encodes 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (CR). Most of the genes were encoded on the majority strand, except NADH dehydrogenase subunit 6 (nad6) and eight tRNAs. All PCGs start with an ATG initiation codon, except for Cytochrome oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 5 (nad5), which both start with GTG codon. The study also found the typical cloverleaf secondary structures in most of the predicted tRNA structures, except for serine (trnS1) which lacks of conventional DHU arm and loop. Both Bayesian and maximum-likelihood phylogenetic inference using 13 concatenated PCGs demonstrated strong support for the monophyly of all 52 Testudines species within their respective families and revealed Batagur trivittata as the nearest neighbor of P. sylhetensis. The mitogenomic phylogeny with other amniotes is congruent with previous research, supporting the sister relationship of Testudines and Archosaurians (birds and crocodilians). Additionally, the mitochondrial Gene Order (GO) analysis indicated plesiomorphy with the typical vertebrate GO in most of the Testudines species.

Analysis of local genome rearrangement improves resolution of ancestral genomic maps in plants

BACKGROUND

Computationally inferred ancestral genomes play an important role in many areas of genome research. We present an improved workflow for the reconstruction from highly diverged genomes such as those of plants.

RESULTS

Our work relies on an established workflow in the reconstruction of ancestral plants, but improves several steps of this process. Instead of using gene annotations for inferring the genome content of the ancestral sequence, we identify genomic markers through a process called genome segmentation. This enables us to reconstruct the ancestral genome from hundreds of thousands of markers rather than the tens of thousands of annotated genes. We also introduce the concept of local genome rearrangement, through which we refine syntenic blocks before they are used in the reconstruction of contiguous ancestral regions. With the enhanced workflow at hand, we reconstruct the ancestral genome of eudicots, a major sub-clade of flowering plants, using whole genome sequences of five modern plants.

CONCLUSIONS

Our reconstructed genome is highly detailed, yet its layout agrees well with that reported in Badouin et al. (2017). Using local genome rearrangement, not only the marker-based, but also the gene-based reconstruction of the eudicot ancestor exhibited increased genome content, evidencing the power of this novel concept.

A Preliminary Study of the Virome of the South American Free-Tailed Bats (Tadarida brasiliensis) and Identification of Two Novel Mammalian Viruses

Bats provide important ecosystem services as pollinators, seed dispersers, and/or insect controllers, but they have also been found harboring different viruses with zoonotic potential. Virome studies in bats distributed in Asia, Africa, Europe, and North America have increased dramatically over the past decade, whereas information on viruses infecting South American species is scarce. We explored the virome of Tadarida brasiliensis, an insectivorous New World bat species inhabiting a maternity colony in Rosario (Argentina), by a metagenomic approach. The analysis of five pooled oral/anal swab samples indicated the presence of 43 different taxonomic viral families infecting a wide range of hosts. By conventional nucleic acid detection techniques and/or bioinformatics approaches, the genomes of two novel viruses were completely covered clustering into the Papillomaviridae (Tadarida brasiliensis papillomavirus type 1, TbraPV1) and Genomoviridae (Tadarida brasiliensis gemykibivirus 1, TbGkyV1) families. TbraPV1 is the first papillomavirus type identified in this host and the prototype of a novel genus. TbGkyV1 is the first genomovirus reported in New World bats and constitutes a new species within the genus Gemykibivirus. Our findings extend the knowledge about oral/anal viromes of a South American bat species and contribute to understand the evolution and genetic diversity of the novel characterized viruses.

De novo sequence of the mitochondrial genome of Tyrophagus putrescentiae (Acari: Sarcoptiformes) including 22 tRNA sequences and the largest non-coding region

In this study, we de novo sequenced and analyzed the circular mitochondrial genome (mitogenome) of Tyrophagus putrescentiae. It was 14,156 bp long and contained a complete set of 37 genes, contrary to the initial published sequences; it included 22 tRNA sequences and the largest non-coding region. The mtDNA gene order of T. putrescentiae was found to be identical to that of Aleuroglyphus ovatus, Caloglyphus berlesei, and Rhizoglyphus robini (all Acaroidea). Most tRNAs of T. putrescentiae lack at least a D-arm or T-arm. Tyrophagus putrescentiae tRNAs also shared considerable structural and sequence similarity with the tRNAs of other reported Acaroidea species that have the full set of tRNAs. The largest non-coding region was located between trnF and trnS1, and it contained a microsatellite-like (AT)_n sequence, short palindromic sequences, and several hairpin loops, as observed in other reported Acaroidea species (excepting Tyrophagus longior).

Bacillus anthracis Responds to Targocil-Induced Envelope Damage through EdsRS Activation of Cardiolipin Synthesis

Bacillus anthracis is a spore-forming bacterium that causes devastating infections and has been used as a bioterror agent. This pathogen can survive hostile environments through the signaling activity of two-component systems, which couple environmental sensing with transcriptional activation to initiate a coordinated response to stress. In this work, we describe the identification of a two-component system, EdsRS, which mediates the B. anthracis response to the antimicrobial compound targocil. Targocil is a cell envelope-targeting compound that is toxic to B. anthracis at high concentrations. Exposure to targocil causes damage to the cellular barrier and activates EdsRS to induce expression of a previously uncharacterized cardiolipin synthase, which we have named ClsT. Both EdsRS and ClsT are required for protection against targocil-dependent damage. Induction of clsT by EdsRS during targocil treatment results in an increase in cardiolipin levels, which protects B. anthracis from envelope damage. Together, these results reveal that a two-component system signaling response to an envelope-targeting antimicrobial induces production of a phospholipid associated with stabilization of the membrane. Cardiolipin is then used to repair envelope damage and promote B. anthracis viability.IMPORTANCE Compromising the integrity of the bacterial cell barrier is a common action of antimicrobials. Targocil is an antimicrobial that is active against the bacterial envelope. We hypothesized that Bacillus anthracis, a potential weapon of bioterror, senses and responds to targocil to alleviate targocil-dependent cell damage. Here, we show that targocil treatment increases the permeability of the cellular envelope and is particularly toxic to B. anthracis spores during outgrowth. In vegetative cells, two-component system signaling through EdsRS is activated by targocil. This results in an increase in the production of cardiolipin via a cardiolipin synthase, ClsT, which restores the loss of barrier function, thereby reducing the effectiveness of targocil. By elucidating the B. anthracis response to targocil, we have uncovered an intrinsic mechanism that this pathogen employs to resist toxicity and have revealed therapeutic targets that are important for bacterial defense against structural damage.

Development of a Honey Bee RNA Virus Vector Based on the Genome of a Deformed Wing Virus

We developed a honey bee RNA-virus vector based on the genome of a picorna-like Deformed wing virus (DWV), the main viral pathogen of the honey bee (Apis mellifera). To test the potential of DWV to be utilized as a vector, the 717 nt sequence coding for the enhanced green fluorescent protein (eGFP), flanked by the peptides targeted by viral protease, was inserted into an infectious cDNA clone of DWV in-frame between the leader protein and the virus structural protein VP2 genes. The in vitro RNA transcripts from egfp-tagged DWV cDNA clones were infectious when injected into honey bee pupae. Stable DWV particles containing genomic RNA of the recovered DWV with egfp inserts were produced, as evidenced by cesium chloride density gradient centrifugation. These particles were infectious to honey bee pupae when injected intra-abdominally. Fluorescent microscopy showed GFP expression in the infected cells and Western blot analysis demonstrated accumulation of free eGFP rather than its fusions with DWV leader protein (LP) and/or viral protein (VP) 2. Analysis of the progeny egfp-tagged DWV showed gradual accumulation of genome deletions for egfp, providing estimates for the rate of loss of a non-essential gene an insect RNA virus genome during natural infection.

Cytogenetics of Blastic Plasmacytoid Dendritic Cell Neoplasm: Chromosomal Rearrangements and DNA Copy-Number Alterations

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a skin-tropic hematopoietic malignancy. Approximately 60% of cases with analyzable karyotyping results show complex karyotypes. Losses are more frequently found than copy-number gains. Recurrently deleted regions include tumor suppressor genes. No specific chromosomal abnormalities have been demonstrated in BPDCN, but genomic rearrangements involving the MYB family genes and MYC were identified. One-third of cases of BPDCN harbor the 8q24 rearrangement, most frequently with 6p21 harboring RUNX2, which is associated with immunoblastoid cytomorphology and MYC expression. MYB rearrangement is detected in 20% of patients with BPDCN. We review copy-number alterations and chromosomal rearrangements.

ORF-based binarized structure network analysis of plasmids (OSNAp), a novel approach to core gene-independent plasmid phylogeny

OBJECTIVES

Systematic comparison of multiple plasmids remains challenging. We aimed to develop a new method for phylogenetic analysis of plasmids, open reading frame (ORF)-based binarized structure network analysis of plasmids (OSNAp).

METHODS

With the OSNAp, the genetic structures of plasmids in a given plasmid group are expressed as binary sequences based on the presence or absence of ORFs regardless of their positions or directions. As a proof-of-concept, ORFs were collected from 101 complete I1 plasmid sequences, and their corresponding binary sequences were generated. A tree was generated using the neighbor-net, an algorithm for constructing phylogenetic networks based on distance between taxa, to visualize the plasmid phylogeny drawn from binary sequences. The results were compared with those of plasmid sequence types (pSTs) defined by plasmid multilocus sequence typing (pMLST).

RESULTS

All I1 plasmids were placed on the phylogenetic tree constructed from the binary sequences. Most plasmids belonging to the same pSTs had Dice indices of ≥0.95 and were placed in the same OSNAp split. On the other hand, pST12 plasmids were distributed on separate splits due to differences in ORFs not used in pMLST, suggesting improved differentiation of the plasmids with OSNAp compared with pMLST.

CONCLUSION

OSNAp is a novel holistic approach to assess relatedness of a population of plasmids in a given plasmid group based on nucleotide sequence data. It provides higher discrimination than pMLST, which may prove useful in tracing bacteria that harbor plasmids of shared origins.