Thirty new loci for age at menarche identified by a meta-analysis of genome-wide association studies.

To identify loci for age at menarche, we performed a meta-analysis of 32 genome-wide association studies in 87,802 women of European descent, with replication in up to 14,731 women. In addition to the known loci at LIN28B (P = 5.4 × 10⁻⁶⁰) and 9q31.2 (P = 2.2 × 10⁻³³), we identified 30 new menarche loci (all P < 5 × 10⁻⁸) and found suggestive evidence for a further 10 loci (P < 1.9 × 10⁻⁶). The new loci included four previously associated with body mass index (in or near FTO, SEC16B, TRA2B and TMEM18), three in or near other genes implicated in energy homeostasis (BSX, CRTC1 and MCHR2) and three in or near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and gene-set enrichment pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing.

The population structure of Trypanosoma cruzi: expanded analysis of 54 strains using eight polymorphic CA-repeat microsatellites

Recently we cloned and sequenced the first eight Trypanosoma cruzi polymorphic microsatellite loci and studied 31 clones and strains to obtain valuable information about the population structure of the parasite. We have now studied 23 further strains, increasing from 11 to 31 the number of strains obtained from patients with chronic Chagas disease. This expanded set of 54 strains and clones analyzed with the eight microsatellites markers confirmed the previously observed diploidy, clonal population organization and very high polymorphism of T. cruzi. Moreover, this new study disclosed two new features of the population genetic structure of T. cruzi. The first was the discovery that, similarly to what we had previously shown for strains isolated from insect vectors, mammals and humans with acute disease, isolates from patients in the chronic phase of Chagas disease could also be multiclonal, albeit at a reduced proportion. Second, when we used parsimony to display the genetic relationship among the clonal lineages in an unrooted Wagner network we observed, like before, a good correlation of the tree topography with the classification in three clusters on the basis of single locus analysis of the ribosomal RNA genes. However, a significant new finding was that now the strains belonging to cluster 2 split in two distant sub-clusters. This observation suggests that the evolutionary history of T. cruzi may be more complex than we previously thought.

Analysis of genetic diversity of Trypanosoma cruzi: an application of riboprinting and gradient gel electrophoresis methods

Analysis of restriction fragment length polymorphism (RFLP) profiles derived from digestion of polymerase chain reaction (PCR) products of the ribosomal 18S from Trypanosoma cruzi yields a typical `riboprint' profile that can vary intraspecifically. A selection of 21 stocks of T. cruzi and three outgroup taxa: T. rangeli, T. conorhini and Leishmania braziliensis were analysed by riboprinting to assess divergence within and between taxa. T. rangeli, T. conorhini and L. braziliensis could be easily differentiated from each other and from T. cruzi. Phenetic analysis of PCR-RFLP profiles indicated that, with one or two exceptions, stocks of T. cruzi could be broadly partitioned into two groups that formally corresponded to T. cruzi I and T. cruzi II respectively. To test if ribosomal 18S sequences were homogeneous within each taxon, gradient gel electrophoresis methods were employed utilising either chemical or temperature gradients. Upon interpretation of the melting profiles of riboprints and a section of the 18S independently amplified by PCR, there would appear to be at least two divergent 18S types present within T. cruzi. Heterogeneity within copies of the ribosomal 18S within a single genome has therefore been demonstrated and interestingly, this dimorphic arrangement was also present in the outgroup taxa. Presumably the ancestral duplicative event that led to the divergent 18S types preceded that of speciation within this group. These divergent 18S paralogues may have, or had, different functional pressures or rates of molecular evolution. Whether or not these divergent types are equally transcriptionally active throughout the life cycle, remain to be assessed.

Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals.

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10(-8)-1.2×10(-43)). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10(-4)). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10(-3), n = 22,044), increased triglycerides (p = 2.6×10(-14), n = 93,440), increased waist-to-hip ratio (p = 1.8×10(-5), n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10(-3), n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10(-13), n = 96,748) and decreased BMI (p = 1.4×10(-4), n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.

Genome-wide association analysis identifies 20 loci that influence adult height.

Adult height is a model polygenic trait, but there has been limited success in identifying the genes underlying its normal variation. To identify genetic variants influencing adult human height, we used genome-wide association data from 13,665 individuals and genotyped 39 variants in an additional 16,482 samples. We identified 20 variants associated with adult height (P < 5 x 10(-7), with 10 reaching P < 1 x 10(-10)). Combined, the 20 SNPs explain approximately 3% of height variation, with a approximately 5 cm difference between the 6.2% of people with 17 or fewer 'tall' alleles compared to the 5.5% with 27 or more 'tall' alleles. The loci we identified implicate genes in Hedgehog signaling (IHH, HHIP, PTCH1), extracellular matrix (EFEMP1, ADAMTSL3, ACAN) and cancer (CDK6, HMGA2, DLEU7) pathways, and provide new insights into human growth and developmental processes. Finally, our results provide insights into the genetic architecture of a classic quantitative trait.

Identification of snails within the Bulinus africanus group from East Africa by multiplex SNaPshotTManalysis of single nucleotide polymorphisms within the cytochrome oxidase subunit I

Identification of populations of Bulinus nasutus and B. globosus from East Africa is unreliable using characters of the shell. In this paper, a molecular method of identification is presented for each species based on DNA sequence variation within the mitochondrial cytochrome oxidase subunit I (COI) as detected by a novel multiplexed SNaPshotTM assay. In total, snails from 7 localities from coastal Kenya were typed using this assay and variation within shell morphology was compared to reference material from Zanzibar. Four locations were found to contain B. nasutus and 2 locations were found to contain B. globosus. A mixed population containing both B. nasutus and B. globosus was found at Kinango. Morphometric variation between samples was considerable and UPGMA cluster analysis failed to differentiate species. The multiplex SNaPshotTM assay is an important development for more precise methods of identification of B. africanus group snails. The assay could be further broadened for identification of other snail intermediate host species.

Genetic and epigenetic analysis of SSAT gene dysregulation in suicidal behavior.

It has recently been proposed that the SSAT gene plays a role in the predisposition to suicidal behavior. SSAT expression was found to be down-regulated in the brain of suicide completers. In addition, a single nucleotide polymorphism (SNP) rs6526342 was associated both with variation in SSAT expression and with suicidal behavior. In this study, we aimed to characterize the relationship between SSAT dysregulation and suicide behavior. To this end, we measured SSAT expression levels in the ventral prefrontal cortex (VPFC) of suicide completers (n = 20) and controls (n = 20) and found them to be significantly down-regulated in suicide victims (P = 0.007). To identify the basis of the regulation of SSAT expression, we performed an association analysis of 309 SNPs with SSAT transcript levels in 53 lymphoblastoid cell lines from the CEPH collection. We then examined the methylation status of the SSAT promoter region in males and females suicide completers and control subjects whose SSAT brain expression had been measured. We found no evidence to support a role for SNPs in controlling the level of SSAT expression. SSAT promoter methylation levels were not different between suicide completers and controls and did not correlate with SSAT expression levels. In addition, we found no indication of a genetic association between suicidal behavior and SNPs located within the SSAT gene. Our study provides new results which show that dysregulation of SSAT expression does play a role in suicide behavior. However, our data do not support any association between rs6526342 and variation in SSAT expression or suicidal behavior.

DIP-STR: Highly Sensitive Markers for the Analysis of Unbalanced Genomic Mixtures.

Samples containing highly unbalanced DNA mixtures from two individuals commonly occur both in forensic mixed stains and in peripheral blood DNA microchimerism induced by pregnancy or following organ transplant. Because of PCR amplification bias, the genetic identification of a DNA that contributes trace amounts to a mixed sample represents a tremendous challenge. This means that standard genetic markers, namely microsatellites, also referred as short tandem repeats (STR), and single-nucleotide polymorphism (SNP) have limited power in addressing common questions of forensic and medical genetics. To address this issue, we developed a molecular marker, named DIP-STR that relies on pairing deletion-insertion polymorphisms (DIP) with STR. This novel analytical approach allows for the unambiguous genotyping of a minor component in the presence of a major component, where DIP-STR genotypes of the minor were successfully procured at ratios up to 1:1,000. The compound nature of this marker generates a high level of polymorphism that is suitable for identity testing. Here, we demonstrate the power of the DIP-STR approach on an initial set of nine markers surveyed in a Swiss population. Finally, we discuss the limitations and potential applications of our new system including preliminary tests on clinical samples and estimates of their performance on simulated DNA mixtures.

Large-scale gene-centric analysis identifies novel variants for coronary artery disease.

Coronary artery disease (CAD) has a significant genetic contribution that is incompletely characterized. To complement genome-wide association (GWA) studies, we conducted a large and systematic candidate gene study of CAD susceptibility, including analysis of many uncommon and functional variants. We examined 49,094 genetic variants in ∼2,100 genes of cardiovascular relevance, using a customised gene array in 15,596 CAD cases and 34,992 controls (11,202 cases and 30,733 controls of European descent; 4,394 cases and 4,259 controls of South Asian origin). We attempted to replicate putative novel associations in an additional 17,121 CAD cases and 40,473 controls. Potential mechanisms through which the novel variants could affect CAD risk were explored through association tests with vascular risk factors and gene expression. We confirmed associations of several previously known CAD susceptibility loci (eg, 9p21.3:p<10(-33); LPA:p<10(-19); 1p13.3:p<10(-17)) as well as three recently discovered loci (COL4A1/COL4A2, ZC3HC1, CYP17A1:p<5×10(-7)). However, we found essentially null results for most previously suggested CAD candidate genes. In our replication study of 24 promising common variants, we identified novel associations of variants in or near LIPA, IL5, TRIB1, and ABCG5/ABCG8, with per-allele odds ratios for CAD risk with each of the novel variants ranging from 1.06-1.09. Associations with variants at LIPA, TRIB1, and ABCG5/ABCG8 were supported by gene expression data or effects on lipid levels. Apart from the previously reported variants in LPA, none of the other ∼4,500 low frequency and functional variants showed a strong effect. Associations in South Asians did not differ appreciably from those in Europeans, except for 9p21.3 (per-allele odds ratio: 1.14 versus 1.27 respectively; P for heterogeneity = 0.003). This large-scale gene-centric analysis has identified several novel genes for CAD that relate to diverse biochemical and cellular functions and clarified the literature with regard to many previously suggested genes.

Comprehensive analysis of RET common and rare variant patientsin a series of Spanish Hirschsprung patients confirms a synergistic effect of both kinds of events

Background. RET is the major gene associated to Hirschsprung disease (HSCR) with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. In the present study, we have performed a comprehensive study of our HSCR series evaluating the involvement of both RET rare variants (RVs) and common variants (CVs) in the context of the disease. Methods. RET mutational screening was performed by dHPLC and direct sequencing for the identification of RVs. In addition Taqman technology was applied for the genotyping of 3 RET CVs previously associated to HSCR, including a variant lying in an enhancer domain within RET intron 1 (rs2435357). Statistical analyses were performed using the SPSS v.17.0 to analyze the distribution of the variants. Results. Our results confirm the strongest association to HSCR for the "enhancer" variant, and demonstrate a significantly higher impact of it in male versus female patients. Integration of the RET RVs and CVs analysis showed that in 91.66% of cases with both kinds of mutational events, the enhancer allele is in trans with the allele bearing the RET RV. Conclusions. A gender effect exists on both the transmission and distribution of rare coding and common HSCR causing mutations. In addition, these RET CVs and RVs seem to act in a synergistic way leading to HSCR phenotype.

Association analysis of urotensin II gene (UTS2) and flanking regions with biochemical parameters related to insulin resistance.

Background. Urotensin II (UII) is a potent vasoconstrictor peptide, which signals through a G-protein coupled receptor (GPCR) known as GPR14 or urotensin receptor (UTR). UII exerts a broad spectrum of actions in several systems such as vascular cell, heart muscle or pancreas, where it inhibits insulin release. Objective. Given the reported role of UII in insulin secretion, we have performed a genetic association analysis of the UTS2 gene and flanking regions with biochemical parameters related to insulin resistance (fasting glucose, glucose 2 hours after a glucose overload, fasting insulin and insulin resistance estimated as HOMA). Results and Conclusions. We have identified several polymorphisms associated with the analysed clinical traits, not only at the UTS2 gene, but also in thePER3 gene, located upstream from UTS2. Our results are compatible with a role for UII in glucose homeostasis and diabetes although we cannot rule out the possibility that PER3 gene may underlie the reported associations.

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.

The T-381C SNP in BNP gene may be modestly associated with type 2 diabetes: an updated meta-analysis in 49 279 subjects.

A recent study reported an association between the brain natriuretic peptide (BNP) promoter T-381C polymorphism (rs198389) and protection against type 2 diabetes (T2D). As replication in several studies is mandatory to confirm genetic results, we analyzed the T-381C polymorphism in seven independent case-control cohorts and in 291 T2D-enriched pedigrees totalling 39 557 subjects of European origin. A meta-analysis of the seven case-control studies (n = 39 040) showed a nominal protective effect [odds ratio (OR) = 0.86 (0.79-0.94), P = 0.0006] of the CC genotype on T2D risk, consistent with the previous study. By combining all available data (n = 49 279), we further confirmed a modest contribution of the BNP T-381C polymorphism for protection against T2D [OR = 0.86 (0.80-0.92), P = 1.4 x 10(-5)]. Potential confounders such as gender, age, obesity status or family history were tested in 4335 T2D and 4179 normoglycemic subjects and they had no influence on T2D risk. This study provides further evidence of a modest contribution of the BNP T-381C polymorphism in protection against T2D and illustrates the difficulty of unambiguously proving modest-sized associations even with large sample sizes.

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.

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.