Identification of hub genes and potential molecular mechanisms related to drug sensitivity in acute myeloid leukemia based on machine learning.

Background: Acute myeloid leukemia (AML) is the most common form of leukemia among adults and is characterized by uncontrolled proliferation and clonal expansion of hematopoietic cells. There has been a significant improvement in the treatment of younger patients, however, prognosis in the elderly AML patients remains poor. Methods: We used computational methods and machine learning (ML) techniques to identify and explore the differential high-risk genes (DHRGs) in AML. The DHRGs were explored through multiple in silico approaches including genomic and functional analysis, survival analysis, immune infiltration, miRNA co-expression and stemness features analyses to reveal their prognostic importance in AML. Furthermore, using different ML algorithms, prognostic models were constructed and validated using the DHRGs. At the end molecular docking studies were performed to identify potential drug candidates targeting the selected DHRGs. Results: We identified a total of 80 DHRGs by comparing the differentially expressed genes derived between AML patients and normal controls and high-risk AML genes identified by Cox regression. Genetic and epigenetic alteration analyses of the DHRGs revealed a significant association of their copy number variations and methylation status with overall survival (OS) of AML patients. Out of the 137 models constructed using different ML algorithms, the combination of Ridge and plsRcox maintained the highest mean C-index and was used to build the final model. When AML patients were classified into low- and high-risk groups based on DHRGs, the low-risk group had significantly longer OS in the AML training and validation cohorts. Furthermore, immune infiltration, miRNA coexpression, stemness feature and hallmark pathway analyses revealed significant differences in the prognosis of the low- and high-risk AML groups. Drug sensitivity and molecular docking studies revealed top 5 drugs, including carboplatin and austocystin-D that may significantly affect the DHRGs in AML. Conclusion: The findings from the current study identified a set of high-risk genes that may be used as prognostic and therapeutic markers for AML patients. In addition, significant use of the ML algorithms in constructing and validating the prognostic models in AML was demonstrated. Although our study used extensive bioinformatics and machine learning methods to identify the hub genes in AML, their experimental validations using knock-out/-in methods would strengthen our findings.

Predicting Prognosis and Immunotherapy Response in Multiple Cancers Based on the Association of PANoptosis-Related Genes with Tumor Heterogeneity.

PANoptosis is a newly recognized inflammatory pathway for programmed cell death (PCD). It participates in regulating the internal environment, homeostasis, and disease process in various complex ways and plays a crucial role in tumor development, but its mechanism of action is still unclear. In this study, we comprehensively analyzed the expression of 14 PANoptosis-related genes (PANRGs) in 28 types of tumors. Most PANRGs are upregulated in tumors, including Z-DNA binding protein 1 (ZBP1), nucleotide-binding oligomerization domain (NOD)-like receptor pyrin domain-containing 3 (NLRP3), caspase (CASP) 1, CASP6, CASP8, PYCARD, FADD, MAP3K7, RNF31, and RBCK1. PANRGs are highly expressed in GBM, LGG, and PAAD, while their levels in ACC are much lower than those in normal tissues. We found that both the CNV and SNV gene sets in BLCA are closely related to survival performance. Subsequently, we conducted clustering and LASSO analysis on each tumor and found that the inhibitory and the stimulating immune checkpoints positively correlate with ZBP1, NLRP3, CASP1, CASP8, and TNFAIP3. The immune infiltration results indicated that KIRC is associated with most infiltrating immune cells. According to the six tumor dryness indicators, PANRGs in LGG show the strongest tumor dryness but have a negative correlation with RNAss. In KIRC, LIHC, and TGCT, most PANRGs play an important role in tumor heterogeneity. Additionally, we analyzed the linear relationship between PANRGs and miRNA and found that MAP3K7 correlates to many miRNAs in most cancers. Finally, we predicted the possible drugs for targeted therapy of the cancers. These data greatly enhance our understanding of the components of cancer and may lead to the discovery of new biomarkers for predicting immunotherapy response and improving the prognosis of cancer patients.

A pleiotropic-epistatic entangelement model of drug response.

Because drug response is multifactorial, graph models are uniquely powerful for comprehending its genetic architecture. We deconstruct drug response into many different and interdependent sub-traits, with each sub-trait controlled by multiple genes that act and interact in a complicated manner. The outcome of drug response is the consequence of multileveled intertwined interactions between pleiotropic effects and epistatic effects. Here, we propose a general statistical physics framework to chart the 3D geometric network that codes how epistasis pleiotropically influences a complete set of sub-traits to shape body-drug interactions. This model can dissect the topological architecture of epistatically induced pleiotropic networks (EiPN) and pleiotropically influenced epistatic networks (PiEN). We analyze and interpret the practical implications of the pleiotropic-epistatic entanglement model for pharmacogenomic studies.

Knockdown of disheveled-associated activator of morphogenesis 2 disrupts cytoskeletal organization and phagocytosis in rat Sertoli cells.

Disheveled-associated activator of morphogenesis 2 (DAAM2) regulates actin polymerization and cell motility. In this study, we investigated the role of DAAM2 in the cytoskeleton and phagocytosis of rat Sertoli cells in vitro and in vivo through siRNA transfection and intratesticular injection. We found that knockdown of DAAM2 significantly attenuated cytoskeletal and tight junction marker expression and reduced the integrity of the Sertoli cell monolayer. In rats, loss of DAAM2 induced disarrangement and deformation of sperms and promoted accumulation of apoptotic sperms in the testis, accompanied by morphological abnormalities in the blood-testis barrier. DAAM2 silencing also reduced the ability of Sertoli cells to engulf apoptotic spermatogenic cells and green fluorescence-labeled beads. RNA sequencing and bioinformatics analysis revealed that phagocytosis and cytoskeleton-related genes and pathways were significantly associated with DAAM2. Our study suggests that DAAM2 may be involved in spermatogenesis possibly by regulating cytoskeleton organization and phagocytosis of Sertoli cells.

Defects of microtubule cytoskeletal organization in NOA human testes.

The importance of actin and microtubule (MT) cytoskeletons in testis function in rodents is known to some extent, but its role in the etiology of azoospermia in humans remains unexplored. Here, we examined if MT cytoskeleton was defective in NOA (non-obstructive azoospermia) testes versus normal human testes based on histopathological, immunofluorescence (IF), and scRNA-Seq transcriptome profiling. Testis biopsy samples from n = 6 normal men versus n = 3 Sertoli cell only (SCO) and n = 3 MA (meiotic arrest) of NOA patients were used for histopathological analysis. IF analysis was also used to examine MT organization across the seminiferous epithelium, investigating the likely involvement of microtubule-associated proteins (MAPs). scRNA-Seq transcriptome profiling datasets from testes of 3 SCO patients versus 3 normal men in public domain in Gene Expression Omnibus (GEO) Sample (GSM) with identifiers were analyzed to examine relevant genes that regulate MT dynamics. NOA testes of MA and SCO patients displayed notable defects in MT organization across the epithelium with extensive truncation, mis-alignments and appeared as collapsed structures near the base of the tubules. These changes are in contrast to MTs in testes of normal men. scRNA-Seq analyses revealed considerable loss of spermatogenesis capacity in SCO testes of NOA patients versus normal men. An array of genes that support MT dynamics displayed considerable changes in expression and in spatial distribution. In summary, defects in MT cytoskeleton were noted in testes of NOA (SCO) patients, possibly mediated by defective spatial expression and/or distribution of MAPs. These changes, in turn, may impede spermatogenesis in SCO testes of NOA patients.

Genetic Architecture of Heterophylly: Single and Multi-Leaf Genome-Wide Association Mapping in Populus euphratica.

Heterophylly is an adaptive strategy used by some plants in response to environmental changes. Due to the lack of representative plants with typical heteromorphic leaves, little is known about the genetic architecture of heterophylly in plants and the genes underlying its control. Here, we investigated the genetic characteristics underlying changes in leaf shape based on the model species, Populus euphratica, which exhibits typical heterophylly. A set of 401,571 single-nucleotide polymorphisms (SNPs) derived from whole-genome sequencing of 860 genotypes were associated with nine leaf traits, which were related to descriptive and shape data using single- and multi-leaf genome-wide association studies (GWAS). Multi-leaf GWAS allows for a more comprehensive understanding of the genetic architecture of heterophylly by considering multiple leaves simultaneously. The single-leaf GWAS detected 140 significant SNPs, whereas the multi-leaf GWAS detected 200 SNP-trait associations. Markers were found across 19 chromosomes, and 21 unique genes were implicated in traits and serve as potential targets for selection. Our results provide novel insights into the genomic architecture of heterophylly, and provide candidate genes for breeding or engineering P. euphratica. Our observations also improve understanding of the intrinsic mechanisms of plant growth, evolution, and adaptation in response to climate change.

BET bromodomain inhibitor JQ1 regulates spermatid development by changing chromatin conformation in mouse spermatogenesis.

As a BET bromodomain inhibitor, JQ1 has been proven have efficacy against a number of different cancers. In terms of male reproduction, JQ1 may be used as a new type of contraceptive, since JQ1 treatment in male mice could lead to germ cell defects and a decrease of sperm motility, moreover, this effect is reversible. However, the mechanism of JQ1 acting on gene regulation in spermatogenesis remains unclear. Here, we performed single-cell RNA sequencing (scRNA-seq) on mouse testes treated with JQ1 or vehicle control to determine the transcriptional regulatory function of JQ1 in spermatogenesis at the single cell resolution. We confirmed that JQ1 treatment could increase the numbers of somatic cells and spermatocytes and decrease the numbers of spermatid cells. Gene Ontology (GO) analysis demonstrated that differentially expressed genes which were down-regulated after JQ1 injection were mainly enriched in "DNA conformation change" biological process in early developmental germ cells and "spermatid development" biological process in spermatid cells. ATAC-seq data further confirmed that JQ1 injection could change the open state of chromatin. In addition, JQ1 could change the numbers of accessible meiotic DNA double-stranded break sites and the types of transcription factor motif that functioned in pachytene spermatocytes and round spermatids. The multi-omics analysis revealed that JQ1 had the ability to regulate gene transcription by changing chromatin conformation in mouse spermatogenesis, which would potentiate the availability of JQ1 in male contraceptive.

Mapping Floral Genetic Architecture in Prunus mume, an Ornamental Woody Plant.

Floral traits are both evolutionarily and economically relevant for ornamental plants. However, their underlying genetic architecture, especially in woody ornamental plants, is still poorly understood. We perform mapping experiments aimed at identifying specific quantitative trait loci (QTLs) that control the size, shape, architecture, color, and timing of flowers in mei (Prunus mume). We find that the narrow region of chromosome 1 (5-15 Mb) contains a number of floral QTLs. Most QTLs detected from this mapping study are annotated to candidate genes that regulate various biological functions toward the floral formation. We identify strong pleiotropic control on different aspects of flower morphology (including shape, petal number, pistil number, petal color, and calyx color) and flower timing, but find different genetic systems that mediate whether a flower produces pistils and how many pistils a flower produces. We find that many floral QTLs display pleiotropic effects on shoot length growth but shoot radial growth, implicating a possible association of floral display with light capture. We conduct a transcriptomic study to characterize the genomic signature of floral QTLs expressed in mei. Our mapping results about the genetic control of floral features make it promising to select superior varieties for mei carrying flowers of ornamental value.

A graph model of combination therapies.

The simultaneous use of multiple medications causes drug-drug interactions (DDI) that impact therapeutic efficacy. Here, we argue that graph theory, in conjunction with game theory and ecosystem theory, can address this issue. We treat the coexistence of multiple drugs as a system in which DDI is modeled by game theory. We develop an ordinary differential equation model to characterize how the concentration of a drug changes as a result of its independent capacity and the dependent influence of other drugs through the metabolic response of the host. We coalesce all drugs into personalized and context-specific networks, which can reveal key DDI determinants of therapeutical efficacy. Our model can quantify drug synergy and antagonism and test the translational success of combination therapies to the clinic.

Peptides from the croceine croaker (Larimichthys crocea) swim bladder attenuate busulfan-induced oligoasthenospermia in mice.

CONTEXT: The swim bladder of the croceine croaker is believed to have a therapeutic effect on various diseases. However, there is no research about its effect on mammalian spermatogenesis. OBJECTIVE: We investigated the swim bladder peptides (SBPs) effect on busulfan-induced oligoasthenospermia in mice. MATERIALS AND METHODS: We first extracted SBP from protein hydrolysate of the croceine croaker swim bladder, and then five groups of ICR male mice were randomly assigned: control, control + SBP 60 mg/kg, busulfan, busulfan + SBP 30 mg/kg and busulfan + SBP 60 mg/kg. Mice received bilateral intratesticular injections of busulfan to establish oligoasthenospermia model. After treatment with SBP for 4 weeks, testis and epididymis were collected from all mice for further analysis. RESULTS: After treatment with SBP 30-60 mg/kg for 4 weeks, epididymal sperm concentration and motility increased by 3.9-9.6- and 1.9-2.4-fold than those of oligoasthenospermia mice induced by busulfan. Meanwhile, histology showed that spermatogenic cells decreased, leading to increased lumen diameters and vacuolization in the busulfan group. These features were reversed by SBP treatment. RNA-sequencing analysis revealed that, compared with the busulfan group, Lin28b and Igf2bp1 expression related to germ cell proliferation, increased with a >1.5-fold change after SBP treatment. Additionally, PGK2 and Cfap69 mRNAs associated with sperm motility, also increased with a >1.5-fold change. Furthermore, these findings were validated by quantitative real-time PCR and Western blotting. DISCUSSION AND CONCLUSIONS: This is the first reported evidence for the therapeutic effect of SBP on oligoasthenospermia. SBP may be a promising drug for oligoasthenospermia in humans.

Nanopore Sequencing and Hi-C Based De Novo Assembly of Trachidermus fasciatus Genome.

Trachidermus fasciatus is a roughskin sculpin fish widespread across the coastal areas of East Asia. Due to environmental destruction and overfishing, the population of this species is under threat. In order to protect this endangered species, it is important to have the genome sequenced. Reference genomes are essential for studying population genetics, domestic farming, and genetic resource protection. However, currently, no reference genome is available for Trachidermus fasciatus, and this has greatly hindered the research on this species. In this study, we integrated nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C methods to thoroughly assemble the Trachidermus fasciatus genome. Our results provided a chromosome-level high-quality genome assembly with a predicted genome size of 542.6 Mbp (2n = 40) and a scaffold N50 of 24.9 Mbp. The BUSCO value for genome assembly completeness was higher than 96%, and the single-base accuracy was 99.997%. Based on EVM-StringTie genome annotation, a total of 19,147 protein-coding genes were identified, including 35,093 mRNA transcripts. In addition, a novel gene-finding strategy named RNR was introduced, and in total, 51 (82) novel genes (transcripts) were identified. Lastly, we present here the first reference genome for Trachidermus fasciatus; this sequence is expected to greatly facilitate future research on this species.

Spatio-temporal landscape of mouse epididymal cells and specific mitochondria-rich segments defined by large-scale single-cell RNA-seq.

Spermatozoa acquire their fertilizing ability and forward motility during epididymal transit, suggesting the importance of the epididymis. Although the cell atlas of the epididymis was reported recently, the heterogeneity of the cells and the gene expression profile in the epididymal tube are still largely unknown. Considering single-cell RNA sequencing results, we thoroughly studied the cell composition, spatio-temporal differences in differentially expressed genes (DEGs) in epididymal segments and mitochondria throughout the epididymis with sufficient cell numbers. In total, 40,623 cells were detected and further clustered into 8 identified cell populations. Focused analyses revealed the subpopulations of principal cells, basal cells, clear/narrow cells, and halo/T cells. Notably, two subtypes of principal cells, the Prc7 and Prc8 subpopulations were enriched as stereocilia-like cells according to GO analysis. Further analysis demonstrated the spatially specific pattern of the DEGs in each cell cluster. Unexpectedly, the abundance of mitochondria and mitochondrial transcription (MT) was found to be higher in the corpus and cauda epididymis than in the caput epididymis by scRNA-seq, immunostaining, and qPCR validation. In addition, the spatio-temporal profile of the DEGs from the P42 and P56 epididymis, including transiting spermatozoa, was depicted. Overall, our study presented the single-cell transcriptome atlas of the mouse epididymis and revealed the novel distribution pattern of mitochondria and key genes that may be linked to sperm functionalities in the first wave and subsequent wave of sperm, providing a roadmap to be emulated in efforts to achieve sperm maturation regulation in the epididymis.

Genetic adaptation of Tibetan poplar (Populus szechuanica var. tibetica) to high altitudes on the Qinghai-Tibetan Plateau.

Plant adaptation to high altitudes has long been a substantial focus of ecological and evolutionary research. However, the genetic mechanisms underlying such adaptation remain poorly understood. Here, we address this issue by sampling, genotyping, and comparing populations of Tibetan poplar, Populus szechuanica var. tibetica, distributed from low (~2,000 m) to high altitudes (~3,000 m) of Sejila Mountain on the Qinghai-Tibet Plateau. Population structure analyses allow clear classification of two groups according to their altitudinal distributions. However, in contrast to the genetic variation within each population, differences between the two populations only explain a small portion of the total genetic variation (3.64%). We identified asymmetrical gene flow from high- to low-altitude populations. Integrating population genomic and landscape genomic analyses, we detected two hotspot regions, one containing four genes associated with altitudinal variation, and the other containing ten genes associated with response to solar radiation. These genes participate in abiotic stress resistance and regulation of reproductive processes. Our results provide insight into the genetic mechanisms underlying high-altitude adaptation in Tibetan poplar.

Two resveratrol analogs, pinosylvin and 4,4'-dihydroxystilbene, improve oligoasthenospermia in a mouse model by attenuating oxidative stress via the Nrf2-ARE pathway.

Two synthesized resveratrol analogs from our laboratory, namely pinosylvin (3,5-dihydroxy-trans-stilbene, PIN) and 4,4'-dihydroxystilbene (DHS), have been carefully evaluated for treatment of oligoasthenospermia. Recent studies have demonstrated that PIN and DHS improved sperm quality in the mouse. However, the mechanism of action of PIN and DHS on oligoasthenospermia remains unknown. Herein, we investigated the mechanistic basis for improvements in sperm parameters by PIN and DHS in a mouse model of oligoasthenospermia induced by treatment with busulfan (BUS) at 6 mg/kg b.w.. Two weeks following busulfan treatment, mice were administered different concentrations of PIN or DHS daily for 2 consecutive weeks. Thereafter, epididymal sperm concentration and motility were determined, and histopathology of the testes was performed. Serum hormone levels including testosterone (T), luteinizing hormone (LH), and follicle stimulating hormone (FSH) were measured using corresponding specific enzyme-linked immunosorbent assay (ELISA) kits. Testicular mRNA expression profiles were determined by RNA sequencing analysis. These findings were validated by quantitative real-time PCR, western blotting and ELISA. Both PIN and DHS improved the epididymal sperm concentration and motility, enhanced testosterone levels, and promoted testicular morphological recovery following BUS treatment. PIN treatment was found to significantly reduce oxidative stress via the nuclear factor erythroid 2-related factor 2 (Nrf2)-antioxidant response element (ARE)-dependent antioxidant, glutathione peroxidase 3. DHS treatment significantly reduced oxidative stress via the Nrf2-ARE-dependent antioxidants glutathione S-transferase theta 2 and glutathione S-transferase omega 2. In summary, PIN and DHS ameliorated oligoasthenospermia in this mouse model by attenuating oxidative stress via the Nrf2-ARE pathway.

Corrigendum: A Computational Model for Inferring QTL Control Networks Underlying Developmental Covariation.

[This corrects the article DOI: 10.3389/fpls.2019.01557.].

A rewiring model of intratumoral interaction networks.

Intratumoral heterogeneity (ITH) has been regarded as a key cause of the failure and resistance of cancer therapy, but how it behaves and functions remains unclear. Advances in single-cell analysis have facilitated the collection of a massive amount of data about genetic and molecular states of individual cancer cells, providing a fuel to dissect the mechanistic organization of ITH at the molecular, metabolic and positional level. Taking advantage of these data, we propose a computational model to rewire up a topological network of cell-cell interdependences and interactions that operate within a tumor mass. The model is grounded on the premise of game theory that each interactive cell (player) strives to maximize its fitness by pursuing a "rational self-interest" strategy, war or peace, in a way that senses and alters other cells to respond properly. By integrating this idea with genome-wide association studies for intratumoral cells, the model is equipped with a capacity to visualize, annotate and quantify how somatic mutations mediate ITH and the network of intratumoral interactions. Taken together, the model provides a topological flow by which cancer cells within a tumor cooperate or compete with each other to downstream pathogenesis. This topological flow can be potentially used as a blueprint for genetically intervening the pattern and strength of cell-cell interactions towards cancer control.

A Regulatory Network for miR156-SPL Module in Arabidopsis thaliana.

Vegetative phase changes in plants describes the transition between juvenile and adult phases of vegetative growth before flowering. It is one of the most fundamental mechanisms for plants to sense developmental signals, presenting a complex process involving many still-unknown determinants. Several studies in annual and perennial plants have identified the conservative roles of miR156 and its targets, SBP/SPL genes, in guiding the switch of plant growth from juvenile to adult phases. Here, we review recent progress in understanding the regulation of miR156 expression and how miR156-SPLs mediated plant age affect other processes in Arabidopsis. Powerful high-throughput sequencing techniques have provided rich data to systematically study the regulatory mechanisms of miR156 regulation network. From this data, we draw an expanded miR156-regulated network that links plant developmental transition and other fundamental biological processes, gaining novel and broad insight into the molecular mechanisms of plant-age-related processes in Arabidopsis.

A Drive to Driven Model of Mapping Intraspecific Interaction Networks.

Community ecology theory suggests that an individual's phenotype is determined by the phenotypes of its coexisting members to the extent at which this process can shape community evolution. Here, we develop a mapping theory to identify interaction quantitative trait loci (QTL) governing inter-individual dependence. We mathematically formulate the decision-making strategy of interacting individuals. We integrate these mathematical descriptors into a statistical procedure, enabling the joint characterization of how QTL drive the strengths of ecological interactions and how the genetic architecture of QTL is driven by ecological networks. In three fish full-sib mapping experiments, we identify a set of genome-wide QTL that control a range of societal behaviors, including mutualism, altruism, aggression, and antagonism, and find that these intraspecific interactions increase the genetic variation of body mass by about 50%. We showcase how the interaction QTL can be used as editors to reconstruct and engineer new social networks for ecological communities.

A dissection model for mapping complex traits.

Many quantitative traits are composites of other traits that contribute differentially to genetic variation. Quantitative trait locus (QTL) mapping of these composite traits can benefit by incorporating the mechanistic process of how their formation is mediated by the underlying components. We propose a dissection model by which to map these interconnected components traits under a joint likelihood setting. The model can test how a composite trait is determined by pleiotropic QTLs for its component traits or jointly by different sets of QTLs each responsible for a different component. The model can visualize the pattern of time-varying genetic effects for individual components and their impacts on composite traits. The dissection model was used to map two composite traits, stemwood volume growth decomposed into its stem height, stem diameter and stem form components for Euramerican poplar adult trees, and total lateral root length constituted by its average lateral root length and lateral root number components for Euphrates poplar seedlings. We found the pattern of how QTLs for different components contribute to phenotypic variation in composite traits. The detailed understanding of the genetic machineries of composite traits will not only help in the design of molecular breeding in plants and animals, but also shed light on the evolutionary processes of quantitative traits under natural selection.

A Computational Model for Inferring QTL Control Networks Underlying Developmental Covariation.

How one trait developmentally varies as a function of others shapes a spectrum of biological phenomena. Despite its importance to trait dissection, the understanding of whether and how genes mediate such developmental covariation is poorly understood. We integrate developmental allometry equations into the functional mapping framework to map specific QTLs that govern the correlated development of different traits. Based on evolutionary game theory, we assemble and contextualize these QTLs into an intricate but organized network coded by bidirectional, signed, and weighted QTL-QTL interactions. We use this approach to map shoot height-diameter allometry QTLs in an ornamental woody species, mei (Prunus mume). We detect "pioneering" QTLs (piQTLs) and "maintaining" QTLs (miQTLs) that determine how shoot height varies with diameter and how shoot diameter varies with height, respectively. The QTL networks inferred can visualize how each piQTL regulates others to promote height growth at a cost of diameter growth, how miQTL regulates others to benefit radial growth at a cost of height growth, and how piQTLs and miQTLs regulate each other to form a pleiotropic web of primary and secondary growth in trees. Our approach provides a unique gateway to explore the genetic architecture of developmental covariation, a widespread phenomenon in nature.