Structural insights into the cause of human RSPH4A primary ciliary dyskinesia.

Cilia and flagella are evolutionarily conserved eukaryotic organelles involved in cell motility and signaling. In humans, mutations in Radial Spoke Head Component 4A (RSPH4A) can lead to primary ciliary dyskinesia (PCD), a life-shortening disease characterized by chronic respiratory tract infections, abnormal organ positioning, and infertility. Despite its importance for human health, the location of RSPH4A in human cilia has not been resolved, and the structural basis of RSPH4A-/- PCD remains elusive. Here, we present the native three-dimensional structure of RSPH4A-/- human respiratory cilia using samples collected noninvasively from a PCD patient. Using cryo-electron tomography (cryo-ET) and subtomogram averaging, we compared the structures of control and RSPH4A-/- cilia, revealing primary defects in two of the three radial spokes (RSs) within the axonemal repeat and secondary (heterogeneous) defects in the central pair complex. Similar to RSPH1-/- cilia, the radial spoke heads of RS1 and RS2, but not RS3, were missing in RSPH4A-/- cilia. However, RSPH4A-/- cilia also exhibited defects within the arch domains adjacent to the RS1 and RS2 heads, which were not observed with RSPH1 loss. Our results provide insight into the underlying structural basis for RSPH4A-/- PCD and highlight the benefits of applying cryo-ET directly to patient samples for molecular structure determination.

INTRODUCTION

Cilia and flagella are highly conserved, microtubule-based organelles that extend from the surfaces of eukaryotic cells and play critical roles in cell motility, sensing, and signaling. In mammals, motile cilia are required for a number of physiological processes, including embryonic development, cerebrospinal fluid circulation, mucus clearance, and sperm motility (Satir and Christensen, 2007; Drummond, 2012). Ciliary dysfunction causes diverse disorders collectively known as ciliopathies, including the severe inherited disorder PCD (OMIM#244400), also known as immotile cilia syndrome (Fliegauf ; Brown and Witman, 2014). PCD is usually an autosomal recessive, heterogeneous genetic disorder characterized by chronic respiratory disease, organ laterality defects, and infertility (Afzelius, 1976). PCD is caused by mutations in genes involved in the biogenesis, assembly, and/or function of motile cilia (Afzelius and Stenram, 2006; Kennedy ; Escudier ; Knowles ). Although the number of known PCD-associated genes/mutations continues to grow (Zariwala ), approximately one-third of PCD cases have no clear genetic basis (Kurkowiak ; Horani and Ferkol, 2018). Furthermore, abnormalities within the ciliary ultrastructure that contribute to PCD clinical phenotypes remain poorly understood.

Most motile cilia are built on a “9+2″ axonemal core with nine peripheral doublet microtubules (DMTs) surrounding two central singlet microtubules known as the central pair complex (CPC; Figure 1A). The axoneme consists of 96-nm repeat units containing regularly spaced protein complexes, including dynein motors and their regulators. Each axonemal repeat contains four identical outer dynein arms (ODAs) and six unique inner dynein arms (IDAs) that are arranged in two rows along the length of the A-tubule (Figure 1, A–C). The dynein motors “walk” along the neighboring B-tubule to generate interdoublet sliding, which is transferred into bending by interdoublet linkers that resist the sliding (Porter and Sale, 2000). Dynein activity is coordinated to generate the undulating motion typical of cilia (Lin and Nicastro, 2018), which requires regulatory signals to be transmitted from the CPC to the radial spokes (RSs) and to two regulatory hubs: the double-headed I1/f dynein and the nexin–dynein regulatory complex (N-DRC) (Figure 1, A–C; Rupp and Porter, 2003; Wirschell , 2009; Heuser ). In humans, each axonemal repeat contains three RSs that are anchored to the A-tubule and project to the axoneme center, where they interact with the CPC (Figure 1, A and B; Warner and Satir, 1974; Goodenough and Heuser, 1985; Smith and Lefebvre, 1997; Mitchell and Sale, 1999; Omoto ; Nicastro ). Cryo–electron tomography (cryo-ET) has revealed RS heterogeneity; that is, RS1 and RS2 share a similar protein composition and structure, whereas RS3 has a different and mostly unknown proteome (Dymek ; Pigino ; Barber ; Lin ). Chlamydomonas reinhardtii, a biflagellated unicellular alga commonly used as a ciliary model system, has a similar RS arrangement but lacks a full-length RS3, containing instead a shortened RS3S stump (Barber ). Full-length RSs are structurally composed of a stalk, an arch, and a head region, which interacts with the CPC projections (Figure 1C; Pigino ; Barber ; Lin ). Though exact mechanisms remain elusive, RSs are involved in the signaling that regulates ciliary motility (Smith and Lefebvre, 1997; Smith and Yang, 2004), and mutations in RS proteins result in severe motility defects or paralyzed flagella (Warner and Satir, 1974; Piperno ; Mitchell and Sale, 1999).

FIGURE 1:

Structural defects in RSPH4A cilia. (A–C) Schematic cross view of a human motile cilium (A) and zoomed-in cross-sectional (B) and longitudinal (C) views of one of the nine doublet microtubules, highlighting major components of the 96-nm repeat. Colored lines in B indicate the locations of the tomographic slices shown in D and E. (D–G) Longitudinal tomographic slices of averaged 96-nm repeats from control HAE (D, E) and RSPH4A (F, G) cilia. White arrows in F denote the missing distal regions of RS1 and RS2 in RSPH4A cilia. RS3, I1/f, and IDAs a, b, c, and e also exhibit lower densities in RSPH4A cilia than the control. Labels: CPC, central pair complex; At and Bt, A- and B-tubules; RS1/2/3, radial spoke1/2/3; a–e and g, inner dynein arm isoforms; I1, I1 dynein complex; N-DRC, nexin dynein regulatory complex; ODA, outer dynein arm; IDA, inner dynein arm. Panels B, D, and E were generated from previously reported HAE control data (Lin ). Scale bar: D–G, 20 nm.

Biallelic pathogenic variants (termed mutations hereafter) in the gene encoding putative RS head protein RSPH4A are a relatively frequent cause of PCD in Eastern European populations (Castleman ; Zietkiewicz ) and in patients of Puerto Rican descent (Daniels ). However, conventional electron microscopy studies of RSPH4A (as well as RSPH1 and RSPH9) mutant cilia have failed to detect clear RS defects, finding instead that up to 60% of the cilia contained heterogeneous defects in the CPC, such as “9+0″ and “8+1″ axonemes (Castleman ; Zietkiewicz ; Daniels ; Kott ; Knowles ; Onoufriadis ). Electron tomography using chemically fixed and plastic-embedded RSPH4A mutant cilia also revealed CPC defects, but no RS defects (Burgoyne ). Examinations of ciliary motility in nasal biopsies from PCD subjects with RSPH4A mutations have reported a range of motility defects, including reduced cilia beat frequency, dyskinesia, and rotational, instead of planar, beating patterns (Castleman ; Daniels ; Frommer ). In contrast, mutating the orthologous RSP4 gene in C. reinhardtii causes clear loss of RS head structures and complete paralysis of ciliary motion (Piperno ; Yang ; Pigino ; Barber ). These studies raise questions about the location and function of RSPH4A in human cilia as well as the underlying basis for PCD caused by RSPH4A mutations.

Here, we use procedures previously established in our laboratories (Lin ) to collect RSPH4A–/– cilia noninvasively from a PCD patient and characterize their structures. We present the native three-dimensional (3D) structure of RSPH4A–/– cilia, identifying primary defects in RS1 and RS2, as well as secondary CPC and IDA defects. By comparing RSPH4A–/– and RSPH1–/– cilia structures from two PCD patients, we find that the RSPH4A mutation causes more severe RS defects, characterized by both shorter RSs and a higher proportion of affected neighboring structures. Our findings suggest that RSPH4A localizes to the radial spoke heads, similar to RSPH1. However, in contrast to RSPH1, RSPH4A also plays a role in stabilizing the arch domain of RS1 and RS2. Together, these results provide insights into PCD etiology in patients with RSPH4A mutations and lay the foundation for structure-based improvements in diagnosis and treatment of PCD.

RESULTS

Cryo-ET reveals structural defects in RSPH4A cilia

Previous work has characterized the structural defects underlying PCD in a patient with a mutation in RSPH1 (Lin ). However, the structural basis for the phenotypically and clinically more severe RSPH4A–/– PCD has not been investigated (Knowles ). To address this, we cultured control HAE cells and respiratory epithelial cells obtained from a patient with homozygous pathogenic mutations in RSPH4A by nasal scrape biopsy as previously reported (Lin ). We then isolated ciliary axonemes from each sample and visualized their 3D structures using cryo-ET and subtomogram averaging. The 96-nm axonemal repeat of the control HAE cilia contained four dimeric ODAs, six monomeric IDAs (a–e and g), the dimeric I1/f dynein, the N-DRC, and three full-length RSs (Figure 1, D–E; Lin ).

In contrast to the control cilia, we found that RSPH4A–/– cilia consistently lacked densities in the distal RS-head regions of RS1 and RS2 (Figure 1F, white arrows). RS3 remained intact in RSPH4A cilia, although its density was lower than in the control (compare Figure 1, D and F). In addition to the RS head defects, we observed lower densities for IDAs a, b, c, e, and I1/f dynein in the averaged axonemal repeat of the mutant than in control cilia (compare Figure 1, E and G). We also observed heterogeneous CPC defects; that is, in most axonemes, the CPC was completely missing (9+0)–sometimes with one DMT moving to the axoneme center instead ((8+1)+0), and occasionally a partial CPC (9+1) was still present (Supplemental Figure S1).

RSPH4A mutation causes primary defects in the distal regions of RS1 and RS2

To determine which abnormalities were primary defects caused by the RSPH4A mutation rather than secondary defects from destabilization of neighboring structures or the axoneme isolation procedure, we next performed classification analyses (Heumann ) using masks focused on the axonemal complexes with observed defects, RS1-3, IDAc, and the I1 dynein, which appear to be missing or reduced in the subtomogram averages (Figure 2; Supplemental Figures S2 and S3). Automated classification separates the axonemal repeats into classes based on structural similarities, determining whether a structure is missing from all axonemal repeats, missing from only a subset of the repeats, or present but positionally flexible. After careful analyses, we report that the RS head-containing distal regions of RS1 and RS2 were missing in 100% of axonemal repeats of RSPH4A cilia (Figure 2, A–D, white arrows). In addition, RS1 and RS2 were further shortened, meaning the arch and distal RS stalk region were also reduced, in 37% and 30% of repeats, respectively (Figure 2, B and D). Similarly shortened radial spokes were also found in the classification analyses from control and RSPH1 PCD cilia, albeit with much lower frequency (RS1: 8% and 11%, respectively; RS2: 6% and 18%, respectively; see Supplemental Figure S2, Class 2). In contrast, the RS3 head was missing in only about one-third of axonemal repeats (Figure 2, E and F), consistent with the RSPH1 phenotype (Supplemental Figure S2). Classification analyses on individual IDAs in RSPH4A cilia revealed a loss of IDA c and I1/f dynein in 84% and 57% of axonemal repeats, respectively (Supplemental Figure S3, B and D). Consistent with the weaker densities of IDA c and I1/f observed in the initial subtomogram averages (Figure 1, E vs. G), the partial loss of the RS3 head, IDA c, and I1/f in the classification analyses indicates that these defects are likely secondary in contrast to the primary 100% loss of distal RS1 and RS2 in RSPH4A cilia. Alternatively, some of the IDA defects in RSPH4A cilia could also be caused by extraction during the sample preparation, as suggested by our previous classification analyses on axonemal structures from wild-type human and Tetrahymena cilia (Lin ). Taken together, these data strongly suggest that loss of the distal regions of RS1 and RS2 are the primary defects in RSPH4A cilia.

FIGURE 2:

RSPH4A cilia exhibit primary defects in RS1 and RS2. (A–F) Longitudinal tomographic slices of two class averages focusing on RS1 (A and B), RS2 (C and D), and RS3 (E and F) from RSPH4A patient cilia. The location of the slice is the same as that shown in Figure 1D. White arrows highlight the missing distal regions of radial spokes, while blue arrows indicate their presence. Note: the distal regions of RS1 and RS2 were missing in all RSPH4A 96-nm repeats, with remnant stalk regions varying in length between classes, whereas RS3 was only absent in a subset of repeats. Scale bar: A–F, 20 nm.

RSPH4A mutation affects the arch domain of RS1 and RS2

Similar to the RSPH4A cilia presented here, RSPH1 cilia exhibit primary defects in the head domains of RS1 and RS2 (Figure 3, A–F, and Supplemental Figure S2; Lin ). This is consistent with the finding that the Chlamydomonas homologues of both RSPH4A and RSPH1, that is, RSP4 and RSP1, were reported as RS head proteins (Yang ; Kohno ). However, PCD patients with mutations in RSPH4A have more severe symptoms than patients with mutations in RSPH1, suggesting potential differences between RSPH1 and RSPH4A localization and function (Knowles ). To better understand the relative positions of RSPH1 and RSPH4A in human cilia, we performed local alignments focusing on each RS to improve the resolution and therefore the comparison of RS structures between control, RSPH1 PCD, and RSPH4A PCD cilia (Figure 3, G and H; Supplemental Movie S1). Isosurface renderings of the axonemal averages clearly reveal the similarities between the primary defects in the 96-nm repeat of RSPH4A–/– and RSPH1–/– cilia, with densities missing in the distal regions of RS1 and RS2 (Figure 3, D–F). However, the remnant structures of the RS1 and RS2 stalks are approximately 2 nm shorter in RSPH4A–/– cilia than in RSPH1–/–, corresponding to a volume of about 100–150 kDa (Figure 3, G and H; see direct comparison in the last column; Supplemental Movie S1). This suggests that RSPH4A localizes to the RS head domain of RS1 and RS2 close to or reaching into the RS arch region between the stalk and head. These differences in the primary defects between RSPH4A–/– and RSPH1–/– cilia may reflect a functional significance of the arch region in ciliary motility.

FIGURE 3:

RS1 and RS2 defects are more severe in RSPH4A cilia than RSPH1 cilia. (A–C) Longitudinal tomographic slices show the averaged 96-nm repeat from control HAE (A), RSPH1 patient (B), and RSPH4A patient cilia (C; the class average with I1 dynein is shown). White arrowheads indicate loss of RS1/RS2 head regions, while black arrows indicate full length RS3 structures. (D–F) Isosurface renderings of the averaged 96-nm repeat from control (D), RSPH1 (E), and RSPH4A cilia (F; the class average with I1 dynein is shown). (G, H) Isosurface renderings of the averaged radial spokes viewed from the longitudinal back (G) and the proximal side (H). RS1, RS2, and RS3 structures are shown in green, blue, and orange, respectively. Difference maps were obtained by superimposing the control HAE or RSPH1 patient structures (gray) over the RSPH4A patient structures (green, blue, or orange), highlighting the structural defects in the RSPH4A patient cilia and the differences between the two PCD patient defects. Images of control and RSPH1 cilia were reproduced from previously reported data (Lin ). Labels: RS1/2/3, radial spoke1/2/3; a–e and g, inner dynein arm isoforms; I1, I1 dynein complex; N-DRC, nexin dynein regulatory complex; ODA, outer dynein arm; IDA, inner dynein arm. Scale bar: A–C, 20 nm.

Movie S1

Animated visualization comparing averaged axonemal repeats from control HAE, RSPH1-/-, and RSPH4A-/- cilia. Each isosurface rendering begins from a back view and rotates clockwise 360 degrees around its longitudinal axis. The first part of the video displays (from left to right) RS1 renderings from control, RSPH1-/-, and RSPH4A-/- cilia (green) and difference maps obtained by superimposing control or RSPH1-/- structures (gray) over the RSPH4A-/- patient structures (green). The second part of the video includes difference maps between RSPH1-/- (gray) and RSPH4A-/- structures of RS1 (green), RS2 (blue), and RS3 (orange). Images of control and RSPH1-/- cilia were generated from previously reported data (Lin, et al., 2014).

DISCUSSION

Diagnosis and treatment of PCD is challenging due to the heterogeneity of underlying genetic defects and their resulting phenotypes. Understanding abnormalities in the ciliary ultrastructure of PCD-associated mutations using conventional methods, which rely on model organisms that do not fully recapitulate human phenotypes, as well as chemical fixation methods that may damage ciliary structures, has been difficult. However, directly studying axonemal structure and function within human tissues presents its own challenges, particularly with preparation and collection of sufficient quantities of high-quality patient sample material. Our previous work has substantially improved methods for collecting and preparing axonemes from human respiratory cells for cryo-ET analysis by optimizing noninvasive cell sampling, in-culture expansion and differentiation, and a gentler method of cilia isolation (Lin ). These advances, combined with the application of robust image processing and classification methods (Heumann ), make it possible to reconstruct high-quality human axonemal structures and to detect and dissect the heterogeneous structural defects observed in patient samples. The latter is critical to distinguish sample variation and preparation artifacts from biologically meaningful defects resulting from genetic mutations.

Here, we show that the primary defect of the RSPH4A mutation is the complete lack of the distal regions, that is, of the entire RS head and at least part of the arch, of RS1 and RS2. In contrast, the distal region of RS3 was reduced in only one-third of the axonemal repeats, demonstrating that RSPH4A is not a RS3 subunit (Figure 2, E and F). This is consistent with previous conclusions that RS3 has a different subunit composition than RS1 and RS2 (Dymek ; Pigino ; Barber ; Lin ). Therefore, the observed RS3 defects are secondary and likely caused by structural destabilization of the RS3 head without its connection to the head of RS2 (Lin ). A similar, partial reduction of the RS3 distal region was observed in RSPH1 cilia that also lack the RS heads of RS1 and RS2 (Supplemental Figure S2L; Lin ). The RS2–RS3 head connection is not only present in control human cilia (Figure 3D), but was also observed in all previously studied organisms with three full-length ciliary radial spokes (e.g., Tetrahymena cilia and sea urchin sperm flagella; Lin ). In addition to the RS3 defect, we observed heterogeneous secondary CPC defects in both RSPH4A and RSPH1 cilia with disrupted RS1 and RS2 heads, suggesting a common underlying mechanism in which RS head interactions with CPC projections are critical for assembly and/or maintenance of the CPC (Oda ).

Although the overall axonemal structure of motile cilia is broadly conserved between species, cryo-ET studies have also revealed species-specific structural differences in major axonemal complexes, making the correlation between structure and axonemal function more challenging. Within the RSs, in particular, the head regions are structurally dichotomous between unicellular protists (e.g., Chlamydomonas, Tetrahymena) and multicellular metazoans, including humans (Lin , 2014; Pigino ). In protists, the RS heads (especially of RS1 and RS2) are considerably larger, with multiple projections that form an extended surface for interactions with CPC projections, whereas the RS heads in metazoans are two small parallel densities resembling ice skates (Lin ). Furthermore, in Chlamydomonas, which is commonly used as a model organism to study cilia structure and motility, the distal half of RS3 is completely absent (Lin ). These observations warrant caution in using Chlamydomonas or other nonhuman species as cellular models for PCD and other ciliopathies. Interestingly, mutating RSP4, the Chlamydomonas ortholog of RSPH4A, results in complete paralysis of Chlamydomonas flagella, whereas RSPH4A mutations in humans result in milder and less penetrant motility defects (Piperno ; Pigino ; Barber ; Knowles ). The increased severity of the RSP4 mutant motility defects in Chlamydomonas compared with human cilia suggests that the presence of a full-length RS3 in human cilia can mitigate some of the deleterious effects of the RSPH4A mutation on ciliary motility. In our analyses of human RSPH4A, RS3 remained intact in approximately two-thirds of axonemal repeats, which may functionally compensate for the loss of RSPH4A in human ciliary motility, though not completely. Interestingly, the angle of the RS3 head compared with its stalk is slightly altered in the RSPH1, but not RSPH4A cilia, compared with the control (compare Figure 2E, Supplemental Figures S2I and S2K). It is unclear whether this angular difference reflects a functional significance. Previous studies demonstrated that the RS3 proteome differs from that of RS1 and RS2 (Dymek ; Pigino ; Barber ; Lin ), but the actual protein composition of the stalk and head region of RS3 and specific protein interactions between RS2 and RS3 are so far unknown. The difference in RS3 head tilt between RSPH1 and RSPH4A cilia could mean that one of these proteins either is present in or interacts with RS3. However, future studies are needed to investigate the RS3 proteome and the mechanisms by which RS3 contributes to ciliary motility.

Classical EM studies of RPSH4A mutant cilia have reported a heterogeneous set of CPC defects in about 50% of the cilia, whereas the remaining cilia appeared structurally intact (Castleman ; Zietkiewicz ; Daniels ; Kott ; Burgoyne ; Knowles ; Onoufriadis ). Despite their apparent structural integrity, however, all RSPH4A mutant cilia exhibited lower mean ciliary beat frequencies as well as abnormal waveforms, with decreased range of motion and lack of coordinated movement (dyskinesis) between adjacent cilia (Daniels ). Our cryo-ET data suggest that these motility defects are in fact due to RS head defects that were previously not visible in chemically fixed EM samples. RSPH4A is either 600 or 716 amino acids in length, depending on isoform, and it contains no known functional domains (Yang ). Although it is possible that RSPH4A regulates ciliary motility through an unknown signaling function, our data strongly support the conclusion that RSPH4A is critically important for the assembly of the head and likely the arch domains of RS1 and RS2, as well as assembly and/or maintenance of a regular CPC complex with two singlet microtubules. The RS heads interact with the CPC projections, and the integrity of RSs, as well as RS-CPC interactions, is critical for stable CPC assembly and function, as they distribute and resist internal mechanical forces during bending that cause local cilia compression (Oda ). Cilia from mice deficient in RSPH4A showed similar CPC defects, and were more sensitive to the destabilizing effects of Taxol (Shinohara ). Interestingly, in the Rsph4a mouse model, the radial spoke heads were missing in all 3 radial spokes (Yoke ). In human RSPH4A cilia, RS1 and RS2 heads were completely missing, but the RS3 heads were missing from only one-third of the axonemal repeats. This is higher than in control cilia (9%) and suggests that the RS2–RS3 head connection is important for RS3 head stability. Although we cannot exclude that RS3 from mouse and human respiratory cilia may contain different proteins, it is more likely that the reduction of RS3 heads in Rsph4a mouse cilia is also a secondary effect caused by the primary RS2 head defect. The observed heterogeneous CPC abnormalities suggest that when the RS heads are not intact, the increased central cylinder space within the axoneme could cause abnormal assembly or arrangement between the CPC and the DMTs in some, but not necessarily all of the patient’s cilia. Given these findings, PCD diagnosis could be improved by directly visualizing the primary mutant defects in the axonemal structure instead of relying on secondary defects in CPC morphology and functional phenotypes, which may be heterogeneous (Supplemental Figure S1; Daniels ; Knowles ).

PCD patients with RSPH1 mutations show milder clinical manifestations than those with RSPH4A mutations, as characterized by a lower incidence of neonatal respiratory distress, a later onset of daily year-round wet cough, a delayed onset of pulmonary symptoms, significantly higher levels of nasal nitric oxide, and better overall lung function (Knowles ). Similar findings have been reported in mice, in which Rsph1 animals exhibit a low level of mucociliary clearance and a less-severe PCD phenotype than those lacking dynein axonemal intermediate chain 1 (Yin ). In contrast, patients with RSPH4A mutations exhibit symptoms similar to those of PCD patients with “classic” dynein arm defects, including a high prevalence of neonatal respiratory distress and low nasal nitric oxide (Noone ). This raises the question: why do RSPH4A mutations lead to poorer PCD outcomes than mutations in RSPH1? Our comparative analyses suggest that RPSH4A cilia exhibit more severe RS defects, with the remnant of the stalk structures of RS1 and RS2 appearing shorter than those of RSPH1 cilia in a larger proportion of axonemal repeats (Figures 2 and 3, G and H; Supplemental Figure S2; Supplemental Movie S1). These data support a model in which RSPH4A is more centrally located within the RS head than RSPH1, potentially reaching into the arch region or interacting strongly with arch proteins, in agreement with recent findings (Zheng ; Grossman-Haham ; Gui ). This is consistent with IHC data that demonstrate loss of RSPH4A results in the loss of RSPH1; while the loss of RSPH1 does not result in the loss of RSPH4A (Frommer ). Of note, a recent study using chemical cross-linking followed by mass spectrometry in Chlamydomonas has established putative interactions between RSPH4A homologue RSP4 and arch protein RSP2 (Poghosyan ). Therefore, arch destabilization due to RSPH4A loss may cause RSPH4A PCD to present a more severe clinical phenotype than RSPH1 PCD (Knowles ). In addition, the differences in the extent and location of RS defects clearly correspond to increased severity and prevalence of CPC defects; classical EM studies of patient cilia showed the axonemal organization deviating from the standard 9+2 arrangement in approximately 50% of RSPH4A mutant cilia compared with 20% of RSPH1 mutant cilia (Daniels ; Knowles ). Our cryo-ET studies of isolated axonemes also show that 77% of RSPH4A cilia completely lacked the CPC without any remnant density in the axoneme center (9+0), whereas only 30% in RSPH1 cilia lacked any density in the center (Supplemental Figure S1B). Although usually the RS heads (and not the arches) are interacting with the CPC, it is possible that the remaining arch region of RS1 and RS2 in RSPH1 cilia provides some support to withstand forces exerted during ciliary beating, thus maintaining more of the ciliary functionality. Alternatively, one could invoke a mechanism by which RSPH4A itself or RS arch proteins that require RPSH4A for proper assembly maintain important signaling functions necessary for proper ciliary motility. Future studies might consider replacing RSPH1 and RSPH4A with altered versions (more flexible or rigid, mutations in key residues) to determine their effects on CPC morphology and overall ciliary motility. By better understanding how mutations within axonemal structures contribute to ciliary disease, we can work toward improving the diagnosis and treatment of PCD and other ciliopathies.

MATERIALS AND METHODS

Human nasal epithelial cell culture

Non-PCD, non–cystic fibrosis control human airway epithelial (HAE) cells were obtained from the UNC Cell and Tissue Culture Core Facility under protocols approved by the Institutional Review Board for the Protection of the Rights of Human Subjects at UNC. All human tissues were obtained with informed consent. Cells were cultured as previously described (Fulcher and Randell, 2013). Briefly, human airway epithelial cells (approximately 106) were plated on human type IV collagen–coated (250 μg, Sigma) Millicell culture inserts 4.2 cm2 in area (Millipore). Then an air–liquid interface was established on confluent cells (approximately 3 days after plating) by removing the apical media and feeding from the basolateral side. Cultures were fed three times per week with ALI media (Fulcher and Randell, 2013; 50:50 mixture of LHC-basal and DMEM-H media) for 2 weeks, after which feedings were reduced to twice per week. Along with medium changes, phosphate-buffered saline (PBS) was used to wash the apical surfaces of the cultures to remove mucus and cell debris. Phase contrast microscopy was used to monitor ciliated-cell differentiation, and cultures became heavily ciliated at the air–liquid interface after approximately 4 weeks.

We identified a PCD patient with a classic PCD phenotype, including neonatal respiratory distress despite term birth, along with daily nasal congestion and year-round daily wet cough from infancy, chronic recurrent sinusitis, severe otitis media and chronic mastoiditis, bronchiectasis, and low nasal nitric oxide (21 nl/min). The patient had normal situs and homozygous pathogenic mutations in the RSPH4A gene (c. 921+3-6delAAGT [p.splice]; Daniels ). Nasal epithelial cells were obtained via nasal scrape biopsy under protocols approved by the Institutional Review Board for the Protection of the Rights of Human subjects at UNC (Knowles ). As previously described (Suprynowicz ), the cells were initially expanded on plastic dishes containing fibroblast-conditioned medium and Rho kinase inhibitor (Y-27632). Then, as described above, the expanded cells were cultured at an air–liquid interface on collagen-coated culture inserts until ciliated cell differentiation occurred (∼4 weeks). Consistent with in vivo observations (Daniels ), ciliary beat in the RSPH4A–/– cultures was dyskinetic, and ciliary beat frequency was reduced (3.6 Hz; n = 6). PCR and Western blot analyses confirmed that no full-length mRNA or protein for RSPH4A could be detected in the cultured RSPH4A–/– cells or cilia, respectively.

Isolation of human ciliary axonemes

Axonemes from respiratory cilia were isolated from human nasal epithelial cells as previously described (Lin ). Briefly, well-differentiated cells were gently washed with PBS to remove mucus and cell debris. The air–liquid interface culture dishes were then filled (both apical and basal compartments) with ice-cold PBS and placed on a bed of ice for approximately 5 min to reduce ciliary beating until it nearly ceased. The PBS was then removed from the apical surface of the culture insert and replaced with 50 μl of ice-cold deciliation buffer (10 mM Tris pH 7.5, 50 mM NaCl, 10 mM CaCl2, 1 mM EDTA, 0.1% Triton X-100, 7 mM β-mercaptoethanol, 1% protease inhibitor cocktail [Sigma P8340]) per 12-mm insert. After 2 min incubation without shaking, the apical solution containing cilia was transferred to a microcentrifuge tube and gently centrifuged at 4°C for 1 min at 500 × g to remove mucus and cellular debris. The cilia-containing supernatant was then centrifuged at 4 °C for 5 min at 5000 × g to collect the ciliary axonemes. Finally, the axonemal pellet was gently dispersed with resuspension buffer (30 mM HEPES pH 7.3, 1 mM EGTA, 5 mM MgSO4, 0.1 mM EDTA, 25 mM NaCl, 1 mM dithiothreitol, 1% protease inhibitor cocktail, and 100 gml-1 soybean trypsin inhibitor [Sigma T9128]) and kept on ice until cryo-sample preparation.

Cryo-sample preparation and cryo–electron tomography

Quantifoil holey carbon grids (R2/2, Quantifoil Micro Tools GmbH, Germany) were glow discharged and loaded into a homemade plunge-freezing device. Isolated axoneme samples (3 μl) and fivefold concentrated 10-nm colloidal gold (1 μl, Sigma-Aldrich) were added to the grid. The grid was blotted with filter paper for 1.5–2.5 s and immediately plunge-frozen in liquid ethane that was cooled by liquid nitrogen. The vitrified samples were then stored in liquid nitrogen until data collection by cryo-ET.

Single-axis tilt series were collected using a Tecnai F30 transmission electron microscope (FEI, Hillsboro, OR) equipped with a postcolumn energy filter (Gatan, Pleasanton, CA). Cryo-samples were transferred into the microscope with a cryo-holder (Gatan) and kept below the devitrification temperature (approximately –170°C). The samples were imaged at 300 kV, with –6 or –8 μm defocus, under low-dose conditions and in the zero-loss mode of the energy filter (slit width 20 eV). Axonemes that appeared well preserved and not obviously compressed by EM inspection were selected for cryo–electron tomography. Tilt series images were recorded automatically over a tilt range of –65° to +65° with 1.5–2.5° angular increments using the SerialEM image acquisition software (Mastronarde, 2005). The cumulative electron dose was limited to 100 e/Å2 for individual tilt series. All images were recorded on a 2k × 2k charge-coupled device camera (Gatan) at 13,500× magnification, resulting in a pixel size of approximately 1 nm.

Image processing

Tilt series images with sufficient and evenly distributed gold fiducial markers were aligned and then reconstructed by weighted back projection into tomograms using the IMOD software (Kremer ). Tomograms with compressed axonemes (caused by insufficient ice thickness during cryo-sample preparation) were excluded from further subtomogram averaging and analysis. Subtomograms of axonemal repeats were extracted from the raw tomograms, aligned, and averaged (with missing wedge compensation) using the PEET software (Nicastro ) to enhance the signal-to-noise ratio and improve resolution (Table 1). Two different alignment strategies were employed to optimize alignment accuracy: 1) for the entire 96-nm axonemal repeat, global alignment with relatively large particle volumes was used (Figures 1, 2, and 3, A–F; Supplemental Figures S1, S2, and S3); 2) to resolve the RSs with greater detail, smaller subtomograms centered on the structure of interest were used for local alignment and then averaged (Figure 3, G and H). Subtomogram averaging of positionally flexible structures or structures that are present in only a subset of averaged particles causes blurring out of the structure signal in the average. Therefore, we used the unsupervised classification tool built into the PEET program (Heumann ) to analyze the heterogeneity of axonemal structures that appeared blurred in the subtomogram averages. Briefly, after alignment, soft-edged masks were applied to focus the classification on the structures of interest, that is, areas with lower electron density. The classification uses principal component analysis and clustering to divide the particles (axonemal repeats) into different classes based on structural similarity. If a structure is positionally flexible, the class averages will show the unblurred structure in different positions in the class averages, whereas reduced structures will result in classes that either contain or lack the structure of interest (Figure 2; Supplemental Figures S2 and S3). The resolution of the resulting averages was estimated in a 20 nm3 subvolume in the 96-nm axonemal repeat using the Fourier shell correlation method and the 0.5 criterion (Harauz and van Heel, 1986). The IMOD software (Kremer ) and the UCSF Chimera package (Pettersen ) were used for visualization of tomographic slices and 3D isosurface renderings, respectively.

TABLE 1:

Image processing information for specimens used in this study.

Specimen	Tomograms included	Averaged repeats	Resolution (nm)a	Used in figure(s)
Control ciliab	12	857	3.44	1,3,S1,S2
RSPH1–/– ciliab	20	3044	3.36	3,S1,S2
RSPH4A–/– cilia	13	1637	4.18	1,2,3,S1,S3

The 0.5 Fourier shell correlation criterion was used to estimate the resolution of the averaged repeats at the site of the doublet microtubule.

Data have been previously published (Lin ).

Click here for additional data file.

Click here for additional data file.