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.

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.

Natural cryptic variation in epigenetic modulation of an embryonic gene regulatory network.

Gene regulatory networks (GRNs) that direct animal embryogenesis must respond to varying environmental and physiological conditions to ensure robust construction of organ systems. While GRNs are evolutionarily modified by natural genomic variation, the roles of epigenetic processes in shaping plasticity of GRN architecture are not well understood. The endoderm GRN in Caenorhabditis elegans is initiated by the maternally supplied SKN-1/Nrf2 bZIP transcription factor; however, the requirement for SKN-1 in endoderm specification varies widely among distinct C. elegans wild isotypes, owing to rapid developmental system drift driven by accumulation of cryptic genetic variants. We report here that heritable epigenetic factors that are stimulated by transient developmental diapause also underlie cryptic variation in the requirement for SKN-1 in endoderm development. This epigenetic memory is inherited from the maternal germline, apparently through a nuclear, rather than cytoplasmic, signal, resulting in a parent-of-origin effect (POE), in which the phenotype of the progeny resembles that of the maternal founder. The occurrence and persistence of POE varies between different parental pairs, perduring for at least 10 generations in one pair. This long-perduring POE requires piwi-interacting RNA (piRNA) function and the germline nuclear RNA interference (RNAi) pathway, as well as MET-2 and SET-32, which direct histone H3K9 trimethylation and drive heritable epigenetic modification. Such nongenetic cryptic variation may provide a resource of additional phenotypic diversity through which adaptation may facilitate evolutionary changes and shape developmental regulatory systems.

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.

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

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

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.

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.

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.

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.

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.

Brain urea increase is an early Huntington's disease pathogenic event observed in a prodromal transgenic sheep model and HD cases.

The neurodegenerative disorder Huntington's disease (HD) is typically characterized by extensive loss of striatal neurons and the midlife onset of debilitating and progressive chorea, dementia, and psychological disturbance. HD is caused by a CAG repeat expansion in the Huntingtin (HTT) gene, translating to an elongated glutamine tract in the huntingtin protein. The pathogenic mechanism resulting in cell dysfunction and death beyond the causative mutation is not well defined. To further delineate the early molecular events in HD, we performed RNA-sequencing (RNA-seq) on striatal tissue from a cohort of 5-y-old OVT73-line sheep expressing a human CAG-expansion HTT cDNA transgene. Our HD OVT73 sheep are a prodromal model and exhibit minimal pathology and no detectable neuronal loss. We identified significantly increased levels of the urea transporter SLC14A1 in the OVT73 striatum, along with other important osmotic regulators. Further investigation revealed elevated levels of the metabolite urea in the OVT73 striatum and cerebellum, consistent with our recently published observation of increased urea in postmortem human brain from HD cases. Extending that finding, we demonstrate that postmortem human brain urea levels are elevated in a larger cohort of HD cases, including those with low-level neuropathology (Vonsattel grade 0/1). This elevation indicates increased protein catabolism, possibly as an alternate energy source given the generalized metabolic defect in HD. Increased urea and ammonia levels due to dysregulation of the urea cycle are known to cause neurologic impairment. Taken together, our findings indicate that aberrant urea metabolism could be the primary biochemical disruption initiating neuropathogenesis in HD.

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.

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.

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.

Identification of an immune modulation locus utilising a bovine mammary gland infection challenge model.

Inflammation of the mammary gland following bacterial infection, commonly known as mastitis, affects all mammalian species. Although the aetiology and epidemiology of mastitis in the dairy cow are well described, the genetic factors mediating resistance to mammary gland infection are not well known, due in part to the difficulty in obtaining robust phenotypic information from sufficiently large numbers of individuals. To address this problem, an experimental mammary gland infection experiment was undertaken, using a Friesian-Jersey cross breed F2 herd. A total of 604 animals received an intramammary infusion of Streptococcus uberis in one gland, and the clinical response over 13 milkings was used for linkage mapping and genome-wide association analysis. A quantitative trait locus (QTL) was detected on bovine chromosome 11 for clinical mastitis status using micro-satellite and Affymetrix 10 K SNP markers, and then exome and genome sequence data used from the six F1 sires of the experimental animals to examine this region in more detail. A total of 485 sequence variants were typed in the QTL interval, and association mapping using these and an additional 37 986 genome-wide markers from the Illumina SNP50 bovine SNP panel revealed association with markers encompassing the interleukin-1 gene cluster locus. This study highlights a region on bovine chromosome 11, consistent with earlier studies, as conferring resistance to experimentally induced mammary gland infection, and newly prioritises the IL1 gene cluster for further analysis in genetic resistance to mastitis.

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.

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.

Early and progressive circadian abnormalities in Huntington's disease sheep are unmasked by social environment.

Insidious changes in behaviour herald the onset of progressive neurodegenerative disorders such as Huntington's disease (HD), sometimes years before overt symptoms are seen. Sleep and circadian disturbances are particularly disruptive symptoms in patients with neurological disorders, but they are difficult to measure in humans. Here we studied circadian behaviour in transgenic HD sheep expressing the full-length human huntingtin protein with an expanded CAG repeat mutation in the juvenile range. Young HD sheep with no other symptoms exhibited circadian behavioural abnormalities that worsened with age. The most obvious change was a disturbed evening behaviour reminiscent of 'sundowning' that is seen in some patients with dementia. There were no structural abnormalities seen with magnetic resonance imaging, even in 5-year-old HD sheep. Interestingly, detection of the circadian abnormalities depended upon their social grouping. Abnormalities emerged in sheep kept in an 'HD-only' flock, whereas the behaviour of HD sheep kept mixed with normal sheep was relatively normal. Sleep-wake abnormalities in HD patients are also likely to be hidden, and may precede overt symptoms by many years. Sleep disruption has deleterious effects, even in normal people. The knock-on effects of sleep-wake disturbance may exacerbate, or even cause symptoms such as irritability and depression that are common in early stage HD patients. HD sheep will be useful models for probing the mechanisms underlying circadian behavioural disorder in HD.