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.

SARS-CoV-2 evolution has increased resistance to monoclonal antibodies and first-generation COVID-19 vaccines: Is there a future therapeutic role for soluble ACE2 receptors for COVID-19?

COVID-19 has caused calamitous health, economic and societal consequences. Although several COVID-19 vaccines have received full authorization for use, global deployment has faced political, financial and logistical challenges. The efficacy of first-generation COVID-19 vaccines is waning and breakthrough infections are allowing ongoing transmission and evolution of SARS-CoV-2. Furthermore, COVID-19 vaccine efficacy relies on a functional immune system. Despite receiving three primary doses and three or more heterologous boosters, some immunocompromised patients may not be adequately protected by COVID-19 vaccines and remain vulnerable to severe disease. The evolution of new SARS-CoV-2 variants has also resulted in the rapid obsolescence of monoclonal antibodies. Convalescent plasma from COVID-19 survivors has produced inconsistent results. New drugs such as Paxlovid (nirmatrelvir/ritonavir) are beyond the reach of low- and middle-income countries. With widespread use of Paxlovid, it is likely nirmatrelvir-resistant clades of SARS-CoV-2 will emerge in the future. There is thus an urgent need for new effective anti-SARS-CoV-2 treatments. The in vitro efficacy of soluble ACE2 against multiple SARS-CoV-2 variants including omicron (B.1.1.529), was recently described using a competitive ELISA assay as a surrogate marker for virus neutralization. This indicates soluble wild-type ACE2 receptors are likely to be resistant to viral evolution. Nasal and inhaled treatment with soluble ACE2 receptors has abrogated severe disease in animal models of COVID-19. There is an urgent need for clinical trials of this new class of antiviral therapeutics, which could complement vaccines and Paxlovid.

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.

Challenges for gene editing in common variable immunodeficiency disorders: Current and future prospects.

The original CRISPR Cas9 gene editing system and subsequent innovations offers unprecedented opportunities to correct severe genetic defects including those causing Primary Immunodeficiencies (PIDs). Common Variable Immunodeficiency Disorders (CVID) are the most frequent symptomatic PID in adults and children. Unlike many other PIDs, patients meeting CVID criteria do not have a definable genetic defect and cannot be considered to have an inborn error of immunity (IEI). Patients with a CVID phenotype carrying a causative mutation are deemed to have a CVID-like disorder consequent to an IEI. Patients from consanguineous families often have highly penetrant early-onset autosomal recessive forms of CVID-like disorders. Individuals from non-consanguineous families may have autosomal dominant CVID-like disorders with variable penetrance and expressivity. This essay explores the potential clinical utility as well as the current limitations and risks of gene editing including collateral genotoxicity. In the immediate future the main application of this technology is likely to be the in vitro investigation of epigenetic and polygenic mechanisms, which are likely to underlie many cases of CVID and CVID-like disorders. In the longer-term, the CRISPR Cas9 system and other gene-based therapies could be utilized to treat CVID-like disorders, where the underlying IEI is known.

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.

Soluble wild-type ACE2 molecules inhibit newer SARS-CoV-2 variants and are a potential antiviral strategy to mitigate disease severity in COVID-19.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus responsible for coronavirus disease of 2019 (COVID-19), has caused havoc around the world. While several COVID-19 vaccines and drugs have been authorized for use, these antiviral drugs remain beyond the reach of most low- and middle-income countries. Rapid viral evolution is reducing the efficacy of vaccines and monoclonal antibodies and contributing to the deaths of some fully vaccinated persons. Others with normal immunity may have chosen not to be vaccinated and remain at risk if they contract the infection. Vaccines may not protect some immunodeficient patients from SARS-CoV-2, who are also at increased risk of chronic COVID-19 infection, a dangerous stalemate between the virus and a suboptimal immune response. Intra-host viral evolution could rapidly lead to the selection and dominance of vaccine and monoclonal antibody-resistant clades of SARS-CoV-2. There is thus an urgent need to develop new treatments for COVID-19. The NZACE2-Patari project, comprising modified soluble angiotensin-converting enzyme 2 (ACE2) molecules, seeks to intercept and block SARS-CoV-2 infection of the respiratory mucosa. In vitro data presented here show that soluble wild-type ACE2 molecules retain the ability to effectively block the Spike (S) glycoprotein of SARS-CoV-2 variants including the ancestral Wuhan, delta (B.1.617.2) and omicron (B.1.1.529) strains. This therapeutic strategy may prove effective if implemented early during the nasal phase of the infection and may act synergistically with other antiviral drugs such as Paxlovid to further mitigate disease severity.

Successful editing and maintenance of lactogenic gene expression in primary bovine mammary epithelial cells.

In vitro investigation of bovine lactation processes is limited by a lack of physiologically representative cell models. This deficiency is most evident through the minimal or absent expression of lactation-specific genes in cultured bovine mammary tissues. Primary bovine mammary epithelial cells (pbMECs) extracted from lactating mammary tissue and grown in culture initially express milk protein transcripts at relatively representative levels. However, expression drops dramatically after only three or four passages, which greatly reduces the utility of primary cells to model and further examine lactogenesis. To investigate the effects of alternate alleles in pbMECs including effects on transcription, we have developed methods to deliver CRISPR-Cas9 gene editing reagents to primary mammary cells, resulting in very high editing efficiencies. We have also found that culturing the cells on an imitation basement membrane composed of Matrigel, results in the restoration of a more representative lactogenic gene expression profile and the cells forming three-dimensional structures in vitro. Here, we present data from four pbMEC lines recovered from pregnant cows and detail the expression profile of five key milk synthesis genes in these MECs grown on Matrigel. Additionally, we describe an optimised method for preferentially selecting CRISPR-Cas9-edited cells conferring a knock-out of DGAT1, using fluorescence-activated cell sorting (FACS). The combination of these techniques facilitates the use of pbMECs as a model to investigate the effects of gene introgressions and genetic variation in lactating mammary tissue.

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.

Genetic control of serum 25(OH)D levels and its association with ethnicity.

BACKGROUND: Identified DNA variants associated with serum 25-hydroxyvitamin vitamin D (25[OH]D) concentration may provide mechanistic insights into the vitamin D metabolic pathway in individuals. Our aim was to further characterise participants and their serum 25(OH)D concentration at baseline using candidate single nucleotide polymorphism (SNP) genotyping. METHODS: 5110 participants, aged 50-84 years, were recruited from the community. Blood samples were collected at baseline to measure serum 25(OH)D by liquid chromatography mass spectrometry and the participants were genotyped for four markers close to or within genes in the vitamin D metabolic pathway known to be associated with differences in 25(OH)D. The markers and their associated genes were rs12785878 (DHCR7), rs10741657 (CYP2R1), rs4588 (DBP) and rs2228570 (VDR). RESULTS: All four markers had significantly different genotype distributions and minor allele frequencies between the four self-determined ethnicities (European/Other, Maori, Pacific, and South Asian). For example, the frequency in each ethnic group of the G allele for the marker rs12785878 was 0.26, 0.71, 0.89, and 0.78 respectively. Using multivariable regression in the full cohort, three out of four markers were significantly associated with baseline concentrations of 25(OH)D (mean differences: 2.9-10.9 nmol/L). Collectively, the four markers explained 8.4% of the variation in 25(OH)D concentrations. CONCLUSION: Significant ethnic variations exist in the distribution of alleles associated with serum 25(OH)D concentration, particularly rs12785878, in a multi-ethnic community sample from New Zealand.

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.