DUF1220-Domain Copy Number Implicated in Human Brain-Size Pathology and Evolution

DUF1220 domains show the largest human-lineage-specific increase in copy number of any protein-coding region in the human genome and map primarily to 1q21, where deletions and reciprocal duplications have been associated with microcephaly and macrocephaly, respectively. Given these findings and the high correlation between DUF1220 copy number and brain size across primate lineages (R² = 0.98; p = 1.8 × 10⁻⁶), DUF1220 sequences represent plausible candidates for underlying 1q21-associated brain-size pathologies. To investigate this possibility, we used specialized bioinformatics tools developed for scoring highly duplicated DUF1220 sequences to implement targeted 1q21 array comparative genomic hybridization on individuals (n = 42) with 1q21-associated microcephaly and macrocephaly. We show that of all the 1q21 genes examined (n = 53), DUF1220 copy number shows the strongest association with brain size among individuals with 1q21-associated microcephaly, particularly with respect to the three evolutionarily conserved DUF1220 clades CON1(p = 0.0079), CON2 (p = 0.0134), and CON3 (p = 0.0116). Interestingly, all 1q21 DUF1220-encoding genes belonging to the NBPF family show significant correlations with frontal-occipital-circumference Z scores in the deletion group. In a similar survey of a nondisease population, we show that DUF1220 copy number exhibits the strongest correlation with brain gray-matter volume (CON1, p = 0.0246; and CON2, p = 0.0334). Notably, only DUF1220 sequences are consistently significant in both disease and nondisease populations. Taken together, these data strongly implicate the loss of DUF1220 copy number in the etiology of 1q21-associated microcephaly and support the view that DUF1220 domains function as general effectors of evolutionary, pathological, and normal variation in brain size.     

Male-Biased Autosomal Effect of 16p13.11 Copy Number Variation in Neurodevelopmental Disorders

Copy number variants (CNVs) at chromosome 16p13.11 have been associated with a range of neurodevelopmental disorders including autism, ADHD, intellectual disability and schizophrenia. Significant sex differences in prevalence, course and severity have been described for a number of these conditions but the biological and environmental factors underlying such sex-specific features remain unclear. We tested the burden and the possible sex-biased effect of CNVs at 16p13.11 in a sample of 10,397 individuals with a range of neurodevelopmental conditions, clinically referred for array comparative genomic hybridisation (aCGH); cases were compared with 11,277 controls. In order to identify candidate phenotype-associated genes, we performed an interval-based analysis and investigated the presence of ohnologs at 16p13.11; finally, we searched the DECIPHER database for previously identified 16p13.11 copy number variants. In the clinical referral series, we identified 46 cases with CNVs of variable size at 16p13.11, including 28 duplications and 18 deletions. Patients were referred for various phenotypes, including developmental delay, autism, speech delay, learning difficulties, behavioural problems, epilepsy, microcephaly and physical dysmorphisms. CNVs at 16p13.11 were also present in 17 controls. Association analysis revealed an excess of CNVs in cases compared with controls (OR = 2.59; p = 0.0005), and a sex-biased effect, with a significant enrichment of CNVs only in the male subgroup of cases (OR = 5.62; p = 0.0002), but not in females (OR = 1.19, p = 0.673). The same pattern of results was also observed in the DECIPHER sample. Interval-based analysis showed a significant enrichment of case CNVs containing interval II (OR = 2.59; p = 0.0005), located in the 0.83 Mb genomic region between 15.49–16.32 Mb, and encompassing the four ohnologs NDE1, MYH11, ABCC1 and ABCC6. Our data confirm that duplications and deletions at 16p13.11 represent incompletely penetrant pathogenic mutations that predispose to a range of neurodevelopmental disorders, and suggest a sex-limited effect on the penetrance of the pathological phenotypes at the 16p13.11 locus.     

The Association between DNA Copy Number Aberrations at Chromosome 5q22 and Gastric Cancer

Background

Gastric cancer is common cancer. Discovering novel genetic biomarkers might help to identify high-risk individuals. Copy number variation (CNV) has recently been shown to influence risk for several cancers. The aim of the present study was sought to test the association between copy number at a variant region and GC.

Methods

A total of 110 gastric cancer patients and 325 healthy volunteers were enrolled in this study. We searched for a CNV and found a CNV (Variation 7468) containing part of the APC gene, the SRP19 gene and the REEP5 gene. We chose four probes targeting at APC-intron8, APC-exon9, SRP19 and REEP5 to interrogate this CNV. Specific Taqman probes labeled by different reporter fluorophores were used in a real-time PCR platform to obtain copy number. Both the original non-integer data and transformed integer data on copy number were used for analyses.

Results

Gastric caner patients had a lower non-integer copy number than controls for the APC-exon9 probe (Adjusted p = 0.026) and SRP19 probe (Adjusted p = 0.002). The analysis of integer copy number yielded a similar pattern although less significant (Adjusted p = 0.07 for APC-exon9 probe and Adjusted p = 0.02 for SRP19 probe).

Conclusions

Losses of a CNV at 5q22, especially in the DNA region surrounding APC-exon 9, may be associated with a higher risk of gastric cancer.     

Higher Fibrinogen Level is Independently Linked with the Presence and Severity of New-Onset Coronary Atherosclerosis among Han Chinese Population

Background

Fibrinogen is a coagulation/inflammatory biomarker strongly associated with atherogenesis. However, no data is currently available regarding the association of fibrinogen level with the presence and severity of new-onset coronary atherosclerosis assessed by Gensini score (GS), particularly in Han Chinese with a large sample size.

Methods and Results

We studied 2288 consecutive, new-onset subjects undergoing coronary angiography with angina-like chest pain. Clinical and laboratory data were collected. Coronary stenotic lesions were considered to be the incidence of coronary atherosclerosis. The severity of coronary stenosis was determined by the GS system. Data indicated that patients with high GS had significantly elevated fibrinogen level (p<0.001). The prevalence and severity of coronary atherosclerosis were dramatically increased according to fibrinogen tertiles. Spearman correlation analysis revealed a positive association between fibrinogen level and GS (r = 0.138, p<0.001). Multivariate logistic regression analysis demonstrated that plasma fibrinogen level was independently associated with high GS (OR = 1.275, 95% CI 1.082–1.502, p = 0.004) after adjusting for potential confounders. Moreover, fibrinogen level was also independently related to the presence of coronary atherosclerosis (fibrinogen tertile 2: OR = 1.192, 95% CI 0.889–1.598, p = 0.241; tertile 3: OR = 2.003, 95% CI 1.383–2.903, p <0.001) and high GS (fibrinogen tertile 2: OR = 1.079, 95% CI 0.833–1.397, p = 0.565; tertile 3: OR = 1.524, 95% CI 1.155–2.011, p = 0.003) in a dose-dependent manner. Receiver-operating characteristic curve analysis showed that the best fibrinogen cut-off value for predicting the severity of coronary stenosis was 3.21 g/L.

Conclusions

Higher fibrinogen level is independently linked with the presence and severity of new-onset coronary atherosclerosis in Han Chinese population.     

Impact of Gene Dosage on Gene Expression,Biological Processes and Survival in Cervical Cancer: A Genome-Wide Follow-Up Study

We investigated the role of tumor copy number (CN)–altered genome (CN-AG) in the carcinogenesis of cervical cancer (CC), especially its effect on gene expression, biological processes, and patient survival. Fifty-nine human papillomavirus 16 (HPV16)-positive CCs were investigated with microarrays–31 for mapping CN-AG and 55 for global gene expression, with 27 CCs in common. Five-year survival was investigated in 55 patients. Deletions and amplifications >2.5 Mb were defined as CN alterations. The %CN-AG varied from 0 to 32.2% (mean = 8.1±8.9). Tumors were classified as low (mean = 0.5±0.6, n = 11), medium (mean = 5.4±2.4, n = 10), or high (mean = 19.2±6.6, n = 10) CN. The highest %CN-AG was found in 3q, which contributed an average of 55% of all CN alterations. Genome-wide, only 5.3% of CN-altered genes were deregulated directly by gene dosage. In contrast, the rate in fully duplicated 3q was twice as high. Amplification of 3q explained 23.2% of deregulated genes in whole tumors (r² = 0.232, p = 0.006; analysis of variance), including genes located in 3q and other chromosomes. A total of 862 genes were deregulated exclusively in high-CN tumors, but only 22.9% were CN altered. This suggests that the remaining genes are not deregulated directly by gene dosage, but by mechanisms induced in trans by CN-altered genes. Anaphase-promoting complex/cyclosome (APC/C)-dependent proteasome proteolysis, glycolysis, and apoptosis were upregulated, whereas cell adhesion and angiogenesis were downregulated exclusively in high-CN tumors. The high %CN-AG and upregulated gene expression profile of APC/C-dependent proteasome proteolysis were associated with poor patient survival (p<0.05, log-rank test). Along with glycolysis, they were linearly associated with FIGO stage (r>0.38, p<0.01, Spearman test). Therefore, inhibition of APC/C-dependent proteasome proteolysis and glycolysis could be useful for CC treatment. However, whether they are indispensable for tumor growth remains to be demonstrated.     

Network Topologies and Convergent Aetiologies Arising from Deletions and Duplications Observed in Individuals with Autism

Autism Spectrum Disorders (ASD) are highly heritable and characterised by impairments in social interaction and communication, and restricted and repetitive behaviours. Considering four sets of de novo copy number variants (CNVs) identified in 181 individuals with autism and exploiting mouse functional genomics and known protein-protein interactions, we identified a large and significantly interconnected interaction network. This network contains 187 genes affected by CNVs drawn from 45% of the patients we considered and 22 genes previously implicated in ASD, of which 192 form a single interconnected cluster. On average, those patients with copy number changed genes from this network possess changes in 3 network genes, suggesting that epistasis mediated through the network is extensive. Correspondingly, genes that are highly connected within the network, and thus whose copy number change is predicted by the network to be more phenotypically consequential, are significantly enriched among patients that possess only a single ASD-associated network copy number changed gene (p = 0.002). Strikingly, deleted or disrupted genes from the network are significantly enriched in GO-annotated positive regulators (2.3-fold enrichment, corrected p = 2×10⁻⁵), whereas duplicated genes are significantly enriched in GO-annotated negative regulators (2.2-fold enrichment, corrected p = 0.005). The direction of copy change is highly informative in the context of the network, providing the means through which perturbations arising from distinct deletions or duplications can yield a common outcome. These findings reveal an extensive ASD-associated molecular network, whose topology indicates ASD-relevant mutational deleteriousness and that mechanistically details how convergent aetiologies can result extensively from CNVs affecting pathways causally implicated in ASD.     

Screening for EGFR Amplifications with a Novel Method and Their Significance for the Outcome of Glioblastoma Patients

Glioblastoma is a highly aggressive tumour of the central nervous system, characterised by poor prognosis irrespective of the applied treatment. The aim of our study was to analyse whether the molecular markers of glioblastoma (i.e. TP53 and IDH1 mutations, CDKN2A deletion, EGFR amplification, chromosome 7 polysomy and EGFRvIII expression) could be associated with distinct prognosis and/or response to the therapy. Moreover, we describe a method which allows for a reliable, as well as time- and cost-effective, screening for EGFR amplification and chromosome 7 polysomy with quantitative Real-Time PCR at DNA level. In the clinical data, only the patient’s age had prognostic significance (continuous: HR = 1.04; p<0.01). At the molecular level, EGFRvIII expression was associated with a better prognosis (HR = 0.37; p = 0.04). Intriguingly, EGFR amplification was associated with a worse outcome in younger patients (HR = 3.75; p<0.01) and in patients treated with radiotherapy (HR = 2.71; p = 0.03). We did not observe any difference between the patients with the amplification treated with radiotherapy and the patients without such a treatment. Next, EGFR amplification was related to a better prognosis in combination with the homozygous CDKN2A deletion (HR = 0.12; p = 0.01), but to a poorer prognosis in combination with chromosome 7 polysomy (HR = 14.88; p = 0.01). Importantly, the results emphasise the necessity to distinguish both mechanisms of the increased EGFR gene copy number (amplification and polysomy). To conclude, although the data presented here require validation in different groups of patients, they strongly advocate the consideration of the patient’s tumour molecular characteristics in the selection of the therapy.     

Variants in the 15q24/25 Locus Associate with Lung Function Decline in Active Smokers

Genetic variation in nicotinic acetylcholine receptor subunit genes (nAChRs) is associated with lung function level and chronic obstructive pulmonary disease (COPD). It is unknown whether these variants also predispose to an accelerated lung function decline. We investigated the association of nAChR susceptibility variants with lung function decline and COPD severity. The rs1051730 and rs8034191 variants were genotyped in a population-based cohort of 1,226 heavy smokers (COPACETIC) and in an independent cohort of 883 heavy smokers, of which 653 with COPD of varying severity (LEUVEN). Participants underwent pulmonary function tests at baseline. Lung function decline was assessed over a median follow-up of 3 years in COPACETIC. Current smokers homozygous for the rs1051730 A-allele or rs8034191 G-allele had significantly greater FEV₁/FVC decline than homozygous carriers of wild-type alleles (3.3% and 4.3%, p = 0.026 and p = 0.009, respectively). In the LEUVEN cohort, rs1051730 AA-carriers and rs8034191 GG-carriers had a two-fold increased risk to suffer from COPD GOLD IV (OR 2.29, 95% confidence interval [CI] = 1.11–4.75; p = 0.025 and OR = 2.42, 95% [CI] = 1.18–4.95; p = 0.016, respectively). The same risk alleles conferred, respectively, a five- and four-fold increased risk to be referred for lung transplantation because of end-stage COPD (OR = 5.0, 95% [CI] = 1.68–14.89; p = 0.004 and OR = 4.06, 95% [CI] = 1.39–11.88; p = 0.010). In Europeans, variants in nAChRs associate with an accelerated lung function decline in current smokers and with clinically relevant COPD.     

TP53 Mutational Status Is a Potential Marker for Risk Stratification in Wilms Tumour with Diffuse Anaplasia

Purpose

The presence of diffuse anaplasia in Wilms tumours (DAWT) is associated with TP53 mutations and poor outcome. As patients receive intensified treatment, we sought to identify whether TP53 mutational status confers additional prognostic information.

Patients and Methods

We studied 40 patients with DAWT with anaplasia in the tissue from which DNA was extracted and analysed for TP53 mutations and 17p loss. The majority of cases were profiled by copy number (n = 32) and gene expression (n = 36) arrays. TP53 mutational status was correlated with patient event-free and overall survival, genomic copy number instability and gene expression profiling.

Results

From the 40 cases, 22 (55%) had TP53 mutations (2 detected only after deep-sequencing), 20 of which also had 17p loss (91%); 18 (45%) cases had no detectable mutation but three had 17p loss. Tumours with TP53 mutations and/or 17p loss (n = 25) had an increased risk of recurrence as a first event (p = 0.03, hazard ratio (HR), 3.89; 95% confidence interval (CI), 1.26–16.0) and death (p = 0.04, HR, 4.95; 95% CI, 1.36–31.7) compared to tumours lacking TP53 abnormalities. DAWT carrying TP53 mutations showed increased copy number alterations compared to those with wild-type, suggesting a more unstable genome (p = 0.03). These tumours showed deregulation of genes associated with cell cycle and DNA repair biological processes.

Conclusion

This study provides evidence that TP53 mutational analysis improves risk stratification in DAWT. This requires validation in an independent cohort before clinical use as a biomarker.     

Whole-Genome Sequencing of Individuals from a Founder Population Identifies Candidate Genes for Asthma

Asthma is a complex genetic disease caused by a combination of genetic and environmental risk factors. We sought to test classes of genetic variants largely missed by genome-wide association studies (GWAS), including copy number variants (CNVs) and low-frequency variants, by performing whole-genome sequencing (WGS) on 16 individuals from asthma-enriched and asthma-depleted families. The samples were obtained from an extended 13-generation Hutterite pedigree with reduced genetic heterogeneity due to a small founding gene pool and reduced environmental heterogeneity as a result of a communal lifestyle. We sequenced each individual to an average depth of 13-fold, generated a comprehensive catalog of genetic variants, and tested the most severe mutations for association with asthma. We identified and validated 1960 CNVs, 19 nonsense or splice-site single nucleotide variants (SNVs), and 18 insertions or deletions that were out of frame. As follow-up, we performed targeted sequencing of 16 genes in 837 cases and 540 controls of Puerto Rican ancestry and found that controls carry a significantly higher burden of mutations in IL27RA (2.0% of controls; 0.23% of cases; nominal p = 0.004; Bonferroni p = 0.21). We also genotyped 593 CNVs in 1199 Hutterite individuals. We identified a nominally significant association (p = 0.03; Odds ratio (OR) = 3.13) between a 6 kbp deletion in an intron of NEDD4L and increased risk of asthma. We genotyped this deletion in an additional 4787 non-Hutterite individuals (nominal p = 0.056; OR = 1.69). NEDD4L is expressed in bronchial epithelial cells, and conditional knockout of this gene in the lung in mice leads to severe inflammation and mucus accumulation. Our study represents one of the early instances of applying WGS to complex disease with a large environmental component and demonstrates how WGS can identify risk variants, including CNVs and low-frequency variants, largely untested in GWAS.     

Dysregulation of Human β-Defensin-2 Protein in Inflammatory Bowel Disease

Background

Human β-defensin-2 (HBD2) is an antimicrobial peptide implicated in the pathogenesis of inflammatory bowel disease (IBD). Low copy number and concomitant low mRNA expression of the HBD2 gene have been implicated in susceptibility to colonic Crohn''s Disease (CD). We investigated the colonic distribution of HBD2 mRNA expression, and the contributions of genetic and environmental factors on HBD2 protein production.

Methodology/Principal Findings

We examined HBD2 mRNA expression at three colonic locations by microarray analysis of biopsies from 151 patients (53 CD, 67 ulcerative colitis [UC], 31 controls). We investigated environmental and genetic influences on HBD2 protein production using ex vivo cultured sigmoid colon biopsies from 69 patients (22 CD, 26 UC, 21 controls) stimulated with lipopolysaccharide (LPS) and/or nicotine for 24 hours. HBD2 and cytokines were measured in culture supernatants. Using DNA samples from these patients, regions in the HBD2 gene promoter were sequenced for NF-κB binding-sites and HBD2 gene copy number was determined. HBD2 mRNA expression was highest in inflamed (vs. uninflamed p = 0.0122) ascending colon in CD and in inflamed (vs. uninflamed p<0.0001) sigmoid colon in UC. HBD2 protein production was increased in inflamed UC biopsies (p = 0.0078). There was no difference in HBD2 protein production from unstimulated biopsies of CD, UC and controls. LPS-induced HBD2 production was significantly increased in CD (p = 0.0375) but not UC (p = 0.2017); this LPS-induced response was augmented by nicotine in UC (p = 0.0308) but not CD (p = 0.6872). Nicotine alone did not affect HBD2 production. HBD2 production correlated with IL8 production in UC (p<0.001) and with IL10 in CD (p<0.05). Variations in the HBD2 promoter and HBD2 gene copy number did not affect HBD2 production.

Significance/Conclusions

Colonic HBD2 was dysregulated at mRNA and protein level in IBD. Inflammatory status and stimulus but not germline variations influenced these changes.     

The Relationship between Endothelial Progenitor Cell Populations and Epicardial and Microvascular Coronary Disease—A Cellular,Angiographic and Physiologic Study

Background

Endothelial progenitor cells (EPCs) are implicated in protection against vascular disease. However, studies using angiography alone have reported conflicting results when relating EPCs to epicardial coronary artery disease (CAD) severity. Moreover, the relationship between different EPC types and the coronary microcirculation is unknown. We therefore investigated the relationship between EPC populations and coronary epicardial and microvascular disease.

Methods

Thirty-three patients with a spectrum of isolated left anterior descending artery disease were studied. The coronary epicardial and microcirculation were physiologically interrogated by measurement of fractional flow reserve (FFR), index of microvascular resistance (IMR) and coronary flow reserve (CFR). Two distinct EPC populations (early EPC and late outgrowth endothelial cells [OECs]) were isolated from these patients and studied ex vivo.

Results

There was a significant inverse relationship between circulating OEC levels and epicardial CAD severity, as assessed by FFR and angiography (r = 0.371, p = 0.04; r = -0.358, p = 0.04; respectively). More severe epicardial CAD was associated with impaired OEC migration and tubulogenesis (r = 0.59, p = 0.005; r = 0.589, p = 0.004; respectively). Patients with significant epicardial CAD (FFR<0.75) had lower OEC levels and function compared to those without hemodynamically significant stenoses (p<0.05). In contrast, no such relationship was seen for early EPC number and function, nor was there a relationship between IMR and EPCs. There was a significant relationship between CFR and OEC function.

Conclusions

EPC populations differ in regards to their associations with CAD severity. The number and function of OECs, but not early EPCs, correlated significantly with epicardial CAD severity. There was no relationship between EPCs and severity of coronary microvascular disease.     

Association of Common Genetic Variants with Lipid Traits in the Indian Population

Genome-wide association studies (GWAS) have been instrumental in identifying novel genetic variants associated with altered plasma lipid levels. However, these quantitative trait loci have not been tested in the Indian population, where there is a poorly understood and growing burden of cardiometabolic disorders. We present the association of six single nucleotide polymorphisms in 1671 sib pairs (3342 subjects) with four lipid traits: total cholesterol, triglycerides, high density lipoprotein cholesterol (HDL-C) and low density lipoprotein cholesterol (LDL-C). We also investigated the interaction effects of gender, location, fat intake and physical activity. Each copy of the risk allele of rs964184 at APOA1 was associated with 1.06 mmol/l increase in triglycerides (SE = 0.049; p = 0.006), rs3764261 at CETP with 1.02 mmol/l increase in both total cholesterol (SE = 0.042; p = 0.017) and HDL-C (SE = 0.041; p = 0.008), rs646776 at CELSR2-PSRC1-SORT1 with 0.96 mmol/l decrease in cholesterol (SE = 0.043; p = 0.0003) and 0.15 mmol/l decrease in LDL-C levels (SE = 0.043; p = 0.0003) and rs2954029 at TRIB1 with 1.02 mmol/l increase in HDL-C (SE = 0.039; p = 0.047). A combined risk score of APOA1 and CETP loci predicted an increase of 1.25 mmol/l in HDL-C level (SE = 0.312; p = 0.0007). Urban location and sex had strong interaction effects on the genetic association of most of the studied loci with lipid traits. To conclude, we validated four genetic variants (identified by GWAS in western populations) associated with lipid traits in the Indian population. The interaction effects found here may explain the sex-specific differences in lipid levels and their heritability. Urbanization appears to influence the nature of the association with GWAS lipid loci in this population. However, these findings will require replication in other Indian populations.     

CXCL10 Gene Promoter Polymorphism -1447A>G Correlates with Plasma CXCL10 Levels and is Associated with Male Susceptibility to Cerebral Malaria

The risk factors for cerebral malaria (CM) and the wide variation in clinical manifestations of malaria are poorly understood. Recent studies indicate that interferon gamma inducible chemokine, CXCL10, is a strong predictor of both human and experimental cerebral malaria. Increased plasma and cerebrospinal fluid levels of CXCL10 were tightly associated with fatal CM in Indian and Ghanaian patients. In the present study, we hypothesized that in a subset of malaria patients, CM susceptibility is associated with variation in CXCL10 expression. We determined whether polymorphisms in the CXCL10 gene promoter region played a role in the clinical status of malaria patients and addressed the genetic basis of CXCL10 expression during malaria infection. Following extensive bioinformatics analyses, two reported single nucleotide polymorphisms in the CXCL10 promoter (−135G>A [rs56061981] and −1447A>G [rs4508917]) were identified among 66 CM and 69 non-CM Indian patients using PCR-restriction fragment length polymorphism assay. Individuals with the −1447(A/G) genotype were susceptible to CM (adjusted odds ratio [AOR] = 2.60, 95% CI = 1.51–5.85, p = 0.021). In addition, individuals with the −1447(A/G) genotype had significantly higher plasma CXCL10 levels than individuals with the −1447(A/A) genotype. Stratifying patients according to gender, the observed association of CM with over expression of CXCL10 were more pronounced in males than in female patients (AOR = 5.47, 95% CI = 1.34–22.29, p = 0.018). Furthermore, −135G>A polymorphism conferred a decreased risk of CM among males (AOR = 0.19, 95% CI = 0.05–0.78, p = 0.021). Polymorphisms in the CXCL10 gene promoter sequence were associated with increased CXCL10 production, which is linked to severity of CM. These results suggest that the −1447A>G polymorphism in CXCL10 gene promoter could be partly responsible for the reported variation underlying severity of CM outcomes particularly in males.     

Age-Related Changes in the Global DNA Methylation Profile of Leukocytes Are Linked to Nutrition but Are Not Associated with the MTHFR C677T Genotype or to Functional Capacities

Global DNA methylation of peripheral blood leukocytes has been recently proposed as a potential biomarker for disease risk. However, the amplitude of the changes in DNA methylation associated with normal aging and the impacts of environmental changes on this variation are still unclear. In this context, we evaluated the association of global DNA methylation with nutritional habits, tobacco smoking, body mass index (BMI), clinical laboratory parameters, polymorphism C677T MTHFR, functional cognition and the daily practice of physical activity in a cancer-free older population. Leukocyte global DNA methylation from 126 older individuals was quantified using a high-throughput ELISA-based method. Global DNA hypomethylation was observed in older individuals when compared to a younger population (p = 0.0469), confirming changes in DNA methylation in the aging process. Furthermore, the methylation profile of elders was correlated with the daily ingestion of carbohydrates (p = 0.0494), lipids (p = 0.0494), vitamin B6 (p = 0.0421), magnesium (p = 0.0302), and also to the serum levels of total protein (p = 0.0004), alpha 2 globulin (p = 0.0013) and albumin (p = 0.0015). No statistically significant difference was observed when global DNA methylation were stratified according to C677T MTHFR genotypes (p = 0.7200), BMI (p = 0.1170), smoking habit (p = 0.4382), physical activity in daily life (p = 0.8492), scored cognitive function (p = 0.7229) or depression state (p = 0.8301). Our data indicate that age-related variations in the global DNA methylation profile of leukocytes might be modulated by the daily intake of carbohydrates, lipids, vitamin B6, and magnesium and be associated with serum protein levels, however it is independent of C677T MTHFR genotype and not correlated with BMI, smoking habit, cognitive function or the routine physical activities.     

Finished sequence and assembly of the DUF1220-rich 1q21 region using a haploid human genome

Background

Although the reference human genome sequence was declared finished in 2003, some regions of the genome remain incomplete due to their complex architecture. One such region, 1q21.1-q21.2, is of increasing interest due to its relevance to human disease and evolution. Elucidation of the exact variants behind these associations has been hampered by the repetitive nature of the region and its incomplete assembly. This region also contains 238 of the 270 human DUF1220 protein domains, which are implicated in human brain evolution and neurodevelopment. Additionally, examinations of this protein domain have been challenging due to the incomplete 1q21 build. To address these problems, a single-haplotype hydatidiform mole BAC library (CHORI-17) was used to produce the first complete sequence of the 1q21.1-q21.2 region.

Results

We found and addressed several inaccuracies in the GRCh37sequence of the 1q21 region on large and small scales, including genomic rearrangements and inversions, and incorrect gene copy number estimates and assemblies. The DUF1220-encoding NBPF genes required the most corrections, with 3 genes removed, 2 genes reassigned to the 1p11.2 region, 8 genes requiring assembly corrections for DUF1220 domains (~91 DUF1220 domains were misassigned), and multiple instances of nucleotide changes that reassigned the domain to a different DUF1220 subtype. These corrections resulted in an overall increase in DUF1220 copy number, yielding a haploid total of 289 copies. Approximately 20 of these new DUF1220 copies were the result of a segmental duplication from 1q21.2 to 1p11.2 that included two NBPF genes. Interestingly, this duplication may have been the catalyst for the evolutionarily important human lineage-specific chromosome 1 pericentric inversion.

Conclusions

Through the hydatidiform mole genome sequencing effort, the 1q21.1-q21.2 region is complete and misassemblies involving inter- and intra-region duplications have been resolved. The availability of this single haploid sequence path will aid in the investigation of many genetic diseases linked to 1q21, including several associated with DUF1220 copy number variations. Finally, the corrected sequence identified a recent segmental duplication that added 20 additional DUF1220 copies to the human genome, and may have facilitated the chromosome 1 pericentric inversion that is among the most notable human-specific genomic landmarks.     

Association of CCR2-CCR5 Haplotypes and CCL3L1 Copy Number with Kawasaki Disease,Coronary Artery Lesions,and IVIG Responses in Japanese Children

Background

The etiology of Kawasaki Disease (KD) is enigmatic, although an infectious cause is suspected. Polymorphisms in CC chemokine receptor 5 (CCR5) and/or its potent ligand CCL3L1 influence KD susceptibility in US, European and Korean populations. However, the influence of these variations on KD susceptibility, coronary artery lesions (CAL) and response to intravenous immunoglobulin (IVIG) in Japanese children, who have the highest incidence of KD, is unknown.

Methodology/Principal Findings

We used unconditional logistic regression analyses to determine the associations of the copy number of the CCL3L1 gene-containing duplication and CCR2-CCR5 haplotypes in 133 Japanese KD cases [33 with CAL and 25 with resistance to IVIG] and 312 Japanese controls without a history of KD. We observed that the deviation from the population average of four CCL3L1 copies (i.e., < or > four copies) was associated with an increased risk of KD and IVIG resistance (adjusted odds ratio (OR)  = 2.25, p = 0.004 and OR = 6.26, p = 0.089, respectively). Heterozygosity for the CCR5 HHF*2 haplotype was associated with a reduced risk of both IVIG resistance (OR = 0.21, p = 0.026) and CAL development (OR = 0.44, p = 0.071).

Conclusions/Significance

The CCL3L1-CCR5 axis may play an important role in KD pathogenesis. In addition to clinical and laboratory parameters, genetic markers may also predict risk of CAL and resistance to IVIG.     

Vitamin D Receptor Gene Polymorphisms Are Associated with Obesity and Inflammosome Activity

To explore the mechanisms underlying the suggested role of the vitamin D/vitamin D receptor (VDR) complex in the pathogenesis of obesity we performed genetic and immunologic analyses in obese and non-obese Saudi individuals without other concomitant chronic diseases. Genomic DNA was genotyped for gene single nucleotide polymorphisms (SNPs) of VDR by allelic discrimination in 402 obese (body mass index –BMI≥30 kg/m2) and 489 non-obese (BMI<30 kg/m2) Saudis. Q-PCR analyses were performed using an ABI Prism 7000 Sequence Detection System. The inflammosome pathway was analysed by PCR, cytokines and plasma lipopolysaccaride (LPS) concentrations with ELISA assays. Results showed that the VDR SNPs rs731236 (G) (TaqI) and rs1544410 (T) (Bsm-I) minor allele polymorphisms are significantly more frequent in obese individuals (p = 0.009, β = 0.086 and p = 0.028, β = 0.072, respectively). VDR haplotypes identified are positively (GTA) (p = 0.008, β = 1.560); or negatively (ACC) (p = 0.044, β = 0.766) associated with obesity and higher BMI scores. The GTA "risk" haplotype was characterized by an up-regulation of inflammosome components, a higher production of proinflammatory cytokines (p<0.05) and a lower VDR expression. Plasma LPS concentration was also increased in GTA obese individuals (p<0.05), suggesting an alteration of gut permeability leading to microbial translocation. Data herein indicate that polymorphisms affecting the vitamin D/VDR axis play a role in obesity that is associated with an ongoing degree of inflammation, possibly resulting from alterations of gut permeability and microbial translocation. These results could help the definition of VDR fingerprints that predict an increased risk of developing obesity and might contribute to the identification of novel therapeutic strategies for this metabolic condition.