NME5 frameshift variant in Alaskan Malamutes with primary ciliary dyskinesia

Recognition and Diagnosis of Underlying Disease Processes in Bacterial Pneumonia

Bacterial pneumonia can present in both dogs and cats as either acute or chronic disease, and clinical signs may suggest respiratory or systemic impairment. Aspiration pneumonia, a common form of disease, can be caused by esophageal or laryngeal disease, vomiting, or altered mentation, such as with CNS disease or general anesthesia. Bacterial pneumonia can also develop due to inhaled airway foreign bodies, reduced local immune function, or impaired clearance of respiratory secretions. Rapid identification of the etiology and associated risk factors for pneumonia can help veterinarians develop a targeted therapeutic plan and prevent recurrence of disease.

Mapping the Most Common Founder Variant in RSPH9 That Causes Primary Ciliary Dyskinesia in Multiple Consanguineous Families of Bedouin Arabs

INTRODUCTION

Primary ciliary dyskinesia (PCD) is a congenital thoracic disorder caused by dysfunction of motile cilia, resulting in insufficient mucociliary clearance of the lungs. The overall aim of this study is to identify causative defective genes in PCD-affected individuals in the Kuwaiti population.

METHODS

A cohort of multiple consanguineous PCD families was identified from Kuwaiti patients and genomic DNA from the family members was isolated using standard procedures. The DNA samples from all affected individuals were analyzed by whole exome sequencing (WES). Transmission electron microscopy (TEM) and immunofluorescent analysis (IF) were performed on samples obtained by nasal brushings to identify specific structural abnormalities within ciliated cells.

RESULTS

Here, we present six multiplex families with 11 patients who all presented with typical PCD symptoms. Ten out of eleven patients inherited a 3 bp homozygous deletion of GAA in RSPH9, whereas the eleventh patients inherited this variant in trans with a frameshift deletion in RSPH9. Genetic results were confirmed by segregation analysis. The in-frame deletion of GAA in RSPH9 has previously been published as pathogenic in both annotated RSPH9 transcript variants (1 and 2). In contrast, the previously unpublished RSPH9 frameshift deletion identified in KU-15.IV2 impacts only RSPH9 transcript variant two. Regarding all 11 PCD individuals analyzed, IF results demonstrated absence of RSPH9 protein and TEM analysis showed the typical findings in RSPH9 mutant individuals.

CONCLUSIONS

We present the largest cohort of PCD individuals affected by the founder in-frame deletion GAA in RSPH9. This founder variant is the most common PCD-causing variant in Bedouin Arabs in Kuwait.

Ordered deployment of distinct ciliary beating machines in growing axonemes of vertebrate multiciliated cells

The beating of motile cilia requires the coordinated action of diverse machineries that include not only the axonemal dynein arms, but also the central apparatus, the radial spokes, and the microtubule inner proteins. These machines exhibit complex radial and proximodistal patterns in mature axonemes, but little is known about the interplay between them during motile ciliogenesis. Here, we describe and quantify the relative rates of axonemal deployment for these diverse cilia beating machineries during the final stages of differentiation of Xenopus epidermal multiciliated cells.

Schmidtea mediterranea as a Model Organism to Study the Molecular Background of Human Motile Ciliopathies

Cilia and flagella are evolutionarily conserved organelles that form protrusions on the surface of many growth-arrested or differentiated eukaryotic cells. Due to the structural and functional differences, cilia can be roughly classified as motile and non-motile (primary). Genetically determined dysfunction of motile cilia is the basis of primary ciliary dyskinesia (PCD), a heterogeneous ciliopathy affecting respiratory airways, fertility, and laterality. In the face of the still incomplete knowledge of PCD genetics and phenotype-genotype relations in PCD and the spectrum of PCD-like diseases, a continuous search for new causative genes is required. The use of model organisms has been a great part of the advances in understanding molecular mechanisms and the genetic basis of human diseases; the PCD spectrum is not different in this respect. The planarian model (Schmidtea mediterranea) has been intensely used to study regeneration processes, and-in the context of cilia-their evolution, assembly, and role in cell signaling. However, relatively little attention has been paid to the use of this simple and accessible model for studying the genetics of PCD and related diseases. The recent rapid development of the available planarian databases with detailed genomic and functional annotations prompted us to review the potential of the S. mediterranea model for studying human motile ciliopathies.

STK36 splice site variant in an Australian Shepherd dog with primary ciliary dyskinesia

Primary ciliary dyskinesia (PCD) represents a group of diseases characterized by impaired movement of cilia and subsequent health problems in diverse organ systems, notably the respiratory tract. Almost 50 candidate genes for PCD are known in humans. In this study, we investigated an Australian Shepherd dog with a history of recurrent respiratory infections and nasal discharge. A transmission electron microscopy investigation led to the diagnosis of PCD with central pair defect, in which the normal 9:2 arrangement of respiratory cilia was altered and reduced to a 9:0 arrangement. Whole genome sequencing data from the affected dog was obtained and searched for variants in PCD candidate genes that were not present in 918 control genomes from different breeds. This revealed a homozygous single base pair exchange at a splice site of STK36, XM_038585732.1:c.2868-1G>A. The mutant allele was absent from 281 additionally genotyped Australian Shepherd dogs. RT-PCR confirmed aberrant splicing in the affected dog with the skipping of exon 20 and the insertion of a cryptic exon, which is predicted to lead to a premature stop codon and truncation of 36% of the STK36 wild-type open reading frame, XP_038441660.1:(p.Met957Profs*11). STK36 variants were previously reported to cause PCD in humans and mice. The knowledge from other species together with the absence of the mutant allele in more than 1000 control dogs suggests STK36:c.2868-1G>A as the most likely candidate variant for PCD in the investigated case.

Bioinformatic analysis and machine learning to identify the diagnostic biomarkers and immune infiltration in adenomyosis

Background: Adenomyosis is a hormone-dependent benign gynecological disease characterized by the invasion of the endometrium into the myometrium. Women with adenomyosis can suffer from abnormal uterine bleeding, severe pelvic pain, and subfertility or infertility, which can interfere with their quality of life. However, effective diagnostic biomarkers for adenomyosis are currently lacking. The aim of this study is to explore the mechanism of adenomyosis by identifying biomarkers and potential therapeutic targets for adenomyosis and analyzing their correlation with immune infiltration in adenomyosis. Methods: Two datasets, GSE78851 and GSE68870, were downloaded and merged for differential expression analysis and functional enrichment analysis using R software. Weighted gene co-expression network analysis (WGCNA), the least absolute shrinkage and selection operator (LASSO), and support vector machine-recursive feature elimination (SVE-RFE) were combined to explore candidate genes. Quantitative reverse transcriptase PCR (qRT-PCR) was conducted to verify the biomarkers and receiver operating characteristic curve analysis was used to assess the diagnostic value of each biomarker. Single-sample Gene Set Enrichment Analysis (ssGSEA) and CIBERSORT were used to explore immune cell infiltration in adenomyosis and the correlation between diagnostic biomarkers and immune cells. Results: A total of 318 genes were differentially expressed. Through the analysis of differentially expressed genes and WGCNA, we obtained 189 adenomyosis-related genes. After utilizing the LASSO and SVM-RFE algorithms, four hub genes, namely, six-transmembrane epithelial antigen of the prostate-1 (STEAP1), translocase of outer mitochondrial membrane 20 (TOMM20), glycosyltransferase eight domain-containing 2 (GLT8D2), and NME/NM23 family member 5 (NME5) expressed in nucleoside-diphosphate kinase, were identified and verified by qRT-PCR. Immune infiltration analysis indicated that T helper 17 cells, CD56dim natural killer cells, monocytes, and memory B-cell may be associated with the occurrence of adenomyosis. There were significant correlations between the diagnostic biomarkers and immune cells. Conclusion: STEAP1, TOMM20, GLT8D2, and NME5 were identified as potential biomarkers and therapeutic targets for adenomyosis. Immune infiltration may contribute to the onset and progression of adenomyosis.

Genomic inbreeding and runs of homozygosity analysis of indigenous cattle populations in southern China

Runs of homozygosity (ROH) are continuous homozygous segments from the common ancestor of parents. Evaluating ROH pattern can help to understand inbreeding level and genetic basis of important traits. In this study, three representative cattle populations including Leiqiong cattle (LQC), Lufeng cattle (LFC) and Hainan cattle (HNC) were genotyped using the Illumina BovineHD SNPs array (770K) to assess ROH pattern at genome wide level. Totally, we identified 26,537 ROH with an average of 153 ROH per individual. The sizes of ROH ranged from 0.5 to 53.26Mb, and the average length was 1.03Mb. The average of F_ROH ranged from 0.10 (LQC) to 0.15 (HNC). Moreover, we identified 34 ROH islands (with frequency > 0.5) across genome. Based on these regions, we observed several breed-specific candidate genes related to adaptive traits. Several common genes related to immunity (TMEM173, MZB1 and SIL1), and heat stress (DNAJC18) were identified in all three populations. Three genes related to immunity (UGP2), development (PURA) and reproduction (VPS54) were detected in both HNC and LQC. Notably, we identified several breed-specific genes related to sperm development (BRDT and SPAG6) and heat stress (TAF7) in HNC, and immunity (CDC23 and NME5) and development (WNT87) in LFC. Our findings provided valuable insights into understanding the genomic homozygosity pattern and promoting the conservation of genetic resources of Chinese indigenous cattle.

The Complex Functions of the NME Family-A Matter of Location and Molecular Activity

The family of NME proteins represents a quite complex group of multifunctional enzymes [...].

Genes causing congenital hydrocephalus: Their chromosomal characteristics of telomere proximity and DNA compositions

Congenital hydrocephalus (CH) is caused by genetic mutations, but whether factors impacting human genetic mutations are disease-specific remains elusive. Given two factors associated with high mutation rates, we reviewed how many disease-susceptible genes match with (i) proximity to telomeres or (ii) high adenine and thymine (A + T) content in human CH as compared to other disorders of the central nervous system (CNS). We extracted genomic information using a genome data viewer. Importantly, 98 of 108 genes causing CH satisfied (i) or (ii), resulting in >90% matching rate. However, such a high accordance no longer sustained as we checked two factors in Alzheimer's disease (AD) and/or familial Parkinson's disease (fPD), resulting in 84% and 59% matching, respectively. A disease-specific matching of telomere proximity or high A + T content predicts causative genes of CH much better than neurodegenerative diseases and other CNS conditions, likely due to sufficient number of known causative genes (n = 108) and precise determination and classification of the genotype and phenotype. Our analysis suggests a need for identifying genetic basis of both factors before human clinical studies, to prioritize putative genes found in preclinical models into the likely (meeting at least one) and more likely candidate (meeting both), which predisposes human genes to mutations.

Functional partitioning of a liquid-like organelle during assembly of axonemal dyneins

Ciliary motility is driven by axonemal dyneins that are assembled in the cytoplasm before deployment to cilia. Motile ciliopathy can result from defects in the dyneins themselves or from defects in factors required for their cytoplasmic pre-assembly. Recent work demonstrates that axonemal dyneins, their specific assembly factors, and broadly-acting chaperones are concentrated in liquid-like organelles in the cytoplasm called DynAPs (Dynein Axonemal Particles). Here, we use in vivo imaging in Xenopus to show that inner dynein arm (IDA) and outer dynein arm (ODA) subunits are partitioned into non-overlapping sub-regions within DynAPs. Using affinity- purification mass-spectrometry of in vivo interaction partners, we also identify novel partners for inner and outer dynein arms. Among these, we identify C16orf71/Daap1 as a novel axonemal dynein regulator. Daap1 interacts with ODA subunits, localizes specifically to the cytoplasm, is enriched in DynAPs, and is required for the deployment of ODAs to axonemes. Our work reveals a new complexity in the structure and function of a cell-type specific liquid-like organelle that is directly relevant to human genetic disease.

Generation of two hiPSC clones (MHHi019-A, MHHi019-B) from a primary ciliary dyskinesia patient carrying a homozygous deletion in the NME5 gene (c.415delA (p.Ile139Tyrfs*8))

Primary ciliary dyskinesia (PCD) is a genetic disorder characterized by defects in motile cilia and is known to occur in about 1 in 20,000 live births (Horani and Ferkol, 2018). Among the many genes associated with PCD, NME5, a gene encoding a protein involved in ciliary function, was recently reported to be involved in PCD (Anderegg et al., 2019; Cho et al., 2020). We have established two human induced pluripotent stem cell clones from a PCD patient carrying a deletion in the NME5 gene (c.415delA (p.Ile139Tyrfs*8)).

Motile ciliopathies

Motile cilia are highly complex hair-like organelles of epithelial cells lining the surface of various organ systems. Genetic mutations (usually with autosomal recessive inheritance) that impair ciliary beating cause a variety of motile ciliopathies, a heterogeneous group of rare disorders. The pathogenetic mechanisms, clinical symptoms and severity of the disease depend on the specific affected genes and the tissues in which they are expressed. Defects in the ependymal cilia can result in hydrocephalus, defects in the cilia in the fallopian tubes or in sperm flagella can cause female and male subfertility, respectively, and malfunctional motile monocilia of the left-right organizer during early embryonic development can lead to laterality defects such as situs inversus and heterotaxy. If mucociliary clearance in the respiratory epithelium is severely impaired, the disorder is referred to as primary ciliary dyskinesia, the most common motile ciliopathy. No single test can confirm a diagnosis of motile ciliopathy, which is based on a combination of tests including nasal nitric oxide measurement, transmission electron microscopy, immunofluorescence and genetic analyses, and high-speed video microscopy. With the exception of azithromycin, there is no evidence-based treatment for primary ciliary dyskinesia; therapies aim at relieving symptoms and reducing the effects of reduced ciliary motility.

A nonsense variant in NME5 causes human primary ciliary dyskinesia with radial spoke defects

Primary ciliary dyskinesia (PCD) is a genetically heterogeneous disorder characterized by defects in the function or structure of motitle cilia. In most cases, causative variants result in axonemal dynein arm anomalies, however, PCD due to radial spoke (RS) and central pair (CP) of microtubules has been rarely reported. To identify the molecular basis of PCD characterized by RS/CP defects, we performed whole exome sequencing in PCD patients with RS/CP defects. We identified a homozygous nonsense variant (c.572G>A; p.Trp191*) in NME5, which encodes a protein component of the RS neck, in one PCD patient with situs solitus. Morpholino knockdown of nme5 in zebrafish embryos resulted in motile cilia defects with phenotypes compatible with ciliopathy. This is the first study to show NME5 as a PCD-causative gene in humans. Our findings indicate that NME5 screening should be considered for PCD patients with RS/CP defects.

Rsph4a is essential for the triplet radial spoke head assembly of the mouse motile cilia

Motile cilia/flagella are essential for swimming and generating extracellular fluid flow in eukaryotes. Motile cilia harbor a 9+2 arrangement consisting of nine doublet microtubules with dynein arms at the periphery and a pair of singlet microtubules at the center (central pair). In the central system, the radial spoke has a T-shaped architecture and regulates the motility and motion pattern of cilia. Recent cryoelectron tomography data reveal three types of radial spokes (RS1, RS2, and RS3) in the 96 nm axoneme repeat unit; however, the molecular composition of the third radial spoke, RS3 is unknown. In human pathology, it is well known mutation of the radial spoke head-related genes causes primary ciliary dyskinesia (PCD) including respiratory defect and infertility. Here, we describe the role of the primary ciliary dyskinesia protein Rsph4a in the mouse motile cilia. Cryoelectron tomography reveals that the mouse trachea cilia harbor three types of radial spoke as with the other vertebrates and that all triplet spoke heads are lacking in the trachea cilia of Rsph4a-deficient mice. Furthermore, observation of ciliary movement and immunofluorescence analysis indicates that Rsph4a contributes to the generation of the planar beating of motile cilia by building the distal architecture of radial spokes in the trachea, the ependymal tissues, and the oviduct. Although detailed mechanism of RSs assembly remains unknown, our results suggest Rsph4a is a generic component of radial spoke heads, and could explain the severe phenotype of human PCD patients with RSPH4A mutation.

Recessive missense LAMP3 variant associated with defect in lamellar body biogenesis and fatal neonatal interstitial lung disease in dogs

Neonatal interstitial lung diseases due to abnormal surfactant biogenesis are rare in humans and have never been reported as a spontaneous disorder in animals. We describe here a novel lung disorder in Airedale Terrier (AT) dogs with clinical symptoms and pathology similar to the most severe neonatal forms of human surfactant deficiency. Lethal hypoxic respiratory distress and failure occurred within the first days or weeks of life in the affected puppies. Transmission electron microscopy of the affected lungs revealed maturation arrest in the formation of lamellar bodies (LBs) in the alveolar epithelial type II (AECII) cells. The secretory organelles were small and contained fewer lamellae, often in combination with small vesicles surrounded by an occasionally disrupted common limiting membrane. A combined approach of genome-wide association study and whole exome sequencing identified a recessive variant, c.1159G>A, p.(E387K), in LAMP3, a limiting membrane protein of the cytoplasmic surfactant organelles in AECII cells. The substitution resides in the LAMP domain adjacent to a conserved disulfide bond. In summary, this study describes a novel interstitial lung disease in dogs, identifies a new candidate gene for human surfactant dysfunction and brings important insights into the essential role of LAMP3 in the process of the LB formation.

A comprehensive biomedical variant catalogue based on whole genome sequences of 582 dogs and eight wolves

The domestic dog serves as an excellent model to investigate the genetic basis of disease. More than 400 heritable traits analogous to human diseases have been described in dogs. To further canine medical genetics research, we established the Dog Biomedical Variant Database Consortium (DBVDC) and present a comprehensive list of functionally annotated genome variants that were identified with whole genome sequencing of 582 dogs from 126 breeds and eight wolves. The genomes used in the study have a minimum coverage of 10× and an average coverage of ~24×. In total, we identified 23 133 692 single-nucleotide variants (SNVs) and 10 048 038 short indels, including 93% undescribed variants. On average, each individual dog genome carried ∼4.1 million single-nucleotide and ~1.4 million short-indel variants with respect to the reference genome assembly. About 2% of the variants were located in coding regions of annotated genes and loci. Variant effect classification showed 247 141 SNVs and 99 562 short indels having moderate or high impact on 11 267 protein-coding genes. On average, each genome contained heterozygous loss-of-function variants in 30 potentially embryonic lethal genes and 97 genes associated with developmental disorders. More than 50 inherited disorders and traits have been unravelled using the DBVDC variant catalogue, enabling genetic testing for breeding and diagnostics. This resource of annotated variants and their corresponding genotype frequencies constitutes a highly useful tool for the identification of potential variants causative for rare inherited disorders in dogs.

AKNA Frameshift Variant in Three Dogs with Recurrent Inflammatory Pulmonary Disease

We investigated three related Rough Collies with recurrent inflammatory pulmonary disease. The clinical symptoms were similar to primary ciliary dyskinesia (PCD). However, the affected dogs did not carry any known pathogenic PCD variants. Pedigree analysis suggested a recessive mode of inheritance. Combined linkage and homozygosity mapping in three cases and seven non-affected family members delineated 19 critical intervals on 10 chromosomes comprising a total of 99 Mb. The genome of one affected dog was sequenced and compared to 601 control genomes. We detected only a single private homozygous protein-changing variant in the critical intervals. The detected variant was a 4 bp deletion, c.2717_2720delACAG, in the AKNA gene encoding the AT-hook transcription factor. It causes a frame-shift introducing a premature stop codon and truncates 37% of the open reading frame, p.(Asp906Alafs*173). We genotyped 88 Rough Collies consisting of family members and unrelated individuals. All three available cases were homozygous for the mutant allele and all 85 non-affected dogs were either homozygous wildtype (n = 67) or heterozygous (n = 18). AKNA modulates inflammatory immune responses. Akna-/- knockout mice die shortly after birth due to systemic autoimmune inflammatory processes including lung inflammation that is accompanied by enhanced leukocyte infiltration and alveolar destruction. The perfect genotype-phenotype association and the comparative functional data strongly suggest that the detected AKNA:c.2717_2720delACAG variant caused the observed severe airway inflammation in the investigated dogs. Our findings enable genetic testing, which can be used to avoid the unintentional breeding of affected puppies.