Mutation discovery by targeted genomic enrichment of multiplexed barcoded samples.

Targeted genomic enrichment followed by next-generation DNA sequencing has dramatically increased efficiency of mutation-discovery efforts. We describe a protocol for genomic enrichment of pooled barcoded samples in a single assay that increases experimental flexibility and efficiency. We screened 770 genes (1.4 megabases) in thirty N-ethyl-N-nitrosourea (ENU)-mutagenized rats and identified known variants at >96% sensitivity as well as new mutations at a false positive rate < 1 in 8 megabases.

A hypomorphic mutation in Lpin1 induces progressively improving neuropathy and lipodystrophy in the rat.

The Lpin1 gene encodes the phosphatidate phosphatase (PAP1) enzyme Lipin 1, which plays a critical role in lipid metabolism. In this study we describe the identification and characterization of a rat model with a mutated Lpin1 gene (Lpin1(1Hubr)), generated by N-ethyl-N-nitrosourea mutagenesis. Lpin1(1Hubr) rats are characterized by hindlimb paralysis and mild lipodystrophy that are detectable from the second postnatal week. Sequencing of Lpin1 identified a point mutation in the 5'-end splice site of intron 18 resulting in mis-splicing, a reading frameshift, and a premature stop codon. As this mutation does not induce nonsense-mediated decay, it allows the production of a truncated Lipin 1 protein lacking PAP1 activity. Lpin1(1Hubr) rats developed hypomyelination and mild lipodystrophy rather than the pronounced demyelination and adipocyte defects characteristic of Lpin1(fld/fld) mice, which carry a null allele for Lpin1. Furthermore, biochemical, histological, and molecular analyses revealed that these lesions improve in older Lpin1(1Hubr) rats as compared with young Lpin1(1Hubr) rats and Lpin1(fld/fld) mice. We observed activation of compensatory biochemical pathways substituting for missing PAP1 activity that, in combination with a possible non-enzymatic Lipin 1 function residing outside of its PAP1 domain, may contribute to the less severe phenotypes observed in Lpin1(1Hubr) rats as compared with Lpin1(fld/fld) mice. Although we are cautious in making a direct parallel between the presented rodent model and human disease, our data may provide new insight into the pathogenicity of recently identified human LPIN1 mutations.

Homozygous and heterozygous p53 knockout rats develop metastasizing sarcomas with high frequency.

The TP53 tumor suppressor gene is mutated in the majority of human cancers. Inactivation of p53 in a variety of animal models results in early-onset tumorigenesis, reflecting the importance of p53 as a gatekeeper tumor suppressor. We generated a mutant Tp53 allele in the rat using a target-selected mutagenesis approach. Here, we report that homozygosity for this allele results in complete loss of p53 function. Homozygous mutant rats predominantly develop sarcomas with an onset of 4 months of age with a high occurrence of pulmonary metastases. Heterozygous rats develop sarcomas starting at 8 months of age. Molecular analysis revealed that these tumors exhibit a loss-of-heterozygosity of the wild-type Tp53 allele. These unique features make this rat highly complementary to other rodent p53 knockout models and a versatile tool for investigating tumorigenesis processes as well as genotoxic studies.

Pmch expression during early development is critical for normal energy homeostasis.

Postnatal development and puberty are times of strong physical maturation and require large quantities of energy. The hypothalamic neuropeptide melanin-concentrating hormone (MCH) regulates nutrient intake and energy homeostasis, but the underlying mechanisms are not completely understood. Here we use a novel rat knockout model in which the MCH precursor Pmch has been inactivated to study the effects of loss of MCH on energy regulation in more detail. Pmch(-/-) rats are lean, hypophagic, osteoporotic, and although endocrine parameters were changed in pmch(-/-) rats, endocrine dynamics were normal, indicating an adaptation to new homeostatic levels rather than disturbed metabolic mechanisms. Detailed body weight growth and feeding behavior analysis revealed that Pmch expression is particularly important during early rat development and puberty, i.e., the first 8 postnatal weeks. Loss of Pmch resulted in a 20% lower set point for body weight that was determined solely during this period and remained unchanged during adulthood. Although the final body weight is diet dependent, the Pmch-deficiency effect was similar for all diets tested in this study. Loss of Pmch affected energy expenditure in both young and adult rats, although these effects seem secondary to the observed hypophagia. Our findings show an important role for Pmch in energy homeostasis determination during early development and indicate that the MCH receptor 1 system is a plausible target for childhood obesity treatment, currently a major health issue in first world countries.

Calorie Restriction Increases the Number of Competing Stem Cells and Decreases Mutation Retention in the Intestine.

Calorie restriction (CR) extends lifespan through several intracellular mechanisms, including increased DNA repair, leading to fewer DNA mutations that cause age-related pathologies. However, it remains unknown how CR acts on mutation retention at the tissue level. Here, we use Cre-mediated DNA recombination of the confetti reporter as proxy for neutral mutations and follow these mutations by intravital microscopy to identify how CR affects retention of mutations in the intestine. We find that CR leads to increased numbers of functional Lgr5+ stem cells that compete for niche occupancy, resulting in slower but stronger stem cell competition. Consequently, stem cells carrying neutral or Apc mutations encounter more wild-type competitors, thus increasing the chance that they get displaced from the niche to get lost over time. Thus, our data show that CR not only affects the acquisition of mutations but also leads to lower retention of mutations in the intestine.

Lack of DNA mismatch repair protein MSH6 in the rat results in hereditary non-polyposis colorectal cancer-like tumorigenesis.

To understand genetic instability in relation to tumorigenesis, experimental animal models have proven very useful. The DNA mismatch repair (MMR) machinery safeguards genomic integrity by repairing mismatches, insertion or deletion loops and responding to genotoxic agents. Here, we describe the functional characterization of a novel rat mutant model in which the MMR gene Msh6 has been genetically inactivated by N-ethyl-N-nitrosourea-driven target-selected mutagenesis. This model shows a robust mutator phenotype that is reflected by microsatellite instability and an increased germ line point mutation frequency. Consequently, these rats develop a spectrum of tumors with a high similarity to atypical hereditary non-polyposis colorectal cancer in humans. The MSH6 knockout rat complements existing models for studying genetic instable tumorigenesis as it provides experimental opportunities that are not available or suboptimal in current models.

Improved generation of rat gene knockouts by target-selected mutagenesis in mismatch repair-deficient animals.

BACKGROUND: The laboratory rat (Rattus norvegicus) is one of the preferred model organisms in physiological and pharmacological research, although the availability of specific genetic models, especially gene knockouts, is limited. N-ethyl-N-nitrosourea (ENU)-driven target-selected mutagenesis is currently the most successful method in rats, although it is still very laborious and expensive. RESULTS: As ENU-induced DNA damage is normally recognized by the mismatch repair (MMR) system, we hypothesized that the effectiveness of the target-selected mutagenesis approach could be improved by using a MMR-deficient genetic background. Indeed, Msh6 knockout rats were found to be more sensitive to ENU treatment and the germ line mutation rate was boosted more than two-fold to 1 mutation per 585 kb. In addition, the molecular mutation spectrum was found to be changed in favor of generating knockout-type alleles by approximately 20%, resulting in an overall increase in efficiency of approximately 2.5 fold. The improved effectiveness was demonstrated by high throughput mutation discovery in 70 Mb of sequence in a set of only 310 mutant F1 rats. This resulted in the identification of 89 mutations of which four introduced a premature stopcodon and 64 resulted in amino acid changes. CONCLUSION: Taken together, we show that the use of a MMR-deficient background considerably improves ENU-driven target-selected mutagenesis in the rat, thereby reducing animal use as well as screening costs. The use of a mismatch repair-deficient genetic background for improving mutagenesis and target-selected knockout efficiency is in principle applicable to any organism of interest.

The Genomic Scrapheap Challenge; Extracting Relevant Data from Unmapped Whole Genome Sequencing Reads, Including Strain Specific Genomic Segments, in Rats.

Unmapped next-generation sequencing reads are typically ignored while they contain biologically relevant information. We systematically analyzed unmapped reads from whole genome sequencing of 33 inbred rat strains. High quality reads were selected and enriched for biologically relevant sequences; similarity-based analysis revealed clustering similar to previously reported phylogenetic trees. Our results demonstrate that on average 20% of all unmapped reads harbor sequences that can be used to improve reference genomes and generate hypotheses on potential genotype-phenotype relationships. Analysis pipelines would benefit from incorporating the described methods and reference genomes would benefit from inclusion of the genomic segments obtained through these efforts.

Generation and characterization of rat liver stem cell lines and their engraftment in a rat model of liver failure.

The rat is an important model for liver regeneration. However, there is no in vitro culture system that can capture the massive proliferation that can be observed after partial hepatectomy in rats. We here describe the generation of rat liver stem cell lines. Rat liver stem cells, which grow as cystic organoids, were characterized by high expression of the stem cell marker Lgr5, by the expression of liver progenitor and duct markers, and by low expression of hepatocyte markers, oval cell markers, and stellate cell markers. Prolonged cultures of rat liver organoids depended on high levels of WNT-signalling and the inhibition of BMP-signaling. Upon transplantation of clonal lines to a Fah(-/-) Il2rg(-/-) rat model of liver failure, the rat liver stem cells engrafted into the host liver where they differentiated into areas with FAH and Albumin positive hepatocytes. Rat liver stem cell lines hold potential as consistent reliable cell sources for pharmacological, toxicological or metabolic studies. In addition, rat liver stem cell lines may contribute to the development of regenerative medicine in liver disease. To our knowledge, the here described liver stem cell lines represent the first organoid culture system in the rat.

Lack of major genome instability in tumors of p53 null rats.

Tumorigenesis is often associated with loss of tumor suppressor genes (such as TP53), genomic instability and telomere lengthening. Previously, we generated and characterized a rat p53 knockout model in which the homozygous rats predominantly develop hemangiosarcomas whereas the heterozygous rats mainly develop osteosarcomas. Using genome-wide analyses, we find that the tumors that arise in the heterozygous and homozygous Tp53C273X mutant animals are also different in their genomic instability profiles. While p53 was fully inactivated in both heterozygous and homozygous knockout rats, tumors from homozygous animals show very limited aneuploidy and low degrees of somatic copy number variation as compared to the tumors from heterozygous animals. In addition, complex structural rearrangements such as chromothripsis and breakage-fusion-bridge cycles were never found in tumors from homozygous animals, while these were readily detectable in tumors from heterozygous animals. Finally, we measured telomere length and telomere lengthening pathway activity and found that tumors of homozygous animals have longer telomeres but do not show clear telomerase or alternative lengthening of telomeres (ALT) activity differences as compared to the tumors from heterozygous animals. Taken together, our results demonstrate that host p53 status in this rat p53 knockout model has a large effect on both tumor type and genomic instability characteristics, where full loss of functional p53 is not the main driver of large-scale structural variations. Our results also suggest that chromothripsis primarily occurs under p53 heterozygous rather than p53 null conditions.

Large-scale identification of coregulated enhancer networks in the adult human brain.

Understanding the complexity of the human brain and its functional diversity remain a major challenge. Distinct anatomical regions are involved in an array of processes, including organismal homeostasis, cognitive functions, and susceptibility to neurological pathologies, many of which define our species. Distal enhancers have emerged as key regulatory elements that acquire histone modifications in a cell- and species-specific manner, thus enforcing specific gene expression programs. Here, we survey the epigenomic landscape of promoters and cis-regulatory elements in 136 regions of the adult human brain. We identify a total of 83,553 promoter-distal H3K27ac-enriched regions showing global characteristics of brain enhancers. We use coregulation of enhancer elements across many distinct regions of the brain to uncover functionally distinct networks at high resolution and link these networks to specific neuroglial functions. Furthermore, we use these data to understand the relevance of noncoding genomic variations previously linked to Parkinson's disease incidence.

Olanzapine-induced changes in glucose metabolism are independent of the melanin-concentrating hormone system.

Atypical antipsychotic drugs such as Olanzapine (Ola) induce weight gain and metabolic changes associated with the development of type 2 diabetes. The mechanisms underlying these undesired side-effects are currently unknown. Chagnon et al. showed that the common allele rs7973796 of the prepro-melanin-concentrating hormone (PMCH) gene is associated with a greater body mass index in Ola-treated schizophrenic patients. As PMCH encodes for the orexigenic neuropeptide melanin-concentrating hormone (MCH), it was hypothesized that MCH is involved in Ola-induced metabolic changes. We have recently reported that the intragastric infusion of Ola results in hyperglycaemia and insulin resistance in male rats. In order to test in vivo the possible involvement of the PMCH gene in the pathogenesis of Ola side-effects, we administered Ola intragastrically in wild-type (WT) and PMCH knock-out (KO) rats. Our results show that glucose and corticosterone levels, as well as endogenous glucose production, are elevated by the infusion of Ola in both WT and KO animals. Thus, the lack of MCH does not seem to affect the acute effects of Ola on glucose metabolism. On the other hand, these effects might be obliterated by compensatory changes in other hypothalamic systems. In addition, possible modulatory effects of the MCH KO on the long term effects of Ola, i.e. increased adiposity, body weight gain, have not been investigated yet.

Pmch-deficiency in rats is associated with normal adipocyte differentiation and lower sympathetic adipose drive.

The orexigenic neuropeptide melanin-concentrating hormone (MCH), a product of Pmch, is an important mediator of energy homeostasis. Pmch-deficient rodents are lean and smaller, characterized by lower food intake, body-, and fat mass. Pmch is expressed in hypothalamic neurons that ultimately are components in the sympathetic nervous system (SNS) drive to white and interscapular brown adipose tissue (WAT, iBAT, respectively). MCH binds to MCH receptor 1 (MCH1R), which is present on adipocytes. Currently it is unknown if Pmch-ablation changes adipocyte differentiation or sympathetic adipose drive. Using Pmch-deficient and wild-type rats on a standard low-fat diet, we analyzed dorsal subcutaneous and perirenal WAT mass and adipocyte morphology (size and number) throughout development, and indices of sympathetic activation in WAT and iBAT during adulthood. Moreover, using an in vitro approach we investigated the ability of MCH to modulate 3T3-L1 adipocyte differentiation. Pmch-deficiency decreased dorsal subcutaneous and perirenal WAT mass by reducing adipocyte size, but not number. In line with this, in vitro 3T3-L1 adipocyte differentiation was unaffected by MCH. Finally, adult Pmch-deficient rats had lower norepinephrine turnover (an index of sympathetic adipose drive) in WAT and iBAT than wild-type rats. Collectively, our data indicate that MCH/MCH1R-pathway does not modify adipocyte differentiation, whereas Pmch-deficiency in laboratory rats lowers adiposity throughout development and sympathetic adipose drive during adulthood.

Systematic biases in DNA copy number originate from isolation procedures.

BACKGROUND: The ability to accurately detect DNA copy number variation in both a sensitive and quantitative manner is important in many research areas. However, genome-wide DNA copy number analyses are complicated by variations in detection signal. RESULTS: While GC content has been used to correct for this, here we show that coverage biases are tissue-specific and independent of the detection method as demonstrated by next-generation sequencing and array CGH. Moreover, we show that DNA isolation stringency affects the degree of equimolar coverage and that the observed biases coincide with chromatin characteristics like gene expression, genomic isochores, and replication timing. CONCLUSION: These results indicate that chromatin organization is a main determinant for differential DNA retrieval. These findings are highly relevant for germline and somatic DNA copy number variation analyses.

Quantitative and qualitative proteome characteristics extracted from in-depth integrated genomics and proteomics analysis.

Quantitative and qualitative protein characteristics are regulated at genomic, transcriptomic, and posttranscriptional levels. Here, we integrated in-depth transcriptome and proteome analyses of liver tissues from two rat strains to unravel the interactions within and between these layers. We obtained peptide evidence for 26,463 rat liver proteins. We validated 1,195 gene predictions, 83 splice events, 126 proteins with nonsynonymous variants, and 20 isoforms with nonsynonymous RNA editing. Quantitative RNA sequencing and proteomics data correlate highly between strains but poorly among each other, indicating extensive nongenetic regulation. Our multilevel analysis identified a genomic variant in the promoter of the most differentially expressed gene Cyp17a1, a previously reported top hit in genome-wide association studies for human hypertension, as a potential contributor to the hypertension phenotype in SHR rats. These results demonstrate the power of and need for integrative analysis for understanding genetic control of molecular dynamics and phenotypic diversity in a system-wide manner.

Combined sequence-based and genetic mapping analysis of complex traits in outbred rats.

Genetic mapping on fully sequenced individuals is transforming understanding of the relationship between molecular variation and variation in complex traits. Here we report a combined sequence and genetic mapping analysis in outbred rats that maps 355 quantitative trait loci for 122 phenotypes. We identify 35 causal genes involved in 31 phenotypes, implicating new genes in models of anxiety, heart disease and multiple sclerosis. The relationship between sequence and genetic variation is unexpectedly complex: at approximately 40% of quantitative trait loci, a single sequence variant cannot account for the phenotypic effect. Using comparable sequence and mapping data from mice, we show that the extent and spatial pattern of variation in inbred rats differ substantially from those of inbred mice and that the genetic variants in orthologous genes rarely contribute to the same phenotype in both species.

Melanocortin receptor 4 deficiency affects body weight regulation, grooming behavior, and substrate preference in the rat.

Obesity is caused by an imbalance between energy intake and expenditure and has become a major health-care problem in western society. The central melanocortin system plays a crucial role in the regulation of feeding and energy expenditure, and functional loss of melanocortin receptor 4 (MC4R) is the most common genetic cause of human obesity. In this study, we present the first functional Mc4r knockout model in the rat, resulting from an N-ethyl-N-nitrosourea mutagenesis-induced point mutation. In vitro observations revealed impaired membrane-binding and subsequent nonfunctionality of the receptor, whereas in vivo observations showed that functional loss of MC4R increased body weight, food intake, white adipose mass, and changed substrate preference. In addition, intracerebroventricular (ICV) administration of Agouti-Related Protein(79-129) (AgRP(79-129)), an MC4R inverse agonist, or Melanotan-II (MTII), an MC4R agonist, did affect feeding behavior in wild-type rats but not in homozygous mutant rats, confirming complete loss of MC4R function in vivo. Finally, ICV administration of MTII induced excessive grooming behavior in wild-type rats, whereas this effect was absent in homozygous mutant rats, indicating that MTII-induced grooming behavior is exclusively regulated via MC4R pathways. Taken together, we expect that the MC4R rat model described here will be a valuable tool for studying monogenic obesity in humans. More specifically, the relative big size and increased cognitive capacity of rats as compared to mice will facilitate complex behavioral studies and detailed mechanistic studies regarding central function of MC4R, both of which ultimately may help to further understand the specific mechanisms that induce obesity during loss of MC4R function.

Chronic loss of melanin-concentrating hormone affects motivational aspects of feeding in the rat.

Current epidemic obesity levels apply great medical and financial pressure to the strenuous economy of obesity-prone cultures, and neuropeptides involved in body weight regulation are regarded as attractive targets for a possible treatment of obesity in humans. The lateral hypothalamus and the nucleus accumbens shell (AcbSh) form a hypothalamic-limbic neuropeptide feeding circuit mediated by Melanin-Concentrating Hormone (MCH). MCH promotes feeding behavior via MCH receptor-1 (MCH1R) in the AcbSh, although this relationship has not been fully characterized. Given the AcbSh mediates reinforcing properties of food, we hypothesized that MCH modulates motivational aspects of feeding.Here we show that chronic loss of the rat MCH-precursor Pmch decreased food intake predominantly via a reduction in meal size during rat development and reduced high-fat food-reinforced operant responding in adult rats. Moreover, acute AcbSh administration of Neuropeptide-GE and Neuropeptide-EI (NEI), both additional neuropeptides derived from Pmch, or chronic intracerebroventricular infusion of NEI, did not affect feeding behavior in adult pmch(+/+) or pmch(-/-) rats. However, acute administration of MCH to the AcbSh of adult pmch(-/-) rats elevated feeding behavior towards wild type levels. Finally, adult pmch(-/-) rats showed increased ex vivo electrically evoked dopamine release and increased limbic dopamine transporter levels, indicating that chronic loss of Pmch in the rat affects the limbic dopamine system.Our findings support the MCH-MCH1R system as an amplifier of consummatory behavior, confirming this system as a possible target for the treatment of obesity. We propose that MCH-mediated signaling in the AcbSh positively mediates motivational aspects of feeding behavior. Thereby it provides a crucial signal by which hypothalamic neural circuits control energy balance and guide limbic brain areas to enhance motivational or incentive-related aspects of food consumption.