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 (P202;=202;3.8×10(-10); odds ratio202;=202;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.

The 16p11.2 locus modulates brain structures common to autism, schizophrenia and obesity

Anatomical structures and mechanisms linking genes to neuropsychiatric disorders are not deciphered. Reciprocal copy number variants at the 16p11.2 BP4-BP5 locus offer a unique opportunity to study the intermediate phenotypes in carriers at high risk for autism spectrum disorder (ASD) or schizophrenia (SZ). We investigated the variation in brain anatomy in 16p11.2 deletion and duplication carriers. Beyond gene dosage effects on global brain metrics, we show that the number of genomic copies negatively correlated to the gray matter volume and white matter tissue properties in cortico-subcortical regions implicated in reward, language and social cognition. Despite the near absence of ASD or SZ diagnoses in our 16p11.2 cohort, the pattern of brain anatomy changes in carriers spatially overlaps with the well-established structural abnormalities in ASD and SZ. Using measures of peripheral mRNA levels, we confirm our genomic copy number findings. This combined molecular, neuroimaging and clinical approach, applied to larger datasets, will help interpret the relative contributions of genes to neuropsychiatric conditions by measuring their effect on local brain anatomy.Molecular Psychiatry advance online publication, 25 November 2014; doi:10.1038/mp.2014.145.

Physical activity attenuates the influence of FTO variants on obesity risk: a meta-analysis of 218,166 adults and 19,268 children.

BACKGROUND: The FTO gene harbors the strongest known susceptibility locus for obesity. While many individual studies have suggested that physical activity (PA) may attenuate the effect of FTO on obesity risk, other studies have not been able to confirm this interaction. To confirm or refute unambiguously whether PA attenuates the association of FTO with obesity risk, we meta-analyzed data from 45 studies of adults (n202;=202;218,166) and nine studies of children and adolescents (n202;=202;19,268). METHODS AND FINDINGS: All studies identified to have data on the FTO rs9939609 variant (or any proxy [r(2)>0.8]) and PA were invited to participate, regardless of ethnicity or age of the participants. PA was standardized by categorizing it into a dichotomous variable (physically inactive versus active) in each study. Overall, 25% of adults and 13% of children were categorized as inactive. Interaction analyses were performed within each study by including the FTO×PA interaction term in an additive model, adjusting for age and sex. Subsequently, random effects meta-analysis was used to pool the interaction terms. In adults, the minor (A-) allele of rs9939609 increased the odds of obesity by 1.23-fold/allele (95% CI 1.20-1.26), but PA attenuated this effect (p(interaction) 202;=202;0.001). More specifically, the minor allele of rs9939609 increased the odds of obesity less in the physically active group (odds ratio 202;=202;1.22/allele, 95% CI 1.19-1.25) than in the inactive group (odds ratio 202;=202;1.30/allele, 95% CI 1.24-1.36). No such interaction was found in children and adolescents. CONCLUSIONS: The association of the FTO risk allele with the odds of obesity is attenuated by 27% in physically active adults, highlighting the importance of PA in particular in those genetically predisposed to obesity.

Immunochip SNP array identifies novel genetic variants conferring susceptibility to candidaemia.

Candidaemia is the fourth most common cause of bloodstream infection, with a high mortality rate of up to 40%. Identification of host genetic factors that confer susceptibility to candidaemia may aid in designing adjunctive immunotherapeutic strategies. Here we hypothesize that variation in immune genes may predispose to candidaemia. We analyse 118,989 single-nucleotide polymorphisms (SNPs) across 186 loci known to be associated with immune-mediated diseases in the largest candidaemia cohort to date of 217 patients of European ancestry and a group of 11,920 controls. We validate the significant associations by comparison with a disease-matched control group. We observe significant association between candidaemia and SNPs in the CD58 (P = 1.97 × 10(-11); odds ratio (OR) = 4.68), LCE4A-C1orf68 (P = 1.98 × 10(-10); OR = 4.25) and TAGAP (P = 1.84 × 10(-8); OR = 2.96) loci. Individuals carrying two or more risk alleles have an increased risk for candidaemia of 19.4-fold compared with individuals carrying no risk allele. We identify three novel genetic risk factors for candidaemia, which we subsequently validate for their role in antifungal host defence.

Disease-causing mutations improving the branch site and polypyrimidine tract: pseudoexon activation of LINE-2 and antisense Alu lacking the poly(T)-tail.

Cryptic exons or pseudoexons are typically activated by point mutations that create GT or AG dinucleotides of new 5' or 3' splice sites in introns, often in repetitive elements. Here we describe two cases of tetrahydrobiopterin deficiency caused by mutations improving the branch point sequence and polypyrimidine tracts of repeat-containing pseudoexons in the PTS gene. In the first case, we demonstrate a novel pathway of antisense Alu exonization, resulting from an intronic deletion that removed the poly(T)-tail of antisense AluSq. The deletion brought a favorable branch point sequence within proximity of the pseudoexon 3' splice site and removed an upstream AG dinucleotide required for the 3' splice site repression on normal alleles. New Alu exons can thus arise in the absence of poly(T)-tails that facilitated inclusion of most transposed elements in mRNAs by serving as polypyrimidine tracts, highlighting extraordinary flexibility of Alu repeats in shaping intron-exon structure. In the other case, a PTS pseudoexon was activated by an A>T substitution 9 nt upstream of its 3' splice site in a LINE-2 sequence, providing the first example of a disease-causing exonization of the most ancient interspersed repeat. These observations expand the spectrum of mutational mechanisms that introduce repetitive sequences in mature transcripts and illustrate the importance of intronic mutations in alternative splicing and phenotypic variability of hereditary disorders.

Caffeine intake and CYP1A2 variants associated with high caffeine intake protect non-smokers from hypertension.

The 15q24.1 locus, including CYP1A2, is associated with blood pressure (BP). The CYP1A2 rs762551 C allele is associated with lower CYP1A2 enzyme activity. CYP1A2 metabolizes caffeine and is induced by smoking. The association of caffeine consumption with hypertension remains controversial. We explored the effects of CYP1A2 variants and CYP1A2 enzyme activity on BP, focusing on caffeine as the potential mediator of CYP1A2 effects. Four observational (n = 16 719) and one quasi-experimental studies (n = 106) including European adults were conducted. Outcome measures were BP, caffeine intake, CYP1A2 activity and polymorphisms rs762551, rs1133323 and rs1378942. CYP1A2 variants were associated with hypertension in non-smokers, but not in smokers (CYP1A2-smoking interaction P = 0.01). Odds ratios (95% CIs) for hypertension for rs762551 CC, CA and AA genotypes were 1 (reference), 0.78 (0.59-1.02) and 0.66 (0.50-0.86), respectively, P = 0.004. Results were similar for the other variants. Higher CYP1A2 activity was linearly associated with lower BP after quitting smoking (P = 0.049 and P = 0.02 for systolic and diastolic BP, respectively), but not while smoking. In non-smokers, the CYP1A2 variants were associated with higher reported caffeine intake, which in turn was associated with lower odds of hypertension and lower BP (P = 0.01). In Mendelian randomization analyses using rs1133323 as instrument, each cup of caffeinated beverage was negatively associated with systolic BP [-9.57 (-16.22, -2.91) mmHg]. The associations of CYP1A2 variants with BP were modified by reported caffeine intake. These observational and quasi-experimental results strongly support a causal role of CYP1A2 in BP control via caffeine intake.

Peroxisomal Delta(3),Delta(2)-enoyl CoA isomerases and evolution of cytosolic paralogues in embryophytes

Delta(3),Delta(2)-enoyl CoA isomerase (ECI) is an enzyme that participates in the degradation of unsaturated fatty acids through the beta-oxidation cycle. Three genes encoding Delta(3),Delta(2)-enoyl CoA isomerases and named AtECI1, AtECI2 and AtECI3 have been identified in Arabidopsis thaliana. When expressed heterologously in Saccharomyces cerevisiae, all three ECI proteins were targeted to the peroxisomes and enabled the yeast Deltaeci1 mutant to degrade 10Z-heptadecenoic acid, demonstrating Delta(3),Delta(2)-enoyl CoA isomerase activity in vivo. Fusion proteins between yellow fluorescent protein and AtECI1 or AtECI2 were targeted to the peroxisomes in onion epidermal cells and Arabidopsis root cells, but a similar fusion protein with AtECI3 remained in the cytosol for both tissues. AtECI3 targeting to peroxisomes in S. cerevisiae was dependent on yeast PEX5, while expression of Arabidopsis PEX5 in yeast failed to target AtECI3 to peroxisomes. AtECI2 and AtECI3 are tandem duplicated genes and show a high level of amino acid conservation, except at the C-terminus; AtECI2 ends with the well conserved peroxisome targeting signal 1 (PTS1) terminal tripeptide PKL, while AtECI3 possesses a divergent HNL terminal tripeptide. Evolutionary analysis of ECI genes in plants revealed several independent duplication events, with duplications occurring in rice and Medicago truncatula, generating homologues with divergent C-termini and no recognizable PTS1. All plant ECI genes analyzed, including AtECI3, are under negative purifying selection, implying functionality of the cytosolic AtECI3. Analysis of the mammalian and fungal genomes failed to identify cytosolic variants of the Delta(3),Delta(2)-enoyl CoA isomerase, indicating that evolution of cytosolic Delta(3),Delta(2)-enoyl CoA isomerases is restricted to the plant kingdom

Orosomucoid system: 17 additional orosomucoid variants and proposal for a new nomenclature.

There are two forms of orosomucoid (ORM) in the sera of most individuals. They are encoded by two separate but closely linked loci, ORM1 and ORM2. A number of variants have been identified in various populations. Duplication and nonexpression are also observed in some populations. Thus, the ORM system is very complicated and its nomenclature is very confusing. In order to propose a new nomenclature, ORM variants detected by several laboratories have been compared and characterized by isoelectric focusing (IEF) followed by immunoprinting. A total of 57 different alleles including 17 new ones were identified. The 27 alleles were assigned to the ORM1 locus, and the others to the ORM2 locus. The designations ORM*F1, ORM1*F2, ORM1*S and ORM2*M were adopted for the four common alleles instead of ORM1*1, ORM1*3, ORM1*2 and ORM2*1 (ORM2*A), respectively. The variants were designated alpha numerically according to their relative mobilities after IEF in a pH gradient of 4.5-5.4 with Triton X-100 and glycerol. For the duplicated genes a prefix is added to a combined name of two alleles, e.g. ORM1*dB9S. Silent alleles were named ORM1*Q0 and ORM2*Q0 conventionally. In addition, the effects of diseases to ORM band patterns after IEF are also discussed.

Les variants al·lomòrfiques de la preposició 'a' en algunes varietats del català: descripció i anàlisi

Aquest treball ofereix, en primer lloc, una descripció lingüística de les formes que presenta la preposició a en quatre varietats del català (barceloní, varietats occidentals del català central, balear i tortosí). En segon lloc, presenta una anàlisi de les dades en què es demostra que la llengua resol els contactes vocàlics on intervé la preposició a mitjançant mecanismes que s'escapen del comportament fonològic general, és a dir processos d'elisió o fusió, això és mitjançant les formes /an/ o /am(b)/. En aquest sentit, es tracten aquestes formes com un cas d'al·lomorfia externa amb ordenació lèxica dels al·lomorfs, en el qual tenen un paper rellevant els conceptes d'especificitat i de pressió paradigmàtica. Finalment, es fa un tractament del fenomen en termes de la teoria de l'optimitat

Relevance of biallelic versus monoallelic TNFRSF13B mutations in distinguishing disease-causing from risk-increasing TNFRSF13B variants in antibody deficiency syndromes.

TNFRSF13B encodes transmembrane activator and calcium modulator and cyclophilin ligand interactor (TACI), a B cell- specific tumor necrosis factor (TNF) receptor superfamily member. Both biallelic and monoallelic TNFRSF13B mutations were identified in patients with common variable immunodeficiency disorders. The genetic complexity and variable clinical presentation of TACI deficiency prompted us to evaluate the genetic, immunologic, and clinical condition in 50 individuals with TNFRSF13B alterations, following screening of 564 unrelated patients with hypogammaglobulinemia. We identified 13 new sequence variants. The most frequent TNFRSF13B variants (C104R and A181E; n=39; 6.9%) were also present in a heterozygous state in 2% of 675 controls. All patients with biallelic mutations had hypogammaglobulinemia and nearly all showed impaired binding to a proliferation-inducing ligand (APRIL). However, the majority (n=41; 82%) of the pa-tients carried monoallelic changes in TNFRSF13B. Presence of a heterozygous mutation was associated with antibody deficiency (P< .001, relative risk 3.6). Heterozygosity for the most common mutation, C104R, was associated with disease (P< .001, relative risk 4.2). Furthermore, heterozygosity for C104R was associated with low numbers of IgD(-)CD27(+) B cells (P= .019), benign lymphoproliferation (P< .001), and autoimmune complications (P= .001). These associations indicate that C104R heterozygosity increases the risk for common variable immunodeficiency disorders and influences clinical presentation.

A 600 kb deletion syndrome at 16p11.2 leads to energy imbalance and neuropsychiatric disorders.

BACKGROUND: The recurrent ~600 kb 16p11.2 BP4-BP5 deletion is among the most frequent known genetic aetiologies of autism spectrum disorder (ASD) and related neurodevelopmental disorders. OBJECTIVE: To define the medical, neuropsychological, and behavioural phenotypes in carriers of this deletion. METHODS: We collected clinical data on 285 deletion carriers and performed detailed evaluations on 72 carriers and 68 intrafamilial non-carrier controls. RESULTS: When compared to intrafamilial controls, full scale intelligence quotient (FSIQ) is two standard deviations lower in carriers, and there is no difference between carriers referred for neurodevelopmental disorders and carriers identified through cascade family testing. Verbal IQ (mean 74) is lower than non-verbal IQ (mean 83) and a majority of carriers require speech therapy. Over 80% of individuals exhibit psychiatric disorders including ASD, which is present in 15% of the paediatric carriers. Increase in head circumference (HC) during infancy is similar to the HC and brain growth patterns observed in idiopathic ASD. Obesity, a major comorbidity present in 50% of the carriers by the age of 7 years, does not correlate with FSIQ or any behavioural trait. Seizures are present in 24% of carriers and occur independently of other symptoms. Malformations are infrequently found, confirming only a few of the previously reported associations. CONCLUSIONS: The 16p11.2 deletion impacts in a quantitative and independent manner FSIQ, behaviour and body mass index, possibly through direct influences on neural circuitry. Although non-specific, these features are clinically significant and reproducible. Lastly, this study demonstrates the necessity of studying large patient cohorts ascertained through multiple methods to characterise the clinical consequences of rare variants involved in common diseases.

Association of adiposity genetic variants with menarche timing in 92,105 women of European descent.

Obesity is of global health concern. There are well-described inverse relationships between female pubertal timing and obesity. Recent genome-wide association studies of age at menarche identified several obesity-related variants. Using data from the ReproGen Consortium, we employed meta-analytical techniques to estimate the associations of 95 a priori and recently identified obesity-related (body mass index (weight (kg)/height (m)(2)), waist circumference, and waist:hip ratio) single-nucleotide polymorphisms (SNPs) with age at menarche in 92,116 women of European descent from 38 studies (1970-2010), in order to estimate associations between genetic variants associated with central or overall adiposity and pubertal timing in girls. Investigators in each study performed a separate analysis of associations between the selected SNPs and age at menarche (ages 9-17 years) using linear regression models and adjusting for birth year, site (as appropriate), and population stratification. Heterogeneity of effect-measure estimates was investigated using meta-regression. Six novel associations of body mass index loci with age at menarche were identified, and 11 adiposity loci previously reported to be associated with age at menarche were confirmed, but none of the central adiposity variants individually showed significant associations. These findings suggest complex genetic relationships between menarche and overall obesity, and to a lesser extent central obesity, in normal processes of growth and development.

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.

Variants in MTNR1B influence fasting glucose levels.

To identify previously unknown genetic loci associated with fasting glucose concentrations, we examined the leading association signals in ten genome-wide association scans involving a total of 36,610 individuals of European descent. Variants in the gene encoding melatonin receptor 1B (MTNR1B) were consistently associated with fasting glucose across all ten studies. The strongest signal was observed at rs10830963, where each G allele (frequency 0.30 in HapMap CEU) was associated with an increase of 0.07 (95% CI = 0.06-0.08) mmol/l in fasting glucose levels (P = 3.2 x 10(-50)) and reduced beta-cell function as measured by homeostasis model assessment (HOMA-B, P = 1.1 x 10(-15)). The same allele was associated with an increased risk of type 2 diabetes (odds ratio = 1.09 (1.05-1.12), per G allele P = 3.3 x 10(-7)) in a meta-analysis of 13 case-control studies totaling 18,236 cases and 64,453 controls. Our analyses also confirm previous associations of fasting glucose with variants at the G6PC2 (rs560887, P = 1.1 x 10(-57)) and GCK (rs4607517, P = 1.0 x 10(-25)) loci.

Discovery of amino acid variants in the human glucose-dependent insulinotropic polypeptide (GIP) receptor: the impact on the pancreatic beta cell responses and functional expression studies in Chinese hamster fibroblast cells.

The two incretins, glucose-dependent insulinotropic polypeptide (GIP) and glucagon-like peptide-1 (GLP-1), are insulinotropic factors released from the small intestine to the blood stream in response to oral glucose ingestion. The insulinotropic effect of GLP-1 is maintained in patients with Type II (non-insulin-dependent) diabetes mellitus, whereas, for unknown reasons, the effect of GIP is diminished or lacking. We defined the exon-intron boundaries of the human GIP receptor, made a mutational analysis of the gene and identified two amino acid substitutions, A207 V and E354Q. In an association study of 227 Caucasian Type II diabetic patients and 224 matched glucose tolerant control subjects, the allelic frequency of the A207 V polymorphism was 1.1% in Type II diabetic patients and 0.7% in control subjects (p = 0.48), whereas the allelic frequency of the codon 354 polymorphism was 24.9% in Type II diabetic patients versus 23.2% in control subjects. Interestingly, the glucose tolerant subjects (6% of the population) who were homozygous for the codon 354 variant had on average a 14% decrease in fasting serum C-peptide concentration (p = 0.01) and an 11% decrease in the same variable 30 min after an oral glucose load (p = 0.03) compared with subjects with the wild-type receptor. Investigation of the function of the two GIP receptor variants in Chinese hamster fibroblasts showed, however, that the GIP-induced cAMP formation and the binding of GIP to cells expressing the variant receptors were not different from the findings in cells expressing the wildtype GIP receptor. In conclusion, amino acid variants in the GIP receptor are not associated with random Type II diabetes in patients of Danish Caucasian origin or with altered GIP binding and GIP-induced cAMP production when stably transfected in Chinese hamster fibroblasts. The finding of an association between homozygosity for the codon 354 variant and reduced fasting and post oral glucose tolerance test (OGTT) serum C-peptide concentrations, however, calls for further investigations and could suggest that GIP even in the fasting state regulates the beta-cell secretory response.