- 1College of Liberal Arts Science, National University of Defense Technology, Changsha, China
- 2Logistics Center, National University of Defense Technology, Changsha, China
Reflectin proteins are natural copolymers consisting of repeated canonical domains. They are located in a biophotonic system called Bragg lamellae and manipulate the dynamic structural coloration of iridocytes. Their biological functions are intriguing, but the underlying mechanism is not fully understood. Reflectin A1, A2, B1, and C were found to present distinguished cyto-/nucleoplasmic localization preferences in the work. Comparable intracellular localization was reproduced by truncated reflectin variants, suggesting a conceivable evolutionary order among reflectin proteins. The size-dependent access of reflectin variants into the nucleus demonstrated a potential model of how reflectins get into Bragg lamellae. Moreover, RfA1 was found to extensively interact with the cytoskeleton, including its binding to actin and enrichment at the microtubule organizing center. This implied that the cytoskeleton system plays a fundamental role during the organization and transportation of reflectin proteins. The findings presented here provide evidence to get an in-depth insight into the evolutionary processes and working mechanisms of reflectins, as well as novel molecular tools to achieve tunable intracellular transportation.
Introduction
Cephalopods (squid, octopus, and cuttlefish) have evolved a remarkable ability to manipulate light and alter their appearance via pigmentary and structural elements. In the uppermost layers of the dermis, chromatophore organs proportionally regulate coloration by stretching pigmentary sacs (Hanlon, 2007; Hanlon et al., 2009) while the structural coloration is produced by two classes of cells called iridocytes and leucophores, distributed on the periphery of iridocytes, Bragg reflectors are periodically stacked lamellae that present iridescence by multilayer interference (Tao et al., 2010; DeMartini et al., 2013a). Leucophores containing granular vesicles are responsible for producing bright white color by unselectively reflecting all incident light (Williams et al., 2019).
Interestingly, both periodically stacked lamellae in iridocytes (Tao et al., 2010; DeMartini et al., 2013a) and intracellular vesicles in leucophores (DeMartini et al., 2013b; Mäthger et al., 2013; Hanlon et al., 2018) are unique and composed of reflectin proteins. Reflectin proteins are exclusively expressed in cephalopods, and their members possess different numbers of canonical reflectin motifs in different sequential locations (Tao et al., 2010; DeMartini et al., 2015; Levenson et al., 2019). Guan et al. (2017) reported that reflectin motifs may be traced to a 24-bp transposon-like DNA fragment from the symbiotic bioluminescent bacterium Vibrio fischeri. Afterward, millions of years of self-replication and translocation of transposons forms a prosperous reflectin family. For example, in Doryteuthis opalescens (one representative of the most recently evolved Loliginid species), distinct reflectins, reflectins A1, A2, B, and C, have been identified and well-characterized (Tao et al., 2010; DeMartini et al., 2015; Levenson et al., 2019).
Numerous in vitro assays suggest that the dynamic condensation, folding, and hierarchical assembly of reflectin proteins play fundamental roles during the formation and regulation of those structural coloration elements in squids (DeMartini et al., 2015; Levenson et al., 2016; Levenson et al., 2019). Recently, Morse and his team have preliminarily defined reflectin proteins as a new group of intrinsically disordered proteins (IDPs) and suggested that transient liquid-liquid phase separation (LLPS) may dominate reflectin assembly (Levenson et al., 2019). The dynamic reflectin assembly properties and the intricate reflectin-based biophotonic systems have already inspired the development of various next-generation tunable photonic (Qin et al., 2013; Phan et al., 2015; Dennis et al., 2017) and electronic platforms and devices (Ordinario et al., 2014; Yu et al., 2014; Phan et al., 2016).
Despite the impressive breakthroughs achieved by the abovementioned references, knowledge about how reflectin molecules are spatially organized to form proteinaceous layers or droplets is still limited. Moreover, the limited existing information is mostly obtained from in-tube assays. It is necessary and interesting to explore reflectin properties in cells.
Recently, two different groups made groundbreaking attempts and used reflectins to modify mammalian cell function. The formation of protein-based photonic architecture endows engineered human cells with tunable transparency-changing and light-scattering capabilities by expressing reflectins in the human embryonic kidney (HEK-293) cells (Chatterjee et al., 2020; Ogawa et al., 2020). This is a milestone in the discovery of novel molecular tools to modify mammalian cells with new features (Tang et al., 2020).
However, a comprehensive understanding of reflectin intracellular behaviors is a prerequisite before exploring their application potential in the synthetic biology field one step further. A comparative analysis is not available to investigate the properties of these reflectins exclusively expressed by cephalopods due to the lack of homogenesis with other well-known proteins.
To overcome this drawback, we engineered HEK-293T cells to express either intact RfA1, RfA2, RfB1, or RfC. These proteins are effective controls to each other, with different self-assembly characteristics (Levenson et al., 2016; Levenson et al., 2019) and membrane-bound properties (Song et al., 2020). After transfection, reflectins are found to phase out from the crowded intracellular milieu, similar to what Atrouli (Chatterjee et al., 2020) and Junko (Ogawa et al., 2020) described. Beyond that, more dramatic is the exclusive distribution of RfA1 particles in the cytoplasm, while spherical RfC droplets only exist in the nucleoplasm. As an intermediate state, RfA2 and RfB1 condensates distributes in both the cytoplasm and nucleoplasm. To the best of our knowledge, this is the first report of discrepant intracellular features of reflectins.
As a result of molecular evolution, the number and localization of canonical motifs distributed throughout reflectins may dominate their functional heterogeneities. Though it is impossible to witness this evolution process that takes a million years, the work took advantage of reflectin programmable sequences and designed a series of RfA1 truncations to validate this hypothesis. By cutting off reflectin motifs (RMs) one by one from RfA1, the constitution and properties of RfA1 variants are more and more similar to the shorter natural reflectins, indicating a possible ancestor-descendant sequence and a tight relationship between motifs repetition and reflectin function.
Further, CoIP-MS was conducted to better understand the molecular mechanism underlying the localization preferences among RfA1 and variants. Molecular partners such as RAN and Nups were identified to facilitate the entrance of RfA1 truncations into nuclei. Except for mass spectrum data, both confocal observation and RNA-Seq analysis confirmed the extensive interplay between RfA1 and the cytoskeleton. Especially, RfA1 was observed to be enriched in the microtubule organizing center (MTOC), significantly affecting spindle organization and cell division at the transcriptomic level. This spatial enrichment can be abolished by nocodazole treatment, as it suppresses microtubule polymerization. Therefore, it is speculated that the enrichment of RfA1 may be the result of its orientated movement towards the microtubule minus-ends. This conjectural MT minus-end-directed movement is likely to recruit RfA1 beneath the cytomembrane in polarized epithelial cells, which should be essential for Bragg lamellae formation in iridocytes.
Briefly, the work presents many unprecedented experimental results to better understand the relations and distinctions among reflectins, the penitential mechanism underlying reflectin orientated intracellular transportation, and the formation of biophotonic structures. Moreover, their intracellular localization properties provide alternative options for manipulating the cytoplasmic/nucleic distribution of proteins of interest.
Results
Selective Intracellular Localization of Reflectin Proteins
RfA1, A2, and B1 possess different numbers and localization of two types of canonical reflectin motifs (see Figure 1A), the regular reflectin motifs {RMs [M/FD(X)5MD(X)5MDX3/4]}, and a N-terminal reflectin motif {RMN [MEPMSRM-(T/S)MDF(H/Q)GR (Y/L) (I/M)DS (M/Q) (G/D)R (I/M)VDP (R/G)]} (see Figure 1B). Specifically, RfC is the shortest reflectin, containing a GMXX motif and RM*. The GMXX motif is a unique region of increased overall hydrophobicity composed of a four-amino acid repeat, and “X” represents less conserved locations within the repeat. Asterisk marked RM* of RfC contains substantial deviations in sequence not observed in any other reflectin motifs (Crookes et al., 2004) (see Figure 1B). The disorder tendencies of reflectins were calculated by PONDR (Xue et al., 2010), DISOclust (Ward et al., 2004; McGuffin, 2008), and ESpritz (Walsh et al., 2012). The possibility of forming protein droplets via liquid-liquid phase separation is also predicted by FuzDrop (Hardenberg et al., 2020). These calculations indicate a high disorder (see Supplementary Figure S1) and LLPS tendency (see Supplementary Figure S2) of four reflectin proteins. Reflectin proteins are introduced into the HEK-293T cells to verify these computational calculations.
FIGURE 1. Transfection of reflectins induces protein condensate formation in selective intracellular regions. (A) Cladogram of RfA1, RfA2, RfB1, and RfC. (B) Schematic reflectins sequences. Reflectin motifs are designated by boxes, while reflectin linkers (RLs) are lines. (C) Confocal images of GFP_Reflectin condensates in the HEK-293T cells. Nuclei and membrane are stained with DAPI (blue) and DiD (red) respectively, scale bars = 10 μM. (D) Cell number statistics and quantification of transfection efficiencies of four reflectins. (E) Distribution ratio statistics of fluorescent intensity in transfected cells and their nuclei. Data are presented as mean values ± SD, for n = 3 independent experiments. p-values (from a two-sided unpaired t test):∗, p < 0.05; ∗∗, p < 0.01; ∗∗∗, p < 0.001; all groups are compared to RfA1.
After transfection and culturing for 12 h or longer, all RfA1, RfA2, RfB1, and RfC are phased out from the cytoplasm or nucleoplasm, which is consistent with computational predictions. However, their difference is more significant. For RfA1, protein condensates are exclusively enriched in the cytoplasm. In a sharp contrast, the perfect spherical RfC droplets are only located in the nucleoplasm. While the distribution of RfA2 and RfB1 seems to stay in a middle state, spherical protein droplets mainly exist in the nucleus, with the minority staying in the cytoplasm (see Figures 1C,E).
According to the report by Guan et al., reflectin genes in cephalopods come from a transposon in symbiotic Vibrio fischeri (Guan et al., 2017). The repetitions of conserved motifs in reflectin sequences are the results of transposon self-replication and translocation. Thus, it is speculated that diverse reflectin proteins with different motifs stay at different evolution stages or represent different evolution directions.
From this point, RfC with only a GMXX motif and RM* may be the earliest appeared ancestor, while RfA1 with a RMN and five RMs is the further evolved offspring protein. Based on this consideration, truncated variants were constructed in the work by cutting off the RMs from the RfA1 one by one, to make their composition comparable to shorter natural reflectins. The dynamic intracellular performance of these RfA1 variants were then investigated.
Intracellular Localization Preferences of RfA1 Variations
Six pairs of primers were designed to clone RfA1 truncations (see Supplementary Table S1). The PCR products responsible for coding six peptides, including RMNto5, RMNto4, RMNto3, RMNto2, RMNto1, and RMN (schematics shown in Figure 2A), were subsequently ligated to vector pEGFP-C1. After transfection, the expressions of GFP-tagged RfA1 truncations were validated by an anti-GFP Western blotting (Figure 2B).
FIGURE 2. Expression of RfA1 variants in the HEK293T cells. (A) Schematics of RfA1 truncated variants. (B) Expression of RfA1 variants validated by Western blotting and indicated by anti-GFP antibodies. The molecular weight of EGFP is 27 kDa and GADPH is 36 kDa. (C) Cyto-/nucleo-localization preferences of truncated RfA1 derivatives. Nuclei and membrane are stained with DAPI (blue) and DiD (red), RfA1-derived peptides are indicated by tandem EGFP (green), scale bars = 10 μM. (D) Cell number statistics and quantification of transfection efficiencies of four reflectins. (E) Distribution ratio statistics of fluorescent intensity in transfected cells and their nuclei. Data are presented as mean values ± SD, for n = 3 independent experiments. p-values (from a two-sided unpaired t test): ∗, p < 0.05; ∗∗, p < 0.01; ∗∗∗, p < 0.001; All groups are compared to RMNto2.
Cells were harvested and RfA1-binding protein was extracted by an immunoprecipitation kit (Sangon) 24 h after transfection. GFP monoclonal antibody was used to validate the expression of GFP-tandem RfA1 truncations, while GADPH monoclonal antibody was employed to adjust the sample-loading quantity. In Figure 2B, the intensities of 36kD GADPH are almost the same among different lanes, implying that the sample loadings are normalized. The numbers labeled above bands are molecular weights calculated by ExPASy ProtParam software, consistent with the sums of GFP and RMNtoX (X = 0, 1, 2, 3, 4, and 5). Moreover, for longer variants RMNto5, RMNto4, and RMNto3, there are significant smears with very large molecular weights. Meanwhile, the lanes of shorter variants are quite clean. It indicates that RfA1 truncations with more motifs tend to tightly interact with other cellular proteins, which is discussed in detail in the next section.
On the other hand, a confocal image system is employed to investigate the intracellular localization of GFP-tagged RfA1 truncations (see Figure 2C). The nuclei and membrane are stained with DAPI (blue) and DiD (red), and RfA1-derived peptides are indicated by tandem EGFP (green). For cells transfected with pEGFP-C1-RMN and pEGFP-C1-RMNto1, though the peptide condensation seems absent, a remarkable enrichment of the green fluorescence signal in the nuclei areas appears. Amazingly, this nucleus-enrichment is executed more thoroughly by RMNto2, along with the recurrence of proteinaceous droplets. With increase in conserved motifs, RMNto3 are found to be distributed in both cyto- and nucleo-plasm. For longer derivatives with more motifs, condensates of RMNto4 and RMNto5 are blocked out from nuclei and located exclusively in cytoplasmic localization, extremely resembling their template RfA1.
Briefly, RMNto2 is most similar to RfC among all variants, which forms condensates exclusively in the nucleoplasm. More significantly, RMNto2 is also a definite boundary: truncations shorter than RMNto2 tend to enrich in the nuclei but cannot form droplets; derivatives longer than RMNto2 escape from the nuclei and accumulate and condensate into droplets in the cytoplasm.
Note smears in PAGE lanes of RMNto3, RMNto4, and RMNto5 in Figure 2B. Those are not sample contaminations or non-specific labeling. Subsequent CoIP-MS results reveal that various cellular components tightly interact with RfA1 variants, which assists their phase separation and diverse cellular localization.
Identification of RfA1/RMNto2 Binding Proteins
A monoclonal antibody of GFP was prepared and a CoIP-MS approach was used to identify protein-protein interactions, thus exploring the mechanism underlying phase separation and distinguished distribution of RfA1 and RMNto2.
Aforementioned, there must be a massive interaction between RfA1 variants and host cell proteins, resulting in smears in Figure 2B. These smears are not meaningless miscellaneous proteins, but hybrid aggregates recruited by RfA1 derivatives. Therefore, valuable information would be lost if only one single band was collected and analyzed. However, the collection of multiple bands as a mixed sample will inevitably bring in disturbing data. Based on these considerations, cells that express RMN are selected as a control group instead of GFP-expressed cells since RMN presents no significant phase separation or selective distribution behavior. As shown in Figure 3A, the IgG and IP bands of RMN, RMNto2, and RfA1 are collected and sent for mass-spectrum analysis. After that, proteins specifically interacting with RMN, RMNto2,and RfA1 are screened with a simple mathematical logic: Unique (RMN) = IP(RMN)-IP(RMN)∩IgG (RMN); Unique (RMNto2) = IP(RMNto2) -IP(RMNto2)∩IgG (RMNto2); Unique (RfA1) = IP(RfA1)-IP(RfA1)∩IgG (RfA1) (see Figures 3B,C).
FIGURE 3. Identification of interaction proteins targeted to RfA1 and RMNto2. (A) Silver staining of IP products, whole IP bands (marked by boxes) are sent to LC-MS. (B,C) PPI summary. Unique (RMN) = IP(RMN)-IP(RMN)∩IgG (RMN); Unique (RMNto2) = IP(RMNto2)-IP(RMNto2)∩IgG (RMNto2); Unique (RfA1) = IP(RfA1)-IP(RfA1)∩IgG (RfA1).
To further identify proteins which are responsible for intracellular localization or the formation of condensates, Unique (RMN), Unique (RMNto2), and Unique (RfA1) are divided into subsets according to the venn diagram in Figure 4A. According to the venn diagram, overlaid sets Set4/6/7 are eliminated as backgrounds. Within Set1 are the proteins combined into whole-length RfA1 and contribute to the phase separation and cytoplasmic localization. Besides, proteins in Set3 promote the importation of RMNto2 into nuclei and its phase out from the nucleoplasm. Co-targets of RMNto2 and RfA1 are sorted into Set 5.
FIGURE 4. GO biological processes annotation of proteins associated with LLPS and cellular localization. (A) Venn diagram showing the results of LC-MS, the graph is divided into Set1∼Set7 based on inclusion-exclusion principle of datasets. (B) Interpretations of different sets in the Venn graph. (C) Gene lists in Set1, Set5 and Set3 are converted to ENTREZ ID, and R package clusterProfiler is used to perform Gene Ontology annotation. LLPS and cellular localization-associated BP terms are showed in the figure, and the thresholds set as log10 (adjusted p-value) < −2.5.
Gene Ontology (GO) annotation is then performed to investigate biological processes associated with interacting proteins of RfA1 and RMNto2. Figure 4C shows the items of interest with an adjusted p-value<0.05 (see Supplementary Table S2 for gene lists and Supplementary Figure S3 for PPI networks).
Taking Set1-i for example, since RfA1 binds to proteins like HNRNPA1 (Matthew Michael et al., 1995; Roy et al., 2014), SRSF1 (Hautbergue et al., 2017), and RAN (Connor et al., 2003; Güttler and Görlich, 2011) and participates in processes including “protein export from nucleus”, RfA1 condensates are exclusively distributed in the cytoplasm. Oppositely, RANBP2 (Zhang et al., 2010) and NUPs (Cautain et al., 2015) in Set3-v facilitate the importation of RNNto2 into nuclei.
One remarkable common point among Set1-ii, Set5-iii, and Set3-iv is the involvement of “ribonucleoprotein complex” or “ribonucleotide metabolic”. The interaction with ribonucleoprotein complexes brings in a massive interplay with proteins and mRNA (Sfakianos et al., 2016; Mukherjee et al., 2019), leading to almost inevitable phase separation. This is consistent with the results in Figure 2C. However, the involvement of the ribonucleoprotein complex makes the situation more elusive. It is equivocal whether RfA1/RNNto2 molecules gather together and are phased out from the surroundings on their own, or they are drawn into the formation processes of ribonucleoprotein complex bodies.
Moreover, PPI information also demonstrates that RfA1 extensively interacts with actin and actin-binding proteins (ABPs) (Supplementary Table S3 for gene lists and Supplementary Figure S4 for PPI networks). These proteins are in charge of various biological processes, including actin polymerization, nucleate assembly of new filaments, and promotion of elongation (Uribe and Jay, 2009; Pollard, 2016). Considering the critical role of actin filaments in forming protrusive structures such as lamellipodia and ruffles (Mejillano et al., 2004; Chhabra and Higgs, 2007; Weiner et al., 2007; Innocenti, 2018), the results presented here highly suggest RfA1’s potential contribution in the establishment and regulation of Bragg lamellae via its collaboration with the actin system. RfA1 can anchor at the inside of the Bragg lamella by cooperating with actin and ABPs, or even participate in the organization of the actin network, which regulates cell morphogenesis.
Verification of Co-localization of RfA1 and Actin System
Actin should be indispensable for either the establishment or regulation of delicate Bragg reflectors in iridocytes. However, so far, except for the Co-IP results presented in this work, no evidence has been issued about the interaction between reflectin proteins and the cytoskeleton. As a preliminary start, we stained the HEK-293T cells with ActinRed 555 (Thermo) to trace the actin networks and verified the potential relationship between RfA1 and the cytoskeleton.
One notable trouble is the ubiquity of the actin cytoskeleton. It is challenging to judge the co-localization of RfA1 and actin (see Supplement Movie. 1mmc1). However, after viewing thousands of cells, an interesting scenario emerges. For plenty of transfected cells, there is a bright green sphere in the center or on the central axis, which are condensates of RfA1 (see Figure 5A for GFP channel, Supplementary Figure S5). These central localized protein condensates are not found in RfA2, RfB1, or RfC transfected cells. Around these spherical condensates, the red signal of stained actin is also significantly stronger than in other areas (see Figure 5A for ActinRed channel). The fluorescence statistics of selected areas in Figure 5B show an almost synchronized fluctuation of ActinRed and GFP signals, indicating their co-localization (see Supplementary Table S4 for raw data).
FIGURE 5. Appearance of spherical RfA1 condensates in the center or on the central axis of transfected cells. (A) RfA1 condensates located in the center or on the central axis of cells. Nuclei and actin are stained with DAPI (blue) and ActinRed 555 (red), RfA1 are indicated by tandem EGFP (green), scale bars = 10 μM. (B) Statistics of fluorescence signal distribution in selected areas. (C) 3D images of a transfected cell, with green RfA1 condensates located in the center, scale bars = 10 μM (see Supplementary Movie S2).
RfA1 Interacts With Microtubules and is Highly Enriched in the Microtubule Organization Center
Except for the fluorescent observation of the actin system, RNA-seq results also strongly suggested the tight relationship between RfA1 and the cytoskeleton. Cells transfected with pEGFP-C1-RfA1 and no-load pEGFP-C1 (set as control) are harvested for sequencing. After Metascape screening (Zhou et al., 2019), 461 differentially expressed unigenes are identified in RfA1 expressed cells (Supplementary Table S5). GO analysis is then performed to annotate differentially expressed genes (DEGs). Figure 6A shows the selected items of interest. Genes evolved in biological processes such as “regulation of microtubule cytoskeleton organization,” “microtubule polymerization or depolymerization,” and “microtubule nucleation” respond to RfA1 expression. Correspondingly, these microtubule-relevant genes are also sorted into cellular component items, including “spindle microtubule,” “spindle pole centrosome,” “mitotic spindle,” and “spindle” (Wittmann et al., 2001; Foley and Kapoor, 2013), which indicates the tight relationship between RfA1 and microtubules.
FIGURE 6. Effects of RfA1 on microtubule and mitosis. (A) GO annotation of DEGs between pEGFP-C1 transfected cells and pEGFP-C1-RfA1 transfected cells. Genes are converted to ENTREZ ID and R package clusterProfiler is used to perform GO annotation. Biological Process, Cellular Component and Molecular Function terms are shown in the figure, and the threshold was set as the p-value< 0.05. (B) RfA1 co-localization with microtubule networks in the center, in interphase cells. Nuclei and microtubule are stained with DAPI (blue) and Tubulin-Tracker Red (Beyotime), RfA1 are indicated by tandem EGFP (green), scale bars = 10 μM. (C) RfA1 located at the spindle center, in late mitosis cells. Nuclei and microtubule are stained with DAPI (blue) and Tubulin-Tracker Red (Beyotime), RfA1 is indicated by tandem EGFP (green), scale bars = 5 μM (see Supplementary Movies S3, S4).
The analysis based on RNA-Seq is further confirmed by confocal observation. After staining the fixed cells with Tubulin-Tracker Red (Beyotime), RfA1 condensates are significantly enriched at MTOC, located around the centrosome (see Figure 6B). For cells staying in late mitosis, RfA1 condensates are no longer located in the central areas but seem to chase the movements of daughter centrosomes (see Figure 6C and Supplementary Movies S3, S4 for Z-axis scanning). This significant enrichment of RfA1 at MOTC can be severely abolished by nocodazole treatment (data not shown). More directly, except for fluorescent colocalization observation, RfA1 is also detected to bind to tubulin subunits TBB2B, TBB4A, and TBB4B (Supplementary Table S6) by CoIP, suggesting physical interaction between RfA1 and microtubules.
Discussion
In summary, the work realized the expression of four natural reflectins (RfA1, RfA2, RfB1, and RfC) and six engineered RfA1 in cells and found their tunable intracellular localization. To the best of our knowledge, this is the first report of discrepant intracellular features of reflectins and their derivatives. Further, CoIP-MS was conducted to better understand the molecular mechanism underlying the localization preferences among RfA1 and variants. Moreover, both confocal observation and RNA-Seq analysis confirmed the extensive interplay between RfA1 and the cytoskeleton, providing evidence to delineate the role of RfA1 during Bragg lamellae establishment (Figure 7A). There are multiple reasons to support the significance of the work.
FIGURE 7. (A) Schematic of the workflow. CDS of squid reflectin proteins are obtained from NCBI. Codons are optimized according to mammal cell expression preference and introduced into HEK-293T cells by lipofection. RfA1 truncated variants are designed and cloned into pEGFP-C1 vectors. Anti-GFP monoclonal antibody is used to pull down GFP-tagged RfA1 variants and their binding proteins. Immunoprecipitates are separated using SDS-PAGE. Gel lanes are digested with trypsin for LC-MS/MS analysis. Colocalization of RfA1 and actin cytoskeleton is confirmed by confocal microscope. RNA-Seq is employed to globally investigate the influence of RfA1 expression on cells. (B). Illustration for hypothetically models. i. Origin and evolution of reflectin proteins. ii. Microtubule minus-end-directed transportation of reflectin. iii. Schematic for RfA1 derivatives as programmable and oriented cargo carriers.
Firstly, the cyto-/nucleoplasmic localization preference presented here could be evidence for the molecular evolution of reflectin proteins. According to Guan’s report, reflectin genes in cephalopods come from a transposon in symbiotic Vibrio fischeri (Guan et al., 2017). Afterward, millions of years of replication and translocation of this transposon forms reflectins with different lengths, motifs, and properties (Tao et al., 2010; DeMartini et al., 2015; Levenson et al., 2017; Levenson et al., 2019). In the work, RfA1 consists of one RMN and five RMs exclusively enriched in cytoplasm, while shorter RfA2, RfB1, and special RfC contains only one GMXX and an RM* tends to stay in nuclei. By truncating RfA1 gradually, shorter peptide derivatives (e.g., RMNto2) localize in nuclei, behaving exactly similar to simpler reflectins. From this point of view, the repetition numbers of RMs (as basic units) accurately determine the disparate properties (e.g., intracellular localization) of RfA1, A2, B1, and C. Being a subordinate element and of evolutionary origin, the 24-bp transposon-like DNA fragment is expected to be the root source of reflectin diversification, which makes our finding a strong clue to support Guan’s evolution hypothesis. Besides, earlier and simpler reflectins with fewer motifs tend to stay in the nuclei since there are abundant molecular tools available in the nuclei to evaluate, to manipulate, or even to edit this passage of alien DNA sequence. As for more complex reflectins such as RfA1, they are excluded from the nuclei and participate in other functional processes in the cytoplasm, like establishing and regulating Bragg lamellae [see Figure 7b(i)].
Secondly, the nucleus entry of short reflectins and centrosome-enrichment of RfA1 delineate their potential organization mode in iridocytes, which lay the material foundation for Bragg lamellae establishment. 1) Seeing from the longitudinal section, the pore at the bottom of Bragg lamellae is approximately equal to the size of the ciliary pore and nuclear pore (see Figure 7b(ii)). Though knowledge about protein transportation between the cytoplasm and Bragg lamellae is absent, the entrance of RfA1 truncations and short reflectins into nucleoplasm is observed. Here, RAN and Nups are identified to mediate the access of RfA1 variants into the nuclei. Presumably, the entry of reflectins into Bragg lamellae follows a similar mechanism, under the guidance of similar components (RAN, Nups, etc.). 2) Both CoIP-MS analysis and the RNA-Seq survey reveal the extensive interactions between RfA1 and the cytoskeleton. Actin and ABPs were identified to be combined with RfA1. Considering the critical role of actin filaments in the formation of protrusive structures such as lamellipodia and ruffles (Mejillano et al., 2004; Chhabra and Higgs, 2007; Weiner et al., 2007; Innocenti, 2018), the results presented here highly suggest RfA1’s potential contribution to the establishment and regulation of Bragg lamellae via its collaboration with the actin system. Moreover, RfA1 is observed to enrich at MTOC and co-localize with centrosomes. This spatial enrichment of RfA1 can be effectively abolished by nocodazole treatment. Since nocodazole suppresses the microtubule polymerization, this phenomenon suggests that the integrity of microtubule may be fundamental to RfA1 orientated transportation. Therefore, it is speculated that the orientated trafficking along the MT networks should be the critical property for RfA1 to work in Bragg lamellae. For the HEK-293T cells, microtubule minus-ends are located at or around the centrosome, guiding the movement of RfA1 to MTOC. However, microtubules are oriented apico-basally in polarized epithelial cells, with their minus-ends directed towards the apical side of the cells (Sugioka and Sawa, 2012; Glotfelty et al., 2014). This is partly verified when we perform new studies on the MDCK cells, which oppose apico-basal paralleled microtubule networks (data not shown). To the best of our knowledge, this is the first experiment to depict the potential mechanism underlying Bragg lamellae formation [see Figure 7B(ii)]. Although there is a huge difference between the mammalian cells and squid cells, the evolutionary conservatism of actins among algae, amoeba, fungi, and animals (Dominguez and Holmes, 2011; Gunning et al., 2015) eliminates this drawback.
Thirdly, reflectins are natural block copolymers and programmable biomaterials. Their intracellular localization properties provide alternative options for the manipulation of the cytoplasmic/nucleic distribution of proteins of interest. If we take GFP as a molecular cargo, reflectin derivatives can be regarded as intelligent vehicles to transport cargoes to pre-selected destinations [see Figure 7B(iii)]. Compared to the commonly-used NLS sequence, our programmable tools also have better scalability and modularity. The distribution of selected proteins can be easily altered by linking programmable tools to proteins (merely changing the numbers of parts), making them present only in the cytoplasm or nuclei. Using NLS, however, leads to a one-way trip into the nucleus. As a quick example, our tools might hypothetically achieve a switch of “cytoplasmic/nucleic distribution” at the protein-protein interaction level by linking different lengths of our parts to drug-/light-induced protein dimerization pairs. Further studies have already started in our research group (unpublished yet).
Methods
Construction of Recombinant pEGFP-C1 Vectors
The nucleotide sequence of D. (Loligo) pealeii reflectin A1 (RfA1) (Genbank: ACZ57764.1) and D. (Loligo) Opalescens reflectin C (Genbank: AIN36559.1) were optimized for human cell expression, and then synthesized and sequencing-identified by Sangon Biotech® (Shanghai, China). Primers (F-GAATTCTAT GAATAGATATTTGAATAGACA; R-GGATCCATACATATGATAATCATAATA ATTT) were designed to introduce EcoR I and BamH I cutting sites, so the modified RfA1 CDS can be constructed into pEGFP-C1 via a standard restriction-enzyme cloning process. As for the truncated RfA1 derivatives, six pairs of primers were coupled. For example, if RMN-F and RM3-R primers were selected, then a nucleotide sequence responsible for the coding of RMN-RL1-RM1-RL2-RM2-RL3-RM3 (simplified as RMNto3 in the work) would be obtained after PCR. Meanwhile, 5′ GCATGGACGAGCTGTACAAG 3′ and 5′ TTATGATCAG- TTATCTAGAT 3’ were added to F-primers and R-primers during the synthesis of the primer, respectively, which enables the sequences to be ligated to pEGFP-C1 by Ready-to-Use Seamless Cloning Kit from Sangon Biotech® (Shanghai, China).
Growth and Transfection of Human Cells
The HEK-293T cells (ATCC®, CRL-3216TM) were cultured on plastic dishes in Dulbecco’s Modified Eagle Medium (DMEM, GibcoTM) supplemented with 10% fetal bovine serum (FBS, GibcoTM) at 37°C under 5% CO2. One day before transfection, cells were seeded at ∼33% of the confluent density for the glass bottom dishes from Cellvis (California, United States) and grown for another 24 h. Then a transfection mixture containing Lipofectamine 3,000 (Thermo Scientific) and recombinant vectors were added to the medium and incubated for ∼16–∼24 h. For CCK-8 test, 1 × 104 cells were seeded into each hole of 96-well plates 1 day before transfection, then transfected with recombinant vectors and incubated for another 24 h. After that, 10 μl CCK-8 solution was added to the wells for a chromogenic reaction of ∼2–∼4 h. OD450 was detected by Multiskan FC (Thermo Scientific).
Fluorescence Microscopy of Stained Cells
Transfected HEK-293T cells grown in Cellvis plastic dishes were first fixed with 4% paraformaldehyde at room temperature for 30 min, then stained with DiD or ActinRed (diluted in 0.5% Triton X-100 PBS) for ∼30 min after PBS rinses. After washing off the fluorescent dye with PBS, fixed cells were embedded in DAPI-Fluoromount (Beyotime, Shanghai, China) and characterized with a Leica TCS SP8 imaging system in fluorescence imaging mode. The resulting images were analyzed with ImageJ (Java 1.8.0_172/1.52b) (Schindelin et al., 2012).
Online Computational Prediction Tools
Disorder prediction of RfA1 and derivatives was conducted with PONDER® (http://www.pondr.com/), DISOclust (http://www.reading.ac.uk/bioinf-/DISOclust/) and ESpritz (http://protein.bio.unipd.it/espritz/). The probability of liquid-liquid phase separation was calculated by FuzPred (http://protdyn-fuzpred.org/).
SDS-PAGE, Western Blotting
Total proteins of cells were extracted by using the immunoprecipitation kit (Sangon). Protein samples were separated by SDS-PAGE and transferred to polyvinylidene fluoride (PVDF) filter membranes (Millipore, United States) for immune-hybridization. After 1 h of blocking in PBST (phosphate-buffered saline containing 0.05% Tween-20 and 2.5% BSA), the membranes were incubated with one of the following primary antibodies at corresponding concentration: anti-GFP mouse monoclonal antibody (Sangon), and anti-GAPDH mouse monoclonal antibody (Sangon). Subsequently, band visualization was performed by electro-chemiluminescence (ECL) and detected by the Digit imaging system (Thermo, Japan).
In-Gel Digestion
For enzymolysis in protein gels (the gel strip, protein solution, and protein freeze-dried powder), a clean blade was used to dig out the strip of interest and cut it into 0.5–0.7 mm cubes. After decolorizing with test stain/silver stain decolorizing solution, the rubber block was washed with 500 μl acetonitrile solution three times until the rubber particles become completely white. Add 500 μl of 10 mM DTT to a water bath at 56°C for 30 min. Use low-speed centrifugation for a reduction reaction. Then add 500 μl decolorizing solution and mix at room temperature for 5–10 min Wash the gel spots once, and centrifuge at low speed to discard the supernatant After quickly adding 500 μl 55 mM IAM, place it in a dark room at room temperature for 30 min. Centrifuge at a low speed for an alkylation reaction. Then add 500 μl of decolorizing solution and mix at room temperature for 5–10 min After washing the gel spots once, centrifuge at low speed, and discard the supernatant; then add 500 μl of acetonitrile until the colloidal particles become completely white, and vacuum dry for 5 min. Add 0.01 μg/μl trypsin according to the gel volume. After an ice bath for 30 min, add an appropriate amount of 25 mM NH4HCO3 into pH8.0 enzymatic hydrolysis buffer, and enzymatically hydrolyze overnight at 37°C. After enzymolysis, add 300 μl extract with sonicate for 10 min. Centrifuge at a low speed, to collect the supernatant. After repeating the process twice, combine the obtained extracts and perform a vacuum dry.
Zip-Tip Desalting
Dissolve the sample with 10–20 μl 0.2% TFA, and centrifuge at 10,000 rpm for 20 min. Then wet the Ziptip 15 times with a wetting solution and equilibrate with a balance solution 10 times. Finally, inhale the sample solution for 10 cycles and blow 8 times with the rinse liquid. Note: To avoid cross-contamination, aliquot 100 μl flushing solution per tube and use a separate tube for each sample. Add 50 μl eluent to a clean EP tube, and pipette it repeatedly to elute the peptide. Finally, drain the sample.
LC-MS/MS Analysis
The peptide samples were diluted to 1 μg/μl on the machine buffet. Set the sample volume to 5 μl and collect the scan mode for 60 min. Scan the peptides with a mass-to-charge ratio of 350–1200 in the sample. The mass spectrometry data was collected using the Triple TOF 5600 + LC/MS system (AB SCIEX, United States). The peptide samples were dissolved in 2% acetonitrile/0.1%formic acid and analyzed using the Triple TOF 5600 plus mass spectrometer coupled with the Eksigent nano LC system (AB SCIEX, United States). The peptide solution was added to the C18 capture column (3 μM, 350 μM × 0.5 mm, AB Sciex, United States) and the C18 analytical column (3 μM, 75 μM × 150) was applied with a 60 min time gradient at a flow rate of 300 nl/min mm(Welch Materials, Inc.) for gradient elution. The two mobile phases were buffer A (2% acetonitrile/0.1%formic acid/98% H2O) and buffer B (98% acetonitrile/0.1% formic acid/2% H2O). For IDA (Information Dependent Acquisition), the MS spectrum was scanned with an anion accumulation time of 250 ms, and the MS spectrum of 30 precursor ions was acquired with an anion accumulation time of 50 ms. Collect MS1 spectrum in the range of 350–1200 m/z, and collect MS2 spectrum in the range of 100–1500 m/z. Set the precursor ion dynamic exclusion time to 15 s.
Genes identified from LC-MS and RNA-Seq were converted to ENTREZ ID and an R package, clusterProfiler, was used to perform gene ontology annotation (Yu et al., 2012). The mass spectrometry proteomics data have been deposited to the Proteome Xchange Consortium via the PRIDE (Perez-Riverol et al., 2022) partner repository with the dataset identifier PXD031350. RNA-Seq data were uploaded to the NCBI GEO database, numbered as GSE186861.
All the data needed to evaluate the conclusions are present in the paper and/or the supplementary information. All other relevant data are available from the authors upon reasonable request.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: Proteome Xchange Consortium via the PRIDE partner repository with the dataset identifier PXD031350; RNA sequencing data NCBI GEO accession: GSE186861.
Author Contributions
JS designed and performed the experiments, analyzed data, and completed the manuscript. CL designed and performed the experiments, analyzed the LC-MS data. LZ assisted in directing graduate students to carry out experiments. BL and LL performed the experiments. ZY analyzed the data. TM designed and plotted the schematic diagrams in Figure 1 and Figure 8. WW and BH designed the experiments and wrote the manuscript.
Funding
We are grateful to the National Natural Science Foundation of China (Grants 31971291) and the Natural Science Foundation of Hunan Province (Grants S2020JJQNJJ1684) for their financial support.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We acknowledge Prof. X. Zhu (SIBCB, CAS) and Prof. L. Qiu (MBL, HNU) for valuable discussion. JS gives his personal thanks to Prof. D. E. Morse and R. Levenson (MCDB, UCSB), for their devoted guidance from 2016 to 2018.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.862011/full#supplementary-material
Reference
Cautain, B., Hill, R., de Pedro, N., and Link, W. (2015). Components and Regulation of Nuclear Transport Processes. FEBS J. 282, 445–462. doi:10.1111/febs.13163
Chatterjee, A., Cerna Sanchez, J. A., Yamauchi, T., Taupin, V., Couvrette, J., and Gorodetsky, A. A. (2020). Cephalopod-inspired Optical Engineering of Human Cells. Nat. Commun. 11, 2708. doi:10.1038/s41467-020-16151-6
Chhabra, E. S., and Higgs, H. N. (2007). The Many Faces of Actin: Matching Assembly Factors with Cellular Structures. Nat. Cell Biol. 9, 1110–1121. doi:10.1038/ncb1007-1110
Connor, M. K., Kotchetkov, R., Cariou, S., Resch, A., Lupetti, R., Beniston, R. G., et al. (2003). CRM1/Ran-Mediated Nuclear Export of p27Kip1Involves a Nuclear Export Signal and Links P27 Export and Proteolysis. Mol. Biol. Cell 14, 201–213. doi:10.1091/mbc.e02-06-0319
Crookes, W. J., Ding, L.-L., Huang, Q. L., Kimbell, J. R., Horwitz, J., and McFall-Ngai, M. J. (2004). Reflectins: the Unusual Proteins of Squid Reflective Tissues. Science 303, 235–238. doi:10.1126/science.1091288
DeMartini, D. G., Krogstad, D. V., and Morse, D. E. (2013). Membrane Invaginations Facilitate Reversible Water Flux Driving Tunable Iridescence in a Dynamic Biophotonic System. Proc. Natl. Acad. Sci. U.S.A. 110, 2552–2556. doi:10.1073/pnas.1217260110
DeMartini, D. G., Ghoshal, A., Pandolfi, E., Weaver, A. T., Baum, M., and Morse, D. E. (2013). Dynamic Biophotonics: Female Squid Exhibit Sexually Dimorphic Tunable Leucophores and Iridocytes. J. Exp. Biol. 216, 3733–3741. doi:10.1242/jeb.090415
DeMartini, D. G., Izumi, M., Weaver, A. T., Pandolfi, E., and Morse, D. E. (2015). Structures, Organization, and Function of Reflectin Proteins in Dynamically Tunable Reflective Cells. J. Biol. Chem. 290, 15238–15249. doi:10.1074/jbc.m115.638254
Dennis, P. B., Singh, K. M., Vasudev, M. C., Naik, R. R., and Crookes-Goodson, W. J. (2017). Research Update: a Minimal Region of Squid Reflectin for Vapor-Induced Light Scattering. Apl. Mater. 5, 120701. doi:10.1063/1.4997199
Dominguez, R., and Holmes, K. C. (2011). Actin Structure and Function. Annu. Rev. Biophys. 40, 169–186. doi:10.1146/annurev-biophys-042910-155359
Foley, E. A., and Kapoor, T. M. (2013). Microtubule Attachment and Spindle Assembly Checkpoint Signalling at the Kinetochore. Nat. Rev. Mol. Cell Biol. 14, 25–37. doi:10.1038/nrm3494
Glotfelty, L. G., Zahs, A., Iancu, C., Shen, L., and Hecht, G. A. (2014). Microtubules Are Required for Efficient Epithelial Tight Junction Homeostasis and Restoration. Am. J. Physiol.-Cell Physiol. 307, C245–C254. doi:10.1152/ajpcell.00336.2013
Guan, Z., Cai, T., Liu, Z., Dou, Y., Hu, X., Zhang, P., et al. (2017). Origin of the Reflectin Gene and Hierarchical Assembly of its Protein. Curr. Biol. 27, 2833–2842. doi:10.1016/j.cub.2017.07.061
Gunning, P. W., Ghoshdastider, U., Whitaker, S., Popp, D., and Robinson, R. C. (2015). The Evolution of Compositionally and Functionally Distinct Actin Filaments. J. Cell Sci. 128, 2009–2019. doi:10.1242/jcs.165563
Güttler, T., and Görlich, D. (2011). Ran-Dependent Nuclear Export Mediators: a Structural Perspective. EMBO J. 30, 3457–3474. doi:10.1038/emboj.2011.287
Hanlon, R. T., Chiao, C.-C., Mäthger, L. M., Barbosa, A., Buresch, K. C., and Chubb, C. (2009). Cephalopod Dynamic Camouflage: Bridging the Continuum between Background Matching and Disruptive Coloration. Phil. Trans. R. Soc. B 364, 429–437. doi:10.1098/rstb.2008.0270
Hanlon, R. T., Mäthger, L. M., Bell, G. R. R., Kuzirian, A. M., and Senft, S. L. (2018). White Reflection from Cuttlefish Skin Leucophores. Bioinspir. Biomim. 13, 035002. doi:10.1088/1748-3190/aaa3a9
Hanlon, R. (2007). Cephalopod Dynamic Camouflage. Curr. Biol. 17, R400–R404. doi:10.1016/j.cub.2007.03.034
Hardenberg, M., Horvath, A., Ambrus, V., Fuxreiter, M., and Vendruscolo, M. (2020). Widespread Occurrence of the Droplet State of Proteins in the Human Proteome. Proc. Natl. Acad. Sci. U.S.A. 117, 33254–33262. doi:10.1073/pnas.2007670117
Hautbergue, G. M., Castelli, L. M., Ferraiuolo, L., Sanchez-Martinez, A., Cooper-Knock, J., Higginbottom, A., et al. (2017). SRSF1-Dependent Nuclear Export Inhibition of C9ORF72 Repeat Transcripts Prevents Neurodegeneration and Associated Motor Deficits. Nat. Commun. 8, 16063. doi:10.1038/ncomms16063
Innocenti, M. (2018). New Insights into the Formation and the Function of Lamellipodia and Ruffles in Mesenchymal Cell Migration. Cell Adh Migr. 12, 401–416. doi:10.1080/19336918.2018.1448352
Levenson, R., Bracken, C., Bush, N., and Morse, D. E. (2016). Cyclable Condensation and Hierarchical Assembly of Metastable Reflectin Proteins, the Drivers of Tunable Biophotonics. J. Biol. Chem. 291, 4058–4068. doi:10.1074/jbc.m115.686014
Levenson, R., DeMartini, D. G., and Morse, D. E. (2017). Molecular Mechanism of Reflectin's Tunable Biophotonic Control: Opportunities and Limitations for New Optoelectronics. Apl. Mater. 5, 104801. doi:10.1063/1.4985758
Levenson, R., Bracken, C., Sharma, C., Santos, J., Arata, C., Malady, B., et al. (2019). Calibration between Trigger and Color: Neutralization of a Genetically Encoded Coulombic Switch and Dynamic Arrest Precisely Tune Reflectin Assembly. J. Biol. Chem. 294, 16804–16815. doi:10.1074/jbc.ra119.010339
Mäthger, L. M., Senft, S. L., Gao, M., Karaveli, S., Bell, G. R. R., Zia, R., et al. (2013). Bright White Scattering from Protein Spheres in Color Changing, Flexible Cuttlefish Skin. Adv. Funct. Mat. 23, 3980–3989. doi:10.1002/adfm.201203705
Matthew Michael, W., Choi, M., and Dreyfuss, G. (1995). A Nuclear Export Signal in hnRNP A1: a Signal-Mediated, Temperature-dependent Nuclear Protein Export Pathway. Cell 83, 415–422. doi:10.1016/0092-8674(95)90119-1
McGuffin, L. J. (2008). Intrinsic Disorder Prediction from the Analysis of Multiple Protein Fold Recognition Models. Bioinformatics 24, 1798–1804. doi:10.1093/bioinformatics/btn326
Mejillano, M. R., Kojima, S.-i., Applewhite, D. A., Gertler, F. B., Svitkina, T. M., and Borisy, G. G. (2004). Lamellipodial Versus Filopodial Mode of the Actin Nanomachinery. Cell 118, 363–373. doi:10.1016/j.cell.2004.07.019
Mukherjee, N., Wessels, H.-H., Lebedeva, S., Sajek, M., Ghanbari, M., Garzia, A., et al. (2019). Deciphering Human Ribonucleoprotein Regulatory Networks. Nucleic Acids Res. 47, 570–581. doi:10.1093/nar/gky1185
Ogawa, J., Iwata, Y., Tonnu, N. U., Gopinath, C., Huang, L., Itoh, S., et al. (2020). Genetic Manipulation of the Optical Refractive Index in Living Cells. NY,US: bioRxiv.
Ordinario, D. D., Phan, L., Walkup IV, W. G., Jocson, J.-M., Karshalev, E., Hüsken, N., et al. (2014). Bulk Protonic Conductivity in a Cephalopod Structural Protein. Nat. Chem. 6, 596–602. doi:10.1038/nchem.1960
Perez-Riverol, Y., Bai, J., Bandla, C., García-Seisdedos, D., Hewapathirana, S., Kamatchinathan, S., et al. (2022). The PRIDE Database Resources in 2022: A Hub for Mass Spectrometry-Based Proteomics Evidences. Nucleic Acids Res. 50 (D1), D543–D552. doi:10.1093/nar/gkab1038
Phan, L., Ordinario, D. D., Karshalev, E., Walkup IV, W. G., Shenk, M. A., and Gorodetsky, A. A. (2015). Infrared Invisibility Stickers Inspired by Cephalopods. J. Mat. Chem. C 3, 6493–6498. doi:10.1039/c5tc00125k
Phan, L., Kautz, R., Leung, E. M., Naughton, K. L., Van Dyke, Y., and Gorodetsky, A. A. (2016). Dynamic Materials Inspired by Cephalopods. Chem. Mat. 28, 6804–6816. doi:10.1021/acs.chemmater.6b01532
Pollard, T. D. (2016). Actin and Actin-Binding Proteins. Cold Spring Harb. Perspect. Biol. 8, 1–17. doi:10.1101/cshperspect.a018226
Qin, G., Dennis, P. B., Zhang, Y., Hu, X., Bressner, J. E., Sun, Z., et al. (2013). Recombinant Reflectin-Based Optical Materials. J. Polym. Sci. B Polym. Phys. 51, 254–264. doi:10.1002/polb.23204
Roy, R., Durie, D., Li, H., Liu, B.-Q., Skehel, J. M., Mauri, F., et al. (2014). hnRNPA1 Couples Nuclear Export and Translation of Specific mRNAs Downstream of FGF-2/S6K2 Signalling. Nucleic Acids Res. 42, 12483–12497. doi:10.1093/nar/gku953
Schindelin, J., Arganda-Carreras, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, T., et al. (2012). Fiji: An Open-Source Platform for Biological-Image Analysis. Nat. Methods 9, 676–682. doi:10.1038/nmeth.2019
Sfakianos, A. P., Whitmarsh, A. J., and Ashe, M. P. (2016). Ribonucleoprotein Bodies Are Phased in. Biochem. Soc. Trans. 44, 1411–1416. doi:10.1042/bst20160117
Song, J., Levenson, R., Santos, J., Velazquez, L., Zhang, F., Fygenson, D., et al. (2020). Reflectin Proteins Bind and Reorganize Synthetic Phospholipid Vesicles. Langmuir 36, 2673–2682. doi:10.1021/acs.langmuir.9b03632
Sugioka, K., and Sawa, H. (2012). Formation and Functions of Asymmetric Microtubule Organization in Polarized Cells. Curr. Opin. Cell Biol. 24, 517–525. doi:10.1016/j.ceb.2012.05.007
Tang, T.-C., An, B., Huang, Y., Vasikaran, S., Wang, Y., Jiang, X., et al. (2020). Materials Design by Synthetic Biology. Nat. Rev. Mater 6, 332–350. doi:10.1038/s41578-020-00265-w
Tao, A. R., DeMartini, D. G., Izumi, M., Sweeney, A. M., Holt, A. L., and Morse, D. E. (2010). The Role of Protein Assembly in Dynamically Tunable Bio-Optical Tissues. Biomaterials 31, 793–801. doi:10.1016/j.biomaterials.2009.10.038
Uribe, R., and Jay, D. (2009). A Review of Actin Binding Proteins: New Perspectives. Mol. Biol. Rep. 36, 121–125. doi:10.1007/s11033-007-9159-2
Walsh, I., Martin, A. J. M., Di Domenico, T., and Tosatto, S. C. E. (2012). ESpritz: Accurate and Fast Prediction of Protein Disorder. Bioinformatics 28, 503–509. doi:10.1093/bioinformatics/btr682
Ward, J. J., Sodhi, J. S., McGuffin, L. J., Buxton, B. F., and Jones, D. T. (2004). Prediction and Functional Analysis of Native Disorder in Proteins from the Three Kingdoms of Life. J. Mol. Biol. 337, 635–645. doi:10.1016/j.jmb.2004.02.002
Weiner, O. D., Marganski, W. A., Wu, L. F., Altschuler, S. J., and Kirschner, M. W. (2007). An Actin-Based Wave Generator Organizes Cell Motility. PLoS Biol. 5, e221. doi:10.1371/journal.pbio.0050221
Williams, T. L., Senft, S. L., Yeo, J., Martín-Martínez, F. J., Kuzirian, A. M., Martin, C. A., et al. (2019). Dynamic Pigmentary and Structural Coloration within Cephalopod Chromatophore Organs. Nat. Commun. 10, 1004. doi:10.1038/s41467-019-08891-x
Wittmann, T., Hyman, A., and Desai, A. (2001). The Spindle: a Dynamic Assembly of Microtubules and Motors. Nat. Cell Biol. 3, E28–E34. doi:10.1038/35050669
Xue, B., Dunbrack, R. L., Williams, R. W., Dunker, A. K., and Uversky, V. N. (2010). PONDR-FIT: A Meta-Predictor of Intrinsically Disordered Amino Acids. Biochim. Biophys. Acta (BBA) - Proteins Proteomics 1804, 996–1010. doi:10.1016/j.bbapap.2010.01.011
Yu, G., Wang, L.-G., Han, Y., and He, Q.-Y. (2012). clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. OMICS A J. Integr. Biol. 16, 284–287. doi:10.1089/omi.2011.0118
Yu, C., Li, Y., Zhang, X., Huang, X., Malyarchuk, V., Wang, S., et al. (2014). Adaptive Optoelectronic Camouflage Systems with Designs Inspired by Cephalopod Skins. Proc. Natl. Acad. Sci. U.S.A. 111, 12998–13003. doi:10.1073/pnas.1410494111
Zhang, R., Mehla, R., and Chauhan, A. (2010). Perturbation of Host Nuclear Membrane Component RanBP2 Impairs the Nuclear Import of Human Immunodeficiency Virus -1 Preintegration Complex (DNA). Plos One 5, e15620. doi:10.1371/journal.pone.0015620
Keywords: reflectins, truncated peptides, intracellular localization, molecular evolution, cytoskeleton interaction
Citation: Song J, Liu C, Li B, Liu L, Zeng L, Ye Z, Mao T, Wu W and Hu B (2022) Tunable Cellular Localization and Extensive Cytoskeleton-Interplay of Reflectins. Front. Cell Dev. Biol. 10:862011. doi: 10.3389/fcell.2022.862011
Received: 25 January 2022; Accepted: 09 May 2022;
Published: 23 June 2022.
Edited by:
Venkaiah Betapudi, United States Department of Homeland Security, United StatesReviewed by:
Marco Saltini, Uppsala University, SwedenYaroslav B. Blume, National Academy of Sciences of Ukraine (NAN Ukraine), Ukraine
Copyright © 2022 Song, Liu, Li, Liu, Zeng, Ye, Mao, Wu and Hu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Junyi Song, c29uZ2p1bnlpQG51ZHQuZWR1LmNu; Biru Hu, aHViaXJ1MDhAbnVkdC5lZHUuY24=
†These authors have contributed equally to this work