Single nuclei RNA-seq reveals a medium spiny neuron glutamate excitotoxicity signature prior to the onset of neuronal death in an ovine Huntington's disease model.

Huntington's disease (HD) is a neurodegenerative genetic disorder caused by an expansion in the CAG repeat tract of the huntingtin (HTT) gene resulting in behavioural, cognitive, and motor defects. Current knowledge of disease pathogenesis remains incomplete, and no disease course-modifying interventions are in clinical use. We have previously reported the development and characterisation of the OVT73 transgenic sheep model of HD. The 73 polyglutamine repeat is somatically stable and therefore likely captures a prodromal phase of the disease with an absence of motor symptomatology even at 5-years of age and no detectable striatal cell loss. To better understand the disease-initiating events we have undertaken a single nuclei transcriptome study of the striatum of an extensively studied cohort of 5-year-old OVT73 HD sheep and age matched wild-type controls. We have identified transcriptional upregulation of genes encoding N-methyl-D-aspartate (NMDA), α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) and kainate receptors in medium spiny neurons, the cell type preferentially lost early in HD. Further, we observed an upregulation of astrocytic glutamate uptake transporters and medium spiny neuron GABAA receptors, which may maintain glutamate homeostasis. Taken together, these observations support the glutamate excitotoxicity hypothesis as an early neurodegeneration cascade-initiating process but the threshold of toxicity may be regulated by several protective mechanisms. Addressing this biochemical defect early may prevent neuronal loss and avoid the more complex secondary consequences precipitated by cell death.

Association between maternal depression symptoms and child telomere length.

The biological mechanisms that explain how adverse early life events influence adult disease risk are poorly understood. One proposed mechanism is via the induction of accelerated biological aging, for which telomere length is considered a biomarker. We aimed to determine if maternal depression pre- and post-partum was associated with telomere length in children at 4 years of age (n = 4299). Mothers completed structured questionnaires assessing depression during pregnancy (Edinburgh Depression Scale), at 9 months (Edinburgh Depression Scale), and at 54 months postpartum (Patient Health Questionnaire 9). Regression methods were used to investigate the relationship between telomere length (DNA from saliva) and maternal depression score recorded at each stage. Significant covariates included in the final model were: maternal age at pregnancy; child sex; child ethnicity; gestational age group, and rurality group. Child telomere length was found to be longer if their mother had a higher depression score at both postpartum time points tested (9 months of age; coefficient 0.003, SE = 0.001, P = 0.01, 54 months of age; coefficient 0.003, SE = 0.002, P = 0.02). Although these findings seem paradoxical, increased telomere length may be an adaptive response to early life stressors. We propose several testable hypotheses for these results and to determine if the positive association between depression and telomere length is a developmental adaptation or an indirect consequence of environmental factors.

A common regulatory haplotype doubles lactoferrin concentration in milk.

BACKGROUND: Bovine lactoferrin (Lf) is an iron absorbing whey protein with antibacterial, antiviral, and antifungal activity. Lactoferrin is economically valuable and has an extremely variable concentration in milk, partly driven by environmental influences such as milking frequency, involution, or mastitis. A significant genetic influence has also been previously observed to regulate lactoferrin content in milk. Here, we conducted genetic mapping of lactoferrin protein concentration in conjunction with RNA-seq, ChIP-seq, and ATAC-seq data to pinpoint candidate causative variants that regulate lactoferrin concentrations in milk. RESULTS: We identified a highly-significant lactoferrin protein quantitative trait locus (pQTL), as well as a cis lactotransferrin (LTF) expression QTL (cis-eQTL) mapping to the LTF locus. Using ChIP-seq and ATAC-seq datasets representing lactating mammary tissue samples, we also report a number of regions where the openness of chromatin is under genetic influence. Several of these also show highly significant QTL with genetic signatures similar to those highlighted through pQTL and eQTL analysis. By performing correlation analysis between these QTL, we revealed an ATAC-seq peak in the putative promotor region of LTF, that highlights a set of 115 high-frequency variants that are potentially responsible for these effects. One of the 115 variants (rs110000337), which maps within the ATAC-seq peak, was predicted to alter binding sites of transcription factors known to be involved in lactation-related pathways. CONCLUSIONS: Here, we report a regulatory haplotype of 115 variants with conspicuously large impacts on milk lactoferrin concentration. These findings could enable the selection of animals for high-producing specialist herds.

ICARUS v3, a massively scalable web server for single-cell RNA-seq analysis of millions of cells.

MOTIVATION: In recent years, improvements in throughput of single-cell RNA-seq have resulted in a significant increase in the number of cells profiled. The generation of single-cell RNA-seq datasets comprising >1 million cells is becoming increasingly common, giving rise to demands for more efficient computational workflows. RESULTS: We present an update to our single-cell RNA-seq analysis web server application, ICARUS (available at https://launch.icarus-scrnaseq.cloud.edu.au) that allows effective analysis of large-scale single-cell RNA-seq datasets. ICARUS v3 utilizes the geometric cell sketching method to subsample cells from the overall dataset for dimensionality reduction and clustering that can be then projected to the large dataset. We then extend this functionality to select a representative subset of cells for downstream data analysis applications including differential expression analysis, gene co-expression network construction, gene regulatory network construction, trajectory analysis, cell-cell communication inference, and cell cluster associations to GWAS traits. We demonstrate analysis of single-cell RNA-seq datasets using ICARUS v3 of 1.3 million cells completed within the hour. AVAILABILITY AND IMPLEMENTATION: ICARUS is available at https://launch.icarus-scrnaseq.cloud.edu.au.

Somatic CAG Repeat Stability in a Transgenic Sheep Model of Huntington's Disease.

Somatic instability of the huntingtin (HTT) CAG repeat mutation modifies age-at-onset of Huntington's disease (HD). Understanding the mechanism and pathogenic consequences of instability may reveal therapeutic targets. Using small-pool PCR we analyzed CAG instability in the OVT73 sheep model which expresses a full-length human cDNA HTT transgene. Analyses of five- and ten-year old sheep revealed the transgene (CAG)69 repeat was remarkably stable in liver, striatum, and other brain tissues. As OVT73 sheep at ten years old have minimal cell death and behavioral changes, our findings support instability of the HTT expanded-CAG repeat as being required for the progression of HD.

A novel 11 base pair deletion in KMT2C resulting in Kleefstra syndrome 2.

BACKGROUND: Haploinsufficiency of the Lysine Methyltransferase 2C (KMT2C) gene results in the autosomal dominant disorder, Kleefstra syndrome 2. It is an extremely rare neurodevelopmental condition, with 14 previous reports describing varied clinical manifestations including dysmorphic features, delayed psychomotor development and delayed growth. METHODS: Here, we describe a female with global developmental delay, attention deficit disorder, dyspraxia, short stature and subtle non-specific dysmorphic features. To identify causative mutations, whole exome sequencing was performed on the proband and her younger brother with discrete clinical presentation. RESULTS: Whole exome sequencing identified a novel de novo heterozygous 11 bp deletion in KMT2C (c.1759_1769del), resulting in a frameshift mutation and early termination of the protein (p.Gln587SerfsTer7). This variant is the second-most N-terminal reported mutation, located 4171 amino acids upstream of the critical enzymatically active SET domain (required for chromatin modification and histone methylation). CONCLUSION: The majority of the other reported mutations are frameshift mutations upstream of the SET domain and are predicted to result in protein truncation. It is thought that truncation of the SET domain, results functionally in an inability to modify chromatin through histone methylation. This report expands the clinical and genetic characterisation of Kleefstra syndrome 2.

From Pathogenesis to Therapeutics: A Review of 150 Years of Huntington's Disease Research.

Huntington's disease (HD) is a debilitating neurodegenerative genetic disorder caused by an expanded polyglutamine-coding (CAG) trinucleotide repeat in the huntingtin (HTT) gene. HD behaves as a highly penetrant dominant disorder likely acting through a toxic gain of function by the mutant huntingtin protein. Widespread cellular degeneration of the medium spiny neurons of the caudate nucleus and putamen are responsible for the onset of symptomology that encompasses motor, cognitive, and behavioural abnormalities. Over the past 150 years of HD research since George Huntington published his description, a plethora of pathogenic mechanisms have been proposed with key themes including excitotoxicity, dopaminergic imbalance, mitochondrial dysfunction, metabolic defects, disruption of proteostasis, transcriptional dysregulation, and neuroinflammation. Despite the identification and characterisation of the causative gene and mutation and significant advances in our understanding of the cellular pathology in recent years, a disease-modifying intervention has not yet been clinically approved. This review includes an overview of Huntington's disease, from its genetic aetiology to clinical presentation and its pathogenic manifestation. An updated view of molecular mechanisms and the latest therapeutic developments will also be discussed.

Severe neonatal onset neuroregression with paroxysmal dystonia and apnoea: Expanding the phenotypic and genotypic spectrum of CARS2-related mitochondrial disease.

Disorders of mitochondrial function are a collectively common group of genetic diseases in which deficits in core mitochondrial translation machinery, including aminoacyl tRNA synthetases, are key players. Biallelic variants in the CARS2 gene (NM_024537.4), which encodes the mitochondrial aminoacyl-tRNA synthetase for cysteine (CARS2, mt-aaRScys; MIM*612800), result in childhood onset epileptic encephalopathy and complex movement disorder with combined oxidative phosphorylation deficiency (MIM#616672). Prior to this report, eight unique pathogenic variants in the CARS2 gene had been reported in seven individuals. Here, we describe a male who presented in the third week of life with apnoea. He rapidly deteriorated with paroxysmal dystonic crises and apnoea resulting in death at 16 weeks. He had no evidence of seizure activity or multisystem disease and had normal brain imaging. Skeletal muscle biopsy revealed a combined disorder of oxidative phosphorylation. Whole-exome sequencing identified biallelic variants in the CARS2 gene: one novel (c.1478T>C, p.Phe493Ser), and one previously reported (c.655G>A, p.Ala219Thr; rs727505361). Northern blot analysis of RNA isolated from the patient's fibroblasts confirmed a clear defect in aminoacylation of the mitochondrial tRNA for cysteine (mt-tRNACys). To our knowledge, this is the earliest reported case of CARS2 deficiency with severe, early onset dystonia and apnoea, without epilepsy.

Delineation of complex gene expression patterns in single cell RNA-seq data with ICARUS v2.0.

Complex biological traits and disease often involve patterns of gene expression that can be characterised and examined. Here we present ICARUS v2.0, an update to our single cell RNA-seq analysis web server with additional tools to investigate gene networks and understand core patterns of gene regulation in relation to biological traits. ICARUS v2.0 enables gene co-expression analysis with MEGENA, transcription factor regulated network identification with SCENIC, trajectory analysis with Monocle3, and characterisation of cell-cell communication with CellChat. Cell cluster gene expression profiles may be examined against Genome Wide Association Studies with MAGMA to find significant associations with GWAS traits. Additionally, differentially expressed genes may be compared against the Drug-Gene Interaction database (DGIdb 4.0) to facilitate drug discovery. ICARUS v2.0 offers a comprehensive toolbox of the latest single cell RNA-seq analysis methodologies packed into an efficient, user friendly, tutorial style web server application (accessible at https://launch.icarus-scrnaseq.cloud.edu.au/) that enables single cell RNA-seq analysis tailored to the user's dataset.

Genomic discrimination in New Zealand health and life insurance. AGenDA: Against Genomic Discrimination in Aotearoa.

Nil.

ICARUS, an interactive web server for single cell RNA-seq analysis.

Here we present ICARUS, a web server to enable users without experience in R to undertake single cell RNA-seq analysis. The focal point of ICARUS is its intuitive tutorial-style user interface, designed to guide logical navigation through the multitude of pre-processing, analysis and visualization steps. ICARUS is easily accessible through a dedicated web server (https://launch.icarus-scrnaseq.cloud.edu.au/) and avoids installation of software on the user's computer. Notable features include the facility to apply quality control thresholds and adjust dimensionality reduction and cell clustering parameters. Data is visualized through 2D/3D UMAP and t-SNE plots and may be curated to remove potential confounders such as cell cycle heterogeneity. ICARUS offers flexible differential expression analysis with user-defined cell groups and gene set enrichment analysis to identify likely affected biological pathways. Eleven organisms including human, dog, mouse, rat, zebrafish, fruit fly, nematode, yeast, cattle, chicken and pig are currently supported. Visualization of multimodal data including those generated by CITE-seq and the 10X Genomics Multiome kit is included. ICARUS incorporates a function to save the current state of analysis avoiding computationally intensive steps during repeat analysis. The complete analysis of a typical single cell RNA-seq dataset by inexperienced users may be achieved in 1-2 h.

Proof of concept for multiplex amplicon sequencing for mutation identification using the MinION nanopore sequencer.

Rapid, cost-effective identification of genetic variants in small candidate genomic regions remains a challenge, particularly for less well equipped or lower throughput laboratories. The application of Oxford Nanopore Technologies' MinION sequencer has the potential to fulfil this requirement. We demonstrate a proof of concept for a multiplexing assay that pools PCR amplicons for MinION sequencing to enable sequencing of multiple templates from multiple individuals, which could be applied to gene-targeted diagnostics. A combined strategy of barcoding and sample pooling was developed for simultaneous multiplex MinION sequencing of 100 PCR amplicons. The amplicons are family-specific, spanning a total of 30 loci in DNA isolated from 82 human neurodevelopmental cases and family members. The target regions were chosen for further interrogation because a potentially disease-causative variant had been identified in affected individuals following Illumina exome sequencing. The pooled MinION sequences were deconvoluted by aligning to custom references using the minimap2 aligner software. Our multiplexing approach produced an interpretable and expected sequence from 29 of the 30 targeted genetic loci. The sequence variant which was not correctly resolved in the MinION sequence was adjacent to a five nucleotide homopolymer. It is already known that homopolymers present a resolution problem with the MinION approach. Interestingly despite equimolar quantities of PCR amplicon pooled for sequencing, significant variation in the depth of coverage (127×-19,626×; mean = 8321×, std err = 452.99) was observed. We observed independent relationships between depth of coverage and target length, and depth of coverage and GC content. These relationships demonstrate biases of the MinION sequencer for longer templates and those with lower GC content. We demonstrate an efficient approach for variant discovery or confirmation from short DNA templates using the MinION sequencing device. With less than 130 × depth of coverage required for accurate genotyping, the methodology described here allows for rapid highly multiplexed targeted sequencing of large numbers of samples in a minimally equipped laboratory with a potential cost as much 200 × less than that from Sanger sequencing.

Screening for phenotypic outliers identifies an unusually low concentration of a ß-lactoglobulin B protein isoform in bovine milk caused by a synonymous SNP.

BACKGROUND: Milk samples from 10,641 dairy cattle were screened by a mass spectrometry method for extreme concentrations of the A or B isoforms of the whey protein, ß-lactoglobulin (BLG), to identify causative genetic variation driving changes in BLG concentration. RESULTS: A cohort of cows, from a single sire family, was identified that produced milk containing a low concentration of the BLG B protein isoform. A genome-wide association study (GWAS) of BLG B protein isoform concentration in milk from AB heterozygous cows, detected a group of highly significant single nucleotide polymorphisms (SNPs) within or close to the BLG gene. Among these was a synonymous G/A variation at position + 78 bp in exon 1 of the BLG gene (chr11:103256256G > A). The effect of the A allele of this SNP (which we named B') on BLG expression was evaluated in a luciferase reporter assay in transfected CHO-K1 and MCF-7 cells. In both cell types, the presence of the B' allele in a plasmid containing the bovine BLG gene from -922 to + 898 bp (relative to the transcription initiation site) resulted in a 60% relative reduction in mRNA expression, compared to the plasmid containing the wild-type B sequence allele. Examination of a mammary RNAseq dataset (n = 391) identified 14 heterozygous carriers of the B' allele which were homozygous for the BLG B protein isoform (BB'). The level of expression of the BLG B' allele was 41.9 ± 1.0% of that of the wild-type BLG B allele. Milk samples from three cows, homozygous for the A allele at chr11:103,256,256 (B'B'), were analysed (HPLC) and showed BLG concentrations of 1.04, 1.26 and 1.83 g/L relative to a mean of 4.84 g/L in milk from 16 herd contemporaries of mixed (A and B) BLG genotypes. The mechanism by which B' downregulates milk BLG concentration remains to be determined. CONCLUSIONS: High-throughput screening and identification of outliers, enabled the discovery of a synonymous G > A mutation in exon 1 of the B allele of the BLG gene (B'), which reduced the milk concentration of ß-lactoglobulin B protein isoform, by more than 50%. Milk from cows carrying the B' allele is expected to have improved processing characteristics, particularly for cheese-making.

A Capra hircus chromosome 19 locus linked to milk production influences mammary conformation.

BACKGROUND: Economically important milk production traits including milk volume, milk fat and protein yield vary considerably across dairy goats in New Zealand. A significant portion of the variation is attributable to genetic variation. Discovery of genetic markers linked to milk production traits can be utilised to drive selection of high-performance animals. A previously reported genome wide association study across dairy goats in New Zealand identified a quantitative trait locus (QTL) located on chromosome 19. The most significantly associated single nucleotide polymorphism (SNP) marker for this locus is located at position 26,610,610 (SNP marker rs268292132). This locus is associated with multiple milk production traits including fat, protein and volume. The predicted effect of selection for the beneficial haplotype would result in an average production increase of 2.2 kg fat, 1.9 kg protein and 73.6 kg milk yield. An outstanding question was whether selection for the beneficial allele would co-select for any negative pleiotropic effects. An adverse relationship between milk production and udder health traits has been reported at this locus. Therefore, a genome wide association study was undertaken looking for loci associated with udder traits. RESULTS: The QTL and production associated marker rs268292132 was identified in this study to also be associated with several goat udder traits including udder depth (UD), fore udder attachment (FUA) and rear udder attachment (RUA). Our study replicates the negative relationship between production and udder traits with the high production allele at position 19:26,610,610 (SNP marker rs268292132) associated with an adverse change in UD, FUA and RUA. CONCLUSIONS: Our study has confirmed the negative relationship between udder traits and production traits in the NZ goat population. We have found that the frequency of the high production allele is relatively high in the NZ goat population, indicating that its effect on udder conformation is not significantly detrimental on animal health. It will however be important to monitor udder conformation as the chromosome 19 locus is progressively implemented for marker assisted selection. It will also be of interest to determine if the gene underlying the production QTL has a direct effect on mammary gland morphology or whether the changes observed are a consequence of the increased milk volume.

Novel PRMT7 mutation in a rare case of dysmorphism and intellectual disability.

Protein arginine N-methyltransferase 7 (PRMT7) encodes an arginine methyltransferase central to a number of fundamental biological processes, mutations in which result in an autosomal recessive developmental disorder characterized by short stature, brachydactyly, intellectual developmental disability and seizures (SBIDDS). To date, fewer than 15 patients with biallelic mutations in PRMT7 have been documented. Here we report brothers from a consanguineous Iraqi family presenting with a developmental disorder characterized by global developmental delay, shortened stature, facial dysmorphisms, brachydactyly, and kidney dysfunction. In both affected brothers, whole genome sequencing (WGS) identified a novel homozygous substitution in PRMT7 (ENST00000339507.5), c.1097 G > A (p.Cys366Tyr), considered to account for the majority of the phenotypic presentation. Rare compound heterozygous mutations in the dysplasia-associated perlecan-encoding HSPG2 gene (ENST00000374695.3) were also found (c.10721-2dupA, p.Ser71Asn and c.212 G > A), potentially accounting for the kidney dysfunction. In addition to expanding the known mutational spectrum of variably expressive PRMT7 mutations alongside potential digenic inheritance with HSPG2, this report underlines the diagnostic utility of a WGS-guided analysis in the detection of rare genetic disorders.

A Multi-Omic Huntington's Disease Transgenic Sheep-Model Database for Investigating Disease Pathogenesis.

BACKGROUND: The pathological mechanism of cellular dysfunction and death in Huntington's disease (HD) is not well defined. Our transgenic HD sheep model (OVT73) was generated to investigate these mechanisms and for therapeutic testing. One particular cohort of animals has undergone focused investigation resulting in a large interrelated multi-omic dataset, with statistically significant changes observed comparing OVT73 and control 'omic' profiles and reported in literature. OBJECTIVE: Here we make this dataset publicly available for the advancement of HD pathogenic mechanism discovery. METHODS: To enable investigation in a user-friendly format, we integrated seven multi-omic datasets from a cohort of 5-year-old OVT73 (n = 6) and control (n = 6) sheep into a single database utilising the programming language R. It includes high-throughput transcriptomic, metabolomic and proteomic data from blood, brain, and other tissues. RESULTS: We present the 'multi-omic' HD sheep database as a queryable web-based platform that can be used by the wider HD research community (https://hdsheep.cer.auckland.ac.nz/). The database is supported with a suite of simple automated statistical analysis functions for rapid exploratory analyses. We present examples of its use that validates the integrity relative to results previously reported. The data may also be downloaded for user determined analysis. CONCLUSION: We propose the use of this online database as a hypothesis generator and method to confirm/refute findings made from patient samples and alternate model systems, to expand our understanding of HD pathogenesis. Importantly, additional tissue samples are available for further investigation of this cohort.

Castration delays epigenetic aging and feminizes DNA methylation at androgen-regulated loci.

In mammals, females generally live longer than males. Nevertheless, the mechanisms underpinning sex-dependent longevity are currently unclear. Epigenetic clocks are powerful biological biomarkers capable of precisely estimating chronological age and identifying novel factors influencing the aging rate using only DNA methylation data. In this study, we developed the first epigenetic clock for domesticated sheep (Ovis aries), which can predict chronological age with a median absolute error of 5.1 months. We have discovered that castrated male sheep have a decelerated aging rate compared to intact males, mediated at least in part by the removal of androgens. Furthermore, we identified several androgen-sensitive CpG dinucleotides that become progressively hypomethylated with age in intact males, but remain stable in castrated males and females. Comparable sex-specific methylation differences in MKLN1 also exist in bat skin and a range of mouse tissues that have high androgen receptor expression, indicating that it may drive androgen-dependent hypomethylation in divergent mammalian species. In characterizing these sites, we identify biologically plausible mechanisms explaining how androgens drive male-accelerated aging.

Testes of DAZL null neonatal sheep lack prospermatogonia but maintain normal somatic cell morphology and marker expression.

Multiplying the germline would increase the number of offspring that can be produced from selected animals, accelerating genetic improvement for livestock breeding. This could be achieved by producing multiple chimaeric animals, each carrying a mix of donor and host germ cells in their gonads. However, such chimaeric germlines would produce offspring from both donor and host genotypes, limiting the rate of genetic improvement. To resolve this problem, we disrupted the RNA-binding protein DAZL and generated germ cell-deficient host animals. Using Cas9-mediated homology-directed repair (HDR), we introduced a DAZL loss-of-function mutation in male ovine fetal fibroblasts. Following manual single cell isolation, 4/48 (8.3%) of donor cell strains were homozygously HDR-edited. Sequence-validated strains were used as nuclear donors for somatic cell cloning to generate three lambs, which died at birth. All DAZL null male neonatal sheep lacked germ cells on histological sections and showed greatly reduced germ cell markers. Somatic cells within their testes were morphologically intact and expressed normal levels of lineage-specific markers, suggesting that the germ cell niche remained intact. This extends the DAZL mutant phenotype beyond mice into agriculturally relevant ruminants, providing a pathway for using absolute germline transmitters in rapid livestock improvement.