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.

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.

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.

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.

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.

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.

Universal DNA methylation age across mammalian tissues

Aging, often considered a result of random cellular damage, can be accurately estimated using DNA methylation profiles, the foundation of pan-tissue epigenetic clocks. Here, we demonstrate the development of universal pan-mammalian clocks, using 11,754 methylation arrays from our Mammalian Methylation Consortium, which encompass 59 tissue types across 185 mammalian species. These predictive models estimate mammalian tissue age with high accuracy (r > 0.96). Age deviations correlate with human mortality risk, mouse somatotropic axis mutations and caloric restriction. We identified specific cytosines with methylation levels that change with age across numerous species. These sites, highly enriched in polycomb repressive complex 2-binding locations, are near genes implicated in mammalian development, cancer, obesity and longevity. Our findings offer new evidence suggesting that aging is evolutionarily conserved and intertwined with developmental processes across all mammals.

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.

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.

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.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00711-z.

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

Nil.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00667-y.

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.

Contaminating DNA in human saliva alters the detection of variants from whole genome sequencing

Cells obtained from human saliva are commonly used as an alternative DNA source when blood is difficult or less convenient to collect. Although DNA extracted from saliva is considered to be of comparable quality to that derived from blood, recent studies have shown that non-human contaminating DNA derived from saliva can confound whole genome sequencing results. The most concerning complication is that non-human reads align to the human reference genome using standard methodology, which can critically affect the resulting variant genotypes identified in a genome. We identified clusters of anomalous variants in saliva DNA derived reads which aligned in an atypical manner. These reads had only short regions of identity to the human reference sequence, flanked by soft clipped sequence. Sequence comparisons of atypically aligning reads from eight human saliva-derived samples to RefSeq genomes revealed the majority to be of bacterial origin (63.46%). To partition the non-human reads during the alignment step, a decoy of the most prevalent bacterial genome sequences was designed and utilised. This reduced the number of atypically aligning reads when trialled on the eight saliva-derived samples by 44% and most importantly prevented the associated anomalous genotype calls. Saliva derived DNA is often contaminated by DNA from other species. This can lead to non-human reads aligning to the human reference genome using current alignment best-practices, impacting variant identification. This problem can be diminished by using a bacterial decoy in the alignment process.

Inconsistencies in histone acetylation patterns among different HD model systems and HD post-mortem brains

Huntington's disease (HD) is a neurodegenerative disorder caused by a CAG trinucleotide repeat expansion in exon 1 of the huntingtin gene. Emerging evidence shows that additional epigenetic factors can modify disease phenotypes. Harnessing the ability of the epigenome to modify the disease for therapeutic purposes is therefore of interest. Epigenome modifiers, such as histone deacetylase inhibitors (HDACi), have improved pathology in a range of HD models. Yet in clinical trials, HDACi have failed to alleviate HD symptoms in patients. This study investigated potential reasons for the lack of translation of the therapeutic benefits of HDACi from lab to clinic. We analysed histone acetylation patterns of immuno-positive nuclei from brain sections and tissue microarrays from post-mortem human control and HD cases alongside several well-established HD models (OVT73 transgenic HD sheep, YAC128 mice, and an in vitro cell model expressing 97Q mutant huntingtin). Significant increases in histone H4 acetylation were observed in post-mortem HD cases, OVT73 transgenic HD sheep and in vitro models; these changes were absent in YAC128 mice. In addition, nuclear labelling for acetyl-histone H4 levels were inversely proportional to mutant huntingtin aggregate load in HD human cortex. Our data raise concerns regarding the utility of HDACi for the treatment of HD when regions of pathology exhibit already elevated histone acetylation patterns and emphasize the importance of searching for alternative epigenetic targets in future therapeutic strategies aiming to rescue HD phenotypes.

A new mechanism for a familiar mutation - bovine DGAT1 K232A modulates gene expression through multi-junction exon splice enhancement

BACKGROUND

The DGAT1 gene encodes an enzyme responsible for catalysing the terminal reaction in mammary triglyceride synthesis, and underpins a well-known pleiotropic quantitative trait locus (QTL) with a large influence on milk composition phenotypes. Since first described over 15 years ago, a protein-coding variant K232A has been assumed as the causative variant underlying these effects, following in-vitro studies that demonstrated differing levels of triglyceride synthesis between the two protein isoforms.

RESULTS

We used a large RNAseq dataset to re-examine the underlying mechanisms of this large milk production QTL, and hereby report novel expression-based functions of the chr14 g.1802265AA > GC variant that encodes the DGAT1 K232A substitution. Using expression QTL (eQTL) mapping, we demonstrate a highly-significant mammary eQTL for DGAT1, where the K232A mutation appears as one of the top associated variants for this effect. By conducting in vitro expression and splicing experiments in bovine mammary cell culture, we further show modulation of splicing efficiency by this mutation, likely through disruption of an exon splice enhancer as a consequence of the allele encoding the 232A variant.

CONCLUSIONS

The relative contributions of the enzymatic and transcription-based mechanisms now attributed to K232A remain unclear; however, these results suggest that transcriptional impacts contribute to the diversity of lactation effects observed at the DGAT1 locus.

Genetic correlations of milk fatty acid contents predicted from milk mid-infrared spectra in New Zealand dairy cattle

The objective of this study was to estimate genetic correlations among milk fatty acid (FA) concentrations in New Zealand dairy cattle. Concentrations of each of the most common FA, expressed as a percentage of the total FA, were determined by gas chromatography on a specific cohort of animals. Using this data set, prediction equations were derived using mid-infrared (MIR) spectroscopy data collected from the same samples. These prediction equations were applied to a large data set of MIR measurements in 34,141 milk samples from 3,445 Holstein-Friesian, 2,935 Jersey, and 3,609 crossbred Holstein-Friesian × Jersey cows, sampled an average of 3.42 times during the 2007-2008 season. Data were analyzed using univariate and bivariate repeatability animal models. Heritability of predicted FA concentration in milk fat ranged from 0.21 to 0.42, indicating that genetic selection could be used to change the FA composition of milk. The de novo synthesized FA (C6:0, C8:0, C10:0, C12:0, and C14:0) showed strong positive genetic correlations with each other, ranging from 0.24 to 0.99. Saturated FA were negatively correlated with unsaturated (-0.93) and polyunsaturated (-0.84) FA. The saturated FA were positively correlated with milk fat yield and fat percentage, whereas the unsaturated FA were negatively associated with fat yield and fat percentage. Our results indicate that bovine milk FA composition can be changed through genetic selection using MIR as a phenotypic proxy.

RBV: Read balance validator, a tool for prioritising copy number variations in germline conditions

The popularisation and decreased cost of genome resequencing has resulted in an increased use in molecular diagnostics. While there are a number of established and high quality bioinfomatic tools for identifying small genetic variants including single nucleotide variants and indels, currently there is no established standard for the detection of copy number variants (CNVs) from sequence data. The requirement for CNV detection from high throughput sequencing has resulted in the development of a large number of software packages. These tools typically utilise the sequence data characteristics: read depth, split reads, read pairs, and assembly-based techniques. However, the additional source of information from read balance (defined as relative proportion of reads of each allele at each position) has been underutilised in the existing applications. Here we present Read Balance Validator (RBV), a bioinformatic tool that uses read balance for prioritisation and validation of putative CNVs. The software simultaneously interrogates nominated regions for the presence of deletions or multiplications, and can differentiate larger CNVs from diploid regions. Additionally, the utility of RBV to test for inheritance of CNVs is demonstrated in this report. RBV is a CNV validation and prioritisation bioinformatic tool for both genome and exome sequencing available as a python package from https://github.com/whitneywhitford/RBV.

Genome-wide association analysis reveals QTL and candidate mutations involved in white spotting in cattle

BACKGROUND

White spotting of the coat is a characteristic trait of various domestic species including cattle and other mammals. It is a hallmark of Holstein-Friesian cattle, and several previous studies have detected genetic loci with major effects for white spotting in animals with Holstein-Friesian ancestry. Here, our aim was to better understand the underlying genetic and molecular mechanisms of white spotting, by conducting the largest mapping study for this trait in cattle, to date.

RESULTS

Using imputed whole-genome sequence data, we conducted a genome-wide association analysis in 2973 mixed-breed cows and bulls. Highly significant quantitative trait loci (QTL) were found on chromosomes 6 and 22, highlighting the well-established coat color genes KIT and MITF as likely responsible for these effects. These results are in broad agreement with previous studies, although we also report a third significant QTL on chromosome 2 that appears to be novel. This signal maps immediately adjacent to the PAX3 gene, which encodes a known transcription factor that controls MITF expression and is the causal locus for white spotting in horses. More detailed examination of these loci revealed a candidate causal mutation in PAX3 (p.Thr424Met), and another candidate mutation (rs209784468) within a conserved element in intron 2 of MITF transcripts expressed in the skin. These analyses also revealed a mechanistic ambiguity at the chromosome 6 locus, where highly dispersed association signals suggested multiple or multiallelic QTL involving KIT and/or other genes in this region.

CONCLUSIONS

Our findings extend those of previous studies that reported KIT as a likely causal gene for white spotting, and report novel associations between candidate causal mutations in both the MITF and PAX3 genes. The sizes of the effects of these QTL are substantial, and could be used to select animals with darker, or conversely whiter, coats depending on the desired characteristics.

Extensive intraspecies cryptic variation in an ancient embryonic gene regulatory network

Innovations in metazoan development arise from evolutionary modification of gene regulatory networks (GRNs). We report widespread cryptic variation in the requirement for two key regulatory inputs, SKN-1/Nrf2 and MOM-2/Wnt, into the C. elegans endoderm GRN. While some natural isolates show a nearly absolute requirement for these two regulators, in others, most embryos differentiate endoderm in their absence. GWAS and analysis of recombinant inbred lines reveal multiple genetic regions underlying this broad phenotypic variation. We observe a reciprocal trend, in which genomic variants, or knockdown of endoderm regulatory genes, that result in a high SKN-1 requirement often show low MOM-2/Wnt requirement and vice-versa, suggesting that cryptic variation in the endoderm GRN may be tuned by opposing requirements for these two key regulatory inputs. These findings reveal that while the downstream components in the endoderm GRN are common across metazoan phylogeny, initiating regulatory inputs are remarkably plastic even within a single species.

Cerebral Vitamin B5 (D-Pantothenic Acid) Deficiency as a Potential Cause of Metabolic Perturbation and Neurodegeneration in Huntington's Disease

Huntington's disease (HD) is a neurodegenerative disorder caused by an expanded CAG repeat in exon 1 of the HTT gene. HD usually manifests in mid-life with loss of GABAergic projection neurons from the striatum accompanied by progressive atrophy of the putamen followed by other brain regions, but linkages between the genetics and neurodegeneration are not understood. We measured metabolic perturbations in HD-human brain in a case-control study, identifying pervasive lowering of vitamin B5, the obligatory precursor of coenzyme A (CoA) that is essential for normal intermediary metabolism. Cerebral pantothenate deficiency is a newly-identified metabolic defect in human HD that could potentially: (i) impair neuronal CoA biosynthesis; (ii) stimulate polyol-pathway activity; (iii) impair glycolysis and tricarboxylic acid cycle activity; and (iv) modify brain-urea metabolism. Pantothenate deficiency could lead to neurodegeneration/dementia in HD that might be preventable by treatment with vitamin B5.

Evaluation of the performance of copy number variant prediction tools for the detection of deletions from whole genome sequencing data

BACKGROUND

Whole genome sequencing (WGS) has increased in popularity and decreased in cost over the past decade, rendering this approach as a viable and sensitive method for variant detection. In addition to its utility for single nucleotide variant detection, WGS data has the potential to detect Copy Number Variants (CNV) to fine resolution. Many CNV detection software packages have been developed exploiting four main types of data: read pair, split read, read depth, and assembly based methods. The aim of this study was to evaluate the efficiency of each of these main approaches in detecting germline deletions.

METHODS

WGS data and high confidence deletion calls for the individual NA12878 from the Genome in a Bottle consortium were the benchmark dataset. The performance of BreakDancer, CNVnator, Delly, FermiKit, and Pindel was assessed by comparing the accuracy and sensitivity of each software package in detecting deletions exceeding 1 kb.

RESULTS

There was considerable variability in the outputs of the different WGS CNV detection programs. The best performance was seen from BreakDancer and Delly, with 92.6% and 96.7% sensitivity, respectively and 34.5% and 68.5% false discovery rate (FDR), respectively. In comparison, Pindel, CNVnator, and FermiKit were less effective with sensitivities of 69.1%, 66.0%, and 15.8%, respectively and FDR of 91.3%, 69.0%, and 31.7%, respectively. Concordance across software packages was poor, with only 27 of the total 612 benchmark deletions identified by all five methodologies.

CONCLUSIONS

The WGS based CNV detection tools evaluated show disparate performance in identifying deletions ≥1 kb, particularly those utilising different input data characteristics. Software that exploits read pair based data had the highest sensitivity, namely BreakDancer and Delly. BreakDancer also had the second lowest false discovery rate. Therefore, in this analysis read pair methods (BreakDancer in particular) were the best performing approaches for the identification of deletions ≥1 kb, balancing accuracy and sensitivity. There is potential for improvement in the detection algorithms, particularly for reducing FDR. This analysis has validated the utility of WGS based CNV detection software to reliably identify deletions, and these findings will be of use when choosing appropriate software for deletion detection, in both research and diagnostic medicine.

Widespread cis-regulation of RNA editing in a large mammal

Post-transcriptional RNA editing may regulate transcript expression and diversity in cells, with potential impacts on various aspects of physiology and environmental adaptation. A small number of recent genome-wide studies in Drosophila, mouse, and human have shown that RNA editing can be genetically modulated, highlighting loci that quantitatively impact editing of transcripts. The potential gene expression and physiological consequences of these RNA-editing quantitative trait loci (edQTL), however, are almost entirely unknown. Here, we present analyses of RNA editing in a large domestic mammal (Bos taurus), where we use whole-genome and high-depth RNA sequencing to discover, characterize, and conduct genetic mapping studies of novel transcript edits. Using a discovery population of nine deeply sequenced cows, we identify 2413 edit sites in the mammary transcriptome, the majority of which are adenosine to inosine edits (98.6%). Most sites are predicted to reside in double-stranded secondary structures (85.1%), and quantification of the rates of editing in an additional 355 cows reveals editing is negatively correlated with gene expression in the majority of cases. Genetic analyses of RNA editing and gene expression highlight 152 cis-regulated edQTL, of which 15 appear to cosegregate with expression QTL effects. Trait association analyses in a separate population of 9989 lactating cows also shows 12 of the cis-edQTL coincide with at least one cosegregating lactation QTL. Together, these results enhance our understanding of RNA-editing dynamics in mammals, and suggest mechanistic links by which loci may impact phenotype through RNA editing mediated processes.

Multiple QTL underlie milk phenotypes at the CSF2RB locus

Background

Over many years, artificial selection has substantially improved milk production by cows. However, the genes that underlie milk production quantitative trait loci (QTL) remain relatively poorly characterised. Here, we investigate a previously reported QTL located at the CSF2RB locus on chromosome 5, for several milk production phenotypes, to better understand its underlying genetic and molecular causes.

Results

Using a population of 29,350 taurine dairy cows, we conducted association analyses for milk yield and composition traits, and identified highly significant QTL for milk yield, milk fat concentration, and milk protein concentration. Strikingly, protein concentration and milk yield appear to show co-located yet genetically distinct QTL. To attempt to understand the molecular mechanisms that might be mediating these effects, gene expression data were used to investigate eQTL for 11 genes in the broader interval. This analysis highlighted genetic impacts on CSF2RB and NCF4 expression that share similar association signatures to those observed for lactation QTL, strongly implicating one or both of these genes as responsible for these effects. Using the same gene expression dataset representing 357 lactating cows, we also identified 38 novel RNA editing sites in the 3′ UTR of CSF2RB transcripts. The extent to which two of these sites were edited also appears to be genetically co-regulated with lactation QTL, highlighting a further layer of regulatory complexity that involves the CSF2RB gene.

Conclusions

This locus presents a diversity of molecular and lactation QTL, likely representing multiple overlapping effects that, at a minimum, highlight the CSF2RB gene as having a causal role in these processes.

Electronic supplementary material

The online version of this article (10.1186/s12711-019-0446-x) contains supplementary material, which is available to authorized users.

Modelling brain dopamine-serotonin vesicular transport disease in Caenorhabditis elegans

Brain dopamine-serotonin vesicular transport disease is a rare disease caused by autosomal recessive mutations in the SLC18A2 gene, which encodes the VMAT2 protein. VMAT2 is a membrane protein responsible for vesicular transport of monoamines, and its disruption negatively affects neurotransmission. This results in a severe neurodevelopmental disorder affecting motor skills and development, and causes muscular hypotonia. The condition was initially described in a consanguineous Saudi Arabian family with affected siblings homozygous for a P387L mutation. We subsequently found a second mutation in a New Zealand family (homozygous P237H), which was later also identified in an Iraqi family. Pramipexole has been shown to have some therapeutic benefit. Transgenic Caenorhabditis elegans were developed to model the P237H and P387L mutations. Investigations into dopamine- and serotonin-related C. elegans phenotypes, including pharyngeal pumping and grazing, showed that both mutations cause significant impairment of these processes when compared with a non-transgenic N2 strain and a transgenic containing the wild-type human SLC18A2 gene. Preliminary experiments investigating the therapeutic effects of serotonin and pramipexole demonstrated that serotonin could successfully restore the pharyngeal pumping phenotype. These analyses provide further support for the role of these mutations in this disease.

Summary: The first Caenorhabditis elegans model to study brain dopamine-serotonin vesicular transport disease, demonstrating impairment of pharyngeal pumping and grazing phenotypes.

De Novo Pathogenic Variants in CACNA1E Cause Developmental and Epileptic Encephalopathy with Contractures, Macrocephaly, and Dyskinesias

Developmental and epileptic encephalopathies (DEEs) are severe neurodevelopmental disorders often beginning in infancy or early childhood that are characterized by intractable seizures, abundant epileptiform activity on EEG, and developmental impairment or regression. CACNA1E is highly expressed in the central nervous system and encodes the α₁-subunit of the voltage-gated Ca_V2.3 channel, which conducts high voltage-activated R-type calcium currents that initiate synaptic transmission. Using next-generation sequencing techniques, we identified de novo CACNA1E variants in 30 individuals with DEE, characterized by refractory infantile-onset seizures, severe hypotonia, and profound developmental impairment, often with congenital contractures, macrocephaly, hyperkinetic movement disorders, and early death. Most of the 14, partially recurring, variants cluster within the cytoplasmic ends of all four S6 segments, which form the presumed Ca_V2.3 channel activation gate. Functional analysis of several S6 variants revealed consistent gain-of-function effects comprising facilitated voltage-dependent activation and slowed inactivation. Another variant located in the domain II S4-S5 linker results in facilitated activation and increased current density. Five participants achieved seizure freedom on the anti-epileptic drug topiramate, which blocks R-type calcium channels. We establish pathogenic variants in CACNA1E as a cause of DEEs and suggest facilitated R-type calcium currents as a disease mechanism for human epilepsy and developmental disorders.

Penetrance and expressivity of the R858H CACNA1C variant in a five-generation pedigree segregating an arrhythmogenic channelopathy

BACKGROUND

Isolated cardiac arrhythmia due to a variant in CACNA1C is of recent knowledge. Most reports have been of singleton cases or of quite small families, and estimates of penetrance and expressivity have been difficult to obtain. We here describe a large pedigree, from which such estimates have been calculated.

METHODS

We studied a five-generation family, in which a CACNA1C variant c.2573G>A p.Arg858His co-segregates with syncope and cardiac arrest, documenting electrocardiographic data and cardiac symptomatology. The reported patients/families from the literature with CACNA1C gene variants were reviewed, and genotype-phenotype correlations are drawn.

RESULTS

The range of phenotype in the studied family is wide, from no apparent effect, through an asymptomatic QT interval prolongation on electrocardiography, to episodes of presyncope and syncope, ventricular fibrillation, and sudden death. QT prolongation showed inconsistent correlation with functional cardiology. Based upon analysis of 28 heterozygous family members, estimates of penetrance and expressivity are derived.

CONCLUSIONS

These estimates of penetrance and expressivity, for this specific variant, may be useful in clinical practice. Review of the literature indicates that individual CACNA1C variants have their own particular genotype-phenotype correlations. We suggest that, at least in respect of the particular variant reported here, "arrhythmogenic channelopathy" may be a more fitting nomenclature than long QT syndrome.

Genomic medicine must reduce, not compound, health inequities: the case for hauora-enhancing genomic resources for New Zealand

Precision medicine seeks to draw on data from both individuals and populations across disparate domains to influence and support diagnosis, management and prevention in healthcare at the level of the individual patient and their family/whānau. Central to this initiative is incorporating the effects of the inherent variation that lies within genomes and can influence health outcomes. Identifying and interpreting such variation requires an accurate, valid and representative dataset to firstly define what variants are present and then assess the potential relevance for the health of a person, their family/whānau and the wider community to which they belong. Globally the variation embedded within genomes differs enormously and has been shaped by the size, constitution, historical origins and evolutionary history of their source populations. Māori, and more broadly Pacific peoples, differ substantially in terms of genomic variation compared to the more closely studied European and Asian populations. In the absence of accurate genomic information from Māori and Pacific populations, the precise interpretation of genomic data and the success and benefits of genomic medicine will be disproportionately less for those Māori and Pacific peoples. In this viewpoint article we, as a group of healthcare professionals, researchers and scientists, present a case for assembling genomic resources that catalogue the characteristics of the genomes of New Zealanders, with an emphasis on peoples of Māori and Polynesian ancestry, as a healthcare imperative. In proposing the creation of these resources, we note that their governance and management must be led by iwi and Māori and Pacific representatives. Assembling a genomic resource must be informed by cultural concepts and values most especially understanding that, at a physical and spiritual level, whakapapa is embodied within the DNA of a person. Therefore DNA and genomic data that connects to whakapapa (genealogy) is considered a taonga (something precious and significant), and its storage, utilisation and interpretation is a culturally significant activity. Furthermore, such resources are not proposed to primarily enable comparisons between those with Māori and broader Pacific ancestries and other Aotearoa peoples but to place an understanding of the genetic contributors to their health outcomes in a valid context. Ongoing oversight and governance of such taonga by Māori and Pacific representatives will maximise hauora (health) while also minimising the risk of misuse of this information.

DNA and RNA-sequence based GWAS highlights membrane-transport genes as key modulators of milk lactose content

Background

Lactose provides an easily-digested energy source for neonates, and is the primary carbohydrate in milk in most species. Bovine lactose is also a key component of many human food products. However, compared to analyses of other milk components, the genetic control of lactose has been little studied. Here we present the first GWAS focussed on analysis of milk lactose traits.

Results

Using a discovery population of 12,000 taurine dairy cattle, we detail 27 QTL for lactose concentration and yield, and subsequently validate the effects of 26 of these loci in a distinct population of 18,000 cows. We next present data implicating causative genes and variants for these QTL. Fine mapping of these regions using imputed, whole genome sequence-resolution genotypes reveals protein-coding candidate causative variants affecting the ABCG2, DGAT1, STAT5B, KCNH4, NPFFR2 and RNF214 genes. Eleven of the remaining QTL appear to be driven by regulatory effects, suggested by the presence of co-locating, co-segregating eQTL discovered using mammary RNA sequence data from a population of 357 lactating cows. Pathway analysis of genes representing all lactose-associated loci shows significant enrichment of genes located in the endoplasmic reticulum, with functions related to ion channel activity mediated through the LRRC8C, P2RX4, KCNJ2 and ANKH genes. A number of the validated QTL are also found to be associated with additional milk volume, fat and protein phenotypes.

Conclusions

Overall, these findings highlight novel candidate genes and variants involved in milk lactose regulation, whose impacts on membrane transport mechanisms reinforce the key osmo-regulatory roles of lactose in milk.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4320-3) contains supplementary material, which is available to authorized users.

Compound heterozygous SLC19A3 mutations further refine the critical promoter region for biotin-thiamine-responsive basal ganglia disease

Mutations in the gene SLC19A3 result in thiamine metabolism dysfunction syndrome 2, also known as biotin-thiamine-responsive basal ganglia disease (BTBGD). This neurometabolic disease typically presents in early childhood with progressive neurodegeneration, including confusion, seizures, and dysphagia, advancing to coma and death. Treatment is possible via supplement of biotin and/or thiamine, with early treatment resulting in significant lifelong improvements. Here we report two siblings who received a refined diagnosis of BTBGD following whole-genome sequencing. Both children inherited compound heterozygous mutations from unaffected parents; a missense single-nucleotide variant (p.G23V) in the first transmembrane domain of the protein, and a 4808-bp deletion in exon 1 encompassing the 5′ UTR and minimal promoter region. This deletion is the smallest promoter deletion reported to date, further defining the minimal promoter region of SLC19A3. Unfortunately, one of the siblings died prior to diagnosis, but the other is showing significant improvement after commencement of therapy. This case demonstrates the power of whole-genome sequencing for the identification of structural variants and subsequent diagnosis of rare neurodevelopmental disorders.

Compound Heterozygous Inheritance of Mutations in Coenzyme Q8A Results in Autosomal Recessive Cerebellar Ataxia and Coenzyme Q10 Deficiency in a Female Sib-Pair

Autosomal recessive ataxias are characterised by a fundamental loss in coordination of gait with associated atrophy of the cerebellum. There is significant clinical and genetic heterogeneity amongst inherited ataxias; however, an early molecular diagnosis is essential with low-risk treatments available for some of these conditions. We describe two female siblings who presented early in life with unsteady gait and cerebellar atrophy. Whole exome sequencing revealed compound heterozygous inheritance of two pathogenic mutations (p.Leu277Pro, c.1506+1G>A) in the coenzyme Q8A gene (COQ8A), a gene central to biosynthesis of coenzyme Q (CoQ). The paternally derived p.Leu277Pro mutation is predicted to disrupt a conserved motif in the substrate-binding pocket of the protein, resulting in inhibition of CoQ₁₀ production. The maternal c.1506+1G>A mutation destroys a canonical splice donor site in exon 12 affecting transcript processing and subsequent protein translation. Mutations in this gene can result in primary coenzyme Q₁₀ deficiency type 4, which is characterized by childhood onset of cerebellar ataxia and exercise intolerance, both of which were observed in this sib-pair. Muscle biopsies revealed unequivocally low levels of CoQ_10, and the siblings were subsequently established on a therapeutic dose of CoQ₁₀ with distinct clinical evidence of improvement after 1 year of treatment. This case emphasises the importance of an early and accurate molecular diagnosis for suspected inherited ataxias, particularly given the availability of approved treatments for some subtypes.

Comparison of Huntington's disease CAG Repeat Length Stability in Human Motor Cortex and Cingulate Gyrus

Huntington's disease is caused by expansion of the CAG repeat in Huntingtin. This repeat has shown tissue-specific instability in mouse models and in a small number of post-mortem human samples. We used small-pool PCR to generate a modified instability index to quantify CAG instability within two brain regions from six human samples where cell loss has been associated with motor and mood symptoms: the motor cortex and cingulate gyrus. The expanded allele demonstrated instability in both regions, with minimal instability in the unexpanded allele. Region-specific differences were not observed, suggesting symptomatology may not be determined by repeat length instability.

Functional confirmation of PLAG1 as the candidate causative gene underlying major pleiotropic effects on body weight and milk characteristics

A major pleiotropic quantitative trait locus (QTL) located at ~25 Mbp on bovine chromosome 14 affects a myriad of growth and developmental traits in Bos taurus and indicus breeds. These QTL have been attributed to two functional variants in the bidirectional promoter of PLAG1 and CHCHD7. Although PLAG1 is a good candidate for mediating these effects, its role remains uncertain given that these variants are also associated with expression of five additional genes at the broader locus. In the current study, we conducted expression QTL (eQTL) mapping of this region using a large, high depth mammary RNAseq dataset representing 375 lactating cows. Here we show that of the seven previously implicated genes, only PLAG1 and LYN are differentially expressed by QTL genotype, and only PLAG1 bears the same association signature of the growth and body weight QTLs. For the first time, we also report significant association of PLAG1 genotype with milk production traits, including milk fat, volume, and protein yield. Collectively, these data strongly suggest PLAG1 as the causative gene underlying this diverse range of traits, and demonstrate new effects for the locus on lactation phenotypes.

Potential molecular consequences of transgene integration: The R6/2 mouse example

Integration of exogenous DNA into a host genome represents an important route to generate animal and cellular models for exploration into human disease and therapeutic development. In most models, little is known concerning structural integrity of the transgene, precise site of integration, or its impact on the host genome. We previously used whole-genome and targeted sequencing approaches to reconstruct transgene structure and integration sites in models of Huntington’s disease, revealing complex structural rearrangements that can result from transgenesis. Here, we demonstrate in the R6/2 mouse, a widely used Huntington’s disease model, that integration of a rearranged transgene with coincident deletion of 5,444 bp of host genome within the gene Gm12695 has striking molecular consequences. Gm12695, the function of which is unknown, is normally expressed at negligible levels in mouse brain, but transgene integration has resulted in cortical expression of a partial fragment (exons 8–11) 3’ to the transgene integration site in R6/2. This transcript shows significant expression among the extensive network of differentially expressed genes associated with this model, including synaptic transmission, cell signalling and transcription. These data illustrate the value of sequence-level resolution of transgene insertions and transcription analysis to inform phenotypic characterization of transgenic models utilized in therapeutic research.

Metabolite mapping reveals severe widespread perturbation of multiple metabolic processes in Huntington's disease human brain

Huntington's disease (HD) is a genetically-mediated neurodegenerative disorder wherein the aetiological defect is a mutation in the Huntington's gene (HTT), which alters the structure of the huntingtin protein (Htt) through lengthening of its polyglutamine tract, thus initiating a cascade that ultimately leads to premature death. However, neurodegeneration typically manifests in HD only in middle age, and mechanisms linking the causative mutation to brain disease are poorly understood. Brain metabolism is severely perturbed in HD, and some studies have indicated a potential role for mutant Htt as a driver of these metabolic aberrations. Here, our objective was to determine the effects of HD on brain metabolism by measuring levels of polar metabolites in regions known to undergo varying degrees of damage. We performed gas-chromatography/mass spectrometry-based metabolomic analyses in a case-control study of eleven brain regions in short post-mortem-delay human tissue from nine well-characterized HD patients and nine matched controls. In each patient, we measured metabolite content in representative tissue-samples from eleven brain regions that display varying degrees of damage in HD, thus identifying the presence and abundance of 63 different metabolites from several molecular classes, including carbohydrates, amino acids, nucleosides, and neurotransmitters. Robust alterations in regional brain-metabolite abundances were observed in HD patients: these included changes in levels of small molecules that play important roles as intermediates in the tricarboxylic-acid and urea cycles, and amino-acid metabolism. Our findings point to widespread disruption of brain metabolism and indicate a complex phenotype beyond the gradient of neuropathologic damage observed in HD brain.

Brain dopamine-serotonin vesicular transport disease presenting as a severe infantile hypotonic parkinsonian disorder

Two male siblings from a consanguineous union presented in early infancy with marked truncal hypotonia, a general paucity of movement, extrapyramidal signs and cognitive delay. By mid-childhood they had made little developmental progress and remained severely hypotonic and bradykinetic. They developed epilepsy and had problems with autonomic dysfunction and oculogyric crises. They had a number of orthopaedic problems secondary to their hypotonia. Cerebrospinal fluid (CSF) neurotransmitters were initially normal, apart from mildly elevated 5-hydroxyindolacetic acid, and the children did not respond favourably to a trial of levodopa-carbidopa. The youngest died from respiratory complications at 10 years of age. Repeat CSF neurotransmitters in the older sibling at eight years of age showed slightly low homovanillic acid and 5-hydroxyindoleacetic acid levels. Whole-exome sequencing revealed a novel mutation homozygous in both children in the monoamine transporter gene SLC18A2 (p.Pro237His), resulting in brain dopamine-serotonin vesicular transport disease. This is the second family to be described with a mutation in this gene. Treatment with the dopamine agonist pramipexole in the surviving child resulted in mild improvements in alertness, communication, and eye movements. This case supports the identification of the causal mutation in the original case, expands the clinical phenotype of brain dopamine-serotonin vesicular transport disease and confirms that pramipexole treatment may lead to symptomatic improvement in affected individuals.

Metabolic disruption identified in the Huntington's disease transgenic sheep model

Huntington’s disease (HD) is a dominantly inherited, progressive neurodegenerative disorder caused by a CAG repeat expansion within exon 1 of HTT, encoding huntingtin. There are no therapies that can delay the progression of this devastating disease. One feature of HD that may play a critical role in its pathogenesis is metabolic disruption. Consequently, we undertook a comparative study of metabolites in our transgenic sheep model of HD (OVT73). This model does not display overt symptoms of HD but has circadian rhythm alterations and molecular changes characteristic of the early phase disease. Quantitative metabolite profiles were generated from the motor cortex, hippocampus, cerebellum and liver tissue of 5 year old transgenic sheep and matched controls by gas chromatography-mass spectrometry. Differentially abundant metabolites were evident in the cerebellum and liver. There was striking tissue-specificity, with predominantly amino acids affected in the transgenic cerebellum and fatty acids in the transgenic liver, which together may indicate a hyper-metabolic state. Furthermore, there were more strong pair-wise correlations of metabolite abundance in transgenic than in wild-type cerebellum and liver, suggesting altered metabolic constraints. Together these differences indicate a metabolic disruption in the sheep model of HD and could provide insight into the presymptomatic human disease.

Phenotypic population screen identifies a new mutation in bovine DGAT1 responsible for unsaturated milk fat

Selective breeding has strongly reduced the genetic diversity in livestock species, and contemporary breeding practices exclude potentially beneficial rare genetic variation from the future gene pool. Here we test whether important traits arising by new mutations can be identified and rescued in highly selected populations. We screened milks from 2.5 million cows to identify an exceptional individual which produced milk with reduced saturated fat content, and improved unsaturated and omega-3 fatty acid concentrations. The milk traits were transmitted dominantly to her offspring, and genetic mapping and genome sequencing revealed a new mutation in a previously unknown splice enhancer of the DGAT1 gene. Homozygous carriers show features of human diarrheal disorders, and may be useful for the development of therapeutic strategies. Our study demonstrates that high-throughput phenotypic screening can uncover rich genetic diversity even in inbred populations, and introduces a novel strategy to develop novel milks with improved nutritional properties.