The liver is the largest solid organ in the body and is critical for metabolic and immune functions. However, little is known about the cells that make up the human liver and its immune microenvironment. Here we report a map of the cellular landscape of the human liver using single-cell RNA sequencing. We provide the transcriptional profiles of 8444 parenchymal and non-parenchymal cells obtained from the fractionation of fresh hepatic tissue from five human livers. Using gene expression patterns, flow cytometry, and immunohistochemical examinations, we identify 20 discrete cell populations of hepatocytes, endothelial cells, cholangiocytes, hepatic stellate cells, B cells, conventional and non-conventional T cells, NK-like cells, and distinct intrahepatic monocyte/macrophage populations. Together, our study presents a comprehensive view of the human liver at single-cell resolution that outlines the characteristics of resident cells in the liver, and in particular provides a map of the human hepatic immune microenvironment.

The critical functions of the human liver are coordinated through the interactions of hepatic parenchymal and non-parenchymal cells. Recent advances in single-cell transcriptional approaches have enabled an examination of the human liver with unprecedented resolution. However, dissociation-related cell perturbation can limit the ability to fully capture the human liver's parenchymal cell fraction, which limits the ability to comprehensively profile this organ. Here, we report the transcriptional landscape of 73,295 cells from the human liver using matched single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq). The addition of snRNA-seq enabled the characterization of interzonal hepatocytes at a single-cell resolution, revealed the presence of rare subtypes of liver mesenchymal cells, and facilitated the detection of cholangiocyte progenitors that had only been observed during in vitro differentiation experiments. However, T and B lymphocytes and natural killer cells were only distinguishable using scRNA-seq, highlighting the importance of applying both technologies to obtain a complete map of tissue-resident cell types. We validated the distinct spatial distribution of the hepatocyte, cholangiocyte, and mesenchymal cell populations by an independent spatial transcriptomics data set and immunohistochemistry. Conclusion: Our study provides a systematic comparison of the transcriptomes captured by scRNA-seq and snRNA-seq and delivers a high-resolution map of the parenchymal cell populations in the healthy human liver.

Despite extensive analysis of pRB phosphorylation in vitro, how this modification influences development and homeostasis in vivo is unclear. Here, we show that homozygous Rb∆K4 and Rb∆K7 knock-in mice, in which either four or all seven phosphorylation sites in the C-terminal region of pRb, respectively, have been abolished by Ser/Thr-to-Ala substitutions, undergo normal embryogenesis and early development, notwithstanding suppressed phosphorylation of additional upstream sites. Whereas Rb∆K4 mice exhibit telomere attrition but no other abnormalities, Rb∆K7 mice are smaller and display additional hallmarks of premature aging including infertility, kyphosis, and diabetes, indicating an accumulative effect of blocking pRb phosphorylation. Diabetes in Rb∆K7 mice is insulin-sensitive and associated with failure of quiescent pancreatic β-cells to re-enter the cell cycle in response to mitogens, resulting in induction of DNA damage response (DDR), senescence-associated secretory phenotype (SASP), and reduced pancreatic islet mass and circulating insulin level. Pre-treatment with the epigenetic regulator vitamin C reduces DDR, increases cell cycle re-entry, improves islet morphology, and attenuates diabetes. These results have direct implications for cell cycle regulation, CDK-inhibitor therapeutics, diabetes, and longevity.

PTEN loss and PIK3CA activation both promote the accumulation of phosphatidylinositol (3, 4, 5)-trisphosphate (PIP3). While these proteins also have distinct biochemical functions, beyond the regulation of PIP3, little is known about the consequences of these differences in vivo. Here, we directly compared cancer signalling in mammary tumors from MMTV-Cre:Ptenf/f and MMTV-Cre:Pik3ca(LSL-H1047R) mice. Using unsupervised hierarchical clustering we found that whereas MMTV-Cre:Pik3ca(LSL-H1047R)-derived tumors fall into two separate groups, designated squamous-likeEx and class14(Ex), MMTV-Cre:Ptenf/f tumors cluster as one group together with PIK3CA(H1047R) class14(Ex), exhibiting a 'luminal' expression profile. Gene Set Enrichment Analysis (GSEA) of Pten(Δ)ˆ† and PIK3CA(H1047R) class14(Ex) tumors revealed very similar profiles of signalling pathways as well as some interesting differences. Analysis of 18 signalling signatures revealed that PI3K signalling is significantly induced whereas EGFR signalling is significantly reduced in Pten(∆) versus PIK3CA(H1047R) tumors. Thus, Pten(∆) and PIK3CA(H1047R) tumors exhibit discernable differences that may impact tumorigenesis and response to therapy.

Melanoma differentiation associated gene-7/interleukin-24 (mda-7/IL-24) encodes a tumor suppressor gene implicated in the growth of various tumor types including breast cancer. We previously demonstrated that recombinant adenovirus-mediated mda-7/IL-24 expression in the mammary glands of carcinogen-treated (methylnitrosourea, MNU) rats suppressed mammary tumor development. Since most MNU-induced tumors in rats contain activating mutations in Ha-ras, which arenot frequently detected in humans, we presently examined the effect of MDA-7/IL-24 on Her2/Neu-induced mammary tumors, in which the RAS pathway is induced. We generated tet-inducible MDA-7/IL-24 transgenic mice and crossed them with Her2/Neu transgenic mice. Triple compound transgenic mice treated with doxycycline exhibited a strong inhibition of tumor development, demonstrating tumor suppressor activity by MDA-7/IL-24 in immune-competent mice. MDA-7/IL-24 induction also inhibited growth of tumors generated following injection of Her2/Neu tumor cells isolated from triple compound transgenic mice that had not been treated with doxycycline, into the mammary fat pads of isogenic FVB mice. Despite initial growth suppression, tumors in triple compound transgenic mice lost mda-7/IL-24 expression and grew, albeit after longer latency, indicating that continuous presence of this cytokine within tumor microenvironment is crucial to sustain tumor inhibitory activity. Mechanistically, MDA-7/IL-24 exerted its tumor suppression effect on HER2+ breast cancer cells, at least in part, through PERP, a member of PMP-22 family with growth arrest and apoptosis-inducing capacity. Overall, our results establish mda-7/IL-24 as a suppressor of mammary tumor development and provide a rationale for using this cytokine in the prevention/treatment of human breast cancer.

Defects in a form of noncanonical autophagy, known as LC3-associated phagocytosis (LAP), lead to increased inflammatory pathology during fungal infection. Although LAP contributes to fungal degradation, the molecular mechanisms underlying LAP-mediated modulation of inflammation are unknown. We describe a mechanism by which inflammation is regulated during LAP through the death-associated protein kinase 1 (DAPK1). The ATF6/C/EBP-β/DAPK1 axis activated by IFN-γ not only mediates LAP to Aspergillus fumigatus but also concomitantly inhibits Nod-like receptor protein 3 (NLRP3) activation and restrains pathogenic inflammation. In mouse models and patient samples of chronic granulomatous disease, which exhibit defective autophagy and increased inflammasome activity, IFN-γ restores reduced DAPK1 activity and dampens fungal growth. Additionally, in a cohort of hematopoietic stem cell-transplanted patients, a genetic DAPK1 deficiency is associated with increased inflammation and heightened aspergillosis susceptibility. Thus, DAPK1 is a potential drugable player in regulating the inflammatory response during fungal clearance initiated by IFN-γ.

The eukaryotic elongation factor-2 kinase, eEF2K, which restricts protein translation elongation, has been identified as a potential therapeutic target for diverse types of malignancies including triple negative breast cancer (TNBC). However, the contexts in which eEF2K inhibition is essential in TNBC and its consequences on the proteome are largely unknown. Here we show that genetic or pharmacological inhibition of eEF2K cooperated with glutamine (Gln) starvation, and synergized with glutaminase (GLS1) inhibitors to suppress growth of diverse TNBC cell lines. eEF2K inhibition also synergized with depletion of eukaryotic translation initiation factor 4E-binding protein 1 (eIF4EBP1; 4EBP1), a suppressor of eukaryotic protein translation initiation factor 4E (eIF4E), to induce c-MYC and Cyclin D1 expression, yet attenuate growth of TNBC cells. Proteomic analysis revealed that whereas eEF2K depletion alone uniquely induced Cyclin Dependent Kinase 1 (CDK1) and 6 (CDK6), combined depletion of eEF2K and 4EBP1 resulted in overlapping effects on the proteome, with the highest impact on the 'Collagen containing extracellular matrix' pathway (e.g. COL1A1), as well as the amino-acid transporter, SLC7A5/LAT1, suggesting a regulatory loop via mTORC1. In addition, combined depletion of eEF2K and 4EBP1 indirectly reduced the levels of IFN-dependent innate immune response-related factors. Thus, eEF2K inhibition triggers cell cycle arrest/death under unfavourable metabolic conditions such as Gln-starvation/GLS1 inhibition or 4EBP1 depletion, uncovering new therapeutic avenues for TNBC and underscoring a pressing need for clinically relevant eEF2K inhibitors.

Bisphosphonates are the most important class of antiresorptive agents used against osteoclast-mediated bone loss, and, more recently, in oncology. These compounds have high affinity for calcium ions (Ca(2+)) and therefore target bone mineral, where they appear to be internalized selectively by bone-resorbing osteoclasts and inhibit osteoclast function. They are extensively used in healthcare, however they are affected by severe side effects; pharmacological properties of bisphosphonates depend on their molecular structure, which is frequently the cause of poor intestinal adsorption and low distribution. In this work we synthesized six novel bisphosphonate compounds having a variably substituted indole moiety to evaluate their extra- and intracellular calcium chelating ability in PE/CA-PJ15 cells. Preliminary in silico and in vitro ADME studies were also performed and the results suggested that the indole moiety plays an important role in cell permeability and metabolism properties.

CDK4/6 inhibitors are effective against cancer cells expressing the tumor suppressor RB1, but not RB1-deficient cells, posing the challenge of how to target RB1 loss. In triple-negative breast cancer (TNBC), RB1 and PTEN are frequently inactivated together with TP53. We performed kinome/phosphatase inhibitor screens on primary mouse Rb/p53-, Pten/p53-, and human RB1/PTEN/TP53-deficient TNBC cell lines and identified CDC25 phosphatase as a common target. Pharmacological or genetic inhibition of CDC25 suppressed growth of RB1-deficient TNBC cells that are resistant to combined CDK4/6 plus CDK2 inhibition. Minimal cooperation was observed in vitro between CDC25 antagonists and CDK1, CDK2, or CDK4/6 inhibitors, but strong synergy with WEE1 inhibition was apparent. In accordance with increased PI3K signaling following long-term CDC25 inhibition, CDC25 and PI3K inhibitors effectively synergized to suppress TNBC growth both in vitro and in xenotransplantation models. These results provide a rationale for the development of CDC25-based therapies for diverse RB1/PTEN/TP53-deficient and -proficient TNBCs.

Pineoblastoma is a rare pediatric cancer induced by germline mutations in the tumor suppressors RB1 or DICER1. Presence of leptomeningeal metastases is indicative of poor prognosis. Here we report that inactivation of Rb plus p53 via a WAP-Cre transgene, commonly used to target the mammary gland during pregnancy, induces metastatic pineoblastoma resembling the human disease with 100% penetrance. A stabilizing mutation rather than deletion of p53 accelerates metastatic dissemination. Deletion of Dicer1 plus p53 via WAP-Cre also predisposes to pineoblastoma, albeit with lower penetrance. In silico analysis predicts tricyclic antidepressants such as nortriptyline as potential therapeutics for both pineoblastoma models. Nortriptyline disrupts the lysosome, leading to accumulation of non-functional autophagosome, cathepsin B release and pineoblastoma cell death. Nortriptyline further synergizes with the antineoplastic drug gemcitabine to effectively suppress pineoblastoma in our preclinical models, offering new modality for this lethal childhood malignancy.

Metastatic breast-cancer is a major cause of death in women worldwide, yet the relationship between oncogenic drivers that promote metastatic versus primary cancer is still contentious. To elucidate this relationship in treatment-naive animals, we hereby describe mammary-specific transposon-mutagenesis screens in female mice together with loss-of-function Rb, which is frequently inactivated in breast-cancer. We report gene-centric common insertion-sites (gCIS) that are enriched in primary-tumors, in metastases or shared by both compartments. Shared-gCIS comprise a major MET-RAS network, whereas metastasis-gCIS form three additional hubs: Rho-signaling, Ubiquitination and RNA-processing. Pathway analysis of four clinical cohorts with paired primary-tumors and metastases reveals similar organization in human breast-cancer with subtype-specific shared-drivers (e.g. RB1-loss, TP53-loss, high MET, RAS, ER), primary-enriched (EGFR, TGFβ and STAT3) and metastasis-enriched (RHO, PI3K) oncogenic signaling. Inhibitors of RB1-deficiency or MET plus RHO-signaling cooperate to block cell migration and drive tumor cell-death. Thus, targeting shared- and metastasis- but not primary-enriched derivers offers a rational avenue to prevent metastatic breast-cancer.

Mammalian cells can utilize hydrogen sulfide (H2S) to support mitochondrial respiration. The aim of our study was to explore the potential role of S-sulfhydration (a H2S-induced posttranslational modification, also known as S-persulfidation) of the mitochondrial inner membrane protein ATP synthase (F1F0 ATP synthase/Complex V) in the regulation of mitochondrial bioenergetics. Using a biotin switch assay, we have detected S-sulfhydration of the α subunit (ATP5A1) of ATP synthase in response to exposure to H2S in vitro. The H2S generator compound NaHS induced S-sulfhydration of ATP5A1 in HepG2 and HEK293 cell lysates in a concentration-dependent manner (50-300μM). The activity of immunocaptured mitochondrial ATP synthase enzyme isolated from HepG2 and HEK293 cells was stimulated by NaHS at low concentrations (10-100nM). Site-directed mutagenesis of ATP5A1 in HEK293 cells demonstrated that cysteine residues at positions 244 and 294 are subject to S-sulfhydration. The double mutant ATP synthase protein (C244S/C294S) showed a significantly reduced enzyme activity compared to control and the single-cysteine-mutated recombinant proteins (C244S or C294S). To determine whether endogenous H2S plays a role in the basal S-sulfhydration of ATP synthase in vivo, we compared liver tissues harvested from wild-type mice and mice deficient in cystathionine-gamma-lyase (CSE, one of the three principal mammalian H2S-producing enzymes). Significantly reduced S-sulfhydration of ATP5A1 was observed in liver homogenates of CSE-/- mice, compared to wild-type mice, suggesting a physiological role for CSE-derived endogenous H2S production in the S-sulfhydration of ATP synthase. Various forms of critical illness (including burn injury) upregulate H2S-producing enzymes and stimulate H2S biosynthesis. In liver tissues collected from mice subjected to burn injury, we detected an increased S-sulfhydration of ATP5A1 at the early time points post-burn. At later time points (when systemic H2S levels decrease) S-sulfhydration of ATP5A1 decreased as well. In conclusion, H2S induces S-sulfhydration of ATP5A1 at C244 and C294. This post-translational modification may be a physiological mechanism to maintain ATP synthase in a physiologically activated state, thereby supporting mitochondrial bioenergetics. The sulfhydration of ATP synthase may be a dynamic process, which may be regulated by endogenous H2S levels under various pathophysiological conditions.

Triple negative breast cancer (TNBC) includes basal-like and claudin-low subtypes for which only chemotherapy and radiation therapy are currently available. The retinoblastoma (RB1) tumor suppressor is frequently lost in human TNBC. Knockdown of RB1 in luminal BC cells was shown to affect response to endocrine, radiation and several antineoplastic drugs. However, the effect of RB1 status on radiation and chemo-sensitivity in TNBC cells and whether RB1 status affects response to divergent or specific treatment are unknown. Using multiple basal-like and claudin-low cell lines, we hereby demonstrate that RB-negative TNBC cell lines are highly sensitive to gamma-irradiation, and moderately more sensitive to doxorubicin and methotrexate compared to RB-positive TNBC cell lines. In contrast, RB1 status did not affect sensitivity of TNBC cells to multiple other drugs including cisplatin (CDDP), 5-fluorouracil, idarubicin, epirubicin, PRIMA-1(met), fludarabine and PD-0332991, some of which are used to treat TNBC patients. Moreover, a non-biased screen of ∼3400 compounds, including FDA-approved drugs, revealed similar sensitivity of RB-proficient and -deficient TNBC cells. Finally, ESA(+)/CD24(-/low)/CD44(+) cancer stem cells from RB-negative TNBC lines were consistently more sensitive to gamma-irradiation than RB-positive lines, whereas the effect of chemotherapy on the cancer stem cell fraction varied irrespective of RB1 expression. Our results suggest that patients carrying RB-deficient TNBCs would benefit from gamma-irradiation as well as doxorubicin and methotrexate therapy, but not necessarily from many other anti-neoplastic drugs.

Triple-negative breast cancer (TNBC) represents an aggressive subtype, for which radiation and chemotherapy are the only options. Here we describe the identification of disulfiram, an FDA-approved drug used to treat alcoholism, as well as the related compound thiram, as the most potent growth inhibitors following high-throughput screens of 3185 compounds against multiple TNBC cell lines. The average IC50 for disulfiram was ~300 nM. Drug affinity responsive target stability (DARTS) analysis identified IQ motif-containing factors IQGAP1 and MYH9 as direct binding targets of disulfiram. Indeed, knockdown of these factors reduced, though did not completely abolish, cell growth. Combination treatment with 4 different drugs commonly used to treat TNBC revealed that disulfiram synergizes most effectively with doxorubicin to inhibit cell growth of TNBC cells. Disulfiram and doxorubicin cooperated to induce cell death as well as cellular senescence, and targeted the ESA(+)/CD24(-/low)/CD44(+) cancer stem cell population. Our results suggest that disulfiram may be repurposed to treat TNBC in combination with doxorubicin.

The tumor suppressors Pten and p53 are frequently lost in breast cancer, yet the consequences of their combined inactivation are poorly understood. Here, we show that mammary-specific deletion of Pten via WAP-Cre, which targets alveolar progenitors, induced tumors with shortened latency compared to those induced by MMTV-Cre, which targets basal/luminal progenitors. Combined Pten-p53 mutations accelerated formation of claudin-low, triple-negative-like breast cancer (TNBC) that exhibited hyper-activated AKT signaling and more mesenchymal features relative to Pten or p53 single-mutant tumors. Twenty-four genes that were significantly and differentially expressed between WAP-Cre:Pten/p53 and MMTV-Cre:Pten/p53 tumors predicted poor survival for claudin-low patients. Kinome screens identified eukaryotic elongation factor-2 kinase (eEF2K) inhibitors as more potent than PI3K/AKT/mTOR inhibitors on both mouse and human Pten/p53-deficient TNBC cells. Sensitivity to eEF2K inhibition correlated with AKT pathway activity. eEF2K monotherapy suppressed growth of Pten/p53-deficient TNBC xenografts in vivo and cooperated with doxorubicin to efficiently kill tumor cells in vitro. Our results identify a prognostic signature for claudin-low patients and provide a rationale for using eEF2K inhibitors for treatment of TNBC with elevated AKT signaling.

Liver sinusoidal endothelial cells (LSECs) form a highly specialized microvasculature that plays a critical role in liver function and disease. To better understand this role, we developed a strategy to generate LSECs from human pluripotent stem cells (hPSCs) by first optimizing the specification of arterial and venous angioblasts and derivative endothelial populations. Induction of a LSEC-like fate by hypoxia, cyclic AMP (cAMP) agonism, and transforming growth factor β (TGF-β) inhibition revealed that venous endothelial cells responded more rapidly and robustly than the arterial cells to upregulate LSEC markers and functions in vitro. Upon intrahepatic transplantation in neonates, venous angioblasts engrafted the liver and generated mature, fenestrated LSECs with scavenger functions and molecular profiles of primary human LSECs. When transplanted into the liver of adult mice, angioblasts efficiently gave rise to mature LSECs with robust factor VIII (FVIII) production. Humanization of the murine liver with hPSC-derived LSECs provides a tractable system for studying the biology of this key liver cell type.

Mild traumatic brain injury (mTBI) results in broad neurological symptoms and an increased risk of being diagnosed with a neurodegenerative disease later in life. While the immediate oxidative stress response and post-mortem pathology of the injured brain has been well studied, it remains unclear how early pathogenic changes may drive persistent symptoms and confer susceptibility to neurodegeneration. In this study we have used a mouse model of repeated mTBI (rmTBI) to identify early gene expression changes at 24 h or 7 days post-injury (7 dpi). At 24 h post-injury, gene expression of rmTBI mice shows activation of the DNA damage response (DDR) towards double strand DNA breaks, altered calcium and cell-cell signalling, and inhibition of cell death pathways. By 7 dpi, rmTBI mice had a gene expression signature consistent with induction of cellular senescence, activation of neurodegenerative processes, and inhibition of the DDR. At both timepoints gliosis, microgliosis, and axonal damage were evident in the absence of any gross lesion, and by 7 dpi rmTBI also mice had elevated levels of IL1β, p21, 53BP1, DNA2, and p53, supportive of DNA damage-induced cellular senescence. These gene expression changes reflect establishment of processes usually linked to brain aging and suggests that cellular senescence occurs early and most likely prior to the accumulation of toxic proteins. These molecular changes were accompanied by spatial learning and memory deficits in the Morris water maze. To conclude, we have identified DNA damage-induced cellular senescence as a repercussion of repeated mild traumatic brain injury which correlates with cognitive impairment. Pathways involved in senescence may represent viable treatment targets of post-concussive syndrome. Senescence has been proposed to promote neurodegeneration and appears as an effective target to prevent long-term complications of mTBI, such as chronic traumatic encephalopathy and other related neurodegenerative pathologies.

Systematically investigating the scores of genes mutated in cancer and discerning disease drivers from inconsequential bystanders is a prerequisite for precision medicine but remains challenging. Here, we developed a somatic CRISPR/Cas9 mutagenesis screen to study 215 recurrent "long-tail" breast cancer genes, which revealed epigenetic regulation as a major tumor-suppressive mechanism. We report that components of the BAP1 and COMPASS-like complexes, including KMT2C/D, KDM6A, BAP1, and ASXL1/2 ("EpiDrivers"), cooperate with PIK3CAH1047R to transform mouse and human breast epithelial cells. Mechanistically, we find that activation of PIK3CAH1047R and concomitant EpiDriver loss triggered an alveolar-like lineage conversion of basal mammary epithelial cells and accelerated formation of luminal-like tumors, suggesting a basal origin for luminal tumors. EpiDriver mutations are found in ∼39% of human breast cancers, and ∼50% of ductal carcinoma in situ express casein, suggesting that lineage infidelity and alveogenic mimicry may significantly contribute to early steps of breast cancer etiology.

Triple-negative breast cancer (TNBC), an aggressive disease comprising several subtypes including basal-like and claudin-low, involves frequent deletions or point mutations in TP53, as well as loss of PTEN. We previously showed that combined deletion of both tumor suppressors in the mouse mammary epithelium invariably induced claudin-low-like TNBC. The effect of p53 mutation plus Pten deletion on mammary tumorigenesis and whether this combination can induce basal-like TNBC in the mouse are unknown.

Requiring water and minerals to grow and to develop its organs, Maize (Zea mays L.) production and distribution is highly rainfall-dependent. Current global climatic changes reveal irregular rainfall patterns and this could represent for maize a stressing condition resulting in yield and productivity loss around the world. It is well known that low water availability leads the plant to adopt a number of metabolic alterations to overcome stress or reduce its effects. In this regard, selenium (Se), a trace element, can help reduce water damage caused by the overproduction of reactive oxygen species (ROS). Here we report the effects of exogenous Se supply on physiological and biochemical processes that may influence yield and quality of maize under drought stress conditions. Plants were grown in soil fertilized by adding 150 mg of Se (sodium selenite). We verified the effects of drought stress and Se treatment. Selenium biofortification proved more beneficial for maize plants when supplied at higher Se concentrations. The increase in proline, K concentrations and nitrogen metabolism in aerial parts of plants grown in Se-rich substrates, seems to prove that Se-biofortification increased plant resistance to water shortage conditions. Moreover, the increase of SeMeSeCys and SeCys2 forms in roots and aerial parts of Se-treated plants suggest resistance strategies to Se similar to those existing in Se-hyperaccumulator species. In addition, epigenetic changes in DNA methylation due to water stress and Se treatment were also investigated using methylation sensitive amplified polymorphism (MSAP). Results suggest that Se may be an activator of particular classes of genes that are involved in tolerance to abiotic stresses. In particular, PSY (phytoene synthase) gene, essential for maintaining leaf carotenoid contents, SDH (sorbitol dehydrogenase), whose activity regulates the level of important osmolytes during drought stress and ADH (alcohol dehydrogenase), whose activity plays a central role in biochemical adaptation to environmental stress. In conclusion, Se-biofortification could help maize plants to cope with drought stress conditions, by inducing a higher drought tolerance.