Search PubMed⌕ Search

Biomedical subjects

Javier Costas

Publications and source records attributed to Javier Costas.

12 recordsLinked to original sources

Testing for Genetic Interactions in Complex Disease With Distance Correlation.

Understanding epistasis (genetic interaction) may shed some light on the genomic basis of common diseases, including disorders of maximum interest due to their high socioeconomic burden, like schizophrenia. Distance correlation is an association measure that characterizes general statistical independence between random variables, not only the linear one. Here, we propose distance correlation as a novel tool for the detection of epistasis from case-control data of single-nucleotide polymorphisms. On the methodological side, we highlight the derivation of the explicit asymptotic null distribution of the test statistic. We show that this is the only way to obtain enough computational speed for the method to be used in practice, in a scenario where the resampling techniques found in the literature are impractical. Our simulations show satisfactory calibration of significance, as well as comparable or better power than existing methodology. We conclude with the application of our technique to a schizophrenia genetics dataset, obtaining biologically sound insights.

Epistasis, Genetic↗

Genetic and epigenetic changes to the glucocorticoid receptor gene (NR3C1) and cognition in major depressive disorder.

INTRODUCTION: Many studies have found that hypothalamic-pituitary-adrenal (HPA) axis abnormalities are related to the pathophysiology of major depressive disorder (MDD) and cognitive functioning. Our aim was to assess the influence of genetic polymorphisms and methylation levels in three different promoter regions throughout the glucocorticoid receptor (GR) gene NR3C1 on cognitive performance in MDD. Plausible interactions with childhood adversity and mediation relationships between genetic and epigenetic variables were explored. MATERIALS AND METHODS: The sample included a total of 64 MDD patients and 82 healthy controls. Child maltreatment and neurocognitive performance were assessed in all participants. HPA negative feedback was analyzed using the dexamethasone suppression test after the administration of 0.25mg of dexamethasone. A total of 23 single-nucleotide polymorphisms were genotyped, and methylation levels at several CpGs in exons 1D, 1F and 1H of the GR gene were measured. RESULTS: Results show that, beyond the influence of other covariables, NR3C1 single-nucleotide polymorphisms and methylation levels predicted performance in executive functioning and working memory tasks. No significant interactions or mediation relationships were detected. CONCLUSIONS: Results suggest that genetic variations and epigenetic regulation of the GR gene are relevant factors influencing cognitive performance in MDD and could emerge as significant biomarkers and therapeutic targets in mood disorders and other stress-related disorders.

Humans↗

Extensive linkage disequilibrium mapping at HTR2A and DRD3 for schizophrenia susceptibility genes in the Galician population.

The serotonin and dopamine neurotransmitter systems are candidate pathways in the development of schizophrenia because of the assumed causal relationship with the observed symptoms as well as effective targeting of the corresponding receptors by antipsychotic drugs. However, genetic association studies have systematically focused on a limited set of genes and single nucleotide polymorphisms (SNPs), including T102C at HTR2A and Ser9Gly at DRD3. Meta-analyses of the associations between these two markers and schizophrenia revealed a true increase in risk, the magnitude of the effect being very low. In the present study we analyzed 260 schizophrenic patients and 354 control subjects from a homogeneous population, the Galician population, using an extensive linkage disequilibrium (LD) mapping approach, genotyping a total of 47 SNPs to test for the existence of additional variants that confer higher risk. We detected nominal significant association with schizophrenia for several haplotype tag SNPs (htSNPs) at HTR2A, although the significance was lost after multiple test corrections. In addition, haplotype analyses involving a sliding window approach, with window size 2 to 4 SNPs, revealed significant differences in frequencies of the DRD3 haplotypes at the 3' half of the gene region. This difference, which remains clearly significant after multiple test corrections (p=0.002, 0.0001, and 0.0025, for window sizes 2, 3, and 4, respectively), was mainly due to over-representation of several rare haplotypes in patients, at the expense of a single common haplotype; this represents interesting evidence of rare haplotypes for susceptibility detected using common htSNPs due to their strong effect.

3' Flanking Region↗

High variability in CYP21A2 mutated alleles in Spanish 21-hydroxylase deficiency patients, six novel mutations and a founder effect.

OBJECTIVE: To detect common as well as rare and novel CYP21A mutations in 21-hydroxylase deficiency patients. To estimate the distribution of mutations and compare them with other European studies. To construct haplotypes linked to a recurrent novel mutation. DESIGN: Genetic analysis by sequencing the entire CYP21A2 gene plus Southern blot. PATIENTS: A total of 138 unrelated Spanish patients: 122 nonclassical forms (NCF) and 16 classical forms (CF) were studied. RESULTS: Among the 266 nonrelated mutated alleles; CYP21A2 deletions/conversions and a spectrum of 27 different mutated alleles were found: 15 different single point mutations, 8 nucleotide deletions in exon 3, 3 mutation clusters in exon 6, 9 alleles with more than one mutation, one 21-nucleotide duplication in exon 10, and one allele with CYP21A2 duplicated and both copies mutated. The most frequent mutation in NCF alleles is V281L (71.8%). Among CFs, the most common is I2 g (20%) and Q318X (16%) and rare alleles (21.9%). Six novel causative mutations were found, four associated with CF: I46+1nt, R444X, P463L and M473_R479dup and two associated with NCF: W302 and D322G. The R444X mutation was found in seven unrelated patients and it appeared only once in an ancestral haplotype. In addition, we found a novel single nucleotide polymorphism with a 31.5% frequency for the rare allele. CONCLUSION: A great diversity of haplotypes with a large spectrum of mutated alleles was found. The frequency of the V281L mutation was the highest reported and the relatively high frequency of R444X was the result of a founder effect.

Adrenal Hyperplasia, Congenital↗

SUMO4 and MAP3K7IP2 single nucleotide polymorphisms and susceptibility to rheumatoid arthritis.

OBJECTIVE: To explore the role of single nuclear polymorphisms (SNP) in 2 candidate genes, SUMO4 and MAP3K7IP2, in susceptibility to rheumatoid arthritis (RA). METHODS: Two cohorts from different Spanish towns totalling 635 patients with RA and 826 controls were studied. Six SNP were genotyped by matrix assisted laser desorption-ionization time-of-flight (MALDI-TOF) with the MassARRAY SNP genotyping system. RESULTS: We found no association with susceptibility to RA for any of the SNP including a previously described functional variant in the SUMO4 gene (163A-->G). RA susceptibility was independent of the haplotypes defined by the 6 SNP and there was also no association with clinical features of RA. Conclusion. SUMO4 and MAP3K7IP2 SNP did not significantly influence predisposition to and features of RA, in contrast to previous genetic and functional evidence that suggested their involvement.

Adaptor Proteins, Signal Transducing↗

Human genome-wide screen of haplotype-like blocks of reduced diversity.

An important proportion of the human genome is organized in regions of high linkage disequilibrium (LD) and low haplotype diversity, referred to as haplotype blocks. Here, we perform a genome-wide screen of haplotype-like blocks presenting just two main haplotypes at a frequency higher than 1%, based on single-nucleotide polymorphism (SNP) frequencies from two populations: African-Americans and Caucasians, using data from the Celera SNP database. These haplotype-like blocks of reduced diversity are more abundant and of longer size in Caucasians, in agreement with population history. Several of the discovered blocks are good candidates for targets of natural selection, such as those blocks containing a cluster of bitter taste receptors or the apolipoprotein L1. In addition, several genes putatively involved in susceptibility to common diseases are included in these haplotype-like blocks of reduced diversity. This fact may present important implications in association studies, leading to a reduction of genotyping efforts.

Black People↗

Relative efficiency of the linkage disequilibrium mapping approach in detecting candidate genes for schizophrenia in different European populations.

Detection of susceptibility genes in indirect association studies depends not only on the degree of linkage disequilibrium between the disease variant and the SNP marker but also on the difference in their allele frequencies. Little is known about how variations in these parameters may affect the power of indirect association studies among related populations. Toward this goal, we genotyped 40 SNPs at four loci in samples from three European populations, Galician, Greek, and Norwegian. We compared the relative efficiency of all pairs of SNPs in detecting each other in each one of the populations. Our results show that a low percentage of marker SNPs may detect association in some populations but be totally ineffective in others. Therefore, these differences have to be an additional factor to consider when a replication study fails to confirm initial associations, especially if the replication is focused on very few markers.

Chromosome Mapping↗

Structural and evolutionary analyses of the Ty3/gypsy group of LTR retrotransposons in the genome of Anopheles gambiae.

The recent availability of the genome of Anopheles gambiae offers an extraordinary opportunity for comparative studies of the diversity of transposable elements (TEs) and their evolutionary dynamics between two related species, taking advantage of the existing information from Drosophila melanogaster. To this goal, we screened the genome of A. gambiae for elements belonging to the Ty3/gypsy group of long-terminal repeat (LTR) retrotransposons. The A. gambiae genome displays a rich diversity of LTR retrotransposons, clearly greater than D. melanogaster. We have characterized in detail 63 families, belonging to five of the nine main lineages of the Ty3/gypsy group. The Mag lineage is the most diverse and abundant, with more than 30 families. In sharp contrast with this finding, a single family belonging to this lineage has been found in D. melanogaster, here reported for the first time in the literature, most probably consisting of old inactive elements. The CsRn1 lineage is also abundant in A. gambiae but almost absent from D. melanogaster. Conversely, the Osvaldo lineage has been detected in Drosophila but not in Anopheles. Comparison of structural characteristics of different families led to the identification of several lineage-specific features such as the primer-binding site (PBS), the gag-pol translational recoding signal (TRS), which is extraordinarily diverse within the Ty3/gypsy retrotransposons of A. gambiae, or the presence/absence of specific amino acid motifs. Interestingly, some of these characteristics, although in general well conserved within lineages, may have evolved independently in particular branches of the phylogenetic tree. We also show evidence of recent activity for around 75% of the families. Nevertheless, almost all families contain a high proportion of degenerate members and solitary LTRs (solo LTRs), indicative of a lower turnover rate of retrotransposons belonging to the Ty3/gypsy group in A. gambiae than in D. melanogaster. Finally, we have detected significant overrepresentations of insertions on the X chromosome versus autosomes and of putatively active insertions on euchromatin versus heterochromatin.

Amino Acid Motifs↗

Genomic characterization of a repetitive motif strongly associated with developmental genes in Drosophila.

BACKGROUND: Non-coding DNA represents a high proportion of all metazoan genomes. Although an undetermined fraction of this DNA may be considered devoid of any function, it also contains important information residing in specific cis-regulatory sequences. RESULTS: We report a 27 bp motif that is overrepresented within the fly genome. This motif does not show any significant similarity with transposon sequences and is strongly associated with genes involved in development and/or signal transduction. The 27 bp motif is preferentially located within introns, and has a tendency to be present in multiple copies around genes. Furthermore, it is often found embedded in known non-coding regulatory regions. The regulatory network defined by this motif is partially shared in D. pseudoobscura. CONCLUSION: We have identified a 27 bp cis-regulatory sequence widely distributed within the Drosophila genome in association with developmental genes. This motif may be very useful towards the annotation of functional regulatory regions within the Drosophila genome and the construction of regulatory networks of Drosophila development.

Animals↗

Turnover of binding sites for transcription factors involved in early Drosophila development.

Despite the importance of cis-regulatory regions in evolution, little is know about their evolutionary dynamics. In this report, we analyze the process of evolution of binding sites for transcription factors using as a model a well characterized system, the Drosophila early developmental enhancers. We compare the sequences of eight enhancer regions for early developmental genes between Drosophila melanogaster and other two species, Drosophila virilis and Drosophila pseudoobscura, searching for the presence/absence of 104 biochemically verified binding sites from D. melanogaster. We also modeled the binding specificity of each binding site by the use of well-defined positional weight matrices (PWMs). The comparisons showed that turnover of binding sites seems to fit a molecular clock, at an approximate rate of 0.94% of gain/loss of binding sites per million years. This intense turnover affects both high and low affinity binding sites at the same extent. Furthermore, the subset of overlapping binding sites is also subjected to this high turnover. Conserved binding sites seem to be constrained to maintain not only location but also the exact sequence at each particular position. Finally, we detected a significant decrease in mean PWM scores for the D. virilis binding sites in the case of Hunchback. Possible explanations for this fact are discussed.

Animals↗

Molecular characterization of the recent intragenomic spread of the murine endogenous retrovirus MuERV-L.

The mouse genome has been subjected to two successive amplification bursts of the murine endogenous retrovirus MuERV-L after the Mus/Rattus split. The main objective of this work is to characterize in detail the intragenomic spread giving rise to these two murine bursts using full-length MuERV-L proviruses taken from public databases. Phylogenetic analyses led to the identification of elements putatively amplifying during each one of the two burst. Likelihood-ratio tests were used to confirm that elements supposedly arisen during the first burst have been evolving under lower selective constrains, as expected for older insertions. The data reported here suggested an evolutionary dynamics for MuERV-L amplification characterized by the existence of multiple elements simultaneously active during each one of the bursts while only one or very few closely related proviruses from the first burst gave rise to the second one. Finally, more than one third of the proviruses present 100% identity between the 5' and 3' LTRs, strongly indicating that MuERV-L is currently active within the mouse genome.

Animals↗

Characterization of the intragenomic spread of the human endogenous retrovirus family HERV-W.

This study examines the intragenomic spread of the human endogenous retrovirus family HERV-W from insertions present within the draft sequence of the human genome. Identification of shared diagnostic differences and phylogenetic analyses revealed the existence of three main subfamilies. The average divergence between sequences for each of the subfamilies suggests that most of the HERV-W elements were inserted within the genome during a short period of evolutionary time. Each one of the subfamilies consists of two types of insertions, the expected proviral sequences and other sequences resembling the structure of processed retrogenes. These HERV-W retrosequences extend from the R region of the 5' long-terminal repeat (LTR) to the R region of the 3' LTR (as viral genomic RNAs), end in poly(A) 3' tails, and are flanked by direct repeats longer than the proviral integrations. Furthermore, several of the HERV-W retrosequences are 5'-truncated at different sites. I suggest the involvement of the L1 machinery in these integrations and discuss the characteristic features of the evolutionary history of HERV-W, with emphasis on the putative impact of HERV-W retrosequence integrations on the mammalian genome.

Base Sequence↗