Gene Set Analysis for improving genetic association studies

Introduction. Genetic epidemiology is focused on the study of the genetic causes that determine health and diseases in populations. To achieve this goal a common strategy is to explore differences in genetic variability between diseased and nondiseased individuals. Usual markers of genetic variability are single nucleotide polymorphisms (SNPs) which are changes in just one base in the genome. The usual statistical approach in genetic epidemiology study is a marginal analysis, where each SNP is analyzed separately for association with the phenotype. Motivation. It has been observed, that for common diseases the single-SNP analysis is not very powerful for detecting genetic causing variants. In this work, we consider Gene Set Analysis (GSA) as an alternative to standard marginal association approaches. GSA aims to assess the overall association of a set of genetic variants with a phenotype and has the potential to detect subtle effects of variants in a gene or a pathway that might be missed when assessed individually. Objective. We present a new optimized implementation of a pair of gene set analysis methodologies for analyze the individual evidence of SNPs in biological pathways. We perform a simulation study for exploring the power of the proposed methodologies in a set of scenarios with different number of causal SNPs under different effect sizes. In addition, we compare the results with the usual single-SNP analysis method. Moreover, we show the advantage of using the proposed gene set approaches in the context of an Alzheimer disease case-control study where we explore the Reelin signal pathway.

Short stature and primary ovarian insufficiency possibly due to chromosomal position effect in a balanced X;1 translocation.

BACKGROUND: Primary ovarian insufficiency (POI) is defined as a primary ovarian defect characterized by absent menarche (primary amenorrhea), a decrease in the initial primordial follicle number, high follicle-stimulating hormone (FSH) levels and hypoestrogenism. Although the etiology of a majority of POI cases is not yet identified, several data suggest that POI has a strong genetic component. Conventional cytogenetic and molecular analyses have identified regions of the X chromosome that are associated with ovarian function, as well as POI candidate genes, such as FMR1 and DIAPH2. Here we describe a 10.5-year-old girl presenting with high FSH and luteinizing hormone (LH) levels, pathologic GH stimulation arginine and clonidine tests, short stature, pterygium, ovarian dysgenesis, hirsutism and POI. RESULTS: Cytogenetic analysis demonstrated a balanced reciprocal translocation between the q arms of chromosomes X and 1, with breakpoints falling in Xq21 and 1q41 bands. Molecular studies did not unravel any chromosome microdeletion/microduplication, and no XIST-mediated inactivation was found on the derivative chromosome 1. Interestingly, through immunofluorescence assays, we found that part of the Xq21q22 trait, translocated to chromosome 1q41, was late replicating and therefore possibly inactivated in 30 % metaphases both in lymphocytes and skin fibroblasts, in addition to a skewed 100 % inactivation of the normal X chromosome. These findings suggest that a dysregulation of gene expression might occur in this region. Two genes mapping to the Xq translocated region, namely DIAPH2 and FMR1, were found overexpressed if compared with controls. CONCLUSIONS: We report a case in which gonadal dysgenesis and POI are associated with over-expression of DIAPH2 gene and of FMR1 gene in wild type form. We hypothesize that this over-expression is possibly due to a phenomenon known as "chromosomal position effect", which accounts for gene expression variations depending on their localization within the nucleus. For the same effect a double mosaic inactivation of genes mapping to the Xq21-q22 region, demonstrated by immunofluorescence assays, may be the cause of a functional Xq partial monosomy leading to most Turner traits of the proband's phenotype.

Buried in the Middle but Guilty: Intronic Mutations in the TCIRG1 Gene Cause Human Autosomal Recessive Osteopetrosis.

Autosomal recessive osteopetrosis (ARO) is a rare genetic bone disease with genotypic and phenotypic heterogeneity, sometimes translating into delayed diagnosis and treatment. In particular, cases of intermediate severity often constitute a diagnostic challenge and represent good candidates for exome sequencing. Here, we describe the tortuous path to identification of the molecular defect in two siblings, in which osteopetrosis diagnosed in early childhood followed a milder course, allowing them to reach the adult age in relatively good conditions with no specific therapy. No clearly pathogenic mutation was identified either with standard amplification and resequencing protocols or with exome sequencing analysis. While evaluating the possible impact of a 3'UTR variant on the TCIRG1 expression, we found a novel single nucleotide change buried in the middle of intron 15 of the TCIRG1 gene, about 150 nucleotides away from the closest canonical splice site. By sequencing a number of independent cDNA clones covering exons 14 to 17, we demonstrated that this mutation reduced splicing efficiency but did not completely abrogate the production of the normal transcript. Prompted by this finding, we sequenced the same genomic region in 33 patients from our unresolved ARO cohort and found three additional novel single nucleotide changes in a similar location and with a predicted disruptive effect on splicing, further confirmed in one of them at the transcript level. Overall, we identified an intronic region in TCIRG1 that seems to be particularly prone to splicing mutations, allowing the production of a small amount of protein sufficient to reduce the severity of the phenotype usually associated with TCIRG1 defects. On this basis, we would recommend including TCIRG1 not only in the molecular work-up of severe infantile osteopetrosis but also in intermediate cases and carefully evaluating the possible effects of intronic changes. © 2015 American Society for Bone and Mineral Research.

Migration and horizontal gene transfer divide microbial genomes into multiple niches.

Horizontal gene transfer is central to microbial evolution, because it enables genetic regions to spread horizontally through diverse communities. However, how gene transfer exerts such a strong effect is not understood. Here we develop an eco-evolutionary model and show how genetic transfer, even when rare, can transform the evolution and ecology of microbes. We recapitulate existing models, which suggest that asexual reproduction will overpower horizontal transfer and greatly limit its effects. We then show that allowing immigration completely changes these predictions. With migration, the rates and impacts of horizontal transfer are greatly increased, and transfer is most frequent for loci under positive natural selection. Our analysis explains how ecologically important loci can sweep through competing strains and species. In this way, microbial genomes can evolve to become ecologically diverse where different genomic regions encode for partially overlapping, but distinct, ecologies. Under these conditions ecological species do not exist, because genes, not species, inhabit niches.

Functional studies of a germ-line polymorphism at codon 47 within the p53 gene.

A rare germ-line polymorphism in codon 47 of the p53 gene replaces the wild-type proline (CCG) with a serine (TCG). Restriction analysis of 101 human samples revealed the frequency of the rare allele to be 0% (n = 69) in Caucasians and 4.7% (3/64, n = 32) among African-Americans. To investigate the consequence of this amino acid substitution, a cDNA construct (p53 mut47ser) containing the mutation was introduced into a lung adenocarcinoma cell line (Calu-6) that does not express p53. A growth suppression similar to that obtained after introduction of a wild-type p53 cDNA construct was observed, in contrast to the result obtained by introduction of p53 mut143ala. Furthermore, expression of neither p53 mut47ser nor wild-type p53 was tolerated by growing cells. In transient expression assays, both mut47ser and wild-type p53 activated the expression of a reporter gene linked to a p53 binding sequence (PG13-CAT) and inhibited the expression of the luciferase gene under the control of the Rous sarcoma virus promoter (RSVluc). In the same assay, mut143ala did not activate the expression of PG13-CAT and produced only a slight inhibitory effect on RSVluc. These findings indicate that the p53 variant with a serine at codon 47 should be considered as a rare germ-line polymorphism that does not alter the growth-suppression activity of p53.

Transition mutation in codon 248 of the p53 tumor suppressor gene induced by reactive oxygen species and a nitric oxide-releasing compound.

Exposing the human bronchial epithelial cell line BEAS-2B to the nitric oxide (NO) donor sodium 1-(N,N-diethylamino)diazen-1-ium-1, 2-diolate (DEA/NO) at an initial concentration of 0.6 mM while generating superoxide ion at the rate of 1 microM/min with the hypoxanthine/xanthine oxidase (HX/XO) system induced C:G-->T:A transition mutations in codon 248 of the p53 gene. This pattern of mutagenicity was not seen by 'fish-restriction fragment length polymorphism/polymerase chain reaction' (fish-RFLP/PCR) on exposure to DEA/NO alone, however, exposure to HX/XO led to various mutations, suggesting that co-generation of NO and superoxide was responsible for inducing the observed point mutation. DEA/NO potentiated the ability of HX/XO to induce lipid peroxidation as well as DNA single- and double-strand breaks under these conditions, while 0.6 mM DEA/NO in the absence of HX/XO had no significant effect on these parameters. The results show that a point mutation seen at high frequency in certain common human tumors can be induced by simultaneous exposure to reactive oxygen species and a NO source.

The effect of learning on the evolution of new courtship behavior: a simulation model

The fact that individuals learn can change the relationship between genotype and phenotype in the population, and thus affect the evolutionary response to selection. Here we ask how male ability to learn from female response affects the evolution of a novel male behavioral courtship trait under pre-existing female preference (sensory drive). We assume a courtship trait which has both a genetic and a learned component, and a two-level female response to males. With individual-based simulations we show that, under this scenario, learning generally increases the strength of selection on the genetic component of the courtship trait, at least when the population genetic mean is still low. As a consequence, learning not only accelerates the evolution of the courtship trait, but also enables it when the trait is costly, which in the absence of learning results in an adaptive valley. Furthermore, learning can enable the evolution of the novel trait in the face of gene flow mediated by immigration of males that show superior attractiveness to females based on another, non-heritable trait. However, rather than increasing monotonically with the speed of learning, the effect of learning on evolution is maximized at intermediate learning rates. This model shows that, at least under some scenarios, the ability to learn can drive the evolution of mating behaviors through a process equivalent to Waddington's genetic assimilation.

The Armc10/SVH gene: Genome context, regulation of mitochondrial dynamics and protection against Aβ-induced mitochondrial fragmentation

Mitochondrial function and dynamics are essential for neurotransmission, neural function and neuronal viability. Recently, we showed that the eutherian-specific Armcx gene cluster (Armcx1-6 genes), located in the X chromosome, encodes for a new family of proteins that localise to mitochondria, regulating mitochondrial trafficking. The Armcx gene cluster evolved by retrotransposition of the Armc10 gene mRNA, which is present in all vertebrates and is considered to be the ancestor gene. Here we investigate the genomic organisation, mitochondrial functions and putative neuroprotective role of the Armc10 ancestor gene. The genomic context of the Armc10 locus shows considerable syntenic conservation among vertebrates, and sequence comparisons and CHIP-data suggest the presence of at least three conserved enhancers. We also show that the Armc10 protein localises to mitochondria and that it is highly expressed in the brain. Furthermore, we show that Armc10 levels regulate mitochondrial trafficking in neurons, but not mitochondrial aggregation, by controlling the number of moving mitochondria. We further demonstrate that the Armc10 protein interacts with the KIF5/Miro1-2/Trak2 trafficking complex. Finally, we show that overexpression of Armc10 in neurons prevents A beta-induced mitochondrial fission and neuronal death. Our data suggest both conserved and differential roles of the Armc10/Armcx gene family in regulating mitochondrial dynamics in neurons, and underscore a protective effect of the Armc10 gene against A beta-induced toxicity. Overall, our findings support a further degree of regulation of mitochondrial dynamics in the brain of more evolved mammals.

The Armc10/SVH gene: Genome context, regulation of mitochondrial dynamics and protection against Aβ-induced mitochondrial fragmentation

Mitochondrial function and dynamics are essential for neurotransmission, neural function and neuronal viability. Recently, we showed that the eutherian-specific Armcx gene cluster (Armcx1-6 genes), located in the X chromosome, encodes for a new family of proteins that localise to mitochondria, regulating mitochondrial trafficking. The Armcx gene cluster evolved by retrotransposition of the Armc10 gene mRNA, which is present in all vertebrates and is considered to be the ancestor gene. Here we investigate the genomic organisation, mitochondrial functions and putative neuroprotective role of the Armc10 ancestor gene. The genomic context of the Armc10 locus shows considerable syntenic conservation among vertebrates, and sequence comparisons and CHIP-data suggest the presence of at least three conserved enhancers. We also show that the Armc10 protein localises to mitochondria and that it is highly expressed in the brain. Furthermore, we show that Armc10 levels regulate mitochondrial trafficking in neurons, but not mitochondrial aggregation, by controlling the number of moving mitochondria. We further demonstrate that the Armc10 protein interacts with the KIF5/Miro1-2/Trak2 trafficking complex. Finally, we show that overexpression of Armc10 in neurons prevents A beta-induced mitochondrial fission and neuronal death. Our data suggest both conserved and differential roles of the Armc10/Armcx gene family in regulating mitochondrial dynamics in neurons, and underscore a protective effect of the Armc10 gene against A beta-induced toxicity. Overall, our findings support a further degree of regulation of mitochondrial dynamics in the brain of more evolved mammals.

Divergence with gene flow and fine-scale phylogeographical structure in the wedge-billed woodcreeper, Glyphorynchus spirurus, a Neotropical rainforest bird.

Determining the relative roles of vicariance and selection in restricting gene flow between populations is of central importance to the evolutionary process of population divergence and speciation. Here we use molecular and morphological data to contrast the effect of isolation (by mountains and geographical distance) with that of ecological factors (altitudinal gradients) in promoting differentiation in the wedge-billed woodcreeper, Glyphorynchus spirurus, a tropical forest bird, in Ecuador. Tarsus length and beak size increased relative to body size with altitude on both sides of the Andes, and were correlated with the amount of moss on tree trunks, suggesting the role of selection in driving adaptive divergence. In contrast, molecular data revealed a considerable degree of admixture along these altitudinal gradients, suggesting that adaptive divergence in morphological traits has occurred in the presence of gene flow. As suggested by mitochondrial DNA sequence data, the Andes act as a barrier to gene flow between ancient subspecific lineages. Genome-wide amplified fragment length polymorphism markers reflected more recent patterns of gene flow and revealed fine-scale patterns of population differentiation that were not detectable with mitochondrial DNA, including the differentiation of isolated coastal populations west of the Andes. Our results support the predominant role of geographical isolation in driving genetic differentiation in G. spirurus, yet suggest the role of selection in driving parallel morphological divergence along ecological gradients.

Biological Characterization of Gene Response to Insulin-Induced Hypoglycemia in Mouse Retina.

Glucose is the most important metabolic substrate of the retina and maintenance of normoglycemia is an essential challenge for diabetic patients. Chronic, exaggerated, glycemic excursions could lead to cardiovascular diseases, nephropathy, neuropathy and retinopathy. We recently showed that hypoglycemia induced retinal cell death in mouse via caspase 3 activation and glutathione (GSH) decrease. Ex vivo experiments in 661W photoreceptor cells confirmed the low-glucose induction of death via superoxide production and activation of caspase 3, which was concomitant with a decrease of GSH content. We evaluate herein retinal gene expression 4 h and 48 h after insulin-induced hypoglycemia. Microarray analysis demonstrated clusters of genes whose expression was modified by hypoglycemia and we discuss the potential implication of those genes in retinal cell death. In addition, we identify by gene set enrichment analysis, three important pathways, including lysosomal function, GSH metabolism and apoptotic pathways. Then we tested the effect of recurrent hypoglycemia (three successive 4h periods of hypoglycemia spaced by 48 h recovery) on retinal cell death. Interestingly, exposure to multiple hypoglycemic events prevented GSH decrease and retinal cell death, or adapted the retina to external stress by restoring GSH level comparable to control situation. We hypothesize that scavenger GSH is a key compound in this apoptotic process, and maintaining "normal" GSH level, as well as a strict glycemic control, represents a therapeutic challenge in order to avoid side effects of diabetes, especially diabetic retinopathy.

DNA Topoisomerase I Gene Copy Number and mRNA Expression Assessed as Predictive Biomarkers for Adjuvant Irinotecan in Stage II/III Colon Cancer.

PURPOSE: Prospective-retrospective assessment of theTOP1gene copy number andTOP1mRNA expression as predictive biomarkers for adjuvant irinotecan in stage II/III colon cancer. EXPERIMENTAL DESIGN: Formalin-fixed, paraffin-embedded tissue microarrays were obtained from an adjuvant colon cancer trial (PETACC3) where patients were randomized to 5-fluorouracil/folinic acid with or without additional irinotecan.TOP1copy number status was analyzed by fluorescencein situhybridization (FISH) using aTOP1/CEN20 dual-probe combination.TOP1mRNA data were available from previous analyses. RESULTS: TOP1FISH and follow-up data were obtained from 534 patients.TOP1gain was identified in 27% using a single-probe enumeration strategy (≥4TOP1signals per cell) and in 31% when defined by aTOP1/CEN20 ratio ≥ 1.5. The effect of additional irinotecan was not dependent onTOP1FISH status.TOP1mRNA data were available from 580 patients with stage III disease. Benefit of irinotecan was restricted to patients characterized byTOP1mRNA expression ≥ third quartile (RFS: HRadjusted, 0.59;P= 0.09; OS: HRadjusted, 0.44;P= 0.03). The treatment byTOP1mRNA interaction was not statistically significant, but in exploratory multivariable fractional polynomial interaction analysis, increasingTOP1mRNA values appeared to be associated with increasing benefit of irinotecan. CONCLUSIONS: In contrast to theTOP1copy number, a trend was demonstrated for a predictive property ofTOP1mRNA expression. On the basis ofTOP1mRNA, it might be possible to identify a subgroup of patients where an irinotecan doublet is a clinically relevant option in the adjuvant setting of colon cancer.Clin Cancer Res; 22(7); 1621-31. ©2015 AACR.

Analysis of a Plant Complex Resistance Gene Locus Underlying Immune-Related Hybrid Incompatibility and Its Occurrence in Nature

Mechanisms underlying speciation in plants include detrimental (incompatible) genetic interactions between parental alleles that incur a fitness cost in hybrids. We reported on recessive hybrid incompatibility between an Arabidopsis thaliana strain from Poland, Landsberg erecta (Ler), and many Central Asian A. thaliana strains. The incompatible interaction is determined by a polymorphic cluster of Toll/interleukin-1 receptor-nucleotide binding-leucine rich repeat (TNL) RPP1 (Recognition of Peronospora parasitica1)-like genes in Ler and alleles of the receptor-like kinase Strubbelig Receptor Family 3 (SRF3) in Central Asian strains Kas-2 or Kond, causing temperature-dependent autoimmunity and loss of growth and reproductive fitness. Here, we genetically dissected the RPP1-like Ler locus to determine contributions of individual RPP1-like Ler (R1R8) genes to the incompatibility. In a neutral background, expression of most RPP1-like Ler genes, except R3, has no effect on growth or pathogen resistance. Incompatibility involves increased R3 expression and engineered R3 overexpression in a neutral background induces dwarfism and sterility. However, no individual RPP1-like Ler gene is sufficient for incompatibility between Ler and Kas-2 or Kond, suggesting that co-action of at least two RPP1-like members underlies this epistatic interaction. We find that the RPP1-like Ler haplotype is frequent and occurs with other Ler RPP1-like alleles in a local population in Gorzów Wielkopolski (Poland). Only Gorzów individuals carrying the RPP1-like Ler haplotype are incompatible with Kas-2 and Kond, whereas other RPP1-like alleles in the population are compatible. Therefore, the RPP1-like Ler haplotype has been maintained in genetically different individuals at a single site, allowing exploration of forces shaping the evolution of RPP1-like genes at local and regional population scales.

Early adversity and 5-HTT-BDNF genes: New evidences of Gene-Environment interactions on depressive symptoms in a general population

Background Adverse childhood experiences have been described as one of the major environmental risk factors for depressive disorder. Similarly, the deleterious impact of early traumatic experiences on depression seems to be moderated by individual genetic variability. Serotonin transporter (5-HTT) and brain-derived neurotrophic factor (BDNF) modulate the effect of childhood adversity on adult depression, although inconsistencies across studies have been found. Moreover, the gene×environment (G×E) interaction concerning the different types of childhood adversity remains poorly understood. The aim of this study was to analyse the putative interaction between the 5-HTT gene (5-HTTLPR polymorphism), the BDNF gene (Val66Met polymorphism) and childhood adversity in accounting for adult depressive symptoms. Method A sample of 534 healthy individuals filled in self-report questionnaires of depressive symptomatology [the Symptom Check List 90 Revised (SCL-90-R)] and different types of childhood adversities [the Childhood Trauma Questionnaire (CTQ)]. The 5-HTTLPR polymorphism (5-HTT gene) and the Val66Met polymorphism (BDNF gene) were genotyped in the whole sample. Results Total childhood adversity (β=0.27, p<0.001), childhood sexual abuse (CSA; β=0.17, p<0.001), childhood emotional abuse (β=0.27, p<0.001) and childhood emotional neglect (β=0.22, p<0.001) had an impact on adult depressive symptoms. CSA had a greater impact on depressive symptoms in Met allele carriers of the BDNF gene than in the Val/Val group (F=5.87, p<0.0001), and in S carriers of the 5-HTTLPR polymorphism (5-HTT gene) (F=5.80, p<0.0001). Conclusions Childhood adversity per se predicted higher levels of adult depressive symptoms. In addition, BDNF Val66Met and 5-HTTLPR polymorphisms seemed to moderate the effect of CSA on adult depressive symptoms.