barCoder: a tool to generate unique, orthogonal genetic tags for qPCR detection

BACKGROUND

Tracking dispersal of microbial populations in the environment requires specific detection methods that discriminate between the target strain and all potential natural and artificial interferents, including previously utilized tester strains. Recent work has shown that genomic insertion of short identification tags, called "barcodes" here, allows detection of chromosomally tagged strains by real-time PCR. Manual design of these barcodes is feasible for small sets, but expansion of the technique to larger pools of distinct and well-functioning assays would be significantly aided by software-guided design.

RESULTS

Here we introduce barCoder, a bioinformatics tool that facilitates the process of creating sets of uniquely identifiable barcoded strains. barCoder utilizes the genomic sequence of the target strain and a set of user-specified PCR parameters to generate a list of suggested barcode "modules" that consist of binding sites for primers and probes, and appropriate spacer sequences. Each module is designed to yield optimal PCR amplification and unique identification. Optimal amplification includes metrics such as ideal melting temperature and G+C content, appropriate spacing, and minimal stem-loop formation; unique identification includes low BLAST hits against the target organism, previously generated barcode modules, and databases (such as NCBI). We tested the ability of our algorithm to suggest appropriate barcodes by generating 12 modules for Bacillus thuringiensis serovar kurstaki-a simulant for the potential biowarfare agent Bacillus anthracis-and three each for other potential target organisms with variable G+C content. Real-time PCR detection assays directed at barcodes were specific and yielded minimal cross-reactivity with a panel of near-neighbor and potential contaminant materials.

CONCLUSIONS

The barCoder algorithm facilitates the generation of synthetically barcoded biological simulants by (a) eliminating the task of creating modules by hand, (b) minimizing optimization of PCR assays, and (c) reducing effort wasted on non-unique barcode modules.

Marine health of the Arabian Gulf: Drivers of pollution and assessment approaches focusing on desalination activities

The Arabian Gulf is one of the most adversely affected marine environments worldwide, which results from combined pollution drivers including climate change, oil and gas activities, and coastal anthropogenic disturbances. Desalination activities are one of the major marine pollution drivers regionally and internationally. Arabian Gulf countries represent a hotspot of desalination activities as they are responsible for nearly 50% of the global desalination capacity. Building desalination plants, up-taking seawater, and discharging untreated brine back into the sea adversely affects the biodiversity of the marine ecosystems. The present review attempted to reveal the potential negative effects of desalination plants on the Gulf's marine environments. We emphasised different conventional and innovative assessment tools used to assess the health of marine environments and evaluate the damage exerted by desalination activity in the Gulf. Finally, we suggested effective management approaches to tackle the issue including the significance of national regulations and regional cooperation.

From Genome-Wide to Environment-Wide: Capturing the Environome

Genome-wide association (GWA) studies have shown that genetic influences on individual differences in affect, behavior, and cognition are driven by thousands of DNA variants, each with very small effect sizes. Here, we propose taking inspiration from GWA studies for understanding and modeling the influence of the environment on complex phenotypes. We argue that the availability of DNA microarrays in genetic research is comparable with the advent of digital technologies in psychological science that enable collecting rich, naturalistic observations in real time of the environome, akin to the genome. These data can capture many thousand environmental elements, which we speculate each influence individual differences in affect, behavior, and cognition with very small effect sizes, akin to findings from GWA studies about DNA variants. We outline how the principles and mechanisms of genetic influences on psychological traits can be applied to improve the understanding and models of the environome.

A chromosome-scale assembly of allotetraploid Brassica juncea (AABB) elucidates comparative architecture of the A and B genomes

Brassica juncea (AABB), commonly referred to as mustard, is a natural allopolyploid of two diploid species-B. rapa (AA) and B. nigra (BB). We report a highly contiguous genome assembly of an oleiferous type of B. juncea variety Varuna, an archetypical Indian gene pool line of mustard, with ~100× PacBio single-molecule real-time (SMRT) long reads providing contigs with an N50 value of >5 Mb. Contigs were corrected for the misassemblies and scaffolded with BioNano optical mapping. We also assembled a draft genome of B. nigra (BB) variety Sangam using Illumina short-read sequencing and Oxford Nanopore long reads and used it to validate the assembly of the B genome of B. juncea. Two different linkage maps of B. juncea, containing a large number of genotyping-by-sequencing markers, were developed and used to anchor scaffolds/contigs to the 18 linkage groups of the species. The resulting chromosome-scale assembly of B. juncea Varuna is a significant improvement over the previous draft assembly of B. juncea Tumida, a vegetable type of mustard. The assembled genome was characterized for transposons, centromeric repeats, gene content and gene block associations. In comparison to the A genome, the B genome contains a significantly higher content of LTR/Gypsy retrotransposons, distinct centromeric repeats and a large number of B. nigra specific gene clusters that break the gene collinearity between the A and the B genomes. The B. juncea Varuna assembly will be of major value to the breeding work on oleiferous types of mustard that are grown extensively in south Asia and elsewhere.

Comparative analysis of long noncoding RNAs in angiosperms and characterization of long noncoding RNAs in response to heat stress in Chinese cabbage

Long noncoding RNAs (lncRNAs) are widely present in different species and play critical roles in response to abiotic stresses. However, the functions of lncRNAs in Chinese cabbage under heat stress remain unknown. Here, we first conducted a global comparative analysis of 247,242 lncRNAs among 37 species. The results indicated that lncRNAs were poorly conserved among different species, and only 960 lncRNAs were homologous to 524 miRNA precursors. We then carried out lncRNA sequencing for a genome-wide analysis of lncRNAs and their target genes in Chinese cabbage at different stages of heat treatment. In total, 18,253 lncRNAs were identified, of which 1229 differentially expressed (DE) lncRNAs were characterized as being heat-responsive. The ceRNA network revealed that 38 lncRNAs, 16 miRNAs, and 167 mRNAs were involved in the heat response in Chinese cabbage. Combined analysis of the cis- and trans-regulated genes indicated that the targets of DE lncRNAs were significantly enriched in the "protein processing in endoplasmic reticulum" and "plant hormone signal transduction" pathways. Furthermore, the majority of HSP and PYL genes involved in these two pathways exhibited similar expression patterns and responded to heat stress rapidly. Based on the networks of DE lncRNA-mRNAs, 29 and 22 lncRNAs were found to interact with HSP and PYL genes, respectively. Finally, the expression of several critical lncRNAs and their targets was verified by qRT-PCR. Overall, we conducted a comparative analysis of lncRNAs among 37 species and performed a comprehensive analysis of lncRNAs in Chinese cabbage. Our findings expand the knowledge of lncRNAs involved in the heat stress response in Chinese cabbage, and the identified lncRNAs provide an abundance of resources for future comparative and functional studies.

Centromeres are dismantled by foundational meiotic proteins Spo11 and Rec8

Meiotic processes are potentially dangerous to genome stability and could be disastrous if activated in proliferative cells. Here we show that two key meiosis-defining proteins, the topoisomerase Spo11 (which forms double-strand breaks) and the meiotic cohesin Rec8, can dismantle centromeres. This dismantlement is normally observable only in mutant cells that lack the telomere bouquet, which provides a nuclear microdomain conducive to centromere reassembly¹; however, overexpression of Spo11 or Rec8 leads to levels of centromere dismantlement that cannot be countered by the bouquet. Specific nucleosome remodelling factors mediate centromere dismantlement by Spo11 and Rec8. Ectopic expression of either protein in proliferating cells leads to the loss of mitotic kinetochores in both fission yeast and human cells. Hence, while centromeric chromatin has been characterized as extraordinarily stable, Spo11 and Rec8 challenge this stability and may jeopardize kinetochores in cancers that express meiotic proteins.

Study of Zuojin Pill in Treating Chronic Atrophic Gastritis by UPLC-Q-TOF/MS Based on Serum and Urine Metabolomics Combined with Network Pharmacology

Zuojin Pill (ZJP) is widely used for the treatment of gastrointestinal diseases, while its specific mechanism has not been systematically investigated. The aim of this study was to explore the therapeutic effects and potential mechanism of ZJP in chronic atrophic gastritis (CAG) through UPLC-Q-TOF/MS-based metabolomics combined with network pharmacology. ZJP and omeprazole significantly reduce contents of IL-1β, IL-6, IL-10, and iNOS and improve pathological characteristics. Metabolomic results indicated that the therapeutic effects of ZJP were mainly related to ten metabolites, namely, choline, L-threonine, hydroxypyruvic acid, creatine, taurine, succinic acid, cis-aconitic acid, citric acid, succinic acid semialdehyde, and uric acid. Pathway analysis showed that the treatment of CAG by ZJP was associated with taurine and hypotaurine metabolism; glyoxylate and dicarboxylate metabolism; glycine, serine, and threonine metabolism; glycerophospholipid metabolism; citrate cycle (TCA cycle), alanine, aspartate, and glutamate metabolism; butanoate metabolism; and purine metabolism. Validation of metabolic markers and key targets of network pharmacology through RT-PCR analysis showed that ZJP significantly downregulated a series of inflammatory markers, such as MAPK1, PKIA, RB1, SCN5A, RXRA, E2F1, PTGS1, IGF2, ADRB1, ADRA1B, PTGS2, and GABRA1. This study was the first to use a combination of metabolomics and network pharmacology to clarify the therapeutic effects of ZJP on CAG and the regulation of multiple metabolic pathways.

Circulating microRNA-15b-5p as a biomarker for asthma-COPD overlap

BACKGROUND

It remains unclear how to characterize different subtypes of asthma and chronic obstructive pulmonary disease (COPD). We previously described serum periostin and chitinase-3-like protein 1 (YKL-40) as useful markers for asthma-COPD overlap (ACO). MicroRNAs (miRNAs) are now recognized as markers for identifying the pathophysiological features in several diseases. This study aimed to identify circulating miRNAs that could discriminate patients with ACO from patients with asthma or COPD.

METHODS

This study included two independent cohorts. First, we screened 84 miRNAs for expression levels in patients with ACO (n = 6) or asthma (n = 6) using a quantitative real-time PCR array. The miRNAs showing at least a 2-fold difference in the discovery phase were analyzed in 30 patients each with asthma, COPD, or ACO in the replication phase. The diagnostic accuracy was evaluated using the area under the receiver operating characteristic curve (AUROC).

RESULTS

Nine miRNAs were identified in the discovery phase. Five of these miRNAs (miR-148a-3p, miR-15b-5p, miR-223-3p, miR-23a-3p, and miR-26b-5p) had lower levels in ACO patients and could discriminate between ACO patients and patients with either asthma or COPD. miR-15b-5p was the most accurate miRNA for the discrimination of patients with ACO (AUROC, 0.71). Moreover, the combined assessment of miR-15b-5p, serum periostin, and YKL-40 (AUROC, 0.80) improved diagnostic accuracy for ACO compared with the combined model of periostin and YKL-40 (AUROC, 0.69).

CONCLUSIONS

Circulating miR-15b-5p is a potential marker for identifying patients with ACO. By elucidating the molecular pathways controlled by miRNAs, we may better understand the pathophysiology of ACO.

Male gametophyte development in flowering plants: A story of quarantine and sacrifice

Over 160 years ago, scientists made the first microscopic observations of angiosperm pollen. Unlike in animals, male meiosis in angiosperms produces a haploid microspore that undergoes one asymmetric division to form a vegetative cell and a generative cell. These two cells have distinct fates: the vegetative cell exits the cell cycle and elongates to form a tip-growing pollen tube; the generative cell divides once more in the pollen grain or within the growing pollen tube to form a pair of sperm cells. The concept that male germ cells are less active than the vegetative cell came from biochemical analyses and pollen structure anatomy early in the last century and is supported by the pollen transcriptome data of the last decade. However, the mechanism of how and when the transcriptional repression in male germ cells occurs is still not fully understood. In this review, we provide a brief account of the cytological and metabolic differentiation between the vegetative cell and male germ cells, with emphasis on the role of temporary callose walls, dynamic nuclear pore density, transcription repression, and histone variants. We further discuss the intercellular movement of small interfering RNA (siRNA) derived from transposable elements (TEs) and reexamine the function of TE expression in male germ cells.

Developmental genetics of maize vegetative shoot architecture

More than 1.1 billion tonnes of maize grain were harvested across 197 million hectares in 2019 (FAOSTAT 2020). The vast global productivity of maize is largely driven by denser planting practices, higher yield potential per area of land, and increased yield potential per plant. Shoot architecture, the three-dimensional structural arrangement of the above-ground plant body, is critical to maize grain yield and biomass. Structure of the shoot is integral to all aspects of modern agronomic practices. Here, the developmental genetics of the maize vegetative shoot is reviewed. Plant architecture is ultimately determined by meristem activity, developmental patterning, and growth. The following topics are discussed: shoot apical meristem, leaf architecture, axillary meristem and shoot branching, and intercalary meristem and stem activity. Where possible, classical and current studies in maize developmental genetics, as well as recent advances leveraged by "-omics" analyses, are highlighted within these sections.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01208-1.

Identification of HIF-dependent alternative splicing in gastrointestinal cancers and characterization of a long, coding isoform of SLC35A3

Intra-tumor hypoxia is a common feature in many solid cancers. Although transcriptional targets of hypoxia-inducible factors (HIFs) have been well characterized, alternative splicing or processing of pre-mRNA transcripts which occurs during hypoxia and subsequent HIF stabilization is much less understood. Here, we identify many HIF-dependent alternative splicing events after whole transcriptome sequencing in pancreatic cancer cells exposed to hypoxia with and without downregulation of the aryl hydrocarbon receptor nuclear translocator (ARNT), a protein required for HIFs to form a transcriptionally active dimer. We correlate the discovered hypoxia-driven events with available sequencing data from pan-cancer TCGA patient cohorts to select a narrow set of putative biologically relevant splice events for experimental validation. We validate a small set of candidate HIF-dependent alternative splicing events in multiple human gastrointestinal cancer cell lines as well as patient-derived human pancreatic cancer organoids. Lastly, we report the discovery of a HIF-dependent mechanism to produce a hypoxia-dependent, long and coding isoform of the UDP-N-acetylglucosamine transporter SLC35A3.

Mi3-GPU: MCMC-based Inverse Ising Inference on GPUs for protein covariation analysis

Inverse Ising inference is a method for inferring the coupling parameters of a Potts/Ising model based on observed site-covariation, which has found important applications in protein physics for detecting interactions between residues in protein families. We introduce Mi3-GPU ("mee-three", for MCMC Inverse Ising Inference) software for solving the inverse Ising problem for protein-sequence datasets with few analytic approximations, by parallel Markov-Chain Monte-Carlo sampling on GPUs. We also provide tools for analysis and preparation of protein-family Multiple Sequence Alignments (MSAs) to account for finite-sampling issues, which are a major source of error or bias in inverse Ising inference. Our method is "generative" in the sense that the inferred model can be used to generate synthetic MSAs whose mutational statistics (marginals) can be verified to match the dataset MSA statistics up to the limits imposed by the effects of finite sampling. Our GPU implementation enables the construction of models which reproduce the covariation patterns of the observed MSA with a precision that is not possible with more approximate methods. The main components of our method are a GPU-optimized algorithm to greatly accelerate MCMC sampling, combined with a multi-step Quasi-Newton parameter-update scheme using a "Zwanzig reweighting" technique. We demonstrate the ability of this software to produce generative models on typical protein family datasets for sequence lengths L ~ 300 with 21 residue types with tens of millions of inferred parameters in short running times.

The emerging role of chromatin remodelers in neurodevelopmental disorders: a developmental perspective

Neurodevelopmental disorders (NDDs), including intellectual disability (ID) and autism spectrum disorders (ASD), are a large group of disorders in which early insults during brain development result in a wide and heterogeneous spectrum of clinical diagnoses. Mutations in genes coding for chromatin remodelers are overrepresented in NDD cohorts, pointing towards epigenetics as a convergent pathogenic pathway between these disorders. In this review we detail the role of NDD-associated chromatin remodelers during the developmental continuum of progenitor expansion, differentiation, cell-type specification, migration and maturation. We discuss how defects in chromatin remodelling during these early developmental time points compound over time and result in impaired brain circuit establishment. In particular, we focus on their role in the three largest cell populations: glutamatergic neurons, GABAergic neurons, and glia cells. An in-depth understanding of the spatiotemporal role of chromatin remodelers during neurodevelopment can contribute to the identification of molecular targets for treatment strategies.

Investigating the relationship between the SNCA gene and cognitive abilities in idiopathic Parkinson's disease using machine learning

Cognitive impairments are prevalent in Parkinson's disease (PD), but the underlying mechanisms of their development are unknown. In this study, we aimed to predict global cognition (GC) in PD with machine learning (ML) using structural neuroimaging, genetics and clinical and demographic characteristics. As a post-hoc analysis, we aimed to explore the connection between novel selected features and GC more precisely and to investigate whether this relationship is specific to GC or is driven by specific cognitive domains. 101 idiopathic PD patients had a cognitive assessment, structural MRI and blood draw. ML was performed on 102 input features including demographics, cortical thickness and subcortical measures, and several genetic variants (APOE, MAPT, SNCA, etc.). Using the combination of RRELIEFF and Support Vector Regression, 11 features were found to be predictive of GC including sex, rs894280, Edinburgh Handedness Inventory, UPDRS-III, education, five cortical thickness measures (R-parahippocampal, L-entorhinal, R-rostral anterior cingulate, L-middle temporal, and R-transverse temporal), and R-caudate volume. The rs894280 of SNCA gene was selected as the most novel finding of ML. Post-hoc analysis revealed a robust association between rs894280 and GC, attention, and visuospatial abilities. This variant indicates a potential role for the SNCA gene in cognitive impairments of idiopathic PD.

Genome-wide profiling of miRNAs during seed development reveals their functional relevance in seed size/weight determination in chickpea

MicroRNAs (miRNAs) are non-coding small RNAs that regulate gene expression at transcriptional and post-transcriptional levels. The role of miRNAs in seed development and seed size/weight determination is poorly understood in legumes. In this study, we profiled miRNAs at seven successive stages of seed development in a small-seeded and a large-seeded chickpea cultivar via small RNA sequencing. In total, 113 known and 243 novel miRNAs were identified. Gene ontology analysis revealed the enrichment of seed/reproductive/post-embryonic development and signaling pathways processes among the miRNA target genes. A large fraction of the target genes exhibited antagonistic correlation with miRNA expression. The sets of co-expressed miRNAs showing differential expression between the two cultivars were recognized. Known transcription factor (TF) encoding genes involved in seed size/weight determination, including SPL, GRF, MYB, ARF, HAIKU1, SHB1, KLUH/CYP78A5, and E2Fb along with novel genes were found to be targeted by the predicted miRNAs. Differential expression analysis revealed higher transcript levels of members of SPL and REVOLUTA TF families and lower expression of their corresponding miRNAs in the large-seeded cultivar. At least 19 miRNAs known to be involved in seed development or differentially expressed between small-seeded and large-seeded cultivars at late-embryogenesis and/or mid-maturation stages were located within known quantitative trait loci (QTLs) associated with seed size/weight determination. Moreover, 41 target genes of these miRNAs were also located within these QTLs. Altogether, we revealed important roles of miRNAs in seed development and identified candidate miRNAs and their target genes that have functional relevance in determining seed size/weight in chickpea.

Antibiotic resistance: Global health crisis and metagenomics

Antibiotic resistance is a global problem which affects human health. The imprudent use of antibiotics (medicine, agriculture, aquaculture, and food industry) has resulted in the broader dissemination of resistance. Urban wastewater & sewage treatment plants act as the hotspot for the widespread of antimicrobial resistance. Natural environment also plays an important role in the dissemination of resistance. Mapping of antibiotic resistance genes (ARGS) in environment is essential for mitigating antimicrobial resistance (AMR) widespread. Therefore, the review article emphasizes on the application of metagenomics for the surveillance of antimicrobial resistance. Metagenomics is the next generation tool which is being used for cataloging the resistome of diverse environments. We summarize the different metagenomic tools that can be used for mining of ARGs and acquired AMR present in the metagenomic data. Also, we recommend application of targeted sequencing/ capture platform for mapping of resistome with higher specificity and selectivity.

Androgen receptor decreases renal cell carcinoma bone metastases via suppressing the osteolytic formation through altering a novel circEXOC7 regulatory axis

BACKGROUND

Clear cell renal cell carcinoma (ccRCC) has gender differences, with the androgen receptor (AR) linked positively with metastasis to the lung. Its linkage to ccRCC bone metastases (RBMs), however, remains unclear.

METHODS

In the current study, five human RCC and five RCC bone metastasis tissues were deeply sequenced using Arraystar human circRNA V2.0 microarray. We conducted gain-of-function screening in vitro and in vivo to elucidate the AR's role in the RBM. Loss/gain-of-function was also implemented to verify the roles of related non-coding RNAs and proteins.

RESULTS

We uncovered that RBM also has a gender difference showing higher AR expression may be linked to fewer RBMs, which might involve suppressing osteolytic formation. Mechanism dissection indicates that AR can decrease the circular RNA EXOC7 (circEXOC7), expression via enhancing transcription of DHX9, a regulatory protein in circRNA biogenesis. The circEXOC7 can sponge/suppress miR-149-3p resulting in suppressing the CSF1 expression by directly binding to the 3'UTR region of CSF1 mRNA. Results from clinical epidemiological surveys also found that AR has a positive correlation with miR-149-3p and a negative correlation with CSF1 in AR-positive ccRCC tissues. Preclinical studies with Balb/c nude mouse model also validated that targeting this newly verified AR/DHX9/circEXOC7/miR-149-3p/CSF1 signaling via altering circEXOC7 or AR could lead to suppressing the RBM progression.

CONCLUSIONS

These data showed that AR/DHX9/circEXOC7/miR-149-3p/CSF1 signaling acts as a valuable feature in the bone metastasis of renal cancer, which may benefit in suppressing the RBM progression.

Novel circular RNA circSOBP governs amoeboid migration through the regulation of the miR-141-3p/MYPT1/p-MLC2 axis in prostate cancer

BACKGROUND

Metastatic prostate cancer is a fatal disease despite multiple new approvals in recent years. Recent studies revealed that circular RNAs (circRNAs) can be involved in cancer metastasis. Defining the role of circRNAs in prostate cancer metastasis and discovering therapeutic targets that block cancer metastasis is of great significance for the treatment of prostate cancer.

METHODS

The circSOBP levels in prostate cancer (PCa) were determined by qRT-PCR. We evaluated the function of circSOBP using a transwell assay and nude mice lung metastasis models. Immunofluorescence assay and electron microscopic assay were applied to determine the phenotypes of prostate cancer cells' migration. We used fluorescence in situ hybridization assay to determine the localization of RNAs. Dual luciferase and rescue assays were applied to verify the interactions between circSOBP, miR-141-3p, MYPT1, and phosphomyosin light chain (p-MLC2).

RESULTS

We observed that circSOBP level was significantly lower in PCa specimens compared with adjacent noncancerous prostate specimens, and was correlated with the grade group of PCa. Overexpression of circSOBP suppressed PCa migration and invasion in vitro and metastasis in vivo. CircSOBP depletion increased migration and invasion and induced amoeboid migration of PCa cells. Mechanistically, circSOBP bound miR-141-3p and regulated the MYPT1/p-MLC2 axis. Moreover, the depletion of MYPT1 reversed the inhibitory effect of circSOBP on the migration and invasion of PCa cells. Complementary intronic Alu elements induced but were not necessary for the formation of circSOBP. The nuclear export of circSOBP was mediated by URH49.

CONCLUSION

Our results suggest that circSOBP suppresses amoeboid migration of PCa cells and inhibits migration and invasion through sponging miR-141-3p and regulating the MYPT1/p-MLC2 axis.

Potential of long non-coding RNAs as a therapeutic target and molecular markers in glioblastoma pathogenesis

Glioblastoma (GB) is by far the most hostile type of malignant tumor that primarily affects the brain and spine, derived from star-shaped glial cells that are astrocytes and oligodendrocytes. Despite of significant efforts in recent years in glioblastoma research, the clinical efficacy of existing medical intervention is still limited and very few potential diagnostic markers are available. Long non-coding RNAs (lncRNAs) that lacks protein-coding capabilities were previously thought to be "junk sequences" in mammalian genomes are quite indispensible epigenetic regulators that can positively or negatively regulate gene expression and nuclear architecture, with significant roles in the initiation and development of tumors. Nevertheless, the precise mechanism of these distortedly expressed lncRNAs in glioblastoma pathogenesis is not yet fully understood. Since the advent of high-throughput sequencing technologies, more and more research have elucidated that lncRNAs are one of the most promising prognostic biomarkers and therapeutic targets for glioblastoma. In this paper, I briefly outlined the existing findings of lncRNAs. And also summarizes the profiles of different lncRNAs that have been broadly classified in glioblastoma research, with emphasis on both their prognostic and therapeutic values.

BRCA1 and RNAi factors promote repair mediated by small RNAs and PALB2-RAD52

Strong connections exist between R-loops (three-stranded structures harbouring an RNA:DNA hybrid and a displaced single-strand DNA), genome instability and human disease^1-5. Indeed, R-loops are favoured in relevant genomic regions as regulators of certain physiological processes through which homeostasis is typically maintained. For example, transcription termination pause sites regulated by R-loops can induce the synthesis of antisense transcripts that enable the formation of local, RNA interference (RNAi)-driven heterochromation⁶. Pause sites are also protected against endogenous single-stranded DNA breaks by BRCA1⁷. Hypotheses about how DNA repair is enacted at pause sites include a role for RNA, which is emerging as a normal, albeit unexplained, regulator of genome integrity⁸. Here we report that a species of single-stranded, DNA-damage-associated small RNA (sdRNA) is generated by a BRCA1-RNAi protein complex. sdRNAs promote DNA repair driven by the PALB2-RAD52 complex at transcriptional termination pause sites that form R-loops and are rich in single-stranded DNA breaks. sdRNA repair operates in both quiescent (G0) and proliferating cells. Thus, sdRNA repair can occur in intact tissue and/or stem cells, and may contribute to tumour suppression mediated by BRCA1.

Complex Linkage Disequilibrium Effects in HLA-DPB1 Expression and Molecular Mismatch Analyses of Transplantation Outcomes

BACKGROUND

HLA molecular mismatch (MM) is a risk factor for de novo donor-specific antibody (dnDSA) development in solid organ transplantation. HLA expression differences have also been associated with adverse outcomes in hematopoietic cell transplantation. We sought to study both MM and expression in assessing dnDSA risk.

METHODS

One hundred three HLA-DP-mismatched solid organ transplantation pairs were retrospectively analyzed. MM was computed using amino acids (aa), eplets, and, supplementarily, Grantham/Epstein scores. DPB1 alleles were classified as rs9277534-A (low-expression) or rs9277534-G (high-expression) linked. To determine the associations between risk factors and dnDSA, logistic regression, linkage disequilibrium (LD), and population-based analyses were performed.

RESULTS

A high-risk AA:GX (recipient:donor) expression combination (X = A or G) demonstrated strong association with HLA-DP dnDSA (P = 0.001). MM was also associated with HLA-DP dnDSA when evaluated by itself (eplet P = 0.007, aa P = 0.003, Grantham P = 0.005, Epstein P = 0.004). When attempting to determine the relative individual effects of the risk factors in multivariable analysis, only AA:GX expression status retained a strong association (relative risk = 18.6, P = 0.007 with eplet; relative risk = 15.8, P = 0.02 with aa), while MM was no longer significant (eplet P = 0.56, aa P = 0.51). Importantly, these risk factors are correlated, due to LD between the expression-tagging single-nucleotide polymorphism and polymorphisms along HLA-DPB1.

CONCLUSIONS

The MM and expression risk factors each appear to be strong predictors of HLA-DP dnDSA and to possess clinical utility; however, these two risk factors are closely correlated. These metrics may represent distinct ways of characterizing a common overlapping dnDSA risk profile, but they are not independent. Further, we demonstrate the importance and detailed implications of LD effects in dnDSA risk assessment and possibly transplantation overall.

A tumor-suppressive circular RNA mediates uncanonical integrin degradation by the proteasome in liver cancer

Circular RNAs (circRNAs) have emerged as important regulators of various cellular processes and have been implicated in cancer. Previously, we reported the discovery of several dysregulated circRNAs including circPABPC1 (polyadenylate-binding protein 1) in human hepatocellular carcinoma (HCC), although their roles in HCC development remained unclear. Here, we show that circPABPC1 is preferentially lost in tumor cells from clinical samples and inhibits both intrahepatic and distant metastases in a mouse xenograft model. This tumor-suppressive function of circPABPC1 can be attributed to its inhibition of cell adhesion and migration through down-regulating a key member of the integrin family, ITGB1 (β₁ integrin). Mass spectrometry and biochemical evidence demonstrate that circPABPC1 directly links ITGB1 to the 26S proteasome for degradation in a ubiquitination-independent manner. Our data have revealed an uncanonical route for integrin turnover and a previously unidentified mode of action for circRNAs in HCC that can be harnessed for anticancer treatment.

Differential effect of anesthetics on mucociliary clearance in vivo in mice

Respiratory mucociliary clearance (MCC) is a key defense mechanism that functions to entrap and transport inhaled pollutants, particulates, and pathogens away from the lungs. Previous work has identified a number of anesthetics to have cilia depressive effects in vitro. Wild-type C57BL/6 J mice received intra-tracheal installation of ^99mTc-Sulfur colloid, and were imaged using a dual-modality SPECT/CT system at 0 and 6 h to measure baseline MCC (n = 8). Mice were challenged for one hour with inhalational 1.5% isoflurane, or intraperitoneal ketamine (100 mg/kg)/xylazine (20 mg/kg), ketamine (0.5 mg/kg)/dexmedetomidine (50 mg/kg), fentanyl (0.2 mg/kg)/1.5% isoflurane, propofol (120 mg/Kg), or fentanyl/midazolam/dexmedetomidine (0.025 mg/kg/2.5 mg/kg/0.25 mg/kg) prior to MCC assessment. The baseline MCC was 6.4%, and was significantly reduced to 3.7% (p = 0.04) and 3.0% (p = 0.01) by ketamine/xylazine and ketamine/dexmedetomidine challenge respectively. Importantly, combinations of drugs containing fentanyl, and propofol in isolation did not significantly depress MCC. Although no change in cilia length or percent ciliation was expected, we tried to correlate ex-vivo tracheal cilia ciliary beat frequency and cilia-generated flow velocities with MCC and found no correlation. Our results indicate that anesthetics containing ketamine (ketamine/xylazine and ketamine/dexmedetomidine) significantly depress MCC, while combinations containing fentanyl (fentanyl/isoflurane, fentanyl/midazolam/dexmedetomidine) and propofol do not. Our method for assessing MCC is reproducible and has utility for studying the effects of other drug combinations.

Asymptotic Distributions of Empirical Interaction Information

Interaction Information is one of the most promising interaction strength measures with many desirable properties. However, its use for interaction detection was hindered by the fact that apart from the simple case of overall independence, asymptotic distribution of its estimate has not been known. In the paper we provide asymptotic distributions of its empirical versions which are needed for formal testing of interactions. We prove that for three-dimensional nominal vector normalized empirical interaction information converges to the normal law unless the distribution coincides with its Kirkwood approximation. In the opposite case the convergence is to the distribution of weighted centred chi square random variables. This case is of special importance as it roughly corresponds to interaction information being zero and the asymptotic distribution can be used for construction of formal tests for interaction detection. The result generalizes result in Han (Inf Control 46(1):26–45 1980) for the case when all coordinate random variables are independent. The derivation relies on studying structure of covariance matrix of asymptotic distribution and its eigenvalues. For the case of 3 × 3 × 2 contingency table corresponding to study of two interacting Single Nucleotide Polymorphisms (SNPs) for prediction of binary outcome, we provide complete description of the asymptotic law and construct approximate critical regions for testing of interactions when two SNPs are possibly dependent. We show in numerical experiments that the test based on the derived asymptotic distribution is easy to implement and yields actual significance levels consistently closer to the nominal ones than the test based on chi square reference distribution.

Expression of DNA repair genes in arctic char (Salvelinus alpinus) from Bjørnøya in the Norwegian Arctic

High levels of organochlorines (OCs) have been measured in arctic char (Salvelinus alpinus) from Lake Ellasjøen on Bjørnøya, Norway (74.30°N, 19.0°E). In a nearby lake, Laksvatn, the OC-levels in arctic char were low. A previous study has shown that char from Ellasjøen had significantly higher levels of DNA double strand breaks (DSBs) than char from Lake Laksvatn. Even though there is increasing evidence of the genotoxic effects of OCs, little is known about the effects of OCs on the DNA repair system. The aim of the present study was to determine if the two main DNA DSB repair mechanisms, homologous recombination (HR) and non-homologous end-joining (NHEJ), are affected by the higher OC and DSB level in char from Ellasjøen. This was analysed by comparing the transcript level of 11 genes involved in DNA DSB repair in char liver samples from Ellasjøen (n = 9) with char from Laksvatn (n = 12). Six of the investigated genes were significantly upregulated in char from Ellasjøen. As the expression of DNA DSB repair genes was increased in the contaminant-exposed char, it is likely that the DNA DSB repair capacity is induced in these individuals. This induction was positively correlated with the DNA DSB and negatively correlated with one or several OCs for four of these genes. However, the strongest predictor variable for DNA repair genes was habitat, indicating genetic differences in repair capacity between populations. As char from Ellasjøen still had significantly higher levels of DSBs compared to char from Laksvatn, it is possible that chronic exposure to OCs and continued production of DSB has caused selective pressure within the population for fixation of adaptive alleles. It is also possible that DSB production was exceeding the repair capacity given the prevailing conditions, or that the OC or DSB level was above the threshold value of inhibition of the DNA repair system resulting in the rate of DNA damage exceeding the rate of repair.

Elevated Polygenic Burden for Autism Spectrum Disorder Is Associated With the Broad Autism Phenotype in Mothers of Individuals With Autism Spectrum Disorder

BACKGROUND

Autism spectrum disorder (ASD) is a multifactorial neurodevelopmental disorder that encompasses a complex and heterogeneous set of traits. Subclinical traits that mirror the core features of ASD, referred to as the broad autism phenotype (BAP), have been documented repeatedly in unaffected relatives and are believed to reflect underlying genetic liability to ASD. The BAP may help inform the etiology of ASD by allowing the stratification of families into more phenotypically and etiologically homogeneous subgroups. This study explores polygenic scores related to the BAP.

METHODS

Phenotypic and genotypic information were obtained from 2614 trios from the Simons Simplex Collection. Polygenic scores of ASD (ASD-PGSs) were generated across the sample to determine the shared genetic overlap between the BAP and ASD. Maternal and paternal ASD-PGSs were explored in relation to BAP traits and their child's ASD symptomatology.

RESULTS

Maternal pragmatic language was related to child's social communicative atypicalities. In fathers, rigid personality was related to increased repetitive behaviors in children. Maternal (but not paternal) ASD-PGSs were related to the pragmatic language and rigid BAP domains.

CONCLUSIONS

Associations emerged between parent and child phenotypes, with more associations emerging in mothers than in fathers. ASD-PGS associations emerged with BAP in mothers only, highlighting the potential for a female protective factor, and implicating the polygenic etiology of ASD-related phenotypes in the BAP.

A long non-coding RNAs expression signature to improve prognostic prediction of Wilms tumor in children

BACKGROUND

Wilms tumor (WT) is the most frequent malignancy of the kidney in children, and a subset of patients remains with a poor prognosis. This study aimed to identify key long non-coding RNAs (lncRNAs) related to prognosis and establish a genomic-clinicopathologic nomogram to predict survival in children with WT.

METHODS

Clinical data of 124 WT patients and the relevant RNA sequencing data including lncRNAs expression signature of primary WT samples were obtained from the Therapeutically Applicable Research to Generate Effective Treatment (TARGET) Data Matrix. Then, lncRNAs associated with overall survival (OS) were identified through univariate Cox, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses. The risk scores of 124 participants were calculated, and survival analyses were performed between low- and high-risk groups. A genomic-clinicopathologic nomogram was then developed and evaluated by time-dependent receiver operating characteristic (ROC) curves, including the area under the curve (AUC), calibration curve, and decision curve analysis. Subsequently, bioinformatics analyses were performed to explore the potential molecular mechanisms that affect the prognosis of WT. The package "DESeq2" was used to identify differentially expressed protein-coding genes (DEPCGs) between groups. Gene Set Enrichment Analysis (GSEA) was applied to explore the differences in pathways enrichment. The analytical tools CIBERSORTx and ESTIMATE were used to investigate the discrepancies of the immune microenvironment.

RESULTS

A total of 10 lncRNAs were selected as independent predictors associated with OS (P<0.05). Participants in the high-risk group had a significantly worse OS and event-free survival (EFS) than those in the low-risk group (P<2E-16 and P=2.03E-04, respectively). The risk score and 3 clinicopathological features (gender, cooperative group protocol, and stage) were identified to construct the nomogram (combined model) (P=5.11E-17). The combined model (1-year AUC: 0.9272, 3-year AUC: 0.9428, 5-year AUC: 0.9259) and risk score model (1-year AUC: 0.9285, 3-year AUC: 0.9399, 5-year AUC: 0.9266) displayed higher predictive accuracy than that of the other models. Subsequently, 105 DEPCGs were identified. The GSEA revealed 4 significant pathways. Analysis with CIBERSORTx demonstrated that monocytes, macrophages M1, activated dendritic cells, and resting mast cells had significant infiltration differences between groups.

CONCLUSIONS

This study constructed a genomic-clinicopathologic nomogram, which might present a novel and efficient method for treating patients with WT.

Beyond monogenetic rare variants: tackling the low rate of genetic diagnoses in predominantly antibody deficiency

Predominantly antibody deficiency (PAD) is the most prevalent form of primary immunodeficiency, and is characterized by broad clinical, immunological and genetic heterogeneity. Utilizing the current gold standard of whole exome sequencing for diagnosis, pathogenic gene variants are only identified in less than 20% of patients. While elucidation of the causal genes underlying PAD has provided many insights into the cellular and molecular mechanisms underpinning disease pathogenesis, many other genes may remain as yet undefined to enable definitive diagnosis, prognostic monitoring and targeted therapy of patients. Considering that many patients display a relatively late onset of disease presentation in their 2nd or 3rd decade of life, it is questionable whether a single genetic lesion underlies disease in all patients. Potentially, combined effects of other gene variants and/or non-genetic factors, including specific infections can drive disease presentation. In this review, we define (1) the clinical and immunological variability of PAD, (2) consider how genetic defects identified in PAD have given insight into B-cell immunobiology, (3) address recent technological advances in genomics and the challenges associated with identifying causal variants, and (4) discuss how functional validation of variants of unknown significance could potentially be translated into increased diagnostic rates, improved prognostic monitoring and personalized medicine for PAD patients. A multidisciplinary approach will be the key to curtailing the early mortality and high morbidity rates in this immune disorder.

TGF-β Signaling in Pancreatic Islet β Cell Development and Function

Pancreatic islet beta cells (β-cells) synthesize and secrete insulin in response to rising glucose levels and thus are a prime target in both major forms of diabetes. Type 1 diabetes ensues due to autoimmune destruction of β-cells. On the other hand, the prevailing insulin resistance and hyperglycemia in type 2 diabetes (T2D) elicits a compensatory response from β-cells that involves increases in β-cell mass and function. However, the sustained metabolic stress results in β-cell failure, characterized by severe β-cell dysfunction and loss of β-cell mass. Dynamic changes to β-cell mass also occur during pancreatic development that involves extensive growth and morphogenesis. These orchestrated events are triggered by multiple signaling pathways, including those representing the transforming growth factor β (TGF-β) superfamily. TGF-β pathway ligands play important roles during endocrine pancreas development, β-cell proliferation, differentiation, and apoptosis. Furthermore, new findings are suggestive of TGF-β's role in regulation of adult β-cell mass and function. Collectively, these findings support the therapeutic utility of targeting TGF-β in diabetes. Summarizing the role of the various TGF-β pathway ligands in β-cell development, growth and function in normal physiology, and during diabetes pathogenesis is the topic of this mini-review.

Emerging roles of circular RNAs in systemic lupus erythematosus

Circular RNAs (circRNAs) are a class of non-coding RNAs with covalently closed single-stranded structures lacking 5'-3' polarity and a polyadenine tail. Over recent years, a growing body of studies have been conducted to explore the roles of circRNAs in human diseases. Systemic lupus erythematosus (SLE) is a severe autoimmune disorder characterized by the presence of autoantibodies and excessive inflammation, which impact multiple organs. Recent advances have begun to shed light on the roles of circRNAs in SLE, providing fresh insights into the pathogenesis of SLE and the latent capacity for translation into clinical applications. Here, we briefly introduce these "star molecules" and summarize their roles in SLE. In addition, we outline the limitations of the current studies and raise prospects for future research.

High expression of CDCA7 predicts poor prognosis for clear cell renal cell carcinoma and explores its associations with immunity

BACKGROUND

Cell division cycle-associated 7 (CDCA7), as a member of the cell division cycle associated family, was reported to be aberrantly expressed in both solid tumors and hematological tumors, suggesting its essential role in promoting tumorigenesis. Hence, we aimed to explore its comprehensive roles of overall survival (OS) in clear cell renal cell carcinoma (ccRCC) and emphasize its associations with immunity.

METHODS

The RNA sequencing data and corresponding clinical information were downloaded from The Cancer Genome Atlas (TCGA) database. Gene set enrichment analysis (GSEA) was adopted to explore CDCA7 associated signaling pathways. Univariate and multivariate Cox regression analyses were carried out to assess independent prognostic factors. Furthermore, roles of CDCA7 in human immunity were also investigated.

RESULTS

Our results suggested that CDCA7 was overexpressed in ccRCC and its elevated expression was related to shorter OS (P < 0.01). Univariate and multivariate Cox regression analyses identified CDCA7 as an independent prognostic factor (both P < 0.05). The prognostic nomogram integrating CDCA7 expression level and clinicopathologic variables was constructed to predict 1-, 3- and 5-year OS. GSEA indicated that high CDCA7 expression was related to the apoptosis pathway, cell cycle pathway, JAK-STAT pathway, NOD like receptor pathway, P53 pathway, T cell receptor pathway and toll like receptor pathway, etc. Moreover, CDCA7 was significantly related to microsatellite instability (MSI, P < 0.001) and tumor mutational burden (TMB, P < 0.001). As for immunity, CDCA7 was remarkably associated with immune infiltration, tumor microenvironment, immune checkpoint molecules and immune pathways.

CONCLUSIONS

CDCA7 could serve as an independent prognostic factor for ccRCC and it was closely related to MSI, TMB, and immunity.

Dynamics of transcription-mediated conversion from euchromatin to facultative heterochromatin at the Xist promoter by Tsix

The fine-scale dynamics from euchromatin (EC) to facultative heterochromatin (fHC) has remained largely unclear. Here, we focus on Xist and its silencing initiator Tsix as a paradigm of transcription-mediated conversion from EC to fHC. In mouse epiblast stem cells, induction of Tsix recapitulates the conversion at the Xist promoter. Investigating the dynamics reveals that the conversion proceeds in a stepwise manner. Initially, a transient opened chromatin structure is observed. In the second step, gene silencing is initiated and dependent on Tsix, which is reversible and accompanied by simultaneous changes in multiple histone modifications. At the last step, maintenance of silencing becomes independent of Tsix and irreversible, which correlates with occupation of the -1 position of the transcription start site by a nucleosome and initiation of DNA methylation introduction. This study highlights the hierarchy of multiple chromatin events upon stepwise gene silencing establishment.

At3g53630 encodes a GUN1-interacting protein under norflurazon treatment

Chloroplasts are semi-autonomous organelles, with more than 95% of their proteins encoded by the nuclear genome. The chloroplast-to-nucleus retrograde signals are critical for the nucleus to coordinate its gene expression for optimizing or repairing chloroplast functions in response to changing environments. In chloroplasts, the pentatricopeptide-repeat protein GENOMES UNCOUPLED 1 (GUN1) is a master switch that senses aberrant physiological states, such as the photooxidative stress induced by norflurazon (NF) treatment, and represses the expression of photosynthesis-associated nuclear genes (PhANGs). However, it is largely unknown how the retrograde signal is transmitted beyond GUN1. In this study, a protein GUN1-INTERACTING PROTEIN 1 (GIP1), encoded by At3g53630, was identified to interact with GUN1 by different approaches. We demonstrated that GIP1 has both cytosol and chloroplast localizations, and its abundance in chloroplasts is enhanced by NF treatment with the presence of GUN1. Our results suggest that GIP1 and GUN1 may function antagonistically in the retrograde signaling pathway.

Pairwise linkage disequilibrium estimation for polyploids

Many tasks in statistical genetics involve pairwise estimation of linkage disequilibrium (LD). The study of LD in diploids is mature. However, in polyploids, the field lacks a comprehensive characterization of LD. Polyploids also exhibit greater levels of genotype uncertainty than diploids, yet no methods currently exist to estimate LD in polyploids in the presence of such genotype uncertainty. Furthermore, most LD estimation methods do not quantify the level of uncertainty in their LD estimates. Our study contains three major contributions. (i) We characterize haplotypic and composite measures of LD in polyploids. These composite measures of LD turn out to be functions of common statistical measures of association. (ii) We derive procedures to estimate haplotypic and composite LD in polyploids in the presence of genotype uncertainty. We do this by estimating LD directly from genotype likelihoods, which may be obtained from many genotyping platforms. (iii) We derive standard errors of all LD estimators that we discuss. We validate our methods on both real and simulated data. Our methods are implemented in the R package ldsep, available on the Comprehensive R Archive Network https://cran.r-project.org/package=ldsep.

Structural Racism in the COVID-19 Pandemic: Moving Forward

The COVID-19 pandemic has taken a substantial human, social and economic toll globally, but its impact on Black/African Americans, Latinx, and American Indian/Alaska Native communities in the U.S. is unconscionable. As the U.S. continues to combat the current COVID-19 cycle and prepares for future pandemics, it will be critical to learn from and rectify past and contemporary wrongs. Drawing on experiences in genomic research and intersecting areas in medical ethics, health disparities, and human rights, this article considers three key COVID-19-related issues: research to identify remedies; testing, contact tracing and surveillance; and lingering health needs and disability. It provides a pathway for the future: community engagement to develop culturally-sensitive responses to the myriad genomic/bioethical dilemmas that arise, and the establishment of a Truth and Reconciliation Commission to transition the country from its contemporary state of segregation in healthcare and health outcomes into an equitable and prosperous society for all.

Cognitive analysis of metabolomics data for systems biology

Cognitive computing is revolutionizing the way big data are processed and integrated, with artificial intelligence (AI) natural language processing (NLP) platforms helping researchers to efficiently search and digest the vast scientific literature. Most available platforms have been developed for biomedical researchers, but new NLP tools are emerging for biologists in other fields and an important example is metabolomics. NLP provides literature-based contextualization of metabolic features that decreases the time and expert-level subject knowledge required during the prioritization, identification and interpretation steps in the metabolomics data analysis pipeline. Here, we describe and demonstrate four workflows that combine metabolomics data with NLP-based literature searches of scientific databases to aid in the analysis of metabolomics data and their biological interpretation. The four procedures can be used in isolation or consecutively, depending on the research questions. The first, used for initial metabolite annotation and prioritization, creates a list of metabolites that would be interesting for follow-up. The second workflow finds literature evidence of the activity of metabolites and metabolic pathways in governing the biological condition on a systems biology level. The third is used to identify candidate biomarkers, and the fourth looks for metabolic conditions or drug-repurposing targets that the two diseases have in common. The protocol can take 1-4 h or more to complete, depending on the processing time of the various software used.

Shortened relative leukocyte telomere length is associated with all-cause mortality in type 2 diabetes- analysis from the Hong Kong Diabetes Register

AIMS

Few studies have investigated the relationship between rLTL and mortality in patients with type 2 diabetes in a large prospective study, particularly in the Asian population. This study investigates the relationship between rLTL and the risk of death in Chinese patients with type 2 diabetes.

METHODS

Consecutive Chinese patients with type 2 diabetes (N = 5349) from the Hong Kong Diabetes Register with stored baseline DNA and available follow-up data were studied. rLTL was measured using a quantitative polymerase chain reaction. Mortality and clinical outcomes were obtained based on ICD-9 codes.

RESULTS

The mean (SD) age of the subjects was 57.5 (13.3) years and mean (SD) follow-up duration was 13.4 (5.5) years. Baseline rLTL was significantly shorter in the 1925 subjects who subsequently died compared with the remaining subjects (4.3 ± 1.2 vs. 4.7 ± 1.2, P < 0.001). Shorter rLTL was associated with a higher risk of mortality (HR: 1.19 (1.14-1.23), P < 0.001), which remained significant after adjusting for traditional risk factors.

CONCLUSIONS

Shorter rLTL was significantly associated with an increased risk of all-cause and CVD mortality in patients with type 2 diabetes, independent of established risk factors. Telomere length may be a useful biomarker for mortality risk in patients with type 2 diabetes.

METTL3-mediated m6A modification regulates cell cycle progression of dental pulp stem cells

Background

Dental pulp stem cells (DPSCs) are a promising cell source in endodontic regeneration and tissue engineering with limited self-renewal and pluripotency capacity. N⁶-methyladenosine (m⁶A) is the most prevalent, reversible internal modification in RNAs associated with stem cell fate determination. In this study, we aim to explore the biological effect of m⁶A methylation in DPSCs.

Methods

m⁶A immunoprecipitation with deep sequencing (m⁶A RIP-seq) demonstrated the features of m⁶A modifications in DPSC transcriptome. Lentiviral vectors were constructed to knockdown or overexpress methyltransferase like 3 (METTL3). Cell morphology, viability, senescence, and apoptosis were analyzed by β-galactosidase, TUNEL staining, and flow cytometry. Bioinformatic analysis combing m⁶A RIP and shMETTL3 RNA-seq functionally enriched overlapped genes and screened target of METTL3. Cell cycle distributions were assayed by flow cytometry, and m⁶A RIP-qPCR was used to confirm METTL3-mediated m⁶A methylation.

Results

Here, m⁶A peak distribution, binding area, and motif in DPSCs were first revealed by m⁶A RIP-seq. We also found a relatively high expression level of METTL3 in immature DPSCs with superior regenerative potential and METTL3 knockdown induced cell apoptosis and senescence. A conjoint analysis of m⁶A RIP and RNA sequencing showed METTL3 depletion associated with cell cycle, mitosis, and alteration of METTL3 resulted in cell cycle arrest. Furthermore, the protein interaction network of differentially expressed genes identified Polo-like kinase 1 (PLK1), a critical cycle modulator, as the target of METTL3-mediated m⁶A methylation in DPSCs.

Conclusions

These results revealed m⁶A methylated hallmarks in DPSCs and a regulatory role of METTL3 in cell cycle control. Our study shed light on therapeutic approaches in vital pulp therapy and served new insight into stem cell-based tissue engineering.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-021-02223-x.

Impact of integrated translational research on clinical exome sequencing

PURPOSE

Exome sequencing often identifies pathogenic genetic variants in patients with undiagnosed diseases. Nevertheless, frequent findings of variants of uncertain significance necessitate additional efforts to establish causality before reaching a conclusive diagnosis. To provide comprehensive genomic testing to patients with undiagnosed disease, we established an Individualized Medicine Clinic, which offered clinical exome testing and included a Translational Omics Program (TOP) that provided variant curation, research activities, or research exome sequencing.

METHODS

From 2012 to 2018, 1101 unselected patients with undiagnosed diseases received exome testing. Outcomes were reviewed to assess impact of the TOP and patient characteristics on diagnostic rates through descriptive and multivariate analyses.

RESULTS

The overall diagnostic yield was 24.9% (274 of 1101 patients), with 174 (15.8% of 1101) diagnosed on the basis of clinical exome sequencing alone. Four hundred twenty-three patients with nondiagnostic or without access to clinical exome sequencing were evaluated by the TOP, with 100 (9% of 1101) patients receiving a diagnosis, accounting for 36.5% of the diagnostic yield. The identification of a genetic diagnosis was influenced by the age at time of testing and the disease phenotype of the patient.

CONCLUSION

Integration of translational research activities into clinical practice of a tertiary medical center can significantly increase the diagnostic yield of patients with undiagnosed disease.

FAM72, Glioblastoma Multiforme (GBM) and Beyond

Simple Summary

Glioblastoma multiforme (GBM) is a serious and aggressive cancer disease that has not allowed scientists to rest for decades. In this review, we consider the new gene pair |-SRGAP2–FAM72-| and discuss its role in the cell cycle and the possibility of defining new therapeutic approaches for the treatment of GBM and other cancers via this gene pair |-SRGAP2–FAM72-|.

Abstract

Neural stem cells (NSCs) offer great potential for regenerative medicine due to their excellent ability to differentiate into various specialized cell types of the brain. In the central nervous system (CNS), NSC renewal and differentiation are under strict control by the regulation of the pivotal SLIT-ROBO Rho GTPase activating protein 2 (SRGAP2)—Family with sequence similarity 72 (FAM72) master gene (i.e., |-SRGAP2–FAM72-|) via a divergent gene transcription activation mechanism. If the gene transcription control unit (i.e., the intergenic region of the two sub-gene units, SRGAP2 and FAM72) gets out of control, NSCs may transform into cancer stem cells and generate brain tumor cells responsible for brain cancer such as glioblastoma multiforme (GBM). Here, we discuss the surveillance of this |-SRGAP2–FAM72-| master gene and its role in GBM, and also in light of FAM72 for diagnosing various types of cancers outside of the CNS.

A draft phased assembly of the diploid Cascade hop (Humulus lupulus) genome

Hop (Humulus lupulus L. var Lupulus) is a diploid, dioecious plant with a history of cultivation spanning more than one thousand years. Hop cones are valued for their use in brewing and contain compounds of therapeutic interest including xanthohumol. Efforts to determine how biochemical pathways responsible for desirable traits are regulated have been challenged by the large (2.8 Gb), repetitive, and heterozygous genome of hop. We present a draft haplotype-phased assembly of the Cascade cultivar genome. Our draft assembly and annotation of the Cascade genome is the most extensive representation of the hop genome to date. PacBio long-read sequences from hop were assembled with FALCON and partially phased with FALCON-Unzip. Comparative analysis of haplotype sequences provides insight into selective pressures that have driven evolution in hop. We discovered genes with greater sequence divergence enriched for stress-response, growth, and flowering functions in the draft phased assembly. With improved resolution of long terminal retrotransposons (LTRs) due to long-read sequencing, we found that hop is over 70% repetitive. We identified a homolog of cannabidiolic acid synthase (CBDAS) that is expressed in multiple tissues. The approaches we developed to analyze the draft phased assembly serve to deepen our understanding of the genomic landscape of hop and may have broader applicability to the study of other large, complex genomes.

LncRNA-SNHG7/miR-29b/DNMT3A axis affects activation, autophagy and proliferation of hepatic stellate cells in liver fibrosis

OBJECTIVES

To determine the relative expression of long non-coding small nucleolar RNA host gene 7 (lncRNA-SNHG7) in fibrotic liver and hepatic stellate cells, and investigate the biological effects and mechanisms of SNHG7 on liver fibrosis.

METHODS

Liver fibrosis model of mice was established, primary hepatic stellate cells (HSCs) were cultured from normal mice and induced to activate by TGF-β. Cell viability, proliferation, and autophagy were detected with MTT, BrdU, and MDC stain, respectively. Liver tissue stained with Masson and Sirius Red. The interaction between SNHG7 and miR-29b was investigated by immunoprecipitation, RNA pull-down and Dual-luciferase reporter gene assay. The effects of SNHG7 on the expression of miR-29b, DNMT3A and liver fibrosis related factors were detected in vitro or in vivo transfection experiments.

RESULTS

SNHG7 was signally increased in liver tissue and HSCs of liver fibrosis model of mice, and inhibition of SNHG7 expression in liver fibrosis mice can reduce liver fibrosis. We also found that SNHG7 could bind to miR-29b in HSCs and inhibit the expression of miR-29b. In TGF-β-stimulated normal HSCs, knockdown of SNHG7 expression after shSNHG7 transfection could restrain DNMT3A and HSCs activation factors α-SMA, Collα1 and autophagy-related factors LC3I/II, Beclin1. However, this shSNHG7 effect was reversed by the inhibition of miR-29b.

CONCLUSION

Inhibition of lncRNA-SNHG7 can inhibit liver fibrosis. This is partly due to SNHG7 acts as a competitive endogenous RNA (ceRNA) to affect the expression of DNMT3A, a downstream target gene of miR-29b, by binding to miR-29b, thereby affecting the activation, autophagy and proliferation of HSCs.

Network Analysis Reveals Proteins Associated with Aortic Dilatation in Mucopolysaccharidoses

Mucopolysaccharidoses are caused by a deficiency of enzymes involved in the degradation of glycosaminoglycans. Heart diseases are a significant cause of morbidity and mortality in MPS patients, even in conditions in which enzyme replacement therapy is available. In this sense, cardiovascular manifestations, such as heart hypertrophy, cardiac function reduction, increased left ventricular chamber, and aortic dilatation, are among the most frequent. However, the downstream events which influence the heart dilatation process are unclear. Here, we employed systems biology tools together with transcriptomic data to investigate new elements that may be involved in aortic dilatation in Mucopolysaccharidoses syndrome. We identified candidate genes involved in biological processes related to inflammatory responses, deposition of collagen, and lipid accumulation in the cardiovascular system that may be involved in aortic dilatation in the Mucopolysaccharidoses I and VII. Furthermore, we investigated the molecular mechanisms of losartan treatment in Mucopolysaccharidoses I mice to underscore how this drug acts to prevent aortic dilation. Our data indicate that the association between the TGF-b signaling pathway, Fos, and Col1a1 proteins can play an essential role in aortic dilation's pathophysiology and its subsequent improvement by losartan treatment.

Wound Healing-related Functions of the p160 Steroid Receptor Coactivator Family

Multicellular organisms have evolved sophisticated mechanisms to recover and maintain original tissue functions following injury. Injury responses require a robust transcriptomic response associated with cellular reprogramming involving complex gene expression programs critical for effective tissue repair following injury. Steroid receptor coactivators (SRCs) are master transcriptional regulators of cell-cell signaling that is integral for embryogenesis, reproduction, normal physiological function, and tissue repair following injury. Effective therapeutic approaches for facilitating improved tissue regeneration and repair will likely involve temporal and combinatorial manipulation of cell-intrinsic and cell-extrinsic factors. Pleiotropic actions of SRCs that are critical for wound healing range from immune regulation and angiogenesis to maintenance of metabolic regulation in diverse organ systems. Recent evidence derived from studies of model organisms during different developmental stages indicates the importance of the interplay of immune cells and stromal cells to wound healing. With SRCs being the master regulators of cell-cell signaling integral to physiologic changes necessary for wound repair, it is becoming clear that therapeutic targeting of SRCs provides a unique opportunity for drug development in wound healing. This review will provide an overview of wound healing-related functions of SRCs with a special focus on cellular and molecular interactions important for limiting tissue damage after injury. Finally, we review recent findings showing stimulation of SRCs following cardiac injury with the SRC small molecule stimulator MCB-613 can promote cardiac protection and inhibit pathologic remodeling after myocardial infarction.