Early adversity and 5-HTT-BDNF genes: new evidences of gene-environment interactions on depressive symptoms in a general population

Adverse childhood experiences have been described as one of the major environmental risk factors for depressive disorder. Likewise, the deleterious impact of early traumatic experiences on depression seems to be moderated by individual genetic variability. Serotonin transporter (5-HTT) and the Brain-Derived Neurotrophic Factor (BDNF) seem to modulate the effect of childhood adversity on adult depression, although inconsistencies across studies have been found. Moreover, the GxE interaction concerning the different types of childhood adversity remains poorly understood. The aim of this study is to analyse the putative interaction between the 5-HTT gene (5-HTTLPR polymorphism), BDNF gene (Val66Met polymorphism) and childhood adversity in accounting for adult depressive symptoms.

Functional and evolutionary analysis of DXL1, a non-essential gene encoding a 1-deoxy-D-xylulose 5-phosphate synthase like protein in Arabidopsis thaliana

The synthesis of 1-deoxy-D-xylulose 5-phosphate (DXP), catalyzed by the enzyme DXP synthase (DXS), represents a key regulatory step of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway for isoprenoid biosynthesis. In plants DXS is encoded by small multigene families that can be classified into, at least, three specialized subfamilies. Arabidopsis thaliana contains three genes encoding proteins with similarity to DXS, including the well-known DXS1/CLA1 gene, which clusters within subfamily I. The remaining proteins, initially named DXS2 and DXS3, have not yet been characterized. Here we report the expression and functional analysis of A. thaliana DXS2. Unexpectedly, the expression of DXS2 failed to rescue Escherichia coli and A. thaliana mutants defective in DXS activity. Coherently, we found that DXS activity was negligible in vitro, being renamed as DXL1 following recent nomenclature recommendation. DXL1 is targeted to plastids as DXS1, but shows a distinct expression pattern. The phenotypic analysis of a DXL1 defective mutant revealed that the function of the encoded protein is not essential for growth and development. Evolutionary analyses indicated that DXL1 emerged from DXS1 through a recent duplication apparently specific of the Brassicaceae lineage. Divergent selective constraints would have affected a significant fraction of sites after diversification of the paralogues. Furthermore, amino acids subjected to divergent selection and likely critical for functional divergence through the acquisition of a novel, although not yet known, biochemical function, were identified. Our results provide with the first evidences of functional specialization at both the regulatory and biochemical level within the plant DXS family.

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.

Síntesi de siRNAs (short interfering RNAs) modificats amb nucleobases 5-metil i 5-propinil pirimidíniques i avaluació de la seva activitat in vitro en cultius cel•lulars i de la seva estabilitat en sèrum humà. Disseny de siRNAs dirigits a la inhibició de l'expressió de mutacions de l'oncogen K-ras.

Projecte de recerca elaborat a partir d’una estada a la Stanford University, EEUU, entre 2007 i 2009. El present projecte es basa 1) en la síntesi de cadenes d'ARN dirigides a la inhibició de l'expressió gènica per un mecanisme d'ARN d'interferència (siRNAs o short interefering RNAs) i 2) en l'avaluació de l'activitat in vitro d'aquests oligonucleòtids en cultius cel•lulars. Concretament, la meva recerca ha estat enfocada principalment a l'estudi de cadenes de siRNA modificades amb nucleobases 5-metil i 5-propinil pirimidíniques. Es tractava d'avaluar l'efecte que exerceixen els factors estèrics en el major groove (solc major) dels siRNAs sobre la seva activitat biològica. En aquest sentit, he dut aterme síntesi de fosforamidits de nucleòsis pirimidínics modificats a la posició C-5 de la nucleobase. A continuació he incorporat aquestes unitats nucleosídiques en cadenes d'ARN emprant un sintetitzador d’ADN/ARN i he estudiat l'estabilitat dels corresponents dúplexs d'ARN mitjançant experiments de desnaturalització tèrmica. Finalment he dut a terme experiments d'inhibició de l'expressió gènica en cèl.lules HeLa per tal d'avaluar l'activitat biològia d'aquests siRNAs modificats. Els resultats d'aquests estudis han posat de manifest que la presència de grups voluminosos com el propinil a l'extrem 5' del dúplex de siRNA (definit per la cadena guia o antisense) influeix de forma molt negativa en la seva activitat biològica. En canvi, grups menys voluminosos com el metil hi influeixen positivament, de manera que algunes de les cadenes sintetitzades han resultat ser més actives que els corresponents siRNAs naturals (wild type siRNAs). A més, aquest tipus de modificació contribueix positivament en l'estabilitat de cadenes de siRNA en sèrum humà. Aquest treball ha estat publicat (Terrazas, M.; Kool, E.T. "Major Groove Modifications Improve siRNA Stability and Biological Activity" Nucleic Acids Res. 2009, in press).

Improving gene annotation using peptide mass spectrometry

Annotation of protein-coding genes is a key goal of genome sequencing projects. In spite of tremendous recent advances in computational gene finding, comprehensive annotation remains a challenge. Peptide mass spectrometry is a powerful tool for researching the dynamic proteome and suggests an attractive approach to discover and validate protein-coding genes. We present algorithms to construct and efficiently search spectra against a genomic database, with no prior knowledge of encoded proteins. By searching a corpus of 18.5 million tandem mass spectra (MS/MS) from human proteomic samples, we validate 39,000 exons and 11,000 introns at the level of translation. We present translation-level evidence for novel or extended exons in 16 genes, confirm translation of 224 hypothetical proteins, and discover or confirm over 40 alternative splicing events. Polymorphisms are efficiently encoded in our database, allowing us to observe variant alleles for 308 coding SNPs. Finally, we demonstrate the use of mass spectrometry to improve automated gene prediction, adding 800 correct exons to our predictions using a simple rescoring strategy. Our results demonstrate that proteomic profiling should play a role in any genome sequencing project.

Prominent use of distal 5’ transcription start sites and discovery of a large number of additional exons in ENCODE regions

This report presents systematic empirical annotation of transcript products from 399 annotated protein-coding loci across the 1% of the human genome targeted by the Encyclopedia of DNA elements (ENCODE) pilot project using a combination of 5' rapid amplification of cDNA ends (RACE) and high-density resolution tiling arrays. We identified previously unannotated and often tissue- or cell-line-specific transcribed fragments (RACEfrags), both 5' distal to the annotated 5' terminus and internal to the annotated gene bounds for the vast majority (81.5%) of the tested genes. Half of the distal RACEfrags span large segments of genomic sequences away from the main portion of the coding transcript and often overlap with the upstream-annotated gene(s). Notably, at least 20% of the resultant novel transcripts have changes in their open reading frames (ORFs), most of them fusing ORFs of adjacent transcripts. A significant fraction of distal RACEfrags show expression levels comparable to those of known exons of the same locus, suggesting that they are not part of very minority splice forms. These results have significant implications concerning (1) our current understanding of the architecture of protein-coding genes; (2) our views on locations of regulatory regions in the genome; and (3) the interpretation of sequence polymorphisms mapping to regions hitherto considered to be "noncoding," ultimately relating to the identification of disease-related sequence alterations.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5′-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Gene finding in the chicken genome

Background: Despite the continuous production of genome sequence for a number of organisms,reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularlytrue for genomes for which there is not a large collection of known gene sequences, such as therecently published chicken genome. We used the chicken sequence to test comparative andhomology-based gene-finding methods followed by experimental validation as an effective genomeannotation method.Results: We performed experimental evaluation by RT-PCR of three different computational genefinders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram wascomputed and each component of it was evaluated. The results showed that de novo comparativemethods can identify up to about 700 chicken genes with no previous evidence of expression, andcan correctly extend about 40% of homology-based predictions at the 5' end.Conclusions: De novo comparative gene prediction followed by experimental verification iseffective at enhancing the annotation of the newly sequenced genomes provided by standardhomology-based methods.

Death and resurrection of the human IRGM gene

Immunity-related GTPases (IRG) play an important role in defense against intracellular pathogens. One member of this gene family in humans, IRGM, has been recently implicated as a risk factor for Crohn's disease. We analyzed the detailed structure of this gene family among primates and showed that most of the IRG gene cluster was deleted early in primate evolution, after the divergence of the anthropoids from prosimians ( about 50 million years ago). Comparative sequence analysis of New World and Old World monkey species shows that the single-copy IRGM gene became pseudogenized as a result of an Alu retrotransposition event in the anthropoid common ancestor that disrupted the open reading frame (ORF). We find that the ORF was reestablished as a part of a polymorphic stop codon in the common ancestor of humans and great apes. Expression analysis suggests that this change occurred in conjunction with the insertion of an endogenous retrovirus, which altered the transcription initiation, splicing, and expression profile of IRGM. These data argue that the gene became pseudogenized and was then resurrected through a series of complex structural events and suggest remarkable functional plasticity where alleles experience diverse evolutionary pressures over time. Such dynamism in structure and evolution may be critical for a gene family locked in an arms race with an ever-changing repertoire of intracellular parasites.

Comparative genomics of the odorant-binding and chemosensory protein gene families across the arthropoda: origin and evolutionary history of the chemosensory system

Chemoreception is a biological process essential for the survival of animals, as it allows the recognition of important volatile cues for the detection of food, egg-laying substrates, mates or predators, among other purposes. Furthermore, its role in pheromone detection may contribute to evolutionary processes such as reproductive isolation and speciation. This key role in several vital biological processes makes chemoreception a particularly interesting system for studying the role of natural selection in molecular adaptation. Two major gene families are involved in the perireceptor events of the chemosensory system: the odorant-binding protein (OBP) and chemosensory protein (CSP) families. Here, we have conducted an exhaustive comparative genomic analysis of these gene families in twenty Arthropoda species. We show that the evolution of the OBP and CSP gene families is highly dynamic, with a high number of gains and losses of genes, pseudogenes and independent origins of subfamilies. Taken together, our data clearly support the birth-and-death model for the evolution of these gene families with an overall high gene-turnover rate. Moreover, we show that the genome organization of the two families is significantly more clustered than expected by chance and, more important, that this pattern appears to be actively maintained across the Drosophila phylogeny. Finally, we suggest the homologous nature of the OBP and CSP gene families, dating back their MRCA (most recent common ancestor) to 380¿420 Mya, and we propose a scenario for the origin and diversification of these families.

Genetic variants of the FADS gene cluster and ELOVL gene family, colostrums LC-PUFA levels, breastfeeding, and child cognition

Introduction: Breastfeeding effects on cognition are attributed to long-chain polyunsaturated fatty acids (LC-PUFAs), but controversy persists. Genetic variation in fatty acid desaturase (FADS) and elongase (ELOVL) enzymes has been overlooked when studying the effects of LC-PUFAs supply on cognition. We aimed to: 1) to determine whether maternal genetic variants in the FADS cluster and ELOVL genes contribute to differences in LC-PUFA levels in colostrum; 2) to analyze whether these maternal variants are related to child cognition; and 3) to assess whether children's variants modify breastfeeding effects on cognition. Methods: Data come from two population-based birth cohorts (n = 400 mother-child pairs from INMA-Sabadell; and n = 340 children from INMA-Menorca). LC-PUFAs were measured in 270 colostrum samples from INMA-Sabadell. Tag SNPs were genotyped both in mothers and children (13 in the FADS cluster, 6 in ELOVL2, and 7 in ELOVL5). Child cognition was assessed at 14 mo and 4 y using the Bayley Scales of Infant Development and the McCarthy Scales of Children"s Abilities, respectively. Results: Children of mothers carrying genetic variants associated with lower FADS1 activity (regulating AA and EPA synthesis), higher FADS2 activity (regulating DHA synthesis), and with higher EPA/AA and DHA/AA ratios in colostrum showed a significant advantage in cognition at 14 mo (3.5 to 5.3 points). Not being breastfed conferred an 8- to 9-point disadvantage in cognition among children GG homozygote for rs174468 (low FADS1 activity) but not among those with the A allele. Moreover, not being breastfed resulted in a disadvantage in cognition (5 to 8 points) among children CC homozygote for rs2397142 (low ELOVL5 activity), but not among those carrying the G allele. Conclusion: Genetically determined maternal supplies of LC-PUFAs during pregnancy and lactation appear to be crucial for child cognition. Breastfeeding effects on cognition are modified by child genetic variation in fatty acid desaturase and elongase enzymes.

Gene-disease network analysis reveals functional modules in mendelian, complex and environmental diseases

AbstractBACKGROUND: Scientists have been trying to understand the molecular mechanisms of diseases to design preventive and therapeutic strategies for a long time. For some diseases, it has become evident that it is not enough to obtain a catalogue of the disease-related genes but to uncover how disruptions of molecular networks in the cell give rise to disease phenotypes. Moreover, with the unprecedented wealth of information available, even obtaining such catalogue is extremely difficult.PRINCIPAL FINDINGS: We developed a comprehensive gene-disease association database by integrating associations from several sources that cover different biomedical aspects of diseases. In particular, we focus on the current knowledge of human genetic diseases including mendelian, complex and environmental diseases. To assess the concept of modularity of human diseases, we performed a systematic study of the emergent properties of human gene-disease networks by means of network topology and functional annotation analysis. The results indicate a highly shared genetic origin of human diseases and show that for most diseases, including mendelian, complex and environmental diseases, functional modules exist. Moreover, a core set of biological pathways is found to be associated with most human diseases. We obtained similar results when studying clusters of diseases, suggesting that related diseases might arise due to dysfunction of common biological processes in the cell.CONCLUSIONS: For the first time, we include mendelian, complex and environmental diseases in an integrated gene-disease association database and show that the concept of modularity applies for all of them. We furthermore provide a functional analysis of disease-related modules providing important new biological insights, which might not be discovered when considering each of the gene-disease association repositories independently. Hence, we present a suitable framework for the study of how genetic and environmental factors, such as drugs, contribute to diseases.AVAILABILITY: The gene-disease networks used in this study and part of the analysis are available at http://ibi.imim.es/DisGeNET/DisGeNETweb.html#Download

Gene Promoter Evolution Targets the Center of the Human Protein Interaction Network

Assessing the contribution of promoters and coding sequences to gene evolution is an important step toward discovering the major genetic determinants of human evolution. Many specific examples have revealed the evolutionary importance of cis-regulatory regions. However, the relative contribution of regulatory and coding regions to the evolutionary process and whether systemic factors differentially influence their evolution remains unclear. To address these questions, we carried out an analysis at the genome scale to identify signatures of positive selection in human proximal promoters. Next, we examined whether genes with positively selected promoters (Prom+ genes) show systemic differences with respect to a set of genes with positively selected protein-coding regions (Cod+ genes). We found that the number of genes in each set was not significantly different (8.1% and 8.5%, respectively). Furthermore, a functional analysis showed that, in both cases, positive selection affects almost all biological processes and only a few genes of each group are located in enriched categories, indicating that promoters and coding regions are not evolutionarily specialized with respect to gene function. On the other hand, we show that the topology of the human protein network has a different influence on the molecular evolution of proximal promoters and coding regions. Notably, Prom+ genes have an unexpectedly high centrality when compared with a reference distribution (P = 0.008, for Eigenvalue centrality). Moreover, the frequency of Prom+ genes increases from the periphery to the center of the protein network (P = 0.02, for the logistic regression coefficient). This means that gene centrality does not constrain the evolution of proximal promoters, unlike the case with coding regions, and further indicates that the evolution of proximal promoters is more efficient in the center of the protein network than in the periphery. These results show that proximal promoters have had a systemic contribution to human evolution by increasing the participation of central genes in the evolutionary process.

Id Gene regulation and function in the prosensory domains of the chicken inner ear: a link between Bmp Signaling and Atoh1

Bone morphogenetic proteins (Bmps) regulate the expression of the proneural gene Atoh1 and the generation of hair cells in the developing inner ear. The present work explored the role of Inhibitor of Differentiation genes (Id1-3) in this process. The results show that Id genes are expressed in the prosensory domains of the otic vesicle, along with Bmp4 and Bmp7. Those domains exhibit high levels of the phosphorylated form of Bmp-responding R-Smads (P-Smad1,5,8), and of Bmp-dependent Smad transcriptional activity as shown by the BRE-tk-EGFP reporter. Increased Bmp signaling induces the expression of Id1-3 along with the inhibition of Atoh1. Conversely, the Bmp antagonist Noggin or the Bmp-receptor inhibitor Dorsomorphin elicit opposite effects, indicating that Bmp signaling is necessary for Id expression and Atoh1 regulation in the otocyst. The forced expression of Id3 is sufficient to reduce Atoh1 expression and to prevent the expression of hair cell differentiation markers. Together, these results suggest that Ids are part of the machinery that mediates the regulation of hair cell differentiation exerted by Bmps. In agreement with that, during hair cell differentiation Bmp4 expression, P-Smad1,5,8 levels and Id expression are downregulated from hair cells. However, Ids are also downregulated from the supporting cells which contrarily to hair cells exhibit high levels of Bmp4 expression, P-Smad1,5,8, and BRE-tk-EGFP activity, suggesting that in these cells Ids escape from Bmp/Smad signaling. The differential regulation of Ids in time and space may underlie the multiple functions of Bmp signaling during sensory organ development.

Nicotine-induced antinociception, rewarding effects and physical dependence are decreased in mice lacking the preproenkephalin gene

It has been shown previously that the endogenous opioid system may be involved in the behavioral effects of nicotine. In the present study, the participation of endogenous enkephalins on nicotine responses has been investigated by using preproenkephalin knock-out mice. Acute nicotine-induced hypolocomotion remained unaffected in these mice. In contrast, antinociception elicited in the tail-immersion and hot-plate tests by acute nicotine administration was reduced in mutant animals. The rewarding properties of nicotine were then investigated using the place-conditioning paradigm. Nicotine induced a conditioned place preference in wild-type animals, but this effect was absent in knock-out mice. Accordingly, in vivo microdialysis studies revealed that the enhancement in dopamine extracellular levels in the nucleus accumbens induced by nicotine was also reduced in preproenkephalin-deficient mice. Finally, the somatic expression of the nicotine withdrawal syndrome precipitated in nicotine-dependent mice by mecamylamine was significantly attenuated in mutant animals. In summary, the present results indicate that endogenous opioid peptides derived from preproenkephalin are involved in the antinociceptive and rewarding properties of nicotine and participate in the expression of physical nicotine dependence.