Family-based association analysis with ordered categorical phenotypes, covariates and interactions

Genetic association analyses of family-based studies with ordered categorical phenotypes are often conducted using methods either for quantitative or for binary traits, which can lead to suboptimal analyses. Here we present an alternative likelihood-based method of analysis for single nucleotide polymorphism (SNP) genotypes and ordered categorical phenotypes in nuclear families of any size. Our approach, which extends our previous work for binary phenotypes, permits straightforward inclusion of covariate, gene-gene and gene-covariate interaction terms in the likelihood, incorporates a simple model for ascertainment and allows for family-specific effects in the hypothesis test. Additionally, our method produces interpretable parameter estimates and valid confidence intervals. We assess the proposed method using simulated data, and apply it to a polymorphism in the c-reactive protein (CRP) gene typed in families collected to investigate human systemic lupus erythematosus. By including sex interactions in the analysis, we show that the polymorphism is associated with anti-nuclear autoantibody (ANA) production in females, while there appears to be no effect in males.

Host niches and defensive extended phenotypes structure parasitoid wasp communities

Oak galls are spectacular extended phenotypes of gallwasp genes in host oak tissues and have evolved complex morphologies that serve, in part, to exclude parasitoid natural enemies. Parasitoids and their insect herbivore hosts have coevolved to produce diverse communities comprising about a third of all animal species. The factors structuring these communities, however, remain poorly understood. An emerging theme in community ecology is the need to consider the effects of host traits, shaped by both natural selection and phylogenetic history, on associated communities of natural enemies. Here we examine the impact of host traits and phylogenetic relatedness on 48 ecologically closed and species-rich communities of parasitoids attacking gall-inducing wasps on oaks. Gallwasps induce the development of spectacular and structurally complex galls whose species- and generation-specific morphologies are the extended phenotypes of gallwasp genes. All the associated natural enemies attack their concealed hosts through gall tissues, and several structural gall traits have been shown to enhance defence against parasitoid attack. Here we explore the significance of these and other host traits in predicting variation in parasitoid community structure across gallwasp species. In particular, we test the "Enemy Hypothesis,'' which predicts that galls with similar morphology will exclude similar sets of parasitoids and therefore have similar parasitoid communities. Having controlled for phylogenetic patterning in host traits and communities, we found significant correlations between parasitoid community structure and several gall structural traits (toughness, hairiness, stickiness), supporting the Enemy Hypothesis. Parasitoid community structure was also consistently predicted by components of the hosts' spatiotemporal niche, particularly host oak taxonomy and gall location (e.g., leaf versus bud versus seed). The combined explanatory power of structural and spatiotemporal traits on community structure can be high, reaching 62% in one analysis. The observed patterns derive mainly from partial niche specialisation of highly generalist parasitoids with broad host ranges (>20 hosts), rather than strict separation of enemies with narrower host ranges, and so may contribute to maintenance of the richness of generalist parasitoids in gallwasp communities. Though evolutionary escape from parasitoids might most effectively be achieved via changes in host oak taxon, extreme conservatism in this trait for gallwasps suggests that selection is more likely to have acted on gall morphology and location. Any escape from parasitoids associated with evolutionary shifts in these traits has probably only been transient, however, due to subsequent recruitment of parasitoid species already attacking other host galls with similar trait combinations.

An intronic growth hormone receptor mutation causing activation of a pseudoexon is associated with a broad spectrum of growth hormone insensitivity phenotypes

Context: Inherited GH insensitivity (GHI) is usually caused by mutations in the GH receptor (GHR). Patients present with short stature associated with high GH and low IGF-I levels and may have midfacial hypoplasia ( typical Laron syndrome facial features). We previously described four mildly affected GHI patients with an intronic mutation in the GHR gene (A.(1) -> G.(1) substitution in intron 6), resulting in the activation of a pseudoexon (6 Psi) and inclusion of 36 amino acids. Objective: The study aimed to analyze the clinical and genetic characteristics of additional GHI patients with the pseudoexon (6 Psi) mutation. Design/Patients: Auxological, biochemical, genetic, and haplotype data from seven patients with severe short stature and biochemical evidence of GHI were assessed. Main Outcome Measures: We assessed genotype-phenotype relationship. Results: One patient belongs to the same extended family, previously reported. She has normal facial features, and her IGF-I levels are in the low-normal range for age. The six unrelated patients, four of whom have typical Laron syndrome facial features, have heights ranging from -3.3 to -6.0 SD and IGF-I levels that vary from normal to undetectable. We hypothesize that the marked difference in biochemical and clinical phenotypes might be caused by variations in the splicing efficiency of the pseudoexon. Conclusions: Activation of the pseudoexon in the GHR gene can lead to a variety of GHI phenotypes. Therefore, screening for the presence of this mutation should be performed in all GHI patients without mutations in the coding exons.

Variation in antibiotic-induced microbial recolonization impacts on the host metabolic phenotypes of rats

The interaction between the gut microbiota and their mammalian host is known to have far-reaching consequences with respect to metabolism and health. We investigated the effects of eight days of oral antibiotic exposure (penicillin and streptomycin sulfate) on gut microbial composition and host metabolic phenotype in male Han-Wistar rats (n = 6) compared to matched controls. Early recolonization was assessed in a third group exposed to antibiotics for four days followed by four days recovery (n = 6). Fluorescence in situ hybridization analysis of the intestinal contents collected at eight days showed a significant reduction in all bacterial groups measured (control, 1010.7 cells/g feces; antibiotic-treated, 108.4). Bacterial suppression reduced the excretion of mammalian-microbial urinary cometabolites including hippurate, phenylpropionic acid, phenylacetylglycine and indoxyl-sulfate whereas taurine, glycine, citrate, 2-oxoglutarate, and fumarate excretion was elevated. While total bacterial counts remained notably lower in the recolonized animals (109.1 cells/g faeces) compared to the controls, two cage-dependent subgroups emerged with Lactobacillus/Enterococcus probe counts dominant in one subgroup. This dichotomous profile manifested in the metabolic phenotypes with subgroup differences in tricarboxylic acid cycle metabolites and indoxyl-sulfate excretion. Fecal short chain fatty acids were diminished in all treated animals. Antibiotic treatment induced a profound effect on the microbiome structure, which was reflected in the metabotype. Moreover, the recolonization process was sensitive to the microenvironment, which may impact on understanding downstream consequences of antibiotic consumption in human populations.

Systemic multicompartmental effects of the gut microbiome on mouse metabolic phenotypes

To characterize the impact of gut microbiota on host metabolism, we investigated the multicompartmental metabolic profiles of a conventional mouse strain (C3H/HeJ) (n=5) and its germ-free (GF) equivalent (n=5). We confirm that the microbiome strongly impacts on the metabolism of bile acids through the enterohepatic cycle and gut metabolism (higher levels of phosphocholine and glycine in GF liver and marked higher levels of bile acids in three gut compartments). Furthermore we demonstrate that (1) well-defined metabolic differences exist in all examined compartments between the metabotypes of GF and conventional mice: bacterial co-metabolic products such as hippurate (urine) and 5-aminovalerate (colon epithelium) were found at reduced concentrations, whereas raffinose was only detected in GF colonic profiles. (2) The microbiome also influences kidney homeostasis with elevated levels of key cell volume regulators (betaine, choline, myo-inositol and so on) observed in GF kidneys. (3) Gut microbiota modulate metabotype expression at both local (gut) and global (biofluids, kidney, liver) system levels and hence influence the responses to a variety of dietary modulation and drug exposures relevant to personalized health-care investigations.

Pharmacometabonomic characterization of xenobiotic and endogenous metabolic phenotypes that account for inter-individual variation in isoniazid-induced toxicological response

An NMR-based pharmacometabonomic approach was applied to investigate inter-animal variation in response to isoniazid (INH; 200 and 400 mg/kg) in male Sprague-Dawley rats, alongside complementary clinical chemistry and histopathological analysis. Marked inter-animal variability in central nervous system (CNS) toxicity was identified following administration of a high dose of INH, which enabled characterization of CNS responders and CNS non-responders. High-resolution post-dose urinary (1)H NMR spectra were modeled both by their xenobiotic and endogenous metabolic information sets, enabling simultaneous identification of the differential metabolic fate of INH and its associated endogenous metabolic consequences in CNS responders and CNS non-responders. A characteristic xenobiotic metabolic profile was observed for CNS responders, which revealed higher urinary levels of pyruvate isonicotinylhydrazone and β-glucosyl isonicotinylhydrazide and lower levels of acetylisoniazid compared to CNS non-responders. This suggested that the capacity for acetylation of INH was lower in CNS responders, leading to increased metabolism via conjugation with pyruvate and glucose. In addition, the endogenous metabolic profile of CNS responders revealed higher urinary levels of lactate and glucose, in comparison to CNS non-responders. Pharmacometabonomic analysis of the pre-dose (1)H NMR urinary spectra identified a metabolic signature that correlated with the development of INH-induced adverse CNS effects and may represent a means of predicting adverse events and acetylation capacity when challenged with high dose INH. Given the widespread use of INH for the treatment of tuberculosis, this pharmacometabonomic screening approach may have translational potential for patient stratification to minimize adverse events.

Association between serotonin 5-HT-2C receptor gene (HTR2C) polymorphisms and obesity- and mental health-related phenotypes in a large population-based cohort

OBJECTIVE: Studies have shown that common single-nucleotide polymorphisms (SNPs) in the serotonin 5-HT-2C receptor (HTR2C) are associated with antipsychotic agent-induced weight gain and the development of behavioural and psychological symptoms. We aimed to analyse whether variation in the HTR2C is associated with obesity- and mental health-related phenotypes in a large population-based cohort. METHOD: Six tagSNPs, which capture all common genetic variation in the HTR2C gene, were genotyped in 4978 men and women from the European Prospective Investigation into Cancer (EPIC)-Norfolk study, an ongoing prospective population-based cohort study in the United Kingdom. To confirm borderline significant associations, the -759C/T SNP (rs3813929) was genotyped in the remaining 16 003 individuals from the EPIC-Norfolk study. We assessed social and psychological circumstances using the Health and Life Experiences Questionnaire. Genmod models were used to test associations between the SNPs and the outcomes. Logistic regression was performed to test for association of SNPs with obesity- and mental health- related phenotypes. RESULTS: Of the six HTR2C SNPs, only the T allele of the -759C/T SNP showed borderline significant associations with higher body mass index (BMI) (0.23 kg m(-2); (95% confidence interval (CI): 0.01-0.44); P=0.051) and increased risk of lifetime major depressive disorder (MDD) (Odds ratio (OR): 1.13 (95% CI: 1.01-1.22), P=0.02). The associations between the -759C/T and BMI and lifetime MDD were independent. As associations only achieved borderline significance, we aimed to validate our findings on the -759C/T SNP in the full EPIC-Norfolk cohort (n=20 981). Although the association with BMI remained borderline significant (beta=0.20 kg m(-2); 95% CI: 0.04-0.44, P=0.09), that with lifetime MDD (OR: 1.01; 95% CI: 0.94-1.09, P=0.73) was not replicated. CONCLUSIONS: Our findings suggest that common HTR2C gene variants are unlikely to have a major role in obesity- and mental health-related traits in the general population.

Population structure and associated phenotypes of Salmonella enterica serovars Derby and Mbandaka overlap with host range

BACKGROUND:The Salmonella enterica serovar Derby is frequently isolated from pigs and turkeys whereas serovar Mbandaka is frequently isolated from cattle, chickens and animal feed in the UK. Through comparative genomics, phenomics and mutant construction we previously suggested possible mechanistic reasons why these serovars demonstrate apparently distinct host ranges. Here, we investigate the genetic and phenotypic diversity of these two serovars in the UK. We produce a phylogenetic reconstruction and perform several biochemical assays on isolates of S. Derby and S. Mbandaka acquired from sites across the UK between the years 2000 and 2010. RESULTS:We show that UK isolates of S. Mbandaka comprise of one clonal lineage which is adapted to proficient utilisation of metabolites found in soya beans under ambient conditions. We also show that this clonal lineage forms a biofilm at 25 °C, suggesting that this serovar maybe well adapted to survival ex vivo, growing in animal feed. Conversely, we show that S. Derby is made of two distinct lineages, L1 and L2. These lineages differ genotypically and phenotypically, being divided by the presence and absence of SPI-23 and the ability to more proficiently invade porcine jejunum derived cell line IPEC-J2. CONCLUSION:The results of this study lend support to the hypothesis that the differences in host ranges of S. Derby and S. Mbandaka are adaptations to pathogenesis, environmental persistence, as well as utilisation of metabolites abundant in their respective host environments.

Distinct behavioral phenotypes in male mice lacking the thyroid hormone receptor alpha1 or beta isoforms

Thyroid hormones influence both neuronal development and anxiety via the thyroid hormone receptors (TRs). The TRs are encoded by two different genes, TRalpha and TRbeta. The loss of TRalpha1 is implicated in increased anxiety in males, possibly via a hippocampal increase in GABAergic activity. We compared both social behaviors and two underlying and related non-social behaviors, state anxiety and responses to acoustic and tactile startle in the gonadally intact TRalpha1 knockout (alpha1KO) and TRbeta (betaKO) male mice to their wild-type counterparts. For the first time, we show an opposing effect of the two TR isoforms, TRalpha1 and TRbeta, in the regulation of state anxiety, with alpha1 knockout animals (alpha1KO) showing higher levels of anxiety and betaKO males showing less anxiety compared to respective wild-type mice. At odds with the increased anxiety in non-social environments, alpha1KO males also show lower levels of responsiveness to acoustic and tactile startle stimuli. Consistent with the data that T4 is inhibitory to lordosis in female mice, we show subtly increased sex behavior in alpha1KO male mice. These behaviors support the idea that TRalpha1 could be inhibitory to ERalpha driven transcription that ultimately impacts ERalpha driven behaviors such as lordosis. The behavioral phenotypes point to novel roles for the TRs, particularly in non-social behaviors such as state anxiety and startle.

Periodic Reordering

For many networks in nature, science and technology, it is possible to order the nodes so that most links are short-range, connecting near-neighbours, and relatively few long-range links, or shortcuts, are present. Given a network as a set of observed links (interactions), the task of finding an ordering of the nodes that reveals such a range-dependent structure is closely related to some sparse matrix reordering problems arising in scientific computation. The spectral, or Fiedler vector, approach for sparse matrix reordering has successfully been applied to biological data sets, revealing useful structures and subpatterns. In this work we argue that a periodic analogue of the standard reordering task is also highly relevant. Here, rather than encouraging nonzeros only to lie close to the diagonal of a suitably ordered adjacency matrix, we also allow them to inhabit the off-diagonal corners. Indeed, for the classic small-world model of Watts & Strogatz (1998, Collective dynamics of ‘small-world’ networks. Nature, 393, 440–442) this type of periodic structure is inherent. We therefore devise and test a new spectral algorithm for periodic reordering. By generalizing the range-dependent random graph class of Grindrod (2002, Range-dependent random graphs and their application to modeling large small-world proteome datasets. Phys. Rev. E, 66, 066702-1–066702-7) to the periodic case, we can also construct a computable likelihood ratio that suggests whether a given network is inherently linear or periodic. Tests on synthetic data show that the new algorithm can detect periodic structure, even in the presence of noise. Further experiments on real biological data sets then show that some networks are better regarded as periodic than linear. Hence, we find both qualitative (reordered networks plots) and quantitative (likelihood ratios) evidence of periodicity in biological networks.

Immunogenetics of drug-induced skin blistering disorders. Part II: Synthesis

The overall immunopathogenesis relevant to a large series of disorders caused by a drug or its associated hyperimmune condition is discussed based upon examining the genetics of severe drug-induced bullous skin problems (sporadic idiosyncratic adverse events including Stevens-Johnson syndrome and Toxic epidermal necrolysis). New results from an exemplar study on shared precipitating and perpetuating inner causes with other related disease phenotypes including aphtous stomatitis, Behcets, erythema multiforme, Hashimoto's thyroiditis, pemphigus, periodic fevers, Sweet's syndrome and drug-induced multisystem hypersensitivity are presented. A call for a collaborative, wider demographic profiling and deeper immunotyping in suggested future work is made.

Functional genomics in zebrafish permits rapid characterization of novel platelet membrane proteins

In this study, we demonstrate the suitability of the vertebrate Danio rerio (zebrafish) for functional screening of novel platelet genes in vivo by reverse genetics. Comparative transcript analysis of platelets and their precursor cell, the megakaryocyte, together with nucleated blood cell elements, endothelial cells, and erythroblasts, identified novel platelet membrane proteins with hitherto unknown roles in thrombus formation. We determined the phenotype induced by antisense morpholino oligonucleotide (MO)–based knockdown of 5 of these genes in a laser-induced arterial thrombosis model. To validate the model, the genes for platelet glycoprotein (GP) IIb and the coagulation protein factor VIII were targeted. MO-injected fish showed normal thrombus initiation but severely impaired thrombus growth, consistent with the mouse knockout phenotypes, and concomitant knockdown of both resulted in spontaneous bleeding. Knockdown of 4 of the 5 novel platelet proteins altered arterial thrombosis, as demonstrated by modified kinetics of thrombus initiation and/or development. We identified a putative role for BAMBI and LRRC32 in promotion and DCBLD2 and ESAM in inhibition of thrombus formation. We conclude that phenotypic analysis of MO-injected zebrafish is a fast and powerful method for initial screening of novel platelet proteins for function in thrombosis.

Developmental programming: Deficits in reproductive hormone dynamics and ovulatory outcomes in prenatal, testosterone-treated sheep

Prenatal testosterone excess leads to neuroendocrine, ovarian, and metabolic disruptions, culminating in reproductive phenotypes mimicking that of women with polycystic ovary syndrome (PCOS). The objective of this study was to determine the consequences of prenatal testosterone treatment on periovulatory hormonal dynamics and ovulatory outcomes. To generate prenatal testosterone-treated females, pregnant sheep were injected intramuscularly (days 30-90 of gestation, term = 147 days) with 100 mg of testosterone-propionate in cottonseed oil semi-weekly. Female offspring born to untreated control females and prenatal testosterone-treated females were then studied during their first two breeding seasons. Sheep were given two injections of prostaglandin F-2alpha 11 days apart, and blood samples were collected at 2-h intervals for 120 h, 10-min intervals for 8 h during the luteal phase (first breeding season only), and daily for an additional 15 days to characterize changes in reproductive hormonal dynamics. During the first breeding season, prenatal testosterone-treated females manifested disruptions in the timing and magnitude of primary gonadotropin surges, luteal defects, and reduced responsiveness to progesterone negative feedback. Disruptions in the periovulatory sequence of events during the second breeding season included: 1) delayed but increased preovulatory estradiol rise, 2) delayed and severely reduced primary gonadotropin surge in prenatal testosterone-treated females having an LH surge, 3) tendency for an amplified secondary FSH surge and a shift in the relative balance of FSH regulatory proteins, and 4) luteal responses that ranged from normal to anovulatory. These outcomes are likely to be of relevance to developmental origin of infertility disorders and suggest that differences in fetal exposure or fetal susceptibility to testosterone may account for the variability in reproductive phenotypes.

Mutational activation of niche-specific genes provides insight into regulatory networks and bacterial function in a complex environment

The genome of the plant-colonizing bacterium Pseudomonas fluorescens SBW25 harbors a subset of genes that are expressed specifically on plant surfaces. The function of these genes is central to the ecological success of SBW25, but their study poses significant challenges because no phenotype is discernable in vitro. Here, we describe a genetic strategy with general utility that combines suppressor analysis with IVET (SPyVET) and provides a means of identifying regulators of niche-specific genes. Central to this strategy are strains carrying operon fusions between plant environment-induced loci (EIL) and promoterless 'dapB. These strains are prototrophic in the plant environment but auxotrophic on laboratory minimal medium. Regulatory elements were identified by transposon mutagenesis and selection for prototrophs on minimal medium. Approximately 106 mutants were screened for each of 27 strains carrying 'dapB fusions to plant EIL and the insertion point for the transposon determined in approximately 2,000 putative regulator mutants. Regulators were functionally characterized and used to provide insight into EIL phenotypes. For one strain carrying a fusion to the cellulose-encoding wss operon, five different regulators were identified including a diguanylate cyclase, the flagella activator, FleQ, and alginate activator, AmrZ (AlgZ). Further rounds of suppressor analysis, possible by virtue of the SPyVET strategy, revealed an additional two regulators including the activator AlgR, and allowed the regulatory connections to be determined.

In vivo selection for metastasis promoting genes in the mouse

Here, we report the identification of a metastasis promoting factor by a forward genetic screen in mice. A retroviral cDNA library was introduced into the nonmetastatic cancer cell line 168FARN, which was then orthotopically transplanted into mouse mammary fat pads, followed by selection for cells that metastasize to the lung. The genes encoding the disulfide isomerase ERp5 and beta-catenin were found to promote breast cancer invasion and metastasis. Disulfide isomerases (thiol isomerases), which catalyze disulfide bond formation, reduction, and isomerization, have not previously been implicated in cancer cell signaling and tumor metastasis. Overexpression of ERp5 promotes both in vitro migration and invasion and in vivo metastasis of breast cancer cells. These effects were shown to involve activation of ErbB2 and phosphoinositicle 3-kinase (PI3K) pathways through dimerization of ErbB2. Activation of ErbB2 and PI3K subsequently stimulates RhoA and beta-catenin, which mediate the migration and invasion of tumor cells. Inhibition of ErbB2 and PI3K reverses the phenotypes induced by ERp5. Finally, ERp5 was shown to be up-regulated in human surgical samples of invasive breast cancers. These data identify a link between disulfide isomerases and tumor development, and provide a mechanism that modulates ErbB2 and PI3K signaling in the promotion of cancer progression.