Genome-wide analysis of rare copy number variations reveals PARK2 as a candidate gene for attention-deficit/hyperactivity disorder

Attention-deficit/hyperactivity disorder (ADHD) is a common, highly heritable neurodevelopmental disorder. Genetic loci have not yet been identified by genome-wide association studies. Rare copy number variations (CNVs), such as chromosomal deletions or duplications, have been implicated in ADHD and other neurodevelopmental disorders. To identify rare (frequency 1%) CNVs that increase the risk of ADHD, we performed a whole-genome CNV analysis based on 489 young ADHD patients and 1285 adult population-based controls and identified one significantly associated CNV region. In tests for a global burden of large (>500 kb) rare CNVs, we observed a nonsignificant (P=0.271) 1.126-fold enriched rate of subjects carrying at least one such CNV in the group of ADHD cases. Locus-specific tests of association were used to assess if there were more rare CNVs in cases compared with controls. Detected CNVs, which were significantly enriched in the ADHD group, were validated by quantitative (q)PCR. Findings were replicated in an independent sample of 386 young patients with ADHD and 781 young population-based healthy controls. We identified rare CNVs within the parkinson protein 2 gene (PARK2) with a significantly higher prevalence in ADHD patients than in controls (P=2.8 × 10(-4) after empirical correction for genome-wide testing). In total, the PARK2 locus (chr 6: 162 659 756-162 767 019) harboured three deletions and nine duplications in the ADHD patients and two deletions and two duplications in the controls. By qPCR analysis, we validated 11 of the 12 CNVs in ADHD patients (P=1.2 × 10(-3) after empirical correction for genome-wide testing). In the replication sample, CNVs at the PARK2 locus were found in four additional ADHD patients and one additional control (P=4.3 × 10(-2)). Our results suggest that copy number variants at the PARK2 locus contribute to the genetic susceptibility of ADHD. Mutations and CNVs in PARK2 are known to be associated with Parkinson disease.Molecular Psychiatry advance online publication, 20 November 2012; doi:10.1038/mp.2012.161.

Non-coding keratin variants associate with liver fibrosis progression in patients with hemochromatosis

Background Keratins 8 and 18 (K8/K18) are intermediate filament proteins that protect the liver from various forms of injury. Exonic K8/K18 variants associate with adverse outcome in acute liver failure and with liver fibrosis progression in patients with chronic hepatitis C infection or primary biliary cirrhosis. Given the association of K8/K18 variants with end-stage liver disease and progression in several chronic liver disorders, we studied the importance of keratin variants in patients with hemochromatosis. Methods The entire K8/K18 exonic regions were analyzed in 162 hemochromatosis patients carrying homozygous C282Y HFE (hemochromatosis gene) mutations. 234 liver-healthy subjects were used as controls. Exonic regions were PCR-amplified and analyzed using denaturing high-performance liquid chromatography and DNA sequencing. Previously-generated transgenic mice overexpressing K8 G62C were studied for their susceptibility to iron overload. Susceptibility to iron toxicity of primary hepatocytes that express K8 wild-type and G62C was also assessed. Results We identified amino-acid-altering keratin heterozygous variants in 10 of 162 hemochromatosis patients (6.2%) and non-coding heterozygous variants in 6 additional patients (3.7%). Two novel K8 variants (Q169E/R275W) were found. K8 R341H was the most common amino-acid altering variant (4 patients), and exclusively associated with an intronic KRT8 IVS7+10delC deletion. Intronic, but not amino-acid-altering variants associated with the development of liver fibrosis. In mice, or ex vivo, the K8 G62C variant did not affect iron-accumulation in response to iron-rich diet or the extent of iron-induced hepatocellular injury. Conclusion In patients with hemochromatosis, intronic but not exonic K8/K18 variants associate with liver fibrosis development.

The Impact of PPARγ Genetic Variants on IBD Susceptibility and IBD Disease Course

PPARγ is a nuclear receptor that regulates numerous pathways including cytokine expression and immune responses and plays an important role in controlling colon inflammation. We aimed at determining the occurring PPARγ SNPs, at predicting the haplotypes, and at determining the frequency outcome in inflammatory bowel disease (IBD) patients in comparison with healthy controls. We determined genetic variants in the coding exons and flanking intronic sequences of the NR1C3 gene in 284 IBD patients and 194 controls and predicted NR1C3 haplotypes via bioinformatic analysis. We investigated whether certain NR1C3 variants are associated with susceptibility to IBD or its disease course. None of the detected 22 NR1C3 variants were associated with IBD. Two variants with allelic frequencies over 1% were included in haplotype/diplotype analyses. None of the NR3C1 haplotypes showed association with IBD development or disease course. We conclude that NR1C3 haplotypes are not related to IBD susceptibility or IBD disease activity.

Molecular characterization of the porcine deleted in malignant brain tumors 1 gene (DMBT1)

The human gene deleted in malignant brain tumors 1 (DMBT1) is considered to play a role in tumorigenesis and pathogen defense. It encodes a protein with multiple scavenger receptor cysteine-rich (SRCR) domains, which are involved in recognition and binding of a broad spectrum of bacterial pathogens. The SRCR domains are encoded by highly homologous repetitive exons, whose number in humans may vary from 8 to 13 due to genetic polymorphism. Here, we characterized the porcine DMBT1 gene on the mRNA and genomic level. We assembled a 4.5 kb porcine DMBT1 cDNA sequence from RT-PCR amplified seminal vesicle RNA. The porcine DMBT1 cDNA contains an open reading frame of 4050 nt. The transcript gives rise to a putative polypeptide of 1349 amino acids with a calculated mass of 147.9 kDa. Compared to human DMBT1, it contains only four N-terminal SRCR domains. Northern blotting revealed transcripts of approximately 4.7 kb in size in the tissues analyzed. Analysis of ESTs suggested the existence of secreted and transmembrane variants. The porcine DMBT1 gene spans about 54 kb on chromosome 14q28-q29. In contrast to the characterized cDNA, the genomic BAC clone only contained 3 exons coding for N-terminal SRCR domains. In different mammalian DMBT1 orthologs large interspecific differences in the number of SRCR exons and utilization of the transmembrane exon exist. Our data suggest that the porcine DMBT1 gene may share with the human DMBT1 gene additional intraspecific variations in the number of SRCR-coding exons.

Novel mutations in the human sucrase-isomaltase gene (SI) that cause congenital carbohydrate malabsorption

Disaccharide intolerance I or congenital sucrase-isomaltase deficiency (CSID) is a disorder leading to maldigestion of disaccharides, which is autosomal recessively inherited. Here we analyzed the sucrase-isomaltase (SI) gene from 11 patients of Hungarian origin with congenital sucrase-isomaltase deficiency. Variants in the SI gene had previously been described in CSID patients, which cause amino acid exchanges that affect the transport, the processing, or the function of the SI protein. None of our patients had known mutations for CSID. Our analyses revealed 43 SI variants in total, 15 within exons and one at a splice site. Eight of the exonic mutations lead to amino acid exchanges, causing hypomorph or null alleles. One new variation affects a splice site, which is also predicted to result in a null allele. All potential pathological alterations were present on one allele only. In six out of the 11 patients the phenotype of CSID could be explained by compound heterozygosity.

Characterization of the canine CLCN3 gene and evaluation as candidate for late-onset NCL

BACKGROUND: The neuronal ceroid lipofuscinoses (NCL) are a heterogenous group of inherited progressive neurodegenerative diseases in different mammalian species. Tibetan Terrier and Polish Owczarek Nizinny (PON) dogs show rare late-onset NCL variants with autosomal recessive inheritance, which can not be explained by mutations of known human NCL genes. These dog breeds represent animal models for human late-onset NCL. In mice the chloride channel 3 gene (Clcn3) encoding an intracellular chloride channel was described to cause a phenotype similar to NCL. RESULTS: Two full-length cDNA splice variants of the canine CLCN3 gene are reported. The current canine whole genome sequence assembly was used for gene structure analyses and revealed 13 coding CLCN3 exons in 52 kb of genomic sequence. Sequence analysis of the coding exons and flanking intron regions of CLCN3 using six NCL-affected Tibetan terrier dogs and an NCL-affected Polish Owczarek Nizinny (PON) dog, as well as eight healthy Tibetan terrier dogs revealed 13 SNPs. No consistent CLCN3 haplotype was associated with NCL. CONCLUSION: For the examined animals we excluded the complete coding region and adjacent intronic regions of canine CLCN3 to harbor disease-causing mutations. Therefore it seems to be unlikely that a mutation in this gene is responsible for the late-onset NCL phenotype in these two dog breeds.

Isolated autosomal dominant growth hormone deficiency: stimulating mutant GH-1 gene expression drives GH-1 splice-site selection, cell proliferation, and apoptosis

The majority of mutations that cause isolated GH deficiency type II (IGHD II) affect splicing of GH-1 transcripts and produce a dominant-negative GH isoform lacking exon 3 resulting in a 17.5-kDa isoform, which further leads to disruption of the GH secretory pathway. A clinical variability in the severity of the IGHD II phenotype depending on the GH-1 gene alteration has been reported, and in vitro and transgenic animal data suggest that the onset and severity of the phenotype relates to the proportion of 17.5-kDa produced. The removal of GH in IGHD creates a positive feedback loop driving more GH expression, which may itself increase 17.5-kDa isoform productions from alternate splice sites in the mutated GH-1 allele. In this study, we aimed to test this idea by comparing the impact of stimulated expression by glucocorticoids on the production of different GH isoforms from wild-type (wt) and mutant GH-1 genes, relying on the glucocorticoid regulatory element within intron 1 in the GH-1 gene. AtT-20 cells were transfected with wt-GH or mutated GH-1 variants (5'IVS-3 + 2-bp T->C; 5'IVS-3 + 6 bp T->C; ISEm1: IVS-3 + 28 G->A) known to cause clinical IGHD II of varying severity. Cells were stimulated with 1 and 10 mum dexamethasone (DEX) for 24 h, after which the relative amounts of GH-1 splice variants were determined by semiquantitative and quantitative (TaqMan) RT-PCR. In the absence of DEX, only around 1% wt-GH-1 transcripts were the 17.5-kDa isoform, whereas the three mutant GH-1 variants produced 29, 39, and 78% of the 17.5-kDa isoform. DEX stimulated total GH-1 gene transcription from all constructs. Notably, however, DEX increased the amount of 17.5-kDa GH isoform relative to the 22- and 20-kDa isoforms produced from the mutated GH-1 variants, but not from wt-GH-1. This DEX-induced enhancement of 17.5-kDa GH isoform production, up to 100% in the most severe case, was completely blocked by the addition of RU486. In other studies, we measured cell proliferation rates, annexin V staining, and DNA fragmentation in cells transfected with the same GH-1 constructs. The results showed that that the 5'IVS-3 + 2-bp GH-1 gene mutation had a more severe impact on those measures than the splice site mutations within 5'IVS-3 + 6 bp or ISE +28, in line with the clinical severity observed with these mutations. Our findings that the proportion of 17.5-kDa produced from mutant GH-1 alleles increases with increased drive for gene expression may help to explain the variable onset progression, and severity observed in IGHD II.

Short communication: duplication in the 5'-flanking region of the beta-lactoglobulin gene is linked to the BLG A allele

beta-Lactoglobulin (beta-LG) is the major whey protein in the milk of cows and other ruminants. It is well established that the predominant genetic variants beta-LG A and B are differentially expressed. Extensive investigation of the genetic variation in the promoter region of the BLG gene revealed the existence of specific haplotypes associated with the A and B variants. However, the genetic basis for the differentially expressed BLG A and B alleles is still elusive. In this study additional genetic variation further upstream in the 5'-flanking region of the BLG gene was identified, including 6 single nucleotide substitutions, a single nucleotide deletion, and a 7-bp duplication. Comparison of DNA sequences showed that the investigated 5'-flanking region is highly conserved between ruminants, and the duplication g.-1885_-1879dupCTCTCGC and the substitution g.-1888A>G are only found in the BLG A and D alleles in cattle. The cytosine at position g.-1957 and the thymines at positions g.-2008 and g.-2049 are only found in BLG B alleles of cattle. It is suggested that the described genetic variability contributes to the differential allelic expression of the BLG gene.

Identification of polymorphisms in the human 11beta-hydroxysteroid dehydrogenase type 2 gene promoter: functional characterization and relevance for salt sensitivity

Reduced activity of 11beta-hydroxysteroid dehydrogenase type 2 (11beta-HSD2) plays a role in essential hypertension and the sensitivity of blood pressure to dietary salt. Nonconservative mutations in the coding region are extremely rare and do not explain the variable 11beta-HSD2 activity. We focused therefore on the 5'-regulatory region and identified and characterized the first promoter polymorphisms. Transfections of variants G-209A and G-126A into SW620 cells reduced promoter activity and affinity for activators nuclear factor 1 (NF1) and Sp1. Chromatin immunoprecipitation revealed Sp1, NF1, and glucocorticoid receptor (GR) binding to the HSD11B2 promoter. Dexamethasone induced expression of mRNA and activity of HSD11B2. GR and/or NF1 overexpression increased endogenous HSD11B2 mRNA and activity. GR complexes cooperated with NF1 to activate HSD11B2, an effect diminished in the presence of the G-209A variant. When compared to salt-resistant subjects (96), salt-sensitive volunteers (54) more frequently had the G-209A variant, higher occurrence of alleles A4/A7 of polymorphic microsatellite marker, and higher urinary ratios of cortisol to cortisone metabolites. First, we conclude that the mechanism of glucocorticoid-induced HSD11B2 expression is mainly mediated by cooperation between GR and NF1 on the HSD11B2 promoter and, second, that the newly identified promoter variants reduce activity and cooperation of cognate transcription factors, resulting in diminished HSD11B2 transcription, an effect favoring salt sensitivity.

Gene duplication: a drive for phenotypic diversity and cause of human disease

Gene duplication is one of the key factors driving genetic innovation, i.e., producing novel genetic variants. Although the contribution of whole-genome and segmental duplications to phenotypic diversity across species is widely appreciated, the phenotypic spectrum and potential pathogenicity of small-scale duplications in individual genomes are less well explored. This review discusses the nature of small-scale duplications and the phenotypes produced by such duplications. Phenotypic variation and disease phenotypes induced by duplications are more diverse and widespread than previously anticipated, and duplications are a major class of disease-related genomic variation. Pathogenic duplications particularly involve dosage-sensitive genes with both similar and dissimilar over- and underexpression phenotypes, and genes encoding proteins with a propensity to aggregate. Phenotypes related to human-specific copy number variation in genes regulating environmental responses and immunity are increasingly recognized. Small genomic duplications containing defense-related genes also contribute to complex common phenotypes.

High variability and non-neutral evolution of the mammalian avpr1a gene

BACKGROUND: The arginine-vasopressin 1a receptor has been identified as a key determinant for social behaviour in Microtus voles, humans and other mammals. Nevertheless, the genetic bases of complex phenotypic traits like differences in social and mating behaviour among species and individuals remain largely unknown. Contrary to previous studies focusing on differences in the promotor region of the gene, we investigate here the level of functional variation in the coding region (exon 1) of this locus. RESULTS: We detected high sequence diversity between higher mammalian taxa as well as between species of the genus Microtus. This includes length variation and radical amino acid changes, as well as the presence of distinct protein variants within individuals. Additionally, negative selection prevails on most parts of the first exon of the arginine-vasopressin receptor 1a (avpr1a) gene but it contains regions with higher rates of change that harbour positively selected sites. Synonymous and non-synonymous substitution rates in the avpr1a gene are not exceptional compared to other genes, but they exceed those found in related hormone receptors with similar functions. DISCUSSION: These results stress the importance of considering variation in the coding sequence of avpr1a in regards to associations with life history traits (e.g. social behaviour, mating system, habitat requirements) of voles, other mammals and humans in particular.

Differential expression of stx2 variants in Shiga toxin-producing Escherichia coli belonging to seropathotypes A and C

Only a subset of Shiga toxin (Stx)-producing Escherichia coli (STEC) are human pathogens, but the characteristics that account for differences in pathogenicity are not well understood. In this study, we investigated the distribution of the stx variants coding for Stx2 and its variants in highly virulent STEC of seropathotype A and low-pathogenic STEC of seropathotype C. We analysed and compared transcription of the corresponding genes, production of Shiga toxins, and stx-phage release in basal as well as in induced conditions. We found that the stx(2) variant was mainly associated with strains of seropathotype A, whereas most of the strains of seropathotype C possessed the stx(2-vhb) variant, which was frequently associated with stx(2), stx(2-vha) or stx(2c). Levels of stx(2) and stx(2)-related mRNA were higher in strains belonging to seropathotype A and in those strains of seropathotype C that express the stx(2) variant than in the remaining strains of seropathotype C. The stx(2-vhb) genes were the least expressed, in basal as well as in induced conditions, and in many cases did not seem to be carried by an inducible prophage. A clear correlation was observed between stx mRNA levels and stx-phage DNA in the culture supernatants, suggesting that most stx(2)-related genes are expressed only when they are carried by a phage. In conclusion, some relationship between stx(2)-related gene expression in vitro and the seropathotype of the STEC strains was observed. A higher expression of the stx(2) gene and a higher release of its product, in basal as well as in induced conditions, was observed in pathogenic strains of seropathotype A. A subset of strains of seropathotype C shows the same characteristics and could be a high risk to human health.

Do common genetic variants in endotoxin signaling pathway contribute to predisposition to alcoholic liver cirrhosis?

BACKGROUND: Tumor necrosis factor-alpha (TNF-alpha) and interleukin-1beta (IL-1beta), produced by endotoxin-activated Kupffer cells, play a key role in the pathogenesis of alcoholic liver cirrhosis (ALC). Alleles TNFA -238A, IL1B -31T and variant IL1RN*2 of repeat polymorphism in the gene encoding the IL-1 receptor antagonist increase production of TNF-alpha and IL-1beta, respectively. Alleles CD14 -159T, TLR4 c.896G and TLR4 c.1196T modify activation of Kupffer cells by endotoxin. We confirmed the published associations between these common variants and genetic predisposition to ALC by means of a large case-control association study conducted on two Central European populations. METHODS: The study population comprised a Czech sample of 198 ALC patients and 370 controls (MONICA project), and a German sample of 173 ALC patients and 331 controls (KORA-Augsburg), and 109 heavy drinkers without liver disease. RESULTS: Single locus analysis revealed no significant difference between patients and controls in all tested loci. Diplotype [IL1RN 2/ 2; IL1B -31T+] was associated with increased risk of ALC in the pilot study, but not in the validation samples. CONCLUSIONS: Although cytokine mediated immune reactions play a role in the pathogenesis of ALC, hereditary susceptibility caused by variants in the corresponding genes is low in Central European populations.

High aldosterone-to-renin variants of CYP11B2 and pregnancy outcome

BACKGROUND: Increased aldosterone concentrations and volume expansion of normal pregnancies are hallmarks of normal pregnancies and blunted in pre-eclampsia. Accordingly, we hypothesized an active mineralocorticoid system to protect from pre-eclampsia. METHODS: In pregnant women (normotensive n = 44; pre-eclamptic n = 48), blood pressure, urinary tetrahydro-aldosterone excretion and activating polymorphisms (SF-1 site and intron 2) of the aldosterone synthase gene (CYP11B2) were determined; 185 non-pregnant normotensive individuals served as control. Amino acid-changing polymorphisms of the DNA- and agonist-binding regions of the mineralocorticoid receptor were evaluated by RT-PCR, SSCP and sequencing. RESULTS: Urinary tetrahydro-aldosterone excretion was reduced in pre-eclampsia as compared to normal pregnancy (P < 0.05). It inversely correlated with blood pressure (r = 0.99, P < 0.04). Homozygosity for activating CYP11B2 polymorphisms was preferably present in normotensive as compared to pre-eclamptic pregnancies, identified (intron 2, P = 0.005; SF-1 site, P = 0.016). Two mutant haplotypes decreased the risk of developing pre-eclampsia (RR 0.16; CI 0.05-0.54; P < 0.001). In contrast, intron 2 wild type predisposed to pre-eclampsia (P < 0.0015). No functional mineralocorticoid receptor mutant has been observed. CONCLUSIONS: High aldosterone availability is associated with lower maternal blood pressure. In line with this observation, gain-of-function variants of the CYP11B2 reduce the risk of developing pre-eclampsia. Mutants of the mineralocorticoid receptor cannot explain the frequent syndrome of pre-eclampsia.

Bacterial phenotype variants in group B streptococcal toxic shock syndrome

We conducted genetic and functional analyses of isolates from a patient with group B streptococcal (GBS) necrotizing fasciitis and toxic shock syndrome. Tissue cultures simultaneously showed colonies with high hemolysis (HH) and low hemolysis (LH). Conversely, the HH and LH variants exhibited low capsule (LC) and high capsule (HC) expression, respectively. Molecular analysis demonstrated that the 2 GBS variants were of the same clonal origin. Genetic analysis found a 3-bp deletion in the covR gene of the HH/LC variant. Functionally, this isolate was associated with an increased growth rate in vitro and with higher interleukin-8 induction. However, in whole blood, opsonophagocytic and intracellular killing assays, the LH/HC phenotype demonstrated higher resistance to host phagocytic killing. In a murine model, LH/HC resulted in higher levels of bacteremia and increased host mortality rate. These findings demonstrate differences in GBS isolates of the same clonal origin but varying phenotypes.