CDK12 loss drives prostate cancer progression, transcription-replication conflicts, and synthetic lethality with paralog CDK13.

Biallelic loss of cyclin-dependent kinase 12 (CDK12) defines a metastatic castration-resistant prostate cancer (mCRPC) subtype. It remains unclear, however, whether CDK12 loss drives prostate cancer (PCa) development or uncovers pharmacologic vulnerabilities. Here, we show Cdk12 ablation in murine prostate epithelium is sufficient to induce preneoplastic lesions with lymphocytic infiltration. In allograft-based CRISPR screening, Cdk12 loss associates positively with Trp53 inactivation but negatively with Pten inactivation. Moreover, concurrent Cdk12/Trp53 ablation promotes proliferation of prostate-derived organoids, while Cdk12 knockout in Pten-null mice abrogates prostate tumor growth. In syngeneic systems, Cdk12/Trp53-null allografts exhibit luminal morphology and immune checkpoint blockade sensitivity. Mechanistically, Cdk12 inactivation mediates genomic instability by inducing transcription-replication conflicts. Strikingly, CDK12-mutant organoids and patient-derived xenografts are sensitive to inhibition or degradation of the paralog kinase, CDK13. We therein establish CDK12 as a bona fide tumor suppressor, mechanistically define how CDK12 inactivation causes genomic instability, and advance a therapeutic strategy for CDK12-mutant mCRPC.

MpPUB9, a U-box E3 ubiquitin ligase, acts as a positive regulator by promoting the turnover of MpEXO70.1 under high salinity in Marchantia polymorpha.

Marchantia polymorpha, occupying a basal position in the monophyletic assemblage of land plants, displays a notable expansion of plant U-box (PUB) proteins compared with those in animals. We elucidated the roles of MpPUB9 in regulating salt stress tolerance in M. polymorpha. MpPUB9 expression was rapidly induced by high salinity and dehydration. MpPUB9 possessed an intact U-box domain in the N-terminus. MpPUB9-Citrine localized to punctate structures and was peripherally associated with microsomal membranes. Phenotypic analyses demonstrate that the hyponastic and epinastic thallus growth phenotypes, which were induced by the overexpression and suppression of MpPUB9, may provoke salt stress-resistant and -susceptible phenotypes, respectively. MpPUB9 was also found to directly interact with the exocyst protein MpEXO70.1, leading to its ubiquitination. Under high-salinity conditions, though the stability of MpPUB9 was dramatically increased, MpEXO70.1 showed slightly faster turnover rates. Transcriptome analyses showed that salt treatment and the overexpression of MpPUB9 co-upregulated the genes related to the modulation of H2O2 and cell wall organization. Overall, our results suggest that MpPUB9 plays a crucial role in the positive regulation of salt stress tolerance, resulting from its interaction with MpEXO70.1 and modulating turnover of the protein under high-salt conditions via the coordination of UPS with autophagy.

Evolution and Spatiotemporal Expression of ankha and ankhb in Zebrafish.

Craniometaphyseal Dysplasia (CMD) is a rare skeletal disorder that can result from mutations in the ANKH gene. This gene encodes progressive anksylosis (ANK), which is responsible for transporting inorganic pyrophosphate (PPi) and ATP from the intracellular to the extracellular environment, where PPi inhibits bone mineralization. When ANK is dysfunctional, as in patients with CMD, the passage of PPi to the extracellular environment is reduced, leading to excess mineralization, particularly in bones of the skull. Zebrafish may serve as a promising model to study the mechanistic basis of CMD. Here, we provide a detailed analysis of the zebrafish Ankh paralogs, Ankha and Ankhb, in terms of their phylogenic relationship with ANK in other vertebrates as well as their spatiotemporal expression patterns during zebrafish development. We found that a closer evolutionary relationship exists between the zebrafish Ankhb protein and its human and other vertebrate counterparts, and stronger promoter activity was predicted for ankhb compared to ankha. Furthermore, we noted distinct temporal expression patterns, with ankha more prominently expressed in early development stages, and both paralogs also being expressed at larval growth stages. Whole-mount in situ hybridization was used to compare the spatial expression patterns of each paralog during bone development, and both showed strong expression in the craniofacial region as well as the notochord and somites. Given the substantial overlap in spatiotemporal expression but only subtle patterning differences, the exact roles of these genes remain speculative. In silico analyses predicted that Ankha and Ankhb have the same function in transporting PPi across the membrane. Nevertheless, this study lays the groundwork for functional analyses of each ankh paralog and highlights the potential of using zebrafish to find possible targeted therapies for CMD.

GABARAP interacts with EGFR - supporting the unique role of this hAtg8 protein during receptor trafficking.

The human Atg8 family member GABARAP is involved in numerous autophagy-related and -unrelated processes. We recently observed that specifically the deficiency of GABARAP enhances epidermal growth factor receptor (EGFR) degradation upon ligand stimulation. Here, we report on two putative LC3-interacting regions (LIRs) within EGFR, the first of which (LIR1) is selected as a GABARAP binding site in silico. Indeed, in vitro interaction studies reveal preferential binding of LIR1 to GABARAP and GABARAPL1. Our X-ray data demonstrate interaction of core LIR1 residues FLPV with both hydrophobic pockets of GABARAP suggesting canonical binding. Although LIR1 occupies the LIR docking site, GABARAP Y49 and L50 appear dispensable in this case. Our data support the hypothesis that GABARAP affects the fate of EGFR at least in part through direct binding.

Varying Selection Pressure for a Na+ Sensing Site in Epithelial Na+ Channel Subunits Reflect Divergent Roles in Na+ Homeostasis.

The epithelial Na+ channel (ENaC) emerged early in vertebrates and has played a role in Na+ and fluid homeostasis throughout vertebrate evolution. We previously showed that proteolytic activation of the channel evolved at the water-to-land transition of vertebrates. Sensitivity to extracellular Na+, known as Na+ self-inhibition, reduces ENaC function when Na+ concentrations are high and is a distinctive feature of the channel. A fourth ENaC subunit, δ, emerged in jawed fishes from an α subunit gene duplication. Here, we analyzed 849 α and δ subunit sequences and found that a key Asp in a postulated Na+ binding site was nearly always present in the α subunit, but frequently lost in the δ subunit (e.g. human). Analysis of site evolution and codon substitution rates provide evidence that the ancestral α subunit had the site and that purifying selection for the site relaxed in the δ subunit after its divergence from the α subunit, coinciding with a loss of δ subunit expression in renal tissues. We also show that the proposed Na+ binding site in the α subunit is a bona fide site by conferring novel function to channels comprising human δ subunits. Together, our findings provide evidence that ENaC Na+ self-inhibition improves fitness through its role in Na+ homeostasis in vertebrates.

Paralog switching facilitates diadromy: ontogenetic, microevolutionary and macroevolutionary evidence.

Identifying how the demands of migration are met at the level of gene expression is critical for understanding migratory physiology and can potentially reveal how migratory forms evolve from nonmigratory forms and vice versa. Among fishes, migration between freshwater and seawater (diadromy) requires considerable osmoregulatory adjustments, powered by the ion pump Na+, K+-ATPase (NKA) in the gills. Paralogs of the catalytic α-subunit of the pump (NKA α1a and α1b) are reciprocally upregulated in fresh- and seawater, a response known as paralog-switching, in gills of some diadromous species. We tested ontogenetic changes in NKA α-subunit paralog expression patterns, comparing pre-migrant and migrant alewife (Alosa pseudoharengus) sampled in their natal freshwater environment and after 24 h in seawater. In comparison to pre-migrants, juvenile out-migrants exhibited stronger paralog switching via greater downregulation of NKA α1a in seawater. We also tested microevolutionary changes in the response, exposing juvenile diadromous and landlocked alewife to freshwater (0 ppt) and seawater (30 ppt) for 2, 5, and 15 days. Diadromous and landlocked alewife exhibited salinity-dependent paralog switching, but levels of NKA α1b transcription were higher and the decrease in NKA α1a was greater after seawater exposure in diadromous alewife. Finally, we placed alewife α-subunit NKA paralogs in a macroevolutionary context. Molecular phylogenies show alewife paralogs originated independently of paralogs in salmonids and other teleosts. This study demonstrated that NKA paralog switching is tied to halohabitat profile and that duplications of the NKA gene provided the substrate for multiple, independent molecular solutions that support a diadromous life history.

Lhx4 surpasses its paralog Lhx3 in promoting the differentiation of spinal V2a interneurons.

Paralog factors are considered to ensure the robustness of biological processes by providing redundant activity in cells where they are co-expressed. However, the specific contribution of each factor is frequently underestimated. In the developing spinal cord, multiple families of transcription factors successively contribute to differentiate an initially homogenous population of neural progenitors into a myriad of neuronal subsets with distinct molecular, morphological, and functional characteristics. The LIM-homeodomain transcription factors Lhx3, Lhx4, Isl1 and Isl2 promote the segregation and differentiation of spinal motor neurons and V2 interneurons. Based on their high sequence identity and their similar distribution, the Lhx3 and Lhx4 paralogs are considered to contribute similarly to these processes. However, the specific contribution of Lhx4 has never been studied. Here, we provide evidence that Lhx3 and Lhx4 are present in the same cell populations during spinal cord development. Similarly to Lhx3, Lhx4 can form multiproteic complexes with Isl1 or Isl2 and the nuclear LIM interactor NLI. Lhx4 can stimulate a V2-specific enhancer more efficiently than Lhx3 and surpasses Lhx3 in promoting the differentiation of V2a interneurons in chicken embryo electroporation experiments. Finally, Lhx4 inactivation in mice results in alterations of differentiation of the V2a subpopulation, but not of motor neuron production, suggesting that Lhx4 plays unique roles in V2a differentiation that are not compensated by the presence of Lhx3. Thus, Lhx4 could be the major LIM-HD factor involved in V2a interneuron differentiation during spinal cord development and should be considered for in vitro differentiation of spinal neuronal populations.

A novel chemical genetic approach reveals paralog-specific role of ERK1/2 in mouse embryonic stem cell fate control.

Introduction: Mouse embryonic stem cell (ESC) self-renewal can be maintained through dual inhibition of GSK3 and MEK kinases. MEK has two highly homologous downstream kinases, extracellular signal-regulated kinase 1 and 2 (ERK1/2). However, the exact roles of ERK1/2 in mouse ESC self-renewal and differentiation remain unclear. Methods: We selectively deleted or inhibited ERK1, ERK2, or both using genetic and chemical genetic approaches combined with small molecule inhibitors. The effects of ERK paralog-specific inhibition on mouse ESC self-renewal and differentiation were then assessed. Results: ERK1/2 were found to be dispensable for mouse ESC survival and self-renewal. The inhibition of both ERK paralogs, in conjunction with GSK3 inhibition, was sufficient to maintain mouse ESC self-renewal. In contrast, selective deletion or inhibition of only one ERK paralog did not mimic the effect of MEK inhibition in promoting mouse ESC self-renewal. Regarding ESC differentiation, inhibition of ERK1/2 prevented mesendoderm differentiation. Additionally, selective inhibition of ERK1, but not ERK2, promoted mesendoderm differentiation. Discussion: These findings suggest that ERK1 and ERK2 have both overlapping and distinct roles in regulating ESC self-renewal and differentiation. This study provides new insights into the molecular mechanisms of ERK1/2 in governing ESC maintenance and lineage commitment, potentially informing future strategies for controlling stem cell fate in research and therapeutic applications.

The TSC22D, WNK, and NRBP gene families exhibit functional buffering and evolved with Metazoa for cell volume regulation.

The ability to sense and respond to osmotic fluctuations is critical for the maintenance of cellular integrity. We used gene co-essentiality analysis to identify an unappreciated relationship between TSC22D2, WNK1, and NRBP1 in regulating cell volume homeostasis. All of these genes have paralogs and are functionally buffered for osmo-sensing and cell volume control. Within seconds of hyperosmotic stress, TSC22D, WNK, and NRBP family members physically associate into biomolecular condensates, a process that is dependent on intrinsically disordered regions (IDRs). A close examination of these protein families across metazoans revealed that TSC22D genes evolved alongside a domain in NRBPs that specifically binds to TSC22D proteins, which we have termed NbrT (NRBP binding region with TSC22D), and this co-evolution is accompanied by rapid IDR length expansion in WNK-family kinases. Our study reveals that TSC22D, WNK, and NRBP genes evolved in metazoans to co-regulate rapid cell volume changes in response to osmolarity.

Mitoferrin2 is a synthetic lethal target for chromosome 8p deleted cancers.

BACKGROUND: Somatic copy number alterations are a hallmark of cancer that offer unique opportunities for therapeutic exploitation. Here, we focused on the identification of specific vulnerabilities for tumors harboring chromosome 8p deletions. METHODS: We developed and applied an integrative analysis of The Cancer Genome Atlas (TCGA), the Cancer Dependency Map (DepMap), and the Cancer Cell Line Encyclopedia to identify chromosome 8p-specific vulnerabilities. We employ orthogonal gene targeting strategies, both in vitro and in vivo, including short hairpin RNA-mediated gene knockdown and CRISPR/Cas9-mediated gene knockout to validate vulnerabilities. RESULTS: We identified SLC25A28 (also known as MFRN2), as a specific vulnerability for tumors harboring chromosome 8p deletions. We demonstrate that vulnerability towards MFRN2 loss is dictated by the expression of its paralog, SLC25A37 (also known as MFRN1), which resides on chromosome 8p. In line with their function as mitochondrial iron transporters, MFRN1/2 paralog protein deficiency profoundly impaired mitochondrial respiration, induced global depletion of iron-sulfur cluster proteins, and resulted in DNA-damage and cell death. MFRN2 depletion in MFRN1-deficient tumors led to impaired growth and even tumor eradication in preclinical mouse xenograft experiments, highlighting its therapeutic potential. CONCLUSIONS: Our data reveal MFRN2 as a therapeutic target of chromosome 8p deleted cancers and nominate MFNR1 as the complimentary biomarker for MFRN2-directed therapies.

OrthoRefine: automated enhancement of prior ortholog identification via synteny.

BACKGROUND: Identifying orthologs continues to be an early and imperative step in genome analysis but remains a challenging problem. While synteny (conservation of gene order) has previously been used independently and in combination with other methods to identify orthologs, applying synteny in ortholog identification has yet to be automated in a user-friendly manner. This desire for automation and ease-of-use led us to develop OrthoRefine, a standalone program that uses synteny to refine ortholog identification. RESULTS: We developed OrthoRefine to improve the detection of orthologous genes by implementing a look-around window approach to detect synteny. We tested OrthoRefine in tandem with OrthoFinder, one of the most used software for identification of orthologs in recent years. We evaluated improvements provided by OrthoRefine in several bacterial and a eukaryotic dataset. OrthoRefine efficiently eliminates paralogs from orthologous groups detected by OrthoFinder. Using synteny increased specificity and functional ortholog identification; additionally, analysis of BLAST e-value, phylogenetics, and operon occurrence further supported using synteny for ortholog identification. A comparison of several window sizes suggested that smaller window sizes (eight genes) were generally the most suitable for identifying orthologs via synteny. However, larger windows (30 genes) performed better in datasets containing less closely related genomes. A typical run of OrthoRefine with ~ 10 bacterial genomes can be completed in a few minutes on a regular desktop PC. CONCLUSION: OrthoRefine is a simple-to-use, standalone tool that automates the application of synteny to improve ortholog detection. OrthoRefine is particularly efficient in eliminating paralogs from orthologous groups delineated by standard methods.

CDK12 Loss Promotes Prostate Cancer Development While Exposing Vulnerabilities to Paralog-Based Synthetic Lethality.

Biallelic loss of cyclin-dependent kinase 12 (CDK12) defines a unique molecular subtype of metastatic castration-resistant prostate cancer (mCRPC). It remains unclear, however, whether CDK12 loss per se is sufficient to drive prostate cancer development-either alone, or in the context of other genetic alterations-and whether CDK12-mutant tumors exhibit sensitivity to specific pharmacotherapies. Here, we demonstrate that tissue-specific Cdk12 ablation is sufficient to induce preneoplastic lesions and robust T cell infiltration in the mouse prostate. Allograft-based CRISPR screening demonstrated that Cdk12 loss is positively associated with Trp53 inactivation but negatively associated with Pten inactivation-akin to what is observed in human mCRPC. Consistent with this, ablation of Cdk12 in prostate organoids with concurrent Trp53 loss promotes their proliferation and ability to form tumors in mice, while Cdk12 knockout in the Pten-null prostate cancer mouse model abrogates tumor growth. Bigenic Cdk12 and Trp53 loss allografts represent a new syngeneic model for the study of androgen receptor (AR)-positive, luminal prostate cancer. Notably, Cdk12/Trp53 loss prostate tumors are sensitive to immune checkpoint blockade. Cdk12-null organoids (either with or without Trp53 co-ablation) and patient-derived xenografts from tumors with CDK12 inactivation are highly sensitive to inhibition or degradation of its paralog kinase, CDK13. Together, these data identify CDK12 as a bona fide tumor suppressor gene with impact on tumor progression and lends support to paralog-based synthetic lethality as a promising strategy for treating CDK12-mutant mCRPC.

NSD2 is a requisite subunit of the AR/FOXA1 neo-enhanceosome in promoting prostate tumorigenesis.

The androgen receptor (AR) is a ligand-responsive transcription factor that binds at enhancers to drive terminal differentiation of the prostatic luminal epithelia. By contrast, in tumors originating from these cells, AR chromatin occupancy is extensively reprogrammed to drive hyper-proliferative, metastatic, or therapy-resistant phenotypes, the molecular mechanisms of which remain poorly understood. Here, we show that the tumor-specific enhancer circuitry of AR is critically reliant on the activity of Nuclear Receptor Binding SET Domain Protein 2 (NSD2), a histone 3 lysine 36 di-methyltransferase. NSD2 expression is abnormally gained in prostate cancer cells and its functional inhibition impairs AR trans-activation potential through partial off-loading from over 40,000 genomic sites, which is greater than 65% of the AR tumor cistrome. The NSD2-dependent AR sites distinctly harbor a chimeric AR-half motif juxtaposed to a FOXA1 element. Similar chimeric motifs of AR are absent at the NSD2-independent AR enhancers and instead contain the canonical palindromic motifs. Meta-analyses of AR cistromes from patient tumors uncovered chimeric AR motifs to exclusively participate in tumor-specific enhancer circuitries, with a minimal role in the physiological activity of AR. Accordingly, NSD2 inactivation attenuated hallmark cancer phenotypes that were fully reinstated upon exogenous NSD2 re-expression. Inactivation of NSD2 also engendered increased dependency on its paralog NSD1, which independently maintained AR and MYC hyper-transcriptional programs in cancer cells. Concordantly, a dual NSD1/2 PROTAC degrader, called LLC0150, was preferentially cytotoxic in AR-dependent prostate cancer as well as NSD2-altered hematologic malignancies. Altogether, we identify NSD2 as a novel subunit of the AR neo-enhanceosome that wires prostate cancer gene expression programs, positioning NSD1/2 as viable paralog co-targets in advanced prostate cancer.

Functional Divergence in Solute Permeability between Ray-Finned Fish-Specific Paralogs of aqp10.

Aquaporin (Aqp) 10 is a member of the aquaglyceroporin subfamily of water channels, and human Aqp10 is permeable to solutes such as glycerol, urea, and boric acid. Tetrapods have a single aqp10 gene, whereas ray-finned fishes have paralogs of this gene through tandem duplication, whole-genome duplication, and subsequent deletion. A previous study on Aqps in the Japanese pufferfish Takifugu rubripes showed that one pufferfish paralog, Aqp10.2b, was permeable to water and glycerol, but not to urea and boric acid. To understand the functional differences of Aqp10s between humans and pufferfish from an evolutionary perspective, we analyzed Aqp10s from an amphibian (Xenopus laevis) and a lobe-finned fish (Protopterus annectens) and Aqp10.1 and Aqp10.2 from several ray-finned fishes (Polypterus senegalus, Lepisosteus oculatus, Danio rerio, and Clupea pallasii). The expression of tetrapod and lobe-finned fish Aqp10s and Aqp10.1-derived Aqps in ray-finned fishes in Xenopus oocytes increased the membrane permeabilities to water, glycerol, urea, and boric acid. In contrast, Aqp10.2-derived Aqps in ray-finned fishes increased water and glycerol permeabilities, whereas those of urea and boric acid were much weaker than those of Aqp10.1-derived Aqps. These results indicate that water, glycerol, urea, and boric acid permeabilities are plesiomorphic activities of Aqp10s and that the ray-finned fish-specific Aqp10.2 paralogs have secondarily reduced or lost urea and boric acid permeability.

Widespread Deviant Patterns of Heterozygosity in Whole-Genome Sequencing Due to Autopolyploidy, Repeated Elements, and Duplication.

Most population genomic tools rely on accurate single nucleotide polymorphism (SNP) calling and filtering to meet their underlying assumptions. However, genomic complexity, resulting from structural variants, paralogous sequences, and repetitive elements, presents significant challenges in assembling contiguous reference genomes. Consequently, short-read resequencing studies can encounter mismapping issues, leading to SNPs that deviate from Mendelian expected patterns of heterozygosity and allelic ratio. In this study, we employed the ngsParalog software to identify such deviant SNPs in whole-genome sequencing (WGS) data with low (1.5×) to intermediate (4.8×) coverage for four species: Arctic Char (Salvelinus alpinus), Lake Whitefish (Coregonus clupeaformis), Atlantic Salmon (Salmo salar), and the American Eel (Anguilla rostrata). The analyses revealed that deviant SNPs accounted for 22% to 62% of all SNPs in salmonid datasets and approximately 11% in the American Eel dataset. These deviant SNPs were particularly concentrated within repetitive elements and genomic regions that had recently undergone rediploidization in salmonids. Additionally, narrow peaks of elevated coverage were ubiquitous along all four reference genomes, encompassed most deviant SNPs, and could be partially associated with transposons and tandem repeats. Including these deviant SNPs in genomic analyses led to highly distorted site frequency spectra, underestimated pairwise FST values, and overestimated nucleotide diversity. Considering the widespread occurrence of deviant SNPs arising from a variety of sources, their important impact in estimating population parameters, and the availability of effective tools to identify them, we propose that excluding deviant SNPs from WGS datasets is required to improve genomic inferences for a wide range of taxa and sequencing depths.

Deep computational phenotyping of genomic variants impacting the SET domain of KMT2C reveal molecular mechanisms for their dysfunction.

Introduction: Kleefstra Syndrome type 2 (KLEFS-2) is a genetic, neurodevelopmental disorder characterized by intellectual disability, infantile hypotonia, severe expressive language delay, and characteristic facial appearance, with a spectrum of other distinct clinical manifestations. Pathogenic mutations in the epigenetic modifier type 2 lysine methyltransferase KMT2C have been identified to be causative in KLEFS-2 individuals. Methods: This work reports a translational genomic study that applies a multidimensional computational approach for deep variant phenotyping, combining conventional genomic analyses, advanced protein bioinformatics, computational biophysics, biochemistry, and biostatistics-based modeling. We use standard variant annotation, paralog annotation analyses, molecular mechanics, and molecular dynamics simulations to evaluate damaging scores and provide potential mechanisms underlying KMT2C variant dysfunction. Results: We integrated data derived from the structure and dynamics of KMT2C to classify variants into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). When compared with controls, these variants show values reflecting alterations in molecular fitness in both structure and dynamics. Discussion: We demonstrate that our 3D models for KMT2C variants suggest distinct mechanisms that lead to their imbalance and are not predictable from sequence alone. Thus, the missense variants studied here cause destabilizing effects on KMT2C function by different biophysical and biochemical mechanisms which we adeptly describe. This new knowledge extends our understanding of how variations in the KMT2C gene cause the dysfunction of its methyltransferase enzyme product, thereby bearing significant biomedical relevance for carriers of KLEFS2-associated genomic mutations.

FgSdhC Paralog Confers Natural Resistance toward SDHI Fungicides in Fusarium graminearum.

Fusarium graminearum exhibited natural resistance to a majority of succinate dehydrogenase inhibitor fungicides (SDHIs) and the molecular mechanisms responsible for the natural resistance were still unknown. Succinate dehydrogenase subunit C (SdhC) is an essential gene for maintaining succinate-ubiquinone oxidoreductase (SQR) function in fungi. In F. graminearum, a paralog of FgSdhC named as FgSdhC1 was identified. Based on RNA-Seq and qRT-PCR assay, we found that the expression level of FgSdhC1 was very low but upregulated by SDHIs treatment. Based on reverse genetics, we demonstrated that FgSdhC1 was an inessential gene in normal growth but was sufficient for maintaining SQR function and conferred natural resistance or reduced sensitivity toward SDHIs. Additionally, we found that the standard F. graminearum isolate PH-1 had high sensitivity to a majority of SDHIs. A single nucleotide variation (C to T) in the FgSdhC1 of isolate PH-1, resulting in a premature termination codon (TAA) replacing the fourth amino acid glutamine (Q), led to the failure of FgSdhC1 to perform functions of conferring nature resistance. These results established that a dispensable paralogous gene determined SDHIs resistance in natural populations of F. graminearum.

Cardiolipin occupancy profiles of YidC paralogs reveal the significance of respective TM2 helix residues in determining paralog-specific phenotypes.

YidC belongs to an evolutionarily conserved family of insertases, YidC/Oxa1/Alb3, in bacteria, mitochondria, and chloroplasts, respectively. Unlike Gram-negative bacteria, Gram-positives including Streptococcus mutans harbor two paralogs of YidC. The mechanism for paralog-specific phenotypes of bacterial YidC1 versus YidC2 has been partially attributed to the differences in their cytoplasmic domains. However, we previously identified a W138R gain-of-function mutation in the YidC1 transmembrane helix 2. YidC1W138R mostly phenocopied YidC2, yet the mechanism remained unknown. Primary sequence comparison of streptococcal YidCs led us to identify and mutate the YidC1W138 analog, YidC2S152 to W/A, which resulted in a loss of YidC2- and acquisition of YidC1-like phenotype. The predicted lipid-facing side chains of YidC1W138/YidC2S152 led us to propose a role for membrane phospholipids in specific-residue dependent phenotypes of S. mutans YidC paralogs. Cardiolipin (CL), a prevalent phospholipid in the S. mutans cytoplasmic membrane during acid stress, is encoded by a single gene, cls. We show a concerted mechanism for cardiolipin and YidC2 under acid stress based on similarly increased promoter activities and similar elimination phenotypes. Using coarse grain molecular dynamics simulations with the Martini2.2 Forcefield, YidC1 and YidC2 wild-type and mutant interactions with CL were assessed in silico. We observed substantially increased CL interaction in dimeric versus monomeric proteins, and variable CL occupancy in YidC1 and YidC2 mutant constructs that mimicked characteristics of the other wild-type paralog. Hence, paralog-specific amino acid- CL interactions contribute to YidC1 and YidC2-associated phenotypes that can be exchanged by point mutation at positions 138 or 152, respectively.

Tackling redundancy: genetic mechanisms underlying paralog compensation in plants.

Gene duplication is a powerful source of biological innovation giving rise to paralogous genes that undergo diverse fates. Redundancy between paralogous genes is an intriguing outcome of duplicate gene evolution, and its maintenance over evolutionary time has long been considered a paradox. Redundancy can also be dubbed 'a geneticist's nightmare': It hinders the predictability of genome editing outcomes and limits our ability to link genotypes to phenotypes. Genetic studies in yeast and plants have suggested that the ability of ancient redundant duplicates to compensate for dosage perturbations resulting from a loss of function depends on the reprogramming of gene expression, a phenomenon known as active compensation. Starting from considerations on the stoichiometric constraints that drive the evolutionary stability of redundancy, this review aims to provide insights into the mechanisms of active compensation between duplicates that could be targeted for breaking paralog dependencies - the next frontier in plant functional studies.

Systematic thermal analysis of the Arabidopsis proteome: Thermal tolerance, organization, and evolution.

Understanding the thermal stability of the plant proteome in the context of the native cellular environment would aid the design of crops with high thermal tolerance, but only limited such data are available. Here, we applied quantitative mass spectrometry to profile the thermal stability of the Arabidopsis proteome and identify thermo-sensitive and thermo-resilient protein networks in Arabidopsis, providing a basis for understanding heat-induced damage. We also show that the similarities of the protein-melting curves can be used as a proxy to evaluate system-wide protein-protein interactions in non-engineered plants and enable the identification of transient interactions exhibited by metabolons in the context of the cellular environment. Finally, we report a systematic comparison of the thermal stability of paralogs in Arabidopsis to aid the investigation and understanding of gene duplication and protein evolution. Taken together, our results could have broad implications for the fields of plant thermal tolerance, plant protein assemblies, and evolution.