Multivariate Statistical Analysis Applied to an IL6 Signal Transduction Model in Hepatocytes

This paper introduces the application of linear multivariate statistical techniques, including partial least squares (PLS), canonical correlation analysis (CCA) and reduced rank regression (RRR), into the area of Systems Biology. This new approach aims to extract the important proteins embedded in complex signal transduction pathway models.The analysis is performed on a model of intracellular signalling along the janus-associated kinases/signal transducers and transcription factors (JAK/STAT) and mitogen activated protein kinases (MAPK) signal transduction pathways in interleukin-6 (IL6) stimulated hepatocytes, which produce signal transducer and activator of transcription factor 3 (STAT3).A region of redundancy within the MAPK pathway that does not affect the STAT3 transcription was identified using CCA. This is the core finding of this analysis and cannot be obtained by inspecting the model by eye. In addition, RRR was found to isolate terms that do not significantly contribute to changes in protein concentrations, while the application of PLS does not provide such a detailed picture by virtue of its construction.This analysis has a similar objective to conventional model reduction techniques with the advantage of maintaining the meaning of the states prior to and after the reduction process. A significant model reduction is performed, with a marginal loss in accuracy, offering a more concise model while maintaining the main influencing factors on the STAT3 transcription.The findings offer a deeper understanding of the reaction terms involved, confirm the relevance of several proteins to the production of Acute Phase Proteins and complement existing findings regarding cross-talk between the two signalling pathways.

Genome-wide linkage analysis for human longevity: Genetics of healthy aging study

Clear evidence exists for heritability of humanlongevity, and much interest is focused on identifying genes associated with longer lives. To identify such longevity alleles, we performed the largest genome-wide linkage scan thus far reported. Linkage analyses included 2118nonagenarian Caucasian sibling pairs that have been enrolled in 15 study centers of 11 European countries as part of the Genetics of Healthy Aging (GEHA) project. In the joint linkage analyses, we observed four regions that show linkage with longevity; chromosome 14q11.2 (LOD = 3.47), chromosome 17q12-q22 (LOD = 2.95), chromosome 19p13.3-p13.11 (LOD = 3.76), and chromosome 19q13.11-q13.32 (LOD = 3.57). To fine map these regions linked to longevity, we performed association analysis using GWAS data in a subgroup of 1228 unrelated nonagenarian and 1907 geographically matched controls. Using a fixed-effect meta-analysis approach, rs4420638 at the TOMM40/ APOE/APOC1 gene locus showed significant association with longevity (P-value = 9.6 × 10). By combined modeling of linkage and association, we showed that association of longevity with APOEe4 and APOEe2 alleles explain the linkage at 19q13.11-q13.32 with P-value = 0.02 and P-value = 1.0 × 10, respectively. In the largest linkage scan thus far performed for human familial longevity, we confirm that the APOE locus is a longevity gene and that additional longevity loci may be identified at 14q11.2, 17q12-q22, and 19p13.3-p13.11. As the latter linkage results are not explained by common variants, we suggest that rare variants play an important role in human familial longevity.

Meta-analysis of 32 genome-wide linkage studies of schizophrenia

A genome scan meta-analysis (GSMA) was carried out on 32 independent genome-wide linkage scan analyses that included 3255 pedigrees with 7413 genotyped cases affected with schizophrenia (SCZ) or related disorders. The primary GSMA divided the autosomes into 120 bins, rank-ordered the bins within each study according to the most positive linkage result in each bin, summed these ranks (weighted for study size) for each bin across studies and determined the empirical probability of a given summed rank (P-SR) by simulation. Suggestive evidence for linkage was observed in two single bins, on chromosomes 5q (142-168 Mb) and 2q (103-134 Mb). Genome-wide evidence for linkage was detected on chromosome 2q (119-152 Mb) when bin boundaries were shifted to the middle of the previous bins. The primary analysis met empirical criteria for 'aggregate' genome-wide significance, indicating that some or all of 10 bins are likely to contain loci linked to SCZ, including regions of chromosomes 1, 2q, 3q, 4q, 5q, 8p and 10q. In a secondary analysis of 22 studies of European-ancestry samples, suggestive evidence for linkage was observed on chromosome 8p (16-33 Mb). Although the newer genome-wide association methodology has greater power to detect weak associations to single common DNA sequence variants, linkage analysis can detect diverse genetic effects that segregate in families, including multiple rare variants within one locus or several weakly associated loci in the same region. Therefore, the regions supported by this meta-analysis deserve close attention in future studies. Molecular Psychiatry (2009) 14, 774-785; doi:10.1038/mp.2008.135; published online 30 December 2008

A genome-wide linkage scan for genes controlling variation in urinary albumin excretion in type ll diabetes

The main hallmark of diabetic nephropathy is elevation in urinary albumin excretion. We performed a genome-wide linkage scan in 63 extended families with multiple members with type II diabetes. Urinary albumin excretion, measured as the albumin-to-creatinine ratio (ACR), was determined in 426 diabetic and 431 nondiabetic relatives who were genotyped for 383 markers. The data were analyzed using variance components linkage analysis. Heritability (h2) of ACR was significant in diabetic (h2=0.23, P=0.0007), and nondiabetic (h2=0.39, P=0.0001) relatives. There was no significant difference in genetic variance of ACR between diabetic and nondiabetic relatives (P=0.16), and the genetic correlation (rG=0.64) for ACR between these two groups was not different from 1 (P=0.12). These results suggested that similar genes contribute to variation in ACR in diabetic and nondiabetic relatives. This hypothesis was supported further by the linkage results.

Genomewide linkage scan of schizophrenia in a large multicenter pedigree sample using single nucleotide polymorphisms

A genomewide linkage scan was carried out in eight clinical samples of informative schizophrenia families. After all quality control checks, the analysis of 707 European-ancestry families included 1615 affected and 1602 unaffected genotyped individuals, and the analysis of all 807 families included 1900 affected and 1839 unaffected individuals. Multipoint linkage analysis with correction for marker-marker linkage disequilibrium was carried out with 5861 single nucleotide polymorphisms (SNPs; Illumina version 4.0 linkage map). Suggestive evidence for linkage ( European families) was observed on chromosomes 8p21, 8q24.1, 9q34 and 12q24.1 in nonparametric and/or parametric analyses. In a logistic regression allele-sharing analysis of linkage allowing for intersite heterogeneity, genomewide significant evidence for linkage was observed on chromosome 10p12. Significant heterogeneity was also observed on chromosome 22q11.1. Evidence for linkage across family sets and analyses was most consistent on chromosome 8p21, with a one-LOD support interval that does not include the candidate gene NRG1, suggesting that one or more other susceptibility loci might exist in the region. In this era of genomewide association and deep resequencing studies, consensus linkage regions deserve continued attention, given that linkage signals can be produced by many types of genomic variation, including any combination of multiple common or rare SNPs or copy number variants in a region. Molecular Psychiatry (2009) 14, 786-795; doi:10.1038/mp.2009.11; published online 17 February 2009

Irish study on high-density schizophrenia families:field methods and power to detect linkage

Large samples of multiplex pedigrees will probably be needed to detect susceptibility loci for schizophrenia by linkage analysis. Standardized ascertainment of such pedigrees from culturally and ethnically homogeneous populations may improve the probability of detection and replication of linkage. The Irish Study of High-Density Schizophrenia Families (ISHDSF) was formed from standardized ascertainment of multiplex schizophrenia families in 39 psychiatric facilities covering over 90% of the population in Ireland and Northern Ireland. We here describe a phenotypic sample and a subset thereof, the linkage sample. Individuals were included in the phenotypic sample if adequate diagnostic information, based on personal interview and/or hospital record, was available. Only individuals with available DNA were included in the linkage sample. Inclusion of a pedigree into the phenotypic sample required at least two first, second, or third degree relatives with non-affective psychosis (NAP), one whom had schizophrenia (S) or poor-outcome schizo-affective disorder (PO-SAD). Entry into the linkage sample required DNA samples on at least two individuals with NAP, of whom at least one had S or PO-SAD. Affection was defined by narrow, intermediate, and broad criteria. The phenotypic sample contained 277 pedigrees and 1,770 individuals and the linkage sample 265 pedigrees and 1,408 individuals. Using the intermediate definition of affection, the phenotypic sample contained 837 affected individuals and 526 affected sibling pairs. Parallel figures for the linkage sample were 700 and 420. Individuals with schizophrenia from these multiplex pedigrees resembled epidemiologically sampled cases with respect to age at onset, gender distribution, and most clinical symptoms, although they were more thought-disordered and had a poorer outcome. Power analyses based on the model of linkage heterogeneity indicated that the ISHDSF should be able to detect a major locus that influences susceptibility to schizophrenia in as few as 20% of families. Compared to first-degree relatives of epidemiologically sampled schizophrenic probands, first-degree relatives of schizophrenic members from the ISHDSF had a similar risk for schizotypal personality disorder, affective illness, alcoholism, and anxiety disorder. With sufficient resources, large-scale ascertainment of multiplex schizophrenia pedigrees is feasible, especially in countries with catchmented psychiatric care and stable populations. Although somewhat more severely ill, schizophrenic members of such pedigrees appear to clinically resemble typical schizophrenic patients. Our ascertainment process for multiplex schizophrenia families did not select for excess familial risk for affective illness or alcoholism. With its large sample ascertained in a standardized manner from a relatively homogeneous population, the ISHDSF provides considerable power to detect susceptibility loci for schizophrenia.

No evidence for linkage or association of neuregulin-1 (NRG1) with disease in the Irish study of high-density schizophrenia families (ISHDSF)

The neuregulin-1 gene (NRG1) at chromosome 8p21-22 has been implicated as a schizophrenia susceptibility gene in Icelandic, Scottish, Irish and mixed UK populations. The shared ancestry between these populations led us to investigate the NRG1 polymorphisms and appropriate marker haplotypes for linkage and/or association to schizophrenia in the Irish study of high-density schizophrenia families (ISHDSF). Neither single-point nor multi-point linkage analysis of NRG1 markers gave evidence for linkage independent of our pre-existing findings telomeric on 8p. Analysis of linkage disequilibrium (LD) across the 252 kb interval encompassing the 7 marker core Icelandic/Scottish NRG1 haplotype revealed two separate regions of modest LD, comprising markers SNP8NRG255133, SNP8NRG249130 and SNP8NRG243177 (telomeric) and microsatellites 478B14-428, 420M9-1395, D8S1810 and 420M9-116I12 (centromeric). From single marker analysis by TRANSMIT and FBAT we found no evidence for association with schizophrenia for any marker. Haplotype analysis for the three SNPs in LD region 1 and, separately, the four microsatellites in LD region 2 (analyzed in overlapping 2-marker windows), showed no evidence for overtransmission of specific haplotypes to affected individuals. We therefore conclude that if NRG1 does contain susceptibility alleles for schizophrenia, they impact quite weakly on risk in the ISHDSF.

Evaluation of predicted network modules in yeast metabolism using NMR-based metabolite profiling

Genome-scale metabolic models promise important insights into cell function. However, the definition of pathways and functional network modules within these models, and in the biochemical literature in general, is often based on intuitive reasoning. Although mathematical methods have been proposed to identify modules, which are defined as groups of reactions with correlated fluxes, there is a need for experimental verification. We show here that multivariate statistical analysis of the NMR-derived intra- and extracellular metabolite profiles of single-gene deletion mutants in specific metabolic pathways in the yeast Saccharomyces cerevisiae identified outliers whose profiles were markedly different from those of the other mutants in their respective pathways. Application of flux coupling analysis to a metabolic model of this yeast showed that the deleted gene in an outlying mutant encoded an enzyme that was not part of the same functional network module as the other enzymes in the pathway. We suggest that metabolomic methods such as this, which do not require any knowledge of how a gene deletion might perturb the metabolic network, provide an empirical method for validating and ultimately refining the predicted network structure.

The finding of reduced estimated glomerular filtration rate is associated with increased mortality in a large UK population

BACKGROUND: CKD as defined by KDIGO/KDOQI has been shown to affect ~ 8.5% of the UK population. The prevalence of CKD in the UK is similar to that in the USA, yet incident dialysis rates are dramatically different. This retrospective cohort study investigates the association between reduced kidney function and mortality in a large UK population. METHODS: All serum creatinine results covering Northern Ireland's 1.7 million population were collected between 1 January 2001 and 31 December 2002. Estimated glomerular filtration rates (eGFR) were calculated for all serum creatinine measurements using four-variable MDRD equation (IDMS aligned). Patients were followed up for both all-cause and cardiovascular mortality data until the end of December 2006. Patients on renal replacement therapy were excluded. Subgroup analysis in the 75 345 subjects enrolled within a parallel primary care study permitted additional survival analysis with adjustment for traditional cardiovascular risk factors. RESULTS: A total of 1 967 827 serum creatinine results from 533 798 patients were collected. During the period of follow-up, 59 980 deaths occurred. In multivariate survival analysis, using eGFR as a time-varying covariate, a graded association between CKD (defined by eGFR) and all-cause mortality was identified. Compared with participants with an eGFR of > 60 mL/min/1.73 m(2), the adjusted hazard ratios (and 95% confidence intervals) for participants with an eGFR of 45-59 mL/min/1.73 m(2) was 1.02 (0.99-1.04), an eGFR of 30-44 mL/min/1.73 m(2) was 1.44 (1.40-1.47), an eGFR of 15-29 mL/min/1.73 m(2) was 2.12 (2.05-2.20) and an eGFR of

Label-free detection of nanomolar unmodified single- and double-stranded dna by using surface-enhanced raman spectroscopy on Ag and Au colloids

Unlabelled single- and double-stranded DNA (ssDNA and dsDNA, respectively) has been detected at concentrations =10-9?M by surface-enhanced Raman spectroscopy. Under appropriate conditions the sequences spontaneously adsorbed to the surface of both Ag and Au colloids through their nucleobases; this allowed highly reproducible spectra with good signal-to-noise ratios to be recorded on completely unmodified samples. This eliminated the need to promote absorption by introducing external linkers, such as thiols. The spectra of model ssDNA sequences contained bands of all the bases present and showed systematic changes when the overall base composition was altered. Initial tests also showed that small but reproducible changes could be detected between oligonucleotides with the same bases arranged in a different order. The spectra of five ssDNA sequences that correspond to different strains of the Escherichia coli bacterium were found to be sufficiently composition-dependent so that they could be differentiated without the need for any advanced multivariate data analysis techniques.

Use of NMR metabolomic plasma profiling methodologies to identify illicit growth-promoting administrations

Detection of growth-promoter use in animal production systems still proves to be an analytical challenge despite years of activity in the field. This study reports on the capability of NMR metabolomic profiling techniques to discriminate between plasma samples obtained from cattle treated with different groups of growth-promoting hormones (dexamethasone, prednisolone, oestradiol) based on recorded metabolite profiles. Two methods of NMR analysis were investigated—a Carr–Purcell–Meiboom–Gill (CPMG)-pulse sequence technique and a conventional 1H NMR method using pre-extracted plasma. Using the CPMG method, 17 distinct metabolites could be identified from the spectra. 1H NMR analysis of extracted plasma facilitated identification of 23 metabolites—six more than the alternative method and all within the aromatic region. Multivariate statistical analysis of acquired data from both forms of NMR analysis separated the plasma metabolite profiles into distinct sample cluster sets representative of the different animal study groups. Samples from both sets of corticosteroid-treated animals—dexamethasone and prednisolone—were found to be clustered relatively closely and had similar alterations to identified metabolite panels. Distinctive metabolite profiles, different from those observed within plasma from corticosteroid-treated animal plasma, were observed in oestradiol-treated animals and samples from these animals formed a cluster spatially isolated from control animal plasma samples. These findings suggest the potential use of NMR methodologies of plasma metabolite analysis as a high-throughput screening technique to aid detection of growth promoter use.

Improved prognosis for congenital nephrotic syndrome of the Finnish type in Irish families

Congenital nephrotic syndrome of the Finnish type is a rare autosomal recessive disease with a high infant mortality without aggressive treatment. The biochemical basis of the disease is not understood fully but the disease locus has been mapped recently to chromosome 19q12-q13.1 in Finnish families. This paper describes the clinical features and outcome of 20 patients in Ireland with congenital nephrotic syndrome of the Finnish type who have presented since 1980. Before 1987, all infants died by the age of 3 years. After the introduction of daily intravenous albumin infusion, nutritional support, elective bilateral nephrectomy, and renal transplantation, mortality in the past decade has fallen to 30%, with no deaths in the past five years. Genetic linkage analysis was performed in six families in whom DNA was available and the locus responsible was mapped to the same region on chromosome 19 as in Finnish families, suggesting that Irish families share the same disease locus.

Comorbidities and current smoking status amongst patients starting renal replacement therapy in England, Wales and Northern Ireland from 2003 to 2008: national and centre-specific analyses

Introduction: The prevalence of comorbidities in incident renal replacement therapy (RRT) patients changes with age and varies between ethnic groups. This study describes these associations and the independent effect of comorbidities on outcomes. Methods: Adult patients starting RRT between 2003 and 2008 in centres reporting to the UK Renal Registry (UKRR) with data on comorbidity (n ¼ 14,909) were included. The UKRR studied the association of comorbidity with patient demographics, treatment modality, haemoglobin, renal function at start of RRT and subsequent listing for kidney transplantation. The relationship between comorbidities and mortality at 90 days and one year after 90 days from start of RRT was explored using Cox regression. Results: Completeness of comorbidity data was 40.0% compared with 54.3% in 2003. Of patients with data, 53.8% had one or more comorbidities. Diabetes mellitus and ischaemic heart disease were the most common conditions seen in 30.1% and 22.7% of patients respectively. Current smoking was recorded for 14.5% of incident RRT patients in the 6-year period. Comorbidities became more common with increasing age in all ethnic groups although the difference between the 65–74 and 75+ age groups was not significant. Within each age group, South Asians and Blacks had lower rates of comorbidity, despite higher rates of diabetes mellitus. In multivariate survival analysis, malignancy and ischaemic/neuropathic ulcers were the strongest independent predictors of poor survival at 1 year after 90 days from the start of RRT. Conclusion: Differences in prevalence of comorbid illnesses in incident RRT patients may reflect variation in access to health care or competing risk prior to commencing treatment. At the same time, smoking rates remained high in this ‘at risk’ population. Further work on this and ways to improve comorbidity reporting should be priorities for 2010–11.