Gene Expression Profiling and Association Studies Implicate the Neuregulin Signaling Pathway in Behçet's Disease Susceptibility

Behçet's disease (BD) is a complex disease with genetic and environmental risk factors implicated in its etiology; however, its pathophysiology is poorly understood. To decipher BD's genetic underpinnings, we combined gene expression profiling with pathway analysis and association studies. We compared the gene expression profiles in peripheral blood mononuclear cells (PBMCs) of 15 patients and 14 matched controls using Affymetrix microarrays and found that the neuregulin signaling pathway was over-represented among the differentially expressed genes. The Epiregulin (EREG), Amphiregulin (AREG), and Neuregulin-1 (NRG1) genes of this pathway stand out as they are also among the top differentially expressed genes. Twelve haplotype tagging SNPs at the EREG-AREG locus and 15 SNPs in NRG1 found associated in at least one published BD genome-wide association study were tested for association with BD in a dataset of 976 Iranian patients and 839 controls. We found a novel association with BD for the rs6845297 SNP located downstream of EREG, and replicated three associations at NRG1 (rs4489285, rs383632, and rs1462891). Multifactor dimensionality reduction analysis indicated the existence of epistatic interactions between EREG and NRG1 variants. EREG-AREG and NRG1, which are members of the epidermal growth factor (EGF) family, seem to modulate BD susceptibility through main effects and gene–gene interactions. These association findings support a role for the EGF/ErbB signaling pathway inBD pathogenesis that warrants further investigation and highlight the importance of combining genetic and genomic approaches to dissect the genetic architecture of complex diseases.

Functional genomic analysis of heat stress in Vitis vinifera

Grapevine (Vitis vinifera) is one of most agro-economically important fruit crops worldwide, with a special relevance in Portugal where over 300 varieties are used for wine production. Due to global warming, temperature stress is currently a serious issue affecting crop production especially in temperate climates. Mobile genetic elements such as retrotransposons have been shown to be involved in environmental stress induced genetic and epigenetic modifications. In this study, sequences related to Grapevine Retrotransposon 1 (Gret1) were utilized to determine heat induced genomic and transcriptomic modifications in Touriga Nacional, a traditional Portuguese grapevine variety. For this purpose, growing canes were treated to 42 oC for four hours and leaf genomic DNA and RNA was utilized for various techniques to observe possible genomic alterations and variation in transcription levels of coding and non-coding sequences between non-treated plants and treated plants immediately after heat stress (HS-0 h) or after a 24 hour recovery period (HS-24 h). Heat stress was found to induce a significant decrease in Gret1 related sequences in HS-24 h leaves, indicating an effect of heat stress on genomic structure. In order to identify putative heat induced DNA modifications, genome wide approaches such as Amplified Fragment Length Polymorphism were utilized. This resulted in the identification of a polymorphic DNA fragment in HS-0 h and HS-24 h leaves whose sequence mapped to a genomic region flanking a house keeping gene (NADH) that is represented in multiple copies in the Vitis vinifera genome. Heat stress was also found to affect the transcript levels of various non-coding and gene coding sequences. Accordingly, quantitative real time PCR results established that Gret1 related sequences are up regulated immediately after heat stress whereas the level of transcript of genes involved in identification and repair of double strand breaks are significantly down regulated in HS-0 h plants. Taken together, the results of this work demonstrated heat stress affects both genomic integrity and transcription levels.

Urothelial bladder cancer progression: lessons learned from the bench

Urothelial bladder carcinoma (UBC) is an intricate malignancy with a variable natural history and clinical behavior. Despite developments in diagnosis/prognosis refinement and treatment modalities, the recurrence rate is high, and progression from non-muscle to muscle invasive UBC commonly leads to metastasis. Moreover, patients with muscle-invasive or extra-vesical disease often fail the standard chemotherapy treatment, and overall survival rates are poor. Thus, UBC remains a challenge in the oncology field, representing an ideal candidate for research on biomarkers that could identify patients at increased risk of recurrence, progression, and chemo-refractoriness. However, progress toward personalized medicine has been hampered by the unique genetic complexity of UBC. Recent genome-wide expression and sequencing studies have brought new insights into its molecular features, pathogenesis and clinical diversity, revealing a landscape where classical pathology is intersected by the novel and heterogeneous molecular groups. Hence, it seems plausible to postulate that only an integrated signature of prognostic/predictive biomarkers inherent in different cancer hallmarks will reach clinical validation. In this review, we have summarized ours and others' research into novel putative biomarkers of progression and chemoresistance that encompass several hallmarks of cancer: tumor neovascularization, invasion and metastasis, and energy metabolism reprogramming of the tumor microenvironment.

Extensive admixture and selective pressure across the Sahel Belt

Genome-wide studies of African populations have the potential to reveal powerful insights into the evolution of our species, as these diverse populations have been exposed to intense selective pressures imposed by infectious diseases, diet, and environmental factors. Within Africa, the Sahel Belt extensively overlaps the geographical center of several endemic infections such as malaria, trypanosomiasis, meningitis, and hemorrhagic fevers. We screened 2.5 million single nucleotide polymorphisms in 161 individuals from 13 Sahelian populations, which together with published data cover Western, Central, and Eastern Sahel, and include both nomadic and sedentary groups. We confirmed the role of this Belt as a main corridor for human migrations across the continent. Strong admixture was observed in both Central and Eastern Sahelian populations, with North Africans and Near Eastern/Arabians, respectively, but it was inexistent in Western Sahelian populations. Genome-wide local ancestry inference in admixed Sahelian populations revealed several candidate regions that were significantly enriched for non-autochthonous haplotypes, and many showed to be under positive selection. The DARC gene region in Arabs and Nubians was enriched for African ancestry, whereas the RAB3GAP1/LCT/MCM6 region in Oromo, the TAS2R gene family in Fulani, and the ALMS1/NAT8 in Turkana and Samburu were enriched for non-African ancestry. Signals of positive selection varied in terms of geographic amplitude. Some genomic regions were selected across the Belt, the most striking example being the malaria-related DARC gene. Others were Western-specific (oxytocin, calcium, and heart pathways), Eastern-specific (lipid pathways), or even population-restricted (TAS2R genes in Fulani, which may reflect sexual selection).

Type 2 diabetes-related variants influence the risk of developing multiple myeloma: results from the IMMEnSE consortium

Type 2 diabetes (T2D) has been suggested to be a risk factor for multiple myeloma (MM), but the relationship between the two traits is still not well understood. The aims of this study were to evaluate whether 58 genome-wide-association-studies (GWAS)-identified common variants for T2D influence the risk of developing MM and to determine whether predictive models built with these variants might help to predict the disease risk. We conducted a case–control study including 1420 MM patients and 1858 controls ascertained through the International Multiple Myeloma (IMMEnSE) consortium. Subjects carrying the KCNQ1rs2237892T allele or the CDKN2A-2Brs2383208G/G, IGF1rs35767T/T and MADDrs7944584T/T genotypes had a significantly increased risk of MM (odds ratio (OR)=1.32–2.13) whereas those carrying the KCNJ11rs5215C, KCNJ11rs5219T and THADArs7578597C alleles or the FTOrs8050136A/A and LTArs1041981C/C genotypes showed a significantly decreased risk of developing the disease (OR=0.76–0.85). Interestingly, a prediction model including those T2D-related variants associated with the risk of MM showed a significantly improved discriminatory ability to predict the disease when compared to a model without genetic information (area under the curve (AUC)=0.645 vs AUC=0.629; P=4.05×10-06). A gender-stratified analysis also revealed a significant gender effect modification for ADAM30rs2641348 and NOTCH2rs10923931 variants (Pinteraction=0.001 and 0.0004, respectively). Men carrying the ADAM30rs2641348C and NOTCH2rs10923931T alleles had a significantly decreased risk of MM whereas an opposite but not significant effect was observed in women (ORM=0.71 and ORM=0.66 vs ORW=1.22 and ORW=1.15, respectively). These results suggest that TD2-related variants may influence the risk of developing MM and their genotyping might help to improve MM risk prediction models.

Resolving the ancestry of Austronesian‑speaking populations

There are two very different interpretations of the prehistory of Island Southeast Asia (ISEA), with genetic evidence invoked in support of both. The "out-of-Taiwan" model proposes a major Late Holocene expansion of Neolithic Austronesian speakers from Taiwan. An alternative, proposing that Late Glacial/postglacial sea-level rises triggered largely autochthonous dispersals, accounts for some otherwise enigmatic genetic patterns, but fails to explain the Austronesian language dispersal. Combining mitochondrial DNA (mtDNA), Y-chromosome and genome-wide data, we performed the most comprehensive analysis of the region to date, obtaining highly consistent results across all three systems and allowing us to reconcile the models. We infer a primarily common ancestry for Taiwan/ISEA populations established before the Neolithic, but also detected clear signals of two minor Late Holocene migrations, probably representing Neolithic input from both Mainland Southeast Asia and South China, via Taiwan. This latter may therefore have mediated the Austronesian language dispersal, implying small-scale migration and language shift rather than large-scale expansion.

A genetic perspective on African prehistory

The various genetic systems (mitochondrial DNA, the Y-chromosome and the genome-wide autosomes) indicate that Africa is the most genetically diverse continent in the world and the most likely place of origin for anatomically modern humans. However, where in Africa modern humans arose and how the current genetic makeup within the continent was shaped is still open to debate. Here, we summarize the debate and focus especially on the maternally inherited mitochondrial DNA (mtDNA) and a recently revised chronology for the African mtDNA tree. We discuss the possible origin of modern humans in southern, eastern or Central Africa; the possibility of a migration from southern to eastern Africa more than 100 ka, carrying lineages within mtDNA haplogroup L0; the evidence for a climate-change-mediated population expansion in eastern Africa involving mtDNA haplogroup L3, leading to the “out-of-Africa” migration around 70–60 ka; the re-population of North Africa from the Near East around 40–30 ka suggested by mtDNA haplogroups U6 and M1; the evidence for population expansions and dispersals across the continent at the onset of the Holocene ; and the impact of the Bantu dispersals in Central, eastern and southern Africa within the last few millennia.

CLOCK-controlled polyphonic regulation of circadian rhythms through canonical and noncanonical E-boxes.

In mammalian circadian clockwork, the CLOCK-BMAL1 complex binds to DNA enhancers of target genes and drives circadian oscillation of transcription. Here we identified 7,978 CLOCK-binding sites in mouse liver by chromatin immunoprecipitation-sequencing (ChIP-Seq), and a newly developed bioinformatics method, motif centrality analysis of ChIP-Seq (MOCCS), revealed a genome-wide distribution of previously unappreciated noncanonical E-boxes targeted by CLOCK. In vitro promoter assays showed that CACGNG, CACGTT, and CATG(T/C)G are functional CLOCK-binding motifs. Furthermore, we extensively revealed rhythmically expressed genes by poly(A)-tailed RNA-Seq and identified 1,629 CLOCK target genes within 11,926 genes expressed in the liver. Our analysis also revealed rhythmically expressed genes that have no apparent CLOCK-binding site, indicating the importance of indirect transcriptional and posttranscriptional regulations. Indirect transcriptional regulation is represented by rhythmic expression of CLOCK-regulated transcription factors, such as Krüppel-like factors (KLFs). Indirect posttranscriptional regulation involves rhythmic microRNAs that were identified by small-RNA-Seq. Collectively, CLOCK-dependent direct transactivation through multiple E-boxes and indirect regulations polyphonically orchestrate dynamic circadian outputs.

IFN stimulated gene expression in the liver is a better predictor of treatment response in chronic hepatitis c than the IL28b (IFN lambda 3) genotype

Background: Therapy of chronic hepatitis C (CHC) with pegIFNa/ribavirin achieves sustained virologic response (SVR) in ~55%. Pre-activation of the endogenous interferon system in the liver is associated non-response (NR). Recently, genome-wide association studies described associations of allelic variants near the IL28B (IFNλ3) gene with treatment response and with spontaneous clearance of the virus. We investigated if the IL28B genotype determines the constitutive expression of IFN stimulated genes (ISGs) in the liver of patients with CHC. Methods: We genotyped 93 patients with CHC for 3 IL28B single nucleotide polymorphisms (SNPs, rs12979860, rs8099917, rs12980275), extracted RNA from their liver biopsies and quantified the expression of IL28B and of 8 previously identified classifier genes which discriminate between SVR and NR (IFI44L, RSAD2, ISG15, IFI22, LAMP3, OAS3, LGALS3BP and HTATIP2). Decision tree ensembles in the form of a random forest classifier were used to calculate the relative predictive power of these different variables in a multivariate analysis. Results: The minor IL28B allele (bad risk for treatment response) was significantly associated with increased expression of ISGs, and, unexpectedly, with decreased expression of IL28B. Stratification of the patients into SVR and NR revealed that ISG expression was conditionally independent from the IL28B genotype, i.e. there was an increased expression of ISGs in NR compared to SVR irrespective of the IL28B genotype. The random forest feature score (RFFS) identified IFI27 (RFFS = 2.93), RSAD2 (1.88) and HTATIP2 (1.50) expression and the HCV genotype (1.62) as the strongest predictors of treatment response. ROC curves of the IL28B SNPs showed an AUC of 0.66 with an error rate (ERR) of 0.38. A classifier with the 3 best classifying genes showed an excellent test performance with an AUC of 0.94 and ERR of 0.15. The addition of IL28B genotype information did not improve the predictive power of the 3-gene classifier. Conclusions: IL28B genotype and hepatic ISG expression are conditionally independent predictors of treatment response in CHC. There is no direct link between altered IFNλ3 expression and pre-activation of the endogenous system in the liver. Hepatic ISG expression is by far the better predictor for treatment response than IL28B genotype.

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.

SUCLG2 identified as both a determinator of CSF Aβ1-42 levels and an attenuator of cognitive decline in Alzheimer's disease.

Cerebrospinal fluid amyloid-beta 1-42 (Aβ1-42) and phosphorylated Tau at position 181 (pTau181) are biomarkers of Alzheimer's disease (AD). We performed an analysis and meta-analysis of genome-wide association study data on Aβ1-42 and pTau181 in AD dementia patients followed by independent replication. An association was found between Aβ1-42 level and a single-nucleotide polymorphism in SUCLG2 (rs62256378) (P = 2.5×10(-12)). An interaction between APOE genotype and rs62256378 was detected (P = 9.5 × 10(-5)), with the strongest effect being observed in APOE-ε4 noncarriers. Clinically, rs62256378 was associated with rate of cognitive decline in AD dementia patients (P = 3.1 × 10(-3)). Functional microglia experiments showed that SUCLG2 was involved in clearance of Aβ1-42.

Genomic, Epigenomic and Transcriptomic Molecular Characterization of Splenic Marginal Zone Lymphoma Reveals Recurrent Abnormalities in miR-182

Splenic marginal zone lymphoma (SMZL) is a low grade B-cell non-Hodgkin's lymphoma. The molecular pathology of this entity remains poorly understood. To characterise this lymphoma at the molecular level, we performed an integrated analysis of 1) genome wide genetic copy number alterations 2) gene expression profiles and 3) epigenetic DNA methylation profiles.We have previously shown that SMZL is characterised by recurrent alterations of chromosomes 7q, 6q, 3q, 9q and 18; however, gene resolution oligonucleotide array comparative genomic hybridisation did not reveal evidence of cryptic amplification or deletion in these regions. The most frequently lost 7q32 region contains a cluster of miRNAs. qRT-PCR revealed that three of these (miR-182/96/183) show underexpression in SMZL, and miR-182 is somatically mutated in >20% of cases of SMZL, as well as in >20% of cases of follicular lymphoma, and between 5-15% of cases of chronic lymphocytic leukaemia, MALT-lymphoma and hairy cell leukaemia. We conclude that miR-182 is a strong candidate novel tumour suppressor miRNA in lymphoma.The overall gene expression signature of SMZL was found to be strongly distinct fromthose of other lymphomas. Functional analysis of gene expression data revealed SMZL to be characterised by abnormalities in B-cell receptor signalling (especially through the CD19/21-PI3K/AKT pathway) and apoptotic pathways. In addition, genes involved in the response to viral infection appeared upregulated. SMZL shows a unique epigenetic profile, but analysis of differentially methylated genes showed few with methylation related transcriptional deregulation, suggesting that DNA methylation abnormalities are not a critical component of the SMZL malignant phenotype.

DQB1 locus alone explains most of the risk and protection in narcolepsy with cataplexy in Europe.

STUDY OBJECTIVE: Prior research has identified five common genetic variants associated with narcolepsy with cataplexy in Caucasian patients. To replicate and/or extend these findings, we have tested HLA-DQB1, the previously identified 5 variants, and 10 other potential variants in a large European sample of narcolepsy with cataplexy subjects. DESIGN: Retrospective case-control study. SETTING: A recent study showed that over 76% of significant genome-wide association variants lie within DNase I hypersensitive sites (DHSs). From our previous GWAS, we identified 30 single nucleotide polymorphisms (SNPs) with P < 10(-4) mapping to DHSs. Ten SNPs tagging these sites, HLADQB1, and all previously reported SNPs significantly associated with narcolepsy were tested for replication. PATIENTS AND PARTICIPANTS: For GWAS, 1,261 narcolepsy patients and 1,422 HLA-DQB1*06:02-matched controls were included. For HLA study, 1,218 patients and 3,541 controls were included. MEASUREMENTS AND RESULTS: None of the top variants within DHSs were replicated. Out of the five previously reported SNPs, only rs2858884 within the HLA region (P < 2x10(-9)) and rs1154155 within the TRA locus (P < 2x10(-8)) replicated. DQB1 typing confirmed that DQB1*06:02 confers an extraordinary risk (odds ratio 251). Four protective alleles (DQB1*06:03, odds ratio 0.17, DQB1*05:01, odds ratio 0.56, DQB1*06:09 odds ratio 0.21, DQB1*02 odds ratio 0.76) were also identified. CONCLUSION: An overwhelming portion of genetic risk for narcolepsy with cataplexy is found at DQB1 locus. Since DQB1*06:02 positive subjects are at 251-fold increase in risk for narcolepsy, and all recent cases of narcolepsy after H1N1 vaccination are positive for this allele, DQB1 genotyping may be relevant to public health policy.

Gene duplication and neofunctionalization: POLR3G and POLR3GL.

RNA polymerase III (Pol III) occurs in two versions, one containing the POLR3G subunit and the other the closely related POLR3GL subunit. It is not clear whether these two Pol III forms have the same function, in particular whether they recognize the same target genes. We show that the POLR3G and POLR3GL genes arose from a DNA-based gene duplication, probably in a common ancestor of vertebrates. POLR3G- as well as POLR3GL-containing Pol III are present in cultured cell lines and in normal mouse liver, although the relative amounts of the two forms vary, with the POLR3G-containing Pol III relatively more abundant in dividing cells. Genome-wide chromatin immunoprecipitations followed by high-throughput sequencing (ChIP-seq) reveal that both forms of Pol III occupy the same target genes, in very constant proportions within one cell line, suggesting that the two forms of Pol III have a similar function with regard to specificity for target genes. In contrast, the POLR3G promoter-not the POLR3GL promoter-binds the transcription factor MYC, as do all other promoters of genes encoding Pol III subunits. Thus, the POLR3G/POLR3GL duplication did not lead to neo-functionalization of the gene product (at least with regard to target gene specificity) but rather to neo-functionalization of the transcription units, which acquired different mechanisms of regulation, thus likely affording greater regulation potential to the cell.

Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution.

Waist-hip ratio (WHR) is a measure of body fat distribution and a predictor of metabolic consequences independent of overall adiposity. WHR is heritable, but few genetic variants influencing this trait have been identified. We conducted a meta-analysis of 32 genome-wide association studies for WHR adjusted for body mass index (comprising up to 77,167 participants), following up 16 loci in an additional 29 studies (comprising up to 113,636 subjects). We identified 13 new loci in or near RSPO3, VEGFA, TBX15-WARS2, NFE2L3, GRB14, DNM3-PIGC, ITPR2-SSPN, LY86, HOXC13, ADAMTS9, ZNRF3-KREMEN1, NISCH-STAB1 and CPEB4 (P = 1.9 × 10⁻⁹ to P = 1.8 × 10⁻⁴⁰) and the known signal at LYPLAL1. Seven of these loci exhibited marked sexual dimorphism, all with a stronger effect on WHR in women than men (P for sex difference = 1.9 × 10⁻³ to P = 1.2 × 10⁻&supl;³). These findings provide evidence for multiple loci that modulate body fat distribution independent of overall adiposity and reveal strong gene-by-sex interactions.