Removal of the proliferation component of gene expression by PCNA adjustment has been addressed in numerous survival prediction studies for breast cancer and all cancers in the TCGA. These studies indicate that widespread co-regulation of proliferation upwardly biases survival prediction when gene selection is performed on a genome-wide basis. In addition, removal of the correlative effects of proliferation does not reduce the random bias associated with survival prediction using random gene selection. Since most cancers become addicted to DNA repair as a result of forced cellular replication, increased oxidation, and repair deficiencies from oncogenic loss or genetic polymorphisms, we pursued an investigation to remove the proliferation component of expression in DNA repair genes to determine survival prediction. This translational hypothesis-driven focus on DNA repair genes is directly amenable to finding new sets of DNA repair genes that could potentially be studied for inhibition therapy. Overall survival (OS) prediction was evaluated in 18 cancers by using normalized RNA-Seq data for 126 DNA repair genes with expression available in TCGA. Transformations for normality and adjustments for age at diagnosis, stage, and PCNA metagene expression were performed for all DNA repair genes. We also analyzed genomic event rates (GER) for somatic mutations, deletions, and amplification in driver genes and DNA repair genes. After performing empirical p-value testing with use of randomly selected gene sets, it was observed that OS could be predicted significantly by sets of DNA repair genes for 61% (11/18) of the cancers. Interestingly, PARP1 was not a significant predictor of survival for any of the 11 cancers. Results from cluster analysis of GERs indicates that the most opportunistic cancers for inhibition therapy may be AML, colorectal, and renal papillary, because of potentially less confounding due to lower GERs for mutations, deletions, and amplifications in DNA repair genes. However, the most opportunistic cancer for inhibition therapy is likely to be AML, since it showed the lowest GERs for mutations, deletions, and amplifications in DNA repair genes. In conclusion, our hypothesis-driven focus to target DNA repair gene expression adjusted for the PCNA metagene as a means of predicting OS in various cancers resulted in statistically significant sets of genes.

Geckos demonstrate a remarkable variability in sex determination systems, but our limited knowledge prohibits accurate conclusions on the evolution of sex determination in this group. Eyelid geckos (Eublepharidae) are of particular interest, as they encompass species with both environmental and genotypic sex determination. We identified for the first time the X-specific gene content in the Yucatán banded gecko, Coleonyx elegans, possessing X1X1X2X2/X1X2Y multiple sex chromosomes by comparative genome coverage analysis between sexes. The X-specific gene content of Coleonyx elegans was revealed to be partially homologous to genomic regions linked to the chicken autosomes 1, 6 and 11. A qPCR-based test was applied to validate a subset of X-specific genes by comparing the difference in gene copy numbers between sexes, and to explore the homology of sex chromosomes across 11 eublepharid, two phyllodactylid and one sphaerodactylid species. Homologous sex chromosomes are shared between Coleonyx elegans and Coleonyx mitratus, two species diverged approximately 34 million years ago, but not with other tested species. As far as we know, the X-specific gene content of Coleonyx elegans / Coleonyx mitratus was never involved in the sex chromosomes of other gecko lineages, indicating that the sex chromosomes in this clade of eublepharid geckos evolved independently.

To guide analysts to select the right tool and parameters in differential gene expression analysis of single-cell RNA sequencing (scRNA-seq) data, we developed a novel simulator that recapitulates the data characteristics of real scRNA-seq datasets while accounting for all the relevant sources of variation in a multi-subject, multi-condition scRNA-seq experiment: the cell-to-cell variation within a subject, the variation across subjects, the variability across cell types, the mean/variance relationship of gene expression across genes, library size effects, group effects, and covariate effects. By applying it to benchmark 12 differential gene expression analysis methods (including cell-level and pseudo-bulk methods) on simulated multi-condition, multi-subject data of the 10x Genomics platform, we demonstrated that methods originating from the negative binomial mixed model such as glmmTMB and NEBULA-HL outperformed other methods. Utilizing NEBULA-HL in a statistical analysis pipeline (https://github.com/interactivereport/scRNAseq_DE) for single cell analysis will enable scientists to better understand cell-type specific transcriptomic response to disease or treatment effects and to discover new drug targets. Further, application to two real datasets showed the outperformance of our differential expression (DE) pipeline, with unified findings of differentially expressed genes (DEG) and a pseudo-time trajectory transcriptomic result. In the end, we made recommendations of filtering strategies of cells and genes based on simulation results to achieve optimal experimental goals.

The Drosophila imaginal disc has been an excellent model for the study of developmental gene regulation. In particular, long non-coding RNAs (lncRNAs) have gained widespread attention in recent years due to their important role in gene regulation. Their specific spatiotemporal expressions further support their role in developmental processes and diseases. In this study, we explored the role of a novel lncRNA in Drosophila leg development by dissecting and dissociating w1118 third-instar larval third leg (L3) discs into single cells and single nuclei, and performing single-cell RNA-sequencing (scRNA-seq) and single-cell assays for transposase-accessible chromatin (scATAC-seq). Single-cell transcriptomics analysis of the L3 discs across three developmental timepoints revealed different cell types and identified lncRNA:CR33938 as a distal specific gene with high expression in late development. This was further validated by fluorescence in-situ hybridization (FISH). The scATAC-seq results reproduced the single-cell transcriptomics landscape and elucidated the distal cell functions at different timepoints. Furthermore, overexpression of lncRNA:CR33938 in the S2 cell line increased the expression of leg development genes, further confirming its important role in development.

The snATAC + snRNA platform allows epigenomic profiling of open chromatin and gene expression with single-cell resolution. The most critical assay step is to isolate high-quality nuclei to proceed with droplet-base single nuclei isolation and barcoding. With the increasing popularity of multiomic profiling in various fields, there is a need for optimized and reliable nuclei isolation methods, mainly for human tissue samples. Herein we compared different nuclei isolation methods for cell suspensions, such as peripheral blood mononuclear cells (PBMC, n=18) and a solid tumor type, ovarian cancer (OC) (n=18), derived from debulking surgery. Nuclei morphology and sequencing output parameters were used to evaluate the quality of preparation. Our results show that NP-40 detergent-based nuclei isolation yields better sequencing results than collagenase tissue dissociation for OC, significantly impacting cell identification and analysis. Given the utility of applying such techniques to frozen samples, we also tested frozen preparation and digestion (n=6). A paired comparison between frozen and fresh samples validated the quality of both specimens. Finally, we demonstrate the reproducibility of scRNA and snATAC + snRNA platform, by comparing the gene expression profiling of PBMC. Our results highlight how the choice of nuclei isolation methods is critical for obtaining quality data in multiomic assays. It also shows that the measurement of expression between scRNA and snRNA is comparable and effective for cell type identification.

Plant life-long organogenesis involves sequential, time and tissue specific expression of developmental genes. This requires activities of Polycomb Group (PcG) and trithorax Group complexes, respectively responsible for repressive Histone 3 trimethylation at lysine 27 (H3K27me3) and activation-related H3K4me3. However, the genome-wide dynamics in histone modifications that occur during developmental processes have remained elusive. Here, we report the distributions of H3K27me3 and H3K4me3 along with transcriptional changes, in a developmental series including Arabidopsis leaf and three stages of flower development. We found that chromatin mark levels are highly dynamic over the time series on nearly half of all Arabidopsis genes. Moreover, during early flower morphogenesis, changes in H3K4me3 prime over changes in H3K27me3 and quantitatively correlate with transcription changes, while H3K27me3 changes occur after prolonged expression changes. Notably, early activation of PcG target genes is dominated by increases in H3K4me3 while H3K27me3 remains present at the locus. Our results reveal H3K4me3 as greater predictor over H3K27me3 for transcription dynamics, unveil unexpected chromatin mechanisms at gene activation and underline the relevance of tissue-specific temporal epigenomics.

Trichomes are a common morphological defense against pests, in particular, type IV glandular trichomes have been associated with resistance against different invertebrates. Cultivated tomatoes usually lack or have a very low density of type IV trichomes. Thus, specific breeding programs to incorporate these natural defences, that are common within the Solanum genus, might improve a more sustainable management. We have identified a S. pimpinellifolium accession with very high density of this type of trichomes. Two F2 mapping populations using two different parents have been developed, characterized and genotyped using a new genotype methodology, K-seq. We have been able to build an ultra-dense genetic map with 147,326 markers with an average distance between markers of 0.2 cM that has allowed us to perform a detailed mapping. We have used two different families and two different approaches, QTL mapping and QTL-seq, to identify several QTLs implicated in the control of trichome type IV developed in this accession on the chromosomes 5, 6, 9 and 11. The QTL located on chromosome 9 is a major QTL that has not been previously reported in S. pimpinellifolium that increases by a factor of 9 the density of trichomes.

Circular RNAs are a special type of RNAs which recently attracted a lot of research interest in studying its formation and function. RNA binding proteins (RBPs) that bind circRNAs are important in these processes but are relatively less studied. CLIP-Seq technology has been invented and applied to profile RBP-RNA interactions on the genome-wide scale. While mRNAs are usually the focus of CLIP-Seq experiments, RBP-circRNA interactions could also be identified through specialized analysis of CLIP-Seq datasets. However, many technical difficulties are involved in this process, such as the usually short read length of CLIP-Seq reads. In this study, we created a pipeline called Clirc specialized for profiling circRNAs in CLIP-Seq data and analyzing the characteristics of RBP- circRNAs interactions. In conclusion, this is one of the first few studies to investigate circRNAs and their binding partners through repurposing CLIP-Seq datasets to our knowledge, and we hope our work will become a valuable resource for future studies into the biogenesis and function of circRNAs. Clirc software is available at https://github.com/Minzhe/Clirc

Natural killer cells are associated with the pathogenesis of ulcerative colitis (UC), but their precise contributions remain unclear. The present study sought to investigate the diagnostic value of activated NK-associated genes (ANAGs) in UC. Bulk RNA-seq and scRNA-seq datasets were obtained from the Gene Expression Omnibus (GEO) and Single Cell Portal (SCP) databases. In the bulk RNA-seq, 92 differentially expressed genes (DEGs) were screened out by the “Batch correction” and “Robust rank aggregation” (RRA) methods. The immune infiltration landscape was estimated by single-sample gene set enrichment analysis (ssGSEA) and CIBERSORT, which revealed a higher abundance of activated NK cells in noninflamed UC tissues. 54 DEGs correlated with activated NK cells were identified as ANAGs. Protein-protein interaction (PPI) analysis and least absolute shrinkage and selection operator (LASSO) regression were utilized to screen out 4 key ANAGs (SELP, TIMP1, MMP7, and ABCG2) and establish an activated NK-associated gene score (ANAG score). The ANAG score demonstrated excellent diagnostic value and was validated in three external datasets. The expression of the 4 key ANAGs was validated in UC patients and healthy controls (HC) samples. Through scRNA-seq data analysis, higher expression levels of SELP, TIMP1, MMP7, and ABCG2 were observed in post-capillary venules, inflammatory fibroblasts, enterocytes, and immature enterocytes. The cell scores based on the ANAGs showed enrichment in endothelial cells and fibroblasts. In conclusion, we established and validated an ANAG score with the ability to precisely diagnose UC. The 4 key ANAGs have the potential to serve as therapeutic targets in UC.

Flood forecasting helps anticipate floods and evacuate people, but due to the access of a large number of iot data acquisition devices, the explosive growth of multidimensional data and the increasingly demanding prediction accuracy, classical parameter models and traditional machine learning algorithms are unable to meet the high efficiency and high precision requirements of prediction tasks. In recent years, deep learning algorithms represented by convolutional neural networks, recurrent neural networks and Informer models have achieved fruitful results in time series prediction tasks. The Informer model is used to predict the flood flow of the reservoir. At the same time, the prediction results are compared with the prediction results of the traditional method and the LSTM model, and how to apply the Informer model in the field of flood prediction to improve the accuracy of flood prediction is studied. The data of 28 floods in the Wan 'an Reservoir control basin from May 2014 to June 2020 were used, with areal rainfall in five subzones and outflow from two reservoirs as inputs and flood processes with different sequence lengths as outputs. The results show that the Informer model has good accuracy and applicability in flood forecasting. In the flood forecasting with sequence length of 4, 5 and 6, Informer has higher prediction accuracy, and the prediction accuracy is better than other models under the same sequence length, but the prediction accuracy will decline to a certain extent with the increase of sequence length. The Informer model stably predicts the flood peak better, and its average flood peak difference and average maximum flood peak difference are the smallest. As the length of the sequence increases, the number of fields with a maximum flood peak difference less than 15% increases, and the maximum flood peak difference decreases. Therefore, the Informer model can be used as one of the better flood forecasting methods, and it provides a new forecasting method and scientific decision-making basis for reservoir flood control.

The exogenous application of phenolic compounds is increasingly recognized as a valuable strategy for promoting growth and mitigating the adverse effects of abiotic stress. However, the biostimulant effect under optimal conditions have not been thoroughly explored. In this study, we investigated the impact of foliar application of flavonoids, specifically Crop Bio Life (CBL), on tomato plants grown under controlled conditions. Our study focused on determining growth parameters, such as cell size, and assessing the concentration of hormones. Principal component analysis (PCA) from all physiological variables were determined. Additionally, we utilized high-throughput mRNA sequencing technology and bioinformatic methodologies to robustly analyze the transcriptomes of tomato leaves regulated by flavonoids. The findings revealed that CBL primarily influenced cells enlargement, leading to increased growth. Furthermore, CBL-treated plants exhibited higher concentrations of the hormone zeatin but lower concentrations of IAA. Moreover, RNAseq analysis indicated that CBL-treated plants required increased mineral transport and water uptake, as evidenced by gene expression patterns. Genes related to pathways such as fatty acid degradation, phenylpropanoid biosynthesis, and ABC transporters showed regulatory mechanisms governing internal flavonoid biosynthesis, transport, and tissue concentration, ultimately resulting in higher flavonoid concentrations in tomato leaves.

This study explored the effects of pBD114 (porcine β-defensin 114) on inflammatory response. The results in RAW264.7 cells showed that 100 μg/mL of pBD114 significantly increased the concentrations of TNFα and IL-10, and 10 μg/mL of pBD114 significantly decreased the concen-tration of TNFα in LPS-induced RAW264.7 cell cultures at 10 ng/mL. These results suggest that pBD114 can exhibit pro-inflammatory or anti-inflammatory activities under certain conditions. RNA-seq analysis was performed to gain further insight into the effects of pBD114 on the in-flammatory response. Among pBD114 promoting RAW264.7 pro-inflammatory response, pBD114 significantly up-regulated 1170 genes and down-regulated 724 genes. KEGG enrichment found that the DEGs were significantly enriched in immune- and signal-transduction-related signaling pathways. PPI and KDA analyses revealed that Bcl10 and Bcl3 were the key genes. In addition, pBD114 significantly up-regulated 12 genes and down-regulated 38 genes in the an-ti-inflammatory response. KEGG enrichment analysis revealed that the DEGs were mainly en-riched in the "Cytokine-cytokine receptor interaction" signaling pathway, and PPI and KDA analyses showed that Stat1 and Csf2 were the key genes. The results of qRT-PCR verified those of RNA-seq. In vivo mice tests also confirmed the pro- or anti-inflammatory activity of pBD114. Although the inflammatory response is a rapid and complex physiological reaction to noxious stimuli, this study found that pBD114 plays an essential role in the inflammatory response mainly by acting on genes related to immunity, signal transduction, signaling molecules, and interactions.

Cervical cancer presents a significant global health concern with high-risk human papilloma-viruses (HPVs) being identified as the cause of this cancer. Although current treatment methods for cervical cancer can eliminate lesions, preventing metastatic spread and minimizing tissue damage remains a major challenge. Therefore, the development of a safer and innovative thera-peutic approach is of utmost importance. Natural products like Fig latex, derived from the Ficus Carica tree, have demonstrated promising anti-cancer properties when tested on cervical cancer cell lines. However, the specific mechanisms by which Fig latex exerts its effects are still un-known. In this study, we conducted RNA-seq analysis to explore how Fig latex may counteract carcinogenesis in HPV-positive cervical cancer cell lines, namely CaSki (HPV type 16-positive) and HeLa (HPV type 18-positive). Our data from this investigation indicates that Fig latex influ-ences the expression of genes associated with the development and progression of cervical can-cer, including pathways related to " Nonsense-Mediated Decay (NMD) "cell cycle regulation" and Transcriptional Regulation by TP53. This selective impact of Fig latex on cancer-related pathways suggests a potential novel therapeutic approach for HPV-related cervical cancer.

Alternative splicing can produce transcripts that affect cancer development and shows potential for cancer diagnosis and treatment. However, intron retention (IR), a type of alternative splicing, has been less systematically studied in cancer biology research. Here, we generated a pan-cancer IR landscape for more than 10,000 samples across 33 cancer types from The Cancer Genome Atlas (TCGA). We characterized differentially retained introns between tumor and normal samples and identified retained introns associated with survival. We discovered 988 differentially retained introns in 14 cancers, some of which demonstrated diagnostic potential in multiple cancer types. We also inferred a large number of prognosis-related introns in 33 cancer types, and the associated genes included well-known cancer hallmarks such as angiogenesis, metastasis, and DNA mutations. Notably, we discovered a novel intron retention event inside 5′UTR of STN1 that is associated with the survival of lung cancer patients. The retained intron reduces translation efficiency by producing upstream open reading frames (uORFs) and thereby inhibits colony formation and cell migration of lung cancer cells. Besides, the IR-based prognostic model achieved good stratification on certain cancers, as illustrated in acute myeloid leukemia. Taken together, we performed a comprehensive IR survey at a pan-cancer level, and the results implied that IR has the potential to be diagnostic and prognostic cancer biomarkers, as well as new drug targets.

The transcriptomic analysis of the dehydration rate of mature rice seeds was conducted to explore candidate genes related to the dehydration rate and provide a theoretical basis for breeding and utilization. In this study, we selected Baghlani Nangarhar, an extremely fast dehydration material, and Saturn, a slow dehydration material, based on the results of the previous studies of screening of 165 germplasm materials for dehydration rate phenotypes. Fast dehydration experiment on these two types of seeds was conducted. Four comparative groups were set up under control and dehydration conditions. The differentially expressed genes (DEGs) were quantified by transcriptome sequencing and quantitative real-time PCR (qRT-PCR). GO and KEGG analyses were carried out. The results showed that in Baghlani Nangarhar, 53 DEGs were screened, of which 33 were up-regulated and 20 were down-regulated. Similarly, in Saturn, 25 DEGs were screened, of which 19 were up-regulated and 6 down-regulated. The results of the GO analysis showed that the sites of action of the differentially expressed genes enriched in the fast dehydration modes were concentrated in the cytoplasm, internal components of the membrane and nucleosomes. They play regulatory roles in catalysis, binding, translocation, transcription, protein folding, degradation and replication. They are involved in adaptive responses to adverse external environments such as reactive oxygen species and high temperature. The KEGG analysis showed that the main metabolic pathways enriched are protein processing in the endoplasmic reticulum, amino acid biosynthesis, and oxidative phosphorylation. The key differentially expressed genes and the most important metabolic pathways in the rapidly and slowly dehydrated materials obtained in this study were protein processing in the endoplasmic reticulum and oxidative phosphorylation metabolism, which are presumed to have important regulatory roles in stress/defense, energy metabolism, protein synthesis/folding, and signal transduction during dehydration and drying of mature seeds. The results of this study could potentially provide a valuable reference for further research on the genes and metabolic pathways related to the dehydration rate of mature rice seeds and provide theoretical guidance for the selection and breeding of new rice germplasm that can be rapidly dehydrated at the mature stage.

Caseous lymphadenitis (CL) is a chronic contagious disease that affects small ruminants and is characterized by the formation of pyogranulomas in lymph nodes and other organs. However, the pathogenesis of this disease and the response of the host genome to infection are not yet fully understood. This study aimed to investigate the whole blood transcriptome and evaluate differential gene expression during the later stages of CL in naturally infected ewes. The study included diseased, serologically positive (EP), exposed, serologically negative (EN) ewes from the same infected flock and healthy ewes (CN) from a different flock. RNA sequencing was performed using the Illumina NextSeq system, and differential gene expression was estimated using DESeq2 and Edge R approaches. The analysis identified 191 annotated differentially expressed genes (DEGs) in the EP group (102 upregulated and 89 downregulated) and 256 DEGs in the EN group (106 upregulated and 150 downregulated). Infection influenced numerous immunoregulatory interactions between lymphoid and nonlymphoid cells in both EP and EN ewes. Immune DEGs were preferentially assigned to antigen presentation through the MHC complex, T-lymphocyte mediated immunity, and extracellular matrix interactions. Furthermore, the EP group showed altered regulation of cytokine and chemokine signalling and activation and recombination of B-cell receptors. Conversely, NF-kappa B signalling, apoptosis, and stress response were the main processes influenced in the EN group. In addition, statistically significant enrichment of the essential immune pathways of Binding and Uptake of Ligands by Scavenger Receptors in EP and p53 signalling in the EN group was found. In conclusion, this study provides new insights into the disease course and host‒pathogen interaction in naturally CL-infected sheep by investigating the blood transcriptome.

Brown-marbled grouper is one of the most important mariculture species in China, which is used as an important crossbreeding parent in grouper industry. Enhancing growth rates is a key target in fish breeding, and gaining insight into the underlying mechanisms responsible for growth differences among individuals can aid in the improvement of grouper growth rates. However, the mechanism behind this difference in growth in this fish is unclear. The difference of transcriptome profiles of muscle tissue between fast- and slow-growing brown-marbled grouper was analyzed by RNA-Seq. 77 significantly up-regulated genes and 92 significantly down-regulated genes were identified in the growth extreme groups. The up-regulated of ghr and tnni2 and the down-regulated of stc2 led to the growth advantages of brown-marbled grouper. The differently expressed genes (DEGs) were used for GO and KEGG enrichment analysis. The results of GO enrichments indicated that the significantly upregulated genes in the fast-growing group were involved in protein folding, actin cytoskeleton, myosin complex, etc. The results of KEGG enrichments indicated that the significantly upregulated genes in the fast-growing group were involved in glycolysis/ gluconeogenesis, adipocytokine signaling pathway, MAPK signaling pathway, carbon metabolism, PI3K-Akt signaling pathway, etc. To analysis the difference gene sets between fast- and slow-growing group, the RNA-seq data were used by gene set enrichment analysis (GSEA). The results showed that the PI3K/AKT/mTOR pathway was up-regulated in the fast-growing group. The up-regulated of this pathway could lead to higher nutrient absorption efficiency and lead to muscle growth in the fast-growing group. These results contribute to understanding of the molecular mechanisms of fast growth and regulative pathways regulating growth in brown-marbled grouper.

Research Highlights: This study identified the cell cycle genes in birch that likely play important roles during plant growth and development. This analysis provides a basis for understanding the regulatory mechanism of various cell cycles in Betula pendula. Background and Objectives: The cell cycle factors not only influence cell cycle progression together, but also regulate accretion, division and differentiation of cells, and then regulate growth and development of plant. In this study, we identified the putative cell cycle genes in B. pendula genome, based on the annotated cell cycle genes in A. thaliana. It could serve as a foundation for further functional studies. Materials and Methods: The transcript abundance was determined for all the cell cycle genes in xylem, root, leaf and flower tissues using RNA-seq technology. Results: We identified 59 cell cycle gene models in the genome of B. pendula, 17 highly expression genes among them. These genes were BpCDKA.1, BpCDKB1.1, BpCDKB2.1, BpCKS1.2, BpCYCB1.1, BpCYCB1.2, BpCYCB2.1, BpCYCD3.1, BpCYCD3.5, BpDEL1, BpDpa2, BpE2Fa, BpE2Fb, BpKRP1, BpKRP2, BpRb1 and BpWEE1. Conclusions: We identified 17 core cell cycle genes in the genome of birch by combining phylogenetic analysis and tissue specific expression data.

RNA sequencing (RNA-Seq) is a complicated protocol, both in the laboratory in generation of data and at the computer in analysis of results. Several decisions during RNA-Seq library construction have important implications for analysis, most notably strandedness during complementary DNA (cDNA) library construction. Here we clarify bioinformatic decisions related to strandedness in both alignment of DNA sequencing reads to reference genomes and subsequent determination of transcript abundance.

Background: Telomere length (TL) shortening process is associated with several known environment and individual determinants. DNA methylation is the most studied epigenetic process and may be associated with TL. We investigated the associations between DNA methylation and TL in peripheral blood. Methods: Methylation wide association study was conducted in 47 women (37.1±8.8 years) with different nutritional status. Association between TL and DNA methylation levels were explored by univariate and multiple linear regression models, corrected by age and Body Mass Index. Corrections for multiple comparisons by Benjamini‐Hochberg test was also performed. WEBGestalt was used to identify pathways that are responsible for regulating TL. Results: We found negative correlations between TL and BMI (r = -0.641; p = 0.001), abdominal circumference (r = -0.622; p = 0.001) and fat mass (r = -0.656; p = 0.001). 44 CpGs sites were associated with TL, independent of age and BMI. The most of these sites were negatively correlated with TL. For the 7 remained sites, DNA hypomethylation were associated with shorter TL. These CpGs were related to nine different pathways, including thermogenesis, cancer, glutamatergic and serotonergic synapse. Conclusion: There is an epigenetic contribution in TL, independent of nutritional status and age. Genes related to TL are involved in important metabolic pathways.

Replication protein A (RPA) is a heterotrimeric protein complex and the main single-stranded DNA (ssDNA) binding protein in eukaryotes. RPA has key functions in most of the DNA-associated metabolic pathways and DNA damage signalling. Its high affinity for ssDNA helps to stabilise ssDNA structures and prevents the DNA sequence from nuclease attacks. RPA consists of multiple DNA-binding domains which are oligonucleotide/oligosaccharide-binding (OB)-folds, responsible for DNA binding, and interactions with proteins. These RPA-ssDNA and RPA-protein interactions are crucial for DNA replication, DNA repair, DNA damage signalling, and the conversation of the genetic information of cells. Proteins such as ATR use RPA to locate to regions of DNA damage for DNA damage signalling. Recruitment of nucleases and DNA exchange factors to sites of double strand breaks are also an important RPA function to ensure effective DNA recombination to correct these DNA lesions. Due to its high affinity to ssDNA, RPA’s removal from ssDNA is of central importance to allow these metabolic pathway to proceed and processes to exchange RPA against downstream factors are established in all eukaryotes. These facetted and multi-layered functions of RPA will be discussed in detail in the review. RPA is also a major player in a variety of human diseases, which will be discussed.

The Jeryak, a hybrid offspring of Gannan yak and Jersey cattle, exhibits evident hybrid advantages over the Gannan yak in terms of production performance and other various aspects. The small non-coding RNAs known as miRNAs exert a significant regulatory influence on gene expression post-transcriptionally. However, the regulatory mechanism of miRNA associated with muscle development in Jeryak remains elusive. To elucidate the regulatory role of miRNAs in orches-trating the skeletal muscle development in Jeryak, we selected dorsal longissimus muscle tissues from Gannan yak and Jeryak for transcriptome sequencing analysis. A total of 230 (DE)miRNAs and 1819 differentially expressed mRNAs (DEMs) were identified in the longest dorsal muscle of Gannan yak and Jeryak. The functional enrichment analysis revealed a significantly enrichment of target genes from differentially expressed (DE)miRNAs and DEMs in signaling pathways associ-ated with muscle growth, such as the Ras signaling pathway and the MAPK signaling pathway. The network of interactions between miRNA and mRNA suggested that some differentially ex-pressed miRNAs, including miR-2478-z, miR-339-x, novel-m0036-3p and novel-m0037-3p, played a pivotal role in facilitating muscle development. These findings help us to deepen our under-standing of the hybrid dominance of Jeryaks and provide a theoretical basis for further research on the regulatory mechanisms of miRNAs associated with Jeryak muscle growth and development.

The A. japonicus industry has expanded significantly, but no research has focused on how to de-termine the age of A. japonicus during farming. Correctly estimating the age of A. japonicus can provide a decision-making basis for the breeding process, and data for the protection of A. japonicus aquatic germplasm resources. DNA methylation levels in the body wall of Apostichopus japonicus at 4 months, 1 year, 2 years, and 3 years old were determined by MethylRAD-Seq, and differentially methylated genes related to age were screened. The results of the study found that 441 and 966 differentially methylated genes were detected at CCGG and CCWGG sites, respectively. As-partate aminotransferase, succinate semialdehyde dehydrogenase, isocitrate dehydrogenase, the histone H2AX, heat shock protein Hsp90, aminopeptidase N, cell division cycle CDC6, Ras GTPase activating protein (RasGAP), slit guidance ligand slit 1, integrin linked kinase ILK, mechanistic target of rapamycin kinase Mtor, protein kinase A Pka, and autophagy-related 3 atg3 these genes may play key roles in the growth and aging process of A. japonicus. This study provided data for identifying the age of A. japonicus.

The formation and development of potato tissues and organs is a complex process regulated by a variety of genes and environmental factors. However, the regulatory mechanisms underlying the growth and development are still unclear. In this study, we used autotetraploid potato JC14 as experimental subject to analyze the transcriptome of root, stem and leaf at seedling, tuber formation and tuber expansion stages to explore the spatio-temporal expression pattern of genes and genetic development characteristics. The results identified thousands of differentially expressed genes and KEGG pathway enrichment analysis showed that these genes were mainly involved in defense response and carbohydrate metabolism pathways. A total of 12 co-expressed Gene modules were identified by Weighted Gene Co-expression Network Analysis (WGCNA), and 4 modules were screened out with the highest correlation with potato stem developmental traits. Core genes in the network were further investigated and functionally annotated by computing the connectivity of genes within the module. The results unveiled number of hub genes in stems at different developmental stages, including carbohydrate metabolism related genes, the defense response related genes, and transcription factors. These findings provide important leads for further understanding of the molecular regulation and genetic mechanisms of potato tissue development.

Iron is the most abundant micronutrient in plant mitochondria and it has a crucial role in biochemical reactions involving electron transfer. It has been described in Oryza sativa that Mitochondrial Iron Transporter (MIT) is an essential gene and that knockdown mutant rice plants have a decreased amount of iron in mitochondria, strongly suggesting that OsMIT is involved in mitochondrial iron uptake. In Arabidopsis thaliana, two genes encode MIT homologues. In this study, we analyzed different AtMIT1 and AtMIT2 mutant alleles, confirming that individually AtMIT1 nor AtMIT2 genes are essential. When we generated crosses between Atmit1 and Atmit2 alleles we were able to isolate homozygous double mutant plants. Interestingly, homozygous double mutant plants were obtained only when mutant alleles of Atmit2 with the T-DNA insertion in the intron region were used for crossings, and in these cases a correctly spliced AtMIT2 mRNA was generated, although at a low level. Atmit1 Atmit2 double homozygous mutant plants, which were knockout for AtMIT1 and knockdown for AtMIT2, were grown and chacterised in iron sufficient conditions. Pleiotropic developmental defects were observed including abnormal seeds, increased number of cotyledons, slow growth rate, pinoid stems, defects in flower structures and reduced seed set. We observed a possible phenomenon of T-DNA suppression in the next generation of Atmit1 Atmit2 double homozygous mutant plants, correlating with an increased splicing of the AtMIT2 intron containing the T-DNA. Molecular analysis of gene expression markers for mitochondrial and oxidative stress showed that Atmit1 Atmit2 double homozygous mutant plants express a degree of mitochondrial perturbation. A RNA-Seq study was performed and we could identify more than 760 genes differentially expressed in Atmit1 Atmit2, including genes involved in iron transport, coumarin metabolism, and hormones metabolism, transport and signaling. Our data suggest that some of the phenotypes observed in Atmit1 Atmit2 double homozygous mutant plants are mediated by defects in auxin homeostasis.

The mosaic disease in maize is caused by Sugarcane mosaic virus (SCMV), a member of the Potyviridae family. The best strategy to cope with viral infections is the use of disease-resistant maize lines. To better understand the resistance response to SCMV, we analyzed differentially expressed genes among a resistant line (CI-RL1), a susceptible line (B73), and the F1 progeny from a cross between both lines using RNA-Seq data. We also analyzed transcript expression pattern clustering to allocate previously reported resistance candidate genes. GO enrichment analysis of biological processes highlighted a strong regulation in ROS detoxification in both the susceptible and resistant lines. The enrichment of cellular components led to the identification of an integral component of the plasma membrane in the RL line. Transcript expression patterns provide evidence of the importance of host translation in virus response, showing the diverse and complex behavior of eIF4E homologs and the presence of eleven eEF1α factors in maize. In addition, we identified two genes putatively implied in long-distance movement: ZmPiezo and ZmPVIP1. Finally, we propose an ABC transporter to be associated with viral resistance.

The major threats to the sustainable supply of forest tree products are adverse climate, pests and diseases. Climate change, exemplified by increased drought, poses a unique threat to global forest health. This is attributed to the unpredictable behavior of forest pathosystems, which can favor fungal pathogens over the host under persistent drought stress conditions in the future. Currently, the effects of drought on tree resistance against pathogens are hypothetical, thus research is needed to identify these correlations. Norway spruce (Picea abies) is one of the most economically important tree species in Europe, and is considered highly vulnerable to changes in climate. Dedicated experiments to investigate how disturbances will affect the Norway spruce - Heterobasidion sp. pathosystem are important, in order to develop different strategies to limit the spread of H. annosum s.l. under the predicted climate change. Here, we report a transcriptional study to compare Norway spruce gene expressions to evaluate the effects of water availability and the infection of Heterobasidion parviporum. We performed inoculation studies of three-year-old saplings in a greenhouse (purchased from a nursery). Norway spruce saplings were treated in either high (+) or low (-) water groups: high water group received double the water amount than the low water group. RNA was extracted and sequenced. Similarly, we quantified gene expression levels of candidate genes in biotic stress and jasmonic acid (JA) signaling pathways using qRT-PCR, through which we discovered a unique preferential defense response of H. parviporum-infected Norway spruce under drought stress at the molecular level. Disturbances related to water availability, especially low water conditions can have negative effects on the tree host and benefit the infection ability of the pathogens in the host. From our RNA-seq analysis, 114 differentially expressed gene regions were identified between high (+) and low (-) water groups under pathogen attack. None of these gene pathways were identified to be differentially expressed from both non-treated and mock-control treatments between high (+) and low (-) water groups. Finally, only four genes were found to be associated with drought in all treatments.

Long non-coding RNA (lncRNA) can regulate several aspects of gene expression, being associated with complex phenotypes in humans and livestock species. In taurine beef cattle, recent evidence points to the involvement of lncRNA in feed efficiency (FE), a proxy for increased productivity and sustainability. Here, we hypothesized specific regulatory roles of lncRNA in FE of indicine cattle. Using RNA-Seq data from liver, muscle, hypothalamus, pituitary and adrenal gland from Nellore bulls with divergent FE, we submitted new transcripts to a series of filters to confidently predict lncRNA. Then, we identified lncRNA that were differentially expressed (DE) and/or key regulators of FE. Finally, we explored lncRNA genomic location and interactions with miRNA and mRNA to infer potential function. We were able to identify 126 relevant lncRNA for FE in Bos indicus, some with high homology to previously identified lncRNA in Bos taurus and some possible specific regulators of FE in indicine cattle. Moreover, lncRNA identified here were linked to previously described mechanisms related to FE in hypothalamus-pituitary-adrenal axis and are expected to help elucidate this complex phenotype. This study contributes to expanding the catalogue of lncRNA, particularly in indicine cattle, and identifies candidates for further studies in animal selection and management.

In December 2019, a novel coronavirus (SARS-CoV-2) was identified in patients with pneumonia (called COVID-19) in Wuhan, Hubei Province, China. SARS-CoV-2 shares high sequence similarity and uses the same cell entry receptor, angiotensin-converting enzyme 2 (ACE2), as does severe acute respiratory syndrome coronavirus (SARS-CoV). Several studies have provided bioinformatic evidence of potential routes for SARS-CoV-2 infection in respiratory, cardiovascular, digestive and urinary systems. However, whether the reproductive system is a potential target of SARS-CoV-2 infection has not been determined. Here, we investigate the expression pattern of ACE2 in adult human testis at the level of single-cell transcriptomes. The results indicate that ACE2 is predominantly enriched in spermatogonia, Leydig and Sertoli cells. Gene ontology analyses indicate that GO categories associated with viral reproduction and transmission are highly enriched in ACE2-positive spermatogonia while male gamete generation related terms are down-regulated. Cell-cell junction and immunity related GO terms are increased in ACE2-positive Leydig and Sertoli cells, but mitochondria and reproduction related GO terms are decreased. These findings provide evidence that human testes are a potential target of SARS-CoV-2 infection which may have significant impact on our understanding of the pathophysiology of this rapidly spreading disease.

Long non-coding RNAs (lncRNAs) have been found to be involved in many biological processes, including the regulation of cell differentiation, but a complete characterization of lncRNA is still lacking. Additionally, there is evidence that lncRNAs interact with ribosomes, raising questions about their functions in cells. Here, we used a developmentally staged protocol to induce cardiogenic commitment of hESCs and then investigated the differential association of lncRNAs with polysomes. Our results identified lncRNAs in both the ribosome-free and polysome-bound fractions during cardiogenesis and showed a very well-defined temporal lncRNA association with polysomes. Clustering of lncRNAs was performed according to the gene expression patterns during the five timepoints analyzed. In addition, differential lncRNA recruitment to polysomes was observed when comparing the differentially expressed lncRNAs in the ribosome-free and polysome-bound fractions or when calculating the polysome-bound vs ribosome-free ratio. The association of lncRNAs with polysomes could represent an additional cytoplasmic role of lncRNAs, e.g., in translational regulation of mRNA expression.

Peanut (Arachis hypogaea) is a crop plant with high economic value, but the epigenetic regulation of its growth and development has only rarely been studied. The peanut histone deacetylase 1 gene (AhHDA1) has been isolated and is known to be ABA- and drought-responsive. In this paper, we investigate the role of AhHDA1 in more detail, focussing on the effect of altered AhHDA1 expression in hairy roots at both the phenotypic and transcriptional levels. Agrobacterium rhizogenes-mediated transformation of A. hypogaea hairy roots was used to analyse how overexpression or RNA interference of AhHDA1 affects this tissue. In both types of transgenic hairy root, RNA sequencing was adopted to identify genes that were differentially expressed, and these genes were assigned to specific metabolic pathways. AhHDA1-overexpressing hairy roots were growth-retarded after 20 d in vitro cultivation, and superoxide anions and hydrogen peroxide accumulated to a greater extent than in control or RNAi groups. Overexpression of AhHDA1 is likely to accelerate flux through various secondary synthetic metabolic pathways in hairy roots, as well as reduce photosynthesis and oxidative phosphorylation. Genes encoding the critical enzymes caffeoyl-CoA O-methyltransferase (Araip.XGB85) and caffeic acid 3-O-methyltransferase (Araip.Z3XZX) in the phenylpropanoid biosynthesis pathway, chalcone synthase (Araip.B8TJ0) and polyketide reductase (Araip.MKZ27) in the flavonoid biosynthesis pathway, and hydroxyisoflavanone synthase (Araip.0P3RJ) and isoflavone 2'-hydroxylase (Araip.S5EJ7) in the isoflavonoid biosynthesis pathway were significantly upregulated by AhHDA1 overexpression, while their expression in AhHDA1-RNAi and control hairy roots remained at a lower level or was unchanged. Our results suggest that alteration of secondary metabolism activities is related to overexpression of AhHDA1, which is mainly reflected in phenylpropanoid, flavonoid and flavonoid biosynthesis. Future studies will focus on the function of AhHDA1 interacting proteins and their action on cell growth and stress responses.

Recently, selection in pigs has been focused on improving the lean meat content in carcasses; this focus has been most evident in breeds constituting a paternal component in breeding. Such sire-breeds are used to improve the meat quantity of cross-breed pig lines. However, even in one breed, a significant variation in the meatiness level can be observed. In the present study, the comprehensive analysis of genes and microRNA expression profiles in porcine muscle tissue was applied to identify the genetic background of meat content. The comparison was performed between whole gene expression and miRNA profiles of muscle tissue collected from two sire-line pig breeds (Piertain, Hampshire). The RNA-seq approach allowed the identification of 627 and 416 differentially expressed genes (DEGs) between pig groups differing in terms of loin weight between Pietrain and Hampshire breeds, respectively. The comparison of miRNA profiles showed differential expression of 57 microRNAs for Hampshire and 34 miRNAs for Pietrain pigs. Next, 43 genes and 18 miRNAs were selected as differentially expressed in both breeds and potentially related to muscle development. According to Gene Ontology analysis, identified DEGs and microRNAs were involved in the regulation of the cell cycle, fatty acid biosynthesis and regulation of the actin cytoskeleton. The most deregulated pathways dependent on muscle mass were the Hippo signalling pathway connected with the TGF-beta signalling pathway and controlling organ size via the regulation of ubiquitin-mediated proteolysis, cell proliferation and apoptosis. The identified target genes were also involved in pathways such as the FoxO signalling pathway, signalling pathways regulating pluripotency of stem cells and the PI3K-Akt signalling pathway. The obtained results indicate molecular mechanisms controlling porcine muscle growth and development. Identified genes (SOX2, SIRT1, KLF4, PAX6 and genes belonging to the transforming growth factor beta superfamily) could be considered candidate genes for determining muscle mass in pigs.

Track-structure Monte Carlo simulations are useful tools to evaluate initial DNA damage induced by irradiation. In the previous study, we have developed a Gean4-DNA-based application to estimate the cell surviving fraction of V79 cells after irradiation, bridging the gap between the initial DNA damage and the DNA-rejoining kinetics by means of the two-lesion kinetics (TLK) model. However, since the DNA repair performance depends on cell line, the same model parameters cannot be used for the different cell lines. Thus, we extended the Geant4-DNA application with an updated TLK model for the evaluation of DNA damage repair performance in HSGc-C5 carcinoma cells which are typically used for evaluating proton/carbon radiation treatment effects. For this evaluation, we also performed experimental measurements for cell surviving fractions and DNA-rejoining kinetics of the HSGc-C5s cells. Concerning fast- and slow-DNA rejoining, the TLK model parameters were adequately optimized with the simulated initial DNA damage. Using the optimized TLK model, the Geant4-DNA simulation is now able to predict cell survival and DNA-rejoining kinetics for HSGc-C5s cells.

Cancer-associated fibroblasts (CAF) form the basis of tumor microenvironment and possess immunomodulatory functions by interacting with other cells surrounding tumor, including T lymphocytes, macrophages, dendritic cells and natural killer cells. Ionizing radiation is a broadly-used method in radiotherapy to target tumors. In mammalian cells, ionizing radiation induces various types of DNA damages and DNA damage response. Being unspecific, radiotherapy affects all the cells in tumor microenvironment, including the tumor itself, CAFs and immune cells. CAFs are extremely radio-resistant and do not initiate apoptosis even at high doses of radiation. However, following radiation, CAFs become senescent and produce a distinct combination of immunoregulatory molecules. Radiosensitivity of immune cells varies depending on the cell type due to inefficient DNA repair in, for example, monocytes and granulocytes. In this minireview, we are summarizing recent findings on the interaction between CAF, ionizing radiation and immune cells in the tumor microenvironment.

Covalent DNA protein crosslinks (DPCs) are common lesions that block replication. We examine here the consequence of DPCs on mutagenesis involving replicational template-switch reactions in Escherichia coli. 5-azacytidine (5azaC) is a potent mutagen for template-switching, dependent on DNA cytosine methylase (Dcm), implicating the trapped Dcm-DNA covalent complex as the initiator for mutagenesis. The leading strand of replication is more mutable than the lagging strand, explained by blocks to the replicative helicase and/or fork regression. We find that template-switch mutagenesis induced by 5-azaC does not require DSB repair via RecABCD. The ability to induce the SOS response is anti-mutagenic by an unknown mechanism. Mutants in recB, but not recA, exhibit high constitutive rates of template-switching and we suggest that RecBCD-mediated DNA degradation prevents template-switching associated with fork regression. A mutation in the DnaB fork helicase also promotes high levels of template-switching. We also find that other DPC-inducers, formaldehyde (a non-specific crosslinker) and ciprofloxacin (a topoisomerase II poison) are also strong mutagens for template-switching. Induction of mutations and genetic rearrangements that occur by template-switching may constitute a previously unrecognized component of the genotoxicity and genetic instability promoted by DPCs.

It is well known that X-ray crystallography is based on X-ray diffraction (XRD) by atoms and molecules. The diffraction pattern arises as a result of scattering of incident radiation, which makes it possible to determine the structure of the scattering substance. With the advent of ultrashort radiation sources, the theory and interpretation of X-ray diffraction analysis remained the same. This work shows that when an attosecond laser pulse is scattered on a DNA molecule, including during its nicking and bending, the pulse duration is an important characteristic of the scattering. In this case, the diffraction pattern changes significantly compared to the previously known scattering theory. The results obtained must be used in XRD theory to study DNA structures, their mutations and damage, since the previously known theory can produce large errors and, therefore, the DNA structure can be “deciphered” incorrectly.

DNA methylation is an epigenetic mechanism for gene expression modulation and can be used as a predictor of future disease risks. A prospective birth cohort study was performed to clarify the effects of neurotoxicants on child development, namely, the Tohoku Study of Child Development, in Japan. This study aimed to evaluate the association of prenatal exposure to five toxic metals—arsenic, cadmium, mercury, lead (Pb), antimony (Sb), and polychlorinated biphenyls (PCBs, N = 166)—with global DNA methylation in umbilical cord blood DNA. DNA methylation markers, 5-methyl-2’-deoxycytidine (mC) and 5-hydroxymethyl-2’-deoxycytidine (hmC), were determined using liquid chromatography-tandem mass spectrometry. The mC content in cord blood DNA was positively correlated with Pb and Sb levels (r = 0.442 and 0.288, respectively) but not with cord blood PCBs. We also observed significant positive correlations among Pb levels, maternal age, and hmC content (r = 0.159 and 0.243, respectively). The multiple regression analysis among the potential predictors demonstrated consistent positive associations between Pb and Sb levels and mC and hmC content. Our results suggest that global DNA methylation is a promising biomarker for prenatal exposure to Pb and Sb.

COVID-19 was an inflammatory pneumonia caused by a respiratory infection with a coronavirus named SARS-CoV-2. Severely ill COVID-19 patients could die of acute pulmonary and systemic inflammation. Thus, there was a theory that long COVID-19 shared many similarities with systemic autoimmune diseases. There was also a theory that the SARS-CoV-2 virus could infect many cell- and tissue-types distributed throughout multiple organ systems of the body. I aim to review the above theories using unbiased high-throughput datasets cataloguing gene expression. Studies reviewed used the technology of next generation sequencing, either as bulk tissue RNA sequencing, or as single-cell RNAs sequencing; ingenious methods for signal deconvolution were then used to identify individual cell types. Datasets reviewed suggested that severe COVID-19 induced expression of genes associated with pro-inflammatory signaling, tissue development or remodeling, scarring, or cell death. The datasets also provided evidence for several inflammatory syndromes associated with COVID-19, such as neuronal COVID, acute respiratory disease syndrome, vascular inflammation, or multisystem inflammatory syndrome. Indeed, inflammatory signal during severe COVID-19 was reported in many cell types, for example in structural cells, in infiltrating immune cells, in the vascular system, or in circulating cells in the blood. There was also frequently a change in proportions of immune cell-types in the blood, epithelial cells in the lungs, or among infiltrating immune cells.

Abstract5′,8-Cyclo-2′-deoxyadenosine (cdA), in the 5′R and 5′Sdiastereomeric forms, are typical non strand-break oxidative DNA lesions, induced by hydroxyl radicals, with emerging importance as a molecular marker. These lesions are exclusively repaired by nucleotide excision repair (NER) mechanism with a low efficiency, thus readily accumulating in the genome. Poly(ADP-ribose) polymerase1 (PARP1) acts as an early responder to DNA damage and plays a key role as a nick sensor in the maintenance of the integrity of the genome by recognizing nicked DNA. So far, it was unknown whether the diastereomeric cdA lesions could induce specific PARP1 binding. Here we provide the first evidence of PARP1 to selectively recognize the diastereomeric lesions 5′S-cdA and 5′R-cdA in vitro as compared to deoxyadenosine in model DNA substrates (23-mers) by using circular dichroism,fluorescence spectroscopy, immunoblotting analysis and gel mobility shift assay. Several features of the recognition of the damaged and undamaged oligonucleotides by PARP1were characterized. Remarkably, PARP1 efficiently binds to both cdA lesions in the double stranded (ds)-oligonucleotides. In particular, PARP1 proved to bind 5′S-cdAwith a higher affinity constant for the 5'S lesion in a model of ds DNA than 5′R-cdA, showing different recognition patterns, also compared with undamaged dA. This new finding highlights the ability of PARP1 to recognize and differentiate the distorted DNA backbone in a biomimetic system caused by different diastereomeric forms of a cdA lesion.

Recent advances in DNA sequencing technology have shown that the human milk microbiota of healthy women varies substantially. The gDNA extraction method may influence the observed variation, biasing the microbiological reconstruction after all. In this study, a genomic DNA extraction method for DNA isolation from human milk samples was standardized and compared with commercial and standard hose make methods. Spectrophotometric measurements, gel electrophoresis, and PCR amplifications were used as criteria for evaluating the quantity, quality, and functionality of the extracted DNA. Furthermore, the standardized method of extracting gDNA from human milk was evaluated for its ability to isolate functional DNA from gram-positive, and gram-negative bacteria and fungi, to improve the reconstruction of microbiological profiles. The novel DNA extraction method increased the quantity and quality of the gDNA extracted compared with commercial and standard house-make protocols. This method even allowed PCR amplification of the V3-V4 regions of the 16S ribosomal gene in all samples, and the ITS-1 region of the fungal 18S ribosomal gene in 95 % of the samples as well. It is concluded that the proposed method provides better performance for the extraction of gDNA from complex samples such as human milk.

DNA methylation, one of the major epigenetic mechanisms, plays critical roles in regulating gene expression, genomic stability and cell lineage commitment. Establishment and maintenance of DNA methylation in mammals is achieved by two groups of DNA methyltransferases: DNMT3A and DNMT3B, which are responsible for installing DNA methylation patterns during gametogenesis and early embryogenesis, and DNMT1, which is essential for propagating DNA methylation patterns during replication. Both groups of DNMTs are multi-modular proteins, containing a large N-terminal regulatory region in addition to the C-terminal methyltransferase domain. Recent structure-function investigations of the individual domains or large fragments of DNMT1 and DNMT3A have revealed the molecular basis for their substrate recognition and specificity, intramolecular domain-domain interactions, as well as their crosstalk with other epigenetic mechanisms. These studies highlight a multifaceted regulation for both DNMT1 and DNMT3A/3B, which is essential for the precise establishment and maintenance of lineage-specific DNA methylation patterns in cells. This review summarizes current understanding of the structure and mechanism of DNMT1- and DNMT3A-mediated DNA methylation, with emphasis on the functional cooperation between the methyltransferase and regulatory domains.

Background: We investigated the tumor immune response in gastric cancer patients receiving third line nivolumab monotherapy to identify immune-related biomarkers for better patient se-lection. Methods: Nineteen patients (10 males, median age 67 years) who received nivolumab as a third or later line therapy were enrolled. We analyzed the tumor immune response in durable clinical benefit (DCB) and non-DCB patients. Pre-treatment and early-on-treatment tumor transcriptomes were examined, and gene expression profiles, immunograms, and T cell receptor (TCR) repertoire were analyzed. Results: DCB was observed in 15.8% of patients, with comparable secondary endpoints (ORR, OS, PFS) to previous trials. Individual immunograms showed no significant changes before and ear-ly-on-treatment, nor consistent alterations among DCB cases. The intratumoral immune response was suppressed by previous treatments in most third or later line nivolumab recipients. TCR repertoire analysis revealed newly emerged clonotypes in on-treatment tumors, but clonal re-placement did not impact efficacy. High T cells/Tregs ratios and a low UV radiation response gene signature were linked to DCB and treatment response. Conclusions: This study emphasizes the tumor immune response's importance in nivolumab ef-ficacy for gastric cancer. High T cells/Tregs ratios and specific gene expression signatures show promise as potential biomarkers for treatment response. Larger cohort validation is crucial to optimize immune checkpoint inhibitors in gastric cancer treatment.

Single-cell RNA sequencing (scRNA-seq) technology has significantly advanced our understanding of the diversity of cells and how this diversity is implicated in diseases. Yet, translating these findings across various scRNA-seq datasets poses challenges due to technical variability and dataset-specific biases. To overcome this, we present a novel approach that employs both an LLM-based framework and explainable machine learning to facilitate generalization across single-cell datasets and identify gene signatures to capture disease-driven transcriptional changes. Our approach uses scBERT, which harnesses shared transcriptomic features among cell types to establish consistent cell-type annotations across multiple scRNA-seq datasets. Additionally, we employ a symbolic regression algorithm to pinpoint highly relevant yet minimally redundant models and features for inferring a cell type’s disease state based on its transcriptomic profile. We ascertain the versatility of these cell-specific gene signatures across datasets, showcasing their resilience as molecular markers to pinpoint and characterize disease-associated cell types. Validation is carried out using four publicly available scRNA-seq datasets from both healthy individuals and those suffering from ulcerative colitis (UC). This demonstrates our approach’s efficacy in bridging disparities specific to different datasets, fostering comparative analyses. Notably, the simplicity and symbolic nature of the retrieved gene signatures facilitate their interpretability, allowing us to elucidate underlying molecular disease mechanisms using these models.

The achievement of reproductive competence in male mammals is dependent on the testis. Goat testis’ development and spermatogenesis involve physiological events with high complexity. In the current work, 6 testes were respectively collected from immature, sexually mature and physically mature Qianbei Ma goats (1, 6 and 12 months old, respectively). RNA-Seq was carried out to reveal changes in testis mRNA expression levels in Qianbei Ma goats at various developmental stages, and gene expression profiling at different ages was established. Totally 18 libraries were established for screening genes and pathways associated with testis development and spermatogenic processes. Totally 9,724 upregulated and 4,153 downregulated genes were identified between immature (I) and sexually mature (S) testes; 7 upregulated and 3 downregulated genes were detected between sexually mature (S) and physically mature (P) testes, and approximately 4% of genes were alternately spliced between the I and S groups. Selected genes were verified by qRT-PCR, in agreement with sequencing data indicating their reliability. Those genes have critical functions in various developmental stages of goat testicular development and spermatogenesis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were carried out to evaluate differentially expressed genes (DEGs). GO analysis suggested DEGs were involved in “reproduction process”, “channel activity” and “cell periphery part” between I and S, and in “ion transport process”, “channel activity” and “transporter complex part” between S and P. KEGG analysis indicated that pathways such as “glycerolipid metabolism”, “steroid hormone biosynthesis” and “MAPK signaling pathway” may be involved in testis development and spermatogenesis. Genes including IGF1, TGFB1, TGFBR1 and EGFR may regulate the development of the testis from immature to sexually mature, which may be key candidate genes for the development of goat testis. These findings provide novel insights into goat testicular development and spermatogenesis.

Rose (Rosa spp.), especially R. hybrida, is one of the most popular ornamental plants in the world and the third largest cut flower crop in Taiwan. Rose mosaic disease (RMD), showing mosaic, line patterns and ringspots on leaves, is a common rose disease caused by complex infection of various viruses. Due to pests and diseases, the rose planting area in Taiwan has been decreasing since 2008; however, no rose virus disease has been reported in the past five decades. In the spring of 2020, rose samples showing RMD-like symptoms were observed at an organic farm in Chiayi, central Taiwan. The virome in the farm was analyzed by RNA-seq. Rose genomic sequences were filtered from the obtained reads. The remaining reads were de novo assembled to generate 294 contigs, 50 of which were annotated as viral sequences corresponding to 10 viruses. Through reverse transcription-polymerase chain reaction validation, a total of seven viruses were detected, including six known rose viruses, namely apple mosaic virus, prunus necrotic ringspot virus, rose partitivirus, apple stem grooving virus, rose spring dwarf-associated virus and rose cryptic virus 1, and a novel ilarvirus. After completing the whole genome sequencing and sequence analysis, the unknown ilarvirus was demonstrated as a new species, tentatively named rose ilarvirus 2. This is the first report of the rose viruses in Taiwan.

The objective of this study was to identify key transcription factors involved in lipid metabolism and immune response related to the differentially expressed genes (DEG) from the liver samples of 35 pig model for metabolic diseases fed diets containing either 1.5 or 3.0% soybean oil (SOY1.5 or SOY3.0). A total of 281 DEG between SOY1.5 and SOY3.0 diets (log2fold-change ≥ 1 or ≤ −1; FDR-corrected p-value < 0.1) were identified, in which 129 were down-regulated and 152 were up-regulated in SOY1.5 group. The functional annotation analysis detected transcription factors linked to lipid homeostasis and immune response, such as RXRA, EGFR, and SREBP2 precursor. These findings demonstrated that key transcription factors related to lipid metabolism could be modulated by dietary inclusion of soybean oil. It could contribute to nutrigenomics research field that aims to elucidate dietary interventions in animal and human health, as well as to drive the food technology and science.

Background: Jerusalem artichoke (Helianthus tuberosus L.) is tolerant to salinity stress and has high economic value. The salt tolerance mechanisms of Jerusalem artichoke are still unclear. Especially in the early stage of Jerusalem artichoke exposure to salt stress, the plant physiology, biochemistry and gene transcription are likely to undergo large changes. Elucidating these changes may be of great significance to understanding the salt tolerance mechanisms of it. Results: We obtained high-quality transcriptome from leaves and roots of Jerusalem artichoke exposed to salinity (300 mM NaCl) for 0 h, 6 h, 12 h, 24 h and 48 h, with 150,129 unigenes and 9023 DEGs (Differentially Expressed Genes). The RNA-seq data were clustered into time-dependent groups (nine clusters each in leaves and roots); gene functions were distributed evenly among the groups convergence. KEGG enrichment analysis showed the genes related to plant hormone signal transduction were enriched in almost all treatment comparisons. Under salt stress, genes belongs to PYL (abscisic acid receptor PYR / PYL family), PP2C (Type 2C protein phosphatases), GH3 (Gretchen Hagen3), ETR (ethylene receptor), EIN2/3 (ethylene-insensitive protein 2/3), JAZ (Genes such as jasmonate ZIM-domain gene) and MYC2 (Transcription factor MYC2) had extremely similar expression patterns. The results of qPCR of 12 randomly selected genes confirmed the accuracy of RNA-seq. Conclusions: Under the impact of high salinity (300mM) environment, Jerusalem artichoke in the seedling stage was difficult to survive for a long time, and the phenotype was severe in the short term. Based on the expression of genes on the time scale, we found that the distribution of gene functions in time is relatively even. Upregulation of the phytohormone signal transduction had a crucial role in the response of Jerusalem artichoke seedlings to salt stress, the genes of abscisic acid, auxin, ethylene, and jasmonic acid had the most obvious change pattern.

A definitive transcriptome atlas for the non-coding expressed elements of pathogenic mycobacteria does not exist. Incomplete lists of non-coding transcripts can be obtained for some of the reference genomes (e.g. Mycobacterium tuberculosis H37Rv) but to what extent these transcripts have homologues in closely related species or even strains is not clear. This has implications for the analysis of transcriptomic data; non-coding parts of the transcriptome are often ignored in the absence of formal, reliable annotation. Here, we review the state of our knowledge of non-coding RNAs in pathogenic mycobacteria, emphasising the disparities in the information included in commonly used databases. We then proceed to review ways of combining computational solutions for predicting the non- coding transcriptome with experiments that can help refine and confirm these predictions.

Single cell RNA-seq (scRNA-seq) profiles conceal temporal and spatial tissue developmental information. De novo reconstruction of single cell temporal trajectory has been fairly addressed, but reverse engineering single cell 3D spatial tissue localization is hitherto landmark based, and de novo spatial reconstruction is a compelling computational open problem. Here we show that a new algorithm - named D-CE - for coalescent embedding of single cell transcriptomic networks can address this open problem. We rely merely on the spatial information encoded in the expression patterns of developmental signal transcription factor (DST) genes, and we find that D-CE of cell-cell association DST-transcriptomic networks reliably reconstructs the Geo-seq or single cell samples’ 3D spatial tissue distribution. Comparison to the novoSpaRC and CSOmap (recent and only available de novo 3D spatial reconstruction methods) on 16 datasets and 681 reconstructions, reveals a significantly distinctive superior performance of D-CE.

Background: DSCAM-AS1 is a cancer-related long noncoding RNA with higher expression levels in Luminal A, B and HER2-positive Breast Cancer (BC), where its expression is strongly dependent on Estrogen Receptor Alpha (ERα). Methods: To decipher its function, DSCAM-AS1 expression was measured by qRT-PCR in tissue samples from 93 BC patients in addition to a meta-analysis of 30 gene expression datasets, together with the evaluation of its association with clinical data. By computational analyses of our RNA-Seq in MCF-7 cells, we investigated the DSCAM-AS1 knock-down effects at both gene and isoform levels. Results: We confirmed DSCAM-AS1 overexpression in high grade Luminal A, B and HER2+ BCs and found a significant correlation with disease relapse. 908 genes were regulated by DSCAM-AS1-silencing, primarily involved in cell cycle and inflammatory response. Noteworthy, the analysis of alternative splicing and isoform regulation revealed 2,085 splicing events regulated by DSCAM-AS1, enriched in differential polyadenylation sites and 3’UTR shortening events. Finally, the DSCAM-AS1-interacting splicing factor hnRNPL was predicted as the most enriched RBP for exon skipping and 3’UTR events. Conclusion: The relevance of DSCAM-AS1 overexpression in BC is confirmed by clinical data and further enhanced by its possible involvement in the regulation of RNA processing, which is emerging as one of the most important dysfunctions in cancer.

Powdery mildew caused by Blumeria graminis f. sp. tritici (Bgt) is one of many severe diseases that threaten bread wheat (Triticum aestivum L.) yield and quality worldwide. The discovery and deployment of powdery mildew resistance genes (Pm) can prevent this disease epidemic in wheat. In a previous study, we transferred the powdery mildew resistance gene Pm57 from Aegilops searsii into common wheat and cytogenetically mapped the gene in a chromosome region with the fraction length (FL) 0.75-0.87, which represents 12% of 2Ss#1 segment on the long arm of chromosome 2Ss#1. In this study, we performed RNA-Seq on three infected and mock-infected wheat-Ae. searsii 2Ss#1 introgression lines with Bgt-isolates inoculation at 0, 12, 24, and 48 hours after inoculation. Then we designed 79 molecular markers based on transcriptome sequences and physically mapped them to Ae. searsii chromosome 2Ss#1- in seven intervals. We used these markers to identify 46 wheat-Ae. searsii 2Ss#1 recombinants induced by ph1b, a deletion mutant of pairing homoelogous (Ph) genes. Analysis of the 46 ph1b-induced 2Ss#1L recombinants with different Bgt-responses using 28 2Ss#1L-specific molecular markers in the interval FL0.72-0.87 where Pm57 is located, and the flanking intervals, we physically mapped Pm57 gene on the long arm of 2Ss#1 in a 5.13 Mb genomic region, which was flanked by markers X67593 (773.72 Mb) and X62492 (778.85 Mb). By comparative synteny analysis of the corresponding region on chromosome 2B in Chinese spring (T. aestivum L.) with other model species we identified ten genes that are putative plant defense-related (R) genes which includes six coiled-coil nucleotide-binding site-leucine-rich repeat (CNL), three nucleotide-binding site-leucine-rich repeat (NL) and a leucine-rich receptor-like repeat (RLP) encoding proteins. This study will lay a foundation for further cloning of Pm57, and benefit the understanding of interactions between resistance genes of wheat and powdery mildew pathogens.

The term ‘Hallmarks of Cancer’ was coined by Hanahan and Weinberg in their influential reviews and they described genome instability as a property of cells enabling cancer development [1, 2]. Accurate DNA replication of genomes is central to diminish genome instability. Here, the understanding of the initiation of DNA synthesis in origins of DNA replication to start leading strand synthesis and the initiation of Okazaki fragment on the lagging strand are crucial to control genome instability. Recent findings have provided new insights into the mechanism of the remodelling of the prime initiation enzyme, DNA polymerase α-primase, during primer synthesis, how the enzyme complex achieves lagging strand synthesis, and how it is linked to replication forks to achieve optimal initiation of Okazaki fragments. Moreover, the central roles of RNA primer synthesis by Pol-prim in multiple genome stability pathways such as replication fork restart and protection of DNA against degradation by exonucleases during double-strand break repair is discussed.

Cells respond to DNA double-strand breaks by initiating DSB repair and ensuing a cell cycle checkpoint. The primary responder to DSB repair is non-homologous end joining, an error prone repair pathway. However, when DSBs are generated after DNA replication in the G2 phase of the cell cycle, a second DSB repair pathway, homologous recombination, can come into action. Both ATM and ATR are important for DSB-induced DSB repair and checkpoint responses. One method of ATM and ATR working together is through the DNA end-resection of DSBs. As a readout and marker of DNA end-resection, RPA is phosphorylated at Ser4/Ser8 of the N-terminus of RPA32 in response to DSBs. Here, the significance of RPA32 Ser4/Ser8 phosphorylation in response to DNA damage, specifically in S phase to G2 phase of the cell cycle is examined. RPA32 Ser4/Ser8 phosphorylation in G2 synchronized cells is necessary for increases in TopBP1 and Rad9 accumulation on chromatin and full activation of the ATR-dependent G2 checkpoint. In addition, our data suggest RPA Ser4/Ser8 phosphorylation modulates ATM-dependent KAP-1 phosphorylation and Rad51 chromatin loading in G2 cells.

DNA damage repair lies at the core of all cells’ survival strategy, including cancerous. Therefore, targeting such repair mechanisms forms the major goal of cancer therapeutics. The mechanism of DNA repair has been tousled with the discovery of multiple kinases. Recent studies on Tousled like Kinases have brought significant clarity on the effectors of these kinases which stands to regulate DSB repair. In addition to their well-established role in the DDR and cell cycle checkpoint mediation after DNA damage or inhibitors of replication, their suspected involvement in the actual DSB repair process has more recently been strengthened by the important finding that TLK1 phosphorylates RAD54 and regulates some of its activities and localization in the cell. Earlier findings of its regulation of RAD9 during checkpoint deactivation as well as defined steps during NHEJ ends processing were earlier hints of its important involvement broadly in DSB repair. All this has opened up new avenues to target cancer cells in combination therapy with genotoxins and TLK inhibitors.

Head and neck squamous cell carcinomas (HNSCCs) develop from the mucosa of the oral cavity, pharynx, and larynx. The methylation levels of Septin 9 (SEPT9) and short stature homeobox 2 (SHOX2) genes in ccfDNA are considered pan-cancer biomarkers and have shown prognostic value in preliminary reports in HNSCC. Liquid biopsy is a non-invasive procedure that collects tumor-derived molecules including circulating cell-free DNA (ccfDNA). Here, we developed a ddPCR-based assay to detect the DNA methylation levels of plasma circulating SEPT9 and SHOX2 in patients with HNSCC. We first set up the assay on commercial methylated and un-methylated DNA. Then, the dynamic changes of the methylation levels of SEPT9 and SHOX2 were quantified in 20 patients during the follow-up. The results highlighted: i) the capability of the ddPCR-based assay in detecting very low copies of methylated molecules; ii) the significant decrease in methylation levels of SEPT9 and SHOX2 in plasma of HNSCC patients at the first time points of follow-up respect to T0; iii) a different trend of longitudinally DNA methylation variations associable to the clinicopathological features of the patients. The absolute quantification of the methylation levels of SEPT9 and SHOX2 in HNSCC may be used for risk stratification and disease monitoring.

Forensic genetic genealogy, a technique leveraging new DNA capabilities and public genetic databases to identify suspects, raises specific considerations in a law enforcement context. Use of this technique requires consideration of its scientific and technical limitations, including the composition of current online datasets, and consideration of its scientific validity. Additionally, forensic genetic genealogy needs to be considered in the relevant legal context to determine the best way in which to make use of its potential to generate investigative leads while minimising its impact on individual privacy. This article presents these issues from an Australian perspective, with the observations and conclusions likely to be applicable to other jurisdictions.

One of creators of quantum mechanics P. Jordan in his work on quantum biology claimed that life's missing laws were the rules of chance and probability of the quantum world. The article presents author’s results of studying frequencies (or probabilities) of nucleotides on so-called epi-chains of long DNA sequences of various eukaryotic and prokaryotic genomes. DNA epi-chains are algorithmically constructed subsequencies of DNA nucleotide sequences. According to the algorithm of construction of any epi-chain of the order n, the epi-chain is such nucleotide subsequence, in which the numerations of adjacent nucleotides differ by natural number n (n = 1, 2, 3, 4,…). Correspondingly each epi-chain of order n ≥ 2 contains n times less nucleotides than the original DNA sequence. The presented results unexpectedly discover that in long single-stranded and double-stranded DNA of any tested genome its DNA epi-chains of different orders n (values n are not too large) have practically identical frequencies (or probabilities) of each kind of nucleotides. These data allow considering DNA as a regular rich set of epi-chains, which can play a certain role in genetic and epigenetic phenomena as the author belives. Appropriate rules of nucleotide frequencies on epi-chains of long DNA sequences are formulated for further their tests on a wider set of genomes. These results testify on existence of long-range coherence in long DNA and remind the Fröhlich's theory of long-range coherence in biological systems. The phenomenological data are discussed from different standpoints: the DNA double helices and helical antennas with circular polarizations of electromagnetic waves; relations with the Fröhlich's theory; numerical analysis of DNA epi-chains under binary representations of nucleotides. Results are useful for developing quantum and algebraic biology.

DNA double-strand breakage is the most lethal damage to chromosomal DNA. It activates a series of cellular DNA damage response pathways, including DNA damage sensing, control of cell cycle arrest and apoptosis, and DNA repair. DNA damage response pathways are regulated by complex signaling machineries. Of the intracellular signaling cascades, diacylglycerol kinase (DGK) phosphorylates diacylglycerol (DG) to generate phosphatidic acid (PA). Because both DG and PA serve as second messengers, DGK activity induces a shift of signaling pathways from DG-mediated to PA-mediated cascades, thereby implicating DGK in the regulation of widely various functions. Reportedly, one member of the DGK family, DGKζ, is intimately involved in the regulation of stress responses through p53 and NF-κB. Stresses such as ischemia and infarction cause DGKζ downregulation. Experimental DGKζ depletion renders cells and mice vulnerable to various stressors such as chemotherapeutic agents and ionizing irradiation. Nevertheless, how DGKζ is involved in DNA repair, a critical event of DNA damage response for survival remains unknown. For this study, we examined how DGKζ depletion affects DNA repair mechanisms. We demonstrated that DGKζ depletion causes attenuation of Akt activation and DNA-PK protein expression upon DNA damage, which might engender downregulated BRCA1 protein synthesis and stability. Results suggest that DGKζ depletion attenuates BRCA1-mediated DNA repair machinery, thereby conferring vulnerability to DNA damage.

Trichothecene mycotoxin T-2 toxin is an important environmental pollutant and poses a global threat to human and animal health. T-2 toxin could induce nephrotoxicity, however, the precise molecular mechanism remains unclear. In this study, the mice were intraperitoneally administrated with a single dose of 2 mg/kg T-2 toxin. The kidney function and ultrastructural observation were also assessed and found T-2 toxin exposure could cause kidney damage. Transmission electron microscopy showed a significant swelling and vacuolization of the mitochondria of renal cell. A total of 1122, 58, and 391 differential expression genes (DEGs) that were predominantly expressed in kidney tissues at 1 d, 3 d, and 7 d after T-2 toxin exposure, respectively. Early transcriptomic changes at 1 d after exposure were down-regulation of DEGs involved in the cell cycle, p53 signaling pathway, and cellular senescence, while up-regulation of DEGs referred to the ribosome pathway. The temporal variations in gene expression pattern in T-2 toxin exposure in mice kidney were presented and cellular metabolism was disturbed at the exposure recovery period with 7 d. In conclusion, this study, for the first time, provided a comprehensive comparative transcriptomic analysis of gene regulation in T-2 toxin exposure-induced nephrotoxicity at different temporal periods and explored the nephrotoxicity mechanism of T-2 toxin at the mRNA level.

The shoot apical meristem (SAM) gives rise to the aerial structure of plants by producing lateral organs and secondary meristems. The SAM is responsible for plant developmental patterns, thus determining plant morphology and, consequently, many agronomic traits such as the number and size of fruits and flowers and kernel yield. Our current understanding of SAM morphology and regulation is based on studies conducted mainly on some angiosperms, including economically important crops such as maize (Zea mays) and rice (Oryza sativa), and the model species Arabidopsis (Arabidopsis thaliana). However, studies in other plant species from the gymnosperm class are scant, making difficult comparative analyses that help us understand SAM regulation in diverse plant species. This limitation prevents deciphering the mechanisms by which evolution gave rise to the multiple plant structures within the plant kingdom, and determined the conserved mechanisms involved in SAM maintenance and operation. This review aims to integrate and analyze the current knowledge of SAM evolution by combining the morphological and molecular information recently reported from the plant kingdom.

Because rice is native to tropical and subtropical regions, it is generally more sensitive to cold stress compared to other cereals. Low temperature stress has become one of the most important conditions that affects the growth, development, and yield of rice. Plant NAC (NAM, ATAF1/2, and CUC) transcription factors are known to play significant regulatory roles in the stress response. In our research, we found that OsNAC050 (LOC_Os03g60080） is mainly expressed in leaves, and low temperature can further up-regulate its expression. OsNAC050 mutants created using CRISPR-Cas9 gene editing technology showed significantly enhanced tolerance to low temperature treatment. Detection of enzyme activities related to the redox pathway also showed that the mutants had stronger viability under low temperature stress. Comparative transcriptome analysis showed that photosynthesis and soluble sugar metabolism were significantly affected in the osnac050 mutant lines, suggesting that OsNAC050 may participate in the above molecular pathways in the response to low temperature stress. The results will enhance our understanding of molecular mechanisms underlying the responses to cold stress in rice and can provide new strategies for engineering cold-tolerance in high-yielding rice varieties.

The emergence of drug resistant pathogenic bacteria is increasingly challenging conventional antibiotics. Plant derived flavonoids are considered as potential alternatives to antibiotics due to their antimicrobial properties. However, the mechanisms by which flavonoids modulate pathogenic microorganisms’ growth are not fully understood. In our previous studies we found that Rutin, a kind of flavonoids showed inhibition against Klebsiella pneumoniae strains. In order to better understand its inhibitory mechanism on strain growth, we used both Rutin and Luteolin as treatments to incubate K. pneumoniae ATCC700603 strain. After incubation for 4 hours, the homogenous and differential effects of two different flavonoids on the growth of K. pneumoniae were evaluated. The integrated metabolomic and transcriptomic analysis was performed. The expression levels of 5,483 genes and the contents of 882 metabolites were measured. The differentially expressed genes (DEGs) and metabolites (DEMs) were screened and analyzed. The correlation between DEGs and DEMs were also studied. Our results showed that Rutin could inhibit the strain growth by changing metabolic pathways and ABC transporters pathways. Our study also revealed FU841_RS17580 and FU841_RS19145 as functional genes that played vital roles in the strain growth.

Systems genetics is key for integrating a large number of variants associated with diseases. Vitamin K (VK) is one of the scarcely studied conditions in lieu of ascertaining either the differentially expressed genes (DEGs) or variants in an individual subpopulation of diseased phenotypes associated with VK, viz. myocardial infarction, renal failure, prostate cancer, thrombosis, thrombocytopenia, coagulation related diseases to name a few. In this work, we have screened characteristic DEGs common to three VK-related diseases, viz. myocardial infarction, renal failure and prostate cancer and asked whether or not any DEGs in addition to pathogenic variants are common to these conditions. We attempt to bridge the gap in finding characteristic biomarkers and discuss the role of long noncoding RNAs (lncRNAs) in the biogenesis of VK deficiencies.

Single-cell RNA sequencing data facilitates investigation of cell heterogeneity and subpopulations as well as differentially abundant states however modern single-cell RNA sequencing datasets are growing in size and complexity requiring advances in the bioinformatic methods that analyze them. Many methods exist for each step of analysis including read alignment, normalization, quality control, batch effect correction, imputation and dimensionality reduction. With so many options to choose from at each step of the analysis, benchmarking and a synthesis of the literature on the methods available is necessary to inform biological researchers on the most optimal workflow for their data. Here, recent key methods of analysis are highlighted with a focus on methods that facilitate identification of cell subpopulations and differentially abundant cell states. With a constantly expanding toolset for each step in single-cell RNA sequencing dataset analysis, biological researchers should stay informed to utilize the most applicable methods for their own analyses.

Microglia activation causes neuroinflammation, which is a hallmark of neurodegenerative disorders, brain injury, and aging. Ladostigil, a bifunctional reagent with antioxidant and anti-inflammatory properties, reduced microglial activation and enhanced brain functioning in elderly rats. In this study, we studied SH-SY5Y, a human neuroblastoma cell line, and tested viability in the presence of hydrogen peroxide and Sin1 (3-morpholinosydnonimine), which generates reactive oxygen and nitrogen species (ROS/RNS). Both stressors caused significant apoptosis and necrotic cell death that was attenuated by ladostigil. Our results from RNA-seq experiments show that long non-coding RNAs (lncRNAs) account for 30% of all transcripts in SH-SY5Y cells treated with Sin1 for 24 hours. Altogether, we identify 94 differently expressed lncRNAs in the presence of Sin1, including MALAT1, a highly expressed lncRNA with anti-inflammatory and anti-apoptotic functions. Additional activities of Sin-1 upregulated lncRNAs include redox homeostasis (e.g., MIAT, GABPB1-AS1), energy metabolism (HAND2-AS1), and neurodegeneration (e.g., MIAT, GABPB1-AS1, NEAT1). Four lncRNAs implicated as enhancers were significantly upregulated in cells exposed to Sin1 and ladostigil. Finally, we show that H2O2 and Sin1 increased the expression of DJ-1, a redox sensor and modulator of Nrf2 (nuclear factor erythroid 2–related factor 2). Nrf2 (NFE2L2 gene) is a major transcription factor regulating antioxidant genes. In the presence of ladostigil, DJ-1 expression is restored to its baseline. The mechanisms governing SH-SY5Y cell survival and homeostasis are highlighted by the beneficial role of ladostigil in the crosstalk involving Nrf2, antioxidant transcription factor DJ-1, and lncRNAs. Stress-dependent induction of lncRNAs represents an underappreciated regulatory level that contributes to cellular homeostasis and the capacity of SH-SY5Y to cope with oxidative stress.

Intratumoral heterogeneity associates with more aggressive disease progression and worse patient outcomes. Understanding the reasons enabling the emergence of such heterogeneity remains incomplete, which restricts our ability to manage it from a therapeutic perspective. Technological advancements such as high-throughput molecular imaging, single-cell omics and spatial transcriptomics now allow recording the patterns of spatiotemporal heterogeneity in a longitudinal manner, thus offering insights into the multi-scale dynamics of its evolution. Here, we review latest technological trends and biological insights from molecular diagnostics as well as spatial transcriptomics, both of which have witnessed a burgeoning growth in recent past in terms of mapping heterogeneity within tumor cell types as well as stromal constitution. We also discuss ongoing challenges, indicating possible ways to integrate insights across these methods to have a systems-level spatiotemporal map of heterogeneity in each tumor, and a more systematic investigation of implications of heterogeneity for the patient outcomes.

Prostate cancer (Pca) is a highly heterogeneous disease and the second more common tumor in males. Molecular and genetic profiles have been used to identify subtypes and guide therapeutic intervention. However, roughly 26% of primary Pca are driven by unknown molecular lesions. We use Principal Component Analysis (PCA) and custom RNAseq-data normalization to identify a gene expression signature which segregates primary PRAD from normal tissues. This Core-Expression Signature (PRAD-CES) includes 33 genes and accounts for 39% of data complexity along the PC1-cancer axis. The PRAD-CES is populated by protein-coding (AMACR, TP63, HPN) and RNA-genes (PCA3, ARLN1) sparsely found in previous studies, validated/predicted biomarkers (HOXC6, TDRD1, DLX1), and/or cancer drivers (PCA3, ARLN1, PCAT-14). Of note, the PRAD-CES also comprises six over-expressed LncRNAs without previous Pca association, four of them potentially modulating driver’s genes TMPRSS2, PRUNE2 and AMACR. Overall, our PCA capture 57% of data complexity within PC1-3. GO enrichment and correlation analysis involving major clinical features (i.e., Gleason Score, AR Score, TMPRSS2-ERG fusion and Tumor Cellularity) suggest that PC2 and PC3 gene signatures might describe more aggressive and inflammation-prone transitional forms of PRAD. Of note, surfaced genes may entail novel prognostic biomarkers and molecular alterations to intervene. Particularly, our work uncovered RNA genes with appealing implications on Pca biology and progression.

Whole genome pooled sequence data of 12 Pakistani Teddy goats is analyzed for positive selection signatures as their breed defining characteristics. Selection imprints left in the Teddy genome are unveiled by genomic differentiation after the successful paired-end alignment of 635,357,043 reads with (ARS1) reference genome assembly. Pooled-heterozygosity ( ) and Tajima’s D (TD) are applied for validation and getting better hits of selection signals, while pairwise F_ST statistics is conducted on Teddy vs. Bezoar (wild goat ancestor) for genomic differentiation. Annotation of regions under positive selection reveals 59 genes underlying production and adaptive traits. score ≥ 5 detected six windows having highest scores on Chr. 29, 9, 25, 15 and 14 that harbor HRASLS5, LACE1 and AXIN1 genes which are candidate for embryonic development, lactation and body height. Secondly, TD value of ≤ -2.2 showed 4 windows with very strong hits on Chr.5 & 9 harbor STIM1 and ADM genes related to body mass and weight. Lastly, F_ST analysis generated three strong signals with threshold ≤ 0.42 on Chr.12 & 5 harbor ITGB1 gene associated with milk production & lactation traits. Other significant selection signatures encompass genes associated with wool production, prolificacy, immunity and coat colors. In brief, this study identified the genes under selection in this Pakistani goat breed that will be helpful to refining future breeding policies and converging required productive traits within and across other goat breeds and to explore full genetic potential of this valued livestock species.

As a major infectious disease in chickens, Newcastle disease causes considerable economic losses in the poultry industry, especially in developing countries where there is limited access to effective vaccination. Therefore, enhancing resistance to the virus in commercial chickens through breeding is a promising way to promote poultry production. In this study, we investigated gene expression changes at 2 and 6 dpi after infection at day21 with a lentogenic Newcastle disease virus in a commercial egg-laying chicken hybrid using RNA sequencing analysis. By comparing NDV challenged and nonchallenged groups, 526 differentially expressed genes (DEGs) (FDR < 0.05) were identified at 2 dpi, and only 36 at 6 dpi. For the DEGs at 2 dpi, IPA analysis predicted inhibition of multiple signaling pathways in response to NDV that regulate immune cell development and activity, neurogenesis and angiogenesis. Upregulation of Interferon Induced Protein with Tetratricopeptide Repeats 5 (IFIT5) in response to NDV was consistent between the current and most previous studies. Sprouty RTK Signaling Antagonist 1 (SPRY1), a DEG in the current study is located in a significant QTL associated with virus load at 6 dpi in the same population. These identified pathways and DEGs provide potential targets to further study breeding strategy to enhance NDV resistance in chickens.

An outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) occurred in China towards the end of 2019, and has spread rapidly ever since. Previous studies showed that some virus could affect the reproductive system and cause long-term complications. Recent studies exploring the source of SARS-CoV-2 using genomic sequencing have revealed that SARS-CoV-2 enters the host cells via the angiotensin-converting enzyme II (ACE2), the receptor that recognizes SARS-CoV. To investigate the expression of ACE2 and to explore the potential risk of infection in the reproductive system, we performed a thorough bioinformatic analysis on data from public databases involving RNA expression, protein expression, and single-cell RNA expression studies. The analyzed data showed high levels of ACE2 mRNA and protein expression in the testis and spermatids and equal levels of ACE2 expression in the uterus and lung. Comprehensive single-cell analysis identified ACE2 expression in the lung, testis, spermatids, and uterus. In conclusion, this study revealed the potential risk associated with the SARS-CoV-2 infection in the reproductive system and predicted that long-term complications might have a significant impact on the prevention and management of COVID-19, the disease caused upon infection with SARS-CoV-2.

Platinum(II) complexes have been found to be effective against cancer cells. Cisplatin curbs cell replication by interacting with the deoxyribonucleic acid (DNA), eventually leading to cell death and reducing cell proliferation. In order to investigate the ability of platinum complexes to affect cancer cells, two examples from the class of polyflurophenylorganoamidoplatinum(II) complexes were synthesised and tested on isolated DNA. The two compounds trans-[N,N’-bis(1,2,3,5,6-pentafluorophenyl)ethane-1,2-diaminato(1-)](2,3,4,5,6-pentafluorobenzoato)(pyridine)platinum(II) (PFB), and trans-[N,N’-bis(1,2,3,5,6-pentafluorophenyl)ethane-1,2-diaminato(1-)](2,4,6-trimethylbenzoato)(pyridine)platinum(II) (TMB) were compared with cisplatin through their reaction with DNA. Attenuated Total Reflection Fourier Transform Infrared (ATR-FTIR) spectroscopy was applied to analyse the interaction of the Pt^(II) complexes with DNA in the hydrated, dehydrated and rehydrated state. These were compared with control DNA in acetone/water (PFB, TMB) and isotonic saline (cisplatin) under the same conditions. Principle Component Analysis (PCA) was applied to compare the ATR-FTIR spectra of the untreated control DNA with spectra of PFB and TMB treated DNA samples. Disruptions in the conformation of DNA treated with the Pt^(II) complexes upon rehydration were mainly observed by monitoring the position of the IR-band around 1711 cm^-1 assigned to the DNA base-stacking vibration. Furthermore, other intensity changes in the phosphodiester bands of DNA at ~1234 cm^-1 and 1225 cm^-1 and shifts in the dianionic phosphodiester vibration at 966 cm^-1 were observed. The isolated double stranded DNA (dsDNA) or single stranded DNA (ssDNA) showed different structural changes when incubated with the studied compounds. PCA confirmed PFB had the most dramatic effect by denaturing both dsDNA and ssDNA. Both compounds, along with cisplatin, induced changes in DNA bands at 1711, 1088, 1051 and 966 cm^-1 indicative of DNA conformation changes. The ability to monitor conformational change with infrared spectroscopy paves the way for a sensor to screen for new anticancer therapeutic agents.

Four 1,10-phenanthroline derivatives (1-4) have been synthesized as potential telomeric DNA binders, three substituted in their chains with thiosemicarbazones (TSCs), the neutral compound (1) and two ionic derivatives (PF6¯ and BF4¯ salts, 2 and 3, respectively), and a fourth corresponding to a 4-phenylthiazole derivative (4). Quadruplex and dsDNA interactions have been preliminarily studied, especially for the neutral derivative 1, by FRET-based DNA melting assays, equilibrium dialysis (both competitive and non-competitive), circular dichroism and viscosity titrations. The TSC derivatives bind and stabilize the telomeric Tel22 quadruplex more efficiently than dsDNA, with an estimated 24-fold selectivity determined by equilibrium dialysis for compound 1. In addition, cytotoxic activity against different tumor cells (PC-3, DU145, HeLa, MCF-7 and HT29), and two non-tumor cell lines (HFF-1 and RWPE-1) has been evaluated. Except for the 4-phenylthiazole derivative, which was inactive, the compounds show moderate cytotoxic properties, with the salts displaying lower IC50 values, compared to the neutral TSC. However, the neutral derivative was the only compound that exhibited a modest selectivity in the case of prostate cells (tumor PC-3 versus healthy RWPE-1). Cell cycle analysis and Annexin V/PI assays revealed that the compounds can produce cell death by apoptosis, an effect that has proven similar to that demonstrated by other known 1,10-phenanthroline G4 ligands endowed with antitumor properties, such as PhenDC3 or PhenQE8.

B chromosomes are supernumerary chromosomes which are found in addition to the normal standard chromosomes (A chromosomes). B chromosomes are well known to accumulate several types of repeats. Although the evolution of B chromosomes has been subject of numerous studies, the mechanisms of accumulation and evolution of repetitive sequences is not fully understood. Recently, new genomic approaches have shed light on the origin and accumulation of different classes of repetitive sequences in the process of B chromosome formation and evolution. Here we discuss the impact of repetitive sequences accumulation on the evolution of plant B chromosomes.

Cholangiocarcinoma is an epithelial malignancy arising in the region between the intrahepatic bile ducts and the ampulla of Vater at the distal end of the common bile duct. The effect of current chemotherapy regimens against cholangiocarcinoma is limited, and the prognosis of patients with cholangiocarcinoma is poor. Aberrant DNA methylation and histone modification induce silencing of tumor suppressor genes and chromosomal instability during carcinogenesis. Studies have shown that the tumor suppressor genes and microRNAs (miRNAs) including MLH1, p14, p16, DAPK, miR-370 and miR-376c are frequently methylated in cholangiocarcinoma. Silencing of these tumor suppressor genes and miRNAs plays critical roles in the initiation and progression of cholangiocarcinoma. In addition, recent studies have demonstrated that DNA methylation inhibitors induce expression of endogenous retroviruses and exert the anti-tumor effect of via an anti-viral immune response. Aberrant DNA methylation of tumor suppressor genes and miRNAs could be a powerful biomarker for diagnosis and treatment of cholangiocarcinoma. Epigenetic therapy with DNA methylation inhibitors hold considerable promise for the treatment of cholangiocarcinoma through re-activation of tumor suppressor genes and miRNAs as well as induction of an anti-viral immune response.

Astronauts in space are subject to a continuous bombardment of ionizing radiation. The Earth's magnetic field and the ISS shield some biologically damaging particles traveling through. Still, travel beyond low-Earth orbit and extravehicular activities are exceedingly more dangerous and there is a concern for the acute and late-occurring adverse health effects befalling astronauts. So, it is vital to consider the current tools and models used to describe and study the organic consequences of ionizing radiation exposure. It is equally important to see where these models may improve. This article reviews the historical development and current state of knowledge of radiation effects impacting astronauts in orbit. We explain the space radiation environment, cellular microenvironment, and how these may be incorporated in radiobiological models to aid in our understanding of the influence space travel may have on astronaut health. The topics discussed in this paper include a review of DNA damage and repair mechanisms and the numerical models that aim to explain the biological effects resulting from ionizing radiation damage. Historically, radiobiological models focused on how radiation damages nuclear DNA, built upon the hypotheses of Crowther and Lea in the 1940s and 1960s, and neglected other sub-cellular targets outside of nuclear DNA.

Environmental DNA (eDNA) sampling has emerged as a powerful approach for monitoring biodiversity in the context of diverse terrestrial and aquatic ecosystems, and has broad potential for both environmental surveys and species management. eDNA studies focusing on plant and animal species have capitalized on analytical methodologies developed originally for metagenomic studies of microorganisms, and on technological improvements in sample collection. However, while eDNA studies have proliferated rapidly, efforts to standardize sampling and analytical methods are still nascent. Critically, fundamental data concerning the distribution and persistence of eDNA in different ecological contexts are lacking, complicating both confidence in species detection and comparability between studies. We propose that key deficits could be addressed by using well-designed synthetic DNA standards or tracers, which are widely used in fields ranging from medical diagnostics to hydrology. In-field application of a standardized mixture of barcoded synthetic DNA fragments of varying lengths and concentrations that are recovered concommitantly with eDNA sampling should provide critical yet missing data on environmental DNA distribution and turnover, as well as an internal molecular standard to enhance both species detection and comparability between eDNA studies.

Methylation of cytosines at CpG sites in DNA, carried out by de novo by DNA methyltransferase Dnmt3a, is a basic epigenetic modification involved in gene regulation and genome stability. Aberrant CpG methylation in gene promoters leads to oncogenesis. In oncogene promoters, CpG sites often colocalize with guanine-rich sequences capable of folding into G-quadruplexes (G4s). Our in vitro study aimed to investigate how parallel G4s formed by a sequence derived from the c-MYC oncogene promoter region affect the activity of the Dnmt3a catalytic domain (Dnmt3a-CD). For this purpose, we designed synthetic oligonucleotide constructs: c-MYC G4-forming oligonucleotide and linear double-stranded DNA containing an embedded stable extrahelical c-MYC G4. The topology and thermal stability of G4 structures in these DNA models were analyzed by physicochemical techniques. We showed that Dnmt3a-CD specifically binds to an oligonucleotide containing c-MYC G4, resulting in inhibition of its methylation activity. c-MYC G4 formation in a double-stranded context significantly reduce Dnmt3a-CD-induced methylation of a CpG site located in close proximity to the quadruplex structure; this effect depends on the distance between the non-canonical structure and the specific CpG site. One would expect DNA hypomethylation near the G4 structure, while regions distant from this non-canonical form would maintain a regular pattern of high methylation levels. We hypothesize that the G4 structure sequesters the Dnmt3a-CD and impedes its proper binding to B-DNA, resulting in hypomethylation and activation of c-MYC transcription.

A large body of evidence indicates that environmental agents can induce alterations in DNA methylation (DNAm) profiles. Radiofrequency electromagnetic fields (RF-EMFs) are radiations emitted by everyday devices which have been classified as “possibly carcinogenic”, however their biological effects are unclear. As aberrant DNAm of genomic repetitive elements (RE) may promote genomic instability, here we sought to determine whether exposure to RF-EMFs could affect DNAm of different classes of RE, such as long interspersed nuclear element1s (LINE-1), Alu short interspersed nuclear element and ribosomal repeats. To this purpose, we analysed DNAm profiles of cervical cancer and neuroblastoma cell lines (HeLa, BE(2)C and SH-SY5Y) exposed to 900 MHz GSM-modulated RF-EMF through an Illumina-based targeted deep bisulfite sequencing approach. Our findings showed that radiofrequency exposure did not affect DNAm of Alu elements in any of the cell lines analysed. Conversely, it influenced DNAm of LINE1 and ribosomal repeats in terms of both average profiles and organization of methylated and un-methylated CpG sites, in different ways in each of the three cell lines studied.

Here, our data provide the first evidence for the existence of a previously unknown receptive system formed by novel DNA- and RNA-based receptors in eukaryotes. This system, named the TR-system, is capable of recognizing and generating a response to different environmental factors and has been shown to orchestrate major vital functions of fungi, mammalian cells, and plants.Recently, we discovered the existence of a similar regulatory system in prokaryotes. These DNA- and RNA-based receptors are localized outside of the membrane forming a type of a network around cells that respond to a variety of chemical, biological, and physical factors and enabled the TR-system to regulate major aspects of eukaryotic cell life as follows: growth, including reproduction and development of multicellular structures; sensitivity to temperature, geomagnetic field, UV, light, and hormones; interaction with viruses; gene expression, recognition and utilization of nutrients. The TR-system was also implicated in cell memory formation and was determined to be responsible for its maintenance and the forgetting of preceding events. This system is the most distant receptive and regulatory system of the cell that regulates interactions with the outer environment and governs the functions of other receptor-mediated signaling pathways.

Here we present 1. a model for amplicon sequencing 2. a definition of the best assignment of a read to a set of reference sequences 3. strategies and structures for indexing reference sequences and computing the best assign- ments of a set of reads efficiently, based on (ultra)metric spaces and their geometry The models, techniques, and ideas are particularly relevant to scenarios akin to 16S taxonomic profiling, where both the number of reference sequences and the read diversity is considerable.

Structural DNA nanotechnology provides an excellent foundation for diverse nanoscale shapes that can be used in various bioapplications and materials research. From all existing DNA assembly techniques, DNA origami has proven to be the most robust one for creating custom nanoshapes. Since its invention in 2006, building from the bottom up using DNA has drastically advanced, and therefore, more and more complex DNA-based systems have become accessible. So far, vast majority of the demonstrated DNA origami frameworks are static by nature, but interestingly, there also exist dynamic DNA origami devices that are increasingly coming into view. In this review, we discuss DNA origami nanostructures that perform controlled translational or rotational movement triggered by predefined DNA strands, various molecular interactions and/or other external stimuli such as light, pH, temperature and electromagnetic fields. The rapid evolution of such dynamic DNA origami tools will undoubtedly have a significant impact on molecular scale precision measurements, targeted drug delivery and diagnostics, but they can also play a role in development of optical/plasmonic sensors, nanophotonic devices and nanorobotics for numerous different tasks.

We present a numerical study in which large-scale bulk simulations of self-assembled DNA constructs have been carried out with a realistic coarse-grained model. The investigation aims at obtaining a precise, albeit numerically demanding, estimate of the free energy for such systems. We then, in turn, use these accurate results to validate a recently proposed theoretical approach that builds on a liquid-state theory, the Wertheim theory, to compute the phase diagram of all-DNA fluids. This hybrid theoretical/numerical approach, based on the lowest order virial expansion and a nearest-neighbor DNA model, can provide, in an undemanding way, a thermodynamic description of DNA associating fluids that is in semi-quantitative agreement with experiments. We show that the predictions of such scheme are as accurate as the ones obtained with more sophisticated methods. We also demonstrate the flexibility of the approach by incorporating non-trivial additional contributions that go beyond the nearest-neighbor model to compute the DNA hybridization free energy.

Many fundamental molecular techniques (PCR, Microarray, Southern and northern hybridization, siRNA, CRISPR/Cas9 etc.) developed so far shows errors. I wish to highlight these molecular techniques are developed on basis of Watson-Crick DNA model, ignoring the concept of parallel stranded DNA. Through this opinion article, I wish to highlight specificity and accuracy of these molecular techniques can be enhanced by considering both parallel and anti parallel hybridization of DNA. Hopefully my views will also solve issue of irreproducibility in life science research.

The interaction between native and invasive plant species plays a crucial role in understanding invasion mechanisms and developing control strategies. This study investigated the differences in photosynthetic and fluorescence parameters, gene expression profiles, and rhizosphere soil microbial composition between the native species Ficus tikoua Bur. and the invasive species Alternanthera philoxeroides. For the first time, it was demonstrated that that F. tikoua Bur. may enhance photoinhibition ability and the expression of genes involved in carbon metabolism, protein processing in the endoplasmic reticulum, and plant-pathogen interactions to maintain growth. Additionally, it recruits beneficial microbes, e.g. Actinomarina (Bacteria), Nitrosotenuis (Archaea) and Laccaria (fungi), to suppress pathogens, thereby preserving soil DMSP and carbon cycling pathways to resist invasion by A. philoxeroides. These findings provide novel insights into the mechanisms underlying plant invasion and suggest potential strategies for invasive plant management. Further research is needed to investigate this relationship and determine the mechanisms underlying the potential resistance of F. tikoua Bur. to plant invasion.

The corneal epithelium (CE) is spread between two domains, the outer vascularized limbus and the avascular cornea proper. Epithelial cells undergo constant migration from the limbus to the vision-critical central cornea. Coordinated with this migration the cells undergo differentiation changes where a pool of unique stem/precursor cells at the limbus yields the mature cells that reach the corneal center. Differentiation is heralded by the expression of the corneal-specific Krt12. Processing data acquired by scRNA-Seq we found that the increase of Krt12 expression occurs in four distinct steps within the limbus plus a single continuous increase in the cornea. Differential gene analysis demonstrated that these domains reflect discreet stages of CE differentiation and yielded extensive information of the genes undergoing down or up regulation in the sequential transition from less to more differentiate conditions. The approach allowed the identification of multiple gene cohorts, including a) the genes that have maximal expression in the most primitive, Krt12-negative cell cohort, which is likely to include the stem/precursor cells; b) the sets of genes that undergo continuous increase or decrease along the whole differentiation path, and c) and the genes showing maximal positive or negative correlation with the changes in Krt12.

Intramuscular fat content (IFC) and meat color are vital indicators affecting pork quality. A significant positive correlation between IFC and redness meat color (a* value) indicates that these two traits are likely to be regulated by shared molecular pathways. To identify candidate genes, hub genes, and signaling pathways that regulate these two traits, we measured the IFC and a* value in 147 hybrid pigs and selected individuls with extreme phenotypes for transcriptome analysis. The results revealed 485 and 394 overlapping differentially expressed genes (DEGs) using DESeq2, limma, and edgeR packages affecting IFC and a* value, respectively. Weighted gene co-expression network analysis (WGCNA) identified four modules significantly correlated with IFC and a* value. Moreover, we integrated functional enrichment analysis results based on DEGs, GSEA, and WGCNA conditions to identify candidate genes, and identified 47 and 53 candidate genes affecting IFC and a* value, respectively. The protein-protein interaction (PPI) network analysis of candidate genes showed that 6 and 14 hub genes affecting IFC and a* value, respectively. These genes mainly participated in various pathways related to lipid metabolism and redox reactions. Notably, four crucial hub genes (MYC, SOX9, CEBPB, and PPAGRC1A) were shared for these two traits. After functional annotation of these four hub genes, we hypothesized that the SOX9/CEBPB/PPARGC1A axis could co-regulate lipid metabolism and myoglobin redox response. Further research on these hub genes, especially the SOX9/CEBPB/PPARGC1A axis, will help to understand the molecular mechanism of the co-regulation of IFC and a* value, which will provide the theoretical basis for improving pork quality.

Injury to the central nervous system (CNS), in most vertebrate animals, results in permanent damage and lack of function, due to their limited regenerative capacities. In contrast, echinoderms can fully regenerate their radial nerve cord (RNC) following transection, with little or no scarring. Investigators have associated the regenerative capacity of some organisms with the stress response and inflammation produced by the injury. Here we explore the gene activation profile of the stressed holothurian CNS. To do this, we performed RNA sequencing on isolated RNC explants submitted to the stress of transection and enzyme dissection and compared them to explants kept in culture for 3 days following dissection. We describe stress-associated genes, including members of heat-shock families, ubiquitin-related pathways, transposons, and apoptosis that were differentially expressed. Surprisingly, the stress response does not induce apoptosis in this system. Other genes associated with stress in other animal models, such as hero proteins and those associated with the integrated stress response, were not found to be differentially expressed either. Our results provide a new viewpoint on the stress response in the nervous system of an organism with an amazing regenerative capacity. This is the first step to deciphering the molecular processes that allow echinoderms to undergo fully functional CNS regeneration while also providing a comparative view for students of the stress response in other organisms.

Numerous proteomic and transcriptomic studies have been carried out to better understand the current multi-variant SARS-CoV-2 virus mechanisms of action and effects. However, they are mostly centered on mRNAs and proteins. The effect of the virus on human post-transcriptional regulatory agents such as microRNAs (miRNAs) involved in the regulation of 60% of human gene activity remains poorly explored. Similar to what we have previously done with other viruses such as Ebola and HIV, in this study we investigated the miRNA profile of lung epithelial cells following infection with SARS-CoV-2. At the 24 and 72 hours post-infection, SARS-CoV-2 did not drastically alter the miRNome. About 90% of the miRNAs remained non-differentially expressed. The results revealed that miR-1246, miR-1290 and miR-4728-5p were the most upregulated over time. miR-196b-5p and miR-196a-5p were the most downregulated at 24 h while at 72 h, miR-3924, miR-30e-5p and miR-145-3p showed the highest level of downregulation. In the top significantly enriched KEGG pathways of genes targeted by differentially expressed miRNAs we found, among others, MAPK, RAS, P13K-Akt and renin secretion signaling pathways. By RT-qPCR, we also showed that SARS-CoV-2 may regulate several predicted host mRNA targets involved in the entry of the virus into host cells (ACE2, TMPRSS2, ADAM17 and FURIN), in renin–angiotensin system (RAS) (Renin, Angiotensinogen, ACE), innate immune response (IL-6, IFN1β, CXCL10, SOCS4) and fundamental cellular processes (AKT, NOTCH, WNT). Finally, we demonstrated by dual luciferase assay a direct interaction between miR-1246 and ACE-2 mRNA. This study highlights the modulatory role of miRNAs in the pathogenesis of SARS-CoV-2.

In Escherichia coli, DNA replication termination is orchestrated by two clusters of Ter sites forming a DNA replication fork trap when bound by Tus proteins. The formation of a ‘locked’ Tus-Ter complex is essential for halting incoming DNA replication forks. However, the absence of replication fork arrest at some Ter sites raised questions about their significance. In this study, we examined the genome-wide distribution of Tus and found that only the six innermost Ter sites (TerA-E and G) were significantly bound by Tus. We also found that a single ectopic insertion of TerB in its non-permissive orientation could not be achieved, advocating against a need for ‘back-up’ Ter sites. Finally, examination of the genomes of a variety of Enterobacterales revealed a new replication fork trap architecture mostly found outside the Enterobacteriaceae family. Taken together, our data enabled the delineation of a narrow ancestral Tus-dependent DNA replication fork trap consisting of only two Ter sites.

Inducible genetically defined mouse models of cancer uniquely facilitate the investigation of early events in cancer progression, however there are valid concerns about the ability of such models to faithfully recapitulate human disease. We developed an inducible mouse model of progressive lung adenocarcinoma (LuAd) that combines sporadic activation of oncogenic KRas^G12D with modest overexpression of c-MYC (KM model). Histological examination revealed a highly reproducible transition from adenoma to locally invasive adenocarcinoma within 6 weeks of oncogene activation. Laser-capture microdissection coupled with RNA-SEQ was employed to determine transcriptional changes associated with tumour progression. Upregulated genes were triaged for relevance to human LuAd using datasets from Oncomine and cBioportal. Selected genes were validated by RNAi screening in human lung cancer cell lines and examined for association with lung cancer patient overall survival using KMplot.com. Depletion of progression-associated genes resulted in pronounced viability and/or cell migration defects in human lung cancer cells. Progression-associated genes moreover exhibited strong associations with overall survival, specifically in human lung adenocarcinoma, but not in squamous cell carcinoma. The KM mouse model faithfully recapitulates key molecular events in human lung cancer and is a useful tool for mechanistic interrogation of LuAd progression.

Plants grown under fluctuating light impact plant developments compared with those grown under non-fluctuating light conditions. However, our knowledge on the underlying regulatory mechanisms is still quite limited, particularly from the transcriptional perspective. In order to investigate the influence of different light intensity distributions on tomato plant development, we designed three fluctuating light intensity distributions with the non-fluctuating light intensity as control and compared the transcriptional differences after five weeks of treatment. We found plant height and aerial/root weight were significantly reduced under all fluctuating light treatments. Transcriptome analysis revealed that the number of up and down regulated genes had a distinct distribution pattern between different treatments and control. The largest difference between the numbers of down and up regulated genes was found between treatment 1 and 3, reaching to a total of 416 genes. The number and type of the top 20 enriched pathways differed between treatments and control. The largest number of genes enriched was involved in the biosynthesis of secondary metabolites. These results provide insights into the transcriptional regulations of tomato under different light intensity distributions.

B chromosomes (B) are supernumerary elements found in many taxonomic groups. Most B chromosomes are rich in heterochromatin and composed of abundant repetitive sequences, especially transposable elements (TEs). Bs origin is generally linked to the A chromosome complement (A). The first report of a B chromosome in African cichlids was on Astatotilapia latifasciata, which can harbor 0, 1 or 2 B chromosomes. Classical cytogenetics studies found high TE content on the species B chromosome. In this study, we aim to understand TE composition and expression on A. latifasciata genome and its relation to the B chromosome. We use bioinformatics analysis to explore TEs genome organization and also their composition on the B chromosome. Bioinformatics findings were validated by fluorescent in situ hybridization (FISH) and real-time PCR (qPCR). A. latifasciata has a TE content similar to other cichlid fishes and several expanded elements on its B chromosome. With RNA sequencing data (RNA-seq) we showed that all major TE classes are transcribed in brain, muscle and male/female gonads. The evaluation of TE expression between B- and B+ individuals showed that few elements have differential expression among groups and expanded B elements were not highly transcribed. Putative silencing mechanisms may the acting on the B chromosome of A. latifasciata to prevent adverse consequences of repeat transcription and mobilization in the genome.

Solanesol is a noncyclic terpene alcohol composed of nine isoprene units and it mainly accumulates in solanaceous plants, especially tobacco (Nicotiana tabacum L.). Here, RNA-seq analyses of tobacco leaves, stems, and roots were used to identify solanesol biosynthesis genes. Six 1-deoxy-d-xylulose 5-phosphate synthase, two 1-deoxy-d-xylulose 5-phosphate reductoisomerase, two 2-C-methyl-d-erythritol 4-phosphate cytidylyltransferase, four 4-diphosphocytidyl-2-C-methyl-d-erythritol kinase, two 2-C-methyl-d-erythritol 2,4-cyclodiphosphate synthase, four 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate synthase, two 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase, six isopentenyl diphosphate isomerase, and two solanesyl diphosphate synthase (SPS) genes were identified to be involved in solanesol biosynthesis. Furthermore, the two N. tabacum SPS (NtSPS1 and NtSPS2), which had two conserved aspartate-rich DDxxD domains, were highly homologous with SPS enzymes from other solanaceous plant species. In addition, the solanesol contents of three organs, and leaves from four growing stages, corresponded with the distribution of chlorophyll. Our findings provide a comprehensive evaluation of the correlation between the expression of different biosynthetic genes and the accumulation of solanesol in tobacco.

Abstract Background In tumors, somatic mutagenesis presumably drives DNA damage response (DDR) via altered regulatory pathways increasing genomic instability and proliferative activity. These considerations led to the standard therapeutic strategy against cancer: the disruption of mutation activated DNA repair pathways of tumors. Purpose Justifying that cancer cells are not enemies to be killed, but rather they are ill human cells having the remnants of physiologic regulatory pathways. Results 1. Genomic instability and cancer development may be originated from the flaw of estrogen signal rather than excessive estrogen signaling. 2. Healthy cells with genomic instability exhibit somatic mutations helping DNA restitution. 3. Somatic mutations in tumor cells aim the restoration of DNA damage rather than further genomic derangement. 4. In tumors, estrogen signal drives the pathways of DNA stabilization leading to apoptotic death. 5. In the peritumoral cellular infiltration, the genomic damage of tumor induces inflammatory cytokine secretion and increased estrogen synthesis. In the inflammatory cells, increased growth factor receptor (GFR) signal confers unliganded activation of estrogen receptors (ERs). 6. In tumor cells responsive to genotoxic therapy, constitutive mutations help the upregulation of estrogen signal and consequential apoptosis. In tumors non responsive to genotoxic therapy, the possibilities for ER activation via either liganded or unliganded pathway, are exhausted leading to farther genomic instability and unrestrained proliferation. Conclusion Understanding the real character and behavior of human tumors at molecular level suggests that we should learn the genome repairing methods of tumors and follow them by supportive therapy rather than provoking additional genomic damages.

The National Institute of Justice reported that current methods for processing sexual assault samples have a high failure rate, with 60 to 80 percent of tested kits unable to produce usable DNA profiles. Even when samples test positive for male DNA, 34 percent of sexual assault kits (SAK) do not yield recovered male DNA after differential extraction. Less than 30% recovery of available sperm DNA contributes to this low success rate. The SpermX™ method (SX) has been shown to recover 80 percent or more of sperm DNA from sexual assault samples. An interlaboratory evaluation compared SX to standard differential extraction (DE) protocols. Mock samples with known ratios of female epithelial cells and sperm cells were processed using both methods. Results revealed that SX consistently provided CODIS up-loadable DNA profiles, even with as few as 25 sperm cells, whereas DE failed to produce usable results. On average, SX yielded a seven-fold increase in the recovery of unshared male alleles compared to DE. In conclusion, SX outperformed DE in recovering higher quantities of male DNA with minimal female carryover in sexual assault-type samples. This improved success rate in obtaining usable DNA profiles can significantly aid in solving sexual assault cases. Keywords: Differential extraction, Sexual Assault Kit (SAK), Sperm DNA, STRmix™, Combined DNA Index System (CODIS), Forensic DNA analysis, Low Sperm Samples, DNA Mixtures, Autosomal STR Analysis, SpermX.

Comprehensive studies into the various causes of this ubiquitous disease have been prompted by the worldwide increase in cancer incidence. There are many potential causes of cancer, but environmental factors are of particular concern now. This abstract explores the complex link between environmental factors and the growing cancer epidemic, shedding light on the role that our immediate environment plays in the emergence of cancers. The environment can contain carcinogenic agents that gradually influence the course of cancer, including air pollutants, water contaminants, workplace exposures, and industrial activities. The ways through which environmental carcinogens may cause DNA damage and unchecked cell proliferation are discussed in this abstract. Environmental justice and the moral necessity to reduce inequities in exposure and health effects are highlighted to underline the disproportionate impact on vulnerable groups. The study presented here highlights the need for preemptive action to reduce environmental risks. Society can reduce the hidden threat posed by environmental carcinogens through strict legislation, sustainable practices, and community-driven initiatives. In order to lessen the cancer burden and secure a healthier future for all, the authors of this abstract argue that there must be a concerted effort to create surroundings that promote health and well-being.

Topoisomerases, common targets for anti-tumor therapeutics, are crucial enzymes for DNA replication, transcription, and many other aspects of DNA metabolism. The potential anti-cancer effect of thiosemicarbazones (TSC) and metal-TSC complexes has been demonstrated by targeting several biological processes, including DNA metabolism. Human topoisomerases were discovered among the molecular targets for TSCs, and metal-chelated TSCs specifically displayed significant inhibition. The processes by which metal-TSCs or TSCs inhibit topoisomerases are still being studied. In this brief review, we summarize the TSCs and metal-TSCs that inhibit various types of human topoisomerases.

Phenotypes are necessary for genomic evaluations and management. Sometimes genomics can be used to measure phenotypes when other methods are difficult or expensive. Prolificacy of bulls used in multiple-bull pastures for commercial beef production is an example. A retrospective study of 79 bulls aged 2-year-old and older used 141 times in 4-5 pastures across 4 years was used to estimate repeatability from variance components. Traits available before each season’s use were tested for predictive ability. Sires were matched to calves using individual genotypes and evaluating exclusions. A lower cost method of measuring prolificacy was simulated for 5 pastures using the bulls’ genotypes and pooled genotypes to estimate average allele frequencies of calves and of cows. Repeatability of prolificacy was 0.62 ± 0.09. A combination of age-class and scrotal circumference accounted for less than 5 % of variation. Simulated estimation of prolificacy by pooling DNA of calves was accurate. Adding pooling of cow DNA or actual genotypes both increased accuracy about the same. Knowing a bull’s prior prolificacy would help predict future prolificacy for management purposes and could be used in genomic evaluations and research with coordination of breeders and commercial beef producers.

The ability to predict physical characteristics from DNA presents significant opportunities for forensic science. Giving scientists an ability to make predictions about the donor of genetic material at a crime scene can then give investigators new intelligence leads for cold cases where DNA evidence has not identified any person of interest. However, the interpretation of this new form of intelligence requires careful analysis. The responses to an online survey, conducted in 2018-19, were used to examine how actors in the criminal justice system assess and interpret different types of DNA evidence and intelligence. The groups of focus for the survey were investigators, legal practitioners and the general public (as potential jurors). Several statistically significant effects were identified based on occupation and whether an individual had prior exposure to new DNA technology. Monitoring how those involved in interpreting reports from different types of DNA evidence and intelligence interpret them helps to ensure that decisions are made based on a sound understanding of their capabilities and limitations and may inform broader training and awareness strategies.

A fully analytical treatment of the base-pair and codon dynamics in double-stranded DNA molecules is introduced, by means of a realistic treatment which considers different mass values for G, A, T, and C nucleotides and takes into account the intrinsic three-dimensional, helicoidal geometry of DNA in terms of a Hamitonian in cylindrical coordinates. Within the framework of the Peyrard-Dauxois-Bishop model we consider the coupling between stretching and stacking radial oscillations as well as the twisting motion of each base pair around the helix axis. By comparing the linearized dynamical equations for the angular and radial variables when going from the bp local scale to the longer triplet codon scale, we report an underlying hierarchical symmetry. The existence of synchronized collective oscillations of the base-pairs and their related codon triplet units are disclosed from the study of their coupled dynamical equations. The possible biological role of these correlated, long-range oscillation effects in double standed DNA molecules containing mirror-symmetric codons of the form XXX, XX’X, X’XX’, YXY, and XYX is discussed in terms of the dynamical equations solutions and their related dispersion relations.