A novel automated image analysis pipeline for quantifying morphological changes to the endoplasmic reticulum in cultured human cells
BMC Bioinformatics volume 22, Article number: 427 (2021)
In mammalian cells the endoplasmic reticulum (ER) comprises a highly complex reticular morphology that is spread throughout the cytoplasm. This organelle is of particular interest to biologists, as its dysfunction is associated with numerous diseases, which often manifest themselves as changes to the structure and organisation of the reticular network. Due to its complex morphology, image analysis methods to quantitatively describe this organelle, and importantly any changes to it, are lacking.
In this work we detail a methodological approach that utilises automated high-content screening microscopy to capture images of cells fluorescently-labelled for various ER markers, followed by their quantitative analysis. We propose that two key metrics, namely the area of dense ER and the area of polygonal regions in between the reticular elements, together provide a basis for measuring the quantities of rough and smooth ER, respectively. We demonstrate that a number of different pharmacological perturbations to the ER can be quantitatively measured and compared in our automated image analysis pipeline. Furthermore, we show that this method can be implemented in both commercial and open-access image analysis software with comparable results.
We propose that this method has the potential to be applied in the context of large-scale genetic and chemical perturbations to assess the organisation of the ER in adherent cell cultures.
One significant challenge in the field of image analysis is the acquisition of images in a manner that permits readily quantifiable information about individual cells in a well or tissue. In this respect, automated analysis offers many advantages over manual assessment of images, as it is objective, reproducible, and significantly less labour intensive, enabling the analysis of large data sets in a fraction of the time . However, the design of logical sequential steps of image analysis, so-called automated analysis pipelines, using commercial and open-access software (e.g. Harmony®, CellProfiler ), can be anything but straightforward depending on the phenotype being studied. In most image analysis pipelines, the nucleus is the primary organelle that is identified, as it is a large, spherical or ovoid, single-copy organelle per cell, from which basic geometrical features such as length (long-axis), width (short-axis), area, perimeter and circularity can readily be extracted. The next step is identification of the cell plasma membrane, which in turn provides information on the location of the cytoplasm. After these steps, most image analysis strategies then focus on identifying the organelles or subcellular structures of interest in the cytoplasm and implementing appropriate analysis and quantification routines . Cellular organelles display a high degree of morphological diversity, and it is at this stage that image analysis pipelines diverge in order that the most appropriate methods to identify them can be applied. For example, although the Golgi complex is generally also a single organelle in mammalian cells, its appearance can be highly heterogeneous. Nevertheless, increasingly sophisticated image analysis routines are now able to report a variety of phenotypes seen in this organelle in response to various perturbations [4, 5]. Similarly, in recent years excellent progress has been made in automated detection and analysis of other organelles displaying punctate and tubular morphologies, such as endosomes . Indeed even profiling of the Golgi and mitochondria in living cells allowing analysis of organelle fission and fusion events has also been shown to be possible [7, 8]. Despite this progress, the highly dynamic and morphologically complex endoplasmic reticulum (ER) network possesses a significant challenge to analyse using automated image analysis methods. As a result, there is a lack of available tools for the quantitative analysis of ER morphology, despite an increasing body of literature linking ER morphology with human disease [9, 10].
The ER is an extensive network of sheets and tubules extending from the nuclear envelope to the plasma membrane. Despite maintaining continuity of the luminal space, the ER is highly compartmentalised, comprising the nuclear envelope (NE), ribosome-studded sheets of “rough” ER (RER), a network of smooth ER (SER) tubules [11, 12], as well as specialised domains termed ER exit sites, which play a key role in the first steps of secretory transport . Each compartment performs essential yet distinct cellular functions, such as protein synthesis and folding at the RER, and calcium storage, lipid metabolism, and mitochondrial division at the SER . ER compartments also differ structurally. The RER comprises two parallel lipid bilayers that extend over many microns with little curvature, which are sealed at curved rims leaving a 50 nm thick lumen. These sheets are usually stacked, connected by regions of twisted membranes and helical edges . By contrast, the SER comprises tubules with diameters ranging from 25 to 90 nm which are stabilised by ER-shaping proteins that create positive curvature [16,17,18].
Within post-mitotic cells, the ER is arranged such that the RER localises proximal to the nucleus and the Golgi apparatus, while the SER tubular network occupies both the perinuclear region and extends towards the cell periphery [19, 20]. However, ER organisation is dynamic and can alter depending on the physiological conditions of the cell. For example, the ER undergoes a complete reorganisation during mitosis by the interconversion of sheets-to-tubules and disassembly of the NE [21, 22]. Altered availability of ER resident proteins, cytoskeletal proteins, or phospholipids can also shift the balance between ER tubules and ER sheets [19, 23,24,25]. Maintaining ER organisation is vital to preserve the essential functions carried out by the ER . Altered ER morphology has been reported in the context of several pathological conditions, for example (1) during viral infection viruses co-opt host cell membranes for their replication and in many cases remodel the ER membrane through interactions with ER proteins creating so-called replication factories (RFs), which promote enrichment of cellular and viral components necessary for viral replication [26, 27]; (2) the presence of protein aggregates, the hallmark of many neurodegenerative diseases, is linked to alteration of ER organisation by increased presence of enlarged peripheral ER sheets [28, 29]; (3) defects in misfolded protein clearance by impaired ER-associated degradation, or accumulation of misfolded proteins in cystic fibrosis cells, both result in dilated ER elements [30, 31]; (4) the cytotoxic effects of several chemicals have been linked to changes on ER structure [32,33,34,35,36,37]; and (5) mutations in ER-shaping proteins, known to cause neurodegenerative disorders such as hereditary spastic paraplegias (HSPs), have been shown to disrupt ER organisation in both in vitro and in vivo models [24, 38,39,40].
Despite the clear importance of ER organisation in cellular function and dysfunction, image analysis tools for the reliable and reproducible quantification of ER conformations are extremely limited. Published reports largely rely on qualitative descriptions of changes in ER structures, while quantification of changes in ER organisation are scarce or absent in most model systems [26, 29, 41,42,43]. Recently a set of tools has been developed which enables the quantification of the ER network in plant cells . Given that the ER is a single copy organelle with vital roles in cell health and disease, there is a pressing need to establish objective tools for the description of this intricate organelle in cultured mammalian cells that will permit researchers to detect changes in its organisation under different conditions. In this study, we describe the design of a high-content platform for automated analysis of ER organisation. We have used this platform to study the effect of genetic and chemical modifiers on the organisation of ER in adherent cell cultures, demonstrating its robustness and applicability for wider use in various cell-based assays.
A high-content and high-resolution imaging strategy to visualise the endoplasmic reticulum
Successful implementation of downstream image analysis relies upon the quality and consistency of the input images, and so the first step was to optimise our imaging protocol in such a way that it would allow for subsequent scale-up. Experiments were therefore designed in 96-well optically clear plates, compatible with a fully automated confocal high-content screening microscope fitted with a 63 ×/1.15 NA water-immersion objective (Fig. 1A). This experimental set-up was designed to ensure that all imaging conditions remained unaltered throughout the experiments, and was compatible with potential upscaling to large-scale genetic perturbation or chemical screening. One major bottleneck in high-content imaging combined with high-resolution microscopy is data accumulation . Aiming to alleviate this problem, we sampled a central area of 1.53 mm2 in each well (Fig. 1B), ensuring that we captured a minimum of 50–100 cells per well, which is a number generally considered to be sufficient for accurate high-content analysis from images acquired with a 63 × objective lens (Fig. 1C). We next considered what would be the most relevant cell model with which to develop an image analysis pipeline for the ER. We initially selected U-2 OS cells, as these cells possess a large and flat morphology, with an extensive RER and SER network spread throughout their cytoplasm (Fig. 1D) . Specifically, we used a U-2 OS cell line stably expressing the ER marker Sec61β-mEmerald . This marker localises to both SER tubules and perinuclear RER, offers a strong signal to noise ratio (Fig. 1D), and its exogenous expression does not affect ER dynamics, as the microtubule binding site of the Sec61β protein is blocked by its fusion to the fluorescent protein mEmerald . These U-2 OS cells were plated into wells of a 96-well plate and imaged the following day using a 63 ×/1.15 NA water-immersion objective, achieving an effective xy resolution of 0.28 μm, which enabled visualisation of peripheral SER tubules and dense ER structures (brighter) in the perinuclear region (Fig. 1E). Although U-2 OS cells are strongly adherent, and therefore their cytoplasm and organelles extend laterally over several microns in almost only 2-dimensions, they still have a small degree of depth that could potentially lead to under-detection of continuous ER tubules when only analysing a single optical slice. In order to overcome this issue, we determined that best practice was to image several z planes to span the depth of the periphery of the cells, and subsequently to use the maximum projection of these slices as the input image for the analysis (Fig. 1F). This approach also assists in eliminating out-of-focus artefacts.
The analysis workflow developed is summarised in Fig. 2A and illustrates the principles followed. The first steps are standard in automated cell analysis and involve cell segmentation (Fig. 2B), and the filtering of the cell population based on nuclear properties to remove dividing and apoptotic cells (Fig. 2C). Given that we were working with a stably expressing cell line, we also included a filter step to remove cells that had either lost the exogenous DNA plasmid or were only weakly expressing it (Fig. 2D). We initially performed the analysis using Harmony® image analysis software, the bespoke analysis software for the automated microscope used to acquire the images, as this software allows easy stitching of multiple fields. However, in order to demonstrate that our analysis approach is not restricted to this software, we applied the same principles and recreated a similar analysis pipeline in CellProfiler  using maximum projection images of individual fields (Additional file 1: Fig. S1). Once all individual cells were segmented and filtered, we next built image analysis steps to quantify the ER organisation, specifically the area of SER polygons and the proportion of dense perinuclear ER (detailed below).
Definition of SER polygon regions
SER tubules form a network in a polygonal array, in which the “gaps” are termed polygonal regions . SER tubules interconnect at three-way junctions, creating the nodes of the SER network. This morphology was clearly visible in our U-2 OS cells expressing Sec61β-mEmerald (Fig. 3A). A previous study reported an increase in the number of three-way junctions per membrane surface upon over-expression of the three-way junction resident protein Lunapark, giving rise to a densely branched network with smaller polygons . Hypothesising that other perturbations to ER function may alter the number of three-way junctions and the density of the SER network, we explored the use of the polygonal region size as a readout of the density of the SER network in the cell periphery. The tubular SER network of polygons extend from the perinuclear RER, through the periphery of the cell to the plasma membrane [12, 48]. In order to quantify the area of these ER polygons, we first defined the cell periphery as a new region using the Harmony image analysis ‘building block’ termed “select cell region”, which performs radial calculations to shrink the cytoplasm, pushing the inner border (nuclear border) towards the cell perimeter depending on the distance between the two, resulting in a new a ring-shaped area which defines our region of interest (ROI) (Fig. 3B, blue, Additional file 2: Appendix 1, Table 1). Other software packages include alternative building blocks to redefine areas, for example in CellProfiler we used “object processing” to “shrink” the object “cell” and identified a “tertiary object” by subtracting this modified area from the object “cytoplasm”, identifying a new object called “ROI” (Additional file 1: Fig. S1C, Additional file 2: Appendix 2). Within the defined ROI, we next used the Sec61β-mEmerald signal to detect ER tubules, and thereby form an ER network ‘mask’ (Fig. 3C, green; Additional file 1: Fig. S1C). By inverting the ER mask, we obtained the regions bounded by the ER tubules and cell perimeter (Fig. 3D, yellow and blue respectively; Additional file 1: Fig. S1C). As the regions between the ER and cell perimeter (Fig. 3D, blue) do not represent true ER polygons, these were excluded by filtering regions that overlapped with the cell perimeter. In this way, we were able to define the regions of the cytoplasm bounded by ER tubules as polygons (Fig. 3D, yellow; Additional file 1: Fig. S1C) that could be analysed by the software, and measurements such as mean area per cell (Fig. 3E) generated as a readout. This pipeline enables us to screen for chemical or genetic modifiers of the tubular SER network.
Definition of perinuclear dense RER structure
The ribosome-studded RER is composed of a series of interconnected stacked sheets typically located in the perinuclear cell area [12, 48]. Although confocal microscopy cannot resolve ER substructures, we were able to clearly observe higher Sec61β-mEmerald intensity in the perinuclear cell area, relative to that seen in the periphery of the cell (Fig. 4A). The differential distribution in intensity of the ER marker becomes clearer on the pseudo-coloured image (Fig. 4B) with high intensity pixels coloured yellow or orange, and low intensity pixels coloured pink, purple or blue. The perinuclear region of high intensity Sec61β fluorescence likely corresponds to the area in which the RER resides, while the peripheral region corresponds to tubular SER, as previously discussed. Applying an intensity threshold, we generated a mask which broadly corresponds to the dense RER (Fig. 4C, yellow; Additional file 1: Fig. S1D). Once a mask of this region was created, we calculated the proportion of pixels of the cell area that overlap with it (Fig. 4D; Additional file 1: Fig. S1D). This analysis allowed us to calculate the proportion of RER within the cells (Fig. 4E), and determine how this sub compartment of the ER is altered upon treatment with chemical or genetic modifiers.
Drug-induced ER reorganisation in U-2 OS cells
Several studies have reported alteration of ER organisation by different chemical compounds [32, 37, 49]. A limitation of all these studies however is that they employed qualitative analysis and did not enable an objective description of the changes observed. By contrast, the analysis pipeline that we have developed permits the unbiased quantification of changes to the ER network organisation in cultured cells upon drug treatment. We first applied our image analysis pipeline to the U-2 OS cells expressing Sec61β-mEmerald, and quantified the steady-state area of the SER polygons and the dense RER (Fig. 5). These values provided us with benchmark values for all subsequent experiments. Interestingly, when we repeated these experiments in another cell line, HeLa Kyoto, transiently expressing the same Sec61β-mEmerald protein, the values we obtained for SER and RER areas were highly comparable (Additional file 1: Fig. S2). This demonstrates that our image analysis strategy is not only restricted to making measurements from U-2 OS cells. To validate that our pipeline can also detect alterations in ER organisation, we treated U-2 OS cells expressing Sec61β-mEmerald with various drugs that have been reported to disrupt the ER. Following treatment with the saturated fatty acid palmitic acid (250 μM for 16 h) or the mycotoxin cytochalasin B (10 μM for 30 min), U-2 OS cells exhibited significantly larger polygon areas compared to untreated or DMSO vehicle controls (Fig. 5A, B) indicating that the SER network had become less dense. By contrast, treatment with the V-ATPase inhibitor bafilomycin A1 (50 nM for 16 h) resulted in significantly smaller SER polygons (Fig. 5A, B, Additional file 1: Fig. S3) and expansion of the dense RER (Fig. 5A, C, Additional file 1: Fig. S3) indicating an increase in the density of the ER network throughout the cell. Interestingly, drugs known to induce the ER stress response, thapsigargin (100 nM for 16 h) and tunicamycin (500 nM for 16 h) did not result in statistically significant alteration to the organisation of the ER network, at least at the resolution afforded by confocal fluorescence microscopy (Fig. 5).
Analysis of ER distribution in U-2 OS cells using various ER markers
In order to validate that our analysis pipeline was compatible with alternative ER markers, we examined several over-expressed fluorescently-tagged ER-resident proteins, antibodies against ER residents, and the live-cell stain ER Tracker™ (Fig. 6). We selected four ER proteins, each with different functions, tagged with yellow fluorescent protein (YFP) from a library of expression constructs previously characterised in a large-scale subcellular localisation project [50, 51]. The conserved protein Golgi Transport 1B (GOLT1B) has been reported to localise to ER exit sites (ERES) and cis Golgi cisternae in yeast , and functions in mediating anterograde transport between the two organelles in rice endosperm cells . Although the human orthologue has not been functionally described, a GOLT1B-YFP fusion protein has also been shown to localise to the ER in HeLa cells , and in our experiments in U-2 OS cells the GOLT1B-YFP protein localised to tubular structures with a high resemblance to SER tubules and dense perinuclear ER (Fig. 6A), suggesting that it is a bona fide ER marker. Similarly, we observed a high signal of the ER lipid-raft associated protein (ERLIN2) , the E3 ubiquitin-protein ligase synoviolin (SYVN1) , and a subunit of the oligosaccharyltransferase complex A (OST-A) Magnesium Transporter 1 (MAGT1)  YFP-tagged proteins on ER structures in U-2 OS cells (Fig. 6A). We adapted our previously described pipeline designed for Sec61β-mEmerald for each of these markers (Additional file 2: Appendix 1, Table 2) and successfully obtained metrics for ER polygon region area and % dense ‘RER’ in proportion to cell area (Fig. 6B, C).
Finally, we wanted to assess the applicability of our method beyond the use of over-expressed fluorescently-tagged proteins. To this end, we first used an immunofluorescence approach, applying our pipeline to immunolabelled ER-resident proteins. Several of the antibodies we tested only gave a low signal-to-noise ratio. We obtained the best results with an antibody against Reep5, although staining in the cell periphery was weak and did not show continuous tubules (Fig. 6D) precluding analysis of the polygon area. However, we were able to determine the proportion of dense ‘RER’ in relation to cell area using the anti-Reep5 antibody, as it provided a strong stain in the perinuclear area (Fig. 6D, E). Similarly, the frequently used ER Tracker dye failed to strongly label SER tubules but enabled analysis of dense ‘RER’ structures using our pipeline (Fig. 6D, E). Together, these results indicate that our analysis pipeline is compatible with a range of ER markers, while they also suggest caution in their choice, considering the advantages and limitations of each individual maker.
In this study we present an automated image analysis pipeline for the quantitative assessment of the ER using high-content fluorescence microscopy. Until now, the majority of studies reporting on ER morphology have been qualitative; i.e. describing the organelle as ‘normal’ or ‘abnormal’. Qualitative studies are limited to the visual detection of gross changes, such as may occur under high levels of stress or cell toxicity, but often fail to detect more subtle, yet physiologically relevant, changes. When quantitative analysis has been undertaken, it has been manually applied , an approach which is so labour intensive that it precludes its application to high-throughput screening as would be required to identify novel targets which could modify disease-associated ER alterations. Here we present two novel quantitative features that describe the ER, namely the SER tubule polygon area and the proportion of dense ‘RER’ in the cell. Together, these features enable key structures of the SER and RER to be reproducibly quantified in adherent mammalian cell lines.
Our experiments revealed that a cell line stably expressing the ER marker Sec61β-mEmerald provided the highest signal-to-noise ratio which facilitated our analysis. We were also able to demonstrate that our analysis pipeline can also be applied to other cell lines (HeLa Kyoto), as well as to cells stained with different classes of ER markers. Labelling of the ER using over-expression of tagged ER-resident proteins is a widely used strategy to assess ER organisation, however altered ER protein expression can itself disrupt ER organisation. For example, overexpression of the ER-shaping protein Lunapark leads to formation of unbranched SER tubules while its depletion leads to expansion of RER . Consistent with this, our method detected changes to SER and RER organisation resulting from overexpression of the ER-resident proteins GOLT1B  and ERLIN2 . Furthermore, our analysis determined that overexpression of the E3 ubiquitin ligase SYVN1 results in high SER polygon area in U-2 OS cells. This independently supports a recent study showing that overexpression of SYVN1 disrupts the SER network via ubiquitination-mediated inhibition of Atlastin-1 . Therefore, our approach provides an opportunity to screen for proteins that regulate ER morphology and therefore potentially ER function.
The major strength of this automated approach is that it enables the rapid screening of genetic or chemical agents that can modify the ER as these may offer novel therapeutic targets for pathological conditions such as the motor neuron disorders HSPs, which are characterised by disrupted ER organisation. Here we conducted a pilot study on several drugs which have been linked to ER structure and function to validate the use of this pipeline to detect alterations in ER organisation. Thapsigargin and tunicamycin are highly efficient at inducing ER stress. Though ER stress and ER organisation have been linked in many studies, recent evidence suggests that treatment with either of these drugs alone is not sufficient to disrupt ER morphology [58, 59]. Likewise, here we report that treatment with thapsigargin or tunicamycin does not significantly alter tubular SER polygon area or the proportion of dense RER within U-2 OS cells. Lipid availability has been linked to ER morphology  and treatment with the saturated free fatty acid palmitate induces ER expansion and stress in several cell types [32, 60,61,62,63]. Here we show that treatment with a low concentration of palmitic acid increases the area of SER polygons. Similarly, actin filaments have a conserved role in maintenance of the ER network [35, 64, 65] and treatment with the actin‐depolymerising drug cytochalasin B alters ER content in developing oocytes (resulting in ER accumulation at the cortical region during oocyte activation) [37, 66]. Consistent with this, our pipeline was able to detect that brief treatment with cytochalasin B induces ER reorganisation with a significant increase in the area of SER polygons. ER structure and function are maintained by selective autophagic turnover, ER-phagy, impairment of which causes abnormal expansion of RER sheets . Here we show that treatment with the autophagy inhibitor Bafilomycin A1 significantly reduces the area of SER polygons and increases the proportion of dense RER.
Together, these findings support the validity of this automated pipeline for the analysis of ER features. This approach can be incorporated into large-scale image-based phenotypic screens, which could be used to identify genetic or chemical agents that modify the ER in normal or disease model cells.
Human bone osteosarcoma epithelial cells (U2-OS) stably expressing Sec61β-mEmerald , generously provided by Craig Blackstone, were cultured with 1 g/L glucose Dulbecco's modified Eagle's medium (DMEM) (Lonza LZBE12-707F) containing 10% heat inactivated foetal calf serum (Sigma-Aldrich F9665), 0.8% L-Glutamine (Thermo Scientific™ 25030024) and 10 μg/ml (or 0.1%) G-418 (Gibco 10131035), in a humidified incubator at 37 °C and 5% CO2. U-2 OS (ATCC® HTB-96™) and HeLa Kyoto (RRIC:CVCL_0042) cell lines used in transient transfection experiments were maintained in the same conditions but without antibiotic.
U-2 OS cells stably expressing Sec61β-mEmerald  were seeded at 4100 cells per well (equivalent to 1.3 × 104 cells per cm2) in CellCarrier Ultra96-well optically clear plates (PerkinElmer). 16-h treatments started 24 h after seeding while 30-min treatments started at 39 h 30 min after seeding. Staggered treatments allowed synchronised fixation and staining of the entire plate. Cells were treated for 16 h with 0.5% DMSO (Sigma-Aldrich D2650), 250 μM Palmitic acid (Sigma-Aldrich P0500), 50 nM Bafilomycin A1 (B1793 Sigma-Aldrich), 100 nM Thapsigargin (Sigma-Aldrich T9033), and 500 nM Tunicamycin (Sigma-Aldrich T7765); and for 30 min with 10 μM Cytochalasin B (Sigma-Aldrich C6762).
Transient transfection of ER-localising proteins
Non-transfected U-2 OS (RRID:CVCL_0042) and HeLa Kyoto (RRID:CVCL_1922) cells were seeded at 3500 cells per well (equivalent to 1.1 × 104 cells per cm2) in CellCarrier Ultra96-well optically clear plates (PerkinElmer). 24 h after seeding, transfection complexes were prepared using 50 ng per well of various DNA constructs, using 0.13 μl of TransIT-X2 (Mirus MIR6003) in a final reaction volume of 20 μl OptiMEM. After 20 min incubation at room temperature, transfection reactions were added to the cells containing growth medium. The YFP-tagged constructs for transfection of the U-2 OS cells were obtained from a library of expression constructs described previously [50, 51]. For experiments not involving transient transfection, U-2 OS cells were instead stained with ER-Tracker™ (ThermoFisher E34251) or immunolabeled with anti-Reep5 antibodies (Proteintech 14643-1-AP) at the end of the culture time. HeLa Kyoto cells were transfected with Sec61β-mEmerald (Addgene plasmid #90992) . The expression time was 24 h after addition of transfection complexes in all cases.
Sample preparation for imaging
Following the various treatments, cells were incubated at 37 °C for 10 min with a pre-warmed solution of 10 μM Hoechst 33342 (Thermo Scientific™ 62249) and 1:1000 dilution of CellMask™ Deep Red plasma membrane stain (Invitrogen™ C10046) in DMEM medium to counterstain the nuclei and plasma membrane, respectively. In the case of ER Tracker™ stained cells, a 1 μM solution was used (ThermoFisher E34251). Next, cells were washed with phosphate-buffered saline (PBS) buffer, pH 7.4 and fixed with 4% (w/v) paraformaldehyde (PFA) for 15 min at room temperature. Subsequently, cells were washed 3 more times with PBS solution and kept in this solution for imaging. In the case of wells immunostained with anti-Reep5 antibodies, extra steps were followed, which comprised permeabilisation of the cell membrane with 0.01% Triton-X (Alfa Aeser A16046) in PBS for 10 min at room temperature; 30-min blocking with 0.5% FCS (Sigma-Aldrich F9665) in PBS; 1-h incubation with anti-Reep5 antibodies (Proteintech 14643-1-AP) in blocking solution; followed by 3 washes with PBS; 1-h incubation with anti-rabbit secondary Alexa 568 antibodies (ThermoFisher A11036) in blocking solution; and 3 final washes with PBS, after which cells were kept in this solution for imaging.
Automated confocal imaging and analysis
Confocal images of 5 optical slices interspaced by 0.5 μm across 7 × 7 fields of view with a 10% overlay were acquired per well on an Opera Phenix High Content Screening System (PerkinElmer) using a 63 × /1.15 NA water-immersion objective with effective xy resolution of 0.28 μm. Images were analysed using Harmony® v.4.8 high-content analysis software (PerkinElmer), for pipeline see Additional file 2: Appendix 1. We used a computer powered by an Intel Xeon CPU at 3.60 GHz that contains 80 GB RAM, and runs a 64-bit operating system/Windows 10 Pro. Approximated computing time per 7 × 7 fields acquired per well was 10 min. A set of sample images (https://osf.io/nru7g/) was exported and used to recreate the analysis on CellProfiler v.3.1.8 , see Additional file 2: Appendix 2 for pipeline. Specific parameters require optimisation, as indicated Additional file 2: Appendixes 1 and 2, but once inputted the pipelines run to completion for the dataset provided. Cell segmentation was achieved using the CellMask™ channel. Throughout this paper, we quantified and analysed RER and SER morphology across populations of cells, analysing ≥ 50 cells per well. This ensured that inherent variability across a population of cultured cells was accounted for through the high numbers of cells analysed between replicate experiments (compare Fig. 5B,C with Additional file 1: Fig. S3). Raw data files containing per object and per cell results for experiments described within this work are available online at: https://osf.io/nru7g/.
Availability of data and materials
The pipelines and sample raw images analysed during this study, in addition to sample raw data are available in the OSF repository, https://osf.io/nru7g/. The other example image data generated during this study are available from the corresponding author on request.
Dulbecco's modified Eagle's medium
ER Exit Sites
ER lipid-raft associated protein
Golgi Transport 1B
Hereditary Spastic Paraplegias
Magnesium Transporter 1
Oligosaccharyltransferase complex A
Region of interest
Human bone osteosarcoma epithelial cells
Bougen-Zhukov N, Loh SY, Lee HK, Loo LH. Large-scale image-based screening and profiling of cellular phenotypes. Cytometry A. 2017;91(2):115–25.
McQuin C, Goodman A, Chernyshev V, Kamentsky L, Cimini BA, Karhohs KW, et al. Cell Profiler 3.0: next-generation image processing for biology. PLoS Biol. 2018;16(7):e2005970.
Caicedo JC, Cooper S, Heigwer F, Warchal S, Qiu P, Molnar C, et al. Data-analysis strategies for image-based cell profiling. Nat Methods. 2017;14(9):849–63.
Hussain S, Le Guezennec X, Yi W, Dong H, Chia J, Yiping K, et al. Digging deep into Golgi phenotypic diversity with unsupervised machine learning. Mol Biol Cell. 2017;28(25):3686–98.
Galea G, Simpson JC. High-content screening and analysis of the golgi complex. Methods Cell Biol. 2013. https://doi.org/10.1016/B978-0-12-417164-0.00017-3.
Lin D, Lin Z, Cao J, Velmurugan R, Sally Ward E, Ober RJ. A two-stage method for automated detection of ring-like endosomes in fluorescent microscopy images. PLoS ONE. 2019;14(6):e0218931.
Iannetti EF, Smeitink JAM, Beyrath J, Willems PHGM, Koopman WJH. Multiplexed high-content analysis of mitochondrial morphofunction using live-cell microscopy. Nat Protoc. 2016;11(9):1693–710.
Yaothak J, Simpson JC, Heffernan LF, Tsai Y-S, Lin C-C. A semi-automated 2D segmentation and classification system to quantify and characterize morphological features of Golgi-derived membrane structures. Am J Biomed Eng. 2020. https://doi.org/10.5923/j.ajbe.20201002.02.
Öztürk Z, O’Kane CJ, Pérez-Moreno JJ. Axonal endoplasmic reticulum dynamics and its roles in neurodegeneration. Front Neurosci. 2020. https://doi.org/10.3389/fnins.2020.00048.
Wilkinson S. ER-phagy: shaping up and destressing the endoplasmic reticulum. FEBS J. 2019;286(14):2645–63.
Westrate LM, Lee JE, Prinz WA, Voeltz GK. Form follows function: the importance of endoplasmic reticulum shape. Annu Rev Biochem. 2015;84:791–811.
Goyal U, Blackstone C. Untangling the web: mechanisms underlying ER network formation. Biochim Biophys Acta - Mol Cell Res. 2013;1833(11):2492–8.
Peotter J, Kasberg W, Pustova I, Audhya A. COPII-mediated trafficking at the ER/ERGIC interface. Traffic. 2019;20(7):491–503.
Fowler PC, Garcia-Pardo ME, Simpson JC, O’Sullivan NC. NeurodegenERation: the central role for ER contacts in neuronal function and axonopathy, lessons from hereditary spastic paraplegias and related diseases. Front Neurosci. 2019. https://doi.org/10.3389/fnins.2019.01051.
Terasaki M, Shemesh T, Kasthuri N, Klemm RW, Schalek R, Hayworth KJ, et al. Stacked endoplasmic reticulum sheets are connected by helicoidal membrane motifs. Cell. 2013;154:285–96. https://doi.org/10.1016/j.cell.2013.06.031.
Terasaki M. Axonal endoplasmic reticulum is very narrow. J Cell Sci. 2018;131(4):jcs210450.
Schroeder LK, Barentine AES, Merta H, Schweighofer S, Zhang Y, Baddeley D, et al. Dynamic nanoscale morphology of the ER surveyed by STED microscopy. J Cell Biol. 2019;218(1):83–96.
Zhang H, Hu J. Shaping the endoplasmic reticulum into a social network. Trends Cell Biol. 2016;26(12):934–43.
Farías GG, Fréal A, Tortosa E, Stucchi R, Pan X, Portegies S, et al. Feedback-driven mechanisms between microtubules and the endoplasmic reticulum instruct neuronal polarity. Neuron. 2019. https://doi.org/10.1016/j.neuron.2019.01.030.
Shibata Y, Voeltz GK, Rapoport TA. Rough sheets and smooth tubules. Cell. 2006;126:435–9. https://doi.org/10.1016/j.cell.2006.07.019.
Liu S, Pellman D. The coordination of nuclear envelope assembly and chromosome segregation in metazoans. Nucleus. 2020;11:35–52.
Puhka M, Joensuu M, Vihinen H, Belevich I, Jokitalo E. Progressive sheet-to-tubule transformation is a general mechanism for endoplasmic reticulum partitioning in dividing mammalian cells. Mol Biol Cell. 2012;23:2424–32.
Terasaki M, Chen LB, Fujiwara K. Microtubules and the endoplasmic reticulum are highly interdependent structures. Cell. 2007;103:1557–68.
Voeltz GK, Prinz WA, Shibata Y, Rist JM, Rapoport TA. A class of membrane proteins shaping the tubular endoplasmic reticulum. Cell. 2006;124(3):573–86.
Ulloa G, Hamati F, Dick A, Fitzgerald J, Mantell J, Verkade P, et al. Lipid species affect morphology of endoplasmic reticulum: a sea urchin oocyte model of reversible manipulation. J Lipid Res. 2019;60(11):1880–91.
Romero-Brey I, Bartenschlager R. Endoplasmic reticulum: the favorite intracellular niche for viral replication and assembly. Viruses. 2016;8(6):160.
Mohd Ropidi MI, Khazali AS, Nor Rashid N, Yusof R. Endoplasmic reticulum: a focal point of Zika virus infection. J Biomed Sci. 2020;27(1):1–13.
Jones JR, Kong L, Hanna MG, Hoffman B, Krencik R, Bradley R, et al. Mutations in GFAP disrupt the distribution and function of organelles in human astrocytes. Cell Rep. 2018;25(4):947–58.
Mookherjee D, Majumder P, Mukherjee R, Chatterjee D, Kaul Z, Das S, et al. Cytosolic aggregates in presence of non-translocated proteins perturb endoplasmic reticulum structure and dynamics. Traffic. 2019;20(12):943–60.
Sun S, Shi G, Han X, Francisco AB, Ji Y, Mendonça N, et al. Sel1L is indispensable for mammalian endoplasmic reticulum-associated degradation, endoplasmic reticulum homeostasis, and survival. Proc Natl Acad Sci USA. 2014;111(5):E582–91.
Lee JE, Yang YM, Liang FX, Gough DJ, Levy DE, Sehgal PB. Nongenomic STAT5-dependent effects on Golgi apparatus and endoplasmic reticulum structure and function. Am J Physiol - Cell Physiol. 2012;302(5):C804–20.
Wikstrom JD, Israeli T, Bachar-Wikstrom E, Swisa A, Ariav Y, Waiss M, et al. AMPK regulates ER morphology and function in stressed pancreatic β-cells via phosphorylation of DRP1. Mol Endocrinol. 2013;27(10):1706–23.
Zheng P, Chen Q, Tian X, Qian N, Chai P, Liu B, et al. DNA damage triggers tubular endoplasmic reticulum extension to promote apoptosis by facilitating ER-mitochondria signaling. Cell Res. 2018;28(8):833–54.
Bik E, Mielniczek N, Jarosz M, Denbigh J, Budzynska R, Baranska M, et al. Tunicamycin induced endoplasmic reticulum changes in endothelial cells investigated: In vitro by confocal Raman imaging. Analyst. 2019;144(22):6561–9.
Joensuu M, Belevich I, Rämö O, Nevzorov I, Vihinen H, Puhka M, et al. ER sheet persistence is coupled to myosin 1c-regulated dynamic actin filament arrays. Mol Biol Cell. 2014;25(7):1111–26.
Haynes LW, Weller RO. The effects of cytochalasin B and colchicine on cell motility and ultrastructure in primary cultures of malignant gliomas. Acta Neuropathol. 1978;44(1):21–30.
Feitosa WB, Lopes E, Visintin JA, Assumpção MEODA. Endoplasmic reticulum distribution during bovine oocyte activation is regulated by protein kinase C via actin filaments. J Cell Physiol. 2020;235(7–8):5823–34.
Blackstone C. Converging cellular themes for the hereditary spastic paraplegias. Curr Opin Neurobiol. 2018;51:139–46.
Shemesh T, Klemm RW, Romano FB, Wang S, Vaughan J, Zhuang X, et al. A model for the generation and interconversion of ER morphologies. Proc Natl Acad Sci U S A. 2014;11(49):E5243–51.
Hu J, Shibata Y, Zhu PP, Voss C, Rismanchi N, Prinz WA, et al. A class of dynamin-like GTPases involved in the generation of the tubular ER network. Cell. 2009;138(3):549–61.
Kerselidou D, Dohai BS, Nelson DR, Daakour S, de Cock N, Al Z, et al. Alternative glycosylation controls endoplasmic reticulum dynamics and tubular extension in mammalian cells. Sci Adv. 2021. https://doi.org/10.1126/sciadv.abe8349.
Gaffke L, Pierzynowska K, Rintz E, Cyske Z, Giecewicz I, Węgrzyn G. Molecular sciences gene expression-related changes in morphologies of organelles and cellular component organization in mucopolysaccharidoses. Int J Mol Sci. 2021. https://doi.org/10.3390/ijms22052766.
Man H, Bian H, Zhang X, Wang C, Huang Z, Yan Y, et al. Hybrid labeling system for dSTORM imaging of endoplasmic reticulum for uncovering ultrastructural transformations under stress conditions. Biosen Bioelectron. 2021;189:113378.
Pain C, Kriechbaumer V, Kittelmann M, Hawes C, Fricker M. Quantitative analysis of plant ER architecture and dynamics. Nat Commun. 2019;10:1–15.
English AR, Zurek N, Voeltz GK. Peripheral ER structure and function. Curr Opin Cell Biol. 2009;21(4):596–602.
Zhu Y, Zhang G, Lin S, Shi J, Zhang H, Hu J. Sec61β facilitates the maintenance of endoplasmic reticulum homeostasis by associating microtubules. Protein Cell. 2018;9:616–28. https://doi.org/10.1007/s13238-017-0492-5.
Griffing LR. Networking in the endoplasmic reticulum. Biochem Soc Trans. 2010. https://doi.org/10.1042/BST0380747.
Lee CA, Blackstone C. ER morphology and endo-lysosomal crosstalk: Functions and disease implications. Biochim Biophys Acta Mol Cell Biol Lipids. 2020;1865(1):158544.
Costantini L, Snapp E. Probing endoplasmic reticulum dynamics using fluorescence imaging and photobleaching techniques. Curr Protoc Cell Biol. 2013. https://doi.org/10.1002/0471143030.cb2107s60.
Mehrle A, Rosenfelder H, Schupp I, del Val C, Arlt D, Hahne F, et al. The LIFEdb database in 2006. Nucleic Acids Res. 2006;34(suppl_1):D415–8.
Simpson JC, Wellenreuther R, Poustka A, Pepperkok R, Wiemann S. Systematic subcellular localization of novel proteins identified by large-scale cDNA sequencing. EMBO Rep. 2000;1(3):287–92.
Lorente-Rodríguez A, Heidtman M, Barlowe C. Multicopy suppressor analysis of thermosensitive YIP1 alleles implicates GOT1 in transport from the ER. J Cell Sci. 2009;122(10):1540–50.
Wang Y, Liu F, Ren Y, Wang Y, Liu X, Long W, et al. GOLGI TRANSPORT 1B regulates protein export from the endoplasmic reticulum in rice endosperm cells. Plant Cell. 2016;28(11):2850–65.
Zhang X, Cai J, Zheng Z, Polin L, Lin Z, Dandekar A, et al. A novel ER-microtubule-binding protein, ERLIN2, stabilizes Cyclin B1 and regulates cell cycle progression. Cell Discov. 2015;1(1):1–18.
Zhao Y, Feng Z, Zou Y, Liu Y. The E3 ubiquitin ligase SYVN1 ubiquitinates atlastins to remodel the endoplasmic reticulum network. iScience. 2020. https://doi.org/10.1016/j.isci.2020.101494.
Ramírez AS, Kowal J, Locher KP. Cryo-electron microscopy structures of human oligosaccharyltransferase complexes OST-A and OST-B. Science. 2019. https://doi.org/10.1126/science.aaz3505.
Gao G, ZhuID C, Liu E, NabiID IR. Reticulon and CLIMP-63 regulate nanodomain organization of peripheral ER tubules. PLoS Biol. 2019. https://doi.org/10.1371/journal.pbio.3000355.
Lu M, Van Tartwijk FW, Lin JQ, Nijenhuis W, Parutto P, Fantham M, et al. The structure and global distribution of the endoplasmic reticulum network are actively regulated by lysosomes. Sci Adv. 2020;6(51):eabc7209.
Fattouh N, Cazevieille C, Landmann F. Wolbachia endosymbionts subvert the endoplasmic reticulum to acquire host membranes without triggering ER stress. PLoS Negl Trop Dis. 2018. https://doi.org/10.1371/journal.pntd.0007218.
Borradaile NM, Han X, Harp JD, Gale SE, Ory DS, Schaffer JE. Disruption of endoplasmic reticulum structure and integrity in lipotoxic cell death. J Lipid Res. 2006;47(12):2726–37.
Peng G, Li L, Liu Y, Pu J, Zhang S, Yu J, et al. Oleate blocks palmitate-induced abnormal lipid distribution, endoplasmic reticulum expansion and stress, and insulin resistance in skeletal muscle. Endocrinology. 2011;152(6):2206–18.
Hurt CM, Björk S, Ho VK, Gilsbach R, Hein L, Angelotti T. REEP1 and REEP2 proteins are preferentially expressed in neuronal and neuronal-like exocytotic tissues. Brain Res. 2014;1545:12–22. https://doi.org/10.1016/j.brainres.2013.12.008.
Rizzo AM, Colombo I, Montorfano G, Zava S, Corsetto PA. Exogenous fatty acids modulate ER lipid composition and metabolism in breast cancer cells. Cells. 2021;10(1):175.
Fehrenbacher KL, Davis D, Wu M, Boldogh I, Pon LA. Endoplasmic reticulum dynamics, inheritance, and cytoskeletal interactions in budding yeast. Mol Biol Cell. 2002;13(3):854–65.
Hamada T, Ueda H, Kawase T, Hara-Nishimura I. Microtubules contribute to tubule elongation and anchoring of endoplasmic reticulum, resulting in high network complexity in arabidopsis. Plant Physiol. 2014;166(4):1869–76.
Prodon F, Sardet C, Nishida H. Cortical and cytoplasmic flows driven by actin microfilaments polarize the cortical ER-mRNA domain along the a–v axis in ascidian oocytes. Dev Biol. 2008;313(2):682–99.
Khaminets A, Heinrich T, Mari M, Grumati P, Huebner AK, Akutsu M, et al. Regulation of endoplasmic reticulum turnover by selective autophagy. Nature. 2015;522(7556):354–8.
Nixon-Abell J, Obara CJ, Weigel AV, Li D, Legant WR, Xu CS, et al. Increased spatiotemporal resolution reveals highly dynamic dense tubular matrices in the peripheral ER. Science. 2016. https://doi.org/10.1126/science.aaf3928.
We thank Craig Blackstone for provision of the U-2 OS Sec61β-mEmerald stable cell line and members of the JCS laboratory for advice on the image analysis pipeline design.
This work was supported by the SFI-HRB-Wellcome Trust Biomedical Research Partnership, the Irish Research Council (IRC) and University College Dublin (UCD). NCO’S was supported by a Seed Award in Science (Grant number 202020/Z/16/Z) and a UCD PhD Advance award. MEG-P is supported by an IRC-Government of Ireland PhD studentship (Grant number GOIPG/2018/3011). This work in the JCS laboratory was supported by a research grant from Science Foundation Ireland (SFI) (Grant number 16/RI/3745). None of the funding bodies played any role in the design of the study, nor the collection, analysis, or interpretation of data nor in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
ER analysis pipeline in Cell Profiler. Figure S2. Pipeline for analysis of ER polygon regions and dense perinuclear ER applied to HeLa Kyoto (RRID:CVCL_1922) cells expressing Sec61β-mEmerald. Figure S3. Bafilomycin-induced ER reorganisation in U-2 OS cells.
Analysis sequence for ER analysis with Harmony: Table 1. Sequential building blocks in Harmony. Table 2. Pipeline settings adjusted for different experiments. Appendix 2. Analysis sequence for ER analysis with CellProfiler used in Fig. S1.
About this article
Cite this article
Garcia-Pardo, M.E., Simpson, J.C. & O’Sullivan, N.C. A novel automated image analysis pipeline for quantifying morphological changes to the endoplasmic reticulum in cultured human cells. BMC Bioinformatics 22, 427 (2021). https://doi.org/10.1186/s12859-021-04334-x
- Endoplasmic reticulum morphology
- ER function
- High-content imaging
- Automated image analysis