Mastigoneme structure reveals insights into the O-linked glycosylation code of native hydroxyproline-rich helices.

Hydroxyproline-rich glycoproteins (HRGPs) are a ubiquitous class of protein in the extracellular matrices and cell walls of plants and algae, yet little is known of their native structures or interactions. Here, we used electron cryomicroscopy (cryo-EM) to determine the structure of the hydroxyproline-rich mastigoneme, an extracellular filament isolated from the cilia of the alga Chlamydomonas reinhardtii. The structure demonstrates that mastigonemes are formed from two HRGPs (a filament of MST1 wrapped around a single copy of MST3) that both have hyperglycosylated poly(hydroxyproline) helices. Within the helices, O-linked glycosylation of the hydroxyproline residues and O-galactosylation of interspersed serine residues create a carbohydrate casing. Analysis of the associated glycans reveals how the pattern of hydroxyproline repetition determines the type and extent of glycosylation. MST3 possesses a PKD2-like transmembrane domain that forms a heteromeric polycystin-like cation channel with PKD2 and SIP, explaining how mastigonemes are tethered to ciliary membranes.

Structural specializations of the sperm tail.

Sperm motility is crucial to reproductive success in sexually reproducing organisms. Impaired sperm movement causes male infertility, which is increasing globally. Sperm are powered by a microtubule-based molecular machine-the axoneme-but it is unclear how axonemal microtubules are ornamented to support motility in diverse fertilization environments. Here, we present high-resolution structures of native axonemal doublet microtubules (DMTs) from sea urchin and bovine sperm, representing external and internal fertilizers. We identify >60 proteins decorating sperm DMTs; at least 15 are sperm associated and 16 are linked to infertility. By comparing DMTs across species and cell types, we define core microtubule inner proteins (MIPs) and analyze evolution of the tektin bundle. We identify conserved axonemal microtubule-associated proteins (MAPs) with unique tubulin-binding modes. Additionally, we identify a testis-specific serine/threonine kinase that links DMTs to outer dense fibers in mammalian sperm. Our study provides structural foundations for understanding sperm evolution, motility, and dysfunction at a molecular level.

Structure of the Decorated Ciliary Doublet Microtubule.

The axoneme of motile cilia is the largest macromolecular machine of eukaryotic cells. In humans, impaired axoneme function causes a range of ciliopathies. Axoneme assembly, structure, and motility require a radially arranged set of doublet microtubules, each decorated in repeating patterns with non-tubulin components. We use single-particle cryo-electron microscopy to visualize and build an atomic model of the repeating structure of a native axonemal doublet microtubule, which reveals the identities, positions, repeat lengths, and interactions of 38 associated proteins, including 33 microtubule inner proteins (MIPs). The structure demonstrates how these proteins establish the unique architecture of doublet microtubules, maintain coherent periodicities along the axoneme, and stabilize the microtubules against the repeated mechanical stress induced by ciliary motility. Our work elucidates the architectural principles that underpin the assembly of this large, repetitive eukaryotic structure and provides a molecular basis for understanding the etiology of human ciliopathies.

Characterizing the Major Structural Variant Alleles of the Human Genome.

In order to provide a comprehensive resource for human structural variants (SVs), we generated long-read sequence data and analyzed SVs for fifteen human genomes. We sequence resolved 99,604 insertions, deletions, and inversions including 2,238 (1.6 Mbp) that are shared among all discovery genomes with an additional 13,053 (6.9 Mbp) present in the majority, indicating minor alleles or errors in the reference. Genotyping in 440 additional genomes confirms the most common SVs in unique euchromatin are now sequence resolved. We report a ninefold SV bias toward the last 5 Mbp of human chromosomes with nearly 55% of all VNTRs (variable number of tandem repeats) mapping to this portion of the genome. We identify SVs affecting coding and noncoding regulatory loci improving annotation and interpretation of functional variation. These data provide the framework to construct a canonical human reference and a resource for developing advanced representations capable of capturing allelic diversity.

Gene dosage of independent dynein arm motor preassembly factors influences cilia assembly in Chlamydomonas reinhardtii.

Motile cilia assembly utilizes over 800 structural and cytoplasmic proteins. Variants in approximately 58 genes cause primary ciliary dyskinesia (PCD) in humans, including the dynein arm (pre)assembly factor (DNAAF) gene DNAAF4. In humans, outer dynein arms (ODAs) and inner dynein arms (IDAs) fail to assemble motile cilia when DNAAF4 function is disrupted. In Chlamydomonas reinhardtii, a ciliated unicellular alga, the DNAAF4 ortholog is called PF23. The pf23-1 mutant assembles short cilia and lacks IDAs, but partially retains ODAs. The cilia of a new null allele (pf23-4) completely lack ODAs and IDAs and are even shorter than cilia from pf23-1. In addition, PF23 plays a role in the cytoplasmic modification of IC138, a protein of the two-headed IDA (I1/f). As most PCD variants in humans are recessive, we sought to test if heterozygosity at two genes affects ciliary function using a second-site non-complementation (SSNC) screening approach. We asked if phenotypes were observed in diploids with pairwise heterozygous combinations of 21 well-characterized ciliary mutant Chlamydomonas strains. Vegetative cultures of single and double heterozygous diploid cells did not show SSNC for motility phenotypes. When protein synthesis is inhibited, wild-type Chlamydomonas cells utilize the pool of cytoplasmic proteins to assemble half-length cilia. In this sensitized assay, 8 double heterozygous diploids with pf23 and other DNAAF mutations show SSNC; they assemble shorter cilia than wild-type. In contrast, double heterozygosity of the other 203 strains showed no effect on ciliary assembly. Immunoblots of diploids heterozygous for pf23 and wdr92 or oda8 show that PF23 is reduced by half in these strains, and that PF23 dosage affects phenotype severity. Reductions in PF23 and another DNAAF in diploids affect the ability to assemble ODAs and IDAs and impedes ciliary assembly. Thus, dosage of multiple DNAAFs is an important factor in cilia assembly and regeneration.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Mapping and characterization of structural variation in 17,795 human genomes.

A key goal of whole-genome sequencing for studies of human genetics is to interrogate all forms of variation, including single-nucleotide variants, small insertion or deletion (indel) variants and structural variants. However, tools and resources for the study of structural variants have lagged behind those for smaller variants. Here we used a scalable pipeline1 to map and characterize structural variants in 17,795 deeply sequenced human genomes. We publicly release site-frequency data to create the largest, to our knowledge, whole-genome-sequencing-based structural variant resource so far. On average, individuals carry 2.9 rare structural variants that alter coding regions; these variants affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Using a computational model, we estimate that structural variants account for 17.2% of rare alleles genome-wide, with predicted deleterious effects that are equivalent to loss-of-function coding alleles; approximately 90% of such structural variants are noncoding deletions (mean 19.1 per genome). We report 158,991 ultra-rare structural variants and show that 2% of individuals carry ultra-rare megabase-scale structural variants, nearly half of which are balanced or complex rearrangements. Finally, we infer the dosage sensitivity of genes and noncoding elements, and reveal trends that relate to element class and conservation. This work will help to guide the analysis and interpretation of structural variants in the era of whole-genome sequencing.

Robust acoustic trapping and perturbation of single-cell microswimmers illuminate three-dimensional swimming and ciliary coordination.

We report a label-free acoustic microfluidic method to confine single, cilia-driven swimming cells in space without limiting their rotational degrees of freedom. Our platform integrates a surface acoustic wave (SAW) actuator and bulk acoustic wave (BAW) trapping array to enable multiplexed analysis with high spatial resolution and trapping forces that are strong enough to hold individual microswimmers. The hybrid BAW/SAW acoustic tweezers employ high-efficiency mode conversion to achieve submicron image resolution while compensating for parasitic system losses to immersion oil in contact with the microfluidic chip. We use the platform to quantify cilia and cell body motion for wildtype biciliate cells, investigating effects of environmental variables like temperature and viscosity on ciliary beating, synchronization, and three-dimensional helical swimming. We confirm and expand upon the existing understanding of these phenomena, for example determining that increasing viscosity promotes asynchronous beating. Motile cilia are subcellular organelles that propel microorganisms or direct fluid and particulate flow. Thus, cilia are critical to cell survival and human health. The unicellular alga Chlamydomonas reinhardtii is widely used to investigate the mechanisms underlying ciliary beating and coordination. However, freely swimming cells are difficult to image with sufficient resolution to capture cilia motion, necessitating that the cell body be held during experiments. Acoustic confinement is a compelling alternative to use of a micropipette, or to magnetic, electrical, and optical trapping that may modify the cells and affect their behavior. Beyond establishing our approach to studying microswimmers, we demonstrate a unique ability to mechanically perturb cells via rapid acoustic positioning.

Rare coding variants in 35 genes associate with circulating lipid levels-A multi-ancestry analysis of 170,000 exomes.

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels.

Author Correction: Exome sequencing of Finnish isolates enhances rare-variant association power.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Exome sequencing of Finnish isolates enhances rare-variant association power.

Exome-sequencing studies have generally been underpowered to identify deleterious alleles with a large effect on complex traits as such alleles are mostly rare. Because the population of northern and eastern Finland has expanded considerably and in isolation following a series of bottlenecks, individuals of these populations have numerous deleterious alleles at a relatively high frequency. Here, using exome sequencing of nearly 20,000 individuals from these regions, we investigate the role of rare coding variants in clinically relevant quantitative cardiometabolic traits. Exome-wide association studies for 64 quantitative traits identified 26 newly associated deleterious alleles. Of these 26 alleles, 19 are either unique to or more than 20 times more frequent in Finnish individuals than in other Europeans and show geographical clustering comparable to Mendelian disease mutations that are characteristic of the Finnish population. We estimate that sequencing studies of populations without this unique history would require hundreds of thousands to millions of participants to achieve comparable association power.

Mutation of CFAP57, a protein required for the asymmetric targeting of a subset of inner dynein arms in Chlamydomonas, causes primary ciliary dyskinesia.

Primary ciliary dyskinesia (PCD) is characterized by chronic airway disease, reduced fertility, and randomization of the left/right body axis. It is caused by defects of motile cilia and sperm flagella. We screened a cohort of affected individuals that lack an obvious axonemal defect for pathogenic variants using whole exome capture, next generation sequencing, and bioinformatic analysis assuming an autosomal recessive trait. We identified one subject with an apparently homozygous nonsense variant [(c.1762C>T), p.(Arg588*)] in the uncharacterized CFAP57 gene. Interestingly, the variant results in the skipping of exon 11 (58 amino acids), which may be due to disruption of an exonic splicing enhancer. In normal human nasal epithelial cells, CFAP57 localizes throughout the ciliary axoneme. Nasal cells from the PCD patient express a shorter, mutant version of CFAP57 and the protein is not incorporated into the axoneme. The missing 58 amino acids include portions of WD repeats that may be important for loading onto the intraflagellar transport (IFT) complexes for transport or docking onto the axoneme. A reduced beat frequency and an alteration in ciliary waveform was observed. Knockdown of CFAP57 in human tracheobronchial epithelial cells (hTECs) recapitulates these findings. Phylogenetic analysis showed that CFAP57 is highly conserved in organisms that assemble motile cilia. CFAP57 is allelic with the BOP2/IDA8/FAP57 gene identified previously in Chlamydomonas reinhardtii. Two independent, insertional fap57 Chlamydomonas mutant strains show reduced swimming velocity and altered waveforms. Tandem mass tag (TMT) mass spectroscopy shows that FAP57 is missing, and the "g" inner dyneins (DHC7 and DHC3) and the "d" inner dynein (DHC2) are reduced, but the FAP57 paralog FBB7 is increased. Together, our data identify a homozygous variant in CFAP57 that causes PCD that is likely due to a defect in the inner dynein arm assembly process.

De Novo Mutations in FOXJ1 Result in a Motile Ciliopathy with Hydrocephalus and Randomization of Left/Right Body Asymmetry.

Hydrocephalus is one of the most prevalent form of developmental central nervous system (CNS) malformations. Cerebrospinal fluid (CSF) flow depends on both heartbeat and body movement. Furthermore, it has been shown that CSF flow within and across brain ventricles depends on cilia motility of the ependymal cells lining the brain ventricles, which play a crucial role to maintain patency of the narrow sites of CSF passage during brain formation in mice. Using whole-exome and whole-genome sequencing, we identified an autosomal-dominant cause of a distinct motile ciliopathy related to defective ciliogenesis of the ependymal cilia in six individuals. Heterozygous de novo mutations in FOXJ1, which encodes a well-known member of the forkhead transcription factors important for ciliogenesis of motile cilia, cause a motile ciliopathy that is characterized by hydrocephalus internus, chronic destructive airway disease, and randomization of left/right body asymmetry. Mutant respiratory epithelial cells are unable to generate a fluid flow and exhibit a reduced number of cilia per cell, as documented by high-speed video microscopy (HVMA), transmission electron microscopy (TEM), and immunofluorescence analysis (IF). TEM and IF demonstrate mislocalized basal bodies. In line with this finding, the focal adhesion protein PTK2 displays aberrant localization in the cytoplasm of the mutant respiratory epithelial cells.

RPGRIP1L helps to establish the ciliary gate for entry of proteins.

Mutations in transition zone genes change the composition of the ciliary proteome. We isolated new mutations in RPGRIP1L (denotated as RPG1 in algae) that affect the localization of the transition zone protein NPHP4 in the model organism Chlamydomonas reinhardtii NPHP4 localization is not affected in multiple new intraflagellar transport (IFT) mutants. We compared the proteome of cilia from wild-type and mutants that affect the transition zone (RPGRIP1L) or IFT (IFT172 and DHC1b) by mass spectrometry. The rpg1-1 mutant cilia show the most dramatic increase in cytoplasmic proteins. These nonciliary proteins function in translation, membrane remodeling, ATP production and as chaperonins. These proteins are excluded in isolated cilia from fla11-1 (IFT172) and fla24-1 (DHC1b). Our data support the idea that RPGRIP1L, but not IFT proteins, acts as part of the gate for cytoplasmic proteins. The rpg1-1 cilia lack only a few proteins, which suggests that RPGRIP1L only has a minor role of in the retention of ciliary proteins. The fla11-1 mutant shows the greatest loss/reduction of proteins, and one-third of these proteins have a transmembrane domain. Hence, IFT172 may play a role in the retention of proteins.

Chlamydomonas DYX1C1/PF23 is essential for axonemal assembly and proper morphology of inner dynein arms.

Cytoplasmic assembly of ciliary dyneins, a process known as preassembly, requires numerous non-dynein proteins, but the identities and functions of these proteins are not fully elucidated. Here, we show that the classical Chlamydomonas motility mutant pf23 is defective in the Chlamydomonas homolog of DYX1C1. The pf23 mutant has a 494 bp deletion in the DYX1C1 gene and expresses a shorter DYX1C1 protein in the cytoplasm. Structural analyses, using cryo-ET, reveal that pf23 axonemes lack most of the inner dynein arms. Spectral counting confirms that DYX1C1 is essential for the assembly of the majority of ciliary inner dynein arms (IDA) as well as a fraction of the outer dynein arms (ODA). A C-terminal truncation of DYX1C1 shows a reduction in a subset of these ciliary IDAs. Sucrose gradients of cytoplasmic extracts show that preassembled ciliary dyneins are reduced compared to wild-type, which suggests an important role in dynein complex stability. The role of PF23/DYX1C1 remains unknown, but we suggest that DYX1C1 could provide a scaffold for macromolecular assembly.

How Does Cilium Length Affect Beating?

The effects of cilium length on the dynamics of cilia motion were investigated by high-speed video microscopy of uniciliated mutants of the swimming alga, Chlamydomonas reinhardtii. Cells with short cilia were obtained by deciliating cells via pH shock and allowing cilia to reassemble for limited times. The frequency of cilia beating was estimated from the motion of the cell body and of the cilium. Key features of the ciliary waveform were quantified from polynomial curves fitted to the cilium in each image frame. Most notably, periodic beating did not emerge until the cilium reached a critical length between 2 and 4 µm. Surprisingly, in cells that exhibited periodic beating, the frequency of beating was similar for all lengths with only a slight decrease in frequency as length increased from 4 µm to the normal length of 10-12 µm. The waveform average curvature (rad/µm) was also conserved as the cilium grew. The mechanical metrics of ciliary propulsion (force, torque, and power) all increased in proportion to length. The mechanical efficiency of beating appeared to be maximal at the normal wild-type length of 10-12 µm. These quantitative features of ciliary behavior illuminate the biophysics of cilia motion and, in future studies, may help distinguish competing hypotheses of the underlying mechanism of oscillation.

MAPINS, a Highly Efficient Detection Method That Identifies Insertional Mutations and Complex DNA Rearrangements.

Insertional mutagenesis, in which a piece of exogenous DNA is integrated randomly into the genomic DNA of the recipient cell, is a useful method to generate new mutants with phenotypes of interest. The unicellular green alga Chlamydomonas reinhardtii is an outstanding model for studying many biological processes. We developed a new computational algorithm, MAPINS (mapping insertions), to efficiently identify insertion sites created by the integration of an APHVIII (aminoglycoside 3'-phosphotransferase VIII) cassette that confers paromomycin resistance. Using whole-genome sequencing data, this method eliminates the need for genomic DNA manipulation and retains all the sequencing information provided by paired-end sequencing. We experimentally verified 38 insertion sites out of 41 sites (93%) identified by MAPINS from 20 paromomycin-resistant strains. Using meiotic analysis of 18 of these strains, we identified insertion sites that completely cosegregate with paromomycin resistance. In six of the seven strains with a mutant phenotype, we demonstrated complete cosegregation of the mutant phenotype and the verified insertion site. In addition, we provide direct evidence of complex rearrangements of genomic DNA in five strains, three of which involve the APHVIII insertion site. We suggest that strains obtained by insertional mutagenesis are more complicated than expected from previous analyses in Chlamydomonas To map the locations of some complex insertions, we designed 49 molecular markers based on differences identified via whole-genome sequencing between wild-type strains CC-124 and CC-125. Overall, MAPINS provides a low-cost, efficient method to characterize insertional mutants in Chlamydomonas.

Acoustic trap-and-release for rapid assessment of cell motility.

Functional cilia and flagella are crucial to the propulsion of physiological fluids, motile cells, and microorganisms. Motility assessment of individual cells allows discrimination of normal from dysfunctional behavior, but cell-scale analysis of individual trajectories to represent a population is laborious and impractical for clinical, industrial, and even research applications. We introduce an assay that quantifies swimming capability as a function of the variation in polar moment of inertia of cells released from an acoustic trap. Acoustic confinement eliminates the need to trace discrete trajectories and enables automated analysis of hundreds of cells in minutes. The approach closely approximates the average speed estimated from the mean squared displacement of individual cells for wild-type Chlamydomonas reinhardtii and two mutants (ida3 and oda5) that display aberrant swimming behaviors. Large-population acoustic trap-and-release rapidly differentiates these cell types based on intrinsic motility, which provides a highly sensitive and efficient alternative to conventional particle tracing.

Functional characterization of biallelic RTTN variants identified in an infant with microcephaly, simplified gyral pattern, pontocerebellar hypoplasia, and seizures.

BACKGROUND: Biallelic deleterious variants in RTTN, which encodes rotatin, are associated with primary microcephaly, polymicrogyria, seizures, intellectual disability, and primordial dwarfism in human infants. METHODS AND RESULTS: We performed exome sequencing of an infant with primary microcephaly, pontocerebellar hypoplasia, and intractable seizures and his healthy, unrelated parents. We cultured the infant's fibroblasts to determine primary ciliary phenotype. RESULTS: We identified biallelic variants in RTTN in the affected infant: a novel missense variant and a rare, intronic variant that results in aberrant transcript splicing. Cultured fibroblasts from the infant demonstrated reduced length and number of primary cilia. CONCLUSION: Biallelic variants in RTTN cause primary microcephaly in infants. Functional characterization of primary cilia length and number can be used to determine pathogenicity of RTTN variants.

A NIMA-Related Kinase Suppresses the Flagellar Instability Associated with the Loss of Multiple Axonemal Structures.

CCDC39 and CCDC40 were first identified as causative mutations in primary ciliary dyskinesia patients; cilia from patients show disorganized microtubules, and they are missing both N-DRC and inner dynein arms proteins. In Chlamydomonas, we used immunoblots and microtubule sliding assays to show that mutants in CCDC40 (PF7) and CCDC39 (PF8) fail to assemble N-DRC, several inner dynein arms, tektin, and CCDC39. Enrichment screens for suppression of pf7; pf8 cells led to the isolation of five independent extragenic suppressors defined by four different mutations in a NIMA-related kinase, CNK11. These alleles partially rescue the flagellar length defect, but not the motility defect. The suppressor does not restore the missing N-DRC and inner dynein arm proteins. In addition, the cnk11 mutations partially suppress the short flagella phenotype of N-DRC and axonemal dynein mutants, but do not suppress the motility defects. The tpg1 mutation in TTLL9, a tubulin polyglutamylase, partially suppresses the length phenotype in the same axonemal dynein mutants. In contrast to cnk11, tpg1 does not suppress the short flagella phenotype of pf7. The polyglutamylated tubulin in the proximal region that remains in the tpg1 mutant is reduced further in the pf7; tpg1 double mutant by immunofluorescence. CCDC40, which is needed for docking multiple other axonemal complexes, is needed for tubulin polyglutamylation in the proximal end of the flagella. The CCDC39 and CCDC40 proteins are likely to be involved in recruiting another tubulin glutamylase(s) to the flagella. Another difference between cnk11-1 and tpg1 mutants is that cnk11-1 cells show a faster turnover rate of tubulin at the flagellar tip than in wild-type flagella and tpg1 flagella show a slower rate. The double mutant shows a turnover rate similar to tpg1, which suggests the faster turnover rate in cnk11-1 flagella requires polyglutamylation. Thus, we hypothesize that many short flagella mutants in Chlamydomonas have increased instability of axonemal microtubules. Both CNK11 and tubulin polyglutamylation play roles in regulating the stability of axonemal microtubules.