Automated Proteomics Workflows for High-Throughput Library Generation and Biomarker Detection Using Data-Independent Acquisition.

Sequential window acquisition of all theoretical mass spectra-mass spectrometry underpinned by advanced bioinformatics offers a framework for comprehensive analysis of proteomes and the discovery of robust biomarkers. However, the lack of a generic sample preparation platform to tackle the heterogeneity of material collected from different sources may be a limiting factor to the broad application of this technique. We have developed universal and fully automated workflows using a robotic sample preparation platform, which enabled in-depth and reproducible proteome coverage and characterization of bovine and ovine specimens representing healthy animals and a model of myocardial infarction. High correlation (R2 = 0.85) between sheep proteomics and transcriptomics datasets validated the developments. The findings suggest that automated workflows can be employed for various clinical applications across different animal species and animal models of health and disease.

Genome architecture and diverged selection shaping pattern of genomic differentiation in wild barley.

Divergent selection of populations in contrasting environments leads to functional genomic divergence. However, the genomic architecture underlying heterogeneous genomic differentiation remains poorly understood. Here, we de novo assembled two high-quality wild barley (Hordeum spontaneum K. Koch) genomes and examined genomic differentiation and gene expression patterns under abiotic stress in two populations. These two populations had a shared ancestry and originated in close geographic proximity but experienced different selective pressures due to their contrasting micro-environments. We identified structural variants that may have played significant roles in affecting genes potentially associated with well-differentiated phenotypes such as flowering time and drought response between two wild barley genomes. Among them, a 29-bp insertion into the promoter region formed a cis-regulatory element in the HvWRKY45 gene, which may contribute to enhanced tolerance to drought. A single SNP mutation in the promoter region may influence HvCO5 expression and be putatively linked to local flowering time adaptation. We also revealed significant genomic differentiation between the two populations with ongoing gene flow. Our results indicate that SNPs and small SVs link to genetic differentiation at the gene level through local adaptation and are maintained through divergent selection. In contrast, large chromosome inversions may have shaped the heterogeneous pattern of genomic differentiation along the chromosomes by suppressing chromosome recombination and gene flow. Our research offers novel insights into the genomic basis underlying local adaptation and provides valuable resources for the genetic improvement of cultivated barley.

First detection and complete genome sequence of a new potexvirus naturally infecting Adenium obesum.

Here, we report the detection and complete genome sequence of a novel potexvirus, tentatively named "Adenium obesum virus X" (AobVX), isolated from Adenium obesum, that was sent for virus screening at Australian Government post-entry quarantine (PEQ) facilities after being imported into Australia from China. The AobVX genome is 6781 nucleotides in length excluding the poly(A) tail and is predicted to encode conserved potexvirus proteins and sequence motifs across five open reading frames. The RNA-dependent RNA polymerase of this virus shares the highest amino acid sequence similarity with that of nerine potexvirus 1 (58.7% identity) and nerine virus X (58.58% identity). This is the first report of a positive-sense single-stranded RNA virus in A. obesum related to members of the genus Potexvirus in the family Alphaflexiviridae.

A chromosome conformation capture ordered sequence of the barley genome.

Cereal grasses of the Triticeae tribe have been the major food source in temperate regions since the dawn of agriculture. Their large genomes are characterized by a high content of repetitive elements and large pericentromeric regions that are virtually devoid of meiotic recombination. Here we present a high-quality reference genome assembly for barley (Hordeum vulgare L.). We use chromosome conformation capture mapping to derive the linear order of sequences across the pericentromeric space and to investigate the spatial organization of chromatin in the nucleus at megabase resolution. The composition of genes and repetitive elements differs between distal and proximal regions. Gene family analyses reveal lineage-specific duplications of genes involved in the transport of nutrients to developing seeds and the mobilization of carbohydrates in grains. We demonstrate the importance of the barley reference sequence for breeding by inspecting the genomic partitioning of sequence variation in modern elite germplasm, highlighting regions vulnerable to genetic erosion.

Complete genome sequence of a novel potyvirus infecting Miscanthus sinensis (silver grass).

Here, we describe the full-length genome sequence of a novel potyvirus, tentatively named "Miscanthus sinensis mosaic virus" (MsiMV), isolated from Miscanthus sinensis (silver grass) held in a post-entry quarantine facility after being imported into Western Australia, Australia. The MsiMV genome is 9604 nucleotides (nt) in length, encoding a 3071-amino-acid (aa) polyprotein with conserved sequence motifs. The MsiMV genome is most closely related to that of sorghum mosaic virus (SrMV), with 74% nt and 78.5% aa sequence identity to the SrMV polyprotein region. Phylogenetic analysis based on the polyprotein grouped MsiMV with SrMV, sugarcane mosaic virus (SCMV), and maize dwarf mosaic virus (MDMV). This is the first report of a novel monopartite ssRNA virus in Miscanthus sinensis related to members of the genus Potyvirus in the family Potyviridae.

IsomiR-eQTL: A Cancer-Specific Expression Quantitative Trait Loci Database of miRNAs and Their Isoforms.

The identification of expression quantitative trait loci (eQTL) is an important component in efforts to understand how genetic variants influence disease risk. MicroRNAs (miRNAs) are short noncoding RNA molecules capable of regulating the expression of several genes simultaneously. Recently, several novel isomers of miRNAs (isomiRs) that differ slightly in length and sequence composition compared to their canonical miRNAs have been reported. Here we present isomiR-eQTL, a user-friendly database designed to help researchers find single nucleotide polymorphisms (SNPs) that can impact miRNA (miR-eQTL) and isomiR expression (isomiR-eQTL) in 30 cancer types. The isomiR-eQTL includes a total of 152,671 miR-eQTLs and 2,390,805 isomiR-eQTLs at a false discovery rate (FDR) of 0.05. It also includes 65,733 miR-eQTLs overlapping known cancer-associated loci identified through genome-wide association studies (GWAS). To the best of our knowledge, this is the first study investigating the impact of SNPs on isomiR expression at the genome-wide level. This database may pave the way for researchers toward finding a model for personalised medicine in which miRNAs, isomiRs, and genotypes are utilised.

Identification of Candidate mRNA Isoforms for Prostate Cancer-Risk SNPs Utilizing Iso-eQTL and sQTL Methods.

Single nucleotide polymorphisms (SNPs) impacting the alternative splicing (AS) process (sQTLs) or isoform expression (iso-eQTL) are implicated as important cancer regulatory elements. To find the sQTL and iso-eQTL, we retrieved prostate cancer (PrCa) tissue RNA-seq and genotype data originating from 385 PrCa European patients from The Cancer Genome Atlas. We conducted RNA-seq analysis with isoform-based and splice event-based approaches. The MatrixEQTL was used to identify PrCa-associated sQTLs and iso-eQTLs. The overlap between sQTL and iso-eQTL with GWAS loci and those that are differentially expressed between cancer and normal tissue were identified. The cis-acting associations (FDR < 0.05) for PrCa-risk SNPs identified 42, 123, and 90 PrCa-associated cassette exons, intron retention, and mRNA isoforms belonging to 25, 95, and 83 genes, respectively; while assessment of trans-acting association (FDR < 0.05) yielded 59, 65, and 196 PrCa-associated cassette exons, intron retention and mRNA isoforms belonging to 35, 55, and 181 genes, respectively. The results suggest that functional PrCa-associated SNPs can play a role in PrCa genesis by making an important contribution to the dysregulation of AS and, consequently, impacting the expression of the mRNA isoforms.

Bovine extracellular vesicles contaminate human extracellular vesicles produced in cell culture conditioned medium when 'exosome-depleted serum' is utilised.

Extracellular vesicles (EVs) are important intercellular communication messengers. Half of the published studies in the field are in vitro cell culture based in which bovine serum in various concentrations and forms is used to facilitate the production of extracellular vesicles. 'Exosome depleted serum' is the type of bovine serum most widely used in the production of human EVs. Herein, we demonstrate that, despite the initial caution raised in 2014 about the persistence of bovine EVs, 'exosome depleted serum' was still used in 46% of publications on human or rodent EVs between 2015 and 2019. Using nanoparticle tracking analysis combined with detergent lysis of vesicles as well as bovine CD9 ELISA, we show that there were approximately 5.33 x 107/mL of bovine EVs remaining in the 'exosome depleted serum'. Importantly, the 'exosome depleted serum' was relatively enriched in small EVs by approximately 2.7-fold relative to the large EVs compared to that in the original serum. Specifically, the percentage of small EVs in total vesicles had increased from the original 48% in the serum before ultracentrifugation to 92% in the 'exosome depleted serum'. Furthermore, the pervasive bovine EVs carried over by the 'exosome depleted serum', even when the lowest concentration (0.5%) was used in cell culture, resulted in a significant contamination of human EVs in cell culture conditioned medium. Our findings indicate that the use 'exosome depleted serum' in cell culture-based studies may introduce artefacts into research examining the function of human and rodent EVs, in particular those involving EV miRNA. Thus, we appeal to the researchers in the EV field to seriously reconsider the practice of using 'exosome depleted serum' in the production of human and other mammalian EVs in vitro.

Development of a CD63 Aptamer for Efficient Cancer Immunochemistry and Immunoaffinity-Based Exosome Isolation.

CD63, a member of transmembrane-4-superfamily of tetraspanin proteins and a highly N-glycosylated type III lysosomal membrane protein, is known to regulate malignancy of various types of cancers such as melanoma and breast cancer and serves as a potential marker for cancer detection. Recently, its important role as a classic exosome marker was also emphasized. In this work, via using a magnetic bead-based competitive SELEX (systematic evolution of ligands by exponential enrichment) procedure and introducing a 0.5M NaCl as elution buffer, we identified two DNA aptamers (CD63-1 and CD63-2) with high affinity and specificity to CD63 protein (Kd = 38.71nM and 78.43, respectively). Furthermore, CD63-1 was found to be efficient in binding CD63 positive cells, including breast cancer MDA-MB-231 cells and CD63-overexpressed HEK293T cells, with a medium binding affinity (Kd~ 100 nM) as assessed by flow cytometry. When immunostaining assay was performed using clinical breast cancer biopsy, the CD63-1 aptamer demonstrated a comparable diagnostic efficacy for CD63 positive breast cancer with commercial antibodies. After developing a magnetic bead-based exosome immunoaffinity separation system using CD63-1 aptamer, it was found that this bead-based system could effectively isolate exosomes from both MDA-MB-231 and HT29 cell culture medium. Importantly, the introduction of the NaCl elution in this work enabled the isolation of native exosomes via a simple 0.5M NaCl incubation step. Based on these results, we firmly believe that the developed aptamers could be useful towards efficient isolation of native state exosomes from clinical samples and various theranostic applications for CD63-positive cancers.

Identification of a novel vitivirus from grapevines in New Zealand.

We report a sequence of a novel vitivirus from Vitis vinifera obtained using two high-throughput sequencing (HTS) strategies on RNA. The initial discovery from small-RNA sequencing was confirmed by HTS of the total RNA and Sanger sequencing. The new virus has a genome structure similar to the one reported for other vitiviruses, with five open reading frames (ORFs) coding for the conserved domains described for members of that genus. Phylogenetic analysis of the complete genome sequence confirmed its affiliation to the genus Vitivirus, with the closest described viruses being grapevine virus E (GVE) and Agave tequilana leaf virus (ATLV). However, the virus we report is distinct and shares only 51% amino acid sequence identity with GVE in the replicase polyprotein and 66.8% amino acid sequence identity with ATLV in the coat protein. This is well below the threshold determined by the ICTV for species demarcation, and we propose that this virus represents a new species. It is provisionally named "grapevine virus G".

Grapevine virus I, a putative new vitivirus detected in co-infection with grapevine virus G in New Zealand.

A novel virus, with characteristics of viruses classified within the genus Vitivirus, was identified from a sample of Vitis vinifera cv. Chardonnay in New Zealand. The virus was detected with high throughput sequencing (small RNA and total RNA) and its sequence was confirmed by Sanger sequencing. Its genome is 7507 nt long (excluding the polyA tail) with an organisation similar to that described for other classifiable members of the genus Vitivirus. The closest relative of the virus is grapevine virus E (GVE) with 65% aa identity in ORF1 (65% nt identity) and 63% aa identity in the coat protein (66% nt identity). The relationship with GVE was confirmed with phylogenetic analysis, showing the new virus branching with GVE, Agave tequilina leaf virus and grapevine virus G (GVG). A limited survey revealed the presence of this virus in multiple plants from the same location where the newly described GVG was discovered, and in most cases both viruses were detected as co-infections. The genetic characteristics of this virus suggest it represents an isolate of a new species within the genus Vitivirus and following the current nomenclature, we propose the name "Grapevine virus I".

De novo assembly of honey bee RNA viral genomes by tapping into the innate insect antiviral response pathway.

Bee pollination is critical for improving productivity of one third of all plants or plant products consumed by humans. The health of honey bees is in decline in many countries worldwide, and RNA viruses together with other biological, environmental and anthropogenic factors have been identified as the main causes. The rapid genetic variation of viruses represents a challenge for diagnosis. Thus, application of deep sequencing methods for detection and analysis of viruses has increased over the last years. In this study, we leverage from the innate Dicer-2 mediated antiviral response against viruses to reconstruct complete viral genomes using virus-derived small interfering RNAs (vsiRNAs). Symptomatic A. mellifera larvae collected from hives free of Colony Collapse Disorder (CCD) and the parasitic Varroa mite (Varroa destructor) were used to generate more than 107 million small RNA reads. We show that de novo assembly of insect viral sequences is less fragmented using only 22 nt long vsiRNAs rather than a combination of 21-22 nt small RNAs. Our results show that A. mellifera larvae activate the RNAi immune response in the presence of Sacbrood virus (SBV). We assembled three SBV genomes from three individual larvae from different hives in a single apiary, with 1-2% nucleotide sequence variability among them. We found 3-4% variability between SBV genomes generated in this study and earlier published Australian variants suggesting the presence of different SBV quasispecies within the country.

An internet-based bioinformatics toolkit for plant biosecurity diagnosis and surveillance of viruses and viroids.

BACKGROUND: Detection and preventing entry of exotic viruses and viroids at the border is critical for protecting plant industries trade worldwide. Existing post entry quarantine screening protocols rely on time-consuming biological indicators and/or molecular assays that require knowledge of infecting viral pathogens. Plants have developed the ability to recognise and respond to viral infections through Dicer-like enzymes that cleave viral sequences into specific small RNA products. Many studies reported the use of a broad range of small RNAs encompassing the product sizes of several Dicer enzymes involved in distinct biological pathways. Here we optimise the assembly of viral sequences by using specific small RNA subsets. RESULTS: We sequenced the small RNA fractions of 21 plants held at quarantine glasshouse facilities in Australia and New Zealand. Benchmarking of several de novo assembler tools yielded SPAdes using a kmer of 19 to produce the best assembly outcomes. We also found that de novo assembly using 21-25 nt small RNAs can result in chimeric assemblies of viral sequences and plant host sequences. Such non-specific assemblies can be resolved by using 21-22 nt or 24 nt small RNAs subsets. Among the 21 selected samples, we identified contigs with sequence similarity to 18 viruses and 3 viroids in 13 samples. Most of the viruses were assembled using only 21-22 nt long virus-derived siRNAs (viRNAs), except for one Citrus endogenous pararetrovirus that was more efficiently assembled using 24 nt long viRNAs. All three viroids found in this study were fully assembled using either 21-22 nt or 24 nt viRNAs. Optimised analysis workflows were customised within the Yabi web-based analytical environment. We present a fully automated viral surveillance and diagnosis web-based bioinformatics toolkit that provides a flexible, user-friendly, robust and scalable interface for the discovery and diagnosis of viral pathogens. CONCLUSIONS: We have implemented an automated viral surveillance and diagnosis (VSD) bioinformatics toolkit that produces improved viruses and viroid sequence assemblies. The VSD toolkit provides several optimised and reusable workflows applicable to distinct viral pathogens. We envisage that this resource will facilitate the surveillance and diagnosis viral pathogens in plants, insects and invertebrates.

Acetylcholinesterase 1 in populations of organophosphate-resistant North American strains of the cattle tick, Rhipicephalus microplus (Acari: Ixodidae).

Rhipicephalus microplus, the cattle fever tick, is a global economic problem to the cattle industry due to direct infestation of cattle and pathogens transmitted during feeding. Cattle fever tick outbreaks continue to occur along the Mexico-US border even though the tick has been eradicated from the USA. The organophosphate (OP) coumaphos targets acetylcholinesterase (AChE) and is the approved acaricide for eradicating cattle fever tick outbreaks. There is evidence for coumaphos resistance developing in cattle ticks in Mexico, and OP-resistant R. microplus ticks were discovered in outbreak populations of Texas in 2005. The molecular basis of coumaphos resistance is not known, and our study was established to gather further information on whether AChE1 is involved in the resistance mechanism. We also sought information on allele diversity in tick populations with different levels of coumaphos resistance. The overarching project goal was to define OP resistance-associated gene mutations such that a DNA-based diagnostic assay could be developed to assist the management of resistance. Three different AChE transcripts have been reported in R. microplus, and supporting genomic and transcriptomic data are available at CattleTickBase. Here, we report the complete R. microplus AChE1 gene ascertained by sequencing a bacterial artificial chromosome clone containing the entire coding region and the flanking 5' and 3' regions. We also report AChE1 sequences of larval ticks from R. microplus strains having different sensitivities to OP. To accomplish this, we sequenced a 669-bp region of the AChE1 gene corresponding to a 223 amino acid region of exon 2 to assess alleles in seven strains of R. microplus with varying OP resistance phenotypes. We identified 72 AChE1 sequence variants, 2 of which are strongly associated with OP-resistant phenotypes. Esterase-like sequences from the R. microplus transcriptome RmiTr Version 1.0 were compared to the available sequence databases to identify other transcripts with similarity to AChE1.

Genome sequencing and analysis of the paclitaxel-producing endophytic fungus Penicillium aurantiogriseum NRRL 62431.

BACKGROUND: Paclitaxel (Taxol™) is an important anticancer drug with a unique mode of action. The biosynthesis of paclitaxel had been considered restricted to the Taxus species until it was discovered in Taxomyces andreanae, an endophytic fungus of T. brevifolia. Subsequently, paclitaxel was found in hazel (Corylus avellana L.) and in several other endophytic fungi. The distribution of paclitaxel in plants and endophytic fungi and the reported sequence homology of key genes in paclitaxel biosynthesis between plant and fungi species raises the question about whether the origin of this pathway in these two physically associated groups could have been facilitated by horizontal gene transfer. RESULTS: The ability of the endophytic fungus of hazel Penicillium aurantiogriseum NRRL 62431 to independently synthesize paclitaxel was established by liquid chromatography-mass spectrometry and proton nuclear magnetic resonance. The genome of Penicillium aurantiogriseum NRRL 62431 was sequenced and gene candidates that may be involved in paclitaxel biosynthesis were identified by comparison with the 13 known paclitaxel biosynthetic genes in Taxus. We found that paclitaxel biosynthetic gene candidates in P. aurantiogriseum NRRL 62431 have evolved independently and that horizontal gene transfer between this endophytic fungus and its plant host is unlikely. CONCLUSIONS: Our findings shed new light on how paclitaxel-producing endophytic fungi synthesize paclitaxel, and will facilitate metabolic engineering for the industrial production of paclitaxel from fungi.

A novel saliva-based miRNA profile to diagnose and predict oral cancer.

Oral cancer (OC) is the most common form of head and neck cancer. Despite the high incidence and unfavourable patient outcomes, currently, there are no biomarkers for the early detection of OC. This study aims to discover, develop, and validate a novel saliva-based microRNA signature for early diagnosis and prediction of OC risk in oral potentially malignant disorders (OPMD). The Cancer Genome Atlas (TCGA) miRNA sequencing data and small RNA sequencing data of saliva samples were used to discover differentially expressed miRNAs. Identified miRNAs were validated in saliva samples of OC (n = 50), OPMD (n = 52), and controls (n = 60) using quantitative real-time PCR. Eight differentially expressed miRNAs (miR-7-5p, miR-10b-5p, miR-182-5p, miR-215-5p, miR-431-5p, miR-486-3p, miR-3614-5p, and miR-4707-3p) were identified in the discovery phase and were validated. The efficiency of our eight-miRNA signature to discriminate OC and controls was: area under curve (AUC): 0.954, sensitivity: 86%, specificity: 90%, positive predictive value (PPV): 87.8% and negative predictive value (NPV): 88.5% whereas between OC and OPMD was: AUC: 0.911, sensitivity: 90%, specificity: 82.7%, PPV: 74.2% and NPV: 89.6%. We have developed a risk probability score to predict the presence or risk of OC in OPMD patients. We established a salivary miRNA signature that can aid in diagnosing and predicting OC, revolutionising the management of patients with OPMD. Together, our results shed new light on the management of OC by salivary miRNAs to the clinical utility of using miRNAs derived from saliva samples.

A prioritization analysis of disease association by data-mining of functional annotation of human genes.

Complex diseases result from contributions of multiple genes that act in concert through pathways. Here we present a method to prioritize novel candidates of disease-susceptibility genes depending on the biological similarities to the known disease-related genes. The extent of disease-susceptibility of a gene is prioritized by analyzing seven features of human genes captured in H-InvDB. Taking rheumatoid arthritis (RA) and prostate cancer (PC) as two examples, we evaluated the efficiency of our method. Highly scored genes obtained included TNFSF12 and OSM as candidate disease genes for RA and PC, respectively. Subsequent characterization of these genes based upon an extensive literature survey reinforced the validity of these highly scored genes as possible disease-susceptibility genes. Our approach, Prioritization ANalysis of Disease Association (PANDA), is an efficient and cost-effective method to narrow down a large set of genes into smaller subsets that are most likely to be involved in the disease pathogenesis.

Does the registry speak your language? A case study of the Global Angelman Syndrome Registry.

Global disease registries are critical to capturing common patient related information on rare illnesses, allowing patients and their families to provide information about their condition in a safe, accessible, and engaging manner that enables researchers to undertake critical research aimed at improving outcomes. Typically, English is the default language of choice for these global digital health platforms. Unfortunately, language barriers can significantly inhibit participation from non-English speaking participants. In addition, there is potential for compromises in data quality and completeness. In contrast, multinational commercial entities provide access to their websites in the local language of the country they are operating in, and often provide multiple options reflecting ethnic diversity. This paper presents a case study of how the Global Angelman Syndrome Registry (GASR) has used a novel approach to enable multiple language translations for its website. Using a "semi-automated language translation" approach, the GASR, which was originally launched in English in September 2016, is now available in several other languages. In 2020, the GASR adopted a novel approach using crowd-sourcing and machine translation tools leading to the availability of the GASR in Spanish, Traditional Chinese, Italian, and Hindi. As a result, enrolments increased by 124% percent for Spain, 67% percent for Latin America, 46% percent for Asia, 24% for Italy, and 43% for India. We describe our approach here, which we believe presents an opportunity for cost-effective and timely translations responsive to changes to the registry and helps build and maintain engagement with global disease communities.

Salivary micro RNAs as biomarkers for oropharyngeal cancer.

BACKGROUND: Despite the rising incidence, particularly of the human papillomavirus (HPV)-associated fraction of oropharyngeal cancer (OPC), there are no early detection methods for OPC. Considering the close association between saliva and head and neck cancers, this study was designed to investigate salivary micro RNA (miRNAs) associated with OPC, especially focusing on HPV-positive OPC. METHODS: Saliva was collected from OPC patients at diagnosis and patients were clinically followed up ≤5 years. Salivary small RNA isolated from HPV-positive OPC patients (N = 6), and HPV-positive (N = 4) and negative controls (N = 6) were analysed by next-generation sequencing to identify dysregulated miRNAs. Discovered miRNAs were validated by quantitative PCR using two different assays in a separate cohort of patients (OPC = 91, controls = 92). The relative expression was calculated considering SNORD-96A as the normalizer. Candidate miRNAs with diagnostic and prognostic potential were evaluated by generalized logistic regression. RESULTS: A panel consisting of nine miRNAs was identified to have the best diagnostic performance to discriminate HPV-positive OPC from HPV-positive controls (AUC- validation-1 = 94.8%, validation-2 = 98%). Further, a panel consisting of six miRNAs were identified to discriminate OPC from controls regardless of the HPV status (AUC- validation-1 = 77.2%, validation-2 = 86.7%). In addition, the downregulation of hsa-miR-7-5p was significantly associated with poor overall survival of OPC patients (HR = 0.638). A panel consisting of nine miRNAs were identified for the prediction of the overall survival of the OPC patients (log-rank test-p = 0.0008). CONCLUSION: This study highlights that salivary miRNAs can play an essential role in the detection and prognostication of OPC.

A novel mutation causing nephronophthisis in the Lewis polycystic kidney rat localises to a conserved RCC1 domain in Nek8.

BACKGROUND: Nephronophthisis (NPHP) as a cause of cystic kidney disease is the most common genetic cause of progressive renal failure in children and young adults. NPHP is characterized by abnormal and/or loss of function of proteins associated with primary cilia. Previously, we characterized an autosomal recessive phenotype of cystic kidney disease in the Lewis Polycystic Kidney (LPK) rat. RESULTS: In this study, quantitative trait locus analysis was used to define a ~1.6 Mbp region on rat chromosome 10q25 harbouring the lpk mutation. Targeted genome capture and next-generation sequencing of this region identified a non-synonymous mutation R650C in the NIMA (never in mitosis gene a)- related kinase 8 ( Nek8) gene. This is a novel Nek8 mutation that occurs within the regulator of chromosome condensation 1 (RCC1)-like region of the protein. Specifically, the R650C substitution is located within a G[QRC]LG repeat motif of the predicted seven bladed beta-propeller structure of the RCC1 domain. The rat Nek8 gene is located in a region syntenic to portions of human chromosome 17 and mouse 11. Scanning electron microscopy confirmed abnormally long cilia on LPK kidney epithelial cells, and fluorescence immunohistochemistry for Nek8 protein revealed altered cilia localisation. CONCLUSIONS: When assessed relative to other Nek8 NPHP mutations, our results indicate the whole propeller structure of the RCC1 domain is important, as the different mutations cause comparable phenotypes. This study establishes the LPK rat as a novel model system for NPHP and further consolidates the link between cystic kidney disease and cilia proteins.