Scaling read aligners to hundreds of threads on general-purpose processors

Protocol for assay of transposase accessible chromatin sequencing in non-model species

The assay for transposase accessible chromatin (ATAC-seq) is a method for mapping genome-wide chromatin accessibility. Coupled with high-throughput sequencing, it enables integrative epigenomics analyses. ATAC-seq requires direct access to cell nuclei, a major challenge in non-model species such as small invertebrates, whose soft tissue is surrounded by a protective exoskeleton. Here, we present modifications of the ATAC-seq protocol for applications in small crustaceans, extending applications to non-model species.

For complete information on the use and execution of this protocol, please refer to Buenrostro et al. (2013).

•

ATAC-seq modified protocol for applications in non-model species

•

Transposase titration identifies thresholds for optimal transposition at lower costs

•

Minimal number of cells/tissue is identified for low input reactions

•

Proof of concept ATAC-seq analysis for the waterflea Daphnia magna

Transcriptional programming drives Ibrutinib-resistance evolution in mantle cell lymphoma

Ibrutinib, a bruton's tyrosine kinase (BTK) inhibitor, provokes robust clinical responses in aggressive mantle cell lymphoma (MCL), yet many patients relapse with lethal Ibrutinib-resistant (IR) disease. Here, using genomic, chemical proteomic, and drug screen profiling, we report that enhancer remodeling-mediated transcriptional activation and adaptive signaling changes drive the aggressive phenotypes of IR. Accordingly, IR MCL cells are vulnerable to inhibitors of the transcriptional machinery and especially so to inhibitors of cyclin-dependent kinase 9 (CDK9), the catalytic subunit of the positive transcription elongation factor b (P-TEFb) of RNA polymerase II (RNAPII). Further, CDK9 inhibition disables reprogrammed signaling circuits and prevents the emergence of IR in MCL. Finally, and importantly, we find that a robust and facile ex vivo image-based functional drug screening platform can predict clinical therapeutic responses of IR MCL and identify vulnerabilities that can be targeted to disable the evolution of IR.

Genome Sequence of Streptomyces sp. Strain HB-N217, Isolated from the Marine Sponge Forcepia sp

The genome sequence of the Forcepia sponge-derived bacterium Streptomyces sp. strain HB-N217 was determined, with approximately 8.25 Mbp and a G+C content of 72.1%. Thirty biosynthetic gene clusters that bear the capability to produce secondary metabolites were predicted. The results will aid marine natural product chemistry and sponge-microbe association studies.

The genome sequence of the Forcepia sponge-derived bacterium Streptomyces sp. strain HB-N217 was determined, with approximately 8.25 Mbp and a G+C content of 72.1%. Thirty biosynthetic gene clusters that bear the capability to produce secondary metabolites were predicted. The results will aid marine natural product chemistry and sponge-microbe association studies.

Portable nanopore analytics: are we there yet?

MOTIVATION

Oxford Nanopore technologies (ONT) add miniaturization and real time to high-throughput sequencing. All available software for ONT data analytics run on cloud/clusters or personal computers. Instead, a linchpin to true portability is software that works on mobile devices of internet connections. Smartphones' and tablets' chipset/memory/operating systems differ from desktop computers, but software can be recompiled. We sought to understand how portable current ONT analysis methods are.

RESULTS

Several tools, from base-calling to genome assembly, were ported and benchmarked on an Android smartphone. Out of 23 programs, 11 succeeded. Recompilation failures included lack of standard headers and unsupported instruction sets. Only DSK, BCALM2 and Kraken were able to process files up to 16 GB, with linearly scaling CPU-times. However, peak CPU temperatures were high. In conclusion, the portability scenario is not favorable. Given the fast market growth, attention of developers to ARM chipsets and Android/iOS is warranted, as well as initiatives to implement mobile-specific libraries.

AVAILABILITY AND IMPLEMENTATION

The source code is freely available at: https://github.com/marco-oliva/portable-nanopore-analytics.

Programmable human histone phosphorylation and gene activation using a CRISPR/Cas9-based chromatin kinase

Histone phosphorylation is a ubiquitous post-translational modification that allows eukaryotic cells to rapidly respond to environmental stimuli. Despite correlative evidence linking histone phosphorylation to changes in gene expression, establishing the causal role of this key epigenomic modification at diverse loci within native chromatin has been hampered by a lack of technologies enabling robust, locus-specific deposition of endogenous histone phosphorylation. To address this technological gap, here we build a programmable chromatin kinase, called dCas9-dMSK1, by directly fusing nuclease-null CRISPR/Cas9 to a hyperactive, truncated variant of the human MSK1 histone kinase. Targeting dCas9-dMSK1 to human promoters results in increased target histone phosphorylation and gene activation and demonstrates that hyperphosphorylation of histone H3 serine 28 (H3S28ph) in particular plays a causal role in the transactivation of human promoters. In addition, we uncover mediators of resistance to the BRAF V600E inhibitor PLX-4720 in human melanoma cells using genome-scale screening with dCas9-dMSK1. Collectively, our findings enable a facile way to reshape human chromatin using CRISPR/Cas9-based epigenome editing and further define the causal link between histone phosphorylation and human gene activation.

Histone phosphorylation is a ubiquitous post-translational modification. Here the authors present a programmable chromatin kinase, dCas9-dMSK1, that enables controlled histone phosphorylation and specific gene activation.

Transcriptional response to the host cell environment of a multidrug-resistant Mycobacterium tuberculosis clonal outbreak Beijing strain reveals its pathogenic features

Tuberculosis is a global public health problem with emergence of multidrug-resistant infections. Previous epidemiological studies of tuberculosis in Thailand have identified a clonal outbreak multidrug-resistant strain of Mycobacterium tuberculosis in the Kanchanaburi province, designated “MKR superspreader”, and this particular strain later was found to also spread to other regions. In this study, we elucidated its biology through RNA-Seq analyses and identified a set of genes involved in cholesterol degradation to be up-regulated in the MKR during the macrophage cell infection, but not in the H37Rv reference strain. We also found that the bacterium up-regulated genes associated with the ESX-1 secretion system during its intracellular growth phase, while the H37Rv did not. All results were confirmed by qRT-PCR. Moreover, we showed that compounds previously shown to inhibit the mycobacterial ESX-1 secretion system and cholesterol utilisation, and FDA-approved drugs known to interfere with the host cholesterol transportation were able to decrease the intracellular survival of the MKR when compared to the untreated control, while not that of the H37Rv. Altogether, our findings suggested that such pathways are important for the MKR’s intracellular growth, and potentially could be targets for the discovery of new drugs against this emerging multidrug-resistant strain of M. tuberculosis.

High-Resolution Mapping of Transcription Initiation in the Asexual Stages of Toxoplasma gondii

Toxoplasma gondii is a common parasite of humans and animals, causing life-threatening disease in the immunocompromized, fetal abnormalities when contracted during gestation, and recurrent ocular lesions in some patients. Central to the prevalence and pathogenicity of this protozoan is its ability to adapt to a broad range of environments, and to differentiate between acute and chronic stages. These processes are underpinned by a major rewiring of gene expression, yet the mechanisms that regulate transcription in this parasite are only partially characterized. Deciphering these mechanisms requires a precise and comprehensive map of transcription start sites (TSSs); however, Toxoplasma TSSs have remained incompletely defined. To address this challenge, we used 5'-end RNA sequencing to genomically assess transcription initiation in both acute and chronic stages of Toxoplasma. Here, we report an in-depth analysis of transcription initiation at promoters, and provide empirically-defined TSSs for 7603 (91%) protein-coding genes, of which only 1840 concur with existing gene models. Comparing data from acute and chronic stages, we identified instances of stage-specific alternative TSSs that putatively generate mRNA isoforms with distinct 5' termini. Analysis of the nucleotide content and nucleosome occupancy around TSSs allowed us to examine the determinants of TSS choice, and outline features of Toxoplasma promoter architecture. We also found pervasive divergent transcription at Toxoplasma promoters, clustered within the nucleosomes of highly-symmetrical phased arrays, underscoring chromatin contributions to transcription initiation. Corroborating previous observations, we asserted that Toxoplasma 5' leaders are among the longest of any eukaryote studied thus far, displaying a median length of approximately 800 nucleotides. Further highlighting the utility of a precise TSS map, we pinpointed motifs associated with transcription initiation, including the binding sites of the master regulator of chronic-stage differentiation, BFD1, and a novel motif with a similar positional arrangement present at 44% of Toxoplasma promoters. This work provides a critical resource for functional genomics in Toxoplasma, and lays down a foundation to study the interactions between genomic sequences and the regulatory factors that control transcription in this parasite.

Genomic signatures of the evolution of defence against its natural enemies in the poisonous and medicinal plant Datura stramonium (Solanaceae)

Tropane alkaloids and terpenoids are widely used in the medicine and pharmaceutic industry and evolved as chemical defenses against herbivores and pathogens in the annual herb Datura stramonium (Solanaceae). Here, we present the first draft genomes of two plants from contrasting environments of D. stramonium. Using these de novo assemblies, along with other previously published genomes from 11 Solanaceae species, we carried out comparative genomic analyses to provide insights on the genome evolution of D. stramonium within the Solanaceae family, and to elucidate adaptive genomic signatures to biotic and abiotic stresses in this plant. We also studied, in detail, the evolution of four genes of D. stramonium-Putrescine N-methyltransferase, Tropinone reductase I, Tropinone reductase II and Hyoscyamine-6S-dioxygenase-involved in the tropane alkaloid biosynthesis. Our analyses revealed that the genomes of D. stramonium show signatures of expansion, physicochemical divergence and/or positive selection on proteins related to the production of tropane alkaloids, terpenoids, and glycoalkaloids as well as on R defensive genes and other important proteins related with biotic and abiotic pressures such as defense against natural enemies and drought.

Integrated analysis of the methylome and transcriptome of chickens with fatty liver hemorrhagic syndrome

Background

DNA methylation, a biochemical modification of cytosine, has an important role in lipid metabolism. Fatty liver hemorrhagic syndrome (FLHS) is a serious disease and is tightly linked to lipid homeostasis. Herein, we compared the methylome and transcriptome of chickens with and without FLHS.

Results

We found genome-wide dysregulated DNA methylation pattern in which regions up- and down-stream of gene body were hypo-methylated in chickens with FLHS. A total of 4155 differentially methylated genes and 1389 differentially expressed genes were identified. Genes were focused when a negative relationship between mRNA expression and DNA methylation in promoter and gene body were detected. Based on pathway enrichment analysis, we found expression of genes related to lipogenesis and oxygenolysis (e.g., PPAR signaling pathway, fatty acid biosynthesis, and fatty acid elongation) to be up-regulated with associated down-regulated DNA methylation. In contrast, genes related to cellular junction and communication pathways (e.g., vascular smooth muscle contraction, phosphatidylinositol signaling system, and gap junction) were inhibited and with associated up-regulation of DNA methylation.

Conclusions

In the current study, we provide a genome-wide scale landscape of DNA methylation and gene expression. The hepatic hypo-methylation feature has been identified with FLHS chickens. By integrated analysis, the results strongly suggest that increased lipid accumulation and hepatocyte rupture are central pathways that are regulated by DNA methylation in chickens with FLHS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07305-3.

Unraveling the emergence and population diversity of Listeria monocytogenes in a newly built meat facility through whole genome sequencing

The food processing environments of a newly opened meat processing facility were sampled in ten visits carried out during its first 1.5 years of activity and analyzed for the presence of Listeria monocytogenes. A total of 18 L. monocytogenes isolates were obtained from 229 samples, and their genomes were sequenced to perform comparative genomic analyses. An increase in the frequency of isolation of L. monocytogenes and in the diversity of sequence types (STs) detected was observed along time. Although the strains isolated belonged to six different STs (ST8, ST9, ST14, ST37, ST121 and ST155), ST9 was the most abundant (8 out of 18 strains). Low (0 and 2) single nucleotide polymorphism (SNP) distances were found between two pairs of ST9 strains isolated in both cases 3 months apart from the same processing room (Lm-1267 and Lm-1705, with a 2 SNPs distance in the core genome; Lm-1265 and Lm-1706, with a 0 SNPs distance), which suggests that these strains may be persistent L. monocytogenes strains in the food processing environment. Most strains showed an in silico attenuated virulence potential either through the truncation of InlA (in 67% of the isolates) or the absence of other virulence factors involved in cell adhesion or invasion. Twelve of the eighteen L. monocytogenes isolates contained a plasmid, which ranged in size from 4 to 87 Kb and harbored stress survival, in addition to heavy metals and biocides resistance determinants. Identical or highly similar plasmids were identified for various sets of L. monocytogenes ST9 isolates, which suggests the clonal expansion and persistence of plasmid-containing ST9 strains in the processing environments of the meat facility. Finally, the analysis of the L. monocytogenes genomes available in the NCBI database, and their associated metadata, evidenced that strains from ST9 are more frequently reported in Europe, linked to foods, particularly to meat and pork products, and less represented among clinical isolates than other L. monocytogenes STs. It also showed that the ST9 strains here isolated were more closely related to the European isolates, which clustered together and separated from ST9 North American isolates.

Longitudinal dynamics of gut bacteriome, mycobiome and virome after fecal microbiota transplantation in graft-versus-host disease

Fecal microbiota transplant (FMT) has emerged as a potential treatment for severe colitis associated with graft-versus-host disease (GvHD) following hematopoietic stem cell transplant. Bacterial engraftment from FMT donor to recipient has been reported, however the fate of fungi and viruses after FMT remains unclear. Here we report longitudinal dynamics of the gut bacteriome, mycobiome and virome in a teenager with GvHD after receiving four doses of FMT at weekly interval. After serial FMTs, the gut bacteriome, mycobiome and virome of the patient differ from compositions before FMT with variable temporal dynamics. Diversity of the gut bacterial community increases after each FMT. Gut fungal community initially shows expansion of several species followed by a decrease in diversity after multiple FMTs. In contrast, gut virome community varies substantially over time with a stable rise in diversity. The bacterium, Corynebacterium jeikeium, and Torque teno viruses, decrease after FMTs in parallel with an increase in the relative abundance of Caudovirales bacteriophages. Collectively, FMT may simultaneously impact on the various components of the gut microbiome with distinct effects.

FHOD formin and SRF promote post-embryonic striated muscle growth through separate pathways in C. elegans

Previous work with cultured cells has shown transcription of muscle genes by serum response factor (SRF) can be stimulated by actin polymerization driven by proteins of the formin family. However, it is not clear if endogenous formins similarly promote SRF-dependent transcription during muscle development in vivo. We tested whether formin activity promotes SRF-dependent transcription in striated muscle in the simple animal model, Caenorhabditis elegans. Our lab has shown FHOD-1 is the only formin that directly promotes sarcomere formation in the worm's striated muscle. We show here FHOD-1 and SRF homolog UNC-120 both support muscle growth and also muscle myosin II heavy chain A expression. However, while a hypomorphic unc-120 allele blunts expression of a set of striated muscle genes, these genes are largely upregulated or unchanged by absence of FHOD-1. Instead, pharmacological inhibition of the proteasome restores myosin protein levels in worms lacking FHOD-1, suggesting elevated proteolysis accounts for their myosin deficit. Interestingly, proteasome inhibition does not restore normal muscle growth to fhod-1(Δ) mutants, suggesting formin contributes to muscle growth by some alternative mechanism. Overall, we find SRF does not depend on formin to promote muscle gene transcription in a simple in vivo system.

Computational Detection of Plant RNA Editing Events

Computers are able to systematically exploit RNA-seq data allowing us to efficiently detect RNA editing sites in a genome-wide scale. This chapter introduces a very flexible computational framework for detecting RNA editing sites in plant organelles. This framework comprises three major steps: RNA-seq data processing, RNA read alignment, and RNA editing site detection. Each step is discussed in sufficient detail to be implemented by the reader. As a study case, the framework will be used with publicly available sequencing data to detect C-to-U RNA editing sites in the coding sequences of the mitochondrial genome of Nicotiana tabacum.

Intraspecies plasmid and genomic variation of Mycobacterium kubicae revealed by the complete genome sequences of two clinical isolates

Mycobacterium kubicae is 1 of nearly 200 species of nontuberculous mycobacteria (NTM), environmental micro-organisms that in some situations can infect humans and cause severe lung, skin and soft tissue infections. Although numerous studies have investigated the genetic variation among prevalent clinical NTM species, including Mycobacterium abscessus and Mycobacterium avium, many of the less common but clinically relevant NTM species, including M. kubicae, still lack complete genomes to serve as a comparative reference. Well-characterized representative genomes for each NTM species are important both for investigating the pathogenic potential of NTM, as well as for use in diagnostic methods, even for species that less frequently cause human disease. Here, we report the complete genomes of two M. kubicae strains, isolated from two unrelated patients. Hybrid short-read and long-read sequencing and assembly, using sequence reads from Illumina and Oxford Nanopore Technologies platforms, were utilized to resolve the chromosome and plasmid sequences of each isolate. The genome of NJH_MKUB1 had 5135 coding sequences (CDSs), a circular chromosome of length 5.3 Mb and two plasmids. The genome of NJH_MKUB2 had 5957 CDSs, a circular chromosome of 6.0 Mb and five plasmids. We compared our completed genomic assemblies to four recently released draft genomes of M. kubicae in order to better understand intraspecies genomic conservation and variability. We also identified genes implicated in drug resistance, virulence and persistence in the M. kubicae chromosome and plasmids. Virulence factors encoded in the genome and in the plasmids of M. kubicae provide a foundation for investigating how opportunistic environmental NTM may cause disease.

Pervasive prophage recombination occurs during evolution of spore-forming Bacilli

Phages are the main source of within-species bacterial diversity and drivers of horizontal gene transfer, but we know little about the mechanisms that drive genetic diversity of these mobile genetic elements (MGEs). Recently, we showed that a sporulation selection regime promotes evolutionary changes within SPβ prophage of Bacillus subtilis, leading to direct antagonistic interactions within the population. Herein, we reveal that under a sporulation selection regime, SPβ recombines with low copy number phi3Ts phage DNA present within the B. subtilis population. Recombination results in a new prophage occupying a different integration site, as well as the spontaneous release of virulent phage hybrids. Analysis of Bacillus sp. strains suggests that SPβ and phi3T belong to a distinct cluster of unusually large phages inserted into sporulation-related genes that are equipped with a spore-related genetic arsenal. Comparison of Bacillus sp. genomes indicates that similar diversification of SPβ-like phages takes place in nature. Our work is a stepping stone toward empirical studies on phage evolution, and understanding the eco-evolutionary relationships between bacteria and their phages. By capturing the first steps of new phage evolution, we reveal striking relationship between survival strategy of bacteria and evolution of their phages.

Insights into gene expression changes under conditions that facilitate horizontal gene transfer (mating) of a model archaeon

Horizontal gene transfer is a means by which bacteria, archaea, and eukaryotes are able to trade DNA within and between species. While there are a variety of mechanisms through which this genetic exchange can take place, one means prevalent in the archaeon Haloferax volcanii involves the transient formation of cytoplasmic bridges between cells and is referred to as mating. This process can result in the exchange of very large fragments of DNA between the participating cells. Genes governing the process of mating, including triggers to initiate mating, mechanisms of cell fusion, and DNA exchange, have yet to be characterized. We used a transcriptomic approach to gain a more detailed knowledge of how mating might transpire. By examining the differential expression of genes expressed in cells harvested from mating conditions on a filter over time and comparing them to those expressed in a shaking culture, we were able to identify genes and pathways potentially associated with mating. These analyses provide new insights into both the mechanisms and barriers of mating in Hfx. volcanii.

Mycobacteria excise DNA damage in 12- or 13-nucleotide-long oligomers by prokaryotic-type dual incisions and performs transcription-coupled repair

In nucleotide excision repair, bulky DNA lesions such as UV-induced cyclobutane pyrimidine dimers are removed from the genome by concerted dual incisions bracketing the lesion, followed by gap filling and ligation. So far, two dual-incision patterns have been discovered: the prokaryotic type, which removes the damage in 11-13-nucleotide-long oligomers, and the eukaryotic type, which removes the damage in 24-32-nucleotide-long oligomers. However, a recent study reported that the UvrC protein of Mycobacterium tuberculosis removes damage in a manner analogous to yeast and humans in a 25-mer oligonucleotide arising from incisions at 15 nt from the 3´ end and 9 nt from the 5´ end flanking the damage. To test this model, we used the in vivo excision assay and the excision repair sequencing genome-wide repair mapping method developed in our laboratory to determine the repair pattern and genome-wide repair map of Mycobacterium smegmatis We find that M. smegmatis, which possesses homologs of the Escherichia coli uvrA, uvrB, and uvrC genes, removes cyclobutane pyrimidine dimers from the genome in a manner identical to the prokaryotic pattern by incising 7 nt 5´ and 3 or 4 nt 3´ to the photoproduct, and performs transcription-coupled repair in a manner similar to E. coli.

AN1-type zinc finger protein 3 (ZFAND3) is a transcriptional regulator that drives Glioblastoma invasion

The infiltrative nature of Glioblastoma (GBM), the most aggressive primary brain tumor, critically prevents complete surgical resection and masks tumor cells behind the blood brain barrier reducing the efficacy of systemic treatment. Here, we use a genome-wide interference screen to determine invasion-essential genes and identify the AN1/A20 zinc finger domain containing protein 3 (ZFAND3) as a crucial driver of GBM invasion. Using patient-derived cellular models, we show that loss of ZFAND3 hampers the invasive capacity of GBM, whereas ZFAND3 overexpression increases motility in cells that were initially not invasive. At the mechanistic level, we find that ZFAND3 activity requires nuclear localization and integral zinc-finger domains. Our findings indicate that ZFAND3 acts within a nuclear protein complex to activate gene transcription and regulates the promoter of invasion-related genes such as COL6A2, FN1, and NRCAM. Further investigation in ZFAND3 function in GBM and other invasive cancers is warranted.

Glioblastomas (GBMs) are highly invasive brain tumours, but the underlying mechanisms of GBM invasion are unclear. Here, the authors perform an RNA interference screen and identify AN1-Type Zinc Finger protein 3 (ZFAND3) as a regulator of GBM invasion, and find that it acts through the transcriptional regulation of invasion-related genes.

Single-cell transcriptome analyses reveal heterogeneity in suspension cultures and clonal markers of CHO-K1 cells

Cell-to-cell variability in cell populations arises from a combination of intrinsic factors and extrinsic factors related to the milieu. However, the heterogeneity of high cell density suspension cultures for therapeutic protein production remains unknown. Here, we illustrate the increasing heterogeneity in the cellular transcriptome of serum-free adapted CHO K1 cells during high cell density suspension culture over time without concomitant changes in the genomic sequence. Cell cycle-dependent subpopulations and cell clusters, which typically appear in other single-cell transcriptome analyses, were not found in these suspension cultures. Our results indicate that cell division changes the intracellular microenvironment and leads to cell cycle-dependent heterogeneity. Whole mitochondrial single-cell genome sequencing showed cell-to-cell mitochondrial genome variation and heteroplasmy within cells. The mitochondrial genome sequencing method developed here is potentially useful for the validation of cell clonality. The culture time-dependent increase in cellular heterogeneity observed in this study did not show any attenuation in this increasing heterogeneity. Future advances in bioengineering such as culture upscaling, prolonged culturing, and complex culture systems will be confronted with the need to assess and control cellular heterogeneity, and the method described here may prove useful for this purpose.

Antimicrobial use and production system shape the fecal, environmental, and slurry resistomes of pig farms

BACKGROUND

The global threat of antimicrobial resistance (AMR) is a One Health problem impacted by antimicrobial use (AMU) for human and livestock applications. Extensive Iberian swine production is based on a more sustainable and eco-friendly management system, providing an excellent opportunity to evaluate how sustained differences in AMU impact the resistome, not only in the animals but also on the farm environment. Here, we evaluate the resistome footprint of an extensive pig farming system, maintained for decades, as compared to that of industrialized intensive pig farming by analyzing 105 fecal, environmental and slurry metagenomes from 38 farms.

RESULTS

Our results evidence a significantly higher abundance of antimicrobial resistance genes (ARGs) on intensive farms and a link between AMU and AMR to certain antimicrobial classes. We observed differences in the resistome across sample types, with a higher richness and dispersion of ARGs within environmental samples than on those from feces or slurry. Indeed, a deeper analysis revealed that differences among the three sample types were defined by taxa-ARGs associations. Interestingly, mobilome analyses revealed that the observed AMR differences between intensive and extensive farms could be linked to differences in the abundance of mobile genetic elements (MGEs). Thus, while there were no differences in the abundance of chromosomal-associated ARGs between intensive and extensive herds, a significantly higher abundance of integrons in the environment and plasmids, regardless of the sample type, was detected on intensive farms.

CONCLUSIONS

Overall, this study shows how AMU, production system, and sample type influence, mainly through MGEs, the profile and dispersion of ARGs in pig production. Video Abstract.

iCLIP analysis of RNA substrates of the archaeal exosome

Background

The archaeal exosome is an exoribonucleolytic multiprotein complex, which degrades single-stranded RNA in 3′ to 5′ direction phosphorolytically. In a reverse reaction, it can add A-rich tails to the 3′-end of RNA. The catalytic center of the exosome is in the aRrp41 subunit of its hexameric core. Its RNA-binding subunits aRrp4 and aDnaG confer poly(A) preference to the complex. The archaeal exosome was intensely characterized in vitro, but still little is known about its interaction with natural substrates in the cell, particularly because analysis of the transcriptome-wide interaction of an exoribonuclease with RNA is challenging.

Results

To determine binding sites of the exosome to RNA on a global scale, we performed individual-nucleotide resolution UV crosslinking and immunoprecipitation (iCLIP) analysis with antibodies directed against aRrp4 and aRrp41 of the chrenarchaeon Sulfolobus solfataricus. A relatively high proportion (17–19%) of the obtained cDNA reads could not be mapped to the genome. Instead, they corresponded to adenine-rich RNA tails, which are post-transcriptionally synthesized by the exosome, and to circular RNAs (circRNAs). We identified novel circRNAs corresponding to 5′ parts of two homologous, transposase-related mRNAs. To detect preferred substrates of the exosome, the iCLIP reads were compared to the transcript abundance using RNA-Seq data. Among the strongly enriched exosome substrates were RNAs antisense to tRNAs, overlapping 3′-UTRs and RNAs containing poly(A) stretches. The majority of the read counts and crosslink sites mapped in mRNAs. Furthermore, unexpected crosslink sites clustering at 5′-ends of RNAs was detected.

Conclusions

In this study, RNA targets of an exoribonuclease were analyzed by iCLIP. The data documents the role of the archaeal exosome as an exoribonuclease and RNA-tailing enzyme interacting with all RNA classes, and underlines its role in mRNA turnover, which is important for adaptation of prokaryotic cells to changing environmental conditions. The clustering of crosslink sites near 5′-ends of genes suggests simultaneous binding of both RNA ends by the S. solfataricus exosome. This may serve to prevent translation of mRNAs dedicated to degradation in 3′-5′ direction.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07200-x.

The Polycomb group protein Ring1 regulates dorsoventral patterning of the mouse telencephalon

Dorsal-ventral patterning of the mammalian telencephalon is fundamental to the formation of distinct functional regions including the neocortex and ganglionic eminence. While Bone morphogenetic protein (BMP), Wnt, and Sonic hedgehog (Shh) signaling are known to determine regional identity along the dorsoventral axis, how the region-specific expression of these morphogens is established remains unclear. Here we show that the Polycomb group (PcG) protein Ring1 contributes to the ventralization of the mouse telencephalon. Deletion of Ring1b or both Ring1a and Ring1b in neuroepithelial cells induces ectopic expression of dorsal genes, including those for BMP and Wnt ligands, as well as attenuated expression of the gene for Shh, a key morphogen for ventralization, in the ventral telencephalon. We observe PcG protein–mediated trimethylation of histone 3 at lysine-27 and binding of Ring1B at BMP and Wnt ligand genes specifically in the ventral region. Furthermore, forced activation of BMP or Wnt signaling represses Shh expression. Our results thus indicate that PcG proteins suppress BMP and Wnt signaling in a region-specific manner and thereby allow proper Shh expression and development of the ventral telencephalon.

NCOMMS-19-38235B Dorsal-ventral patterning of the mammalian telencephalon is fundamental to the formation of distinct functional regions. Here, the authors find that PcG proteins suppress BMP and Wnt signaling in a region-specific manner, allowing for proper Shh expression and development of the ventral telencephalon.

Arioc: High-concurrency short-read alignment on multiple GPUs

In large DNA sequence repositories, archival data storage is often coupled with computers that provide 40 or more CPU threads and multiple GPU (general-purpose graphics processing unit) devices. This presents an opportunity for DNA sequence alignment software to exploit high-concurrency hardware to generate short-read alignments at high speed. Arioc, a GPU-accelerated short-read aligner, can compute WGS (whole-genome sequencing) alignments ten times faster than comparable CPU-only alignment software. When two or more GPUs are available, Arioc's speed increases proportionately because the software executes concurrently on each available GPU device. We have adapted Arioc to recent multi-GPU hardware architectures that support high-bandwidth peer-to-peer memory accesses among multiple GPUs. By modifying Arioc's implementation to exploit this GPU memory architecture we obtained a further 1.8x-2.9x increase in overall alignment speeds. With this additional acceleration, Arioc computes two million short-read alignments per second in a four-GPU system; it can align the reads from a human WGS sequencer run–over 500 million 150nt paired-end reads–in less than 15 minutes. As WGS data accumulates exponentially and high-concurrency computational resources become widespread, Arioc addresses a growing need for timely computation in the short-read data analysis toolchain.

Integrative Omics Analyses Reveal Epigenetic Memory in Diabetic Renal Cells Regulating Genes Associated With Kidney Dysfunction

Diabetic kidney disease (DKD) is a major complication of diabetes and the leading cause of end-stage renal failure. Epigenetics has been associated with metabolic memory in which prior periods of hyperglycemia enhance the future risk of developing DKD despite subsequent glycemic control. To understand the mechanistic role of such epigenetic memory in human DKD and to identify new therapeutic targets, we profiled gene expression, DNA methylation, and chromatin accessibility in kidney proximal tubule epithelial cells (PTECs) derived from subjects with and without type 2 diabetes (T2D). T2D-PTECs displayed persistent gene expression and epigenetic changes with and without transforming growth factor-β1 treatment, even after culturing in vitro under similar conditions as nondiabetic PTECs, signified by deregulation of fibrotic and transport-associated genes (TAGs). Motif analysis of differential DNA methylation and chromatin accessibility regions associated with genes differentially regulated in T2D revealed enrichment for SMAD3, HNF4A, and CTCF transcription factor binding sites. Furthermore, the downregulation of several TAGs in T2D (including CLDN10, CLDN14, CLDN16, SLC16A2, and SLC16A5) was associated with promoter hypermethylation, decreased chromatin accessibility, and reduced enrichment of HNF4A, histone H3-lysine-27-acetylation, and CTCF. Together, these integrative analyses reveal epigenetic memory underlying the deregulation of key target genes in T2D-PTECs that may contribute to sustained renal dysfunction in DKD.

MicroRNA signatures of perioperative myocardial injury after elective noncardiac surgery: a prospective observational mechanistic cohort study

BACKGROUND

Elevated plasma or serum troponin, indicating perioperative myocardial injury (PMI), is common after noncardiac surgery. However, underlying mechanisms remain unclear. Acute coronary syndrome (ACS) is associated with the early appearance of circulating microRNAs, which regulate post-translational gene expression. We hypothesised that if PMI and ACS share pathophysiological mechanisms, common microRNA signatures should be evident.

METHODS

We performed a nested case control study of samples obtained before and after noncardiac surgery from patients enrolled in two prospective observational studies of PMI (postoperative troponin I/T>99th centile). In cohort one, serum microRNAs were compared between patients with or without PMI, matched for age, gender, and comorbidity. Real-time polymerase chain reaction quantified (qRT-PCR) relative microRNA expression (cycle quantification [Cq] threshold <37) before and after surgery for microRNA signatures associated with ACS, blinded to PMI. In cohort two, we analysed (EdgeR) microRNA from plasma extracellular vesicles using next-generation sequencing (Illumina HiSeq 500). microRNA-messenger RNA-function pathway analysis was performed (DIANA miRPath v3.0/TopGO).

RESULTS

MicroRNAs were detectable in all 59 patients (median age 67 yr [61-75]; 42% male), who had similar clinical characteristics independent of developing PMI. In cohort one, serum microRNA expression increased after surgery (mean fold-change) hsa-miR-1-3p: 3.99 (95% confidence interval [CI: 1.95-8.19]; hsa-miR-133-3p: 5.67 [95% CI: 2.94-10.91]; P<0.001). These changes were not associated with PMI. Bioinformatic analysis of differentially expressed microRNAs from cohorts one (n=48) and two (n=11) identified pathways associated with adrenergic stress and calcium dysregulation, rather than ischaemia.

CONCLUSIONS

Circulating microRNAs associated with cardiac ischaemia were universally elevated in patients after surgery, independent of development of myocardial injury.

Dynamically evolving novel overlapping gene as a factor in the SARS-CoV-2 pandemic

Understanding the emergence of novel viruses requires an accurate and comprehensive annotation of their genomes. Overlapping genes (OLGs) are common in viruses and have been associated with pandemics but are still widely overlooked. We identify and characterize ORF3d, a novel OLG in SARS-CoV-2 that is also present in Guangxi pangolin-CoVs but not other closely related pangolin-CoVs or bat-CoVs. We then document evidence of ORF3d translation, characterize its protein sequence, and conduct an evolutionary analysis at three levels: between taxa (21 members of Severe acute respiratory syndrome-related coronavirus), between human hosts (3978 SARS-CoV-2 consensus sequences), and within human hosts (401 deeply sequenced SARS-CoV-2 samples). ORF3d has been independently identified and shown to elicit a strong antibody response in COVID-19 patients. However, it has been misclassified as the unrelated gene ORF3b, leading to confusion. Our results liken ORF3d to other accessory genes in emerging viruses and highlight the importance of OLGs.

Lifestyle adaptations of Rhizobium from rhizosphere to symbiosis

By analyzing successive lifestyle stages of a model Rhizobium-legume symbiosis using mariner-based transposon insertion sequencing (INSeq), we have defined the genes required for rhizosphere growth, root colonization, bacterial infection, N₂-fixing bacteroids, and release from legume (pea) nodules. While only 27 genes are annotated as nif and fix in Rhizobium leguminosarum, we show 603 genetic regions (593 genes, 5 transfer RNAs, and 5 RNA features) are required for the competitive ability to nodulate pea and fix N₂ Of these, 146 are common to rhizosphere growth through to bacteroids. This large number of genes, defined as rhizosphere-progressive, highlights how critical successful competition in the rhizosphere is to subsequent infection and nodulation. As expected, there is also a large group (211) specific for nodule bacteria and bacteroid function. Nodule infection and bacteroid formation require genes for motility, cell envelope restructuring, nodulation signaling, N₂ fixation, and metabolic adaptation. Metabolic adaptation includes urea, erythritol and aldehyde metabolism, glycogen synthesis, dicarboxylate metabolism, and glutamine synthesis (GlnII). There are 17 separate lifestyle adaptations specific to rhizosphere growth and 23 to root colonization, distinct from infection and nodule formation. These results dramatically highlight the importance of competition at multiple stages of a Rhizobium-legume symbiosis.

A Metabolic Pathway for Activation of Dietary Glucosinolates by a Human Gut Symbiont

Consumption of glucosinolates, pro-drug-like metabolites abundant in Brassica vegetables, has been associated with decreased risk of certain cancers. Gut microbiota have the ability to metabolize glucosinolates, generating chemopreventive isothiocyanates. Here, we identify a genetic and biochemical basis for activation of glucosinolates to isothiocyanates by Bacteroides thetaiotaomicron, a prominent gut commensal species. Using a genome-wide transposon insertion screen, we identified an operon required for glucosinolate metabolism in B. thetaiotaomicron. Expression of BT2159-BT2156 in a non-metabolizing relative, Bacteroides fragilis, resulted in gain of glucosinolate metabolism. We show that isothiocyanate formation requires the action of BT2158 and either BT2156 or BT2157 in vitro. Monocolonization of mice with mutant BtΔ2157 showed reduced isothiocyanate production in the gastrointestinal tract. These data provide insight into the mechanisms by which a common gut bacterium processes an important dietary nutrient.

Within-individual phenotypic plasticity in flowers fosters pollination niche shift

Phenotypic plasticity, the ability of a genotype of producing different phenotypes when exposed to different environments, may impact ecological interactions. We study here how within-individual plasticity in Moricandia arvensis flowers modifies its pollination niche. During spring, this plant produces large, cross-shaped, UV-reflecting lilac flowers attracting mostly long-tongued large bees. However, unlike most co-occurring species, M. arvensis keeps flowering during the hot, dry summer due to its plasticity in key vegetative traits. Changes in temperature and photoperiod in summer trigger changes in gene expression and the production of small, rounded, UV-absorbing white flowers that attract a different assemblage of generalist pollinators. This shift in pollination niche potentially allows successful reproduction in harsh conditions, facilitating M. arvensis to face anthropogenic perturbations and climate change.

Enhancer evolutionary co-option through shared chromatin accessibility input

Form diversity is fueled by changes in the expression of genes that build organisms. New expression often results from the emergence of new DNA switches, known as transcriptional enhancers. Many enhancers are thought to appear through the recycling of older enhancers, a process called evolutionary co-option. Enhancer co-option is difficult to assess, and the molecular mechanisms explaining its prevalence are elusive. Using state-of-the-art quantification and analyses, we reveal that the sequences of an ancestral and a derived enhancer overlap extensively. They contain specific binding sites for regulators imparting spatial activities. We found that the two enhancers also share a site facilitating access to chromatin in a region where they overlap.

The diversity of forms in multicellular organisms originates largely from the spatial redeployment of developmental genes [S. B. Carroll, Cell 134, 25–36 (2008)]. Several scenarios can explain the emergence of cis-regulatory elements that govern novel aspects of a gene expression pattern [M. Rebeiz, M. Tsiantis, Curr. Opin. Genet. Dev. 45, 115–123 (2017)]. One scenario, enhancer co-option, holds that a DNA sequence producing an ancestral regulatory activity also becomes the template for a new regulatory activity, sharing regulatory information. While enhancer co-option might fuel morphological diversification, it has rarely been documented [W. J. Glassford et al., Dev. Cell 34, 520–531 (2015)]. Moreover, if two regulatory activities are borne from the same sequence, their modularity, considered a defining feature of enhancers [J. Banerji, L. Olson, W. Schaffner, Cell 33, 729–740 (1983)], might be affected by pleiotropy. Sequence overlap may thereby play a determinant role in enhancer function and evolution. Here, we investigated this problem with two regulatory activities of the Drosophila gene yellow, the novel spot enhancer and the ancestral wing blade enhancer. We used precise and comprehensive quantification of each activity in Drosophila wings to systematically map their sequences along the locus. We show that the spot enhancer has co-opted the sequences of the wing blade enhancer. We also identified a pleiotropic site necessary for DNA accessibility of a shared regulatory region. While the evolutionary steps leading to the derived activity are still unknown, such pleiotropy suggests that enhancer accessibility could be one of the molecular mechanisms seeding evolutionary co-option.

Phlox Species Show Quantitative and Qualitative Resistance to a Population of Powdery Mildew Isolates from the Eastern United States

Ornamental plants in the genus Phlox are extensively planted in landscapes and home gardens around the world. A major limitation to a more widespread use of these plants is their susceptibility to powdery mildew (PM). In this study, we used multilocus sequence typing (MLST) analysis to gain insights into the population diversity of 32 Phlox PM pathogen (Golovinomyces magnicellulatus and Podosphaera sp.) isolates collected from the eastern United States and relate it to the ability to overcome host resistance. Low genetic diversity and a lack of structure were found within our population. Whole genome comparison of two isolates was used to support low genetic diversity evidence found with the MLST analysis. Recombination was suggested by the incongruences observed in the six phylogenetic trees generated from the housekeeping genes TEF-1α, CSI, ITS, IGS, H3, and TUB. Contrasting with low genetic diversity, we found high phenotypic diversity when using 10 of the 32 isolates to evaluate host resistance in four different Phlox species (P. paniculata 'Dunbar Creek', P. amoena OPGC 3598, P. glaberrima OPGC 3594, and P. subulata OPGC 4185) using in vitro bioassays. We observed quantitative and qualitative resistance in all Phlox species and a consistent low disease severity in our control, P. paniculata 'Dunbar Creek'. Taken together, the results generated in this study constitute a robust screening of popular Phlox germplasm that can be incorporated into breeding programs for PM resistance and provides significant information on the evolution of PM pathogens.

Chd4 choreographs self-antigen expression for central immune tolerance

Autoreactive T cells are eliminated in the thymus to prevent autoimmunity by promiscuous expression of tissue-restricted self-antigens in medullary thymic epithelial cells. This expression is dependent on the transcription factor Fezf2, as well as the transcriptional regulator Aire, but the entire picture of the transcriptional program has been obscure. Here, we found that the chromatin remodeler Chd4, also called Mi-2β, plays a key role in the self-antigen expression in medullary thymic epithelial cells. To maximize the diversity of self-antigen expression, Fezf2 and Aire utilized completely distinct transcriptional mechanisms, both of which were under the control of Chd4. Chd4 organized the promoter regions of Fezf2-dependent genes, while contributing to the Aire-mediated induction of self-antigens via super-enhancers. Mice deficient in Chd4 specifically in thymic epithelial cells exhibited autoimmune phenotypes, including T cell infiltration. Thus, Chd4 plays a critical role in integrating Fezf2- and Aire-mediated gene induction to establish central immune tolerance.

A Transposon Story: From TE Content to TE Dynamic Invasion of Drosophila Genomes Using the Single-Molecule Sequencing Technology from Oxford Nanopore

Transposable elements (TEs) are the main components of genomes. However, due to their repetitive nature, they are very difficult to study using data obtained with short-read sequencing technologies. Here, we describe an efficient pipeline to accurately recover TE insertion (TEI) sites and sequences from long reads obtained by Oxford Nanopore Technology (ONT) sequencing. With this pipeline, we could precisely describe the landscapes of the most recent TEIs in wild-type strains of Drosophila melanogaster and Drosophila simulans. Their comparison suggests that this subset of TE sequences is more similar than previously thought in these two species. The chromosome assemblies obtained using this pipeline also allowed recovering piRNA cluster sequences, which was impossible using short-read sequencing. Finally, we used our pipeline to analyze ONT sequencing data from a D. melanogaster unstable line in which LTR transposition was derepressed for 73 successive generations. We could rely on single reads to identify new insertions with intact target site duplications. Moreover, the detailed analysis of TEIs in the wild-type strains and the unstable line did not support the trap model claiming that piRNA clusters are hotspots of TE insertions.

Heterogeneity of the rice microbial community of the Chinese centuries-old Honghe Hani rice terraces system

The Honghe Hani rice terraces system (HHRTS) is a traditional rice cultivation system where Hani people cultivate remarkably diverse rice varieties. Recent introductions of modern rice varieties to the HHRTS have significantly increased the severity of rice diseases within the terraces. Here, we determine the impacts of these recent introductions on the composition of the rice-associated microbial communities. We confirm that the HHRTS contains a range of both traditional HHRTS landraces and introduced modern rice varieties and find differences between the microbial communities of these two groups. However, this introduction of modern rice varieties has not strongly impacted the overall diversity of the HHRTS rice microbial community. Furthermore, we find that the rice varieties (i.e. groups of closely related genotypes) have significantly structured the rice microbial community composition (accounting for 15%-22% of the variance) and that the core microbial community of HHRTS rice plants represents less than 3.3% of all the microbial taxa identified. Collectively, our study suggests a highly diverse HHRTS rice holobiont (host with its associated microbes) where the diversity of rice hosts mirrors the diversity of their microbial communities. Further studies will be needed to better determine how such changes might impact the sustainability of the HHRTS.

Analysis of H3K4me3 and H3K27me3 bivalent promotors in HER2+ breast cancer cell lines reveals variations depending on estrogen receptor status and significantly correlates with gene expression

BACKGROUND

The role of histone modifications is poorly characterized in breast cancer, especially within the major subtypes. While epigenetic modifications may enhance the adaptability of a cell to both therapy and the surrounding environment, the mechanisms by which this is accomplished remains unclear. In this study we focus on the HER2 subtype and investigate two histone trimethylations that occur on the histone 3; the trimethylation located at lysine 4 (H3K4me3) found in active promoters and the trimethylation located at lysine 27 (H3K27me3) that correlates with gene repression. A bivalency state is the result of the co-presence of these two marks at the same promoter.

METHODS

In this study we investigated the relationship between these histone modifications in promoter regions and their proximal gene expression in HER2+ breast cancer cell lines. In addition, we assessed these patterns with respect to the presence or absence of the estrogen receptor (ER). To do this, we utilized ChIP-seq and matching RNA-seq from publicly available data for the AU565, SKBR3, MB361 and UACC812 cell lines. In order to visualize these relationships, we used KEGG pathway enrichment analysis, and Kaplan-Meyer plots.

RESULTS

We found that the correlation between the three types of promoter trimethylation statuses (H3K4me3, H3K27me3 or both) and the expression of the proximal genes was highly significant overall, while roughly a third of all genes are regulated by this phenomenon. We also show that there are several pathways related to cancer progression and invasion that are associated with the bivalent status of the gene promoters, and that there are specific differences between ER+ and ER- HER2+ breast cancer cell lines. These specific differences that are differentially trimethylated are also shown to be differentially expressed in patient samples. One of these genes, HIF1AN, significantly correlates with patient outcome.

CONCLUSIONS

This study highlights the importance of looking at epigenetic markings at a subtype specific level by characterizing the relationship between the bivalent promoters and gene expression. This provides a deeper insight into a mechanism that could lead to future targets for treatment and prognosis, along with oncogenesis and response to therapy of HER2+ breast cancer patients.

Different Interspecies Electron Transfer Patterns during Mesophilic and Thermophilic Syntrophic Propionate Degradation in Chemostats

Propionate is one of the major intermediates in anaerobic digestion of organic waste to CO₂ and CH₄. In methanogenic environments, propionate is degraded through a mutualistic interaction between symbiotic propionate oxidizers and methanogens. Although temperature heavily influences the microbial ecology and performance of methanogenic processes, its effect on syntrophic interaction during propionate degradation remains poorly understood. In this study, metagenomics and metatranscriptomics were employed to compare mesophilic and thermophilic propionate degradation communities. Mesophilic propionate degradation involved multiple syntrophic organisms (Syntrophobacter, Smithella, and Syntrophomonas), pathways, interactions, and preference toward formate-based electron transfer to methanogenic partners (i.e., Methanoculleus). In thermophilic propionate degradation, one syntrophic organism predominated (Pelotomaculum), interspecies H₂ transfer played a major role, and phylogenetically and metabolically diverse H₂-oxidizing methanogens were present (i.e., Methanoculleus, Methanothermobacter, and Methanomassiliicoccus). This study showed that microbial interactions, metabolic pathways, and niche diversity are distinct between mesophilic and thermophilic microbial communities responsible for syntrophic propionate degradation.

IRF4 instructs effector Treg differentiation and immune suppression in human cancer

The molecular mechanisms responsible for the high immunosuppressive capacity of CD4+ Tregs in tumors are not well known. High-dimensional single-cell profiling of T cells from chemotherapy-naive individuals with non-small-cell lung cancer identified the transcription factor IRF4 as specifically expressed by a subset of intratumoral CD4+ effector Tregs with superior suppressive activity. In contrast to the IRF4- counterparts, IRF4+ Tregs expressed a vast array of suppressive molecules, and their presence correlated with multiple exhausted subpopulations of T cells. Integration of transcriptomic and epigenomic data revealed that IRF4, either alone or in combination with its partner BATF, directly controlled a molecular program responsible for immunosuppression in tumors. Accordingly, deletion of Irf4 exclusively in Tregs resulted in delayed tumor growth in mice while the abundance of IRF4+ Tregs correlated with poor prognosis in patients with multiple human cancers. Thus, a common mechanism underlies immunosuppression in the tumor microenvironment irrespective of the tumor type.

A new strain group of common carp: The genetic differences and admixture events between Cyprinus carpio breeds

Common carp (Cyprinus carpio) has an outstanding economic importance in freshwater aquaculture due to its high adaptive capacity to both food and environment. In fact, it is the third most farmed fish species worldwide according to the Food and Agriculture Organization. More than four million tons of common carp are produced annually in aquaculture, and more than a hundred thousand tons are caught from the wild. Historically, the common carp was also the first fish species to be domesticated in ancient China, and now, there is a huge variety of domestic carp strains worldwide. In the present study, we used double digestion restriction site-associated DNA sequencing to genotype several European common carp strains and showed that they are divided into two distinct groups. One of them includes central European common carp strains as well as Ponto-Caspian wild common carp populations, whereas the other group contains several common carp strains that originated in the Soviet Union, mostly as cold-resistant strains. We believe that breeding with wild Amur carp and subsequent selection of the hybrids for resistance to adverse environmental conditions was the attribute of the second group. We assessed the contribution of wild Amur carp inheritance to the common carp strains and discovered discriminating genes, which differed in allele frequencies between groups. Taken together, our results improve our current understanding of the genetic variability of common carp, namely the structure of natural and artificial carp populations, and the contribution of wild carp traits to domestic strains.

Gene expression trajectories during male and female reproductive development in balsam poplar (Populus balsamifera L.)

Plant reproductive development from the first appearance of reproductively committed axes through to floral maturation requires massive and rapid remarshalling of gene expression. In dioecious species such as poplar this is further complicated by divergent male and female developmental programs. We used seven time points in male and female balsam poplar (Populus balsamifera L.) buds and catkins representing the full annual flowering cycle, to elucidate the effects of time and sex on gene expression during reproductive development. Time (developmental stage) is dominant in patterning gene expression with the effect of sex nested within this. Here, we find (1) evidence for five successive waves of alterations to the chromatin landscape which may be important in setting the overall reproductive trajectory, regardless of sex. (2) Each individual developmental stage is further characterized by marked sex-differential gene expression. (3) Consistent sexually differentiated gene expression regardless of developmental stage reveal candidates for high-level regulators of sex and include the female-specific poplar ARR17 homologue. There is also consistent male-biased expression of the MADS-box genes PISTILLATA and APETALA3. Our work provides insights into expression trajectories shaping reproductive development, its potential underlying mechanisms, and sex-specific translation of the genome information into reproductive structures in balsam poplar.

Drosophila histone locus body assembly and function involves multiple interactions

The histone locus body (HLB) assembles at replication-dependent (RD) histone loci and concentrates factors required for RD histone mRNA biosynthesis. The Drosophila melanogaster genome has a single locus comprised of ∼100 copies of a tandemly arrayed 5-kB repeat unit containing one copy of each of the 5 RD histone genes. To determine sequence elements required for D. melanogaster HLB formation and histone gene expression, we used transgenic gene arrays containing 12 copies of the histone repeat unit that functionally complement loss of the ∼200 endogenous RD histone genes. A 12x histone gene array in which all H3-H4 promoters were replaced with H2a-H2b promoters (12x^PR) does not form an HLB or express high levels of RD histone mRNA in the presence of the endogenous histone genes. In contrast, this same transgenic array is active in HLB assembly and RD histone gene expression in the absence of the endogenous RD histone genes and rescues the lethality caused by homozygous deletion of the RD histone locus. The HLB formed in the absence of endogenous RD histone genes on the mutant 12x array contains all known factors present in the wild-type HLB including CLAMP, which normally binds to GAGA repeats in the H3-H4 promoter. These data suggest that multiple protein–protein and/or protein–DNA interactions contribute to HLB formation, and that the large number of endogenous RD histone gene copies sequester available factor(s) from attenuated transgenic arrays, thereby preventing HLB formation and gene expression on these arrays.

Chromatin accessibility dynamics and single cell RNA-Seq reveal new regulators of regeneration in neural progenitors

Vertebrate appendage regeneration requires precisely coordinated remodeling of the transcriptional landscape to enable the growth and differentiation of new tissue, a process executed over multiple days and across dozens of cell types. The heterogeneity of tissues and temporally-sensitive fate decisions involved has made it difficult to articulate the gene regulatory programs enabling regeneration of individual cell types. To better understand how a regenerative program is fulfilled by neural progenitor cells (NPCs) of the spinal cord, we analyzed pax6-expressing NPCs isolated from regenerating Xenopus tropicalis tails. By intersecting chromatin accessibility data with single-cell transcriptomics, we find that NPCs place an early priority on neuronal differentiation. Late in regeneration, the priority returns to proliferation. Our analyses identify Pbx3 and Meis1 as critical regulators of tail regeneration and axon organization. Overall, we use transcriptional regulatory dynamics to present a new model for cell fate decisions and their regulators in NPCs during regeneration.

Genomic Evidence for Formate Metabolism by Chloroflexi as the Key to Unlocking Deep Carbon in Lost City Microbial Ecosystems

Primitive forms of life may have originated around hydrothermal vents at the bottom of the ancient ocean. The Lost City hydrothermal vent field, fueled by just rock and water, provides an analog for not only primitive ecosystems but also potential extraterrestrial rock-powered ecosystems. The microscopic life covering the towering chimney structures at the Lost City has been previously documented, yet little is known about the carbon cycling in this ecosystem. These results provide a better understanding of how carbon from the deep subsurface can fuel rich microbial ecosystems on the seafloor.

The Lost City hydrothermal field on the Mid-Atlantic Ridge supports dense microbial life on the lofty calcium carbonate chimney structures. The vent field is fueled by chemical reactions between the ultramafic rock under the chimneys and ambient seawater. These serpentinization reactions provide reducing power (as hydrogen gas) and organic compounds that can serve as microbial food; the most abundant of these are methane and formate. Previous studies have characterized the interior of the chimneys as a single-species biofilm inhabited by the Lost City Methanosarcinales, but they also indicated that this methanogen is unable to metabolize formate. The new metagenomic results presented here indicate that carbon cycling in these Lost City chimney biofilms could depend on the metabolism of formate by Chloroflexi populations. Additionally, we present evidence for metabolically diverse, formate-utilizing Sulfurovum populations and new genomic and phylogenetic insights into the unique Lost City Methanosarcinales.

IMPORTANCE Primitive forms of life may have originated around hydrothermal vents at the bottom of the ancient ocean. The Lost City hydrothermal vent field, fueled by just rock and water, provides an analog for not only primitive ecosystems but also potential extraterrestrial rock-powered ecosystems. The microscopic life covering the towering chimney structures at the Lost City has been previously documented, yet little is known about the carbon cycling in this ecosystem. These results provide a better understanding of how carbon from the deep subsurface can fuel rich microbial ecosystems on the seafloor.

An MXD1-derived repressor peptide identifies noncoding mediators of MYC-driven cell proliferation

MYC controls the transcription of large numbers of long noncoding RNAs (lncRNAs). Since MYC is a ubiquitous oncoprotein, some of these lncRNAs probably play a significant role in cancer. We applied CRISPR interference (CRISPRi) to the identification of MYC-regulated lncRNAs that are required for MYC-driven cell proliferation in the P493-6 and RAMOS human lymphoid cell lines. We identified 320 noncoding loci that play positive roles in cell growth. Transcriptional repression of any one of these lncRNAs reduces the proliferative capacity of the cells. Selected hits were validated by RT-qPCR and in CRISPRi competition assays with individual GFP-expressing sgRNA constructs. We also showed binding of MYC to the promoter of two candidate genes by chromatin immunoprecipitation. In the course of our studies, we discovered that the repressor domain SID (SIN3-interacting domain) derived from the MXD1 protein is highly effective in P493-6 and RAMOS cells in terms of the number of guides depleted in library screening and the extent of the induced transcriptional repression. In the cell lines used, SID is superior to the KRAB repressor domain, which serves routinely as a transcriptional repressor domain in CRISPRi. The SID transcriptional repressor domain is effective as a fusion to the MS2 aptamer binding protein MCP, allowing the construction of a doxycycline-regulatable CRISPRi system that allows controlled repression of targeted genes and will facilitate the functional analysis of growth-promoting lncRNAs.

A Nearly Complete Genome of Ciona intestinalis Type A (C. robusta) Reveals the Contribution of Inversion to Chromosomal Evolution in the Genus Ciona

Since its initial publication in 2002, the genome of Ciona intestinalis type A (Ciona robusta), the first genome sequence of an invertebrate chordate, has provided a valuable resource for a wide range of biological studies, including developmental biology, evolutionary biology, and neuroscience. The genome assembly was updated in 2008, and it included 68% of the sequence information in 14 pairs of chromosomes. However, a more contiguous genome is required for analyses of higher order genomic structure and of chromosomal evolution. Here, we provide a new genome assembly for an inbred line of this animal, constructed with short and long sequencing reads and Hi-C data. In this latest assembly, over 95% of the 123 Mb of sequence data was included in the chromosomes. Short sequencing reads predicted a genome size of 114-120 Mb; therefore, it is likely that the current assembly contains almost the entire genome, although this estimate of genome size was smaller than previous estimates. Remapping of the Hi-C data onto the new assembly revealed a large inversion in the genome of the inbred line. Moreover, a comparison of this genome assembly with that of Ciona savignyi, a different species in the same genus, revealed many chromosomal inversions between these two Ciona species, suggesting that such inversions have occurred frequently and have contributed to chromosomal evolution of Ciona species. Thus, the present assembly greatly improves an essential resource for genome-wide studies of ascidians.

The genomic survey of Tc1-like elements in the silkworm microsporidia Nosema bombycis

BACKGROUND

Microsporidia Nosema bombycis is the destructive pathogen in the production of sericulture. The Tc1/mariner elements belong to important component of DNA transposon.

METHODS

The genomic data of N. bombycis and related Nosema species were screened to identify the Tc1-like elements and analyzed the phylogenetic relationship, based on bioinformational analysis. High-throughput data of transcriptomes and small RNAs were used to evaluate the expressed level and potential rasiRNAs for the Tc1-like elements of N. bombycis.

RESULTS

Twelve complete Tc1-like elements belonging to DD34,E clade is confirmed in the whole genome of N. bombycis, and divided into two branches. Six of them are sole in N. bombycis and thereby would be the molecular marker to differentiate this species from others Nosema spp. Most of the elements have the transcriptional active and are the source of sRNAs.

CONCLUSION

Abundant Tc1-like elements in N. bombycis reflect the expansion of transposons for this genomic characters, comparing with others Nosema spp. The finding of distribution, phylogeny and potential functional activity for Tc1Nbs in N. bombycis will help understanding the role of the DNA transposon in genomic evolution of microsporidia.

Transcriptome Analysis of Apple Leaves Infected by the Rust Fungus Gymnosporangium yamadae at Two Sporulation Stages

Apple rust disease caused by Gymnosporangium yamadae is one of the major threats to apple orchards. In this study, dual RNA-seq analysis was conducted to simultaneously monitor gene expression profiles of G. yamadae and infected apple leaves during the formation of rust spermogonia and aecia. The molecular mechanisms underlying this compatible interaction at 10 and 30 days postinoculation (dpi) indicate a significant reaction from the host plant and comprise detoxication pathways at the earliest stage and the induction of secondary metabolism pathways at 30 dpi. Such host reactions have been previously reported in other rust pathosystems and may represent a general reaction to rust infection. G. yamadae transcript profiling indicates a conserved genetic program in spermogonia and aecia that is shared with other rust fungi, whereas secretome prediction reveals the presence of specific secreted candidate effector proteins expressed during apple infection. Unexpectedly, the survey of fungal unigenes in the transcriptome assemblies of inoculated and mock-inoculated apple leaves reveals that G. yamadae infection may modify the fungal community composition in the apple phyllosphere at 30 dpi. Collectively, our results provide novel insights into the compatible apple-G. yamadae interaction and advance the knowledge of this heteroecious demicyclic rust fungus.

Xenosiderophore Utilization Promotes Bacteroides thetaiotaomicron Resilience during Colitis

During short-lived perturbations, such as inflammation, the gut microbiota exhibits resilience and reverts to its original configuration. Although microbial access to the micronutrient iron is decreased during colitis, pathogens can scavenge iron by using siderophores. How commensal bacteria acquire iron during gut inflammation is incompletely understood. Curiously, the human commensal Bacteroides thetaiotaomicron does not produce siderophores but grows under iron-limiting conditions using enterobacterial siderophores. Using RNA-seq, we identify B. thetaiotaomicron genes that were upregulated during Salmonella-induced gut inflammation and were predicted to be involved in iron uptake. Mutants in the xusABC locus (BT2063-2065) were defective for xenosiderophore-mediated iron uptake in vitro. In the normal mouse gut, the XusABC system was dispensable, while a xusA mutant colonized poorly during colitis. This work identifies xenosiderophore utilization as a critical mechanism for B. thetaiotaomicron to sustain colonization during inflammation and suggests a mechanism of how interphylum iron metabolism contributes to gut microbiota resilience.

Host-Specific Evolutionary and Transmission Dynamics Shape the Functional Diversification of Staphylococcus epidermidis in Human Skin

Metagenomic inferences of bacterial strain diversity and infectious disease transmission studies largely assume a dominant, within-individual haplotype. We hypothesize that within-individual bacterial population diversity is critical for homeostasis of a healthy microbiome and infection risk. We characterized the evolutionary trajectory and functional distribution of Staphylococcus epidermidis-a keystone skin microbe and opportunistic pathogen. Analyzing 1,482 S. epidermidis genomes from 5 healthy individuals, we found that skin S. epidermidis isolates coalesce into multiple founder lineages rather than a single colonizer. Transmission events, natural selection, and pervasive horizontal gene transfer result in population admixture within skin sites and dissemination of antibiotic resistance genes within-individual. We provide experimental evidence for how admixture can modulate virulence and metabolism. Leveraging data on the contextual microbiome, we assess how interspecies interactions can shape genetic diversity and mobile gene elements. Our study provides insights into how within-individual evolution of human skin microbes shapes their functional diversification.

Cabbage looper (Trichoplusia ni Hübner) labial glands contain unique bacterial flora in contrast with their alimentary canal, mandibular glands, and Malpighian tubules

In recent years, several studies have examined the gut microbiome of lepidopteran larvae and how factors such as host plant affect it, and in turn, how gut bacteria affect host plant responses to herbivory. In addition, other studies have detailed how secretions of the labial (salivary) glands can alter host plant defense responses. We examined the gut microbiome of the cabbage looper (Trichoplusia ni) feeding on collards (Brassica oleracea) and separately analyzed the microbiomes of various organs that open directly into the alimentary canal, including the labial glands, mandibular glands, and the Malpighian tubules. In this study, the gut microbiome of T. ni was found to be generally consistent with those of other lepidopteran larvae in prior studies. The greatest diversity of bacteria appeared in the Firmicutes, Actinobacteria, Proteobacteria, and Bacteriodetes. Well‐represented genera included Staphylococcus, Streptococcus, Corynebacterium, Pseudomonas, Diaphorobacter, Methylobacterium, Flavobacterium, and Cloacibacterium. Across all organs, two amplicon sequence variants (ASVs) associated with the genera Diaphorobacter and Cloacibacterium appeared to be most abundant. In terms of the most prevalent ASVs, the alimentary canal, Malpighian tubules, and mandibular glands appeared to have similar complements of bacteria, with relatively few significant differences evident. However, aside from the Diaphorobacter and Cloacibacterium ASVs common to all the organs, the labial glands appeared to possess a distinctive complement of bacteria which was absent or poorly represented in the other organs. Among these were representatives of the Pseudomonas, Flavobacterium, Caulobacterium, Anaerococcus, and Methylobacterium. These results suggest that the labial glands present bacteria with different selective pressures than those occurring in the mandibular gland, Malpighian tubules and the alimentary canal. Given the documented effects that labial gland secretions and the gut microbiome can exert on host plant defenses, the effects exerted by the bacteria inhabiting the labial glands themselves deserve further study.