mmquant: how to count multi-mapping reads?

A fusarioid fungus forms mutualistic interactions with poplar trees that resemble ectomycorrhizal symbiosis

Fusarium species, recognised as global priority pathogens, frequently induce severe diseases in crops; however, certain species exhibit alternative symbiotic lifestyles and are either non-pathogenic or endophytic. In this study, we characterised the mutualistic relationship between the eFp isolate of F.pseudograminearum and five poplar species, resulting in formation root structures reminiscent of ectomycorrhizal (ECM) symbiosis. This functional symbiosis is evidenced by enhanced plant growth, reciprocal nutrient exchange, improved nitrogen and phosphorus uptake and upregulation of root sugar transporter gene expression (PtSweet1). Comparative and population genomics confirmed that eFp maintains a structurally similar genome, but exhibits significant divergence from ten conspecific pathogenic isolates. Notably, eFp enhanced the growth of diverse plant lineages (Oryza, Arabidopsis, Pinus and non-vascular liverworts), indicating a near-complete loss of virulence. Although this specialised symbiosis has only been established in vitro, it holds significant value in elucidating the evolutionary track from endophytic to mycorrhizal associations.

The Convergent Evolution of Hummingbird Pollination Results in Repeated Floral Scent Loss Through Gene Downregulation

The repeated evolution of the same trait in different lineages provides powerful natural experiments to study the phenotypic and genotypic predictability of trait gain and loss. A fascinating example is the repeated evolution of hummingbird pollination in plant lineages in the Americas, a widespread and often unidirectional phenomenon. The spiral gingers in the genus Costus are ancestrally bee pollinated, and hummingbird pollination has evolved multiple times independently in the tropical Americas. These pollinator transitions are accompanied by predictable morphological and color changes, but the changes in floral scent have not been described. In this study, we describe the floral scent composition of 30 species of Costus sampled across the phylogeny to understand how floral scent has evolved across the genus with respect to pollinator transitions. We then combine transcriptomics and genomics to identify gene expression differences and gene family evolution associated with pollinator transitions. We show that hummingbird-pollinated species have mostly lost their floral scent, whereas bee-pollinated species exhibit either floral scent maintenance or, in some cases, gains of more diverse scent profiles. We find the floral scent loss appears to be due to gene downregulation rather than pseudogenization. The remarkable consistency of scent loss in hummingbird-pollinated species highlights the shared strong selection pressures experienced by these lineages. Even species with more recent transitions from bee to hummingbird pollination exhibit scent loss, highlighting the rapid breakdown of scent emission following pollinator transitions. This research highlights the capacity for rapid changes when selection pressures are strong through downregulation of floral scent genes.

Accurate allocation of multimapped reads enables regulatory element analysis at repeats

Transposable elements (TEs) and other repetitive regions have been shown to contain gene regulatory elements, including transcription factor binding sites. However, regulatory elements harbored by repeats have proven difficult to characterize using short-read sequencing assays such as ChIP-seq or ATAC-seq. Most regulatory genomics analysis pipelines discard "multimapped" reads that align equally well to multiple genomic locations. Because multimapped reads arise predominantly from repeats, current analysis pipelines fail to detect a substantial portion of regulatory events that occur in repetitive regions. To address this shortcoming, we developed Allo, a new approach to allocate multimapped reads in an efficient, accurate, and user-friendly manner. Allo combines probabilistic mapping of multimapped reads with a convolutional neural network that recognizes the read distribution features of potential peaks, offering enhanced accuracy in multimapping read assignment. Allo also provides read-level output in the form of a corrected alignment file, making it compatible with existing regulatory genomics analysis pipelines and downstream peak-finders. In a demonstration application on CTCF ChIP-seq data, we show that Allo results in the discovery of thousands of new CTCF peaks. Many of these peaks contain the expected cognate motif and/or serve as TAD boundaries. We additionally apply Allo to a diverse collection of ENCODE ChIP-seq data sets, resulting in multiple previously unidentified interactions between transcription factors and repetitive element families. Finally, we show that Allo may be particularly beneficial in identifying ChIP-seq peaks at centromeres, near segmentally duplicated genes, and in younger TEs, enabling new regulatory analyses in these regions.

Intraspecific diploidization of a halophyte root fungus drives heterosis

How organisms respond to environmental stress is a key topic in evolutionary biology. This study focused on the genomic evolution of Laburnicola rhizohalophila, a dark-septate endophytic fungus from roots of a halophyte. Chromosome-level assemblies were generated from five representative isolates from structured subpopulations. The data revealed significant genomic plasticity resulting from chromosomal polymorphisms created by fusion and fission events, known as dysploidy. Analyses of genomic features, phylogenomics, and macrosynteny have provided clear evidence for the origin of intraspecific diploid-like hybrids. Notably, one diploid phenotype stood out as an outlier and exhibited a conditional fitness advantage when exposed to a range of abiotic stresses compared with its parents. By comparing the gene expression patterns in each hybrid parent triad under the four growth conditions, the mechanisms underlying growth vigor were corroborated through an analysis of transgressively upregulated genes enriched in membrane glycerolipid biosynthesis and transmembrane transporter activity. In vitro assays suggested increased membrane integrity and lipid accumulation, as well as decreased malondialdehyde production under optimal salt conditions (0.3 M NaCl) in the hybrid. These attributes have been implicated in salinity tolerance. This study supports the notion that hybridization-induced genome doubling leads to the emergence of phenotypic innovations in an extremophilic endophyte.

Evolution of sex-biased genes in Drosophila species with neo-sex chromosomes: Potential contribution to reducing the sexual conflict

An advantage of sex chromosomes may be the potential to reduce sexual conflict because they provide a basis for selection to operate separately on females and males. However, evaluating the relationship between sex chromosomes and sexual conflict is challenging owing to the difficulty in measuring sexual conflict and substantial divergence between species with and without sex chromosomes. We therefore examined sex-biased gene expression as a proxy for sexual conflict in three sets of Drosophila species with and without young sex chromosomes, the so-called neo-sex chromosomes. In all sets, we detected more sex-biased genes in the species with neo-sex chromosomes than in the species without neo-sex chromosomes in larvae, pupae, and adult somatic tissues but not in gonads. In particular, many unbiased genes became either female- or male-biased after linkage to the neo-sex chromosomes in larvae, despite the low sexual dimorphism. For example, genes involved in metabolism, a key determinant for the rate of development in many animals, were enriched in the genes that acquired sex-biased expression on the neo-sex chromosomes at the larval stage. These genes may be targets of sexually antagonistic selection (i.e., large size and rapid development are selected for in females but selected against in males). These results indicate that acquiring neo-sex chromosomes may have contributed to a reduction in sexual conflict, particularly at the larval stage, in Drosophila..

Multi-omics analysis reveals COVID-19 vaccine induced attenuation of inflammatory responses during breakthrough disease

The immune mechanisms mediating COVID-19 vaccine attenuation of COVID-19 remain undescribed. We conducted comprehensive analyses detailing immune responses to SARS-CoV-2 virus in blood post-vaccination with ChAdOx1 nCoV-19 or a placebo. Samples from randomised placebo-controlled trials (NCT04324606 and NCT04400838) were taken at baseline, onset of COVID-19-like symptoms, and 7 days later, confirming COVID-19 using nucleic amplification test (NAAT test) via real-time PCR (RT-PCR). Serum cytokines were measured with multiplexed immunoassays. The transcriptome was analysed with long, short and small RNA sequencing. We found attenuation of RNA inflammatory signatures in ChAdOx1 nCoV-19 compared with placebo vaccinees and reduced levels of serum proteins associated with COVID-19 severity. KREMEN1, a putative alternative SARS-CoV-2 receptor, was downregulated in placebo compared with ChAdOx1 nCoV-19 vaccinees. Vaccination ameliorates reductions in cell counts across leukocyte populations and platelets noted at COVID-19 onset, without inducing potentially deleterious Th2-skewed immune responses. Multi-omics integration links a global reduction in miRNA expression at COVID-19 onset to increased pro-inflammatory responses at the mRNA level. This study reveals insights into the role of COVID-19 vaccines in mitigating disease severity by abrogating pro-inflammatory responses associated with severe COVID-19, affirming vaccine-mediated benefit in breakthrough infection, and highlighting the importance of clinically relevant endpoints in vaccine evaluation.

Selective gene expression maintains human tRNA anticodon pools during differentiation

Transfer RNAs are essential for translating genetic information into proteins. The human genome contains hundreds of predicted tRNA genes, many in multiple copies. How their expression is regulated to control tRNA repertoires is unknown. Here we combined quantitative tRNA profiling and chromatin immunoprecipitation with sequencing to measure tRNA expression following the differentiation of human induced pluripotent stem cells into neuronal and cardiac cells. We find that tRNA transcript levels vary substantially, whereas tRNA anticodon pools, which govern decoding rates, are more stable among cell types. Mechanistically, RNA polymerase III transcribes a wide range of tRNA genes in human induced pluripotent stem cells but on differentiation becomes constrained to a subset we define as housekeeping tRNAs. This shift is mediated by decreased mTORC1 signalling, which activates the RNA polymerase III repressor MAF1. Our data explain how tRNA anticodon pools are buffered to maintain decoding speed across cell types and reveal that mTORC1 drives selective tRNA expression during differentiation.

The transcriptomics profiling of blood CD4 and CD8 T-cells in narcolepsy type I

Background

Narcolepsy Type I (NT1) is a rare, life-long sleep disorder arising as a consequence of the extensive destruction of orexin-producing hypothalamic neurons. The mechanisms involved in the destruction of orexin neurons are not yet elucidated but the association of narcolepsy with environmental triggers and genetic susceptibility (strong association with the HLA, TCRs and other immunologically-relevant loci) implicates an immuno-pathological process. Several studies in animal models and on human samples have suggested that T-cells are the main pathogenic culprits.

Methods

RNA sequencing was performed on four CD4 and CD8 T-cell subsets (naive, effector, effector memory and central memory) sorted by flow cytometry from peripheral blood mononuclear cells (PBMCs) of NT1 patients and HLA-matched healthy donors as well as (age- and sex-) matched individuals suffering from other sleep disorders (OSD). The RNAseq analysis was conducted by comparing the transcriptome of NT1 patients to that of healthy donors and other sleep disorder patients (collectively referred to as the non-narcolepsy controls) in order to identify NT1-specific genes and pathways.

Results

We determined NT1-specific differentially expressed genes, several of which are involved in tubulin arrangement found in CD4 (TBCB, CCT5, EML4, TPGS1, TPGS2) and CD8 (TTLL7) T cell subsets, which play a role in the immune synapse formation and TCR signaling. Furthermore, we identified genes (GZMB, LTB in CD4 T-cells and NLRP3, TRADD, IL6, CXCR1, FOXO3, FOXP3 in CD8 T-cells) and pathways involved in various aspects of inflammation and inflammatory response. More specifically, the inflammatory profile was identified in the "naive" subset of CD4 and CD8 T-cell.

Conclusion

We identified NT1-specific differentially expressed genes, providing a cell-type and subset specific catalog describing their functions in T-cells as well as their potential involvement in NT1. Several genes and pathways identified are involved in the formation of the immune synapse and TCR activation as well as inflammation and the inflammatory response. An inflammatory transcriptomic profile was detected in both "naive" CD4 and CD8 T-cell subsets suggesting their possible involvement in the development or progression of the narcoleptic process.

Methylated guanosine and uridine modifications in S. cerevisiae mRNAs modulate translation elongation

Chemical modifications to protein encoding messenger RNAs (mRNAs) influence their localization, translation, and stability within cells. Over 15 different types of mRNA modifications have been observed by sequencing and liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) approaches. While LC-MS/MS is arguably the most essential tool available for studying analogous protein post-translational modifications, the high-throughput discovery and quantitative characterization of mRNA modifications by LC-MS/MS has been hampered by the difficulty of obtaining sufficient quantities of pure mRNA and limited sensitivities for modified nucleosides. We have overcome these challenges by improving the mRNA purification and LC-MS/MS pipelines. The methodologies we developed result in no detectable non-coding RNA modifications signals in our purified mRNA samples, quantify 50 ribonucleosides in a single analysis, and provide the lowest limit of detection reported for ribonucleoside modification LC-MS/MS analyses. These advancements enabled the detection and quantification of 13 S. cerevisiae mRNA ribonucleoside modifications and reveal the presence of four new S. cerevisiae mRNA modifications at low to moderate levels (1-methyguanosine, N2-methylguanosine, N2,N2-dimethylguanosine, and 5-methyluridine). We identified four enzymes that incorporate these modifications into S. cerevisiae mRNAs (Trm10, Trm11, Trm1, and Trm2, respectively), though our results suggest that guanosine and uridine nucleobases are also non-enzymatically methylated at low levels. Regardless of whether they are incorporated in a programmed manner or as the result of RNA damage, we reasoned that the ribosome will encounter the modifications that we detect in cells. To evaluate this possibility, we used a reconstituted translation system to investigate the consequences of modifications on translation elongation. Our findings demonstrate that the introduction of 1-methyguanosine, N2-methylguanosine and 5-methyluridine into mRNA codons impedes amino acid addition in a position dependent manner. This work expands the repertoire of nucleoside modifications that the ribosome must decode in S. cerevisiae. Additionally, it highlights the challenge of predicting the effect of discrete modified mRNA sites on translation de novo because individual modifications influence translation differently depending on mRNA sequence context.

A Genome-Wide View of the Transcriptome Dynamics of Fresh-Cut Potato Tubers

Fresh fruits and vegetable products are easily perishable during postharvest handling due to enzymatic browning reactions. This phenomenon has contributed to a significant loss of food. To reveal the physiological changes in fresh-cut potato tubers at the molecular level, a transcriptome analysis of potato tubers after cutting was carried out. A total of 10,872, 10,449, and 11,880 differentially expressed genes (DEGs) were identified at 4 h, 12 h and 24 h after cutting, respectively. More than 87.5% of these DEGs were classified into the categories of biological process (BP) and molecular function (MF) based on Gene Ontology (GO) analysis. There was a difference in the response to cutting at different stages after the cutting of potato tubers. The genes related to the phenol and fatty biosynthesis pathways, which are responsible for enzymatic browning and wound healing in potato tubers, were significantly enriched at 0-24 h after cutting. Most genes related to the enzymatic browning of potato tubers were up-regulated in response to cut-wounding. Plant hormone biosynthesis, signal molecular biosynthesis and transduction-related genes, such as gibberelin (GA), cytokinin (CK), ethylene (ET), auxin (IAA), jasmonic acid (JA), salicylic (SA), and Respiratory burst oxidase (Rboh) significantly changed at the early stage after cutting. In addition, the transcription factors involved in the wound response were the most abundant at the early stage after cutting. The transcription factor with the greatest response to injury was MYB, followed by AP2-EREBP, C3H and WRKY. This study revealed the physiological changes at the molecular level of fresh-cut potato tubers after cutting. This information is needed for developing a better approach to enhancing the postharvest shelf life of fresh processed potato and the breeding of potato plants that are resistant to enzymatic browning.

High-quality read-based phasing of cystic fibrosis cohort informs genetic understanding of disease modification

Phasing of heterozygous alleles is critical for interpretation of cis-effects of disease-relevant variation. We sequenced 477 individuals with cystic fibrosis (CF) using linked-read sequencing, which display an average phase block N50 of 4.39 Mb. We use these samples to construct a graph representation of CFTR haplotypes, demonstrating its utility for understanding complex CF alleles. These are visualized in a Web app, CFTbaRcodes, that enables interactive exploration of CFTR haplotypes present in this cohort. We perform fine-mapping and phasing of the chr7q35 trypsinogen locus associated with CF meconium ileus, an intestinal obstruction at birth associated with more severe CF outcomes and pancreatic disease. A 20-kb deletion polymorphism and a PRSS2 missense variant p.Thr8Ile (rs62473563) are shown to independently contribute to meconium ileus risk (p = 0.0028, p = 0.011, respectively) and are PRSS2 pancreas eQTLs (p = 9.5 × 10⁻⁷ and p = 1.4 × 10⁻⁴, respectively), suggesting the mechanism by which these polymorphisms contribute to CF. The phase information from linked reads provides a putative causal explanation for variation at a CF-relevant locus, which also has implications for the genetic basis of non-CF pancreatitis, to which this locus has been reported to contribute.

srnaMapper: an optimal mapping tool for sRNA-Seq reads

Background

Sequencing is the key method to study the impact of short RNAs, which include micro RNAs, tRNA-derived RNAs, and piwi-interacting RNA, among others. The first step to make use of these reads is to map them to a genome. Existing mapping tools have been developed for long RNAs in mind, and, so far, no tool has been conceived for short RNAs. However, short RNAs have several distinctive features which make them different from messenger RNAs: they are shorter, they are often redundant, they can be produced by duplicated loci, and they may be edited at their ends.

Results

In this work, we present a new tool, srnaMapper, that exhaustively maps these reads with all these features in mind, and is most efficient when applied to reads no longer than 50 base pairs. We show, on several datasets, that srnaMapper is very efficient considering computation time and edition error handling: it retrieves all the hits, with arbitrary number of errors, in time comparable with non-exhaustive tools.

Environmental response in gene expression and DNA methylation reveals factors influencing the adaptive potential of Arabidopsis lyrata

Understanding what factors influence plastic and genetic variation is valuable for predicting how organisms respond to changes in the selective environment. Here, using gene expression and DNA methylation as molecular phenotypes, we study environmentally induced variation among Arabidopsis lyrata plants grown at lowland and alpine field sites. Our results show that gene expression is highly plastic, as many more genes are differentially expressed between the field sites than between populations. These environmentally responsive genes evolve under strong selective constraint - the strength of purifying selection on the coding sequence is high, while the rate of adaptive evolution is low. We find, however, that positive selection on cis-regulatory variants has likely contributed to the maintenance of genetically variable environmental responses, but such variants segregate only between distantly related populations. In contrast to gene expression, DNA methylation at genic regions is largely insensitive to the environment, and plastic methylation changes are not associated with differential gene expression. Besides genes, we detect environmental effects at transposable elements (TEs): TEs at the high-altitude field site have higher expression and methylation levels, suggestive of a broad-scale TE activation. Compared to the lowland population, plants native to the alpine environment harbor an excess of recent TE insertions, and we observe that specific TE families are enriched within environmentally responsive genes. Our findings provide insight into selective forces shaping plastic and genetic variation. We also highlight how plastic responses at TEs can rapidly create novel heritable variation in stressful conditions.

Rta is the principal activator of Epstein-Barr virus epithelial lytic transcription

The transition from latent Epstein-Barr virus (EBV) infection to lytic viral replication is mediated by the viral transcription factors Rta and Zta. Although both are required for virion production, dissecting the specific roles played by Rta and Zta is challenging because they induce each other's expression. To circumvent this, we constructed an EBV mutant deleted for the genes encoding Rta and Zta (BRLF1 and BZLF1, respectively) in the Akata strain BACmid. This mutant, termed EBVΔRZ, was used to infect several epithelial cell lines, including telomerase-immortalized normal oral keratinocytes, a highly physiologic model of EBV epithelial cell infection. Using RNA-seq, we determined the gene expression induced by each viral transactivator. Surprisingly, Zta alone only induced expression of the lytic origin transcripts BHLF1 and LF3. In contrast, Rta activated the majority of EBV early gene transcripts. As expected, Zta and Rta were both required for expression of late gene transcripts. Zta also cooperated with Rta to enhance a subset of early gene transcripts (Rtasynergy transcripts) that Zta was unable to activate when expressed alone. Interestingly, Rta and Zta each cooperatively enhanced the other's binding to EBV early gene promoters, but this effect was not restricted to promoters where synergy was observed. We demonstrate that Zta did not affect Rtasynergy transcript stability, but increased Rtasynergy gene transcription despite having no effect on their transcription when expressed alone. Our results suggest that, at least in epithelial cells, Rta is the dominant transactivator and that Zta functions primarily to support DNA replication and co-activate a subset of early promoters with Rta. This closely parallels the arrangement in KSHV where ORF50 (Rta homolog) is the principal activator of lytic transcription and K8 (Zta homolog) is required for DNA replication at oriLyt.

Alpinetin promotes hair regeneration via activating hair follicle stem cells

Background

Alopecia affects millions of individuals globally, with hair loss becoming more common among young people.  Various traditional Chinese medicines (TCM) have been used clinically for treating alopecia, however, the effective compounds and underlying mechanism are less known. We sought to investigate the effect of Alpinetin (AP), a compound extracted from Fabaceae and Zingiberaceae herbs, in hair regeneration.

Methods

Animal model for hair regeneration was mimicked by depilation in C57BL/6J mice. The mice were then topically treated with 3 mg/ml AP, minoxidil as positive control (PC), or solvent ethanol as vehicle control (VC) on the dorsal skin. Skin color changes which reflected the hair growth stages were monitored and pictured, along with H&E staining and hair shaft length measurement. RNA-seq analysis combined with immunofluorescence staining and qPCR analysis were used for mechanism study. Meanwhile, Gli1^CreERT2; R26R^tdTomato and Lgr5^{EGFP−CreERT2}; R26R^tdTomato transgenic mice were used to monitor the activation and proliferation of Gli1+ and Lgr5+ HFSCs after treatment. Furthermore, the toxicity of AP was tested in keratinocytes and fibroblasts from both human and mouse skin to assess the safety.

Results

When compared to minoxidil-treated and vehicle-treated control mice, topical application of AP promoted anagen initiation and delayed catagen entry, resulting in a longer anagen phase and hair shaft length. Mechanistically, RNA-seq analysis combined with immunofluorescence staining of Lef1 suggested that Lgr5+ HFSCs in lower bulge were activated by AP via Wnt signaling. Other HFSCs, including K15+, Lef1+, and Gli1+ cells, were also promoted into proliferating upon AP treatment. In addition, AP inhibited cleaved caspase 3-dependent apoptosis at the late anagen stage to postpone regression of hair follicles. More importantly, AP showed no cytotoxicity in keratinocytes and fibroblasts from both human and mouse skin.

Conclusion

This study clarified the effect of AP in promoting hair regeneration by activating HFSCs via Wnt signaling. Our findings may contribute to the development of a new generation of pilatory that is more efficient and less cytotoxic for treating alopecia.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13020-022-00619-2.

Accurate Quantification of Overlapping Herpesvirus Transcripts from RNA Sequencing Data

Herpesviruses employ extensive bidirectional transcription of overlapping genes to overcome length constraints on their gene product repertoire. As a consequence, many lytic transcripts cannot be measured individually by reverse transcription-quantitative PCR (RT-qPCR) or conventional RNA sequencing (RNA-seq) analysis. A. G. Bruce, S. Barcy, T. DiMaio, E. Gan, et al. (Pathogens 6:11, 2017, https://doi.org/10.3390/pathogens6010011) have proposed an approximation method using unique coding sequences (UCDS) to estimate lytic gene abundance from Kaposi's sarcoma-associated herpesvirus (KSHV) RNA-seq data. Although UCDS has been widely employed, its accuracy, to our knowledge, has never been rigorously validated for any herpesvirus. In this study, we use cap analysis of gene expression sequencing (CAGE-seq) as a gold-standard to determine the accuracy of UCDS for estimating Epstein-Barr virus (EBV) lytic gene expression levels from RNA-seq data. We also introduce the Unique TranScript (UTS) method, which, like UCDS, estimates transcript abundance from changes in mean RNA-seq read depth. UTS is distinguished by its use of empirically determined 5' and 3' transcript ends rather than coding sequence annotations. Compared to conventional read assignment, both UCDS and UTS improved the accuracy of quantitation of overlapping genes, with UTS giving the most-accurate results. The UTS method discards fewer reads and may be advantageous for experiments with less sequencing depth. UTS is compatible with any aligner and, unlike isoform-aware alignment methods, can be implemented on a laptop computer. Our findings demonstrate that the accuracy achieved by complex and expensive techniques such as CAGE-seq can be approximated using conventional short-read RNA-seq data when read assignment methods address transcript overlap. Although our study focuses on EBV transcription, the UTS method should be applicable across all herpesviruses as well as to other genomes with extensively overlapping transcriptomes. IMPORTANCE Many viruses employ extensively overlapping transcript structures. This complexity makes it difficult to quantify gene expression by using conventional methods, including RNA-seq. Although high-throughput techniques that overcome these limitations exist, they are complex, expensive, and scarce in the herpesvirus literature relative to short-read RNA-seq. Here, using Epstein-Barr virus (EBV) as a model, we demonstrate that conventional RNA-seq analysis methods fail to accurately quantify the abundances of many overlapping transcripts. We further show that the previously described Unique CoDing Sequence (UCDS) method and our Unique TranScript (UTS) method greatly improve the accuracy of EBV lytic gene measurements obtained from RNA-seq data. The UTS method has the advantages of discarding fewer reads and being implementable on a laptop computer. Although this study focuses on EBV, the UCDS and UTS methods should be applicable across herpesviruses and for other viruses that make extensive use of overlapping transcription.

MGcount: a total RNA-seq quantification tool to address multi-mapping and multi-overlapping alignments ambiguity in non-coding transcripts

BACKGROUND

Total-RNA sequencing (total-RNA-seq) allows the simultaneous study of both the coding and the non-coding transcriptome. Yet, computational pipelines have traditionally focused on particular biotypes, making assumptions that are not fullfilled by total-RNA-seq datasets. Transcripts from distinct RNA biotypes vary in length, biogenesis, and function, can overlap in a genomic region, and may be present in the genome with a high copy number. Consequently, reads from total-RNA-seq libraries may cause ambiguous genomic alignments, demanding for flexible quantification approaches.

RESULTS

Here we present Multi-Graph count (MGcount), a total-RNA-seq quantification tool combining two strategies for handling ambiguous alignments. First, MGcount assigns reads hierarchically to small-RNA and long-RNA features to account for length disparity when transcripts overlap in the same genomic position. Next, MGcount aggregates RNA products with similar sequences where reads systematically multi-map using a graph-based approach. MGcount outputs a transcriptomic count matrix compatible with RNA-sequencing downstream analysis pipelines, with both bulk and single-cell resolution, and the graphs that model repeated transcript structures for different biotypes. The software can be used as a python module or as a single-file executable program.

CONCLUSIONS

MGcount is a flexible total-RNA-seq quantification tool that successfully integrates reads that align to multiple genomic locations or that overlap with multiple gene features. Its approach is suitable for the simultaneous estimation of protein-coding, long non-coding and small non-coding transcript concentration, in both precursor and processed forms. Both source code and compiled software are available at https://github.com/hitaandrea/MGcount .

Integrated physiological, transcriptome and metabolome analyses of the hepatopancreas of the female swimming crab Portunus trituberculatus under ammonia exposure

Ammonia is a common environmental pollutant in aquatic ecosystem and is also a significant concern in closed aquaculture systems. The threat of ammonia has been increasing with rising anthropogenic activities including intensified aquaculture. In this study, we aimed to investigate ammonia toxicity and metabolism mechanisms in the hepatopancreas, a major organ for Vitellogenin (Vtg) synthesis and defending against ammonia stress, of female swimming crab Portunus trituberculatus which is an important fishery and aquaculture species, by integrating physiological, transcriptome and metabolome analyses. The results revealed that ammonia exposure (10 mg/L, an environmentally relevant concentration) resulted in a remarkable reduction in vtg expression and depression of multiple signaling pathways for reproductive regulators including methyl farnesoate, ecdysone and neuroparsin, demonstrating for the first time that ammonia impairs swimming crab female reproduction. In addition, a number of important genes and metabolites in glycolysis, the Krebs cycle, fatty acid β-oxidation and synthesis were significantly downregulated, indicating that changes in ammonia levels lead to a general depression of energy metabolism in hepatopancreas. After ammonia exposure, an increased level of urea and a reduction of amino acid catabolism were observed in hepatopancreas, suggesting that urea cycle was utilized to biotransform ammonia, and amino acid catabolism was decreased to reduce endogenous ammonia generation. Furthermore, antioxidant systems were altered following ammonia exposure, which was accompanied by proteins and lipid oxidations, as well as cellular apoptosis. These results indicate that ammonia leads to metabolic suppression, oxidative stress and apoptosis in P. trituberculatus hepatopancreas. The findings improve the understanding for the mechanisms of ammonia toxicity and metabolism in P. trituberculatus, and provide valuable information for assessing potential ecological risk of environmental ammonia and improving aquaculture management.

Identification and Validation of Reference Genes in Clostridium beijerinckii NRRL B-598 for RT-qPCR Using RNA-Seq Data

Gene expression analysis through reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) depends on correct data normalization by reference genes with stable expression. Although Clostridium beijerinckii NRRL B-598 is a promising Gram-positive bacterium for the industrial production of biobutanol, validated reference genes have not yet been reported. In this study, we selected 160 genes with stable expression based on an RNA sequencing (RNA-Seq) data analysis, and among them, seven genes (zmp, rpoB1, rsmB, greA, rpoB2, topB2, and rimO) were selected for experimental validation by RT-qPCR and gene ontology (GO) enrichment analysis. According to statistical analyses, zmp and greA were the most stable and suitable reference genes for RT-qPCR normalization. Furthermore, our methodology can be useful for selection of the reference genes in other strains of C. beijerinckii and it also suggests that the RNA-Seq data can be used for the initial selection of novel reference genes, however, their validation is required.

Slowpoke functions in circadian output cells to regulate rest:activity rhythms

The circadian system produces ~24-hr oscillations in behavioral and physiological processes to ensure that they occur at optimal times of day and in the correct temporal order. At its core, the circadian system is composed of dedicated central clock neurons that keep time through a cell-autonomous molecular clock. To produce rhythmic behaviors, time-of-day information generated by clock neurons must be transmitted across output pathways to regulate the downstream neuronal populations that control the relevant behaviors. An understanding of the manner through which the circadian system enacts behavioral rhythms therefore requires the identification of the cells and molecules that make up the output pathways. To that end, we recently characterized the Drosophila pars intercerebralis (PI) as a major circadian output center that lies downstream of central clock neurons in a circuit controlling rest:activity rhythms. We have conducted single-cell RNA sequencing (scRNAseq) to identify potential circadian output genes expressed by PI cells, and used cell-specific RNA interference (RNAi) to knock down expression of ~40 of these candidate genes selectively within subsets of PI cells. We demonstrate that knockdown of the slowpoke (slo) potassium channel in PI cells reliably decreases circadian rest:activity rhythm strength. Interestingly, slo mutants have previously been shown to have aberrant rest:activity rhythms, in part due to a necessary function of slo within central clock cells. However, rescue of slo in all clock cells does not fully reestablish behavioral rhythms, indicating that expression in non-clock neurons is also necessary. Our results demonstrate that slo exerts its effects in multiple components of the circadian circuit, including PI output cells in addition to clock neurons, and we hypothesize that it does so by contributing to the generation of daily neuronal activity rhythms that allow for the propagation of circadian information throughout output circuits.

Quantitative analysis of the splice variants expressed by the major hepatitis B virus genotypes

Hepatitis B virus (HBV) is a major human pathogen that causes liver diseases. The main HBV RNAs are unspliced transcripts that encode the key viral proteins. Recent studies have shown that some of the HBV spliced transcript isoforms are predictive of liver cancer, yet the roles of these spliced transcripts remain elusive. Furthermore, there are nine major HBV genotypes common in different regions of the world, these genotypes may express different spliced transcript isoforms. To systematically study the HBV splice variants, we transfected human hepatoma cells, Huh7, with four HBV genotypes (A2, B2, C2 and D3), followed by deep RNA-sequencing. We found that 13-28 % of HBV RNAs were splice variants, which were reproducibly detected across independent biological replicates. These comprised 6 novel and 10 previously identified splice variants. In particular, a novel, singly spliced transcript was detected in genotypes A2 and D3 at high levels. The biological relevance of these splice variants was supported by their identification in HBV-positive liver biopsy and serum samples, and in HBV-infected primary human hepatocytes. Interestingly the levels of HBV splice variants varied across the genotypes, but the spliced pregenomic RNA SP1 and SP9 were the two most abundant splice variants. Counterintuitively, these singly spliced SP1 and SP9 variants had a suboptimal 5' splice site, supporting the idea that splicing of HBV RNAs is tightly controlled by the viral post-transcriptional regulatory RNA element.

Probably Correct: Rescuing Repeats with Short and Long Reads

Ever since the introduction of high-throughput sequencing following the human genome project, assembling short reads into a reference of sufficient quality posed a significant problem as a large portion of the human genome-estimated 50-69%-is repetitive. As a result, a sizable proportion of sequencing reads is multi-mapping, i.e., without a unique placement in the genome. The two key parameters for whether or not a read is multi-mapping are the read length and genome complexity. Long reads are now able to span difficult, heterochromatic regions, including full centromeres, and characterize chromosomes from "telomere to telomere". Moreover, identical reads or repeat arrays can be differentiated based on their epigenetic marks, such as methylation patterns, aiding in the assembly process. This is despite the fact that long reads still contain a modest percentage of sequencing errors, disorienting the aligners and assemblers both in accuracy and speed. Here, I review the proposed and implemented solutions to the repeat resolution and the multi-mapping read problem, as well as the downstream consequences of reference choice, repeat masking, and proper representation of sex chromosomes. I also consider the forthcoming challenges and solutions with regards to long reads, where we expect the shift from the problem of repeat localization within a single individual to the problem of repeat positioning within pangenomes.

Transcriptomic profiling of Brassica napus responses to Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic bacterial pathogen of plants. Unlike the well-characterized plant defense responses to highly adapted bacterial phytopathogens, little is known about plant response to P. aeruginosa infection. In this study, we examined the Brassica napus (canola) tissue-specific response to P. aeruginosa infection using RNA sequencing. Transcriptomic analysis of canola seedlings over a 5 day P. aeruginosa infection revealed that many molecular processes involved in plant innate immunity were up-regulated, whereas photosynthesis was down-regulated. Phytohormones control many vital biological processes within plants, including growth and development, senescence, seed setting, fruit ripening, and innate immunity. The three main phytohormones involved in plant innate immunity are salicylic acid (SA), jasmonic acid (JA), and ethylene (ET). Many bacterial pathogens have evolved multiple strategies to manipulate these hormone responses in order to infect plants successfully. Interestingly, gene expression within all three phytohormone (SA, JA, and ET) signaling pathways was up-regulated in response to P. aeruginosa infection. This study identified a unique plant hormone response to the opportunistic bacterial pathogen P. aeruginosa infection.

A clinical and multi‑omics study of Van der Woude syndrome in three generations of a Chinese family

Previous studies have suggested that pathogenic variants in interferon regulatoryse factor 6 (IRF6) can account for almost 70% of familial Van der Woude Syndrome (VWS) cases. However, gene modifiers that account for the phenotypic variability of IRF6 in the context of VWS remain poorly characterized. The aim of this study was to report a family with VWS with variable expressivity and to identify the genetic cause. A 4-month-old boy initially presented with cleft palate and bilateral lower lip pits. Examination of his family history identified similar, albeit milder, clinical features in another four family members, including bilateral lower lip pits and/or hypodontia. Peripheral blood samples of eight members in this three-generation family were subsequently collected, and whole-exome sequencing was performed to detect pathogenic variants. A heterozygous missense IRF6 variant with a c.1198C>T change in exon 9 (resulting in an R400W change at the amino acid level) was detected in five affected subjects, but not in the other three unaffected subjects. Moreover, subsequent structural analysis was indicative of damaged stability to the structure in the mutant IRF protein. Whole-transcriptome sequencing, expression analysis and Gene Ontology enrichment analysis were conducted on two groups of patients with phenotypic diversity from the same family. These analyses identified significant differentially expressed genes and enriched pathways in these two groups. Altogether, these findings provide insight into the mechanism underlying the variable expressivity of VWS.

Handling multi-mapped reads in RNA-seq

Many eukaryotic genomes harbour large numbers of duplicated sequences, of diverse biotypes, resulting from several mechanisms including recombination, whole genome duplication and retro-transposition. Such repeated sequences complicate gene/transcript quantification during RNA-seq analysis due to reads mapping to more than one locus, sometimes involving genes embedded in other genes. Genes of different biotypes have dissimilar levels of sequence duplication, with long-noncoding RNAs and messenger RNAs sharing less sequence similarity to other genes than biotypes encoding shorter RNAs. Many strategies have been elaborated to handle these multi-mapped reads, resulting in increased accuracy in gene/transcript quantification, although separate tools are typically used to estimate the abundance of short and long genes due to their dissimilar characteristics. This review discusses the mechanisms leading to sequence duplication, the biotypes affected, the computational strategies employed to deal with multi-mapped reads and the challenges that still remain to be overcome.

Multimapping confounds ribosome profiling analysis: A case-study of the Hsp90 molecular chaperone

Ribosome profiling (Ribo-seq) can potentially provide detailed information about ribosome position on transcripts and estimates of protein translation levels in vivo. Hsp90 chaperones, which play a critical role in stress tolerance, have characteristic patterns of differential expression under nonstressed and heat shock conditions. By analyzing published Ribo-seq data for the Hsp90 chaperones in S. cerevisiae, we find wide-ranging artifacts originating from "multimapping" reads (reads that cannot be uniquely assigned to one position), which constitute ~25% of typical S. cerevisiae Ribo-seq datasets and ~80% of the reads from HEK293 cells. Estimates of Hsp90 protein production as determined by Ribo-seq are reproducible but not robust, with inferred expression levels that can change 10-fold depending on how multimapping reads are processed. The differential expression of Hsp90 chaperones under nonstressed and heat shock conditions creates artificial peaks and valleys in their ribosome profiles that give a false impression of regulated translational pausing. Indeed, we find that multimapping can even create an appearance of reproducibility to the shape of the Hsp90 ribosome profiles from biological replicates. Adding further complexity, this artificial reproducibility is dependent on the computational method used to construct the ribosome profile. Given the ubiquity of multimapping reads in Ribo-seq experiments and the complexity of artifacts associated with multimapping, we developed a publicly available computational tool to identify transcripts most at risk for multimapping artifacts. In doing so, we identify biological pathways that are enriched in multimapping transcripts, meaning that particular biological pathways will be highly susceptible to multimapping artifacts.

Identification and Quantification of Modified Nucleosides in Saccharomyces cerevisiae mRNAs

Post-transcriptional modifications to messenger RNAs (mRNAs) have the potential to alter the biological function of this important class of biomolecules. The study of mRNA modifications is a rapidly emerging field, and the full complement of chemical modifications in mRNAs is not yet established. We sought to identify and quantify the modifications present in yeast mRNAs using an ultra-high performance liquid chromatography tandem mass spectrometry method to detect 40 nucleoside variations in parallel. We observe six modified nucleosides with high confidence in highly purified mRNA samples (N7-methylguanosine, N6-methyladenosine, 2'-O-methylguanosine, 2'-O-methylcytidine, N4-acetylcytidine, and 5-formylcytidine) and identify the yeast protein responsible for N4-acetylcytidine incorporation in mRNAs (Rra1). In addition, we find that mRNA modification levels change in response to heat shock, glucose starvation, and/or oxidative stress. This work expands the repertoire of potential chemical modifications in mRNAs and highlights the value of integrating mass spectrometry tools in the mRNA modification discovery and characterization pipeline.

Sequencing XMET genes to promote genotype-guided risk assessment and precision medicine

High-throughput next generation sequencing (NGS) is a shotgun approach applied in a parallel fashion by which the genome is fragmented and sequenced through small pieces and then analyzed either by aligning to a known reference genome or by de novo assembly without reference genome. This technology has led researchers to conduct an explosion of sequencing related projects in multidisciplinary fields of science. However, due to the limitations of sequencing-based chemistry, length of sequencing reads and the complexity of genes, it is difficult to determine the sequences of some portions of the human genome, leaving gaps in genomic data that frustrate further analysis. Particularly, some complex genes are difficult to be accurately sequenced or mapped because they contain high GC-content and/or low complexity regions, and complicated pseudogenes, such as the genes encoding xenobiotic metabolizing enzymes and transporters (XMETs). The genetic variants in XMET genes are critical to predicate inter-individual variability in drug efficacy, drug safety and susceptibility to environmental toxicity. We summarized and discussed challenges, wet-lab methods, and bioinformatics algorithms in sequencing "complex" XMET genes, which may provide insightful information in the application of NGS technology for implementation in toxicogenomics and pharmacogenomics.

Chronic circadian misalignment results in reduced longevity and large-scale changes in gene expression in Drosophila

BACKGROUND

Circadian clocks are found in nearly all organisms, from bacteria to mammals, and ensure that behavioral and physiological processes occur at optimal times of day and in the correct temporal order. It is becoming increasingly clear that chronic circadian misalignment (CCM), such as occurs in shift workers or as a result of aberrant sleeping and eating schedules common to modern society, has profound metabolic and cognitive consequences, but the proximate mechanisms connecting CCM with reduced organismal health are unknown. Furthermore, it has been difficult to disentangle whether the health effects are directly induced by misalignment or are secondary to the alterations in sleep and activity levels that commonly occur with CCM. Here, we investigated the consequences of CCM in the powerful model system of the fruit fly, Drosophila melanogaster. We subjected flies to daily 4-h phase delays in the light-dark schedule and used the Drosophila Activity Monitoring (DAM) system to continuously track locomotor activity and sleep while simultaneously monitoring fly lifespan.

RESULTS

Consistent with previous results, we find that exposing flies to CCM leads to a ~ 15% reduction in median lifespan in both male and female flies. Importantly, we demonstrate that the reduced longevity occurs independent of changes in overall sleep or activity. To uncover potential molecular mechanisms of CCM-induced reduction in lifespan, we conducted whole body RNA-sequencing to assess differences in gene transcription between control and misaligned flies. CCM caused progressive, large-scale changes in gene expression characterized by upregulation of genes involved in response to toxic substances, aging and oxidative stress, and downregulation of genes involved in regulation of development and differentiation, gene expression and biosynthesis.

CONCLUSIONS

Many of these gene expression changes mimic those that occur during natural aging, consistent with the idea that CCM results in premature organismal decline, however, we found that genes involved in lipid metabolism are overrepresented among those that are differentially regulated by CCM and aging. This category of genes is also among the earliest to exhibit CCM-induced changes in expression, thus highlighting altered lipid metabolism as a potentially important mediator of the negative health consequences of CCM.

Transcription profiling of butanol producer Clostridium beijerinckii NRRL B-598 using RNA-Seq

BACKGROUND

Thinning supplies of natural resources increase attention to sustainable microbial production of bio-based fuels. The strain Clostridium beijerinckii NRRL B-598 is a relatively well-described butanol producer regarding its genotype and phenotype under various conditions. However, a link between these two levels, lying in the description of the gene regulation mechanisms, is missing for this strain, due to the lack of transcriptomic data.

RESULTS

In this paper, we present a transcription profile of the strain over the whole fermentation using an RNA-Seq dataset covering six time-points with the current highest dynamic range among solventogenic clostridia. We investigated the accuracy of the genome sequence and particular genome elements, including pseudogenes and prophages. While some pseudogenes were highly expressed, all three identified prophages remained silent. Furthermore, we identified major changes in the transcriptional activity of genes using differential expression analysis between adjacent time-points. We identified functional groups of these significantly regulated genes and together with fermentation and cultivation kinetics captured using liquid chromatography and flow cytometry, we identified basic changes in the metabolism of the strain during fermentation. Interestingly, C. beijerinckii NRRL B-598 demonstrated different behavior in comparison with the closely related strain C. beijerinckii NCIMB 8052 in the latter phases of cultivation.

CONCLUSIONS

We provided a complex analysis of the C. beijerinckii NRRL B-598 fermentation profile using several technologies, including RNA-Seq. We described the changes in the global metabolism of the strain and confirmed the uniqueness of its behavior. The whole experiment demonstrated a good reproducibility. Therefore, we will be able to repeat the experiment under selected conditions in order to investigate particular metabolic changes and signaling pathways suitable for following targeted engineering.