AlexandrusPS: A User-Friendly Pipeline for the Automated Detection of Orthologous Gene Clusters and Subsequent Positive Selection Analysis.

The detection of adaptive selection in a system approach considering all protein-coding genes allows for the identification of mechanisms and pathways that enabled adaptation to different environments. Currently, available programs for the estimation of positive selection signals can be divided into two groups. They are either easy to apply but can analyze only one gene family at a time, restricting system analysis; or they can handle larger cohorts of gene families, but require considerable prerequisite data such as orthology associations, codon alignments, phylogenetic trees, and proper configuration files. All these steps require extensive computational expertise, restricting this endeavor to specialists. Here, we introduce AlexandrusPS, a high-throughput pipeline that overcomes technical challenges when conducting transcriptome-wide positive selection analyses on large sets of nucleotide and protein sequences. The pipeline streamlines 1) the execution of an accurate orthology prediction as a precondition for positive selection analysis, 2) preparing and organizing configuration files for CodeML, 3) performing positive selection analysis using CodeML, and 4) generating an output that is easy to interpret, including all maximum likelihood and log-likelihood test results. The only input needed from the user is the CDS and peptide FASTA files of proteins of interest. The pipeline is provided in a Docker image, requiring no program or module installation, enabling the application of the pipeline in any computing environment. AlexandrusPS and its documentation are available via GitHub (https://github.com/alejocn5/AlexandrusPS).

DNA damage repair proteins across the Tree of Life.

Genome maintenance is orchestrated by a highly regulated DNA damage response with specific DNA repair pathways. Here, we investigate the phylogenetic diversity in the recognition and repair of three well-established DNA lesions, primarily repaired by base excision repair (BER) and ribonucleotide excision repair (RER): (1) 8-oxoguanine, (2) abasic site, and (3) incorporated ribonucleotide in DNA in 11 species: Escherichia coli, Bacillus subtilis, Halobacterium salinarum, Trypanosoma brucei, Tetrahymena thermophila, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Caenorhabditis elegans, Homo sapiens, Arabidopsis thaliana, and Zea mays. Using quantitative mass spectrometry, we identified 337 binding proteins across these species. Of these proteins, 99 were previously characterized to be involved in DNA repair. Through orthology, network, and domain analysis, we linked 44 previously unconnected proteins to DNA repair. Our study presents a resource for future study of the crosstalk and evolutionary conservation of DNA damage repair across all domains of life.

Rational targeting of a NuRD subcomplex guided by comprehensive in situ mutagenesis.

Developmental silencing of fetal globins serves as both a paradigm of spatiotemporal gene regulation and an opportunity for therapeutic intervention of ß-hemoglobinopathy. The nucleosome remodeling and deacetylase (NuRD) chromatin complex participates in γ-globin repression. We used pooled CRISPR screening to disrupt NuRD protein coding sequences comprehensively in human adult erythroid precursors. Essential for fetal hemoglobin (HbF) control is a non-redundant subcomplex of NuRD protein family paralogs, whose composition we corroborated by affinity chromatography and proximity labeling mass spectrometry proteomics. Mapping top functional guide RNAs identified key protein interfaces where in-frame alleles resulted in loss-of-function due to destabilization or altered function of subunits. We ascertained mutations of CHD4 that dissociate its requirement for cell fitness from HbF repression in both primary human erythroid precursors and transgenic mice. Finally we demonstrated that sequestering CHD4 from NuRD phenocopied these mutations. These results indicate a generalizable approach to discover protein complex features amenable to rational biochemical targeting.

CRISPRO: identification of functional protein coding sequences based on genome editing dense mutagenesis.

CRISPR/Cas9 pooled screening permits parallel evaluation of comprehensive guide RNA libraries to systematically perturb protein coding sequences in situ and correlate with functional readouts. For the analysis and visualization of the resulting datasets, we develop CRISPRO, a computational pipeline that maps functional scores associated with guide RNAs to genomes, transcripts, and protein coordinates and structures. No currently available tool has similar functionality. The ensuing genotype-phenotype linear and three-dimensional maps raise hypotheses about structure-function relationships at discrete protein regions. Machine learning based on CRISPRO features improves prediction of guide RNA efficacy. The CRISPRO tool is freely available at gitlab.com/bauerlab/crispro .

Genome-wide CRISPR-Cas9 Screen Identifies Leukemia-Specific Dependence on a Pre-mRNA Metabolic Pathway Regulated by DCPS.

To identify novel targets for acute myeloid leukemia (AML) therapy, we performed genome-wide CRISPR-Cas9 screening using AML cell lines, followed by a second screen in vivo. Here, we show that the mRNA decapping enzyme scavenger (DCPS) gene is essential for AML cell survival. The DCPS enzyme interacted with components of pre-mRNA metabolic pathways, including spliceosomes, as revealed by mass spectrometry. RG3039, a DCPS inhibitor originally developed to treat spinal muscular atrophy, exhibited anti-leukemic activity via inducing pre-mRNA mis-splicing. Humans harboring germline biallelic DCPS loss-of-function mutations do not exhibit aberrant hematologic phenotypes, indicating that DCPS is dispensable for human hematopoiesis. Our findings shed light on a pre-mRNA metabolic pathway and identify DCPS as a target for AML therapy.

Aberrant hepatic lipid storage and metabolism in canine portosystemic shunts.

Non-alcoholic fatty liver disease (NAFLD) is a poorly understood multifactorial pandemic disorder. One of the hallmarks of NAFLD, hepatic steatosis, is a common feature in canine congenital portosystemic shunts. The aim of this study was to gain detailed insight into the pathogenesis of steatosis in this large animal model. Hepatic lipid accumulation, gene-expression analysis and HPLC-MS of neutral lipids and phospholipids in extrahepatic (EHPSS) and intrahepatic portosystemic shunts (IHPSS) was compared to healthy control dogs. Liver organoids of diseased dogs and healthy control dogs were incubated with palmitic- and oleic-acid, and lipid accumulation was quantified using LD540. In histological slides of shunt livers, a 12-fold increase of lipid content was detected compared to the control dogs (EHPSS P<0.01; IHPSS P = 0.042). Involvement of lipid-related genes to steatosis in portosystemic shunting was corroborated using gene-expression profiling. Lipid analysis demonstrated different triglyceride composition and a shift towards short chain and omega-3 fatty acids in shunt versus healthy dogs, with no difference in lipid species composition between shunt types. All organoids showed a similar increase in triacylglycerols after free fatty acids enrichment. This study demonstrates that steatosis is probably secondary to canine portosystemic shunts. Unravelling the pathogenesis of this hepatic steatosis might contribute to a better understanding of steatosis in NAFLD.