Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 x 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 x 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 x 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Significant reproductive skew in the facultatively polygynous ant Pheidole pallidula

Reproductive skew - the extent to which reproduction is unevenly shared between individuals in a social group - varies greatly between and within animal species. In this study, we investigated how queens share parentage in polygynous (multiple queen) colonies of the Mediterranean ant Pheidole pallidula. We used highly polymorphic microsatellites markers to determine parentage of gynes (new queens), males and workers in P. pallidula field colonies. The comparison of the genotypes of young and adult workers revealed a very low queen turnover (less than 2%). The first main finding of the study of reproductive skew in these colonies was that there was a significant departure from equal contribution of queens to gyne, male and worker production. Reproductive skew was greater for male production than for queen and worker production. There was no relationship between the magnitude of the reproductive skew and the number of reproductive queens per colony, their relatedness and the overall colony productivity, some of the factors predicted to influence the extent of reproductive skew. Finally, our study revealed for the first time a trade-off in the relative contribution of nestmate queens to gyne and worker production. The queens contributing more to gyne production contributed significantly less to worker production.

Genetic variants in IGF-I, IGF-II, IGFBP-3, and adiponectin genes and colon cancer risk in African Americans and Whites.

PURPOSE: Evaluating genetic susceptibility may clarify effects of known environmental factors and also identify individuals at high risk. We evaluated the association of four insulin-related pathway gene polymorphisms in insulin-like growth factor-1 (IGF-I) (CA)( n ) repeat, insulin-like growth factor-2 (IGF-II) (rs680), insulin-like growth factor-binding protein-3 (IGFBP-3) (rs2854744), and adiponectin (APM1 rs1501299) with colon cancer risk, as well as relationships with circulating IGF-I, IGF-II, IGFBP-3, and C-peptide in a population-based study. METHODS: Participants were African Americans (231 cases and 306 controls) and Whites (297 cases, 530 controls). Consenting subjects provided blood specimens and lifestyle/diet information. Genotyping for all genes except IGF-I was performed by the 5'-exonuclease (Taqman) assay. The IGF-I (CA)(n) repeat was assayed by PCR and fragment analysis. Circulating proteins were measured by enzyme immunoassays. Odds ratios (ORs) and 95 % confidence intervals (CIs) were calculated by logistic regression. RESULTS: The IGF-I (CA)( 19 ) repeat was higher in White controls (50 %) than African American controls (31 %). Whites homozygous for the IGF-I (CA)(19) repeat had a nearly twofold increase in risk of colon cancer (OR = 1.77; 95 % CI = 1.15-2.73), but not African Americans (OR = 0.73, 95 % CI 0.50-1.51). We observed an inverse association between the IGF-II Apa1 A-variant and colon cancer risk (OR = 0.49, 95 % CI 0.28-0.88) in Whites only. Carrying the IGFBP-3 variant alleles was associated with lower IGFBP-3 protein levels, a difference most pronounced in Whites (p-trend <0.05). CONCLUSIONS: These results support an association between insulin pathway-related genes and elevated colon cancer risk in Whites but not in African Americans.

The missense of smell: functional variability in the human odorant receptor repertoire.

Humans have ~400 intact odorant receptors, but each individual has a unique set of genetic variations that lead to variation in olfactory perception. We used a heterologous assay to determine how often genetic polymorphisms in odorant receptors alter receptor function. We identified agonists for 18 odorant receptors and found that 63% of the odorant receptors we examined had polymorphisms that altered in vitro function. On average, two individuals have functional differences at over 30% of their odorant receptor alleles. To show that these in vitro results are relevant to olfactory perception, we verified that variations in OR10G4 genotype explain over 15% of the observed variation in perceived intensity and over 10% of the observed variation in perceived valence for the high-affinity in vitro agonist guaiacol but do not explain phenotype variation for the lower-affinity agonists vanillin and ethyl vanillin.

Effect of interleukin-6 polymorphisms on human longevity:a systematic review and meta-analysis

Several studies have assessed changes in frequency of -174 interleukin (IL)-6 single nucleotide polymorphism (SNP) with age. If IL-6 tracks with disability and age-related diseases, then there should be reduction, in the oldest old, of the frequency of homozygous GG subjects, who produce higher IL-6 levels. However, discordant results have been obtained. To explore the relationship between this polymorphism and longevity, we analyzed individual data on long-living subjects and controls from eight case-control studies conducted in Europeans, using meta-analysis. There was no significant difference in the IL-6 genotype between the oldest old and controls (Odds Ratio [OR]=0.96; 95% C.I.: 0.77-1.20; p=0.71), but there was significant between-study heterogeneity (I2=55.5%). In a subgroup analyses when male centenarians from the three Italian studies were included, the frequency of the IL-6 -174 GG genotype was significantly lower than the other genotypes (OR=0.49; 95% C.I.: 0.31-0.80; p=0.004), with no evidence of heterogeneity (I2=0%). Our data supports a negative association between the GG genotype of IL-6 SNP and longevity in Italian centenarians, with males who carry the genotype being two times less likely to reach extreme old age compared with subjects carrying CC or CG genotypes. These findings were not replicated in other European groups suggesting a possible interaction between genetics, sex and environment in reaching longevity.

Catechol-O-methyltransferase and the clinical features of psychosis

A functional polymorphism (Val-158-Met) at the Catechol-O-methyltransferase (COMT) locus has been identified as a potential etiological factor in schizophrenia. Yet the association has not been convincingly replicated across independent samples. We hypothesized that phenotypic heterogeneity might be diluting the COMT effect. To clarify the putative association, we performed an exploratory analysis to test for association between COMT and five psychosis symptom scales. These were derived through factor analysis of the Operational Criteria Checklist for Psychiatric Illness. Our sample was the Irish Study of High Density Schizophrenia Families, a large collection consisting of 268 multiplex families. This sample has previously shown a small but significant effect of the COMT Val allele in conferring risk for schizophrenia. We tested for preferential transmission of COMT alleles from parent to affected offspring (n = 749) for each of the five factor-derived scales (negative symptoms, delusions, hallucinations, mania, and depression). Significant overtransmission of the Val allele was found for mania (P <0.05) and depression (P = 0.01) scales. Examination of odds ratios (ORs) revealed a heterogeneous effect of COMT, whereby it had no effect on Negative Symptoms, but largest impact on Depression (OR = 1.4). These results suggest a modest affective vulnerability conferred by this allele in psychosis, but will require replication.

Genetic examination of SETD7 and SUV39H1/H2 methyltransferases and the risk of diabetes complications in patients with type 1 diabetes

Hyperglycemia plays a pivotal role in the development and progression of vascular complications, which are the major sources of morbidity and mortality in diabetes. Furthermore, these vascular complications often persist and progress despite improved glucose control, possibly as a result of prior episodes of hyperglycemia. Epigenetic modifications mediated by histone methyltransferases are associated with gene-activating events that promote enhanced expression of key proinflammatory molecules implicated in vascular injury. In this study, we investigated genetic polymorphisms of the SETD7, SUV39H1, and SUV39H2 methyltransferases as predictors of risk for micro- and macrovascular complications in type 1 diabetes.

The HIF-1alpha C1772T polymorphism may be associated with susceptibility to clinically localised prostate cancer but not with elevated expression of hypoxic biomarkers

We investigated the role of the C1772T polymorphisms in exon 12 of the Hypoxia-inducible factor-1 alpha (HIF-1alpha) gene C1772T genotype in prostate cancer (PCa) and amplification of the hypoxic response. We identified the heterozygous germline CT genotype as an increased risk factor for clinically localised prostate cancer (Odds ratio = 6.2; p < 0.0001). While immunostaining intensity for HIF-1alpha and VEGF was significantly enhanced in 75% of PCa specimens when compared to matched benign specimens (p < 0.0001), the CT genotype did not modulate the kinetics of HIF-1alpha protein expression in hypoxia in vitro, and was not associated with enhanced expression of hypoxic biomarkers. This study provides the first evidence of an increased risk for clinically localised prostate cancer in men carrying the C1772T HIF-1alpha gene polymorphism. Although our results did not suggest an association between expression of hypoxic biomarkers and genotype status, the correlation may merit further investigation.

Variation in DNA repair genes XRCC3, XRCC4, XRCC5 and susceptibility to myeloma

Cytogenetic analysis in myeloma reveals marked chromosomal instability. Both widespread genomic alterations and evidence of aberrant class switch recombination, the physiological process that regulates maturation of the antibody response, implicate the DNA repair pathway in disease pathogenesis. We therefore assessed 27 SNPs in three genes (XRCC3, XRCC4 and XRCC5) central to DNA repair in patients with myeloma and controls from the EpiLymph study and from an Irish hospital registry (n = 306 cases, 263 controls). For the haplotype-tagging SNP (htSNP) rs963248 in XRCC4, Allele A was significantly more frequent in cases than in controls (86.4 versus 80.8%; odds ratio 1.51; 95% confidence interval 1.10-2.08; P = 0.0133), as was the AA genotype (74 versus 65%) (P = 0.026). Haplotype analysis was performed using Unphased for rs963248 in combination with additional SNPs in XRCC4. The strongest evidence of association came from the A-T haplotype from rs963248-rs2891980 (P = 0.008). For XRCC5, the genotype GG from rs1051685 was detected in 10 cases from different national populations but in only one control (P = 0.015). This SNP is located in the 3'-UTR of XRCC5. Overall, these data provide support for the hypothesis that common variation in the genes encoding DNA repair proteins contributes to susceptibility to myeloma.

Prevalence of the factor VR506Q mutation in two Irish control populations:use of a novel nested polymerase chain reaction approach

The prevalence of factor V (FV) Leiden among normal populations has primarily been determined using blood donors. This control group is carefully selected and therefore may not accurately reflect the true prevalence within the population. We assessed the prevalence of FV Leiden within the Irish population using Guthrie card samples randomly selected from all newborns. We compared this result with the prevalence of FV Leiden within blood donors. A novel nested polymerase chain reaction (PCR) method for FV Leiden was developed for analysis of the Guthrie card samples. There was no significant difference between the allele frequency within the Guthrie card samples and blood donors (2.07% vs. 2.35%, P = 0.66)

Prevalence of the prothrombin G20210A mutation in the Irish populations:use of a novel polymerase chain reaction approach

The prothrombin G20210A polymorphism is associated with a threefold-increased risk of venous thrombosis. There is considerable variation in the reported prevalence of this polymorphism within normal populations, ranging from 0 to 6.5%. The prevalence within the Irish population has not been determined. A restriction fragment length polymorphism (RFLP)-based assay is commonly used for the detection of the prothrombin 20210A allele. This assay does not include a control restriction digest fragment and, consequently, failure of the enzyme activity or lack of addition of enzyme to the sample cannot be distinguished from wild-type prothrombin. We developed a RFLP-based assay, which incorporates an invariant digest site, resulting in the generation of a control digest fragment. Furthermore, we developed a nested polymerase chain reaction (PCR) method for the amplification and digestion of poor-quality or low-concentration DNA. In the Irish population studied, five of 385 (1.29%) were heterozygous and one patient was homozygous for the prothrombin 20210A polymorphism. This is the first reported data on an Irish or Celtic population and suggests that the allele frequency is similar to Anglo-Saxon populations. The nested PCR method successfully amplified and digested 100/100 (100%) of the archived samples; none of these samples could be analyzed by the standard single-round PCR method. In conclusion, nested PCR should be considered in the analysis of archived samples. Single-round PCR is appropriate for recently collected samples; however, an invariant control digest site should be incorporated in RFLP-based assays to validate the integrity of the digestion enzyme and limit the risk of false-negative results.

Plastin Polymorphisms Predict Gender- and Stage-Specific Colon Cancer Recurrence after Adjuvant Chemotherapy

Tumor recurrence after curative resection remains a major problem in patients with locally advanced colorectal cancer treated with adjuvant chemotherapy. Genetic single-nucleotide polymorphisms (SNP) may serve as useful molecular markers to predict clinical outcomes in these patients and identify targets for future drug development. Recent in vitro and in vivo studies have demonstrated that the plastin genes PLS3 and LCP1 are overexpressed in colon cancer cells and play an important role in tumor cell invasion, adhesion, and migration. Hence, we hypothesized that functional genetic variations of plastin may have direct effects on the progression and prognosis of locally advanced colorectal cancer. We tested whether functional tagging polymorphisms of PLS3 and LCP1 predict time to tumor recurrence (TTR) in 732 patients (training set, 234; validation set, 498) with stage II/III colorectal cancer. The PLS3 rs11342 and LCP1 rs4941543 polymorphisms were associated with a significantly increased risk for recurrence in the training set. PLS3 rs6643869 showed a consistent association with TTR in the training and validation set, when stratified by gender and tumor location. Female patients with the PLS3 rs6643869 AA genotype had the shortest median TTR compared with those with any G allele in the training set [1.7 vs. 9.4 years; HR, 2.84; 95% confidence interval (CI), 1.32-6.1; P = 0.005] and validation set (3.3 vs. 13.7 years; HR, 2.07; 95% CI, 1.09-3.91; P = 0.021). Our findings suggest that several SNPs of the PLS3 and LCP1 genes could serve as gender- and/or stage-specific molecular predictors of tumor recurrence in stage II/III patients with colorectal cancer as well as potential therapeutic targets.

SNP discovery using Next Generation Transcriptomic Sequencing in Atlantic herring (Clupea harengus)

The introduction of Next Generation Sequencing (NGS) has revolutionised population genetics, providing studies of non-model species with unprecedented genomic coverage, allowing evolutionary biologists to address questions previously far beyond the reach of available resources. Furthermore, the simple mutation model of Single Nucleotide Polymorphisms (SNPs) permits cost-effective high-throughput genotyping in thousands of individuals simultaneously. Genomic resources are scarce for the Atlantic herring (Clupea harengus), a small pelagic species that sustains high revenue fisheries. This paper details the development of 578 SNPs using a combined NGS and high-throughput genotyping approach. Eight individuals covering the species distribution in the eastern Atlantic were bar-coded and multiplexed into a single cDNA library and sequenced using the 454 GS FLX platform. SNP discovery was performed by de novo sequence clustering and contig assembly, followed by the mapping of reads against consensus contig sequences. Selection of candidate SNPs for genotyping was conducted using an in silico approach. SNP validation and genotyping were performed simultaneously using an Illumina 1,536 GoldenGate assay. Although the conversion rate of candidate SNPs in the genotyping assay cannot be predicted in advance, this approach has the potential to maximise cost and time efficiencies by avoiding expensive and time-consuming laboratory stages of SNP validation. Additionally, the in silico approach leads to lower ascertainment bias in the resulting SNP panel as marker selection is based only on the ability to design primers and the predicted presence of intron-exon boundaries. Consequently SNPs with a wider spectrum of minor allele frequencies (MAFs) will be genotyped in the final panel. The genomic resources presented here represent a valuable multi-purpose resource for developing informative marker panels for population discrimination, microarray development and for population genomic studies in the wild.

Novel tools for conservation genomics: comparing two high-throughput approaches for SNP discovery in the transcriptome of the European hake

The growing accessibility to genomic resources using next-generation sequencing (NGS) technologies has revolutionized the application of molecular genetic tools to ecology and evolutionary studies in non-model organisms. Here we present the case study of the European hake (Merluccius merluccius), one of the most important demersal resources of European fisheries. Two sequencing platforms, the Roche 454 FLX (454) and the Illumina Genome Analyzer (GAII), were used for Single Nucleotide Polymorphisms (SNPs) discovery in the hake muscle transcriptome. De novo transcriptome assembly into unique contigs, annotation, and in silico SNP detection were carried out in parallel for 454 and GAII sequence data. High-throughput genotyping using the Illumina GoldenGate assay was performed for validating 1,536 putative SNPs. Validation results were analysed to compare the performances of 454 and GAII methods and to evaluate the role of several variables (e.g. sequencing depth, intron-exon structure, sequence quality and annotation). Despite well-known differences in sequence length and throughput, the two approaches showed similar assay conversion rates (approximately 43%) and percentages of polymorphic loci (67.5% and 63.3% for GAII and 454, respectively). Both NGS platforms therefore demonstrated to be suitable for large scale identification of SNPs in transcribed regions of non-model species, although the lack of a reference genome profoundly affects the genotyping success rate. The overall efficiency, however, can be improved using strict quality and filtering criteria for SNP selection (sequence quality, intron-exon structure, target region score).

β-Arrestin2 influences the response to methadone in opioid-dependent patients.

β-Arrestin2 (ARRB2) is a component of the G-protein-coupled receptor complex and is involved in μ-opioid and dopamine D(2) receptor signaling, two central processes in methadone signal transduction. We analyzed 238 patients in methadone maintenance treatment (MMT) and identified a haplotype block (rs34230287, rs3786047, rs1045280 and rs2036657) spanning almost the entire ARRB2 locus. Although none of these single nucleotide polymorphisms (SNPs) leads to a change in amino-acid sequence, we found that for all the SNPs analyzed, with exception of rs34230287, homozygosity for the variant allele confers a nonresponding phenotype (n=73; rs1045280C and rs2036657G: OR=3.1, 95% CI=1.5-6.3, P=0.004; rs3786047A: OR=2.5, 95% CI=1.2-5.1, P=0.02) also illustrated by a 12-fold shorter period of negative urine screening (P=0.01). The ARRB2 genotype may thus contribute to the interindividual variability in the response to MMT and help to predict response to treatment.