Assessment of LD matrix measures for the analysis of biological pathway association.

Complex diseases will have multiple functional sites, and it will be invaluable to understand the cross-locus interaction in terms of linkage disequilibrium (LD) between those sites (epistasis) in addition to the haplotype-LD effects. We investigated the statistical properties of a class of matrix-based statistics to assess this epistasis. These statistical methods include two LD contrast tests (Zaykin et al., 2006) and partial least squares regression (Wang et al., 2008). To estimate Type 1 error rates and power, we simulated multiple two-variant disease models using the SIMLA software package. SIMLA allows for the joint action of up to two disease genes in the simulated data with all possible multiplicative interaction effects between them. Our goal was to detect an interaction between multiple disease-causing variants by means of their linkage disequilibrium (LD) patterns with other markers. We measured the effects of marginal disease effect size, haplotype LD, disease prevalence and minor allele frequency have on cross-locus interaction (epistasis). In the setting of strong allele effects and strong interaction, the correlation between the two disease genes was weak (r=0.2). In a complex system with multiple correlations (both marginal and interaction), it was difficult to determine the source of a significant result. Despite these complications, the partial least squares and modified LD contrast methods maintained adequate power to detect the epistatic effects; however, for many of the analyses we often could not separate interaction from a strong marginal effect. While we did not exhaust the entire parameter space of possible models, we do provide guidance on the effects that population parameters have on cross-locus interaction.

Gait and behavior in an IL1β-mediated model of rat knee arthritis and effects of an IL1 antagonist.

Interleukin-1 beta (IL1β) is a proinflammatory cytokine that mediates arthritic pathologies. Our objectives were to evaluate pain and limb dysfunction resulting from IL1β over-expression in the rat knee and to investigate the ability of local IL1 receptor antagonist (IL1Ra) delivery to reverse-associated pathology. IL1β over-expression was induced in the right knees of 30 Wistar rats via intra-articular injection of rat fibroblasts retrovirally infected with human IL1β cDNA. A subset of animals received a 30 µl intra-articular injection of saline or human IL1Ra on day 1 after cell delivery (0.65 µg/µl hIL1Ra, n = 7 per group). Joint swelling, gait, and sensitivity were investigated over 1 week. On day 8, animals were sacrificed and joints were collected for histological evaluation. Joint inflammation and elevated levels of endogenous IL1β were observed in knees receiving IL1β-infected fibroblasts. Asymmetric gaits favoring the affected limb and heightened mechanical sensitivity (allodynia) reflected a unilateral pathology. Histopathology revealed cartilage loss on the femoral groove and condyle of affected joints. Intra-articular IL1Ra injection failed to restore gait and sensitivity to preoperative levels and did not reduce cartilage degeneration observed in histopathology. Joint swelling and degeneration subsequent to IL1β over-expression is associated limb hypersensitivity and gait compensation. Intra-articular IL1Ra delivery did not result in marked improvement for this model; this may be driven by rapid clearance of administered IL1Ra from the joint space. These results motivate work to further investigate the behavioral consequences of monoarticular arthritis and sustained release drug delivery strategies for the joint space.

Binding of MetJ repressor to specific and nonspecific DNA and effect of S-adenosylmethionine on these interactions.

We have used analytical ultracentrifugation to characterize the binding of the methionine repressor protein, MetJ, to synthetic oligonucleotides containing zero to five specific recognition sites, called metboxes. For all lengths of DNA studied, MetJ binds more tightly to repeats of the consensus sequence than to naturally occurring metboxes, which exhibit a variable number of deviations from the consensus. Strong cooperative binding occurs only in the presence of two or more tandem metboxes, which facilitate protein-protein contacts between adjacent MetJ dimers, but weak affinity is detected even with DNA containing zero or one metbox. The affinity of MetJ for all of the DNA sequences studied is enhanced by the addition of SAM, the known cofactor for MetJ in the cell. This effect extends to oligos containing zero or one metbox, both of which bind two MetJ dimers. In the presence of a large excess concentration of metbox DNA, the effect of cooperativity is to favor populations of DNA oligos bound by two or more MetJ dimers rather than a stochastic redistribution of the repressor onto all available metboxes. These results illustrate the dynamic range of binding affinity and repressor assembly that MetJ can exhibit with DNA and the effect of the corepressor SAM on binding to both specific and nonspecific DNA.

One sound or two? Object-related negativity indexes echo perception.

The ability to isolate a single sound source among concurrent sources and reverberant energy is necessary for understanding the auditory world. The precedence effect describes a related experimental finding, that when presented with identical sounds from two locations with a short onset asynchrony (on the order of milliseconds), listeners report a single source with a location dominated by the lead sound. Single-cell recordings in multiple animal models have indicated that there are low-level mechanisms that may contribute to the precedence effect, yet psychophysical studies in humans have provided evidence that top-down cognitive processes have a great deal of influence on the perception of simulated echoes. In the present study, event-related potentials evoked by click pairs at and around listeners' echo thresholds indicate that perception of the lead and lag sound as individual sources elicits a negativity between 100 and 250 msec, previously termed the object-related negativity (ORN). Even for physically identical stimuli, the ORN is evident when listeners report hearing, as compared with not hearing, a second sound source. These results define a neural mechanism related to the conscious perception of multiple auditory objects.

Inhibition of the anaphase-promoting complex by the Xnf7 ubiquitin ligase.

Degradation of specific protein substrates by the anaphase-promoting complex/cyclosome (APC) is critical for mitotic exit. We have identified the protein Xenopus nuclear factor 7 (Xnf7) as a novel APC inhibitor able to regulate the timing of exit from mitosis. Immunodepletion of Xnf7 from Xenopus laevis egg extracts accelerated the degradation of APC substrates cyclin B1, cyclin B2, and securin upon release from cytostatic factor arrest, whereas excess Xnf7 inhibited APC activity. Interestingly, Xnf7 exhibited intrinsic ubiquitin ligase activity, and this activity was required for APC inhibition. Unlike other reported APC inhibitors, Xnf7 did not associate with Cdc20, but rather bound directly to core subunits of the APC. Furthermore, Xnf7 was required for spindle assembly checkpoint function in egg extracts. These data suggest that Xnf7 is an APC inhibitor able to link spindle status to the APC through direct association with APC core components.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Association between DNA damage response and repair genes and risk of invasive serous ovarian cancer.

BACKGROUND: We analyzed the association between 53 genes related to DNA repair and p53-mediated damage response and serous ovarian cancer risk using case-control data from the North Carolina Ovarian Cancer Study (NCOCS), a population-based, case-control study. METHODS/PRINCIPAL FINDINGS: The analysis was restricted to 364 invasive serous ovarian cancer cases and 761 controls of white, non-Hispanic race. Statistical analysis was two staged: a screen using marginal Bayes factors (BFs) for 484 SNPs and a modeling stage in which we calculated multivariate adjusted posterior probabilities of association for 77 SNPs that passed the screen. These probabilities were conditional on subject age at diagnosis/interview, batch, a DNA quality metric and genotypes of other SNPs and allowed for uncertainty in the genetic parameterizations of the SNPs and number of associated SNPs. Six SNPs had Bayes factors greater than 10 in favor of an association with invasive serous ovarian cancer. These included rs5762746 (median OR(odds ratio)(per allele) = 0.66; 95% credible interval (CI) = 0.44-1.00) and rs6005835 (median OR(per allele) = 0.69; 95% CI = 0.53-0.91) in CHEK2, rs2078486 (median OR(per allele) = 1.65; 95% CI = 1.21-2.25) and rs12951053 (median OR(per allele) = 1.65; 95% CI = 1.20-2.26) in TP53, rs411697 (median OR (rare homozygote) = 0.53; 95% CI = 0.35 - 0.79) in BACH1 and rs10131 (median OR( rare homozygote) = not estimable) in LIG4. The six most highly associated SNPs are either predicted to be functionally significant or are in LD with such a variant. The variants in TP53 were confirmed to be associated in a large follow-up study. CONCLUSIONS/SIGNIFICANCE: Based on our findings, further follow-up of the DNA repair and response pathways in a larger dataset is warranted to confirm these results.

Sequencing three crocodilian genomes to illuminate the evolution of archosaurs and amniotes.

The International Crocodilian Genomes Working Group (ICGWG) will sequence and assemble the American alligator (Alligator mississippiensis), saltwater crocodile (Crocodylus porosus) and Indian gharial (Gavialis gangeticus) genomes. The status of these projects and our planned analyses are described.

Finding our way through phenotypes.

Despite a large and multifaceted effort to understand the vast landscape of phenotypic data, their current form inhibits productive data analysis. The lack of a community-wide, consensus-based, human- and machine-interpretable language for describing phenotypes and their genomic and environmental contexts is perhaps the most pressing scientific bottleneck to integration across many key fields in biology, including genomics, systems biology, development, medicine, evolution, ecology, and systematics. Here we survey the current phenomics landscape, including data resources and handling, and the progress that has been made to accurately capture relevant data descriptions for phenotypes. We present an example of the kind of integration across domains that computable phenotypes would enable, and we call upon the broader biology community, publishers, and relevant funding agencies to support efforts to surmount today's data barriers and facilitate analytical reproducibility.

Chlamydia trachomatis Infection Leads to Defined Alterations to the Lipid Droplet Proteome in Epithelial Cells.

The obligate intracellular bacterium Chlamydia trachomatis is a major human pathogen and a main cause of genital and ocular diseases. During its intracellular cycle, C. trachomatis replicates inside a membrane-bound vacuole termed an "inclusion". Acquisition of lipids (and other nutrients) from the host cell is a critical step in chlamydial replication. Lipid droplets (LD) are ubiquitous, ER-derived neutral lipid-rich storage organelles surrounded by a phospholipids monolayer and associated proteins. Previous studies have shown that LDs accumulate at the periphery of, and eventually translocate into, the chlamydial inclusion. These observations point out to Chlamydia-mediated manipulation of LDs in infected cells, which may impact the function and thereby the protein composition of these organelles. By means of a label-free quantitative mass spectrometry approach we found that the LD proteome is modified in the context of C. trachomatis infection. We determined that LDs isolated from C. trachomatis-infected cells were enriched in proteins related to lipid metabolism, biosynthesis and LD-specific functions. Interestingly, consistent with the observation that LDs intimately associate with the inclusion, a subset of inclusion membrane proteins co-purified with LD protein extracts. Finally, genetic ablation of LDs negatively affected generation of C. trachomatis infectious progeny, consistent with a role for LD biogenesis in optimal chlamydial growth.