Rare genomic structural variants in complex disease: lessons from the replication of associations with obesity.

The limited ability of common variants to account for the genetic contribution to complex disease has prompted searches for rare variants of large effect, to partly explain the 'missing heritability'. Analyses of genome-wide genotyping data have identified genomic structural variants (GSVs) as a source of such rare causal variants. Recent studies have reported multiple GSV loci associated with risk of obesity. We attempted to replicate these associations by similar analysis of two familial-obesity case-control cohorts and a population cohort, and detected GSVs at 11 out of 18 loci, at frequencies similar to those previously reported. Based on their reported frequencies and effect sizes (OR≥25), we had sufficient statistical power to detect the large majority (80%) of genuine associations at these loci. However, only one obesity association was replicated. Deletion of a 220 kb region on chromosome 16p11.2 has a carrier population frequency of 2×10(-4) (95% confidence interval [9.6×10(-5)-3.1×10(-4)]); accounts overall for 0.5% [0.19%-0.82%] of severe childhood obesity cases (P = 3.8×10(-10); odds ratio = 25.0 [9.9-60.6]); and results in a mean body mass index (BMI) increase of 5.8 kg.m(-2) [1.8-10.3] in adults from the general population. We also attempted replication using BMI as a quantitative trait in our population cohort; associations with BMI at or near nominal significance were detected at two further loci near KIF2B and within FOXP2, but these did not survive correction for multiple testing. These findings emphasise several issues of importance when conducting rare GSV association, including the need for careful cohort selection and replication strategy, accurate GSV identification, and appropriate correction for multiple testing and/or control of false discovery rate. Moreover, they highlight the potential difficulty in replicating rare CNV associations across different populations. Nevertheless, we show that such studies are potentially valuable for the identification of variants making an appreciable contribution to complex disease.

Effects of ciliary muscle plasmid electrotransfer of TNF-alpha soluble receptor variants in experimental uveitis.

Intraocular inflammation has been recognized as a major factor leading to blindness. Because tumor necrosis factor-alpha (TNF-alpha) enhances intraocular cytotoxic events, systemic anti-TNF therapies have been introduced in the treatment of severe intraocular inflammation, but frequent re-injections are needed and are associated with severe side effects. We have devised a local intraocular nonviral gene therapy to deliver effective and sustained anti-TNF therapy in inflamed eyes. In this study, we show that transfection of the ciliary muscle by plasmids encoding for three different variants of the p55 TNF-alpha soluble receptor, using electrotransfer, resulted in sustained intraocular secretion of the encoded proteins, without any detection in the serum. In the eye, even the shorter monomeric variant resulted in efficient neutralization of TNF-alpha in a rat experimental model of endotoxin-induced uveitis, as long as 3 months after transfection. A subsequent downregulation of interleukin (IL)-6 and iNOS and upregulation of IL-10 expression was observed together with a decreased rolling of inflammatory cells in anterior segment vessels and reduced infiltration within the ocular tissues. Our results indicate that using a nonviral gene therapy strategy, the local self-production of monomeric TNF-alpha soluble receptors induces a local immunomodulation enabling the control of intraocular inflammation.

The fine-scale architecture of structural variants in 17 mouse genomes.

BACKGROUND: Accurate catalogs of structural variants (SVs) in mammalian genomes are necessary to elucidate the potential mechanisms that drive SV formation and to assess their functional impact. Next generation sequencing methods for SV detection are an advance on array-based methods, but are almost exclusively limited to four basic types: deletions, insertions, inversions and copy number gains. RESULTS: By visual inspection of 100 Mbp of genome to which next generation sequence data from 17 inbred mouse strains had been aligned, we identify and interpret 21 paired-end mapping patterns, which we validate by PCR. These paired-end mapping patterns reveal a greater diversity and complexity in SVs than previously recognized. In addition, Sanger-based sequence analysis of 4,176 breakpoints at 261 SV sites reveal additional complexity at approximately a quarter of structural variants analyzed. We find micro-deletions and micro-insertions at SV breakpoints, ranging from 1 to 107 bp, and SNPs that extend breakpoint micro-homology and may catalyze SV formation. CONCLUSIONS: An integrative approach using experimental analyses to train computational SV calling is essential for the accurate resolution of the architecture of SVs. We find considerable complexity in SV formation; about a quarter of SVs in the mouse are composed of a complex mixture of deletion, insertion, inversion and copy number gain. Computational methods can be adapted to identify most paired-end mapping patterns.

Periphere okklusive Vaskulopathie bei einer Patientin mit kombinierter Präeklampsie und Alpha-Thalassämie minor [Peripheral occlusive retinal vasculopathy in a woman with combined preeclampsia and alpha thalassemia minor].

Normalerweise eine Störung der ersten Schwangerschaft, ist die Präeklampsie charakterisiert durch eine arterielle Hypertonie (> 140 mmHg systolisch oder > 90 mmHg diastolisch), die in der Regel nach der 20. Schwangerschaftswoche auftritt und von einer Proteinurie begleitet wird [1]. Die Präeklampsie wird als ,,schwer" bezeichnet, wenn sie mit einer wesentlichen Erhöhung des Blutdrucks (> 160 mmHg systolisch oder > 110 mmHg diastolisch), schwerer Proteinurie, Oligurie, Lungenödem, abdominalen Schmerzen, Leberfunktionsstörungen, Thrombozytopenie und visuellen oder zerebralen Symptomen einhergeht. Eine Eklampsie wiederum ist durch die Entwicklung von tonisch-klonischen Anfällen bei einer präeklamptischen Patientin charakterisiert. Bei der Alpha-Thalassämie tritt ein Defekt von 2 oder mehr der 4 Alpha-Globin-Gene auf. Von einer Alpha-Thalassämie minor spricht man, wenn 2 Alpha-Ketten-Gene deletiert sind. Sie tritt häufig bei Menschen aus Afrika, Südostasien, dem westindischen und mediterranen Raum auf. Die Alpha-Thalassämie minor verursacht eine milde bis moderate mikrozytäre Anämie. Wir berichten über eine Patientin mit peripherer okklusiver Vaskulopathie im Rahmen einer kombinierten Präeklampsie und Alpha-Thalassämie minor.

Identification and validation of copy number variants using SNP genotyping arrays from a large clinical cohort.

BACKGROUND: Genotypes obtained with commercial SNP arrays have been extensively used in many large case-control or population-based cohorts for SNP-based genome-wide association studies for a multitude of traits. Yet, these genotypes capture only a small fraction of the variance of the studied traits. Genomic structural variants (GSV) such as Copy Number Variation (CNV) may account for part of the missing heritability, but their comprehensive detection requires either next-generation arrays or sequencing. Sophisticated algorithms that infer CNVs by combining the intensities from SNP-probes for the two alleles can already be used to extract a partial view of such GSV from existing data sets. RESULTS: Here we present several advances to facilitate the latter approach. First, we introduce a novel CNV detection method based on a Gaussian Mixture Model. Second, we propose a new algorithm, PCA merge, for combining copy-number profiles from many individuals into consensus regions. We applied both our new methods as well as existing ones to data from 5612 individuals from the CoLaus study who were genotyped on Affymetrix 500K arrays. We developed a number of procedures in order to evaluate the performance of the different methods. This includes comparison with previously published CNVs as well as using a replication sample of 239 individuals, genotyped with Illumina 550K arrays. We also established a new evaluation procedure that employs the fact that related individuals are expected to share their CNVs more frequently than randomly selected individuals. The ability to detect both rare and common CNVs provides a valuable resource that will facilitate association studies exploring potential phenotypic associations with CNVs. CONCLUSION: Our new methodologies for CNV detection and their evaluation will help in extracting additional information from the large amount of SNP-genotyping data on various cohorts and use this to explore structural variants and their impact on complex traits.

Structural variation and its effect on expression.

Structural variation, whether it is caused by copy number variants or present in a balanced form, such as reciprocal translocations and inversions, can have a profound and dramatic effect on the expression of genes mapping within and close to the rearrangement, as well as affecting others genome wide. These effects can be caused by altering the copy number of one or more genes or regulatory elements (dosage effect) or from physical disruption of links between regulatory elements and their associated gene or genes, resulting in perturbation of expression. Similarly, large-scale structural variants can result in genome-wide expression changes by altering the positions that chromosomes occupy within the nucleus, potentially disrupting not only local cis interactions, but also trans interactions that occur throughout the genome. Structural variation is, therefore, a significant factor in the study of gene expression and is discussed here in more detail.

Copy number variation and chromatin structure

The functional consequences of structural variation in the human genome range from adaptation, to phenotypic variation, to predisposition to diseases. Copy number variation (CNV) was shown to influence the phenotype by modifying, in a somewhat dose-dependent manner, the expression of genes that map within them, as well as that of genes located on their flanks. To assess the possible mechanism(s) behind this neighboring effect, we compared histone modification status of cell lines from patients affected by Williams-Beuren, Williams-Beuren region duplication, Smith-Magenis or DiGeorge Syndrome and control individuals using a high-throughput version of chromatin immuno-precipitation method (ChIP), called ChlP-seq. We monitored monomethylation of lysine K20 on histone H4 and trimethylation of lysine K27 on histone H3, as proxies for open and condensed chromatin, respectively. Consistent with the changes in expression levels observed for multiple genes mapping on the entire length of chromosomes affected by structural variants, we also detected regions with modified histone status between samples, up- and downstream from the critical regions, up to the end of the rearranged chromosome. We also gauged the intrachromosomal interactions of these cell lines utilizing chromosome conformation capture (4C-seq) technique. We observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together, possibly forming an interacting cluster with each other and the WBSCR. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We conclude, that large genomic rearrangements can lead to changes in the state of the chromatin spreading far away from the critical region, thus possibly affecting expression globally and as a result modifying the phenotype of the patients. - Les conséquences fonctionnelles des variations structurelles dans le génome humain sont vastes, allant de l'adaptation, en passant par les variations phénotypiques, aux prédispositions à certaines maladies. Il a été démontré que les variations du nombre de copies (CNV) influencent le phénotype en modifiant, d'une manière plus ou moins dose-dépendante, l'expression des gènes se situant à l'intérieur de ces régions, mais également celle des gènes se trouvant dans les régions flanquantes. Afin d'étudier les mécanismes possibles sous-jacents à cet effet de voisinage, nous avons comparé les états de modification des histones dans des lignées cellulaires dérivées de patients atteints du syndrome de Williams-Beuren, de la duplication de la région Williams-Beuren, du syndrome de Smith-Magenis ou du syndrome de Di- George et d'individus contrôles en utilisant une version haut-débit de la méthode d'immunoprécipitation de la chromatine (ChIP), appelée ChIP-seq. Nous avons suivi la mono-méthylation de la lysine K20 sur l'histone H4 et la tri-méthylation de la lysine K27 sur l'histone H3, marqueurs respectifs de la chromatine ouverte et fermée. En accord avec les changements de niveaux d'expression observés pour de multiples gènes tout le long des chromosomes affectés par les CNVs, nous avons aussi détecté des régions présentant des modifications d'histones entre les échantillons, situées de part et d'autre des régions critiques, jusqu'aux extrémités du chromosome réarrangé. Nous avons aussi évalué les interactions intra-chromosomiques ayant lieu dans ces cellules par l'utilisation de la technique de capture de conformation des chromosomes (4C-seq). Nous avons observé qu'un groupe de gènes flanquants la région critique du syndrome de Williams-Beuren (WBSCR) forment souvent une boucle, constituant un groupe d'interactions privilégiées entre ces gènes et la WBSCR. La délétion de la WBSCR perturbe l'expression de ce groupe de gènes flanquants, mais également les interactions à grande échelle entre eux et la région réarrangée. Nous en concluons que les larges réarrangements génomiques peuvent aboutir à des changements de l'état de la chromatine pouvant s'étendre bien plus loin que la région critique, affectant donc potentiellement l'expression de manière globale et ainsi modifiant le phénotype des patients.

A comparative phenotypic and genomic analysis of C57BL/6J and C57BL/6N mouse strains.

BACKGROUND: The mouse inbred line C57BL/6J is widely used in mouse genetics and its genome has been incorporated into many genetic reference populations. More recently large initiatives such as the International Knockout Mouse Consortium (IKMC) are using the C57BL/6N mouse strain to generate null alleles for all mouse genes. Hence both strains are now widely used in mouse genetics studies. Here we perform a comprehensive genomic and phenotypic analysis of the two strains to identify differences that may influence their underlying genetic mechanisms. RESULTS: We undertake genome sequence comparisons of C57BL/6J and C57BL/6N to identify SNPs, indels and structural variants, with a focus on identifying all coding variants. We annotate 34 SNPs and 2 indels that distinguish C57BL/6J and C57BL/6N coding sequences, as well as 15 structural variants that overlap a gene. In parallel we assess the comparative phenotypes of the two inbred lines utilizing the EMPReSSslim phenotyping pipeline, a broad based assessment encompassing diverse biological systems. We perform additional secondary phenotyping assessments to explore other phenotype domains and to elaborate phenotype differences identified in the primary assessment. We uncover significant phenotypic differences between the two lines, replicated across multiple centers, in a number of physiological, biochemical and behavioral systems. CONCLUSIONS: Comparison of C57BL/6J and C57BL/6N demonstrates a range of phenotypic differences that have the potential to impact upon penetrance and expressivity of mutational effects in these strains. Moreover, the sequence variants we identify provide a set of candidate genes for the phenotypic differences observed between the two strains.

Copy number variation of KIR genes influences HIV-1 control.

A genome-wide screen for large structural variants showed that a copy number variant (CNV) in the region encoding killer cell immunoglobulin-like receptors (KIR) associates with HIV-1 control as measured by plasma viral load at set point in individuals of European ancestry. This CNV encompasses the KIR3DL1-KIR3DS1 locus, encoding receptors that interact with specific HLA-Bw4 molecules to regulate the activation of lymphocyte subsets including natural killer (NK) cells. We quantified the number of copies of KIR3DS1 and KIR3DL1 in a large HIV-1 positive cohort, and showed that an increase in KIR3DS1 count associates with a lower viral set point if its putative ligand is present (p = 0.00028), as does an increase in KIR3DL1 count in the presence of KIR3DS1 and appropriate ligands for both receptors (p = 0.0015). We further provide functional data that demonstrate that NK cells from individuals with multiple copies of KIR3DL1, in the presence of KIR3DS1 and the appropriate ligands, inhibit HIV-1 replication more robustly, and associated with a significant expansion in the frequency of KIR3DS1+, but not KIR3DL1+, NK cells in their peripheral blood. Our results suggest that the relative amounts of these activating and inhibitory KIR play a role in regulating the peripheral expansion of highly antiviral KIR3DS1+ NK cells, which may determine differences in HIV-1 control following infection.

Isoelectric focusing of alpha-1 acid glycoprotein (orosomucoid) in immobilized pH-gradients with 8M urea: detection of its desialylated variants using an alkaline phosphatase-linked secondary antibody system.

A method allowing a clear separation of the different variants of desialylated alpha 1-acid glycoprotein (orosomucoid) has been developed using isoelectric focusing in immobilized pH gradients, supplemented with 8 M urea and 2% v/v 2-mercaptoethanol. Immunoblotting with two antibody-steps afforded high sensitivity and permitted the detection of about 700 pg of alpha 1-acid glycoprotein in a 20 microL plasma sample diluted 1:28 672. A one year old bloodstrain, kept at room temperature, could easily be phenotyped.

Binding of d-methadone, l-methadone, and dl-methadone to proteins in plasma of healthy volunteers: role of the variants of alpha 1-acid glycoprotein.

The plasma concentrations of alpha 1-acid glycoprotein (AAG), albumin, triglycerides, cholesterol, and total proteins, as well as the plasma binding of racemic, d-methadone, and l-methadone were measured in 45 healthy subjects. The AAG phenotypes and the concentrations of AAG variants were also determined. The measured free fractions for racemic, d-methadone, and l-methadone were, respectively, 12.7% +/- 3.3%, 10.0% +/- 2.9%, and 14.2% +/- 3.2% (mean +/- SD). A significant correlation was obtained between the binding ratio (B/F) for dl-methadone and the total AAG concentration (r = 0.724; p less than 0.001). A multiple stepwise regression analysis showed that AAG was the main explanatory variable for the binding of the racemate. When concentrations of AAG variants were considered, a significant correlation was obtained between the binding ratio of dl-methadone and orosomucoid2 A concentration (r = 0.715; p less than 0.001), a weak correlation between dl-methadone and orosomucoid1 S concentration (r = 0.494; p less than 0.001), and no correlation between dl-methadone and orosomucoid1 F1 concentration (r = 0.049; not significant). Similar findings were obtained with the enantiomers. This study shows the importance of considering not only total AAG but also concentrations of AAG variants when measuring the binding of methadone and possibly of other drugs in plasma.

Binding of amitriptyline to alpha 1-acid glycoprotein and its variants.

Binding studies have been performed between amitriptyline and i) native alpha 1-acid glycoprotein (AAG); ii) its desialylated form; iii) its two variants, S-AAG and F-AAG; and iv) a mixture of S-AAG and F-AAG. Scatchard analysis revealed the presence of two classes of binding sites on AAG. For native AAG, the first class (of high affinity) has an association constant (Ka1) of 1.5 x 10(6) L mol-1 and a number of binding sites per mole of protein (n1) of 0.25, while the second class (of low affinity) has a Ka2 of 3.2 x 10(4) L mol-1 and a n2 of 0.94. Similar data were found for desialylated AAG. S-AAG and F-AAG do not differ in their association constants measured with amitriptyline, but in their number of binding sites per mole of protein (n): S-AAG: n1 = 0.56, n2 = 0.52; F-AAG: n1 = 0.17, n2 = 0.71. These results confirm those of a previous study, in which a higher affinity of S-AAG towards various basic drugs in comparison with F-AAG has been found.

Use of transgenic mice for the characterization of human alpha 1-acid glycoprotein (orosomucoid) variants.

Sera from transgenic mice (TM) carrying human genes of alpha 1-acid glycoprotein (orosomucoid or ORM) have been analyzed by isoelectrofocusing and subsequent immunoblotting with antihuman ORM antibodies. With this technique it is possible to reveal selectively the human protein secreted in the TM sera. Orosomucoid bands present in TM sera have been compared with those of the most common human ORM phenotypes to correlate the products of specific genes to previously identified genetic variants. In this paper, we report the identification of the genes encoding for variants ORM1 F1 and ORM2 A, which are genes AGP-A and AGP-B/B' respectively. The nucleotide sequences of these genes are known; therefore a direct correlation between variants and specific amino acid sequences can be established.

Selectivity in the binding of psychotropic drugs to the variants of alpha-1 acid glycoprotein.

The S- and F-forms of alpha-1 acid glycoprotein (AAG) variants have been isolated by isoelectric focusing with immobilines from commercially available AAG. In equilibrium dialysis experiments using a multicompartmental system, a higher affinity for various basic drugs has been found with S- in comparison with F-AAG: Amitriptyline, nortriptyline, imipramine, desipramine, trimipramine, methadone, thioridazine, clomipramine, desmethylclomipramine, and maprotiline. The selectivity (binding to S- vs. F-AAG) is the most pronounced for methadone and the lowest for thioridazine, while it is absent for the acidic drug mephenytoin.