George Nicholson

Alternative splicing events are a late feature of pathology in a mouse model of spinal muscular atrophy.

Spinal muscular atrophy is a severe motor neuron disease caused by inactivating mutations in the SMN1 gene leading to reduced levels of full-length functional SMN protein. SMN is a critical mediator of spliceosomal protein assembly, and complete loss or drastic reduction in protein leads to loss of cell viability. However, the reason for selective motor neuron degeneration when SMN is reduced to levels which are tolerated by all other cell types is not currently understood. Widespread splicing abnormalities have recently been reported at end-stage in a mouse model of SMA, leading to the proposition that disruption of efficient splicing is the primary mechanism of motor neuron death. However, it remains unclear whether splicing abnormalities are present during early stages of the disease, which would be a requirement for a direct role in disease pathogenesis. We performed exon-array analysis of RNA from SMN deficient mouse spinal cord at 3 time points, pre-symptomatic (P1), early symptomatic (P7), and late-symptomatic (P13). Compared to littermate control mice, SMA mice showed a time-dependent increase in the number of exons showing differential expression, with minimal differences between genotypes at P1 and P7, but substantial variation in late-symptomatic (P13) mice. Gene ontology analysis revealed differences in pathways associated with neuronal development as well as cellular injury. Validation of selected targets by RT–PCR confirmed the array findings and was in keeping with a shift between physiologically occurring mRNA isoforms. We conclude that the majority of splicing changes occur late in SMA and may represent a secondary effect of cell injury, though we cannot rule out significant early changes in a small number of transcripts crucial to motor neuron survival.

Human metabolic profiles are stably controlled by genetic and environmental variation.

1H Nuclear Magnetic Resonance spectroscopy (1H NMR) is increasingly used to measure metabolite concentrations in sets of biological samples for top‐down systems biology and molecular epidemiology. For such purposes, knowledge of the sources of human variation in metabolite concentrations is valuable, but currently sparse. We conducted and analysed a study to create such a resource. In our unique design, identical and non‐identical twin pairs donated plasma and urine samples longitudinally. We acquired 1H NMR spectra on the samples, and statistically decomposed variation in metabolite concentration into familial (genetic and common‐environmental), individual‐environmental, and longitudinally unstable components. We estimate that stable variation, comprising familial and individual‐environmental factors, accounts on average for 60% (plasma) and 47% (urine) of biological variation in 1H NMR‐detectable metabolite concentrations. Clinically predictive metabolic variation is likely nested within this stable component, so our results have implications for the effective design of biomarker‐discovery studies. We provide a power‐calculation method which reveals that sample sizes of a few thousand should offer sufficient statistical precision to detect 1H NMR‐based biomarkers quantifying predisposition to disease.

Adult Onset Global Loss of the Fto Gene Alters Body Composition and Metabolism in the Mouse

The strongest BMI–associated GWAS locus in humans is the FTO gene. Rodent studies demonstrate a role for FTO in energy homeostasis and body composition. The phenotypes observed in loss of expression studies are complex with perinatal lethality, stunted growth from weaning, and significant alterations in body composition. Thus understanding how and where Fto regulates food intake, energy expenditure, and body composition is a challenge. To address this we generated a series of mice with distinct temporal and spatial loss of Fto expression. Global germline loss of Fto resulted in high perinatal lethality and a reduction in body length, fat mass, and lean mass. When ratio corrected for lean mass, mice had a significant increase in energy expenditure, but more appropriate multiple linear regression normalisation showed no difference in energy expenditure. Global deletion of Fto after the in utero and perinatal period, at 6 weeks of age, removed the high lethality of germline loss. However, there was a reduction in weight by 9 weeks, primarily as loss of lean mass. Over the subsequent 10 weeks, weight converged, driven by an increase in fat mass. There was a switch to a lower RER with no overall change in food intake or energy expenditure. To test if the phenotype can be explained by loss of Fto in the mediobasal hypothalamus, we sterotactically injected adeno-associated viral vectors encoding Cre recombinase to cause regional deletion. We observed a small reduction in food intake and weight gain with no effect on energy expenditure or body composition. Thus, although hypothalamic Fto can impact feeding, the effect of loss of Fto on body composition is brought about by its actions at sites elsewhere. Our data suggest that Fto may have a critical role in the control of lean mass, independent of its effect on food intake.

A reassessment of genetic diversity in Icelanders: Strong evidence from multiple loci for relative homogeneity caused by genetic drift

There has been some controversy in the literature concerning whether Icelanders are genetically homogenous or heterogeneous relative to other European populations. We reassess this question in the light of large data sets spanning 83 autosomal SNP loci, 14 serogenetic loci, 6622 Y‐chromosomes and 3214 sequences from mtDNA hypervariable segments 1 and 2 (HVS1 and HVS2). Our results strongly support the hypothesis that genetic drift, with a consequent loss of variation, has had a greater impact on Icelanders than most other Europeans. We also analyse 7245 HVS1 sequences from 25 European populations. In line with other studies, we observe a deficit of rare HVS1 haplotypes and an excess of intermediate frequency haplotypes in Icelanders compared to most European populations, with some measures of genetic diversity indicating relative heterogeneity and others indicating relative homogeneity of Icelanders. Simulations indicate that genetic drift, and not admixture (as proposed by Árnason, 2003 ) is the most likely cause of the atypical Icelandic HVS1 frequency spectrum. These simulations reveal that gene diversity (heterozygosity) and mean pairwise differences are largely insensitive to events in recent population history, while statistics based on the number of haplotypes or segregating sites are much more sensitive. Overall, our analyses strongly indicate that the Icelandic gene pool is less heterogeneous than those of most other European populations.

The presence of methylation quantitative trait loci indicates a direct genetic influence on the level of DNA methylation in adipose tissue.

Genetic variants that associate with DNA methylation at CpG sites (methylation quantitative trait loci, meQTLs) offer a potential biological mechanism of action for disease associated SNPs. We investigated whether meQTLs exist in abdominal subcutaneous adipose tissue (SAT) and if CpG methylation associates with metabolic syndrome (MetSyn) phenotypes. We profiled 27,718 genomic regions in abdominal SAT samples of 38 unrelated individuals using differential methylation hybridization (DMH) together with genotypes at 5,227,243 SNPs and expression of 17,209 mRNA transcripts. Validation and replication of significant meQTLs was pursued in an independent cohort of 181 female twins. We find that, at 5% false discovery rate, methylation levels of 149 DMH regions associate with at least one SNP in a ±500 kilobase cis-region in our primary study. We sought to validate 19 of these in the replication study and find that five of these significantly associate with the corresponding meQTL SNPs from the primary study. We find that none of the 149 meQTL top SNPs is a significant expression quantitative trait locus in our expression data, but we observed association between expression levels of two mRNA transcripts and cis-methylation status. Our results indicate that DNA CpG methylation in abdominal SAT is partly under genetic control. This study provides a starting point for future investigations of DNA methylation in adipose tissue.

Distinct Developmental Profile of Lower-Body Adipose Tissue Defines Resistance Against Obesity-Associated Metabolic Complications

Upper- and lower-body fat depots exhibit opposing associations with obesity-related metabolic disease. We defined the relationship between DEXA-quantified fat depots and diabetes/cardiovascular risk factors in a healthy population-based cohort (n = 3,399). Gynoid fat mass correlated negatively with insulin resistance after total fat mass adjustment, whereas the opposite was seen for abdominal fat. Paired transcriptomic analysis of gluteal subcutaneous adipose tissue (GSAT) and abdominal subcutaneous adipose tissue (ASAT) was performed across the BMI spectrum (n = 49; 21.4–45.5 kg/m2). In both depots, energy-generating metabolic genes were negatively associated and inflammatory genes were positively associated with obesity. However, associations were significantly weaker in GSAT. At the systemic level, arteriovenous release of the proinflammatory cytokine interleukin-6 (n = 34) was lower from GSAT than ASAT. Isolated preadipocytes retained a depot-specific transcriptional “memory” of embryonic developmental genes and exhibited differential promoter DNA methylation of selected genes (HOTAIR, TBX5) between GSAT and ASAT. Short hairpin RNA–mediated silencing identified TBX5 as a regulator of preadipocyte proliferation and adipogenic differentiation in ASAT. In conclusion, intrinsic differences in the expression of developmental genes in regional adipocytes provide a mechanistic basis for diversity in adipose tissue (AT) function. The less inflammatory nature of lower-body AT offers insight into the opposing metabolic disease risk associations between upper- and lower-body obesity.

Coexpression network analysis in abdominal and gluteal adipose tissue reveals regulatory genetic loci for metabolic syndrome and related phenotypes.

Metabolic Syndrome (MetS) is highly prevalent and has considerable public health impact, but its underlying genetic factors remain elusive. To identify gene networks involved in MetS, we conducted whole-genome expression and genotype profiling on abdominal (ABD) and gluteal (GLU) adipose tissue, and whole blood (WB), from 29 MetS cases and 44 controls. Co-expression network analysis for each tissue independently identified nine, six, and zero MetS–associated modules of coexpressed genes in ABD, GLU, and WB, respectively. Of 8,992 probesets expressed in ABD or GLU, 685 (7.6%) were expressed in ABD and 51 (0.6%) in GLU only. Differential eigengene network analysis of 8,256 shared probesets detected 22 shared modules with high preservation across adipose depots (DABD-GLU = 0.89), seven of which were associated with MetS (FDR P<0.01). The strongest associated module, significantly enriched for immune response–related processes, contained 94/620 (15%) genes with inter-depot differences. In an independent cohort of 145/141 twins with ABD and WB longitudinal expression data, median variability in ABD due to familiality was greater for MetS–associated versus un-associated modules (ABD: 0.48 versus 0.18, P = 0.08; GLU: 0.54 versus 0.20, P = 7.8×10−4). Cis-eQTL analysis of probesets associated with MetS (FDR P<0.01) and/or inter-depot differences (FDR P<0.01) provided evidence for 32 eQTLs. Corresponding eSNPs were tested for association with MetS–related phenotypes in two GWAS of >100,000 individuals; rs10282458, affecting expression of RARRES2 (encoding chemerin), was associated with body mass index (BMI) (P = 6.0×10−4); and rs2395185, affecting inter-depot differences of HLA-DRB1 expression, was associated with high-density lipoprotein (P = 8.7×10−4) and BMI–adjusted waist-to-hip ratio (P = 2.4×10−4). Since many genes and their interactions influence complex traits such as MetS, integrated analysis of genotypes and coexpression networks across multiple tissues relevant to clinical traits is an efficient strategy to identify novel associations.

Inbreeding avoidance under different null models of random mating in the great tit

1. In populations where inbreeding causes a substantial decrease in fitness, selection is expected to favour the evolution of inbreeding avoidance behaviours. Elsewhere we have documented substantial inbreeding depression and the importance of dispersal in avoiding inbreeding in a long-term population study of the great tit Parus major in Wytham (UK). In this study, we ask whether individuals from this population actively avoid mating with kin. 2. We generated four contrasting models of random mate choice that assumed varying levels of mate availability in each year of the data set. This allowed us to compare observed and simulated distributions and frequencies of inbreeding coefficients from 41 years of breeding data. 3. We found no evidence that birds avoid mating with related partners. Our results show that birds breed more often with relatives than expected under null models of mate choice that lack population structure, but not when compared to scenarios where birds were mated with their nearest neighbours. Pedigree-derived F(IS) values were positive for all scenarios of random mating, confirming the lack of inbreeding avoidance in this population. 4. These results imply the existence of spatial genetic structure where related individuals occur closer together than nonrelated individuals while breeding, and suggest that the relatedness between breeding individuals of the opposite sex decreases with distance. Thus, while dispersal from the natal site decreases the number of relatives around an individual, it does not completely homogenize genetic structure. 5. We show that brother-sister pairs are observed more often than under any scenario of random mating, suggesting that not only birds do not avoid mating with kin, but also that the apparently maladaptive choice of mating with a sibling is made more often than expected. 6. Our results provide no evidence to suggest that individuals actively avoid kin. In fact, some types of inbreeding occur more often than expected, despite the substantial fitness costs. The observed lack of inbreeding avoidance is in agreement with other studies of non-cooperatively breeding passerine birds, although the higher than expected frequency of sibling mating remains a puzzling result.

MicroRNA Expression in Abdominal and Gluteal Adipose Tissue Is Associated with mRNA Expression Levels and Partly Genetically Driven

To understand how miRNAs contribute to the molecular phenotype of adipose tissues and related traits, we performed global miRNA expression profiling in subcutaneous abdominal and gluteal adipose tissue of 70 human subjects and characterised which miRNAs were differentially expressed between these tissues. We found that 12% of the miRNAs were significantly differentially expressed between abdominal and gluteal adipose tissue (FDR adjusted p<0.05) in the primary study, of which 59 replicated in a follow-up study of 40 additional subjects. Further, 14 miRNAs were found to be associated with metabolic syndrome case-control status in abdominal tissue and three of these replicated (primary study: FDR adjusted p<0.05, replication: p<0.05 and directionally consistent effect). Genome-wide genotyping was performed in the 70 subjects to enable miRNA expression quantitative trait loci (eQTL) analysis. Candidate miRNA eQTLs were followed-up in the additional 40 subjects and six significant, independent cis-located miRNA eQTLs (primary study: p<0.001; replication: p<0.05 and directionally consistent effect) were identified. Finally, global mRNA expression profiling was performed in both tissues to enable association analysis between miRNA and target mRNA expression levels. We find 22% miRNAs in abdominal and 9% miRNAs in gluteal adipose tissue with expression levels significantly associated with the expression of corresponding target mRNAs (FDR adjusted p<0.05). Taken together, our results indicate a clear difference in the miRNA molecular phenotypic profile of abdominal and gluteal adipose tissue, that the expressions of some miRNAs are influenced by cis-located genetic variants and that miRNAs are associated with expression levels of their predicted mRNA targets.

Comparative protein profiling of serum and plasma using an antibody suspension bead array approach

In the pursuit towards a systematic analysis of human diseases, array‐based approaches within antibody proteomics offer high‐throughput strategies to discover protein biomarkers in serum and plasma. To investigate the influence of sample preparation on such discovery attempts, we report on a systematic effort to compare serum and plasma protein profiles determined with an antibody suspension bead array. The intensity levels were used to define protein profiles and no significant differences between serum and plasma were observed for 79% of the 174 antibodies (targeting 156 proteins). By excluding 36 antibodies giving rise to differential intensity levels, cluster analysis revealed donor‐specific rather than preparation‐dependent grouping. With a cohort from a clinically relevant medical condition, the metabolic syndrome, the influence of the sample type on a multiplexed biomarker discovery approach was further investigated. Independent comparisons of protein profiles in serum and plasma revealed an antibody targeting ADAMTSL‐4, a protein that would qualify to be studied further in association with the condition. In general, the preparation type had an impact on the results of the applied antibody suspension bead array, and while the technical variability was equal, plasma offered a greater biological variability and allowed to give rise to more discoveries than serum.