Fractal and complex network analyses of protein molecular dynamics

Based on protein molecular dynamics, we investigate the fractal properties of energy, pressure and volume time series using the multifractal detrended fluctuation analysis (MF-DFA) and the topological and fractal properties of their converted horizontal visibility graphs (HVGs). The energy parameters of protein dynamics we considered are bonded potential, angle potential, dihedral potential, improper potential, kinetic energy, Van der Waals potential, electrostatic potential, total energy and potential energy. The shape of the h(q)h(q) curves from MF-DFA indicates that these time series are multifractal. The numerical values of the exponent h(2)h(2) of MF-DFA show that the series of total energy and potential energy are non-stationary and anti-persistent; the other time series are stationary and persistent apart from series of pressure (with H≈0.5H≈0.5 indicating the absence of long-range correlation). The degree distributions of their converted HVGs show that these networks are exponential. The results of fractal analysis show that fractality exists in these converted HVGs. For each energy, pressure or volume parameter, it is found that the values of h(2)h(2) of MF-DFA on the time series, exponent λλ of the exponential degree distribution and fractal dimension dBdB of their converted HVGs do not change much for different proteins (indicating some universality). We also found that after taking average over all proteins, there is a linear relationship between 〈h(2)〉〈h(2)〉 (from MF-DFA on time series) and 〈dB〉〈dB〉 of the converted HVGs for different energy, pressure and volume.

Nutritional status and the foraging behaviour of Bactrocera tryoni with particular reference to protein bait spray

Poisoned protein baits comprise a recognized method for controlling tephritid fruit flies in the form of a ‘lure-and-kill’ technique. However, little is known about how a fly's internal protein and carbohydrate levels (i.e. nutritional status) might influence the efficacy of this control. In the present study, the relationships between the internal levels of protein (as measured by total body nitrogen) and carbohydrate (as measured by total body carbon) of the fruit fly Bactrocera tryoni (Froggatt) (Diptera: Tephritidae) are investigated, as well as its foraging behaviours in response to protein, fruit and cue-lure (a male-specific attractant) baits. Small cage behavioural experiments are conducted using flies from cultures of different nutritional status and wild flies sampled from the field during the fruiting cycle of a guava crop. For female flies, increasing total body nitrogen is correlated with decreased protein foraging and increased oviposition activity; increasing total body carbon levels generate the same behavioural changes except that the oviposition response is not significant. For males, there are no significant correlations between changes in total body nitrogen and total body carbon and protein or cue-lure foraging. For wild flies from the guava orchard, almost all of them are sexually mature when entering the crop and, over the entire season, total body nitrogen and total body carbon levels are such that protein hunger is unlikely for most flies. The results infer strongly that the requirements of wild, sexually mature flies for protein are minimal and that flies can readily gain sufficient nutrients from wild sources for their physiological needs. The results offer a mechanistic explanation for the poor response of male and mature female fruit flies to protein bait spray.

Secondary structure element alignment kernel method for prediction of protein structural classes

In this paper, we aim at predicting protein structural classes for low-homology data sets based on predicted secondary structures. We propose a new and simple kernel method, named as SSEAKSVM, to predict protein structural classes. The secondary structures of all protein sequences are obtained by using the tool PSIPRED and then a linear kernel on the basis of secondary structure element alignment scores is constructed for training a support vector machine classifier without parameter adjusting. Our method SSEAKSVM was evaluated on two low-homology datasets 25PDB and 1189 with sequence homology being 25% and 40%, respectively. The jackknife test is used to test and compare our method with other existing methods. The overall accuracies on these two data sets are 86.3% and 84.5%, respectively, which are higher than those obtained by other existing methods. Especially, our method achieves higher accuracies (88.1% and 88.5%) for differentiating the α + β class and the α/β class compared to other methods. This suggests that our method is valuable to predict protein structural classes particularly for low-homology protein sequences. The source code of the method in this paper can be downloaded at http://math.xtu.edu.cn/myphp/math/research/source/SSEAK_source_code.rar.

Genetic determinants of heel bone properties: Genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n 5 14 260), velocity of sound (VOS; n 5 15 514) and BMD (n 5 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n 5 11 452) and new genotyping in 15 cohorts (de novo n 5 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 3 108) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 3 1014). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 3 106 also had the expected direction of association with any fracture (P < 0.05), including threeSNPswithP < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, thisGWAstudy reveals the effect of several genescommon to central DXA-derivedBMDand heel ultrasound/DXAmeasures and points to anewgenetic locus with potential implications for better understanding of osteoporosis pathophysiology.

The function of the RNA-binding protein TEL1 in moss reveals ancient regulatory mechanisms of shoot development

The shoot represents the basic body plan in land plants. It consists of a repeated structure composed of stems and leaves. Whereas vascular plants generate a shoot in their diploid phase, non-vascular plants such as mosses form a shoot (called the gametophore) in their haploid generation. The evolution of regulatory mechanisms or genetic networks used in the development of these two kinds of shoots is unclear. TERMINAL EAR1-like genes have been involved in diploid shoot development in vascular plants. Here, we show that disruption of PpTEL1 from the moss Physcomitrella patens, causes reduced protonema growth and gametophore initiation, as well as defects in gametophore development. Leafy shoots formed on ΔTEL1 mutants exhibit shorter stems with more leaves per shoot, suggesting an accelerated leaf initiation (shortened plastochron), a phenotype shared with the Poaceae vascular plants TE1 and PLA2/LHD2 mutants. Moreover, the positive correlation between plastochron length and leaf size observed in ΔTEL1 mutants suggests a conserved compensatory mechanism correlating leaf growth and leaf initiation rate that would minimize overall changes in plant biomass. The RNA-binding protein encoded by PpTEL1 contains two N-terminus RNA-recognition motifs, and a third C-terminus non-canonical RRM, specific to TEL proteins. Removal of the PpTEL1 C-terminus (including this third RRM) or only 16–18 amino acids within it seriously impairs PpTEL1 function, suggesting a critical role for this third RRM. These results show a conserved function of the RNA-binding PpTEL1 protein in the regulation of shoot development, from early ancestors to vascular plants, that depends on the third TEL-specific RRM.

Meta-analysis of genome-wide studies identifies MEF2C SNPs associated with bone mineral density at forearm

Background Forearm fractures affect 1.7 million individuals worldwide each year and most occur earlier in life than hip fractures. While the heritability of forearm bone mineral density (BMD) and fracture is high, their genetic determinants are largely unknown. Aim To identify genetic variants associated with forearm BMD and forearm fractures. Methods BMD at distal radius, measured by dualenergy x-ray absorptiometry, was tested for association with common genetic variants. We conducted a metaanalysis of genome-wide association studies for BMD in 5866 subjects of European descent and then selected the variants for replication in 715 Mexican American samples. Gene-based association was carried out to supplement the single-nucleotide polymorphism (SNP) association test. We then tested the BMD-associated SNPs for association with forearm fracture in 2023 cases and 3740 controls. Results We found that five SNPs in the introns of MEF2C were associated with forearm BMD at a genome-wide significance level (p<5×10-8) in meta-analysis (lead SNP, rs11951031[T] -0.20 SDs per allele, p=9.01×10-9). The gene-based association test suggested an association between MEF2C and forearm BMD ( p=0.003). The association between MEF2C variants and risk of fracture did not achieve statistical significance (SNP rs12521522[A]: OR=1.14 (95% CI 0.92 to 1.35), p=0.14). Meta-analysis also revealed two genome-wide suggestive loci at CTNNA2 and 6q23.2. Conclusions These findings demonstrate that variants at MEF2C were associated with forearm BMD, implicating this gene in the determination of BMD at forearm.

Genome-wide association study of intraocular pressure identifies the GLCCI1/ICA1 region as a glaucoma susceptibility locus

To discover quantitative trait loci for intraocular pressure, amajor risk factor for glaucoma and the only modifiable one,weperformed agenome-wide association studyonadiscoverycohort of2175individualsfromSydney, Australia. We found a novel association between intraocular pressure and a common variant at 7p21 near to GLCCI1 and ICA1. The findings in this region were confirmed through two UK replication cohorts totalling 4866 individuals (rs59072263, Pcombined = 1.10 × 10-8). A copy of the G allele at this SNP is associated with an increase in mean IOP of 0.45 mmHg (95%CI = 0.30-0.61 mmHg). These results lend support to the implication of vesicle trafficking and glucocorticoid inducibility pathways in the determination of intraocular pressure and in the pathogenesis of primary open-angle glaucoma.

Mosaic PPM1D mutations are associated with predisposition to breast and ovarian cancer

Improved sequencing technologies offer unprecedented opportunities for investigating the role of rare genetic variation in common disease. However, there are considerable challenges with respect to study design, data analysis and replication. Using pooled next-generation sequencing of 507 genes implicated in the repair of DNA in 1,150 samples, an analytical strategy focused on protein-truncating variants (PTVs) and a large-scale sequencing case-control replication experiment in 13,642 individuals, here we show that rare PTVs in the p53-inducible protein phosphatase PPM1D are associated with predisposition to breast cancer and ovarian cancer. PPM1D PTV mutations were present in 25 out of 7,781 cases versus 1 out of 5,861 controls (P = 1.12 × 10-5), including 18 mutations in 6,912 individuals with breast cancer (P = 2.42 × 10-4) and 12 mutations in 1,121 individuals with ovarian cancer (P = 3.10 × 10-9). Notably, all of the identified PPM1D PTVs were mosaic in lymphocyte DNA and clustered within a 370-base-pair region in the final exon of the gene, carboxy-terminal to the phosphatase catalytic domain. Functional studies demonstrate that the mutations result in enhanced suppression of p53 in response to ionizing radiation exposure, suggesting that the mutant alleles encode hyperactive PPM1D isoforms. Thus, although the mutations cause premature protein truncation, they do not result in the simple loss-of-function effect typically associated with this class of variant, but instead probably have a gain-of-function effect. Our results have implications for the detection and management of breast and ovarian cancer risk. More generally, these data provide new insights into the role of rare and of mosaic genetic variants in common conditions, and the use of sequencing in their identification.

Genome-wide association analysis identifies 13 new risk loci for schizophrenia

Schizophrenia is an idiopathic mental disorder with a heritable component and a substantial public health impact. We conducted a multi-stage genome-wide association study (GWAS) for schizophrenia beginning with a Swedish national sample (5,001 cases and 6,243 controls) followed by meta-Analysis with previous schizophrenia GWAS (8,832 cases and 12,067 controls) and finally by replication of SNPs in 168 genomic regions in independent samples (7,413 cases, 19,762 controls and 581 parent-offspring trios). We identified 22 loci associated at genome-wide significance; 13 of these are new, and 1 was previously implicated in bipolar disorder. Examination of candidate genes at these loci suggests the involvement of neuronal calcium signaling. We estimate that 8,300 independent, mostly common SNPs (95% credible interval of 6,300-10,200 SNPs) contribute to risk for schizophrenia and that these collectively account for at least 32% of the variance in liability. Common genetic variation has an important role in the etiology of schizophrenia, and larger studies will allow more detailed understanding of this disorder.

High-resolution SNP microarray investigation of copy number variations on chromosome 18 in a control cohort

Copy number variations (CNVs) as described in the healthy population are purported to contribute significantly to genetic heterogeneity. Recent studies have described CNVs using lymphoblastoid cell lines or by application of specifically developed algorithms to interrogate previously described data. However, the full extent of CNVs remains unclear. Using high-density SNP array, we have undertaken a comprehensive investigation of chromosome 18 for CNV discovery and characterisation of distribution and association with chromosome architecture. We identified 399 CNVs, of which loss represents 98%, 58% are less than 2.5 kb in size and 71% are intergenic. Intronic deletions account for the majority of copy number changes with gene involvement. Furthermore, one-third of CNVs do not have putative breakpoints within repetitive sequences. We conclude that replicative processes, mediated either by repetitive elements or microhomology, account for the majority of CNVs in the healthy population. Genomic instability involving the formation of a non-B structure is demonstrated in one region.

Common variants at the MHC locus and at chromosome 16q24.1 predispose to Barrett's Esophagus

Barrett's esophagus is an increasingly common disease that is strongly associated with reflux of stomach acid and usually a hiatus hernia, and it strongly predisposes to esophageal adenocarcinoma (EAC), a tumor with a very poor prognosis. We report the first genome-wide association study on Barrett's esophagus, comprising 1,852 UK cases and 5,172 UK controls in the discovery stage and 5,986 cases and 12,825 controls in the replication stage. Variants at two loci were associated with disease risk: chromosome 6p21, rs9257809 (P combined = 4.09 × 10-9; odds ratio (OR) = 1.21, 95% confidence interval (CI) =1.13-1.28), within the major histocompatibility complex locus, and chromosome 16q24, rs9936833 (P combined = 2.74 × 10-10; OR = 1.14, 95% CI = 1.10-1.19), for which the closest protein-coding gene is FOXF1, which is implicated in esophageal development and structure. We found evidence that many common variants of small effect contribute to genetic susceptibility to Barrett's esophagus and that SNP alleles predisposing to obesity also increase risk for Barrett's esophagus. © 2012 Nature America, Inc. All rights reserved.

Fine mapping and replication of genetic risk Loci in primary sclerosing cholangitis

Background and aims. Primary sclerosing cholangitis (PSC) is a chronic cholestatic liver disease characterized by progressive inflammation and fibrosis of the bile ducts eventually leading to biliary cirrhosis. Recent genetic studies in PSC have identified associations at 2q13, 2q35, 3p21, 4q27, 13q31 and suggestive association at 10p15. The aim of this study was to further characterize and refine the genetic architecture of PSC. Methods. We analyzed previously reported associated SNPs at four of these non-HLA loci and 59 SNPs tagging the IL-2/IL-21 (4q27) and IL2RA (10p15) loci in 992 UK PSC cases and 5162 healthy UK controls. Results. The most associated SNPs identified were rs3197999 (3p21 (MST1), p = 1.9 × 10 -6, OR A vs G = 1.28, 95% CI (1.16-1.42)); rs4147359 (10p15 (IL2RA), p = 2.6 × 10 -4, OR A vs G = 1.20, 95% CI (1.09-1.33)) and rs12511287 (4q27 (IL-2/IL-21), p = 3.0 × 10 -4, OR A vs T = 1.21, 95% CI (1.09-1.35)). In addition, we performed a meta-analysis for selected SNPs using published summary statistics from recent studies. We observed genome-wide significance for rs3197999 (3p21 (MST1), P combined = 3.8 × 10 -12) and rs4147359 (10p15 (IL2RA), P combined = 1.5 × 10 -8). Conclusion. We have for the first time confirmed the association of PSC with genetic variants at 10p15 (IL2RA) locus at genome-wide significance and replicated the associations at MST1 and IL-2/IL-21 loci in a large homogeneous UK population. These results strongly implicate the role of IL-2/IL2RA pathway in PSC and provide further confirmation of MST1 association. © Informa Healthcare.

Common variants near ATM are associated with glycemic response to metformin in type 2 diabetes

Metformin is the most commonly used pharmacological therapy for type 2 diabetes. We report a genome-wide association study for glycemic response to metformin in 1,024 Scottish individuals with type 2 diabetes with replication in two cohorts including 1,783 Scottish individuals and 1,113 individuals from the UK Prospective Diabetes Study. In a combined meta-analysis, we identified a SNP, rs11212617, associated with treatment success (n = 3,920, P = 2.9 P×-9, odds ratio = 1.35, 95% CI 1.22-1.49) at a locus containing ATM, the ataxia telangiectasia mutated gene. In a rat hepatoma cell line, inhibition of ATM with KU-55933 attenuated the phosphorylation and activation of AMP-activated protein kinase in response to metformin. We conclude that ATM, a gene known to be involved in DNA repair and cell cycle control, plays a role in the effect of metformin upstream of AMP-activated protein kinase, and variation in this gene alters glycemic response to metformin. © 2011 Nature America, Inc. All rights reserved.

Genome-wide association study identifies susceptibility loci for open angle glaucoma at TMCO1 and CDKN2B-AS1

We report a genome-wide association study for open-angle glaucoma (OAG) blindness using a discovery cohort of 590 individuals with severe visual field loss (cases) and 3,956 controls. We identified associated loci at TMCO1 (rs4656461[G] odds ratio (OR) = 1.68, P = 6.1 × 10-10) and CDKN2B-AS1 (rs4977756[A] OR = 1.50, P = 4.7 × 10-9). We replicated these associations in an independent cohort of cases with advanced OAG (rs4656461 P = 0.010; rs4977756 P = 0.042) and two additional cohorts of less severe OAG (rs4656461 combined discovery and replication P = 6.00 × 10-14, OR = 1.51, 95% CI 1.35-1.68; rs4977756 combined P = 1.35 × 10-14, OR = 1.39, 95% CI 1.28-1.51). We show retinal expression of genes at both loci in human ocular tissues. We also show that CDKN2A and CDKN2B are upregulated in the retina of a rat model of glaucoma. © 2011 Nature America, Inc. All rights reserved.

Elucidating the chromosome 9 association with AS; CARD9 is a candidate gene

Ankylosing spondylitis (AS) is polygenic with contributions from the immunologically relevant genes HLA-B27, ERAP1 and IL23R. A recent genome-wide association screen (GWAS) identified associations (P0.005) with the non-synonymous single-nucleotide polymorphisms (nsSNPs), rs4077515 and rs3812571, in caspase recruitment domain-containing protein 9 (CARD9) and small nuclear RNA-activating complex polypeptide 4 (SNAPC4) on chromosome 9q that had previously been linked to AS. We replicated these associations in a study of 730 AS patients compared with 2879 historic disease controls (rs4077515 P0.0004, odds ratio (OR)1.2, 95% confidence interval (CI)1.1-1.4; rs3812571 P0.0003, OR1.2, 95% CI1.1-1.4). Meta-analysis revealed strong associations of both SNPs with AS, rs4077515 P0.000005, OR1.2, 95% CI1.1-1.3 and rs3812571 P0.000006, OR1.2, 95% CI1.1-1.3. We then typed 1604 AS cases and 1020 controls for 13 tagging SNPs; 6 showed at least nominal association, 5 of which were in CARD9. We imputed genotypes for 13 additional SNPs but none was more strongly associated with AS than the tagging SNPs. Finally, interrogation of an mRNA expression database revealed that the SNPs most strongly associated with AS (or in strong linkage disequilibrium) were those most associated with CARD9 expression. CARD9 is a plausible candidate for AS given its central role in the innate immune response.