Common occurrence of hotspots of single strand DNA breaks at transcriptional start sites.

BACKGROUND: We recently developed two high-resolution methods for genome-wide mapping of two prominent types of DNA damage, single-strand DNA breaks (SSBs) and abasic (AP) sites and found highly complex and non-random patterns of these lesions in mammalian genomes. One salient feature of SSB and AP sites was the existence of single-nucleotide hotspots for both lesions. RESULTS: In this work, we show that SSB hotspots are enriched in the immediate vicinity of transcriptional start sites (TSSs) in multiple normal mammalian tissues, however the magnitude of enrichment varies significantly with tissue type and appears to be limited to a subset of genes. SSB hotspots around TSSs are enriched on the template strand and associate with higher expression of the corresponding genes. Interestingly, SSB hotspots appear to be at least in part generated by the base-excision repair (BER) pathway from the AP sites. CONCLUSIONS: Our results highlight complex relationship between DNA damage and regulation of gene expression and suggest an exciting possibility that SSBs at TSSs might function as sensors of DNA damage to activate genes important for DNA damage response.

Tenascin C as a novel zinc finger protein 750 target regulating the immunogenicity via DNA damage in lung squamous cell carcinoma.

Modulation of DNA damage repair in lung squamous cell carcinoma (LUSC) can result in the generation of neoantigens and heightened immunogenicity. Therefore, understanding DNA damage repair mechanisms holds significant clinical relevance for identifying targets for immunotherapy and devising therapeutic strategies. Our research has unveiled that the tumor suppressor zinc finger protein 750 (ZNF750) in LUSC binds to the promoter region of tenascin C (TNC), leading to reduced TNC expression. This modulation may impact the malignant behavior of tumor cells and is associated with patient prognosis. Additionally, single-cell RNA sequencing (scRNA-seq) of LUSC tissues has demonstrated an inverse correlation between ZNF750/TNC expression levels and immunogenicity. Manipulation of the ZNF750-TNC axis in vitro within LUSC cells has shown differential sensitivity to CD8+ cells, underscoring its pivotal role in regulating cellular immunogenicity. Further transcriptome sequencing analysis, DNA damage repair assay, and single-strand break analyses have revealed the involvement of the ZNF750-TNC axis in determining the preference for homologous recombination (HR) repair or non-homologous end joining (NHEJ) repair of DNA damage. with involvement of the Hippo/ERK signaling pathway. In summary, this study sheds light on the ZNF750-TNC axis's role in DNA damage repair regulation in LUSC, laying a groundwork for future translational research in immune cell therapy for LUSC.

Evidence for widespread existence of functional novel and non-canonical human transcripts.

BACKGROUND: Fraction of functional sequence in the human genome remains a key unresolved question in Biology and the subject of vigorous debate. While a plethora of studies have connected a significant fraction of human DNA to various biochemical processes, the classical definition of function requires evidence of effects on cellular or organismal fitness that such studies do not provide. Although multiple high-throughput reverse genetics screens have been developed to address this issue, they are limited to annotated genomic elements and suffer from non-specific effects, arguing for a strong need to develop additional functional genomics approaches. RESULTS: In this work, we established a high-throughput lentivirus-based insertional mutagenesis strategy as a forward genetics screen tool in aneuploid cells. Application of this approach to human cell lines in multiple phenotypic screens suggested the presence of many yet uncharacterized functional elements in the human genome, represented at least in part by novel exons of known and novel genes. The novel transcripts containing these exons can be massively, up to thousands-fold, induced by specific stresses, and at least some can represent bi-cistronic protein-coding mRNAs. CONCLUSIONS: Altogether, these results argue that many unannotated and non-canonical human transcripts, including those that appear as aberrant splice products, have biological relevance under specific biological conditions.

Lessons from discovery of true ADAR RNA editing sites in a human cell line.

BACKGROUND: Conversion or editing of adenosine (A) into inosine (I) catalyzed by specialized cellular enzymes represents one of the most common post-transcriptional RNA modifications with emerging connection to disease. A-to-I conversions can happen at specific sites and lead to increase in proteome diversity and changes in RNA stability, splicing, and regulation. Such sites can be detected as adenine-to-guanine sequence changes by next-generation RNA sequencing which resulted in millions reported sites from multiple genome-wide surveys. Nonetheless, the lack of extensive independent validation in such endeavors, which is critical considering the relatively high error rate of next-generation sequencing, leads to lingering questions about the validity of the current compendiums of the editing sites and conclusions based on them. RESULTS: Strikingly, we found that the current analytical methods suffer from very high false positive rates and that a significant fraction of sites in the public databases cannot be validated. In this work, we present potential solutions to these problems and provide a comprehensive and extensively validated list of A-to-I editing sites in a human cancer cell line. Our findings demonstrate that most of true A-to-I editing sites in a human cancer cell line are located in the non-coding transcripts, the so-called RNA 'dark matter'. On the other hand, many ADAR editing events occurring in exons of human protein-coding mRNAs, including those that can recode the transcriptome, represent false positives and need to be interpreted with caution. Nonetheless, yet undiscovered authentic ADAR sites that increase the diversity of human proteome exist and warrant further identification. CONCLUSIONS: Accurate identification of human ADAR sites remains a challenging problem, particularly for the sites in exons of protein-coding mRNAs. As a result, genome-wide surveys of ADAR editome must still be accompanied by extensive Sanger validation efforts. However, given the vast number of unknown human ADAR sites, there is a need for further developments of the analytical techniques, potentially those that are based on deep learning solutions, in order to provide a quick and reliable identification of the editome in any sample.

Glucosinolate O-methyltransferase mediated callus formation and affected ROS homeostasis in Arabidopsis thaliana.

Auxin-induced callus formation was largely dependent on the function of Lateral Organ Boundaries Domain (LBD) family transcription factors. We previously revealed that two IGMT (Indole glucosinolate oxy-methyl transferase) genes, IGMT2 and IGMT3, may be involved in the callus formation process as potential target genes of LBD29. Overexpression of the IGMT genes induces spontaneous callus formation. However, the details of the IGMT involvement in callus formation process were not well studied. IGMT1-4, but not IGMT5, are targeted and induced by LBD29 during the early stage of callus formation. Cell membrane and nucleus localized IGMT3 was mainly expressed in the elongation and maturation zones tissues of the primary root and lateral root, which could be further accumulated after CIM treatment. The igmts quadruple mutant, which obtained by CRISPR/Cas9 technology, exhibits a phenotype of attenuated callus formation. Enhanced indole glucosinolate anabolic pathway caused by IGMT1-4 overexpression promotes callus formation. In addition, the IGMT genes were involved in the reactive oxygen species homeostasis, which could be responsible for its role on callus formation. This study provides novel insights into the role of IGMTs gene-mediated callus formation. Activation of the Indole glucosinolate anabolic pathway is an inducing factor for plant callus initiation. Supplementary Information: The online version contains supplementary material available at 10.1007/s12298-023-01409-2.

Hovlinc is a recently evolved class of ribozyme found in human lncRNA.

Although naturally occurring catalytic RNA molecules-ribozymes-have attracted a great deal of research interest, very few have been identified in humans. Here, we developed a genome-wide approach to discovering self-cleaving ribozymes and identified a naturally occurring ribozyme in humans. The secondary structure and biochemical properties of this ribozyme indicate that it belongs to an unidentified class of small, self-cleaving ribozymes. The sequence of the ribozyme exhibits a clear evolutionary path, from its appearance between ~130 and ~65 million years ago (Ma), to acquiring self-cleavage activity very recently, ~13-10 Ma, in the common ancestors of humans, chimpanzees and gorillas. The ribozyme appears to be functional in vivo and is embedded within a long noncoding RNA belonging to a class of very long intergenic noncoding RNAs. The presence of a catalytic RNA enzyme in lncRNA creates the possibility that these transcripts could function by carrying catalytic RNA domains.

Evidence for Existence of Multiple Functional Human Small RNAs Derived from Transcripts of Protein-Coding Genes.

The human genome encodes a multitude of different noncoding transcripts that have been traditionally separated on the basis of their lengths into long (>200 nt) or small (<200 nt) noncoding RNAs. The functions, mechanisms of action, and biological relevance of the vast majority of both long and short noncoding transcripts remain unknown. However, according to the functional understanding of the known classes of long and small noncoding RNAs (sncRNAs) that have been shown to play crucial roles in multiple biological processes, it is generally assumed that many unannotated long and small transcripts participate in important cellular functions as well. Nevertheless, direct evidence of functionality is lacking for most noncoding transcripts, especially for sncRNAs that are often dismissed as stable degradation products of longer RNAs. Here, we developed a high-throughput assay to test the functionality of sncRNAs by overexpressing them in human cells. Surprisingly, we found that a significant fraction (>40%) of unannotated sncRNAs appear to have biological relevance. Furthermore, contrary to the expectation, the potentially functional transcripts are not highly abundant and can be derived from protein-coding mRNAs. These results strongly suggest that the small noncoding transcriptome can harbor multiple functional transcripts that warrant future studies.

Genome-Wide Profiling of Endogenous Single-Stranded DNA Using the SSiNGLe-P1 Method.

Endogenous single-stranded DNA (essDNA) can form in a mammalian genome as the result of a variety of molecular processes and can both play important roles inside the cell as well as have detrimental consequences to genome integrity, much of which remains to be fully understood. Here, we established the SSiNGLe-P1 approach based on limited digestion by P1 endonuclease for high-throughput genome-wide identification of essDNA regions. We applied this method to profile essDNA in both human mitochondrial and nuclear genomes. In the mitochondrial genome, the profiles of essDNA provide new evidence to support the strand-displacement model of mitochondrial DNA replication. In the nuclear genome, essDNA regions were found to be enriched in certain types of functional genomic elements, particularly, the origins of DNA replication, R-loops, and to a lesser degree, in promoters. Furthermore, interestingly, many of the essDNA regions identified by SSiNGLe-P1 have not been annotated and thus could represent yet unknown functional elements.

Standing genetic variation as the predominant source for adaptation of a songbird.

What kind of genetic variation contributes the most to adaptation is a fundamental question in evolutionary biology. By resequencing genomes of 80 individuals, we inferred the origin of genomic variants associated with a complex adaptive syndrome involving multiple quantitative traits, namely, adaptation between high and low altitudes, in the vinous-throated parrotbill (Sinosuthora webbiana) in Taiwan. By comparing these variants with those in the Asian mainland population, we revealed standing variation in 24 noncoding genomic regions to be the predominant genetic source of adaptation. Parrotbills at both high and low altitudes exhibited signatures of recent selection, suggesting that not only the front but also the trailing edges of postglacial expanding populations could be subjected to environmental stresses. This study verifies and quantifies the importance of standing variation in adaptation in a cohort of genes, illustrating that the evolutionary potential of a population depends significantly on its preexisting genetic diversity. These findings provide important context for understanding adaptation and conservation of species in the Anthropocene.

Very long intergenic non-coding (vlinc) RNAs directly regulate multiple genes in cis and trans.

BACKGROUND: The majority of the human genome is transcribed in the form of long non-coding (lnc) RNAs. While these transcripts have attracted considerable interest, their molecular mechanisms of function and biological significance remain controversial. One of the main reasons behind this lies in the significant challenges posed by lncRNAs requiring the development of novel methods and concepts to unravel their functionality. Existing methods often lack cross-validation and independent confirmation by different methodologies and therefore leave significant ambiguity as to the authenticity of the outcomes. Nonetheless, despite all the caveats, it appears that lncRNAs may function, at least in part, by regulating other genes via chromatin interactions. Therefore, the function of a lncRNA could be inferred from the function of genes it regulates. In this work, we present a genome-wide functional annotation strategy for lncRNAs based on identification of their regulatory networks via the integration of three distinct types of approaches: co-expression analysis, mapping of lncRNA-chromatin interactions, and assaying molecular effects of lncRNA knockdowns obtained using an inducible and highly specific CRISPR/Cas13 system. RESULTS: We applied the strategy to annotate 407 very long intergenic non-coding (vlinc) RNAs belonging to a novel widespread subclass of lncRNAs. We show that vlincRNAs indeed appear to regulate multiple genes encoding proteins predominantly involved in RNA- and development-related functions, cell cycle, and cellular adhesion via a mechanism involving proximity between vlincRNAs and their targets in the nucleus. A typical vlincRNAs can be both a positive and negative regulator and regulate multiple genes both in trans and cis. Finally, we show vlincRNAs and their regulatory networks potentially represent novel components of DNA damage response and are functionally important for the ability of cancer cells to survive genotoxic stress. CONCLUSIONS: This study provides strong evidence for the regulatory role of the vlincRNA class of lncRNAs and a potentially important role played by these transcripts in the hidden layer of RNA-based regulation in complex biological systems.

Strategies to Annotate and Characterize Long Noncoding RNAs: Advantages and Pitfalls.

The past decade has seen an explosion of interest in long noncoding RNAs (lncRNAs). However, despite the massive volume of scientific data implicating these transcripts in a plethora of molecular and cellular processes, a great deal of controversy surrounds these RNAs. One of the main reasons for this lies in the multiple unique features of lncRNAs which limit the available methods used to characterize them. Combined with their vast numbers and inadequate classification, comprehensive annotation of these transcripts becomes a daunting task. The solution to this complex challenge likely lies in deep understanding of the strengths and weaknesses of each computational and empirical approach, and integration of multiple strategies to reduce noise, authenticate the results, and classify lncRNAs. We review here both the advantages and caveats of strategies commonly used for functional characterization and annotation of lncRNAs in the context of emerging conceptual guidelines for their application.

Graphene oxide exhibited positive effects on the growth of Aloe vera L.

There is increasing evidence for graphene associated plant growth promotion, however, the chronic effects of soil-applied graphene remain largely unexplored. The present study investigated the morphological, physiological and biochemical responses of graphene oxide (GO) on Aloe vera L. over the concentration range of 0-100 mg/L for four months. Our results demonstrated that GO, with the best efficiency at 50 mg/L, could enhance the photosynthetic capacity of leaves, increase the yield and morphological characters of root and leaf, improve the nutrient (protein and amino acid) contents of leaf, without reducing the content of the main bioactive compound aloin. Compared with leaves, the effect of GO on root growth was more obvious. Although the electrolyte leakage and MDA content were raised at high concentrations, GO treatment did not increase the root antioxidant enzymes activity or decrease the root vigor, which excluding typical stress response. Furthermore, injection experiments showed that the GO in vivo did not change the plant growth state obviously. Taken together, our study revealed the role of GO in promoting Aloe vera growth by stimulating root growth and photosynthesis, which would provide theory basis for GO application in agriculture and forestry. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s12298-021-00979-3.

Very-long-chain fatty acids restrict regeneration capacity by confining pericycle competence for callus formation in Arabidopsis.

The already differentiated organs in plants have a remarkable capacity to regenerate new individuals under culture conditions. Plant in vitro regeneration practically starts with the induction of a pluripotent cell mass, the callus, from detached organs on auxin-rich callus-inducing medium (CIM), which is generally required for subsequent regeneration of new bodies. Recent studies show that CIM-induced callus formation occurs from the pericycle or pericycle-like cells through a root developmental pathway, whereas the signals involved in governing callus-forming capacity of pericycle cells remain unknown. Here we report that very-long-chain fatty acids (VLCFAs) play a critical role in confining the pericycle competence for callus formation and thus the regeneration capacity of Arabidopsis By genetic screening, we identified the callus formation-related 1 (cfr1) mutant, which bypasses the inhibition of callus-forming capacity in roots by solitary-root (slr/iaa14). We show that CFR1 encodes 3-ketoacyl-CoA synthase 1 (KCS1), which catalyzes a rate-limiting step of VLCFA biosynthesis. Our biochemical and genetic analyses demonstrate that VLCFAs restrict the pericycle competence for callus formation, at least in part, by regulating the transcription of Aberrant Lateral Root Formation 4 (ALF4). Moreover, we provide evidence that VLCFAs act as cell layer signals to mediate the pericycle competence for callus formation. Taken together, our results identify VLCFAs or their derivatives as the confining signals for mediating the pericycle competence for callus formation and thus the regeneration capacity of plant organs.

Genome-Wide Identification of Arabidopsis LBD29 Target Genes Reveals the Molecular Events behind Auxin-Induced Cell Reprogramming during Callus Formation.

Auxin-induced callus formation represents an important cell reprogramming process during in vitro regeneration of plants, in which the pericycle or pericycle-like cells within plant organs are reprogrammed into the pluripotent cell mass termed callus that is generally required for subsequent regeneration of root or shoot. However, the molecular events behind cell reprogramming during auxin-induced callus formation are largely elusive. We previously identified that auxin-induced LATERAL ORGAN BOUNDARIES DOMAIN (LBD) transcription factors act as the master regulators to trigger auxin-induced callus formation. Here, by ChIP-seq (chromatin immunoprecipitation-based sequencing) and RNA sequencing approaches, we identified the potential LBD29 target genes at the genome-wide level and outlined the molecular events of LBD-triggered cell reprogramming during callus formation. We showed that LBD29 preferentially bound to the G-box (CACGTG) and TGGGC[C/T] motifs and potentially targeted >350 genes, among which the genes related to methylation, reactive oxygen species (ROS) metabolism, cell wall hydrolysis and lipid metabolism were rapidly activated, while most of the light-responsive genes were suppressed by LBD29. Further examination of a few representative genes validated that they were targeted by LBD29 and participated in the regulation of cell reprogramming during callus formation. Our data not only outline a framework of the early molecular events behind auxin-induced cell reprogramming of callus formation, but also provide a valuable resource for identification of genes that regulate cell fate switch during in vitro regeneration of plants.

Single nucleus genome sequencing reveals high similarity among nuclei of an endomycorrhizal fungus.

Nuclei of arbuscular endomycorrhizal fungi have been described as highly diverse due to their asexual nature and absence of a single cell stage with only one nucleus. This has raised fundamental questions concerning speciation, selection and transmission of the genetic make-up to next generations. Although this concept has become textbook knowledge, it is only based on studying a few loci, including 45S rDNA. To provide a more comprehensive insight into the genetic makeup of arbuscular endomycorrhizal fungi, we applied de novo genome sequencing of individual nuclei of Rhizophagus irregularis. This revealed a surprisingly low level of polymorphism between nuclei. In contrast, within a nucleus, the 45S rDNA repeat unit turned out to be highly diverged. This finding demystifies a long-lasting hypothesis on the complex genetic makeup of arbuscular endomycorrhizal fungi. Subsequent genome assembly resulted in the first draft reference genome sequence of an arbuscular endomycorrhizal fungus. Its length is 141 Mbps, representing over 27,000 protein-coding gene models. We used the genomic sequence to reinvestigate the phylogenetic relationships of Rhizophagus irregularis with other fungal phyla. This unambiguously demonstrated that Glomeromycota are more closely related to Mucoromycotina than to its postulated sister Dikarya.

Positive selection and functional divergence of R2R3-MYB paralogous genes expressed in inflorescence buds of Scutellaria species (Labiatae).

Anthocyanin is the main pigment forming floral diversity. Several transcription factors that regulate the expression of anthocyanin biosynthetic genes belong to the R2R3-MYB family. Here we examined the transcriptomes of inflorescence buds of Scutellaria species (skullcaps), identified the expression R2R3-MYBs, and detected the genetic signatures of positive selection for adaptive divergence across the rapidly evolving skullcaps. In the inflorescence buds, seven R2R3-MYBs were identified. MYB11 and MYB16 were detected to be positively selected. The signature of positive selection on MYB genes indicated that species diversification could be affected by transcriptional regulation, rather than at the translational level. When comparing among the background lineages of Arabidopsis, tomato, rice, and Amborella, heterogeneous evolutionary rates were detected among MYB paralogs, especially between MYB13 and MYB19. Significantly different evolutionary rates were also evidenced by type-I functional divergence between MYB13 and MYB19, and the accelerated evolutionary rates in MYB19, implied the acquisition of novel functions. Another paralogous pair, MYB2/7 and MYB11, revealed significant radical amino acid changes, indicating divergence in the regulation of different anthocyanin-biosynthetic enzymes. Our findings not only showed that Scutellaria R2R3-MYBs are functionally divergent and positively selected, but also indicated the adaptive relevance of regulatory genes in floral diversification.

Effects of Graphene-Based Nanomaterials on Microorganisms and Soil Microbial Communities.

The past decades have witnessed intensive research on the biological effects of graphene-based nanomaterials (GBNs) and the application of GBNs in different fields. The published literature shows that GBNs exhibit inhibitory effects on almost all microorganisms under pure culture conditions, and that this inhibitory effect is influenced by the microbial species, the GBN's physicochemical properties, the GBN's concentration, treatment time, and experimental surroundings. In addition, microorganisms exist in the soil in the form of microbial communities. Considering the complex interactions between different soil components, different microbial communities, and GBNs in the soil environment, the effects of GBNs on soil microbial communities are undoubtedly intertwined. Since bacteria and fungi are major players in terrestrial biogeochemistry, this review focuses on the antibacterial and antifungal performance of GBNs, their antimicrobial mechanisms and influencing factors, as well as the impact of this effect on soil microbial communities. This review will provide a better understanding of the effects of GBNs on microorganisms at both the individual and population scales, thus providing an ecologically safe reference for the release of GBNs to different soil environments.

Genome-wide profiles of DNA damage represent highly accurate predictors of mammalian age.

The identification of novel age-related biomarkers represents an area of intense research interest. Despite multiple studies associating DNA damage with aging, there is a glaring paucity of DNA damage-based biomarkers of age, mainly due to the lack of precise methods for genome-wide surveys of different types of DNA damage. Recently, we developed two techniques for genome-wide mapping of the most prevalent types of DNA damage, single-strand breaks and abasic sites, with nucleotide-level resolution. Herein, we explored the potential of genomic patterns of DNA damage identified by these methods as a source of novel age-related biomarkers using mice as a model system. Strikingly, we found that models based on genomic patterns of either DNA lesion could accurately predict age with higher precision than the commonly used transcriptome analysis. Interestingly, the informative patterns were limited to relatively few genes and the DNA damage levels were positively or negatively correlated with age. These findings show that previously unexplored high-resolution genomic patterns of DNA damage contain useful information that can contribute significantly to both practical applications and basic science.

Methods to Analyze the Non-Coding RNA Interactome-Recent Advances and Challenges.

Most of the human genome is transcribed to generate a multitude of non-coding RNAs. However, while these transcripts have generated an immense amount of scientific interest, their biological function remains a subject of an intense debate. Understanding mechanisms of action of non-coding RNAs is a key to addressing the issue of biological relevance of these transcripts. Based on some well-understood non-coding RNAs that function inside the cell by interacting with other molecules, it is generally believed many other non-coding transcripts could also function in a similar fashion. Therefore, development of methods that can map RNA interactome is the key to understanding functionality of the extensive cellular non-coding transcriptome. Here, we review the vast progress that has been made in the past decade in technologies that can map RNA interactions with different sites in DNA, proteins or other RNA molecules; the general approaches used to validate the existence of novel interactions; and the challenges posed by interpreting the data obtained using the interactome mapping methods.

Complex genomic patterns of abasic sites in mammalian DNA revealed by a high-resolution SSiNGLe-AP method.

DNA damage plays a critical role in biology and diseases; however, how different types of DNA lesions affect cellular functions is far from clear mostly due to the paucity of high-resolution methods that can map their locations in complex genomes, such as those of mammals. Here, we present the development and validation of SSiNGLe-AP method, which can map a common type of DNA damage, abasic (AP) sites, in a genome-wide and high-resolution manner. We apply this method to six different tissues of mice with different ages and human cancer cell lines. We find a nonrandom distribution of AP sites in the mammalian genome that exhibits dynamic enrichment at specific genomic locations, including single-nucleotide hotspots, and is significantly influenced by gene expression, age and tissue type in particular. Overall, these results suggest that we are only starting to understand the true complexities in the genomic patterns of DNA damage.