Genome-wide meta-analysis of common variant differences between men and women.

The male-to-female sex ratio at birth is constant across world populations with an average of 1.06 (106 male to 100 female live births) for populations of European descent. The sex ratio is considered to be affected by numerous biological and environmental factors and to have a heritable component. The aim of this study was to investigate the presence of common allele modest effects at autosomal and chromosome X variants that could explain the observed sex ratio at birth. We conducted a large-scale genome-wide association scan (GWAS) meta-analysis across 51 studies, comprising overall 114 863 individuals (61 094 women and 53 769 men) of European ancestry and 2 623 828 common (minor allele frequency >0.05) single-nucleotide polymorphisms (SNPs). Allele frequencies were compared between men and women for directly-typed and imputed variants within each study. Forward-time simulations for unlinked, neutral, autosomal, common loci were performed under the demographic model for European populations with a fixed sex ratio and a random mating scheme to assess the probability of detecting significant allele frequency differences. We do not detect any genome-wide significant (P < 5 × 10(-8)) common SNP differences between men and women in this well-powered meta-analysis. The simulated data provided results entirely consistent with these findings. This large-scale investigation across ~115 000 individuals shows no detectable contribution from common genetic variants to the observed skew in the sex ratio. The absence of sex-specific differences is useful in guiding genetic association study design, for example when using mixed controls for sex-biased traits.

Analysis of TNFAIP3, a feedback inhibitor of nuclear factor-κB and the neighbor intergenic 6q23 region in rheumatoid arthritis susceptibility

INTRODUCTION Genome-wide association studies of rheumatoid arthritis (RA) have identified an association of the disease with a 6q23 region devoid of genes. TNFAIP3, an RA candidate gene, flanks this region, and polymorphisms in both the TNFAIP3 gene and the intergenic region are associated with systemic lupus erythematosus. We hypothesized that there is a similar association with RA, including polymorphisms in TNFAIP3 and the intergenic region. METHODS To test this hypothesis, we selected tag-single nucleotide polymorphisms (SNPs) in both loci. They were analyzed in 1,651 patients with RA and 1,619 control individuals of Spanish ancestry. RESULTS Weak evidence of association was found both in the 6q23 intergenic region and in the TNFAIP3 locus. The rs582757 SNP and a common haplotype in the TNFAIP3 locus exhibited association with RA. In the intergenic region, two SNPs were associated, namely rs609438 and rs13207033. The latter was only associated in patients with anti-citrullinated peptide antibodies. Overall, statistical association was best explained by the interdependent contribution of SNPs from the two loci TNFAIP3 and the 6q23 intergenic region. CONCLUSIONS Our data are consistent with the hypothesis that several RA genetic factors exist in the 6q23 region, including polymorphisms in the TNFAIP3 gene, like that previously described for systemic lupus erythematosus.

Genetic analysis of caveolin-1 and eNOS genes in colorectal cancer.

Caveolae are involved in physical compartmentalization between different groups of signaling events. Its main component, CAV1, modulates different pathways in cellular physiology. The emerging evidence pointing to the role of CAV1 in cancer led us to study whether different alleles of this gene are associated with colorectal cancer (CRC). Since one of the most characterized enzymes regulated by CAV1 is eNOS, we decided to include both genes in this study. We analyzed five SNPs in 360 unrelated CRC patients and 550 controls from the general population. Two of these SNPs were located within eNOS and three within the CAV1 gene. Although haplotype distribution was not associated with CRC, haplotype TiA (CAV1) was associated with familiar forms of CRC (p<0.05). This was especially evident in CRC antecedents and nuclear forms of CRC. If both CG (eNOS) and TiA (CAV1) haplotypes were taken together, this association increased in significance. Thus, we propose that CAV1, either alone or together with eNOS alleles, might modify CRC heritability.

Sequence analysis of coding DNA fragments of pfcrt and pfmdr-1 genes in Plasmodium falciparum isolates from Odisha, India

The global emergence and spread of malaria parasites resistant to antimalarial drugs is the major problem in malaria control. The genetic basis of the parasite's resistance to the antimalarial drug chloroquine (CQ) is well-documented, allowing for the analysis of field isolates of malaria parasites to address evolutionary questions concerning the origin and spread of CQ-resistance. Here, we present DNA sequence analyses of both the second exon of the Plasmodium falciparum CQ-resistance transporter (pfcrt) gene and the 5' end of the P. falciparum multidrug-resistance 1 (pfmdr-1) gene in 40 P. falciparum field isolates collected from eight different localities of Odisha, India. First, we genotyped the samples for the pfcrt K76T and pfmdr-1 N86Y mutations in these two genes, which are the mutations primarily implicated in CQ-resistance. We further analyzed amino acid changes in codons 72-76 of the pfcrt haplotypes. Interestingly, both the K76T and N86Y mutations were found to co-exist in 32 out of the total 40 isolates, which were of either the CVIET or SVMNT haplotype, while the remaining eight isolates were of the CVMNK haplotype. In total, eight nonsynonymous single nucleotide polymorphisms (SNPs) were observed, six in the pfcrt gene and two in the pfmdr-1 gene. One poorly studied SNP in the pfcrt gene (A97T) was found at a high frequency in many P. falciparum samples. Using population genetics to analyze these two gene fragments, we revealed comparatively higher nucleotide diversity in the pfcrt gene than in the pfmdr-1 gene. Furthermore, linkage disequilibrium was found to be tight between closely spaced SNPs of the pfcrt gene. Finally, both the pfcrt and the pfmdr-1 genes were found to evolve under the standard neutral model of molecular evolution.

Polymorphism analysis of the CTLA-4 gene in paracoccidioidomycosis patients

The CTLA-4 protein is expressed in activated T cells and plays an essential role in the immune response through its regulatory effect on T cell activation. Polymorphisms of the CTLA-4 gene have been correlated with autoimmune, neoplastic and infectious illnesses. This work aimed to verify possible associations between single nucleotide polymorphisms (SNPs) in CTLA-4, -318C/T in the promoter and +49A/G in exon 1 and paracoccidioidomycosis (PCM) caused by Paracoccidioides brasiliensis. For this purpose, 66 chronic form PCM patients and 76 healthy controls had their allele, genotype and haplotype frequencies determined. The genetic admixture structure of the patients and controls was evaluated to eliminate ancestral bias. The comparison of frequencies indicated no significant differences between patients and controls that could link the SNPs to PCM. Groups were admixture matched with no difference observed in population ancestry inference, indicating that the absence of association between CTLA-4 polymorphisms and PCM could not be attributed to ancestral bias. This study showed that there was no association between the CTLA-4 SNPs -318 and +49 and the resistance or susceptibility to PCM.

Comprehensive analysis of NRG1 common and rare variants in Hirschsprung patients.

Hirschsprung disease (HSCR, OMIM 142623) is a developmental disorder characterized by the absence of ganglion cells along variable lengths of the distal gastrointestinal tract, which results in tonic contraction of the aganglionic gut segment and functional intestinal obstruction. The RET proto-oncogene is the major gene for HSCR with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. Many other genes have been described to be associated with the pathology, as NRG1 gene (8p12), encoding neuregulin 1, which is implicated in the development of the enteric nervous system (ENS), and seems to contribute by both common and rare variants. Here we present the results of a comprehensive analysis of the NRG1 gene in the context of the disease in a series of 207 Spanish HSCR patients, by both mutational screening of its coding sequence and evaluation of 3 common tag SNPs as low penetrance susceptibility factors, finding some potentially damaging variants which we have functionally characterized. All of them were found to be associated with a significant reduction of the normal NRG1 protein levels. The fact that those mutations analyzed alter NRG1 protein would suggest that they would be related with HSCR disease not only in Chinese but also in a Caucasian population, which reinforces the implication of NRG1 gene in this pathology.

Validated SNPs for eGFR and their associations with albuminuria.

Albuminuria and reduced glomerular filtration rate are manifestations of chronic kidney disease (CKD) that predict end-stage renal disease, acute kidney injury, cardiovascular disease and death. We hypothesized that SNPs identified in association with the estimated glomerular filtration rate (eGFR) would also be associated with albuminuria. Within the CKDGen Consortium cohort (n= 31 580, European ancestry), we tested 16 eGFR-associated SNPs for association with the urinary albumin-to-creatinine ratio (UACR) and albuminuria [UACR >25 mg/g (women); 17 mg/g (men)]. In parallel, within the CARe Renal Consortium (n= 5569, African ancestry), we tested seven eGFR-associated SNPs for association with the UACR. We used a Bonferroni-corrected P-value of 0.003 (0.05/16) in CKDGen and 0.007 (0.05/7) in CARe. We also assessed whether the 16 eGFR SNPs were associated with the UACR in aggregate using a beta-weighted genotype score. In the CKDGen Consortium, the minor A allele of rs17319721 in the SHROOM3 gene, known to be associated with a lower eGFR, was associated with lower ln(UACR) levels (beta = -0.034, P-value = 0.0002). No additional eGFR-associated SNPs met the Bonferroni-corrected P-value threshold of 0.003 for either UACR or albuminuria. In the CARe Renal Consortium, there were no associations between SNPs and UACR with a P< 0.007. Although we found the genotype score to be associated with albuminuria (P= 0.0006), this result was driven almost entirely by the known SHROOM3 variant, rs17319721. Removal of rs17319721 resulted in a P-value 0.03, indicating a weak residual aggregate signal. No alleles, previously demonstrated to be associated with a lower eGFR, were associated with the UACR or albuminuria, suggesting that there may be distinct genetic components for these traits.

Fine mapping and functional analysis of the multiple sclerosis risk gene CD6.

CD6 has recently been identified and validated as risk gene for multiple sclerosis (MS), based on the association of a single nucleotide polymorphism (SNP), rs17824933, located in intron 1. CD6 is a cell surface scavenger receptor involved in T-cell activation and proliferation, as well as in thymocyte differentiation. In this study, we performed a haptag SNP screen of the CD6 gene locus using a total of thirteen tagging SNPs, of which three were non-synonymous SNPs, and replicated the recently reported GWAS SNP rs650258 in a Spanish-Basque collection of 814 controls and 823 cases. Validation of the six most strongly associated SNPs was performed in an independent collection of 2265 MS patients and 2600 healthy controls. We identified association of haplotypes composed of two non-synonymous SNPs [rs11230563 (R225W) and rs2074225 (A257V)] in the 2(nd) SRCR domain with susceptibility to MS (P max(T) permutation = 1×10(-4)). The effect of these haplotypes on CD6 surface expression and cytokine secretion was also tested. The analysis showed significantly different CD6 expression patterns in the distinct cell subsets, i.e. - CD4(+) naïve cells, P = 0.0001; CD8(+) naïve cells, P<0.0001; CD4(+) and CD8(+) central memory cells, P = 0.01 and 0.05, respectively; and natural killer T (NKT) cells, P = 0.02; with the protective haplotype (RA) showing higher expression of CD6. However, no significant changes were observed in natural killer (NK) cells, effector memory and terminally differentiated effector memory T cells. Our findings reveal that this new MS-associated CD6 risk haplotype significantly modifies expression of CD6 on CD4(+) and CD8(+) T cells.

Analysis of the association between CD40 and CD40 ligand polymorphisms and systemic sclerosis.

INTRODUCTION The aim of the present study was to investigate the possible role of CD40 and CD40 ligand (CD40LG) genes in the susceptibility and phenotype expression of systemic sclerosis (SSc). METHODS In total, 2,670 SSc patients and 3,245 healthy individuals from four European populations (Spain, Germany, The Netherlands, and Italy) were included in the study. Five single-nucleotide polymorphisms (SNPs) of CD40 (rs1883832, rs4810485, rs1535045) and CD40LG (rs3092952, rs3092920) were genotyped by using a predesigned TaqMan allele-discrimination assay technology. Meta-analysis was assessed to determine whether an association exists between the genetic variants and SSc or its main clinical subtypes. RESULTS No evidence of association between CD40 and CD40LG genes variants and susceptibility to SSc was observed. Similarly, no significant statistical differences were observed when SSc patients were stratified by the clinical subtypes, the serologic features, and pulmonary fibrosis. CONCLUSIONS Our results do not suggest an important role of CD40 and CD40LG gene polymorphisms in the susceptibility to or clinical expression of SSc.

Silver nanoparticle production by the fungus Fusarium oxysporum: nanoparticle characterisation and analysis of antifungal activity against pathogenic yeasts

The microbial synthesis of nanoparticles is a green chemistry approach that combines nanotechnology and microbial biotechnology. The aim of this study was to obtain silver nanoparticles (SNPs) using aqueous extract from the filamentous fungus Fusarium oxysporum as an alternative to chemical procedures and to evaluate its antifungal activity. SNPs production increased in a concentration-dependent way up to 1 mM silver nitrate until 30 days of reaction. Monodispersed and spherical SNPs were predominantly produced. After 60 days, it was possible to observe degenerated SNPs with in additional needle morphology. The SNPs showed a high antifungal activity against Candida and Cryptococcus , with minimum inhibitory concentration values ≤ 1.68 µg/mL for both genera. Morphological alterations of Cryptococcus neoformans treated with SNPs were observed such as disruption of the cell wall and cytoplasmic membrane and lost of the cytoplasm content. This work revealed that SNPs can be easily produced by F. oxysporum aqueous extracts and may be a feasible, low-cost, environmentally friendly method for generating stable and uniformly sized SNPs. Finally, we have demonstrated that these SNPs are active against pathogenic fungi, such as Candida and Cryptococcus .

FGF21 signalling pathway and metabolic traits - genetic association analysis.

Fibroblast growth factor 21 (FGF21) is a novel master regulator of metabolic profile. The biological actions of FGF21 are elicited upon its klotho beta (KLB)-facilitated binding to FGF receptor 1 (FGFR1), FGFR2 and FGFR3. We hypothesised that common polymorphisms in the FGF21 signalling pathway may be associated with metabolic risk. At the screening stage, we examined associations between 63 common single-nucleotide polymorphisms (SNPs) in five genes of this pathway (FGF21, KLB, FGFR1, FGFR2, FGFR3) and four metabolic phenotypes (LDL cholesterol - LDL-C, HDL-cholesterol - HDL-C, triglycerides and body mass index) in 629 individuals from Silesian Hypertension Study (SHS). Replication analyses were performed in 5478 unrelated individuals of the Swiss CoLaus cohort (imputed genotypes) and in 3030 directly genotyped individuals of the German Myocardial Infarction Family Study (GerMIFS). Of 54 SNPs that met quality control criteria after genotyping in SHS, 4 (rs4733946 and rs7012413 in FGFR1; rs2071616 in FGFR2 and rs7670903 in KLB) showed suggestive association with LDL-C (P=0.0006, P=0.0013, P=0.0055, P=0.011, respectively) and 1 (rs2608819 in KLB) was associated with body mass index (P=0.011); all with false discovery rate q<0.5. Of these, only one FGFR2 polymorphism (rs2071616) showed replicated association with LDL-C in both CoLaus (P=0.009) and men from GerMIFS (P=0.017). The direction of allelic effect of rs2071616 upon LDL-C was consistent in all examined populations. These data show that common genetic variations in FGFR2 may be associated with LDL-C in subjects of white European ancestry.

A genomic analysis identifies a novel component in the genetic structure of sub-Saharan African populations

Studies of large sets of SNP data have proven to be a powerful tool in the analysis of the genetic structure of human populations. In this work, we analyze genotyping data for 2,841 SNPs in 12 Sub-Saharan African populations, including a previously unsampled region of south-eastern Africa (Mozambique). We show that robust results in a world-wide perspective can be obtained when analyzing only 1,000 SNPs. Our main results both confirm the results of previous studies, and show new and interesting features in Sub-Saharan African genetic complexity. There is a strong differentiation of Nilo-Saharans, much beyond what would be expected by geography. Hunter-gatherer populations (Khoisan and Pygmies) show a clear distinctiveness with very intrinsic Pygmy (and not only Khoisan) genetic features. Populations of the West Africa present an unexpected similarity among them, possibly the result of a population expansion. Finally, we find a strong differentiation of the south-eastern Bantu population from Mozambique, which suggests an assimilation of a pre-Bantu substrate by Bantu speakers in the region.

From SNPs to pathways: integration of functional effect of sequence variations on models of cell signalling pathways

Background: Single nucleotide polymorphisms (SNPs) are the most frequent type of sequence variation between individuals, and represent a promising tool for finding genetic determinants of complex diseases and understanding the differences in drug response. In this regard, it is of particular interest to study the effect of non-synonymous SNPs in the context of biological networks such as cell signalling pathways. UniProt provides curated information about the functional and phenotypic effects of sequence variation, including SNPs, as well as on mutations of protein sequences. However, no strategy has been developed to integrate this information with biological networks, with the ultimate goal of studying the impact of the functional effect of SNPs in the structure and dynamics of biological networks. Results: First, we identified the different challenges posed by the integration of the phenotypic effect of sequence variants and mutations with biological networks. Second, we developed a strategy for the combination of data extracted from public resources, such as UniProt, NCBI dbSNP, Reactome and BioModels. We generated attribute files containing phenotypic and genotypic annotations to the nodes of biological networks, which can be imported into network visualization tools such as Cytoscape. These resources allow the mapping and visualization of mutations and natural variations of human proteins and their phenotypic effect on biological networks (e.g. signalling pathways, protein-protein interaction networks, dynamic models). Finally, an example on the use of the sequence variation data in the dynamics of a network model is presented. Conclusion: In this paper we present a general strategy for the integration of pathway and sequence variation data for visualization, analysis and modelling purposes, including the study of the functional impact of protein sequence variations on the dynamics of signalling pathways. This is of particular interest when the SNP or mutation is known to be associated to disease. We expect that this approach will help in the study of the functional impact of disease-associated SNPs on the behaviour of cell signalling pathways, which ultimately will lead to a better understanding of the mechanisms underlying complex diseases.

Association between variants of the leptin receptor gene (LEPR) and overweight: a systematic review and an analysis of the CoLaus study.

BACKGROUND: Three non-synonymous single nucleotide polymorphisms (Q223R, K109R and K656N) of the leptin receptor gene (LEPR) have been tested for association with obesity-related outcomes in multiple studies, showing inconclusive results. We performed a systematic review and meta-analysis on the association of the three LEPR variants with BMI. In addition, we analysed 15 SNPs within the LEPR gene in the CoLaus study, assessing the interaction of the variants with sex. METHODOLOGY/PRINCIPAL FINDINGS: We searched electronic databases, including population-based studies that investigated the association between LEPR variants Q223R, K109R and K656N and obesity- related phenotypes in healthy, unrelated subjects. We furthermore performed meta-analyses of the genotype and allele frequencies in case-control studies. Results were stratified by SNP and by potential effect modifiers. CoLaus data were analysed by logistic and linear regressions and tested for interaction with sex. The meta-analysis of published data did not show an overall association between any of the tested LEPR variants and overweight. However, the choice of a BMI cut-off value to distinguish cases from controls was crucial to explain heterogeneity in Q223R. Differences in allele frequencies across ethnic groups are compatible with natural selection of derived alleles in Q223R and K109R and of the ancient allele in K656N in Asians. In CoLaus, the rs10128072, rs3790438 and rs3790437 variants showed interaction with sex for their association with overweight, waist circumference and fat mass in linear regressions. CONCLUSIONS: Our systematic review and analysis of primary data from the CoLaus study did not show an overall association between LEPR SNPs and overweight. Most studies were underpowered to detect small effect sizes. A potential effect modification by sex, population stratification, as well as the role of natural selection should be addressed in future genetic association studies.

Causal relationship between obesity and vitamin D status: bi-directional Mendelian randomization analysis of multiple cohorts.

BACKGROUND: Obesity is associated with vitamin D deficiency, and both are areas of active public health concern. We explored the causality and direction of the relationship between body mass index (BMI) and 25-hydroxyvitamin D [25(OH)D] using genetic markers as instrumental variables (IVs) in bi-directional Mendelian randomization (MR) analysis. METHODS AND FINDINGS: We used information from 21 adult cohorts (up to 42,024 participants) with 12 BMI-related SNPs (combined in an allelic score) to produce an instrument for BMI and four SNPs associated with 25(OH)D (combined in two allelic scores, separately for genes encoding its synthesis or metabolism) as an instrument for vitamin D. Regression estimates for the IVs (allele scores) were generated within-study and pooled by meta-analysis to generate summary effects. Associations between vitamin D scores and BMI were confirmed in the Genetic Investigation of Anthropometric Traits (GIANT) consortium (n = 123,864). Each 1 kg/m(2) higher BMI was associated with 1.15% lower 25(OH)D (p = 6.52×10⁻²⁷). The BMI allele score was associated both with BMI (p = 6.30×10⁻⁶²) and 25(OH)D (-0.06% [95% CI -0.10 to -0.02], p = 0.004) in the cohorts that underwent meta-analysis. The two vitamin D allele scores were strongly associated with 25(OH)D (p≤8.07×10⁻⁵⁷ for both scores) but not with BMI (synthesis score, p = 0.88; metabolism score, p = 0.08) in the meta-analysis. A 10% higher genetically instrumented BMI was associated with 4.2% lower 25(OH)D concentrations (IV ratio: -4.2 [95% CI -7.1 to -1.3], p = 0.005). No association was seen for genetically instrumented 25(OH)D with BMI, a finding that was confirmed using data from the GIANT consortium (p≥0.57 for both vitamin D scores). CONCLUSIONS: On the basis of a bi-directional genetic approach that limits confounding, our study suggests that a higher BMI leads to lower 25(OH)D, while any effects of lower 25(OH)D increasing BMI are likely to be small. Population level interventions to reduce BMI are expected to decrease the prevalence of vitamin D deficiency.