The transcriptional signature associated with human motile cilia

Brief biology and pathophysiology of Tekt bundles

Tektins, a family of microtubule-stabilizing proteins, are critical for cilia and flagella assembly in mammals. They maintain doublet microtubule stability and ciliary/flagellar motility. Loss of Tekt1-5 causes microtubule instability, impaired motility, and diseases like infertility, retinal degeneration, Mainzer-Saldino syndrome, and diabetic nephropathy. Pathophysiological stimuli regulate Tektin expression through transcriptional, posttranscriptional, translational, and posttranslational modifications. This review summarizes the latest findings on Tektin functions and their role in diseases.

The absence of both RIBC1 and RIBC2 induces decreased sperm motility and litter size in male mice

BACKGROUND

RIBC1 (RIB43A domain with coiled-coils 1) and RIBC2 (RIB43A domain with coiled-coils 2) are homolog proteins of RIB43a which is localized to microtubules in the cilia and flagella of unicellular organisms. Cryo-electron microscopy and artificial intelligence studies showed that RIBC1 and RIBC2 are microtubule inner proteins (MIPs) localized in the inner lumen of the doublet microtubules (DMTs) in mouse sperm flagella. However, the function of RIBC1 and RIBC2 in mammalian reproduction and sperm flagella is still unknown.

OBJECTIVE

To clarify the function of RIBC1 and RIBC2 in mouse spermatozoa.

MATERIALS AND METHODS

We generated Ribc1 knockout (KO), Ribc2 KO, and Ribc1 and Ribc2 double-knockout (Ribc1/2 DKO) mice using the CRISPR/Cas9 system and analyzed their phenotypes.

RESULTS

We revealed that the loss of either RIBC1 or RIBC2 alone did not affect male fertility, but the absence of both caused a decrease in pup numbers. Sperm motility analysis showed that Ribc1 KO spermatozoa had reduced velocity, but Ribc2 KO sperm velocities were comparable with WT mice. However, Ribc1/2 DKO sperm velocities were significantly lower than those from Ribc1 KO mice. No structural abnormalities in the axonemal structure at the transmission electron microscope (TEM) level and no abnormalities in the flagellar waveform pattern were observed in Ribc1/2 DKO spermatozoa.

DISCUSSION AND CONCLUSION

Both RIBC1 and RIBC2 are not significant for maintaining the axonemal structure in mouse spermatozoa, but both proteins function cooperatively in sperm motility. This result may indicate that minor structural changes due to RIBC protein absence are not detected at the TEM level, and RIBC2 function depends on RIBC1 in sperm motility. We think that reduced litter size in Ribc1/2 DKO mice is caused by reduced sperm motility due to minor structural abnormalities caused by the loss of two RIBC proteins.

Identification of Candidate Genes and eQTLs Related to Porcine Reproductive Function

Expression quantitative trait locus (eQTL) mapping is an effective tool for identifying genetic variations that regulate gene expression. An increasing number of studies suggested that SNPs associated with complex traits in farm animals are considered as expression quantitative trait loci. Identifying eQTLs associated with gene expression levels in the endometrium helps to unravel the regulatory mechanisms of genes related to reproductive functions in this tissue and provides molecular markers for the genetic improvement of high-fertility sow breeding. In this study, 218 RNA-seq data from pig endometrial tissue were used for eQTL analysis to identify genetic variants regulating gene expression. Additionally, weighted gene co-expression network analysis (WGCNA) was performed to identify hub genes involved in reproductive functions. The eQTL analysis identified 34,876 significant cis-eQTLs regulating the expression of 5632 genes (FDR ≤ 0.05), and 90 hub genes were identified by WGCNA analysis. By integrating eQTL and WGCNA results, 14 candidate genes and 16 fine-mapped cis-eQTLs were identified, including FRK, ARMC3, SLC35F3, TMEM72, FFAR4, SOWAHA, PSPH, FMO5, HPN, FUT2, RAP1GAP, C6orf52, SEL1L3, and CLGN, which were involved in the physiological processes of reproduction in sows through hormone regulation, cell adhesion, and amino acid and lipid metabolism. These eQTLs regulate the high expression of candidate genes in the endometrium, thereby affecting reproductive-related physiological functions. These findings enhance our understanding of the genetic basis of reproductive traits and provide valuable genetic markers for marker-assisted selection (MAS), which can be applied to improve sow fecundity and optimize breeding strategies for high reproductive performance.

Single-cell Transcriptomics Unravel Stereocilia Degeneration as a Key Contributor to Age-related Vestibular Dysfunction in Mice and Humans

Age-related vestibular dysfunction (ARVD) is a prevalent, debilitating condition in the elderly. The etiology and molecular mechanisms are poorly understood. We focused on mechanosensitive hair cells (HCs) as they are particularly vulnerable to aging. Using single-cell RNA-seq transcriptomes of young and old mouse vestibular HCs, we show that aging HCs display both universal molecular blueprints, such as genomic instability, mitochondrial dysfunction, and impaired proteostasis, and cell type-specific aging signatures associated with deterioration of hair bundles and mechanotransduction. These signatures are also observed in aged human vestibular HCs, suggesting shared mechanisms. Importantly, morphological and functional analysis revealed that bundle degeneration and vestibular functional decline precede HC loss, highlighting the deterioration of mechanotransduction as a key contributor to ARVD. Furthermore, molecular and cellular changes associated with aging signatures are less pronounced in vestibular HCs than in cochlear HCs, underscoring the different pace of aging between the two mammalian inner ear sensory epithelia.

Characterization of a primary cellular airway model for inhalative drug delivery in comparison with the established permanent cell lines CaLu3 and RPMI 2650

Purpose

For optimization of respiratory drug delivery, the selection of suitable in vitro cell models plays an important role in predicting the efficacy and safety of (bio)pharmaceutics and pharmaceutical formulations. Therefore, an in-depth comparison of different primary and permanent in vitro cellular airway models was performed with a focus on selecting a suitable model for inhalative antibodies.

Methods

Primary cells isolated from the porcine trachea were compared with the established human cell lines CaLu3 and RPMI 2650. The in vitro models were characterized for different epithelial markers by real-time quantitative polymerase chain reaction, which provides insight into the cellular composition of each model. For a few selected markers, the results from RT-qPCR were confirmed via immunofluorescence. Barrier integrity was assessed by transepithelial electrical resistance measurements and FITC-dextran permeability.

Results

Primary cell models retain key features of the respiratory epithelium, e.g., the formation of a tight epithelial barrier, mucin production, and the presence of club/basal cells. Furthermore, the expression of Fc receptors in the primary cell models closely resembles that in respiratory mucosal tissue, an essential parameter to consider when developing therapeutic antibodies for inhalation.

Conclusion

The study underlines the importance of selecting wisely appropriate in vitro models. Despite the greater effort and variability in cultivating primary airway cells, they are far superior to permanent cells and a suitable model for drug development.

Supplementary Information

The online version contains supplementary material available at 10.1007/s44164-024-00079-y.

Primary Cilia as a Tumor Marker in Pituitary Neuroendocrine Tumors

Pituitary neuroendocrine tumors (PitNETs) account for approximately 15% of all intracranial neoplasms. Although they usually appear to be benign, some tumors display worse behavior, displaying rapid growth, invasion, refractoriness to treatment, and recurrence. Increasing evidence supports the role of primary cilia (PC) in regulating cancer development. Here, we showed that PC are significantly increased in PitNETs and are associated with increased tumor invasion and recurrence. Serial electron micrographs of PITNETs demonstrated different ciliation phenotypes (dot-like versus normal-like cilia) that represented PC at different stages of ciliogenesis. Molecular findings demonstrated that 123 ciliary-associated genes (eg, doublecortin domain containing protein 2, Sintaxin-3, and centriolar coiled-coil protein 110) were dysregulated in PitNETs, representing the upregulation of markers at different stages of intracellular ciliogenesis. Our results demonstrate, for the first time, that ciliogenesis is increased in PitNETs, suggesting that this process might be used as a potential target for therapy in the future.

Sperm DNA methylation defects in a new mouse model of the 5,10-methylenetetrahydrofolate reductase 677C>T variant and correction with moderate dose folic acid supplementation

5,10-Methylenetetrahydrofolate reductase (MTHFR) is an enzyme that plays a key role in providing methyl groups for DNA methylation, including during spermatogenesis. A common genetic variant in humans (MTHFR 677C>T) results in reduced enzyme activity and has been linked to various disorders, including male infertility. A new animal model has been created by reproducing the human equivalent of the polymorphism in mice using CRISPR/Cas9. Biochemical parameters in the Mthfr 677TT mice recapitulate alterations found in MTHFR 677TT men. Our aims were to characterize the sperm DNA methylome of the Mthfr 677CC and TT mice on a control diet (2 mg folic acid/kg diet) and assess the effects of folic acid supplementation (10 mg/kg diet) on the sperm DNA methylome. Body and reproductive organ weights, testicular sperm counts, and histology were examined. DNA methylation in sperm was assessed using bisulfite pyrosequencing and whole-genome bisulfite sequencing (WGBS). Reproductive parameters and locus-specific imprinted gene methylation were unaffected by genotype or diet. Using WGBS, sperm from 677TT mice had 360 differentially methylated tiles as compared to 677CC mice, predominantly hypomethylation (60% of tiles). Folic acid supplementation mostly caused hypermethylation in sperm of males of both genotypes and was found to partially correct the DNA methylation alterations in sperm associated with the TT genotype. The new mouse model will be useful in understanding the role of MTHFR deficiency in male fertility and in designing folate supplementation regimens for the clinic.

Single-cell analysis reveals insights into epithelial abnormalities in ovarian endometriosis

Ovarian endometriosis is characterized by the growth of endometrial tissue within the ovary, causing infertility and chronic pain. However, its pathophysiology remains unclear. Utilizing high-precision single-cell RNA sequencing, we profile the normal, eutopic, and ectopic endometrium from 34 individuals across proliferative and secretory phases. We observe an increased proportion of ciliated cells in both eutopic and ectopic endometrium, characterized by a diminished expression of estrogen sulfotransferase, which likely confers apoptosis resistance. After translocating to ectopic lesions, endometrial epithelium upregulates nicotinamide N-methyltransferase expression that inhibits apoptosis by promoting deacetylation and subsequent nuclear exclusion of transcription factor forkhead box protein O1, thereby leading to the downregulation of the apoptotic gene BIM. Moreover, epithelial cells in ectopic lesions elevate HLA class II complex expression, which stimulates CD4+ T cells and consequently contributes to chronic inflammation. Altogether, our study provides a comprehensive atlas of ovarian endometriosis and highlights potential therapeutic targets for modulating apoptosis and inflammation.

Convergent differentiation of multiciliated cells

Multiciliated cells (MCCs) are epithelial cells that control body fluid flow and contribute to the clearance of pathogenic microbes and other particles from the airways, egg transport in oviducts, and circulation of cerebrospinal fluid in the central nervous system. Although MCCs have shared functions to control fluid flow via coordinated motility of multiple ciliary structures, they are found in multiple mammalian tissues originating from distinct germ layers and differentiate via distinct developmental pathways. To understand the similarities and differences of MCCs in multiple tissues, we investigated single-cell transcriptome data of nasal epithelial cells, bronchial tubes, fallopian tubes, and ependymal cells in the subventricular zone from humans and mice by cross-species data integration. Expression of cilia-associated genes was indistinguishable between these MCCs, although cell populations had unique properties by the species and tissue, demonstrating that they share the same final differentiation status for ciliary functions. We further analyzed the final differentiation step of MCCs from their distinctive progenitors and confirmed their convergent gene set expression for ciliogenesis at the final step. These results may provide new insight into understanding ciliogenesis during the developmental process.

Single-cell analysis of shared signatures and transcriptional diversity during zebrafish development

During development, animals generate distinct cell populations with specific identities, functions, and morphologies. We mapped transcriptionally distinct populations across 489,686 cells from 62 stages during wild-type zebrafish embryogenesis and early larval development (3-120 h post-fertilization). Using these data, we identified the limited catalog of gene expression programs reused across multiple tissues and their cell-type-specific adaptations. We also determined the duration each transcriptional state is present during development and identify unexpected long-term cycling populations. Focused clustering and transcriptional trajectory analyses of non-skeletal muscle and endoderm identified transcriptional profiles and candidate transcriptional regulators of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and recently discovered best4+ cells. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

Cellular heterogeneity of the developing worker honey bee (Apis mellifera) pupa: a single cell transcriptomics analysis

It is estimated that animals pollinate 87.5% of flowering plants worldwide and that managed honey bees (Apis mellifera) account for 30-50% of this ecosystem service to agriculture. In addition to their important role as pollinators, honey bees are well-established insect models for studying learning and memory, behavior, caste differentiation, epigenetic mechanisms, olfactory biology, sex determination, and eusociality. Despite their importance to agriculture, knowledge of honey bee biology lags behind many other livestock species. In this study, we have used scRNA-Seq to map cell types to different developmental stages of the worker honey bee (prepupa at day 11 and pupa at day 15) and sought to determine their gene expression signatures. To identify cell-type populations, we examined the cell-to-cell network based on the similarity of the single-cells transcriptomic profiles. Grouping similar cells together we identified 63 different cell clusters of which 17 clusters were identifiable at both stages. To determine genes associated with specific cell populations or with a particular biological process involved in honey bee development, we used gene coexpression analysis. We combined this analysis with literature mining, the honey bee protein atlas, and gene ontology analysis to determine cell cluster identity. Of the cell clusters identified, 17 were related to the nervous system and sensory organs, 7 to the fat body, 19 to the cuticle, 5 to muscle, 4 to compound eye, 2 to midgut, 2 to hemocytes, and 1 to malpighian tubule/pericardial nephrocyte. To our knowledge, this is the first whole single-cell atlas of honey bees at any stage of development and demonstrates the potential for further work to investigate their biology at the cellular level.

VWA3A-derived ependyma promoter drives increased therapeutic protein secretion into the CSF

Recombinant adeno-associated viral vectors (rAAVs) are a promising strategy to treat neurodegenerative diseases because of their ability to infect non-dividing cells and confer long-term transgene expression. Despite an ever-growing library of capsid variants, widespread delivery of AAVs in the adult central nervous system remains a challenge. We have previously demonstrated successful distribution of secreted proteins by infection of the ependyma, a layer of post-mitotic epithelial cells lining the ventricles of the brain and central column of the spinal cord, and subsequent protein delivery via the cerebrospinal fluid (CSF). Here we define a functional ependyma promoter to enhance expression from this cell type. Using RNA sequencing on human autopsy samples, we identified disease- and age-independent ependyma gene signatures. Associated promoters were cloned and screened as libraries in mouse and rhesus macaque to reveal cross-species function of a human DNA-derived von Willebrand factor domain containing 3A (VWA3A) promoter. When tested in mice, our VWA3A promoter drove strong, ependyma-localized expression of eGFP and increased secreted ApoE protein levels in the CSF by 2-12× over the ubiquitous iCAG promoter.

Concurrent loss of ciliary genes WDR93 and CFAP46 in phylogenetically distant birds

The respiratory system is the primary route of infection for many contagious pathogens. Mucociliary clearance of inhaled pathogens is an important innate defence mechanism sustained by the rhythmic movement of epithelial cilia. To counter host defences, viral pathogens target epithelial cells and cilia. For instance, the avian influenza virus that targets ciliated cells modulates the expression of WDR93, a central ciliary apparatus C1d projection component. Lineage-specific prevalence of such host defence genes results in differential susceptibility. In this study, the comparative analysis of approximately 500 vertebrate genomes from seven taxonomic classes spanning 73 orders confirms the widespread conservation of WDR93 across these different vertebrate groups. However, we established loss of the WDR93 in landfowl, geese and other phylogenetically independent bird species due to gene-disrupting changes. The lack of WDR93 transcripts in species with gene loss in contrast to its expression in species with an intact gene confirms gene loss. Notably, species with WDR93 loss have concurrently lost another C1d component, CFAP46, through large segmental deletions. Understanding the consequences of such gene loss may provide insight into their role in host-pathogen interactions and benefit global pathogen surveillance efforts by prioritizing species missing host defence genes and identifying putative zoonotic reservoirs.

Congenital heart defects caused by FOXJ1

FOXJ1 is expressed in ciliated cells of the airways, testis, oviduct, central nervous system and the embryonic left-right organizer. Ablation or targeted mutation of Foxj1 in mice, zebrafish and frogs results in loss of ciliary motility and/or reduced length and number of motile cilia, affecting the establishment of the left-right axis. In humans, heterozygous pathogenic variants in FOXJ1 cause ciliopathy leading to situs inversus, obstructive hydrocephalus and chronic airway disease. Here, we report a novel truncating FOXJ1 variant (c.784_799dup; p.Glu267Glyfs*12) identified by clinical exome sequencing from a patient with isolated congenital heart defects (CHD) which included atrial and ventricular septal defects, double outlet right ventricle (DORV) and transposition of the great arteries. Functional experiments show that FOXJ1 c.784_799dup; p.Glu267Glyfs*12, unlike FOXJ1, fails to induce ectopic cilia in frog epidermis in vivo or to activate the ADGB promoter, a downstream target of FOXJ1 in cilia, in transactivation assays in vitro. Variant analysis of patients with heterotaxy or heterotaxy-related CHD indicates that pathogenic variants in FOXJ1 are an infrequent cause of heterotaxy. Finally, we characterize embryonic-stage CHD in Foxj1 loss-of-function mice, demonstrating randomized heart looping. Abnormal heart looping includes reversed looping (dextrocardia), ventral looping and no looping/single ventricle hearts. Complex CHDs revealed by histological analysis include atrioventricular septal defects, DORV, single ventricle defects as well as abnormal position of the great arteries. These results indicate that pathogenic variants in FOXJ1 can cause isolated CHD.

LRRC6 regulates biogenesis of motile cilia by aiding FOXJ1 translocation into the nucleus

BACKGROUND

LRRC6 is an assembly factor for dynein arms in the cytoplasm of motile ciliated cells, and when mutated, dynein arm components remained in the cytoplasm. Here, we demonstrate the role of LRRC6 in the active nuclear translocation of FOXJ1, a master regulator for cilia-associated gene transcription.

METHODS

We generated Lrrc6 knockout (KO) mice, and we investigated the role of LRRC6 on ciliopathy development by using proteomic, transcriptomic, and immunofluorescence analysis. Experiments on mouse basal cell organoids confirmed the biological relevance of our findings.

RESULTS

The absence of LRRC6 in multi-ciliated cells hinders the assembly of ODA and IDA components of cilia; in this study, we showed that the overall expression of proteins related to cilia decreased as well. Expression of cilia-related transcripts, specifically ODA and IDA components, dynein axonemal assembly factors, radial spokes, and central apparatus was lower in Lrrc6 KO mice than in wild-type mice. We demonstrated that FOXJ1 was present in the cytoplasm and translocated into the nucleus when LRRC6 was expressed and that this process was blocked by INI-43, an importin α inhibitor.

CONCLUSIONS

Taken together, these results hinted at the LRRC6 transcriptional regulation of cilia-related genes via the nuclear translocation of FOXJ1. Video Abstract.

Identification of a heterogeneous and dynamic ciliome during embryonic development and cell differentiation

Primary cilia are nearly ubiquitous organelles that transduce molecular and mechanical signals. Although the basic structure of the cilium and the cadre of genes that contribute to ciliary formation and function (the ciliome) are believed to be evolutionarily conserved, the presentation of ciliopathies with narrow, tissue-specific phenotypes and distinct molecular readouts suggests that an unappreciated heterogeneity exists within this organelle. Here, we provide a searchable transcriptomic resource for a curated primary ciliome, detailing various subgroups of differentially expressed genes within the ciliome that display tissue and temporal specificity. Genes within the differentially expressed ciliome exhibited a lower level of functional constraint across species, suggesting organism and cell-specific function adaptation. The biological relevance of ciliary heterogeneity was functionally validated by using Cas9 gene-editing to disrupt ciliary genes that displayed dynamic gene expression profiles during osteogenic differentiation of multipotent neural crest cells. Collectively, this novel primary cilia-focused resource will allow researchers to explore longstanding questions related to how tissue and cell-type specific functions and ciliary heterogeneity may contribute to the range of phenotypes associated with ciliopathies.

Ckb and Ybx2 interact with Ribc2 and are necessary for the ciliary beating of multi-cilia

BACKGROUND

Motile cilia in a vertebrate are important to sustaining activities of life. Fluid flow on the apical surface of several tissues, including bronchial epithelium, ependymal epithelium, and fallopian tubules is generated by the ciliary beating of motile cilia. Multi-ciliated cells in ependymal tissue are responsible for the circulation of cerebrospinal fluid (CSF), which is essential for the development and homeostasis of the central nervous system, and airway tissues are protected from external contaminants by cilia-driven mucosal flow over the top of the airway epithelium.

OBJECTIVE

A previous study reported that reduction of Ribc2 protein leads to disruption of ciliary beating in multi-ciliated cells. However, knowledge regarding the molecular function of Ribc2 is limited, thus currently available information is also limited. Therefore, we evaluated the importance of proteins involved in the interaction with Ribc2 in the process of ciliary beating.

METHODS

Immunoprecipitation and mass spectrometry analysis was performed for the discovery of proteins involved in the interaction with Ribc2. Expression of the target gene was inhibited by injection of antisense morpholinos and measurement of the fluid flow on the embryonic epidermis of Xenopus was performed using fluorescent beads for examination of the ciliary beating of multi cilia. In addition, the flag-tagged protein was expressed by injection of mRNA and the changes in protein localization in the cilia were measured by immunostaining and western blot analysis for analysis of the molecular interaction between Ribc2 and Ribc2 binding proteins in multi-cilia.

RESULTS

The IP/MS analysis identified Ckb and Ybx2 as Ribc2 binding proteins and our results showed that localization of both Ckb and Ybx2 occurs at the axoneme of multi-cilia on the embryonic epithelium of Xenopus laevis. In addition, our findings confirmed that knock-down of Ckb or Ybx2 resulted in abnormal ciliary beating and reduction of cilia-driven fluid flow on multi-cilia of Xenopus laevis. In addition, significantly decreased localization of Ckb or Ybx2 in the ciliary axoneme was observed in Ribc2-depleted multi-cilia.

CONCLUSION

Ckb and Ybx2 are involved in the interaction with Ribc2 and are necessary for the ciliary beating of multi-cilia.

Variable phenotypes and penetrance between and within different zebrafish ciliary transition zone mutants

Meckel syndrome, nephronophthisis, Joubert syndrome and Bardet-Biedl syndrome are caused by mutations in proteins that localize to the ciliary transition zone (TZ). The phenotypically distinct syndromes suggest that these TZ proteins have differing functions. However, mutations in a single TZ gene can result in multiple syndromes, suggesting that the phenotype is influenced by modifier genes. We performed a comprehensive analysis of ten zebrafish TZ mutants, including mks1, tmem216, tmem67, rpgrip1l, cc2d2a, b9d2, cep290, tctn1, nphp1 and nphp4, as well as mutants in ift88 and ift172. Our data indicate that variations in phenotypes exist between different TZ mutants, supporting different tissue-specific functions of these TZ genes. Further, we observed phenotypic variations within progeny of a single TZ mutant, reminiscent of multiple disease syndromes being associated with mutations in one gene. In some mutants, the dynamics of the phenotype became complex with transitory phenotypes that are corrected over time. We also demonstrated that multiple-guide-derived CRISPR/Cas9 F0 'crispant' embryos recapitulate zygotic null phenotypes, and rapidly identified ciliary phenotypes in 11 cilia-associated gene candidates (ankfn1, ccdc65, cfap57, fhad1, nme7, pacrg, saxo2, c1orf194, ttc26, zmynd12 and cfap52).

Is IIIG9 a New Protein with Exclusive Ciliary Function? Analysis of Its Potential Role in Cancer and Other Pathologies

The identification of new proteins that regulate the function of one of the main cellular phosphatases, protein phosphatase 1 (PP1), is essential to find possible pharmacological targets to alter phosphatase function in various cellular processes, including the initiation and development of multiple diseases. IIIG9 is a regulatory subunit of PP1 initially identified in highly polarized ciliated cells. In addition to its ciliary location in ependymal cells, we recently showed that IIIG9 has extraciliary functions that regulate the integrity of adherens junctions. In this review, we perform a detailed analysis of the expression, localization, and function of IIIG9 in adult and developing normal brains. In addition, we provide a 3D model of IIIG9 protein structure for the first time, verifying that the classic structural and conformational characteristics of the PP1 regulatory subunits are maintained. Our review is especially focused on finding evidence linking IIIG9 dysfunction with the course of some pathologies, such as ciliopathies, drug dependence, diseases based on neurological development, and the development of specific high-malignancy and -frequency brain tumors in the pediatric population. Finally, we propose that IIIG9 is a relevant regulator of PP1 function in physiological and pathological processes in the CNS.

Clustering of Genetic Anomalies of Cilia Outer Dynein Arm and Central Apparatus in Patients with Transposition of the Great Arteries

Transposition of the great arteries (TGA) is a congenital heart defect with a complex pathogenesis that has not been fully elucidated. In this study, we performed whole-exome sequencing (WES) in isolated TGA-diagnosed patients and analyzed genes of motile and non-motile cilia ciliogenesis and ciliary trafficking, as well as genes previously associated with this heart malformation. Deleterious missense and splicing variants of genes DNAH9, DNAH11, and ODAD4 of cilia outer dynein arm and central apparatus, HYDIN, were found in our TGA patients. Remarkable, there is a clustering of deleterious genetic variants in cilia genes, suggesting it could be an oligogenic disease. Our data evidence the genetic diversity and etiological complexity of TGA and point out that population allele determination and genetic aggregation studies are required to improve genetic counseling.

Pathogenic Variants in Cardiomyopathy Disorder Genes Underlie Pediatric Myocarditis—Further Impact of Heterozygous Immune Disorder Gene Variants?

Myocarditis is an inflammatory disease of the heart. Pediatric myocarditis with the dilated cardiomyopathy (DCM) phenotype may be caused by likely pathogenic or pathogenic genetic variants [(L)P] in cardiomyopathy (CMP) genes. Systematic analysis of immune disorder gene defects has not been performed so far. We analyzed 12 patients with biopsy-proven myocarditis and the DCM phenotype together with their parents using whole-exome sequencing (WES). The WES data were filtered for rare pathogenic variants in CMP (n = 89) and immune disorder genes (n = 631). Twelve children with a median age of 2.9 (1.0–6.8) years had a mean left ventricular ejection fraction of 28% (22–32%) and myocarditis was confirmed by endomyocardial biopsy. Patients with primary immunodeficiency were excluded from the study. Four patients underwent implantation of a ventricular assist device and subsequent heart transplantation. Genetic analysis of the 12 families revealed an (L)P variant in the CMP gene in 8/12 index patients explaining DCM. Screening of recessive immune disorder genes identified a heterozygous (L)P variant in 3/12 index patients. This study supports the genetic impact of CMP genes for pediatric myocarditis with the DCM phenotype. Piloting the idea that additional immune-related genetic defects promote myocarditis suggests that the presence of heterozygous variants in these genes needs further investigation. Altered cilium function might play an additional role in inducing inflammation in the context of CMP.

Graphia: A platform for the graph-based visualisation and analysis of high dimensional data

Graphia is an open-source platform created for the graph-based analysis of the huge amounts of quantitative and qualitative data currently being generated from the study of genomes, genes, proteins metabolites and cells. Core to Graphia’s functionality is support for the calculation of correlation matrices from any tabular matrix of continuous or discrete values, whereupon the software is designed to rapidly visualise the often very large graphs that result in 2D or 3D space. Following graph construction, an extensive range of measurement algorithms, routines for graph transformation, and options for the visualisation of node and edge attributes are available, for graph exploration and analysis. Combined, these provide a powerful solution for the interpretation of high-dimensional data from many sources, or data already in the form of a network or equivalent adjacency matrix. Several use cases of Graphia are described, to showcase its wide range of applications in the analysis biological data. Graphia runs on all major desktop operating systems, is extensible through the deployment of plugins and is freely available to download from https://graphia.app/.

Graphia is a new visual analytics platform specifically created for the network-based analysis of large and complex data, such as that generated in huge amounts by modern biological analyses. It works in a data agnostic, hypothesis-free manner to generate correlation networks from any table of numerical or discrete values, thereafter providing a means to rapidly visualise the often very large networks that result, in either 2D or 3D space. Following network construction, the tool offers an extensive range of analysis algorithms, routines for network transformation, and options for the visualisation of metadata. This provides a powerful analysis solution for the exploration and interpretation of high-dimensional data from any source, as well as any data already defined as a network. Several use cases of Graphia are described to showcase its wide range of applications in the analysis biological data. Graphia is open source and free to all.

Single-Cell RNA-Seq Analysis Reveals Lung Epithelial Cell Type-Specific Responses to HDM and Regulation by Tet1

Tet1 protects against house dust mite (HDM)-induced lung inflammation in mice and alters the lung methylome and transcriptome. In order to explore the role of Tet1 in individual lung epithelial cell types in HDM-induced inflammation, we established a model of HDM-induced lung inflammation in Tet1 knockout and littermate wild-type mice, then studied EpCAM+ lung epithelial cells using single-cell RNA-seq analysis. We identified eight EpCAM+ lung epithelial cell types, among which AT2 cells were the most abundant. HDM challenge altered the relative abundance of epithelial cell types and resulted in cell type-specific transcriptomic changes. Bulk and cell type-specific analysis also showed that loss of Tet1 led to the altered expression of genes linked to augmented HDM-induced lung inflammation, including alarms, detoxification enzymes, oxidative stress response genes, and tissue repair genes. The transcriptomic regulation was accompanied by alterations in TF activities. Trajectory analysis supports that HDM may enhance the differentiation of AP and BAS cells into AT2 cells, independent of Tet1. Collectively, our data showed that lung epithelial cells had common and unique transcriptomic signatures of allergic lung inflammation. Tet1 deletion altered transcriptomic networks in various lung epithelial cells, which may promote allergen-induced lung inflammation.

Age-Related Differences in Structure and Function of Nasal Epithelial Cultures From Healthy Children and Elderly People

The nasal epithelium represents the first line of defense against inhaled pathogens, allergens, and irritants and plays a key role in the pathogenesis of a spectrum of acute and chronic airways diseases. Despite age-dependent clinical phenotypes triggered by these noxious stimuli, little is known about how aging affects the structure and function of the airway epithelium that is crucial for lung homeostasis and host defense. The aim of this study was therefore to determine age-related differences in structural and functional properties of primary nasal epithelial cultures from healthy children and non-smoking elderly people. To achieve this goal, highly differentiated nasal epithelial cultures were established from nasal brushes at air–liquid interface and used to study epithelial cell type composition, mucin (MUC5AC and MUC5B) expression, and ion transport properties. Furthermore, we determined age-dependent molecular signatures using global proteomic analysis. We found lower numeric densities of ciliated cells and higher levels of MUC5AC expression in cultures from children vs. elderly people. Bioelectric studies showed no differences in basal ion transport properties, ENaC-mediated sodium absorption, or CFTR-mediated chloride transport, but detected decreased calcium-activated TMEM16A-mediated chloride secretory responses in cultures from children vs. elderly people. Proteome analysis identified distinct age-dependent molecular signatures associated with ciliation and mucin biosynthesis, as well as other pathways implicated in aging. Our data identified intrinsic, age-related differences in structure and function of the nasal epithelium and provide a basis for further studies on the role of these findings in age-dependent airways disease phenotypes observed with a spectrum of respiratory infections and other noxious stimuli.

Sex-biased dynamics of three-spined stickleback (Gasterosteus aculeatus) gene expression patterns

The study of the differences between sexes presents an excellent model to unravel how phenotypic variation is achieved from a similar genetic background. Sticklebacks are of particular interest since evidence of a heteromorphic chromosome pair has not always been detected. The present study investigated sex-biased mRNA and small non-coding RNA (sncRNA) expression patterns in the brain, adipose tissues, and gonads of the three-spined stickleback. The sncRNA analysis indicated that regulatory functions occurred mainly in the gonads. Alleged miRNA-mRNA interactions were established and a mapping bias of differential expressed transcripts towards chromosome 19 was observed. Key players previously shown to control sex determination and differentiation in other fish species but also genes like gapdh were among the transcripts identified. This is the first report in the three-spined stickleback demonstrating tissue-specific expression comprising both mRNA and sncRNA between sexes, emphasizing the importance of mRNA-miRNA interactions as well as new presumed genes not yet identified to have gender-specific roles.

Cilium Expression Score Predicts Glioma Survival

The accurate classification, prognostication, and treatment of gliomas has been hindered by an existing cellular, genomic, and transcriptomic heterogeneity within individual tumors and their microenvironments. Traditional clustering is limited in its ability to distinguish heterogeneity in gliomas because the clusters are required to be exclusive and exhaustive. In contrast, biclustering can identify groups of co-regulated genes with respect to a subset of samples and vice versa. In this study, we analyzed 1,798 normal and tumor brain samples using an unsupervised biclustering approach. We identified co-regulated gene expression profiles that were linked to proximally located brain regions and detected upregulated genes in subsets of gliomas, associated with their histologic grade and clinical outcome. In particular, we present a cilium-associated signature that when upregulated in tumors is predictive of poor survival. We also introduce a risk score based on expression of 12 cilium-associated genes which is reproducibly informative of survival independent of other prognostic biomarkers. These results highlight the role of cilia in development and progression of gliomas and suggest potential therapeutic vulnerabilities for these highly aggressive tumors.

Kinase-anchoring proteins in ciliary signal transduction

Historically, the diffusion of chemical signals through the cell was thought to occur within a cytoplasmic soup bounded by the plasma membrane. This theory was predicated on the notion that all regulatory enzymes are soluble and moved with a Brownian motion. Although enzyme compartmentalization was initially rebuffed by biochemists as a 'last refuge of a scoundrel', signal relay through macromolecular complexes is now accepted as a fundamental tenet of the burgeoning field of spatial biology. A-Kinase anchoring proteins (AKAPs) are prototypic enzyme-organizing elements that position clusters of regulatory proteins at defined subcellular locations. In parallel, the primary cilium has gained recognition as a subcellular mechanosensory organelle that amplifies second messenger signals pertaining to metazoan development. This article highlights advances in our understanding of AKAP signaling within the primary cilium and how defective ciliary function contributes to an increasing number of diseases known as ciliopathies.

Spatiotemporal transcriptional dynamics of the cycling mouse oviduct

Female fertility in mammals requires iterative remodeling of the entire adult female reproductive tract across the menstrual/estrous cycle. However, while transcriptome dynamics across the estrous cycle have been reported in human and bovine models, no global analysis of gene expression across the estrous cycle has yet been reported for the mouse. Here, we examined the cellular composition and global transcriptional dynamics of the mouse oviduct along the anteroposterior axis and across the estrous cycle. We observed robust patterns of differential gene expression along the anteroposterior axis, but we found surprisingly few changes in gene expression across the estrous cycle. Notable gene expression differences along the anteroposterior axis included a surprising enrichment for genes related to embryonic development, such as Hox and Wnt genes. The relatively stable transcriptional dynamics across the estrous cycle differ markedly from other mammals, leading us to speculate that this is an evolutionarily derived state that may reflect the extremely rapid five-day mouse estrous cycle. This dataset fills a critical gap by providing an important genomic resource for a highly tractable genetic model of mammalian female reproduction.

Autosomal recessive loci contribute significantly to quantitative variation of male fertility in a dairy cattle population

BACKGROUND

Cattle are ideally suited to investigate the genetics of male fertility. Semen from individual bulls is used for thousands of artificial inseminations for which the fertilization success is monitored. Results from the breeding soundness examination and repeated observations of semen quality complement the fertility evaluation for each bull.

RESULTS

In a cohort of 3881 Brown Swiss bulls that had genotypes at 683,609 SNPs, we reveal four novel recessive QTL for male fertility on BTA1, 18, 25, and 26 using haplotype-based association testing. A QTL for bull fertility on BTA1 is also associated with sperm head shape anomalies. All other QTL are not associated with any of the semen quality traits investigated. We perform complementary fine-mapping approaches using publicly available transcriptomes as well as whole-genome sequencing data of 125 Brown Swiss bulls to reveal candidate causal variants. We show that missense or nonsense variants in SPATA16, VWA3A, ENSBTAG00000006717 and ENSBTAG00000019919 are in linkage disequilibrium with the QTL. Using whole-genome sequence data, we detect strong association (P = 4.83 × 10^- 12) of a missense variant (p.Ile193Met) in SPATA16 with male fertility. However, non-coding variants exhibit stronger association at all QTL suggesting that variants in regulatory regions contribute to variation in bull fertility.

CONCLUSION

Our findings in a dairy cattle population provide evidence that recessive variants may contribute substantially to quantitative variation in male fertility in mammals. Detecting causal variants that underpin variation in male fertility remains difficult because the most strongly associated variants reside in poorly annotated non-coding regions.

Limitations and opportunities in the pharmacotherapy of ciliopathies

Ciliopathies are a family of rather diverse conditions, which have been grouped based on the finding of altered or dysfunctional cilia, potentially motile, small cellular antennae extending from the surface of postmitotic cells. Cilia-related disorders include embryonically arising conditions such as Joubert, Usher or Kartagener syndrome, but also afflictions with a postnatal or even adult onset phenotype, i.e. autosomal dominant polycystic kidney disease. The majority of ciliopathies are syndromic rather than affecting only a single organ due to cilia being found on almost any cell in the human body. Overall ciliopathies are considered rare diseases. Despite that, pharmacological research and the strive to help these patients has led to enormous therapeutic advances in the last decade. In this review we discuss new treatment options for certain ciliopathies, give an outlook on promising future therapeutic strategies, but also highlight the limitations in the development of therapeutic approaches of ciliopathies.

CiliOPD: a ciliopathy-associated COPD endotype

The pathophysiology of chronic obstructive pulmonary disease (COPD) relies on airway remodelling and inflammation. Alterations of mucociliary clearance are a major hallmark of COPD caused by structural and functional cilia abnormalities. Using transcriptomic databases of whole lung tissues and isolated small airway epithelial cells (SAEC), we comparatively analysed cilia-associated and ciliopathy-associated gene signatures from a set of 495 genes in 7 datasets including 538 non-COPD and 508 COPD patients. This bio-informatics approach unveils yet undescribed cilia and ciliopathy genes associated with COPD including NEK6 and PROM2 that may contribute to the pathology, and suggests a COPD endotype exhibiting ciliopathy features (CiliOPD).

Motile cilia and airway disease

A finely regulated system of airway epithelial development governs the differentiation of motile ciliated cells of the human respiratory tract, conferring the body's mucociliary clearance defence system. Human cilia dysfunction can arise through genetic mutations and this is a cause of debilitating disease morbidities that confer a greatly reduced quality of life. The inherited human motile ciliopathy disorder, primary ciliary dyskinesia (PCD), can arise from mutations in genes affecting various aspects of motile cilia structure and function through deficient production, transport and assembly of cilia motility components or through defective multiciliogenesis. Our understanding about the development of the respiratory epithelium, motile cilia biology and the implications for human pathology has expanded greatly over the past 20 years since isolation of the first PCD gene, rising to now nearly 50 genes. Systems level insights about cilia motility in health and disease have been made possible through intensive molecular and omics (genomics, transcriptomics, proteomics) research, applied in ciliate organisms and in animal and human disease modelling. Here, we review ciliated airway development and the genetic stratification that underlies PCD, for which the underlying genotype can increasingly be connected to biological mechanism and disease prognostics. Progress in this field can facilitate clinical translation of research advances, with potential for great medical impact, e.g. through improvements in ciliopathy disease diagnosis, management, family counselling and by enhancing the potential for future genetically tailored approaches to disease therapeutics.