Evolution-based screening enables genome-wide prioritization and discovery of DNA repair genes

Co-evolution of gene transfer agents and their alphaproteobacterial hosts

Gene transfer agents (GTAs) are enigmatic elements that resemble small viruses and are known to be produced during nutritional stress by some bacteria and archaea. The production of GTAs is regulated by quorum sensing, under which a small fraction of the population acts as GTA producers, while the rest becomes GTA recipients. In contrast to canonical viruses, GTAs cannot propagate themselves because they package pieces of the producing cell's genome. In alphaproteobacteria, GTAs are mostly vertically inherited and reside in their hosts' genomes for hundreds of millions of years. While GTAs' ability to transfer genetic material within a population and their long-term preservation suggest an increased fitness of GTA-producing microbes, the associated benefits and type of selection that maintains GTAs are poorly understood. By comparing rates of evolutionary change in GTA genes to the rates in gene families abundantly present across 293 alphaproteobacterial genomes, we detected 59 gene families that likely co-evolve with GTA genes. These gene families are predominantly involved in stress response, DNA repair, and biofilm formation. We hypothesize that biofilm formation enables the physical proximity of GTA-producing cells, limiting GTA-derived benefits only to a group of closely related cells. We further conjecture that the population structure of biofilm-forming sub-populations ensures that the trait of GTA production is maintained despite the inevitable rise of "cheating" genotypes. Because release of GTA particles kills the producing cell, maintenance of GTAs is an exciting example of social evolution in a microbial population.IMPORTANCEGene transfer agents (GTAs) are viruses domesticated by some archaea and bacteria as vehicles for carrying pieces of the host genome. Produced under certain environmental conditions, GTA particles can deliver DNA to neighboring, closely related cells. The function of GTAs remains uncertain. While making GTAs is suicidal for a cell, GTA-encoding genes are widespread in genomes of alphaproteobacteria. Such GTA persistence implies functional benefits but raises questions about how selection maintains this lethal trait. By showing that GTA genes co-evolve with genes involved in stress response, DNA repair, and biofilm formation, we provide support for the hypothesis that GTAs facilitate DNA exchange during the stress conditions and present a model for how GTAs persist in biofilm-forming bacterial populations despite being lethal.

Mechanical stretch promotes apoptosis and impedes ciliogenesis of primary human airway basal stem cells

BACKGROUND

Airway basal stem cells (ABSCs) have self-renewal and differentiation abilities. Although an abnormal mechanical environment related to chronic airway disease (CAD) can cause ABSC dysfunction, it remains unclear how mechanical stretch regulates the behavior and structure of ABSCs. Here, we explored the effect of mechanical stretch on primary human ABSCs.

METHODS

Primary human ABSCs were isolated from healthy volunteers. A Flexcell FX-5000 Tension system was used to mimic the pathological airway mechanical stretch conditions of patients with CAD. ABSCs were stretched for 12, 24, or 48 h with 20% elongation. We first performed bulk RNA sequencing to identify the most predominantly changed genes and pathways. Next, apoptosis of stretched ABSCs was detected with Annexin V-FITC/PI staining and a caspase 3 activity assay. Proliferation of stretched ABSCs was assessed by measuring MKI67 mRNA expression and cell cycle dynamics. Immunofluorescence and hematoxylin-eosin staining were used to demonstrate the differentiation state of ABSCs at the air-liquid interface.

RESULTS

Compared with unstretched control cells, apoptosis and caspase 3 activation of ABSCs stretched for 48 h were significantly increased (p < 0.0001; p < 0.0001, respectively), and MKI67 mRNA levels were decreased (p < 0.0001). In addition, a significant increase in the G0/G1 population (20.2%, p < 0.001) and a significant decrease in S-phase cells (21.1%, p < 0.0001) were observed. The ratio of Krt5+ ABSCs was significantly higher (32.38% vs. 48.71%, p = 0.0037) following stretching, while the ratio of Ac-tub+ cells was significantly lower (37.64% vs. 21.29%, p < 0.001). Moreover, compared with the control, the expression of NKX2-1 was upregulated significantly after stretching (14.06% vs. 39.51%, p < 0.0001). RNA sequencing showed 285 differentially expressed genes, among which 140 were upregulated and 145 were downregulated, revealing that DDIAS, BIRC5, TGFBI, and NKX2-1 may be involved in the function of primary human ABSCs during mechanical stretch. There was no apparent difference between stretching ABSCs for 24 and 48 h compared with the control.

CONCLUSIONS

Pathological stretching induces apoptosis of ABSCs, inhibits their proliferation, and disrupts cilia cell differentiation. These features may be related to abnormal regeneration and repair observed after airway epithelium injury in patients with CAD.

Genomic Instability Evolutionary Footprints on Human Health: Driving Forces or Side Effects?

In this work, we propose a comprehensive perspective on genomic instability comprising not only the accumulation of mutations but also telomeric shortening, epigenetic alterations and other mechanisms that could contribute to genomic information conservation or corruption. First, we present mechanisms playing a role in genomic instability across the kingdoms of life. Then, we explore the impact of genomic instability on the human being across its evolutionary history and on present-day human health, with a particular focus on aging and complex disorders. Finally, we discuss the role of non-coding RNAs, highlighting future approaches for a better living and an expanded healthy lifespan.

DDIAS, DNA damage-induced apoptosis suppressor, is a potential therapeutic target in cancer

Increasing evidence indicates that DNA damage-induced apoptosis suppressor (DDIAS) is an oncogenic protein that is highly expressed in a variety of cancers, including colorectal cancer, lung cancer, breast cancer, and hepatocellular carcinoma (HCC). The discovery of DDIAS as a novel therapeutic target and its role in human cancer biology is fascinating and noteworthy. Recent studies have shown that DDIAS is involved in tumorigenesis, metastasis, DNA repair and synthesis, and drug resistance and that it plays multiple roles with distinct binding partners in several human cancers. This review focuses on the function of DDIAS and its regulatory proteins in human cancer as potential targets for cancer therapy, as well as the development and future prospects of DDIAS inhibitors.

MAIA, Fc receptor-like 3, supersedes JUNO as IZUMO1 receptor during human fertilization

Gamete fusion is a critical event of mammalian fertilization. A random one-bead one-compound combinatorial peptide library represented synthetic human egg mimics and identified a previously unidentified ligand as Fc receptor-like 3, named MAIA after the mythological goddess intertwined with JUNO. This immunoglobulin super family receptor was expressed on human oolemma and played a major role during sperm-egg adhesion and fusion. MAIA forms a highly stable interaction with the known IZUMO1/JUNO sperm-egg complex, permitting specific gamete fusion. The complexity of the MAIA isotype may offer a cryptic sexual selection mechanism to avoid genetic incompatibility and achieve favorable fitness outcomes.

PSPC1 is a potential prognostic marker for hormone-dependent breast cancer patients and modulates RNA processing of ESR1 and SCFD2

Breast cancer is the most common cancer type among women worldwide. The majority of breast cancer expresses estrogen receptor (ER) and endocrine therapy is a standard treatment of ER-positive breast cancer. However, development of the therapy resistance is still a major challenge and thus new therapeutic approaches are needed. Here we show that an RNA-binding protein, PSPC1, play a crucial role in ER-positive breast cancer growth through post-transcriptional gene regulation. We showed that siRNA-mediated PSPC1 silencing suppressed the proliferation of ER-positive breast cancer cells. Strong immunoreactivity (IR) of PSPC1 was correlated with poor prognosis for ER-positive breast cancer patients. Using immunoprecipitation, RNA-immunoprecipitation (RIP) and quantitative PCR (qPCR) experiments, we showed that PSPC1 interacted with PSF and was involved in post-transcriptional regulation of PSF target genes, ESR1 and SCFD2. Strong SCFD2 IR was correlated with poor prognosis for ER-positive breast cancer patients and combinations of PSPC1, PSF, and SCFD2 IRs were potent prognostic factors. Moreover, we identified DDIAS and MYBL1 as SCFD2 downstream target genes using microarray analysis, and finally showed that SCFD2 silencing suppressed tamoxifen-resistant breast tumor growth in vivo. These results indicated that PSPC1 and SCFD2 axis could be a promising target in the clinical management of the disease.

An orthologous gene coevolution network provides insight into eukaryotic cellular and genomic structure and function

The evolutionary rates of functionally related genes often covary. We present a gene coevolution network inferred from examining nearly 3 million orthologous gene pairs from 332 budding yeast species spanning ~400 million years of evolution. Network modules provide insight into cellular and genomic structure and function. Examination of the phenotypic impact of network perturbation using deletion mutant data from the baker's yeast Saccharomyces cerevisiae, which were obtained from previously published studies, suggests that fitness in diverse environments is affected by orthologous gene neighborhood and connectivity. Mapping the network onto the chromosomes of S. cerevisiae and Candida albicans revealed that coevolving orthologous genes are not physically clustered in either species; rather, they are often located on different chromosomes or far apart on the same chromosome. The coevolution network captures the hierarchy of cellular structure and function, provides a roadmap for genotype-to-phenotype discovery, and portrays the genome as a linked ensemble of genes.

Novel ACE2 protein interactions relevant to COVID-19 predicted by evolutionary rate correlations

Angiotensin-converting enzyme 2 (ACE2) is the cell receptor that the coronavirus SARS-CoV-2 binds to and uses to enter and infect human cells. COVID-19, the pandemic disease caused by the coronavirus, involves diverse pathologies beyond those of a respiratory disease, including micro-thrombosis (micro-clotting), cytokine storms, and inflammatory responses affecting many organ systems. Longer-term chronic illness can persist for many months, often well after the pathogen is no longer detected. A better understanding of the proteins that ACE2 interacts with can reveal information relevant to these disease manifestations and possible avenues for treatment. We have undertaken an approach to predict candidate ACE2 interacting proteins which uses evolutionary inference to identify a set of mammalian proteins that “coevolve” with ACE2. The approach, called evolutionary rate correlation (ERC), detects proteins that show highly correlated evolutionary rates during mammalian evolution. Such proteins are candidates for biological interactions with the ACE2 receptor. The approach has uncovered a number of key ACE2 protein interactions of potential relevance to COVID-19 pathologies. Some proteins have previously been reported to be associated with severe COVID-19, but are not currently known to interact with ACE2, while additional predicted novel ACE2 interactors are of potential relevance to the disease. Using reciprocal rankings of protein ERCs, we have identified strongly interconnected ACE2 associated protein networks relevant to COVID-19 pathologies. ACE2 has clear connections to coagulation pathway proteins, such as Coagulation Factor V and fibrinogen components FGA, FGB, and FGG, the latter possibly mediated through ACE2 connections to Clusterin (which clears misfolded extracellular proteins) and GPR141 (whose functions are relatively unknown). ACE2 also connects to proteins involved in cytokine signaling and immune response (e.g. XCR1, IFNAR2 and TLR8), and to Androgen Receptor (AR). The ERC prescreening approach has elucidated possible functions for relatively uncharacterized proteins and possible new functions for well-characterized ones. Suggestions are made for the validation of ERC-predicted ACE2 protein interactions. We propose that ACE2 has novel protein interactions that are disrupted during SARS-CoV-2 infection, contributing to the spectrum of COVID-19 pathologies.

Evolutionary rate covariation identifies SLC30A9 (ZnT9) as a mitochondrial zinc transporter

Recent advances in genome sequencing have led to the identification of new ion and metabolite transporters, many of which have not been characterized. Due to the variety of subcellular localizations, cargo and transport mechanisms, such characterization is a daunting task, and predictive approaches focused on the functional context of transporters are very much needed. Here we present a case for identifying a transporter localization using evolutionary rate covariation (ERC), a computational approach based on pairwise correlations of amino acid sequence evolutionary rates across the mammalian phylogeny. As a case study, we find that poorly characterized transporter SLC30A9 (ZnT9) coevolves with several components of the mitochondrial oxidative phosphorylation chain, suggesting mitochondrial localization. We confirmed this computational finding experimentally using recombinant human SLC30A9. SLC30A9 loss caused zinc mishandling in the mitochondria, suggesting that under normal conditions it acts as a zinc exporter. We therefore propose that ERC can be used to predict the functional context of novel transporters and other poorly characterized proteins.

PhyKIT: a broadly applicable UNIX shell toolkit for processing and analyzing phylogenomic data

MOTIVATION

Diverse disciplines in biology process and analyze multiple sequence alignments (MSAs) and phylogenetic trees to evaluate their information content, infer evolutionary events and processes, and predict gene function. However, automated processing of MSAs and trees remains a challenge due to the lack of a unified toolkit. To fill this gap, we introduce PhyKIT, a toolkit for the UNIX shell environment with 30 functions that process MSAs and trees, including but not limited to estimation of mutation rate, evaluation of sequence composition biases, calculation of the degree of violation of a molecular clock, and collapsing bipartitions (internal branches) with low support.

RESULTS

To demonstrate the utility of PhyKIT, we detail three use cases: (1) summarizing information content in MSAs and phylogenetic trees for diagnosing potential biases in sequence or tree data; (2) evaluating gene-gene covariation of evolutionary rates to identify functional relationships, including novel ones, among genes; and (3) identify lack of resolution events or polytomies in phylogenetic trees, which are suggestive of rapid radiation events or lack of data. We anticipate PhyKIT will be useful for processing, examining, and deriving biological meaning from increasingly large phylogenomic datasets.

AVAILABILITY

PhyKIT is freely available on GitHub (https://github.com/JLSteenwyk/PhyKIT), PyPi (https://pypi.org/project/phykit/), and the Anaconda Cloud (https://anaconda.org/JLSteenwyk/phykit) under the MIT license with extensive documentation and user tutorials (https://jlsteenwyk.com/PhyKIT).

SUPPLEMENTARY INFORMATION

Supplementary data are available on figshare (doi: 10.6084/m9.figshare.13118600) and are available at Bioinformatics online.

Regulation of ddb2 expression in blind cavefish and zebrafish reveals plasticity in the control of sunlight-induced DNA damage repair

We have gained considerable insight into the mechanisms which recognize and repair DNA damage, but how they adapt to extreme environmental challenges remains poorly understood. Cavefish have proven to be fascinating models for exploring the evolution of DNA repair in the complete absence of UV-induced DNA damage and light. We have previously revealed that the Somalian cavefish Phreatichthys andruzzii, lacks photoreactivation repair via the loss of light, UV and ROS-induced photolyase gene transcription mediated by D-box enhancer elements. Here, we explore whether other systems repairing UV-induced DNA damage have been similarly affected in this cavefish model. By performing a comparative study using P. andruzzii and the surface-dwelling zebrafish, we provide evidence for a conservation of sunlight-regulated Nucleotide Excision Repair (NER). Specifically, the expression of the ddb2 gene which encodes a key NER recognition factor is robustly induced following exposure to light, UV and oxidative stress in both species. As in the case of the photolyase genes, D-boxes in the ddb2 promoter are sufficient to induce transcription in zebrafish. Interestingly, despite the loss of D-box-regulated photolyase gene expression in P. andruzzii, the D-box is required for ddb2 induction by visible light and oxidative stress in cavefish. However, in the cavefish ddb2 gene this D-box-mediated induction requires cooperation with an adjacent, highly conserved E2F element. Furthermore, while in zebrafish UV-induced ddb2 expression results from transcriptional activation accompanied by stabilization of the ddb2 mRNA, in P. andruzzii UV induces ddb2 expression exclusively via an increase in mRNA stability. Thus, we reveal plasticity in the transcriptional and post transcriptional mechanisms regulating the repair of sunlight-induced DNA damage under long-term environmental challenges.

Formation and repair of unavoidable, endogenous interstrand cross-links in cellular DNA

Genome integrity is essential for life and, as a result, DNA repair systems evolved to remove unavoidable DNA lesions from cellular DNA. Many forms of life possess the capacity to remove interstrand DNA cross-links (ICLs) from their genome but the identity of the naturally-occurring, endogenous substrates that drove the evolution and retention of these DNA repair systems across a wide range of life forms remains uncertain. In this review, we describe more than a dozen chemical processes by which endogenous ICLs plausibly can be introduced into cellular DNA. The majority involve DNA degradation processes that introduce aldehyde residues into the double helix or reactions of DNA with endogenous low molecular weight aldehyde metabolites. A smaller number of the cross-linking processes involve reactions of DNA radicals generated by oxidation.

A Drosophila screen identifies NKCC1 as a modifier of NGLY1 deficiency

N-Glycanase 1 (NGLY1) is a cytoplasmic deglycosylating enzyme. Loss-of-function mutations in the NGLY1 gene cause NGLY1 deficiency, which is characterized by developmental delay, seizures, and a lack of sweat and tears. To model the phenotypic variability observed among patients, we crossed a Drosophila model of NGLY1 deficiency onto a panel of genetically diverse strains. The resulting progeny showed a phenotypic spectrum from 0 to 100% lethality. Association analysis on the lethality phenotype, as well as an evolutionary rate covariation analysis, generated lists of modifying genes, providing insight into NGLY1 function and disease. The top association hit was Ncc69 (human NKCC1/2), a conserved ion transporter. Analyses in NGLY1-/- mouse cells demonstrated that NKCC1 has an altered average molecular weight and reduced function. The misregulation of this ion transporter may explain the observed defects in secretory epithelium function in NGLY1 deficiency patients.

MCM8IP activates the MCM8-9 helicase to promote DNA synthesis and homologous recombination upon DNA damage

Homologous recombination (HR) mediates the error-free repair of DNA double-strand breaks to maintain genomic stability. Here we characterize C17orf53/MCM8IP, an OB-fold containing protein that binds ssDNA, as a DNA repair factor involved in HR. MCM8IP-deficient cells exhibit HR defects, especially in long-tract gene conversion, occurring downstream of RAD51 loading, consistent with a role for MCM8IP in HR-dependent DNA synthesis. Moreover, loss of MCM8IP confers cellular sensitivity to crosslinking agents and PARP inhibition. Importantly, we report that MCM8IP directly associates with MCM8-9, a helicase complex mutated in primary ovarian insufficiency, and RPA1. We additionally show that the interactions of MCM8IP with MCM8-9 and RPA facilitate HR and promote replication fork progression and cellular viability in response to treatment with crosslinking agents. Mechanistically, MCM8IP stimulates the helicase activity of MCM8-9. Collectively, our work identifies MCM8IP as a key regulator of MCM8-9-dependent DNA synthesis during DNA recombination and replication.