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

T Cells With Activated STAT4 Drive the High-Risk Rejection State to Renal Allograft Failure After Kidney Transplantation

In kidney transplantation, deteriorated progression of rejection is considered to be a leading course of postoperative mortality. However, the conventional histologic diagnosis is limited in reading the rejection status at the molecular level, thereby triggering mismatched pathogenesis with clinical phenotypes. Here, by applying uniform manifold approximation and projection and Leiden algorithms to 2,611 publicly available microarray datasets of renal transplantation, we uncovered six rejection states with corresponding signature genes and revealed a high-risk (HR) state that was essential in promoting allograft loss. By identifying cell populations from single-cell RNA sequencing data that were associated with the six rejection states, we identified a T-cell population to be the pathogenesis-triggering cells associated with the HR rejection state. Additionally, by constructing gene regulatory networks, we identified that activated STAT4, as a core transcription factor that was regulated by PTPN6 in T cells, was closely linked to poor allograft function and prognosis. Taken together, our study provides a novel strategy to help with the precise diagnosis of kidney allograft rejection progression, which is powerful in investigating the underlying molecular pathogenesis, and therefore, for further clinical intervention.

Tracking the turnover of SARS-CoV-2 VOCs Gamma to Delta in a Brazilian state (Minas Gerais) with a high-vaccination status

The emergence and global dissemination of Severe Acute Respiratory Syndrome virus 2 (SARS-CoV-2) variants of concern (VOCs) have been described as the main factor driving the Coronavirus Disease 2019 pandemic. In Brazil, the Gamma variant dominated the epidemiological scenario during the first period of 2021. Many Brazilian regions detected the Delta variant after its first description and documented its spread. To monitor the introduction and spread of VOC Delta, we performed Polymerase Chain Reaction (PCR) genotyping and genome sequencing in ten regional sentinel units from June to October 2021 in the State of Minas Gerais (MG). We documented the introduction and spread of Delta, comprising 70 per cent of the cases 8 weeks later. Comparing the viral loads of the Gamma and Delta dominance periods, we provide additional evidence that the latter is more transmissible. The spread and dominance of Delta did not culminate in the increase in cases and deaths, suggesting that the vaccination may have restrained the epidemic growth. Analysis of 224 novel Delta genomes revealed that Rio de Janeiro state was the primary source for disseminating this variant in the state of MG. We present the establishment of Delta, providing evidence of its enhanced transmissibility and showing that this variant shift did not aggravate the epidemiological scenario in a high immunity setting.

Multilocus Sequence Typing and Virulence Potential of Vibrio parahaemolyticus Strains Isolated from Aquatic Bird Feces

Vibrio parahaemolyticus is a Gram-negative, foodborne pathogenic bacterium that causes human gastroenteritis. This organism is ubiquitously present in the marine environment. Detection of V. parahaemolyticus in aquatic birds has been previously reported; however, the characterization of isolates of this bacterium recovered from these birds remains limited. The present study isolated and characterized V. parahaemolyticus from aquatic bird feces at the Bangpu Recreation Center (Samut Prakan province, Thailand) from 2016 to 2017, using multilocus sequence typing (MLST) and genome analysis. The results showed that V. parahaemolyticus was present in 34.9% (76/218) of the collected bird fecal samples. Among the ldh-positive V. parahaemolyticus isolates (n = 308), 1% (3/308) were positive for tdh, 1.3% (4/308) were positive for trh, and 0.3% (1/308) were positive for both tdh and trh. In turn, the MLST analysis revealed that 49 selected V. parahaemolyticus isolates resolved to 36 STs, 26 of which were novel (72.2%). Moreover, a total of 10 identified STs were identical to globally reported pathogenic strains (ST1309, ST1919, ST491, ST799, and ST2516) and environmental strains (ST1879, ST985, ST288, ST1925, and ST260). The genome analysis of isolates possessing tdh and/or trh (ST985, ST1923, ST1924, ST1929 and ST2516) demonstrated that the organization of the T3SS2α and T3SS2β genes in bird fecal isolates were almost identical to those of human clinical strains posing public health concerns of pathogen dissemination in the recreational area. The results of this study suggest that aquatic birds are natural reservoirs of new strains with high genetic diversity and are alternative sources of potentially pathogenic V. parahaemolyticus in the marine environment.

IMPORTANCE To our knowledge, infection of foodborne bacterium V. parahamolyticus occurs via the consumption of undercooked seafood contaminated with pathogenic strains. Aquatic bird is a neglectable source that can transmit V. parahaemolyticus along coastal areas. This study reported the detection of potentially pathogenic V. parahamolyticus harboring virulence genes from aquatic bird feces at the recreational center situated near the Gulf of Thailand. These strains shared identical genetic profile to the clinical isolates that previously reported in many countries. Furthermore, the strains from aquatic birds showed extremely high genetic diversity. Our research pointed out that the aquatic bird is possibly involved in the evolution of novel strains of V. parahaemolyticus and play a role in dissimilation of the potentially pathogenic strains across geographical distance.

Arabinoxylan and Pectin Metabolism in Crohn’s Disease Microbiota: An In Silico Study

Inflammatory bowel disease is a chronic disorder including ulcerative colitis and Crohn’s disease (CD). Gut dysbiosis is often associated with CD, and metagenomics allows a better understanding of the microbial communities involved. The objective of this study was to reconstruct in silico carbohydrate metabolic capabilities from metagenome-assembled genomes (MAGs) obtained from healthy and CD individuals. This computational method was developed as a mean to aid rationally designed prebiotic interventions to rebalance CD dysbiosis, with a focus on metabolism of emergent prebiotics derived from arabinoxylan and pectin. Up to 1196 and 1577 MAGs were recovered from CD and healthy people, respectively. MAGs of Akkermansia muciniphila, Barnesiella viscericola DSM 18177 and Paraprevotella xylaniphila YIT 11841 showed a wide range of unique and specific enzymes acting on arabinoxylan and pectin. These glycosidases were also found in MAGs recovered from CD patients. Interestingly, these arabinoxylan and pectin degraders are predicted to exhibit metabolic interactions with other gut microbes reduced in CD. Thus, administration of arabinoxylan and pectin may ameliorate dysbiosis in CD by promoting species with key metabolic functions, capable of cross-feeding other beneficial species. These computational methods may be of special interest for the rational design of prebiotic ingredients targeting at CD.

Blood-derived lncRNAs as biomarkers for cancer diagnosis: the Good, the Bad and the Beauty

Cancer ranks as one of the deadliest diseases worldwide. The high mortality rate associated with cancer is partially due to the lack of reliable early detection methods and/or inaccurate diagnostic tools such as certain protein biomarkers. Cell-free nucleic acids (cfNA) such as circulating long noncoding RNAs (lncRNAs) have been proposed as a new class of potential biomarkers for cancer diagnosis. The reported correlation between the presence of tumors and abnormal levels of lncRNAs in the blood of cancer patients has notably triggered a worldwide interest among clinicians and oncologists who have been actively investigating their potentials as reliable cancer biomarkers. In this report, we review the progress achieved ("the Good") and challenges encountered ("the Bad") in the development of circulating lncRNAs as potential biomarkers for early cancer diagnosis. We report and discuss the diagnostic performance of more than 50 different circulating lncRNAs and emphasize their numerous potential clinical applications ("the Beauty") including therapeutic targets and agents, on top of diagnostic and prognostic capabilities. This review also summarizes the best methods of investigation and provides useful guidelines for clinicians and scientists who desire conducting their own clinical studies on circulating lncRNAs in cancer patients via RT-qPCR or Next Generation Sequencing (NGS).

PIMGAVir and Vir-MinION: Two Viral Metagenomic Pipelines for Complete Baseline Analysis of 2nd and 3rd Generation Data

The taxonomic classification of viral sequences is frequently used for the rapid identification of pathogens, which is a key point for when a viral outbreak occurs. Both Oxford Nanopore Technologies (ONT) MinION and the Illumina (NGS) technology provide efficient methods to detect viral pathogens. Despite the availability of many strategies and software, matching them can be a very tedious and time-consuming task. As a result, we developed PIMGAVir and Vir-MinION, two metagenomics pipelines that automatically provide the user with a complete baseline analysis. The PIMGAVir and Vir-MinION pipelines work on 2nd and 3rd generation data, respectively, and provide the user with a taxonomic classification of the reads through three strategies: assembly-based, read-based, and clustering-based. The pipelines supply the scientist with comprehensive results in graphical and textual format for future analyses. Finally, the pipelines equip the user with a stand-alone platform with dedicated and various viral databases, which is a requirement for working in field conditions without internet connection.

Comparative Genomics of Fusarium circinatum Isolates Used to Screen Southern Pines for Pitch Canker Resistance

Pitch canker, caused by the fungal pathogen Fusarium circinatum, is a global disease affecting many Pinus spp. Often fatal, this disease causes significant mortality in both commercially grown and natural pine forests and is an issue of current and growing concern. F. circinatum isolates collected from three locations in the U.S. state of Florida were shown to be virulent on both slash and loblolly pine, with two of the isolates causing equivalent and significantly larger lesions than those caused by the third isolate during pathogenicity trials. In addition, significant genetic variation in lesion length in the pedigreed slash pine population was evident and rankings of parents for lesion length were similar across isolates. Experimental data demonstrate that both host and pathogen genetics contribute to disease severity. High-quality genomic assemblies of all three isolates were created and compared for structural differences and gene content. No major structural differences were observed among the isolates; however, missing or altered genes do contribute to genomic variation in the pathogen population. This work evaluates in planta virulence among three isolates of F. circinatum, provides genomic resources to facilitate study of this organism, and details comparative genomic methods that may be used to explore the pathogen's contribution to disease development.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Fine-tuned KDM1A alternative splicing regulates human cardiomyogenesis through an enzymatic-independent mechanism

The histone demethylase KDM1A is a multi-faceted regulator of vital developmental processes, including mesodermal and cardiac tube formation during gastrulation. However, it is unknown whether the fine-tuning of KDM1A splicing isoforms, already shown to regulate neuronal maturation, is crucial for the specification and maintenance of cell identity during cardiogenesis. Here, we discovered a temporal modulation of ubKDM1A and KDM1A+2a during human and mice fetal cardiac development and evaluated their impact on the regulation of cardiac differentiation. We revealed a severely impaired cardiac differentiation in KDM1A^−/− hESCs that can be rescued by re-expressing ubKDM1A or catalytically impaired ubKDM1A-K661A, but not by KDM1A+2a or KDM1A+2a-K661A. Conversely, KDM1A+2a^−/− hESCs give rise to functional cardiac cells, displaying increased beating amplitude and frequency and enhanced expression of critical cardiogenic markers. Our findings prove the existence of a divergent scaffolding role of KDM1A splice variants, independent of their enzymatic activity, during hESC differentiation into cardiac cells.

•

ubKDM1A and KDM1A+2a isoforms are fine-tuned during fetal cardiac development

•

Depletion of KDM1A isoforms impairs hESC differentiation into cardiac cells

•

KDM1A+2a ablation enhances the expression of key cardiac markers

•

KDM1A isoforms exhibit enzymatic-independent divergent roles during cardiogenesis

Decoding microbial genomes to understand their functional roles in human complex diseases

Complex diseases such as cardiovascular disease (CVD), obesity, inflammatory bowel disease (IBD), kidney disease, type 2 diabetes (T2D), and cancer have become a major burden to public health and affect more than 20% of the population worldwide. The etiology of complex diseases is not yet clear, but they are traditionally thought to be caused by genetics and environmental factors (e.g., dietary habits), and by their interactions. Besides this, increasing pieces of evidence now highlight that the intestinal microbiota may contribute substantially to the health and disease of the human host via their metabolic molecules. Therefore, decoding the microbial genomes has been an important strategy to shed light on their functional potential. In this review, we summarize the roles of the gut microbiome in complex diseases from its functional perspective. We further introduce artificial tools in decoding microbial genomes to profile their functionalities. Finally, state-of-the-art techniques have been highlighted which may contribute to a mechanistic understanding of the gut microbiome in human complex diseases and promote the development of the gut microbiome-based personalized medicine.

Eukaryotic Initiation Factor 5A2 Regulates Expression of Antiviral Genes

Translation factors are essential for regulation of protein synthesis. The eukaryotic translation initiation factor 5A (eIF5A) family is made up of two paralogues - eIF5A1 and eIF5A2 - which display high sequence homology but distinct tissue tropism. While eIF5A1 directly binds to the ribosome and regulates translation initiation, elongation, and termination, the molecular function of eIF5A2 remains poorly understood. Here, we engineer an eIF5A2 knockout allele in the SW480 colon cancer cell line. Using ribosome profiling and RNA-Sequencing, we reveal that eIF5A2 is functionally distinct from eIF5A1 and does not regulate transcript-specific or global protein synthesis. Instead, eIF5A2 knockout leads to decreased intrinsic antiviral gene expression, including members of the IFITM and APOBEC3 family. Furthermore, cells lacking eIF5A2 display increased permissiveness to virus infection. Our results uncover eIF5A2 as a factor involved regulating the antiviral transcriptome, and reveal an example of how gene duplications of translation factors can result in proteins with distinct functions.

γδ T cells license immature B cells to produce a broad range of polyreactive antibodies

Immature autoreactive B cells are present in all healthy individuals, but it is unclear which signals are required for their maturation into antibody-producing cells. Inducible depletion of γδ T cells show that direct interaction between γδ T cells and immature B cells in the spleen support an "innate" transition to mature B cells with a broad range of antigen specificities. IL-4 production of γδ T cells and cell-to-cell contact via CD30L support B cell maturation and induce genes of the unfolded protein response and mTORC1 signaling. Eight days after in vivo depletion of γδ T cells, increased numbers of B cells are already stuck in the transitional phase and express increased levels of IgD and CD21. Absence of γδ T cells leads also to reduced levels of serum anti-nuclear autoantibodies, making γδ T cells an attractive target to treat autoimmunity.

Diversity and distribution of sulfur metabolic genes in the human gut microbiome and their association with colorectal cancer

BACKGROUND

Recent evidence implicates microbial sulfidogenesis as a potential trigger of colorectal cancer (CRC), highlighting the need for comprehensive knowledge of sulfur metabolism within the human gut. Microbial sulfidogenesis produces genotoxic hydrogen sulfide (H₂S) in the human colon using inorganic (sulfate) and organic (taurine/cysteine/methionine) substrates; however, the majority of studies have focused on sulfate reduction using dissimilatory sulfite reductases (Dsr).

RESULTS

Here, we show that genes for microbial sulfur metabolism are more abundant and diverse than previously observed and are statistically associated with CRC. Using ~ 17,000 bacterial genomes from publicly available stool metagenomes, we studied the diversity of sulfur metabolic genes in 667 participants across different health statuses: healthy, adenoma, and carcinoma. Sulfidogenic genes were harbored by 142 bacterial genera and both organic and inorganic sulfidogenic genes were associated with carcinoma. Significantly, the anaerobic sulfite reductase (asr) genes were twice as abundant as dsr, demonstrating that Asr is likely a more important contributor to sulfate reduction in the human gut than Dsr. We identified twelve potential pathways for reductive taurine metabolism and discovered novel genera harboring these pathways. Finally, the prevalence of metabolic genes for organic sulfur indicates that these understudied substrates may be the most abundant source of microbially derived H₂S.

CONCLUSIONS

Our findings significantly expand knowledge of microbial sulfur metabolism in the human gut. We show that genes for microbial sulfur metabolism in the human gut are more prevalent than previously known, irrespective of health status (i.e., in both healthy and diseased states). Our results significantly increase the diversity of pathways and bacteria that are associated with microbial sulfur metabolism in the human gut. Overall, our results have implications for understanding the role of the human gut microbiome and its potential contributions to the pathogenesis of CRC. Video abstract.

Regulatory variants at 2q33.1 confer schizophrenia risk by modulating distal gene TYW5 expression

Genome-wide association studies have shown that genetic variants at 2q33.1 are strongly associated with schizophrenia. However, potential causal variants in this locus and their roles in schizophrenia remain unknown. Here, we identified two functional variants (rs796364 and rs281759) that disrupt CTCF, RAD21 and FOXP2 binding at 2q33.1. We systematically investigated the regulatory mechanisms of these two variants with serial experiments, including reporter gene assays and electrophoretic mobility shift assay. Intriguingly, these two single nucleotide polymorphisms physically interacted with TYW5 and showed the most significant associations with TYW5 expression in human brain. Consistently, CRISPR-Cas9-mediated genome editing confirmed the regulatory effect of the two single nucleotide polymorphisms on TYW5 expression. Additionally, expression analysis indicated that TYW5 was significantly upregulated in brains of schizophrenia cases compared with controls, suggesting that rs796364 and rs281759 might confer schizophrenia risk by modulating TYW5 expression. We over-expressed TYW5 in mouse neural stem cells and rat primary neurons to mimic its upregulation in schizophrenia and found significant alterations in the proliferation and differentiation of neural stem cells, as well as dendritic spine density following TYW5 overexpression, indicating its important roles in neurodevelopment and spine morphogenesis. Furthermore, we independently confirmed the association between rs796364 and schizophrenia in a Chinese cohort of 8202 subjects. Finally, transcriptome analysis revealed that TYW5 affected schizophrenia-associated pathways. These lines of evidence consistently revealed that rs796364 and rs281759 might contribute to schizophrenia risk by regulating the expression of TYW5, a gene whose expression dysregulation affects two important schizophrenia pathophysiological processes (i.e. neurodevelopment and dendritic spine formation).

Performance optimization in DNA short-read alignment

SUMMARY

Over the past decade, short-read sequence alignment has become a mature technology. Optimized algorithms, careful software engineering and high-speed hardware have contributed to greatly increased throughput and accuracy. With these improvements, many opportunities for performance optimization have emerged. In this review, we examine three general-purpose short-read alignment tools-BWA-MEM, Bowtie 2 and Arioc-with a focus on performance optimization. We analyze the performance-related behavior of the algorithms and heuristics each tool implements, with the goal of arriving at practical methods of improving processing speed and accuracy. We indicate where an aligner's default behavior may result in suboptimal performance, explore the effects of computational constraints such as end-to-end mapping and alignment scoring threshold, and discuss sources of imprecision in the computation of alignment scores and mapping quality. With this perspective, we describe an approach to tuning short-read aligner performance to meet specific data-analysis and throughput requirements while avoiding potential inaccuracies in subsequent analysis of alignment results. Finally, we illustrate how this approach avoids easily overlooked pitfalls and leads to verifiable improvements in alignment speed and accuracy.

CONTACT

richard.wilton@jhu.edu.

SUPPLEMENTARY INFORMATION

Appendices referenced in this article are available at Bioinformatics online.

Time-series transcriptome comparison reveals the gene regulation network under salt stress in soybean (Glycine max) roots

BACKGROUND

Soil salinity is a primary factor limiting soybean (Glycine max) productivity. Breeding soybean for tolerance to high salt conditions is therefore critical for increasing yield. To explore the molecular mechanism of soybean responses to salt stress, we performed a comparative transcriptome time-series analysis of root samples collected from two soybean cultivars with contrasting salt sensitivity.

RESULTS

The salt-tolerant cultivar 'Qi Huang No.34' (QH34) showed more differential expression of genes than the salt-sensitive cultivar 'Dong Nong No.50' (DN50). We identified 17,477 genes responsive to salt stress, of which 6644 exhibited distinct expression differences between the two soybean cultivars. We constructed the corresponding co-expression network and performed Gene Ontology term and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. The results suggested that phytohormone signaling, oxidoreduction, phenylpropanoid biosynthesis, the mitogen-activated protein kinase pathway and ribosome metabolism may play crucial roles in response to salt stress.

CONCLUSIONS

Our comparative analysis offers a comprehensive understanding of the genes involved in responding to salt stress and maintaining cell homeostasis in soybean. The regulatory gene networks constructed here also provide valuable molecular resources for future functional studies and breeding of soybean with improved tolerance to salinity.

Despite its sequence identity with canonical H4, Drosophila H4r product is enriched at specific chromatin regions

Histone variants are different from their canonical counterparts in structure and are encoded by solitary genes with unique regulation to fulfill tissue or differentiation specific functions. A single H4 variant gene (His4r or H4r) that is located outside of the histone cluster and gives rise to a polyA tailed messenger RNA via replication-independent expression is preserved in Drosophila strains despite that its protein product is identical with canonical H4. In order to reveal information on the possible role of this alternative H4 we epitope tagged endogenous H4r and studied its spatial and temporal expression, and revealed its genome-wide localization to chromatin at the nucleosomal level. RNA and immunohistochemistry analysis of H4r expressed under its cognate regulation indicate expression of the gene throughout zygotic and larval development and presence of the protein product is evident already in the pronuclei of fertilized eggs. In the developing nervous system a slight disequibrium in H4r distribution is observable, cholinergic neurons are the most abundant among H4r-expressing cells. ChIP-seq experiments revealed H4r association with regulatory regions of genes involved in cellular stress response. The data presented here indicate that H4r has a variant histone function.

Mapping bacterial diversity and metabolic functionality of the human respiratory tract microbiome

Background

The Human Respiratory Tract (HRT) is colonized by various microbial taxa, known as HRT microbiota, in a manner that is indicative of mutualistic interaction between such microorganisms and their host.

Aim

To investigate the microbial composition of the HRT and its possible correlation with the different compartments of the respiratory tract.

Methods

In the current study, we performed an in-depth meta-analysis of 849 HRT samples from public shotgun metagenomic datasets obtained through several distinct collection methods.

Results

The statistical robustness provided by this meta-analysis allowed the identification of 13 possible HRT-specific Community State Types (CSTs), which appear to be specific to each anatomical region of the respiratory tract. Furthermore, functional characterization of the metagenomic datasets revealed specific microbial metabolic features correlating with the different compartments of the respiratory tract.

Conclusion

The meta-analysis here performed suggested that the variable presence of certain bacterial species seems to be linked to a location-related abundance gradient in the HRT and seems to be characterized by a specific microbial metabolic capability.

Wt1 transcription factor impairs cardiomyocyte specification and drives a phenotypic switch from myocardium to epicardium

During development, the heart grows by addition of progenitor cells to the poles of the primordial heart tube. In the zebrafish, Wilms tumor 1 transcription factor a (wt1a) and b (wt1b) genes are expressed in the pericardium, at the venous pole of the heart. From this pericardial layer, the proepicardium emerges. Proepicardial cells are subsequently transferred to the myocardial surface and form the epicardium, covering the myocardium. We found that while wt1a and wt1b expression is maintained in proepicardial cells, it is downregulated in pericardial cells that contribute cardiomyocytes to the developing heart. Sustained wt1b expression in cardiomyocytes reduced chromatin accessibility of specific genomic loci. Strikingly, a subset of wt1a- and wt1b-expressing cardiomyocytes changed their cell-adhesion properties, delaminated from the myocardium and upregulated epicardial gene expression. Thus, wt1a and wt1b act as a break for cardiomyocyte differentiation, and ectopic wt1a and wt1b expression in cardiomyocytes can lead to their transdifferentiation into epicardial-like cells.

Phase variation in the glycosyltransferase genes of Pasteurella multocida associated with outbreaks of fowl cholera on free-range layer farms

Fowl cholera caused by Pasteurella multocida has re-emerged in Australian poultry production since the increasing adoption of free-range production systems. Currently, autogenous killed whole-cell vaccines prepared from the isolates previously obtained from each farm are the main preventative measures used. In this study, we use whole-genome sequencing and phylogenomic analysis to investigate outbreak dynamics, as well as monitoring and comparing the variations in the lipopolysaccharide (LPS) outer core biosynthesis loci of the outbreak and vaccine strains. In total, 73 isolates from two different free-range layer farms were included. Our genomic analysis revealed that all investigated isolates within the two farms (layer A and layer B) carried LPS type L3, albeit with a high degree of genetic diversity between them. Additionally, the isolates belonged to five different sequence types (STs), with isolates belonging to ST9 and ST20 being the most prevalent. The isolates carried ST-specific mutations within their LPS type L3 outer core biosynthesis loci, including frameshift mutations in the outer core heptosyltransferase gene (htpE) (ST7 and ST274) or galactosyltransferase gene (gatG) (ST20). The ST9 isolates could be separated into three groups based on their LPS outer core biosynthesis loci sequences, with evidence for potential phase variation mechanisms identified. The potential phase variation mechanisms included a tandem repeat insertion in natC and a single base deletion in a homopolymer region of gatG. Importantly, our results demonstrated that two of the three ST9 groups shared identical rep-PCR (repetitive extragenic palindromic PCR) patterns, while carrying differences in their LPS outer core biosynthesis loci region. In addition, we found that ST9 isolates either with or without the natC tandem repeat insertion were both associated with a single outbreak, which would indicate the importance of screening more than one isolate within an outbreak. Our results strongly suggest the need for a metagenomics culture-independent approach, as well as a genetic typing scheme for LPS, to ensure an appropriate vaccine strain with a matching predicted LPS structure is used.

Taxonomic distribution and evolutionary analysis of the equol biosynthesis gene cluster

Background

Equol, an isoflavonoid metabolite with possible health benefits in humans, is known to be produced by some human gut bacteria. While the genes encoding the equol production pathway have been characterized in a few bacterial strains, a systematic analysis of the equol production pathway is currently lacking.

Results

This study presents an analysis of the taxonomic distribution and evolutionary history of the gene cluster encoding the equol production pathway. A survey for equol gene clusters within the Genome Taxonomy Database bacterial genomes and human gut metagenomes resulted in the identification of a highly conserved gene cluster found in nine bacterial species from the Eggerthellaceae family. The identified gene clusters from human gut metagenomes revealed potential variations in the equol gene cluster organization and gene content within the equol-producing Eggerthellaceae clades. Subsequent analysis showed that in addition to the four genes directly involved in equol production, multiple other genes were consistently found in the equol gene clusters. These genes were predicted to encode a putative electron transport complex and hydrogenase maturase system, suggesting potential roles for them in the equol production pathway. Analysis of the gene clusters and a phylogenetic reconstruction of a putative NAD kinase gene provided evidence of the recent transfer of the equol gene cluster from a basal Eggerthellaceae species to Slackia_A equolifaciens, Enteroscipio sp000270285, and Lactococcus garvieae 20–92.

Conclusions

This analysis demonstrates that the highly conserved equol gene cluster is taxonomically restricted to the Eggerthellaceae family of bacteria and provides evidence of the role of horizontal gene transfer in the evolutionary history of these genes. These results provide a foundation for future studies of equol production in the human gut and future efforts related to bioengineering and the use of equol-producing bacteria as probiotics.

Analysis and Performance Assessment of the Whole Genome Bisulfite Sequencing Data Workflow: Currently Available Tools and a Practical Guide to Advance DNA Methylation Studies

DNA methylation is associated with transcriptional repression, genomic imprinting, stem cell differentiation, embryonic development, and inflammation. Aberrant DNA methylation can indicate disease states, including cancer and neurological disorders. Therefore, the prevalence and location of 5-methylcytosine in the human genome is a topic of interest. Whole-genome bisulfite sequencing (WGBS) is a high-throughput method for analyzing DNA methylation. This technique involves library preparation, alignment, and quality control. Advancements in epigenetic technology have led to an increase in DNA methylation studies. This review compares the detailed experimental methodology of WGBS using accessible and up-to-date analysis tools. Practical codes for WGBS data processing are included as a general guide to assist progress in DNA methylation studies through a comprehensive case study.

Integrative Comparison of Burrows-Wheeler Transform-Based Mapping Algorithm with de Bruijn Graph for Identification of Lung/Liver Cancer-Specific Gene

Cancers of the lung and liver are the top 10 leading causes of cancer death worldwide. Thus, it is essential to identify the genes specifically expressed in these two cancer types to develop new therapeutics. Although many messenger RNA (mRNA) sequencing data related to these cancer cells are available due to the advancement of next-generation sequencing (NGS) technologies, optimized data processing methods need to be developed to identify the novel cancer-specific genes. Here, we conducted an analytical comparison between Bowtie2, a Burrows-Wheeler transform-based alignment tool, and Kallisto, which adopts pseudo alignment based on a transcriptome de Bruijn graph using mRNA sequencing data on normal cells and lung/liver cancer tissues. Before using cancer data, simulated mRNA sequencing reads were generated, and the high Transcripts Per Million (TPM) values were compared. mRNA sequencing reads data on lung/liver cancer cells were also extracted and quantified. While Kallisto could directly give the output in TPM values, Bowtie2 provided the counts. Thus, TPM values were calculated by processing the Sequence Alignment Map (SAM) file in R using package Rsubread and subsequently in python. The analysis of the simulated sequencing data revealed that Kallisto could detect more transcripts and had a higher overlap over Bowtie2. The evaluation of these two data processing methods using the known lung cancer biomarkers concludes that in standard settings without any dedicated quality control, Kallisto is more effective at producing faster and more accurate results than Bowtie2. Such conclusions were also drawn and confirmed with the known biomarkers specific to liver cancer.

PRDM paralogs antagonistically balance Wnt/β-catenin activity during craniofacial chondrocyte differentiation

Cranial neural crest cell (NCC)-derived chondrocyte precursors undergo a dynamic differentiation and maturation process to establish a scaffold for subsequent bone formation, alterations in which contribute to congenital birth defects. Here, we demonstrate that transcription factor and histone methyltransferase proteins Prdm3 and Prdm16 control the differentiation switch of cranial NCCs to craniofacial cartilage. Loss of either paralog results in hypoplastic and disorganized chondrocytes due to impaired cellular orientation and polarity. We show that these proteins regulate cartilage differentiation by controlling the timing of Wnt/β-catenin activity in strikingly different ways: Prdm3 represses whereas Prdm16 activates global gene expression, although both act by regulating Wnt enhanceosome activity and chromatin accessibility. Finally, we show that manipulating Wnt/β-catenin signaling pharmacologically or generating prdm3-/-;prdm16-/- double mutants rescues craniofacial cartilage defects. Our findings reveal upstream regulatory roles for Prdm3 and Prdm16 in cranial NCCs to control Wnt/β-catenin transcriptional activity during chondrocyte differentiation to ensure proper development of the craniofacial skeleton.

Whole-Genome Sequence of the Wolbachia Strain wTcon, an Endosymbiont of the Confused Flour Beetle, Tribolium confusum

Up to 60% of insects are infected with symbiont intracellular alphaproteobacteria of the genus Wolbachia, which are often able to manipulate their host's reproduction. Here, we report the annotated draft genome sequence of strain wTcon from the confused flour beetle, Tribolium confusum, based on long- and short-read sequence data. The assembled genome is located on 12 contigs with a total size of 1,418,452 bp.

RdRp-based sensitive taxonomic classification of RNA viruses for metagenomic data

With advances in library construction protocols and next-generation sequencing technologies, viral metagenomic sequencing has become the major source for novel virus discovery. Conducting taxonomic classification for metagenomic data is an important means to characterize the viral composition in the underlying samples. However, RNA viruses are abundant and highly diverse, jeopardizing the sensitivity of comparison-based classification methods. To improve the sensitivity of read-level taxonomic classification, we developed an RNA-dependent RNA polymerase (RdRp) gene-based read classification tool RdRpBin. It combines alignment-based strategy with machine learning models in order to fully exploit the sequence properties of RdRp. We tested our method and compared its performance with the state-of-the-art tools on the simulated and real sequencing data. RdRpBin competes favorably with all. In particular, when the query RNA viruses share low sequence similarity with the known viruses (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim 0.4$\end{document}), our tool can still maintain a higher F-score than the state-of-the-art tools. The experimental results on real data also showed that RdRpBin can classify more RNA viral reads with a relatively low false-positive rate. Thus, RdRpBin can be utilized to classify novel and diverged RNA viruses.

The gut virome in two indigenous populations from Malaysia

The human gut contains a complex microbiota dominated by bacteriophages but also containing other viruses and bacteria and fungi. There are a growing number of techniques for the extraction, sequencing, and analysis of the virome but currently no standardized protocols. This study established an effective workflow for virome analysis to investigate the virome of stool samples from two understudied ethnic groups from Malaysia: the Jakun and Jehai Orang Asli. By using the virome extraction and analysis workflow with the Oxford Nanopore Technology, long-read sequencing successfully captured close to full-length viral genomes. The virome composition of the two indigenous Malaysian communities were remarkably different from those found in other parts of the world. Additionally, plant viruses found in the viromes of these individuals were attributed to traditional food-seeking methods. This study establishes a human gut virome workflow and extends insights into the healthy human gut virome, laying the groundwork for comparative studies.

Genome diversity of domesticated Acinetobacter baumannii ATCC 19606T strains

Acinetobacter baumannii has emerged as an important opportunistic pathogen worldwide, being responsible for large outbreaks for nosocomial infections, primarily in intensive care units. A. baumannii ATCC 19606^T is the species type strain, and a reference organism in many laboratories due to its low virulence, amenability to genetic manipulation and extensive antibiotic susceptibility. We wondered if frequent propagation of A. baumannii ATCC 19606^T in different laboratories may have driven micro- and macro-evolutionary events that could determine inter-laboratory differences of genome-based data. By combining Illumina MiSeq, MinION and Sanger technologies, we generated a high-quality whole-genome sequence of A. baumannii ATCC 19606^T, then performed a comparative genome analysis between A. baumannii ATCC 19606^T strains from several research laboratories and a reference collection. Differences between publicly available ATCC 19606^T genome sequences were observed, including SNPs, macro- and micro-deletions, and the uneven presence of a 52 kb prophage belonging to genus Vieuvirus. Two plasmids, pMAC and p1ATCC19606, were invariably detected in all tested strains. The presence of a putative replicase, a replication origin containing four 22-mer direct repeats, and a toxin-antitoxin system implicated in plasmid stability were predicted by in silico analysis of p1ATCC19606, and experimentally confirmed. This work refines the sequence, structure and functional annotation of the A. baumannii ATCC 19606^T genome, and highlights some remarkable differences between domesticated strains, likely resulting from genetic drift.

The virulome of Streptomyces scabiei in response to cello-oligosaccharide elicitors

The development of spots or lesions symptomatic of common scab on root and tuber crops is caused by few pathogenic Streptomyces with Streptomyces scabiei 87–22 as the model species. Thaxtomin phytotoxins are the primary virulence determinants, mainly acting by impairing cellulose synthesis, and their production in S. scabiei is in turn boosted by cello-oligosaccharides released from host plants. In this work we aimed to determine which molecules and which biosynthetic gene clusters (BGCs) of the specialized metabolism of S. scabiei 87–22 show a production and/or a transcriptional response to cello-oligosaccharides. Comparative metabolomic analyses revealed that molecules of the virulome of S. scabiei induced by cellobiose and cellotriose include (i) thaxtomin and concanamycin phytotoxins, (ii) desferrioxamines, scabichelin and turgichelin siderophores in order to acquire iron essential for housekeeping functions, (iii) ectoine for protection against osmotic shock once inside the host, and (iv) bottromycin and concanamycin antimicrobials possibly to prevent other microorganisms from colonizing the same niche. Importantly, both cello-oligosaccharides reduced the production of the spore germination inhibitors germicidins thereby giving the ‘green light’ to escape dormancy and trigger the onset of the pathogenic lifestyle. For most metabolites - either with induced or reduced production - cellotriose was revealed to be a slightly stronger elicitor compared to cellobiose, supporting an earlier hypothesis which suggested the trisaccharide was the real trigger for virulence released from the plant cell wall through the action of thaxtomins. Interestingly, except for thaxtomins, none of these BGCs’ expression seems to be under direct control of the cellulose utilization repressor CebR suggesting the existence of a yet unknown mechanism for switching on the virulome. Finally, a transcriptomic analysis revealed nine additional cryptic BGCs that have their expression awakened by cello-oligosaccharides, suggesting that other and yet to be discovered metabolites could be part of the virulome of S. scabiei.

regCNN: identifying Drosophila genome-wide cis-regulatory modules via integrating the local patterns in epigenetic marks and transcription factor binding motifs

Transcription regulation in metazoa is controlled by the binding events of transcription factors (TFs) or regulatory proteins on specific modular DNA regulatory sequences called cis-regulatory modules (CRMs). Understanding the distributions of CRMs on a genomic scale is essential for constructing the metazoan transcriptional regulatory networks that help diagnose genetic disorders. While traditional reporter-assay CRM identification approaches can provide an in-depth understanding of functions of some CRM, these methods are usually cost-inefficient and low-throughput. It is generally believed that by integrating diverse genomic data, reliable CRM predictions can be made. Hence, researchers often first resort to computational algorithms for genome-wide CRM screening before specific experiments. However, current existing in silico methods for searching potential CRMs were restricted by low sensitivity, poor prediction accuracy, or high computation time from TFBS composition combinatorial complexity. To overcome these obstacles, we designed a novel CRM identification pipeline called regCNN by considering the base-by-base local patterns in TF binding motifs and epigenetic profiles. On the test set, regCNN shows an accuracy/auROC of 84.5%/92.5% in CRM identification. And by further considering local patterns in epigenetic profiles and TF binding motifs, it can accomplish 4.7% (92.5%–87.8%) improvement in the auROC value over the average value-based pure multi-layer perceptron model. We also demonstrated that regCNN outperforms all currently available tools by at least 11.3% in auROC values. Finally, regCNN is verified to be robust against its resizing window hyperparameter in dealing with the variable lengths of CRMs. The model of regCNN can be downloaded athttp://cobisHSS0.im.nuk.edu.tw/regCNN/.

Exposure to soil environments during earlier life stages is distinguishable in the gut microbiome of adult mice

Environmental exposure during earlier life stages can govern the assembly and development of gut microbiota, yet it is insufficiently understood. In this study, ex-germ-free mice were cohoused with distinct soil-microbiota (from desert, steppe, and forest) beddings within 60 days after birth and subsequently transferred to new soil beddings from 60 to 90th day. Using metagenomic shotgun sequencing, firstly, we found soil microbes from natural environments (birthplace) greatly influenced the gut community assembly in the housing experiment. About 27% microbial species and 12% functional components that associated with birthplaces at Day 60 were still discriminatory of birthplaces after transferring mice to new environments. Moreover, prior soil-exposure types are associated with the magnitude of temporal microbiome change due to environmental shifts. The appropriate soil-exposure (e.g., steppe) might help mice gut microbiome adapt to changing environments or host development. Our study demonstrated the continuous soil-exposure history earlier is associated with the gut microbiome individuality and development later.

Inverting the model of genomics data sharing with the NHGRI Genomic Data Science Analysis, Visualization, and Informatics Lab-space

The NHGRI Genomic Data Science Analysis, Visualization, and Informatics Lab-space (AnVIL; https://anvilproject.org) was developed to address a widespread community need for a unified computing environment for genomics data storage, management, and analysis. In this perspective, we present AnVIL, describe its ecosystem and interoperability with other platforms, and highlight how this platform and associated initiatives contribute to improved genomic data sharing efforts. The AnVIL is a federated cloud platform designed to manage and store genomics and related data, enable population-scale analysis, and facilitate collaboration through the sharing of data, code, and analysis results. By inverting the traditional model of data sharing, the AnVIL eliminates the need for data movement while also adding security measures for active threat detection and monitoring and provides scalable, shared computing resources for any researcher. We describe the core data management and analysis components of the AnVIL, which currently consists of Terra, Gen3, Galaxy, RStudio/Bioconductor, Dockstore, and Jupyter, and describe several flagship genomics datasets available within the AnVIL. We continue to extend and innovate the AnVIL ecosystem by implementing new capabilities, including mechanisms for interoperability and responsible data sharing, while streamlining access management. The AnVIL opens many new opportunities for analysis, collaboration, and data sharing that are needed to drive research and to make discoveries through the joint analysis of hundreds of thousands to millions of genomes along with associated clinical and molecular data types.

Spatiotemporal expression profile of novel and known small RNAs throughout rice plant development focussing on seed tissues

Background

Small RNAs (sRNAs) regulate numerous plant processes directly related to yield, such as disease resistance and plant growth. To exploit this yield-regulating potential of sRNAs, the sRNA profile of one of the world’s most important staple crops – rice – was investigated throughout plant development using next-generation sequencing.

Results

Root and leaves were investigated at both the vegetative and generative phase, and early-life sRNA expression was characterized in the embryo and endosperm. This led to the identification of 49,505 novel sRNAs and 5581 tRNA-derived sRNAs (tsRNAs). In all tissues, 24 nt small interfering RNAs (siRNAs) were highly expressed and associated with euchromatic, but not heterochromatic transposable elements. Twenty-one nt siRNAs deriving from genic regions in the endosperm were exceptionally highly expressed, mimicking previously reported expression levels of 24 nt siRNAs in younger endosperm samples. In rice embryos, sRNA content was highly diverse while tsRNAs were underrepresented, possibly due to snoRNA activity. Publicly available mRNA expression and DNA methylation profiles were used to identify putative siRNA targets in embryo and endosperm. These include multiple genes related to the plant hormones gibberellic acid and ethylene, and to seed phytoalexin and iron content.

Conclusions

This work introduces multiple sRNAs as potential regulators of rice yield and quality, identifying them as possible targets for the continuous search to optimize rice production.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08264-z.

Patterns of piRNA Regulation in Drosophila Revealed through Transposable Element Clade Inference

Transposable elements (TEs) are self-replicating "genetic parasites" ubiquitous to eukaryotic genomes. In addition to conflict between TEs and their host genomes, TEs of the same family are in competition with each other. They compete for the same genomic niches while experiencing the same regime of copy-number selection. This suggests that competition among TEs may favor the emergence of new variants that can outcompete their ancestral forms. To investigate the sequence evolution of TEs, we developed a method to infer clades: collections of TEs that share SNP variants and represent distinct TE family lineages. We applied this method to a panel of 85 Drosophila melanogaster genomes and found that the genetic variation of several TE families shows significant population structure that arises from the population-specific expansions of single clades. We used population genetic theory to classify these clades into younger versus older clades and found that younger clades are associated with a greater abundance of sense and antisense piRNAs per copy than older ones. Further, we find that the abundance of younger, but not older clades, is positively correlated with antisense piRNA production, suggesting a general pattern where hosts preferentially produce antisense piRNAs from recently active TE variants. Together these findings suggest a pattern whereby new TE variants arise by mutation and then increase in copy number, followed by the host producing antisense piRNAs that may be used to silence these emerging variants.

A drug discovery approach based on comparative transcriptomics between two toxin-secreting marine annelids: Glycera alba and Hediste diversicolor

While Glycera alba secretes neurotoxins, Hediste diversicolor may secrete fewer toxins with a broader action. Transcriptomics and human interactome-directed analysis unraveled promising candidates for biomedical applications from either annelid.

Myeloid cyclooxygenase-2/prostaglandin E2/E-type prostanoid receptor 4 promotes transcription factor MafB-dependent inflammatory resolution in acute kidney injury

Following acute injury to the kidney, macrophages play an important role in recovery of functional and structural integrity, but organ fibrosis and progressive functional decline occur with incomplete recovery. Pro-resolving macrophages are characterized by increased cyclooxygenase 2 (COX-2) expression and this expression was selectively increased in kidney macrophages following injury and myeloid-specific COX-2 deletion inhibited recovery. Deletion of the myeloid prostaglandin E2 (PGE2) receptor, E-type prostanoid receptor 4 (EP4), mimicked effects seen with myeloid COX-2-/- deletion. PGE2-mediated EP4 activation induced expression of the transcription factor MafB in kidney macrophages, which upregulated anti-inflammatory genes and suppressed pro-inflammatory genes. Myeloid Mafb deletion recapitulated the effects seen with either myeloid COX-2 or EP4 deletion following acute kidney injury, with delayed recovery, persistent presence of pro-inflammatory kidney macrophages, and increased kidney fibrosis. Thus, our studies identified a previously unknown mechanism by which prostaglandins modulate macrophage phenotype following acute organ injury and provide new insight into mechanisms underlying detrimental kidney effects of non-steroidal anti-inflammatory drugs that inhibit cyclooxygenase activity.

Identification and Molecular Characterization of Novel Mycoviruses in Saccharomyces and Non-Saccharomyces Yeasts of Oenological Interest

Wine yeasts can be natural hosts for dsRNA, ssRNA viruses and retrotransposon elements. In this study, high-throughput RNA sequencing combined with bioinformatic analyses unveiled the virome associated to 16 Saccharomyces cerevisiae and 8 non-Saccharomyces strains of oenological interest. Results showed the presence of six viruses and two satellite dsRNAs from four different families, two of which-Partitiviridae and Mitoviridae-were not reported before in yeasts, as well as two ORFan contigs of viral origin. According to phylogenetic analysis, four new putative mycoviruses distributed in Totivirus, Cryspovirus, and Mitovirus genera were identified. The majority of commercial S. cerevisiae strains were confirmed to be the host for helper L-A type totiviruses and satellite M dsRNAs associated with the killer phenotype, both in single and mixed infections with L-BC totiviruses, and two viral sequences belonging to a new cryspovirus putative species discovered here for the first time. Moreover, single infection by a narnavirus 20S-related sequence was also found in one S. cerevisiae strain. Considering the non-Saccharomyces yeasts, Starmerella bacillaris hosted four RNAs of viral origin-two clustering in Totivirus and Mitovirus genera, and two ORFans with putative satellite behavior. This study confirmed the infection of wine yeasts by viruses associated with useful technological characteristics and demonstrated the presence of complex mixed infections with unpredictable biological effects.

Distinct gene expression pattern of RUNX1 mutations coordinated by target repression and promoter hypermethylation in acute myeloid leukemia

Runt-related transcription factor 1 (RUNX1) is an essential regulator of normal hematopoiesis. Its dysfunction, caused by either fusions or mutations, is frequently reported in acute myeloid leukemia (AML). However, RUNX1 mutations have been largely under-explored compared with RUNX1 fusions mainly due to their elusive genetic characteristics. Here, based on 1741 patients with AML, we report a unique expression pattern associated with RUNX1 mutations in AML. This expression pattern was coordinated by target repression and promoter hypermethylation. We first reanalyzed a joint AML cohort that consisted of three public cohorts and found that RUNX1 mutations were mainly distributed in the Runt domain and almost mutually exclusive with NPM1 mutations. Then, based on RNA-seq data from The Cancer Genome Atlas AML cohort, we developed a 300-gene signature that significantly distinguished the patients with RUNX1 mutations from those with other AML subtypes. Furthermore, we explored the mechanisms underlying this signature from the transcriptional and epigenetic levels. Using chromatin immunoprecipitation sequencing data, we found that RUNX1 target genes tended to be repressed in patients with RUNX1 mutations. Through the integration of DNA methylation array data, we illustrated that hypermethylation on the promoter regions of RUNX1-regulated genes also contributed to dysregulation in RUNX1-mutated AML. This study revealed the distinct gene expression pattern of RUNX1 mutations and the underlying mechanisms in AML development.

The Mouse Gastrointestinal Bacteria Catalogue enables translation between the mouse and human gut microbiotas via functional mapping

Human health and disease have increasingly been shown to be impacted by the gut microbiota, and mouse models are essential for investigating these effects. However, the compositions of human and mouse gut microbiotas are distinct, limiting translation of microbiota research between these hosts. To address this, we constructed the Mouse Gastrointestinal Bacteria Catalogue (MGBC), a repository of 26,640 high-quality mouse microbiota-derived bacterial genomes. This catalog enables species-level analyses for mapping functions of interest and identifying functionally equivalent taxa between the microbiotas of humans and mice. We have complemented this with a publicly deposited collection of 223 bacterial isolates, including 62 previously uncultured species, to facilitate experimental investigation of individual commensal bacteria functions in vitro and in vivo. Together, these resources provide the ability to identify and test functionally equivalent members of the host-specific gut microbiotas of humans and mice and support the informed use of mouse models in human microbiota research.

The role of epithelial progesterone receptor isoforms in embryo implantation

The loss of uterine epithelial progesterone receptor (PGR) is crucial for successful embryo implantation in both humans and mice. The two major isoforms PGRA and PGRB have divergent functions under both physiological and pathological conditions. The present study compares phenotypes and gene signatures of PGRA and PGRB in uterine epithelium using uterine epithelial-specific constitutively expressed PGRA or PGRB mouse models. The cistrome and transcriptome analysis reveals substantial overlap between epithelial PGRA and PGRB, and both disrupt embryo implantation through FOXO1 pathways. Constitutive epithelial PGRA and PGRB expression impairs ESR1 occupancy at the promoter of Lif leading to reduced Lif transcription and further exaggerates SGK1 expression leading to enhanced PI3K-SGK1 activities, and both contribute to the decline of nuclear FOXO1 expression. Our study demonstrates that PGRA and PGRB in the uterine epithelium act on a similar set of target genes and commonly regulate the LIF-SGK1-FOXO1 signaling pathway for embryo implantation.

UBR7 acts as a histone chaperone for post-nucleosomal histone H3

Histone chaperones modulate the stability of histones beginning from histone synthesis, through incorporation into DNA, and during recycling during transcription and replication. Following histone removal from DNA, chaperones regulate histone storage and degradation. Here, we demonstrate that UBR7 is a histone H3.1 chaperone that modulates the supply of pre-existing post-nucleosomal histone complexes. We demonstrate that UBR7 binds to post-nucleosomal H3K4me3 and H3K9me3 histones via its UBR box and PHD. UBR7 binds to the non-nucleosomal histone chaperone NASP. In the absence of UBR7, the pool of NASP-bound post-nucleosomal histones accumulate and chromatin is depleted of H3K4me3-modified histones. We propose that the interaction of UBR7 with NASP and histones opposes the histone storage functions of NASP and that UBR7 promotes reincorporation of post-nucleosomal H3 complexes.

Secondary metabolite biosynthetic diversity in Arctic Ocean metagenomes

Polyketide synthases (PKSs) and non-ribosomal peptide synthetases (NRPSs) are mega enzymes responsible for the biosynthesis of a large fraction of natural products (NPs). Molecular markers for biosynthetic genes, such as the ketosynthase (KS) domain of PKSs, have been used to assess the diversity and distribution of biosynthetic genes in complex microbial communities. More recently, metagenomic studies have complemented and enhanced this approach by allowing the recovery of complete biosynthetic gene clusters (BGCs) from environmental DNA. In this study, the distribution and diversity of biosynthetic genes and clusters from Arctic Ocean samples (NICE-2015 expedition), was assessed using PCR-based strategies coupled with high-throughput sequencing and metagenomic analysis. In total, 149 KS domain OTU sequences were recovered, 36 % of which could not be assigned to any known BGC. In addition, 74 bacterial metagenome-assembled genomes were recovered, from which 179 BGCs were extracted. A network analysis identified potential new NP families, including non-ribosomal peptides and polyketides. Complete or near-complete BGCs were recovered, which will enable future heterologous expression efforts to uncover the respective NPs. Our study represents the first report of biosynthetic diversity assessed for Arctic Ocean metagenomes and highlights the potential of Arctic Ocean planktonic microbiomes for the discovery of novel secondary metabolites. The strategy employed in this study will enable future bioprospection, by identifying promising samples for bacterial isolation efforts, while providing also full-length BGCs for heterologous expression.

Novel Tyrosine Kinase-Mediated Phosphorylation With Dual Specificity Plays a Key Role in the Modulation of Streptococcus pyogenes Physiology and Virulence

Streptococcus pyogenes (Group A Streptococcus, GAS) genomes do not contain a gene encoding a typical bacterial-type tyrosine kinase (BY-kinase) but contain an orphan gene-encoding protein Tyr-phosphatase (SP-PTP). Hence, the importance of Tyr-phosphorylation is underappreciated and not recognized for its role in GAS pathophysiology and pathogenesis. The fact that SP-PTP dephosphorylates Abl-tyrosine kinase-phosphorylated myelin basic protein (MBP), and SP-STK (S. pyogenes Ser/Thr kinase) also autophosphorylates its Tyr¹⁰¹-residue prompted us to identify a putative tyrosine kinase and Tyr-phosphorylation in GAS. Upon a genome-wide search of kinases possessing a classical Walker motif, we identified a non-canonical tyrosine kinase M5005_Spy_1476, a ∼17 kDa protein (153 aa) (SP-TyK). The purified recombinant SP-TyK autophosphorylated in the presence of ATP. In vitro and in vivo phosphoproteomic analyses revealed two key phosphorylated tyrosine residues located within the catalytic domain of SP-TyK. An isogenic mutant lacking SP-TyK derived from the M1T1 strain showed a retarded growth pattern. It displayed defective cell division and long chains with multiple parallel septa, often resulting in aggregates. Transcriptomic analysis of the mutant revealed 287 differentially expressed genes responsible for GAS pathophysiology and pathogenesis. SP-TyK also phosphorylated GAS CovR, WalR, SP-STP, and SDH/GAPDH proteins with dual specificity targeting their Tyr/Ser/Thr residues as revealed by biochemical and mass-spectrometric-based phosphoproteomic analyses. SP-TyK-phosphorylated CovR bound to PcovR efficiently. The mutant displayed sustained release of IL-6 compared to TNF-α during co-culturing with A549 lung cell lines, attenuation in mice sepsis model, and significantly reduced ability to adhere to and invade A549 lung cells and form biofilms on abiotic surfaces. SP-TyK, thus, plays a critical role in fine-tuning the regulation of key cellular functions essential for GAS pathophysiology and pathogenesis through post-translational modifications and hence, may serve as a promising target for future therapeutic developments.

Recurrent integration of human papillomavirus genomes at transcriptional regulatory hubs

Oncogenic human papillomavirus (HPV) genomes are often integrated into host chromosomes in HPV-associated cancers. HPV genomes are integrated either as a single copy or as tandem repeats of viral DNA interspersed with, or without, host DNA. Integration occurs frequently in common fragile sites susceptible to tandem repeat formation and the flanking or interspersed host DNA often contains transcriptional enhancer elements. When co-amplified with the viral genome, these enhancers can form super-enhancer-like elements that drive high viral oncogene expression. Here we compiled highly curated datasets of HPV integration sites in cervical (CESC) and head and neck squamous cell carcinoma (HNSCC) cancers, and assessed the number of breakpoints, viral transcriptional activity, and host genome copy number at each insertion site. Tumors frequently contained multiple distinct HPV integration sites but often only one “driver” site that expressed viral RNA. As common fragile sites and active enhancer elements are cell-type-specific, we mapped these regions in cervical cell lines using FANCD2 and Brd4/H3K27ac ChIP-seq, respectively. Large enhancer clusters, or super-enhancers, were also defined using the Brd4/H3K27ac ChIP-seq dataset. HPV integration breakpoints were enriched at both FANCD2-associated fragile sites and enhancer-rich regions, and frequently showed adjacent focal DNA amplification in CESC samples. We identified recurrent integration “hotspots” that were enriched for super-enhancers, some of which function as regulatory hubs for cell-identity genes. We propose that during persistent infection, extrachromosomal HPV minichromosomes associate with these transcriptional epicenters and accidental integration could promote viral oncogene expression and carcinogenesis.

Kr-h1 maintains distinct caste-specific neurotranscriptomes in response to socially regulated hormones

Behavioral plasticity is key to animal survival. Harpegnathos saltator ants can switch between worker and queen-like status (gamergate) depending on the outcome of social conflicts, providing an opportunity to study how distinct behavioral states are achieved in adult brains. Using social and molecular manipulations in live ants and ant neuronal cultures, we show that ecdysone and juvenile hormone drive molecular and functional differences in the brains of workers and gamergates and direct the transcriptional repressor Kr-h1 to different target genes. Depletion of Kr-h1 in the brain caused de-repression of "socially inappropriate" genes: gamergate genes were upregulated in workers, whereas worker genes were upregulated in gamergates. At the phenotypic level, loss of Kr-h1 resulted in the emergence of worker-specific behaviors in gamergates and gamergate-specific traits in workers. We conclude that Kr-h1 is a transcription factor that maintains distinct brain states established in response to socially regulated hormones.

Genome-Wide DNA Methylation Profile in Jejunum Reveals the Potential Genes Associated With Paratuberculosis in Dairy Cattle

Paratuberculosis in cattle causes substantial economic losses to the dairy industry. Exploring functional genes and corresponding regulatory pathways related to resistance or susceptibility to paratuberculosis is essential to the breeding of disease resistance in cattle. Co-analysis of genome-wide DNA methylation and transcriptome profiles is a critically important approach to understand potential regulatory mechanism underlying the development of diseases. In this study, we characterized the profiles of DNA methylation of jejunum from nine Holstein cows in clinical, subclinical, and healthy groups using whole-genome bisulfite sequencing (WGBS). The average methylation level in functional regions was 29.95% in the promoter, 29.65% in the 5’ untranslated region (UTR), 68.24% in exons, 71.55% in introns, and 72.81% in the 3’ UTR. A total of 3,911, 4,336, and 4,094 differentially methylated genes (DMGs) were detected in clinical vs. subclinical, clinical vs. healthy, and subclinical vs. healthy comparative group, respectively. Gene ontology (GO) and analysis based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that these DMGs were significantly enriched in specific biological processes related to immune response, such as Th1 and Th2 cell differentiation, wnt, TNF, MAPK, ECM-receptor interaction, cellular senescence, calcium, and chemokine signaling pathways (q value <0.05). The integration of information about DMGs, differentially expressed genes (DEGs), and biological functions suggested nine genes CALCRL, TNC, GATA4, CD44, TGM3, CXCL9, CXCL10, PPARG, and NFATC1 as promising candidates related to resistance/susceptibility to Mycobacterium avium subspecies paratuberculosis (MAP). This study reports on the high-resolution DNA methylation landscapes of the jejunum methylome across three conditions (clinical, subclinical, and healthy) in dairy cows. Our investigations integrated different sources of information about DMGs, DEGs, and pathways, enabling us to find nine functional genes that might have potential application in resisting paratuberculosis in dairy cattle.

Doxorubicin induces wide-spread transcriptional changes in the myocardium of hearts distinguishing between mice with preserved and impaired cardiac function

AIMS

Doxorubicin (DOX) is an important drug for the treatment of various tumor entities. However, the occurrence of heart failure limits its application. This study investigated differential gene expression profiles in the left and right ventricles of DOX treated mice with either preserved or impaired myocardial function. We provide new mechanistic insights into the pathophysiology of DOX-induced heart failure and have discovered pathways that counteract DOX-induced cardiotoxicity.

MAIN METHODS

We used in total 48 male mice and applied a chronic low dose DOX administration (5 mg/kg per injection, in total 20 mg/kg over 4 weeks) to induce heart failure. Echocardiographic parameters were evaluated one week after the final dose and mice were separated according to functional parameters into doxorubicin responding and non-responding animals. Post mortem, measurements of reactive oxygen species (ROS) and gene expression profiling was performed in separated right and left hearts.

KEY FINDINGS

We detected significant ROS production in the left heart of the mice in response to DOX treatment, although interestingly, not in the right heart. We found that transcriptional changes differ between right and left heart correlating with the occurrence of myocardial dysfunction.

SIGNIFICANCE

Doxorubicin induces changes in gene expression in the entire heart of animals without necessarily impairing cardiac function. We identified a set of transcripts that are associated with DOX cardiotoxicity. These might represent promising targets to ameliorate DOX-induced heart failure. Moreover, our results emphasize that parameters of left and right heart function should be evaluated during standardized echocardiography in patients undergoing DOX therapy.

A cis-regulatory element promoting increased transcription at low temperature in cultured ectothermic Drosophila cells

Background

Temperature change affects the myriad of concurrent cellular processes in a non-uniform, disruptive manner. While endothermic organisms minimize the challenge of ambient temperature variation by keeping the core body temperature constant, cells of many ectothermic species maintain homeostatic function within a considerable temperature range. The cellular mechanisms enabling temperature acclimation in ectotherms are still poorly understood. At the transcriptional level, the heat shock response has been analyzed extensively. The opposite, the response to sub-optimal temperature, has received lesser attention in particular in animal species. The tissue specificity of transcriptional responses to cool temperature has not been addressed and it is not clear whether a prominent general response occurs. Cis-regulatory elements (CREs), which mediate increased transcription at cool temperature, and responsible transcription factors are largely unknown.

Results

The ectotherm Drosophila melanogaster with a presumed temperature optimum around 25 °C was used for transcriptomic analyses of effects of temperatures at the lower end of the readily tolerated range (14–29 °C). Comparative analyses with adult flies and cell culture lines indicated a striking degree of cell-type specificity in the transcriptional response to cool. To identify potential cis-regulatory elements (CREs) for transcriptional upregulation at cool temperature, we analyzed temperature effects on DNA accessibility in chromatin of S2R+ cells. Candidate cis-regulatory elements (CREs) were evaluated with a novel reporter assay for accurate assessment of their temperature-dependency. Robust transcriptional upregulation at low temperature could be demonstrated for a fragment from the pastrel gene, which expresses more transcript and protein at reduced temperatures. This CRE is controlled by the JAK/STAT signaling pathway and antagonizing activities of the transcription factors Pointed and Ets97D.

Conclusion

Beyond a rich data resource for future analyses of transcriptional control within the readily tolerated range of an ectothermic animal, a novel reporter assay permitting quantitative characterization of CRE temperature dependence was developed. Our identification and functional dissection of the pst_E1 enhancer demonstrate the utility of resources and assay. The functional characterization of this CoolUp enhancer provides initial mechanistic insights into transcriptional upregulation induced by a shift to temperatures at the lower end of the readily tolerated range.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08057-4.