RESUMEN

The nearly neutral theory of molecular evolution posits variation among species in the effectiveness of selection. In an idealized model, the census population size determines both this minimum magnitude of the selection coefficient required for deleterious variants to be reliably purged, and the amount of neutral diversity. Empirically, an 'effective population size' is often estimated from the amount of putatively neutral genetic diversity and is assumed to also capture a species' effectiveness of selection. A potentially more direct measure of the effectiveness of selection is the degree to which selection maintains preferred codons. However, past metrics that compare codon bias across species are confounded by among-species variation in %GC content and/or amino acid composition. Here, we propose a new Codon Adaptation Index of Species (CAIS), based on Kullback-Leibler divergence, that corrects for both confounders. We demonstrate the use of CAIS correlations, as well as the Effective Number of Codons, to show that the protein domains of more highly adapted vertebrate species evolve higher intrinsic structural disorder.

Evolution is the process through which populations change over time, starting with mutations in the genetic sequence of an organism. Many of these mutations harm the survival and reproduction of an organism, but only by a very small amount. Some species, especially those with large populations, can purge these slightly harmful mutations more effectively than other species. This fact has been used by the 'drift barrier theory' to explain various profound differences amongst species, including differences in biological complexity. In this theory, the effectiveness of eliminating slightly harmful mutations is specified by an 'effective' population size, which depends on factors beyond just the number of individuals in the population. Effective population size is normally calculated from the amount of time a 'neutral' mutation (one with no effect at all) stays in the population before becoming lost or taking over. Estimating this time requires both representative data for genetic diversity and knowledge of the mutation rate. A major limitation is that these data are unavailable for most species. A second limitation is that a brief, temporary reduction in the number of individuals has an oversized impact on the metric, relative to its impact on the number of slighly harmful mutations accumulated. Weibel, Wheeler et al. developed a new metric to more directly determine how effectively a species purges slightly harmful mutations. Their approach is based on the fact that the genetic code has 'synonymous' sequences. These sequences code for the same amino acid building block, with one of these sequences being only slightly preferred over others. The metric by Weibel, Wheeler et al. quantifies the proportion of the genome from which less preferred synonymous sequences have been effectively purged. It judges a population to have a higher effective population size when the usage of synonymous sequences departs further from the usage predicted from mutational processes. The researchers expected that natural selection would favour 'ordered' proteins with robust three-dimensional structures, i.e., that species with a higher effective population size would tend to have more ordered versions of a protein. Instead, they found the opposite: species with a higher effective population size tend to have more disordered versions of the same protein. This changes our view of how natural selection acts on proteins. Why species are so different remains a fundamental question in biology. Weibel, Wheeler et al. provide a useful tool for future applications of drift barrier theory to a broad range of ways that species differ.

Asunto(s)

RESUMEN

The m6A epitranscriptomic mark is the most abundant and widespread internal RNA chemical modification, which through the control of RNA acts as an important factor of eukaryote reproduction, growth, morphogenesis and stress response. The main m6A readers constitute a super family of proteins with hundreds of members that share a so-called YTH RNA binding domain. The majority of YTH proteins carry no obvious additional domain except for an Intrinsically Disordered Region (IDR). In Arabidopsis thaliana IDRs are important for the functional specialization among the different YTH proteins, known as Evolutionarily Conserved C-Terminal region, ECT 1 to 12. Here by studying the ECT2 protein and using an in vitro biochemical characterization, we show that full-length ECT2 and its YTH domain alone have a distinct ability to bind m6A, conversely to previously characterized YTH readers. We identify peptide regions outside of ECT2 YTH domain, in the N-terminal IDR, that regulate its binding to m6A-methylated RNA. Furthermore, we show that the selectivity of ECT2 binding for m6A is enhanced by a high uridine content within its neighbouring sequence, where ECT2 N-terminal IDR is believed to contact the target RNA in vivo. Finally, we also identify small structural elements, located next to ECT2 YTH domain and conserved in a large set of YTH proteins, that enhance its binding to m6A-methylated RNA. We propose from these findings that some of these regulatory regions are not limited to ECT2 or YTH readers of flowering plants but may be widespread among eukaryotic YTH readers.

Asunto(s)

RESUMEN

In eukaryotes, the origin recognition complex (ORC) faciliates the assembly of pre-replicative complex (pre-RC) at origin DNA for replication licensing. Here we show that the N-terminal intrinsically disordered region (IDR) of the yeast Orc2 subunit is crucial for this process. Removing a segment (residues 176-200) from Orc2-IDR or mutating a key isoleucine (194) significantly inhibits replication initiation across the genome. These Orc2-IDR mutants are capable of assembling the ORC-Cdc6-Cdt1-Mcm2-7 intermediate, which exhibits impaired ATP hydrolysis and fails to be convered into the subsequent Mcm2-7-ORC complex and pre-RC. These defects can be partially rescued by the Orc2-IDR peptide. Moreover, the phosphorylation of this Orc2-IDR region by S cyclin-dependent kinase blocks its binding to Mcm2-7 complex, causing a defective pre-RC assembly. Our findings provide important insights into the multifaceted roles of ORC in supporting origin licensing during the G1 phase and its regulation to restrict origin firing within the S phase.

Asunto(s)

RESUMEN

Intrinsically disordered proteins (IDPs) perform a broad range of functions in biology, suggesting that the ability to design IDPs could help expand the repertoire of proteins with novel functions. Computational design of IDPs with specific conformational properties has, however, been difficult because of their substantial dynamics and structural complexity. We describe a general algorithm for designing IDPs with specific structural properties. We demonstrate the power of the algorithm by generating variants of naturally occurring IDPs that differ in compaction, long-range contacts, and propensity to phase separate. We experimentally tested and validated our designs and analyzed the sequence features that determine conformations. We show how our results are captured by a machine learning model, enabling us to speed up the algorithm. Our work expands the toolbox for computational protein design and will facilitate the design of proteins whose functions exploit the many properties afforded by protein disorder.

Asunto(s)

RESUMEN

Tuftelin Interacting Protein 11 (TFIP11) was identified as a critical human spliceosome assembly regulator, interacting with multiple proteins and localising in membrane-less organelles. However, a lack of structural information on TFIP11 limits the rationalisation of its biological role. TFIP11 is predicted as an intrinsically disordered protein (IDP), and more specifically concerning its N-terminal (N-TER) region. IDPs lack a defined tertiary structure, existing as a dynamic conformational ensemble, favouring protein-protein and protein-RNA interactions. IDPs are involved in liquid-liquid phase separation (LLPS), driving the formation of subnuclear compartments. Combining disorder prediction, molecular dynamics, and spectroscopy methods, this contribution shows the first evidence TFIP11 N-TER is a polyampholytic IDP, exhibiting a structural duality with the coexistence of ordered and disordered assemblies, depending on the ionic strength. Increasing the salt concentration enhances the protein conformational flexibility, presenting a more globule-like shape, and a fuzzier unstructured arrangement that could favour LLPS and protein-RNA interaction. The most charged and hydrophilic regions are the most impacted, including the G-Patch domain essential to TFIP11 function. This study gives a better understanding of the salt-dependent conformational behaviour of the N-TER TFIP11, supporting the hypothesis of the formation of different types of protein assembly, in line with its multiple biological roles.

Asunto(s)

RESUMEN

Complementary studies by Zhao et al.1 and Chen et al.2 reveal how an intrinsically disordered region in MED13 controls mutually exclusive binding of RNA Polymerase II and CDK8 kinase module to Mediator, switching Mediator and transcription activation on and off.

Asunto(s)

RESUMEN

In this work, we explored the intrinsic disorder status of the three members of the synuclein family of proteins-α-, ß-, and γ-synucleins-and showed that although all three human synucleins are highly disordered, the highest levels of disorder are observed in γ-synuclein. Our analysis of the peculiarities of the amino acid sequences and modeled 3D structures of the human synuclein family members revealed that the pathological mutations A30P, E46K, H50Q, A53T, and A53E associated with the early onset of Parkinson's disease caused some increase in the local disorder propensity of human α-synuclein. A comparative sequence-based analysis of the synuclein proteins from various evolutionary distant species and evaluation of their levels of intrinsic disorder using a set of commonly used bioinformatics tools revealed that, irrespective of their origin, all members of the synuclein family analyzed in this study were predicted to be highly disordered proteins, indicating that their intrinsically disordered nature represents an evolutionary conserved and therefore functionally important feature. A detailed functional disorder analysis of the proteins in the interactomes of the human synuclein family members utilizing a set of commonly used disorder analysis tools showed that the human α-synuclein interactome has relatively higher levels of intrinsic disorder as compared with the interactomes of human ß- and γ- synucleins and revealed that, relative to the ß- and γ-synuclein interactomes, α-synuclein interactors are involved in a much broader spectrum of highly diversified functional pathways. Although proteins interacting with three human synucleins were characterized by highly diversified functionalities, this analysis also revealed that the interactors of three human synucleins were involved in three common functional pathways, such as the synaptic vesicle cycle, serotonergic synapse, and retrograde endocannabinoid signaling. Taken together, these observations highlight the critical importance of the intrinsic disorder of human synucleins and their interactors in various neuronal processes.

Asunto(s)

RESUMEN

Temperature can impact every reaction essential to a cell. For organisms that cannot regulate their own temperature, adapting to temperatures that fluctuate unpredictably and on variable timescales is a major challenge. Extremes in the magnitude and frequency of temperature changes are increasing across the planet, raising questions as to how the biosphere will respond. To examine mechanisms of adaptation to temperature, we collected wild isolates from different climates of the fungus Ashbya gossypii, which has a compact genome of only ∼4,600 genes. We found control of the nuclear division cycle and polarized morphogenesis, both critical processes for fungal growth, were temperature sensitive and varied among the isolates. The phenotypes were associated with naturally varying sequences within the glutamine-rich region (QRR) IDR of an RNA-binding protein called Whi3. This protein regulates both nuclear division and polarized growth via its ability to form biomolecular condensates. In cells and in cell-free reconstitution assays, we found that temperature tunes the properties of Whi3-based condensates. Exchanging Whi3 sequences between isolates was sufficient to rescue temperature-sensitive phenotypes, and specifically, a heptad repeat sequence within the QRR confers temperature-sensitive behavior. Together, these data demonstrate that sequence variation in the size and composition of an IDR can promote cell adaptation to growth at specific temperature ranges. These data demonstrate the power of IDRs as tuning knobs for rapid adaptation to environmental fluctuations.

Asunto(s)

RESUMEN

An intrinsically disordered protein (IDP) or region (IDR) lacks or has little protein structure but still maintains function. This lack of structure creates flexibility and fluidity, allowing multiple protein conformations and potentially transient interactions with more than one partner. Caliciviruses are positive-sense ssRNA viruses, containing a relatively small genome of 7.6-8.6 kb and have a broad host range. Many viral proteins are known to contain IDRs, which benefit smaller viral genomes by expanding the functional proteome through the multifunctional nature of the IDR. The percentage of intrinsically disordered residues within the total proteome for each calicivirus type species can range between 8 and 23%, and IDRs have been experimentally identified in NS1-2, VPg and RdRP proteins. The IDRs within a protein are not well conserved across the genera, and whether this correlates to different activities or increased tolerance to mutations, driving virus adaptation to new selection pressures, is unknown. The function of norovirus NS1-2 has not yet been fully elucidated but includes involvement in host cell tropism, the promotion of viral spread and the suppression of host interferon-λ responses. These functions and the presence of host cell-like linear motifs that interact with host cell caspases and VAPA/B are all found or affected by the disordered region of norovirus NS1-2. The IDRs of calicivirus VPg are involved in viral transcription and translation, RNA binding, nucleotidylylation and cell cycle arrest, and the N-terminal IDR within the human norovirus RdRP could potentially drive liquid-liquid phase separation. This review identifies and summarises the IDRs of proteins within the Caliciviridae family and their importance during viral replication and subsequent host interactions.

Asunto(s)

RESUMEN

Several mammalian genes have originated from the domestication of retrotransposons, selfish mobile elements related to retroviruses. Some of the proteins encoded by these genes have maintained virus-like features; including self-processing, capsid structure formation, and the generation of different isoforms through -1 programmed ribosomal frameshifting. Using quantitative approaches in molecular evolution and biophysical analyses, we studied 28 retrotransposon-derived genes, with a focus on the evolution of virus-like features. By analyzing the rate of synonymous substitutions, we show that the -1 programmed ribosomal frameshifting mechanism in three of these genes (PEG10, PNMA3, and PNMA5) is conserved across mammals and originates alternative proteins. These genes were targets of positive selection in primates, and one of the positively selected sites affects a B-cell epitope on the spike domain of the PNMA5 capsid, a finding reminiscent of observations in infectious viruses. More generally, we found that retrotransposon-derived proteins vary in their intrinsically disordered region content and this is directly associated with their evolutionary rates. Most positively selected sites in these proteins are located in intrinsically disordered regions and some of them impact protein posttranslational modifications, such as autocleavage and phosphorylation. Detailed analyses of the biophysical properties of intrinsically disordered regions showed that positive selection preferentially targeted regions with lower conformational entropy. Furthermore, positive selection introduces variation in binary sequence patterns across orthologues, as well as in chain compaction. Our results shed light on the evolutionary trajectories of a unique class of mammalian genes and suggest a novel approach to study how intrinsically disordered region biophysical characteristics are affected by evolution.

Asunto(s)

RESUMEN

All cyanobacteria and some chemoautotrophic bacteria fix CO2 into sugars using specialized proteinaceous compartments called carboxysomes. Carboxysomes enclose the enzymes Rubisco and carbonic anhydrase inside a layer of shell proteins to increase the CO2 concentration for efficient carbon fixation by Rubisco. In the ⍺-carboxysome lineage, a disordered and highly repetitive protein named CsoS2 is essential for carboxysome formation and function. Without it, the bacteria require high CO2 to grow. How does a protein predicted to be lacking structure serve as the architectural scaffold for such a vital cellular compartment? In this study, we identify key residues present in the repeats of CsoS2, VTG and Y, which are necessary for building functional ⍺-carboxysomes in vivo. These highly conserved and repetitive residues contribute to the multivalent binding interaction and phase separation behavior between CsoS2 and shell proteins. We also demonstrate 3-component reconstitution of CsoS2, Rubisco, and shell proteins into spherical condensates and show the utility of reconstitution as a biochemical tool to study carboxysome biogenesis. The precise self-assembly of thousands of proteins is crucial for carboxysome formation, and understanding this process could enable their use in alternative biological hosts or industrial processes as effective tools to fix carbon.

Asunto(s)

RESUMEN

Human ZC3H11A is an RNA-binding zinc finger protein involved in mRNA export and required for the efficient growth of human nuclear replicating viruses. Its biochemical properties are largely unknown so our goal has been to produce the protein in a pure and stable form suitable for its characterization. This has been challenging since the protein is large (810 amino acids) and with only the N-terminal zinc finger domain (amino acids 1-86) being well structured, the remainder is intrinsically disordered. Our production strategies have encompassed recombinant expression of full-length, truncated and mutated ZC3H11A variants with varying purification tags and fusion proteins in several expression systems, with or without co-expression of chaperones and putative interaction partners. A range of purification schemes have been explored. Initially, only truncated ZC3H11A encompassing the zinc finger domain could successfully be produced in a stable form. It required recombinant expression in insect cells since expression in E. coli gave a protein that aggregated. To reduce problematic nucleic acid contaminations, Cys8, located in one of the zinc fingers, was substituted by Ala and Ser. Interestingly, this did not affect nucleic acid binding, but the full-length protein was stabilised while the truncated version was insoluble. Ultimately, we discovered that when using alkaline buffers (pH 9) for purification, full-length ZC3H11A expressed in Sf9 insect cells was obtained in a stable and >90 % pure form, and as a mixture of monomers, dimers, tetramers and hexamers. Many of the challenges experienced are consistent with its predicted structure and unusual charge distribution.

Asunto(s)

RESUMEN

The abundant nuclear protein hnRNP U interacts with a broad array of RNAs along with DNA and protein to regulate nuclear chromatin architecture. The RNA-binding activity is achieved via a disordered ∼100 residue C-terminal RNA-binding domain (RBD) containing two distinct RGG/RG motifs. Although the RNA-binding capabilities of RGG/RG motifs have been widely reported, less is known about hnRNP U's RNA-binding selectivity. Furthermore, while it is well established that hnRNP U binds numerous nuclear RNAs, it remains unknown whether it selectively recognizes sequence or structural motifs in target RNAs. To address this question, we performed equilibrium binding assays using fluorescence anisotropy (FA) and electrophoretic mobility shift assays (EMSAs) to quantitatively assess the ability of human hnRNP U RBD to interact with segments of cellular RNAs identified from eCLIP data. These RNAs often, but not exclusively, contain poly-uridine or 5'-AGGGAG sequence motifs. Detailed binding analysis of several target RNAs reveal that the hnRNP U RBD binds RNA in a promiscuous manner with high affinity for a broad range of structured RNAs, but with little preference for any distinct sequence motif. In contrast, the isolated RGG/RG of hnRNP U motif exhibits a strong preference for G-quadruplexes, similar to that observed for other RGG motif bearing peptides. These data reveal that the hnRNP U RBD attenuates the RNA binding selectivity of its core RGG motifs to achieve an extensive RNA interactome. We propose that a critical role of RGG/RG motifs in RNA biology is to alter binding affinity or selectivity of adjacent RNA-binding domains.

Asunto(s)

RESUMEN

Whole genome and exome sequencing are reporting on hundreds of thousands of missense mutations. Taking a pan-disease approach, we explored how mutations in intrinsically disordered regions (IDRs) break or generate protein interactions mediated by short linear motifs. We created a peptide-phage display library tiling ~57,000 peptides from the IDRs of the human proteome overlapping 12,301 single nucleotide variants associated with diverse phenotypes including cancer, metabolic diseases and neurological diseases. By screening 80 human proteins, we identified 366 mutation-modulated interactions, with half of the mutations diminishing binding, and half enhancing binding or creating novel interaction interfaces. The effects of the mutations were confirmed by affinity measurements. In cellular assays, the effects of motif-disruptive mutations were validated, including loss of a nuclear localisation signal in the cell division control protein CDC45 by a mutation associated with Meier-Gorlin syndrome. The study provides insights into how disease-associated mutations may perturb and rewire the motif-based interactome.

Asunto(s)

RESUMEN

This study has explored the involvement of Intrinsically Disordered Proteins (IDPs) in cyanobacterial stress response. IDPs possess distinct physicochemical properties, which allow them to execute diverse functions. Anabaena PCC 7120, the model photosynthetic, nitrogen-fixing cyanobacterium encodes 688 proteins (11 % of the total proteome) with at least one intrinsically disordered region (IDR). Of these, 130 proteins that showed >30 % overall disorder were designated as IDPs. Physico-chemical analysis, showed these IDPs to adopt shapes ranging from 'globular' to 'tadpole-like'. Upon exposure to NaCl, 41 IDP-encoding genes were found to be differentially expressed. Surprisingly, most of these were induced, indicating the importance of IDP-accumulation in overcoming salt stress. Subsequently, six IDPs were identified to be induced by multiple stresses (salt, ammonium and selenite). Interestingly, the presence of these 6-multiple stress-induced IDPs was conserved in filamentous cyanobacteria. Utilizing the experimental proteomic data of Anabaena, these 6 IDPs were found to interact with many proteins involved in diverse pathways, underscoring their physiological importance as protein hubs. This study lays the framework for IDP-related research in Anabaena by (a) identifying, as well as physiochemically characterizing, all the disordered proteins and (b) uncovering a subset of IDPs that are likely to be critical in adaptation to environmental stresses.

Asunto(s)

RESUMEN

Germ granules are biomolecular condensates present in most animal germ cells. One function of germ granules is to help maintain germ cell totipotency by organizing mRNA regulatory machinery, including small RNA-based gene regulatory pathways. The C. elegans germ granule is compartmentalized into multiple subcompartments whose biological functions are largely unknown. Here, we identify an uncharted subcompartment of the C. elegans germ granule, which we term the E granule. The E granule is nonrandomly positioned within the germ granule. We identify five proteins that localize to the E granule, including the RNA-dependent RNA polymerase (RdRP) EGO-1, the Dicer-related helicase DRH-3, the Tudor domain-containing protein EKL-1, and two intrinsically disordered proteins, EGC-1 and ELLI-1. Localization of EGO-1 to the E granule enables synthesis of a specialized class of 22G RNAs, which derive exclusively from 5' regions of a subset of germline-expressed mRNAs. Defects in E granule assembly elicit disordered production of endogenous siRNAs, which disturbs fertility and the RNAi response. Our results define a distinct subcompartment of the C. elegans germ granule and suggest that one function of germ granule compartmentalization is to facilitate the localized production of specialized classes of small regulatory RNAs.

Asunto(s)

RESUMEN

The early phases of clathrin mediated endocytosis are organized through a highly complex interaction network mediated by clathrin associated sorting proteins (CLASPs) that comprise long intrinsically disordered regions (IDRs). AP180 is a CLASP exclusively expressed in neurons and comprises a long IDR of around 600 residues, whose function remains partially elusive. Using NMR spectroscopy, we discovered an extended and strong interaction site within AP180 with the major adaptor protein AP2, and describe its binding dynamics at atomic resolution. We find that the 70 residue-long site determines the overall interaction between AP180 and AP2 in a dynamic equilibrium between its bound and unbound states, while weaker binding sites contribute to the overall affinity at much higher concentrations of AP2. Our data suggest that this particular interaction site might play a central role in recruitment of adaptors to the clathrin coated pit, whereas more transient and promiscuous interactions allow reshaping of the interaction network until cargo uptake inside a coated vesicle.

Asunto(s)

RESUMEN

The essential Mediator (MED) coactivator complex plays a well-understood role in regulation of basal transcription in all eukaryotes, but the mechanism underlying its role in activator-dependent transcription remains unknown. We investigated modulation of metazoan MED interaction with RNA polymerase II (RNA Pol II) by antagonistic effects of the MED26 subunit and the CDK8 kinase module (CKM). Biochemical analysis of CKM-MED showed that the CKM blocks binding of the RNA Pol II carboxy-terminal domain (CTD), preventing RNA Pol II interaction. This restriction is eliminated by nuclear receptor (NR) binding to CKM-MED, which enables CTD binding in a MED26-dependent manner. Cryoelectron microscopy (cryo-EM) and crosslinking-mass spectrometry (XL-MS) revealed that the structural basis for modulation of CTD interaction with MED relates to a large intrinsically disordered region (IDR) in CKM subunit MED13 that blocks MED26 and CTD interaction with MED but is repositioned upon NR binding. Hence, NRs can control transcription initiation by priming CKM-MED for MED26-dependent RNA Pol II interaction.

Asunto(s)

RESUMEN

G protein-coupled receptors (GPCRs) control intracellular signaling cascades via agonist-dependent coupling to intracellular transducers including heterotrimeric G proteins, GPCR kinases (GRKs), and arrestins. In addition to their critical interactions with the transmembrane core of active GPCRs, all three classes of transducers have also been reported to interact with receptor C-terminal domains (CTDs). An underexplored aspect of GPCR CTDs is their possible role as lipid sensors given their proximity to the membrane. CTD-membrane interactions have the potential to control the accessibility of key regulatory CTD residues to downstream effectors and transducers. Here, we report that the CTDs of two closely related family C GPCRs, metabotropic glutamate receptor 2 (mGluR2) and mGluR3, bind to membranes and that this interaction can regulate receptor function. We first characterize CTD structure with NMR spectroscopy, revealing lipid composition-dependent modes of membrane binding. Using molecular dynamics simulations and structure-guided mutagenesis, we then identify key conserved residues and cancer-associated mutations that modulate CTD-membrane binding. Finally, we provide evidence that mGluR3 transducer coupling is controlled by CTD-membrane interactions in live cells, which may be subject to regulation by CTD phosphorylation and changes in membrane composition. This work reveals an additional mechanism of GPCR modulation, suggesting that CTD-membrane binding may be a general regulatory mode throughout the broad GPCR superfamily.

Asunto(s)

RESUMEN

Transcription factors (TFs) control specificity and activity of gene transcription, but whether a relationship between these two features exists is unclear. Here we provide evidence for an evolutionary trade-off between the activity and specificity in human TFs encoded as submaximal dispersion of aromatic residues in their intrinsically disordered protein regions. We identified approximately 500 human TFs that encode short periodic blocks of aromatic residues in their intrinsically disordered regions, resembling imperfect prion-like sequences. Mutation of periodic aromatic residues reduced transcriptional activity, whereas increasing the aromatic dispersion of multiple human TFs enhanced transcriptional activity and reprogramming efficiency, promoted liquid-liquid phase separation in vitro and more promiscuous DNA binding in cells. Together with recent work on enhancer elements, these results suggest an important evolutionary role of suboptimal features in transcriptional control. We propose that rational engineering of amino acid features that alter phase separation may be a strategy to optimize TF-dependent processes, including cellular reprogramming.

Asunto(s)