Enrico Petretto

Copy number polymorphism in Fcgr3 predisposes to glomerulonephritis in rats and humans

Identification of the genes underlying complex phenotypes and the definition of the evolutionary forces that have shaped eukaryotic genomes are among the current challenges in molecular genetics. Variation in gene copy number is increasingly recognized as a source of inter-individual differences in genome sequence and has been proposed as a driving force for genome evolution and phenotypic variation. Here we show that copy number variation of the orthologous rat and human Fcgr3 genes is a determinant of susceptibility to immunologically mediated glomerulonephritis. Positional cloning identified loss of the newly described, rat-specific Fcgr3 paralogue, Fcgr3-related sequence (Fcgr3-rs), as a determinant of macrophage overactivity and glomerulonephritis in Wistar Kyoto rats. In humans, low copy number of FCGR3B, an orthologue of rat Fcgr3, was associated with glomerulonephritis in the autoimmune disease systemic lupus erythematosus. The finding that gene copy number polymorphism predisposes to immunologically mediated renal disease in two mammalian species provides direct evidence for the importance of genome plasticity in the evolution of genetically complex phenotypes, including susceptibility to common human disease.

Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease

Integration of genome-wide expression profiling with linkage analysis is a new approach to identifying genes underlying complex traits. We applied this approach to the regulation of gene expression in the BXH/HXB panel of rat recombinant inbred strains, one of the largest available rodent recombinant inbred panels and a leading resource for genetic analysis of the highly prevalent metabolic syndrome. In two tissues important to the pathogenesis of the metabolic syndrome, we mapped cis- and trans-regulatory control elements for expression of thousands of genes across the genome. Many of the most highly linked expression quantitative trait loci are regulated in cis, are inherited essentially as monogenic traits and are good candidate genes for previously mapped physiological quantitative trait loci in the rat. By comparative mapping we generated a data set of 73 candidate genes for hypertension that merit testing in human populations. Mining of this publicly available data set is expected to lead to new insights into the genes and regulatory pathways underlying the extensive range of metabolic and cardiovascular disease phenotypes that segregate in these recombinant inbred strains.

Repeatability of published microarray gene expression analyses

Given the complexity of microarray-based gene expression studies, guidelines encourage transparent design and public data availability. Several journals require public data deposition and several public databases exist. However, not all data are publicly available, and even when available, it is unknown whether the published results are reproducible by independent scientists. Here we evaluated the replication of data analyses in 18 articles on microarray-based gene expression profiling published in Nature Genetics in 2005–2006. One table or figure from each article was independently evaluated by two teams of analysts. We reproduced two analyses in principle and six partially or with some discrepancies; ten could not be reproduced. The main reason for failure to reproduce was data unavailability, and discrepancies were mostly due to incomplete data annotation or specification of data processing and analysis. Repeatability of published microarray studies is apparently limited. More strict publication rules enforcing public data availability and explicit description of data processing and analysis should be considered.

A trans-acting locus regulates an anti-viral expression network and type 1 diabetes risk

Combined analyses of gene networks and DNA sequence variation can provide new insights into the aetiology of common diseases that may not be apparent from genome-wide association studies alone. Recent advances in rat genomics are facilitating systems-genetics approaches. Here we report the use of integrated genome-wide approaches across seven rat tissues to identify gene networks and the loci underlying their regulation. We defined an interferon regulatory factor 7 (IRF7)-driven inflammatory network (IDIN) enriched for viral response genes, which represents a molecular biomarker for macrophages and which was regulated in multiple tissues by a locus on rat chromosome 15q25. We show that Epstein–Barr virus induced gene 2 (Ebi2, also known as Gpr183), which lies at this locus and controls B lymphocyte migration, is expressed in macrophages and regulates the IDIN. The human orthologous locus on chromosome 13q32 controlled the human equivalent of the IDIN, which was conserved in monocytes. IDIN genes were more likely to associate with susceptibility to type 1 diabetes (T1D)—a macrophage-associated autoimmune disease—than randomly selected immune response genes (P = 8.85 × 10−6). The human locus controlling the IDIN was associated with the risk of T1D at single nucleotide polymorphism rs9585056 (P = 7.0 × 10−10; odds ratio, 1.15), which was one of five single nucleotide polymorphisms in this region associated with EBI2 (GPR183) expression. These data implicate IRF7 network genes and their regulatory locus in the pathogenesis of T1D.

Integrated genomic approaches implicate osteoglycin (Ogn) in the regulation of left ventricular mass

Left ventricular mass (LVM) and cardiac gene expression are complex traits regulated by factors both intrinsic and extrinsic to the heart. To dissect the major determinants of LVM, we combined expression quantitative trait locus and quantitative trait transcript (QTT) analyses of the cardiac transcriptome in the rat. Using these methods and in vitro functional assays, we identified osteoglycin (Ogn) as a major candidate regulator of rat LVM, with increased Ogn protein expression associated with elevated LVM. We also applied genome-wide QTT analysis to the human heart and observed that, out of ∼22,000 transcripts, OGN transcript abundance had the highest correlation with LVM. We further confirmed a role for Ogn in the in vivo regulation of LVM in Ogn knockout mice. Taken together, these data implicate Ogn as a key regulator of LVM in rats, mice and humans, and suggest that Ogn modifies the hypertrophic response to extrinsic factors such as hypertension and aortic stenosis.

Genome Sequencing Reveals Loci under Artificial Selection that Underlie Disease Phenotypes in the Laboratory Rat

Summary Large numbers of inbred laboratory rat strains have been developed for a range of complex disease phenotypes. To gain insights into the evolutionary pressures underlying selection for these phenotypes, we sequenced the genomes of 27 rat strains, including 11 models of hypertension, diabetes, and insulin resistance, along with their respective control strains. Altogether, we identified more than 13 million single-nucleotide variants, indels, and structural variants across these rat strains. Analysis of strain-specific selective sweeps and gene clusters implicated genes and pathways involved in cation transport, angiotensin production, and regulators of oxidative stress in the development of cardiovascular disease phenotypes in rats. Many of the rat loci that we identified overlap with previously mapped loci for related traits in humans, indicating the presence of shared pathways underlying these phenotypes in rats and humans. These data represent a step change in resources available for evolutionary analysis of complex traits in disease models. PaperClip

Gene copy number variation and common human disease

Fanciulli M, Petretto E, Aitman TJ. Gene copy number variation and common human disease.

Endonuclease G is a novel determinant of cardiac hypertrophy and mitochondrial function.

Left ventricular mass (LVM) is a highly heritable trait and an independent risk factor for all-cause mortality. So far, genome-wide association studies have not identified the genetic factors that underlie LVM variation, and the regulatory mechanisms for blood-pressure-independent cardiac hypertrophy remain poorly understood. Unbiased systems genetics approaches in the rat now provide a powerful complementary tool to genome-wide association studies, and we applied integrative genomics to dissect a highly replicated, blood-pressure-independent LVM locus on rat chromosome 3p. Here we identified endonuclease G (Endog), which previously was implicated in apoptosis but not hypertrophy, as the gene at the locus, and we found a loss-of-function mutation in Endog that is associated with increased LVM and impaired cardiac function. Inhibition of Endog in cultured cardiomyocytes resulted in an increase in cell size and hypertrophic biomarkers in the absence of pro-hypertrophic stimulation. Genome-wide network analysis unexpectedly implicated ENDOG in fundamental mitochondrial processes that are unrelated to apoptosis. We showed direct regulation of ENDOG by ERR-α and PGC1α (which are master regulators of mitochondrial and cardiac function), interaction of ENDOG with the mitochondrial genome and ENDOG-mediated regulation of mitochondrial mass. At baseline, the Endog-deleted mouse heart had depleted mitochondria, mitochondrial dysfunction and elevated levels of reactive oxygen species, which were associated with enlarged and steatotic cardiomyocytes. Our study has further established the link between mitochondrial dysfunction, reactive oxygen species and heart disease and has uncovered a role for Endog in maladaptive cardiac hypertrophy.

The genome sequence of the spontaneously hypertensive rat: Analysis and functional significance

The spontaneously hypertensive rat (SHR) is the most widely studied animal model of hypertension. Scores of SHR quantitative loci (QTLs) have been mapped for hypertension and other phenotypes. We have sequenced the SHR/OlaIpcv genome at 10.7-fold coverage by paired-end sequencing on the Illumina platform. We identified 3.6 million high-quality single nucleotide polymorphisms (SNPs) between the SHR/OlaIpcv and Brown Norway (BN) reference genome, with a high rate of validation (sensitivity 96.3%-98.0% and specificity 99%-100%). We also identified 343,243 short indels between the SHR/OlaIpcv and reference genomes. These SNPs and indels resulted in 161 gain or loss of stop codons and 629 frameshifts compared with the BN reference sequence. We also identified 13,438 larger deletions that result in complete or partial absence of 107 genes in the SHR/OlaIpcv genome compared with the BN reference and 588 copy number variants (CNVs) that overlap with the gene regions of 688 genes. Genomic regions containing genes whose expression had been previously mapped as cis-regulated expression quantitative trait loci (eQTLs) were significantly enriched with SNPs, short indels, and larger deletions, suggesting that some of these variants have functional effects on gene expression. Genes that were affected by major alterations in their coding sequence were highly enriched for genes related to ion transport, transport, and plasma membrane localization, providing insights into the likely molecular and cellular basis of hypertension and other phenotypes specific to the SHR strain. This near complete catalog of genomic differences between two extensively studied rat strains provides the starting point for complete elucidation, at the molecular level, of the physiological and pathophysiological phenotypic differences between individuals from these strains.

A new essential hypertension susceptibility locus on chromosome 2p24-p25, detected by genomewide search.

Essential hypertension (EH) is a complex disorder that results from the interaction of a number of susceptibility genes and environmental factors. We studied an isolated Sardinian village (Talana) in which the prevalence of hypertension is comparable to that in most Western populations. Talana exhibits features, such as slow demographic growth, high inbreeding, a low number of founders, stable lifestyle and culture, and accurate genealogical records, that make it suitable for the study of complex disorders. Clinical assessment of the entire adult population (N= approximately 1,000) identified approximately 100 hypertensive subjects. For our study, we selected the individuals with the most-severe EH (i.e., diastolic blood pressure >100 mm Hg), belonging to a single deep-rooted pedigree (12 generations), whose common ancestors lived in the 17th century. We performed a three-stage genomewide search using 36 affected individuals, by means of parametric linkage and allele-sharing approaches. LOD scores >1 were observed on chromosomes 1, 2, 13, 15, 17, and 19 (stage I). The most striking result was found in a 7.57-cM region on chromosome 2p24-p25. All five nonparametric linkage statistics estimated by the SimWalk2 program lie above the significance threshold of P<.008 for the whole region. Similar significance was obtained for 2p24-25 when parametric linkage (LOD score 1.99) and linkage disequilibrium mapping (P=.00006) were used, suggesting that a hypertension-susceptibility locus is located between D2S2278 and D2S168. This finding is strengthened by a recent report of linkage with marker D2S168 in a hypertensive sib-pair sample from China.