Prediction of alternatively skipped exons and splicing enhancers from exon junction arrays

<p>Abstract</p> <p>Background</p> <p>Alternative splicing of exons in a pre-mRNA transcript is an important mechanism which contributes to protein diversity in human. Arrays for detecting alternative splicing are available using several different probe designs, including those based on exon-junctions. In this work, we introduce a new method for predicting alternatively skipped exons from exon-junction arrays. Predictions based on our method are compared against controls and their sequences are analyzed to identify motifs important for regulating alternative splicing.</p> <p>Results</p> <p>Our comparison of several alternative methods shows that an exon-skipping score based on neighboring junctions best discriminates between positive and negative controls. Sequence analysis of our predicted exons confirms the presence of known splicing regulatory sequences. In addition, we also derive a set of development-related alternatively spliced genes based on fetal versus adult tissue comparisons and find that our predictions are consistent with their functional annotations. <it>Ab initio </it>motif finding algorithms are applied to identify several motifs that may be relevant for splicing during development.</p> <p>Conclusion</p> <p>This work describes a new method for analyzing exon-junction arrays, identifies sequence motifs that are specific for alternative and constitutive splicing and suggests a role for several known splicing factors and their motifs in developmental regulation.</p

Conserved Amino Acid Sequence Features in the α Subunits of MoFe, VFe, and FeFe Nitrogenases

BACKGROUND:This study examines the structural features and phylogeny of the alpha subunits of 69 full-length NifD (MoFe subunit), VnfD (VFe subunit), and AnfD (FeFe subunit) sequences. METHODOLOGY/PRINCIPAL FINDINGS:The analyses of this set of sequences included BLAST scores, multiple sequence alignment, examination of patterns of covariant residues, phylogenetic analysis and comparison of the sequences flanking the conserved Cys and His residues that attach the FeMo cofactor to NifD and that are also conserved in the alternative nitrogenases. The results show that NifD nitrogenases fall into two distinct groups. Group I includes NifD sequences from many genera within Bacteria, including all nitrogen-fixing aerobes examined, as well as strict anaerobes and some facultative anaerobes, but no archaeal sequences. In contrast, Group II NifD sequences were limited to a small number of archaeal and bacterial sequences from strict anaerobes. The VnfD and AnfD sequences fall into two separate groups, more closely related to Group II NifD than to Group I NifD. The pattern of perfectly conserved residues, distributed along the full length of the Group I and II NifD, VnfD, and AnfD, confirms unambiguously that these polypeptides are derived from a common ancestral sequence. CONCLUSIONS/SIGNIFICANCE:There is no indication of a relationship between the patterns of covariant residues specific to each of the four groups discussed above that would give indications of an evolutionary pathway leading from one type of nitrogenase to another. Rather the totality of the data, along with the phylogenetic analysis, is consistent with a radiation of Group I and II NifDs, VnfD and AnfD from a common ancestral sequence. All the data presented here strongly support the suggestion made by some earlier investigators that the nitrogenase family had already evolved in the last common ancestor of the Archaea and Bacteria

Bayesian Inference of Networks Across Multiple Sample Groups and Data Types

In this paper, we develop a graphical modeling framework for the inference of networks across multiple sample groups and data types. In medical studies, this setting arises whenever a set of subjects, which may be heterogeneous due to differing disease stage or subtype, is profiled across multiple platforms, such as metabolomics, proteomics, or transcriptomics data. Our proposed Bayesian hierarchical model first links the network structures within each platform using a Markov random field prior to relate edge selection across sample groups, and then links the network similarity parameters across platforms. This enables joint estimation in a flexible manner, as we make no assumptions on the directionality of influence across the data types or the extent of network similarity across the sample groups and platforms. In addition, our model formulation allows the number of variables and number of subjects to differ across the data types, and only requires that we have data for the same set of groups. We illustrate the proposed approach through both simulation studies and an application to gene expression levels and metabolite abundances on subjects with varying severity levels of Chronic Obstructive Pulmonary Disease (COPD)

Multiple Amino Acid Sequence Alignment Nitrogenase Component 1: Insights into Phylogenetics and Structure-Function Relationships

Amino acid residues critical for a protein's structure-function are retained by natural selection and these residues are identified by the level of variance in co-aligned homologous protein sequences. The relevant residues in the nitrogen fixation Component 1 α- and β-subunits were identified by the alignment of 95 protein sequences. Proteins were included from species encompassing multiple microbial phyla and diverse ecological niches as well as the nitrogen fixation genotypes, anf, nif, and vnf, which encode proteins associated with cofactors differing at one metal site. After adjusting for differences in sequence length, insertions, and deletions, the remaining >85% of the sequence co-aligned the subunits from the three genotypes. Six Groups, designated Anf, Vnf , and Nif I-IV, were assigned based upon genetic origin, sequence adjustments, and conserved residues. Both subunits subdivided into the same groups. Invariant and single variant residues were identified and were defined as “core” for nitrogenase function. Three species in Group Nif-III, Candidatus Desulforudis audaxviator, Desulfotomaculum kuznetsovii, and Thermodesulfatator indicus, were found to have a seleno-cysteine that replaces one cysteinyl ligand of the 8Fe:7S, P-cluster. Subsets of invariant residues, limited to individual groups, were identified; these unique residues help identify the gene of origin (anf, nif, or vnf) yet should not be considered diagnostic of the metal content of associated cofactors. Fourteen of the 19 residues that compose the cofactor pocket are invariant or single variant; the other five residues are highly variable but do not correlate with the putative metal content of the cofactor. The variable residues are clustered on one side of the cofactor, away from other functional centers in the three dimensional structure. Many of the invariant and single variant residues were not previously recognized as potentially critical and their identification provides the bases for new analyses of the three-dimensional structure and for mutagenesis studies

Model based heritability scores for high-throughput sequencing data

Supplementary materials. (PDF 1370 KB

Longitudinal changes in DNA methylation during the onset of islet autoimmunity differentiate between reversion versus progression of islet autoimmunity

BackgroundType 1 diabetes (T1D) is preceded by a heterogenous pre-clinical phase, islet autoimmunity (IA). We aimed to identify pre vs. post-IA seroconversion (SV) changes in DNAm that differed across three IA progression phenotypes, those who lose autoantibodies (reverters), progress to clinical T1D (progressors), or maintain autoantibody levels (maintainers).MethodsThis epigenome-wide association study (EWAS) included longitudinal DNAm measurements in blood (Illumina 450K and EPIC) from participants in Diabetes Autoimmunity Study in the Young (DAISY) who developed IA, one or more islet autoantibodies on at least two consecutive visits. We compared reverters - individuals who sero-reverted, negative for all autoantibodies on at least two consecutive visits and did not develop T1D (n=41); maintainers - continued to test positive for autoantibodies but did not develop T1D (n=60); progressors - developed clinical T1D (n=42). DNAm data were measured before (pre-SV visit) and after IA (post-SV visit). Linear mixed models were used to test for differences in pre- vs post-SV changes in DNAm across the three groups. Linear mixed models were also used to test for group differences in average DNAm. Cell proportions, age, and sex were adjusted for in all models. Median follow-up across all participants was 15.5 yrs. (interquartile range (IQR): 10.8-18.7).ResultsThe median age at the pre-SV visit was 2.2 yrs. (IQR: 0.8-5.3) in progressors, compared to 6.0 yrs. (IQR: 1.3-8.4) in reverters, and 5.7 yrs. (IQR: 1.4-9.7) in maintainers. Median time between the visits was similar in reverters 1.4 yrs. (IQR: 1-1.9), maintainers 1.3 yrs. (IQR: 1.0-2.0), and progressors 1.8 yrs. (IQR: 1.0-2.0). Changes in DNAm, pre- vs post-SV, differed across the groups at one site (cg16066195) and 11 regions. Average DNAm (mean of pre- and post-SV) differed across 22 regions.ConclusionDifferentially changing DNAm regions were located in genomic areas related to beta cell function, immune cell differentiation, and immune cell function

Differentially methylated regions interrogated for metastable epialleles associate with offspring adiposity

Aim: Assess if cord blood differentially methylated regions (DMRs) representing human metastable epialleles (MEs) associate with offspring adiposity in 588 maternal-infant dyads from the Colorado Health Start Study.Materials &amp; methods: DNA methylation was assessed via the Illumina 450K array (~439,500 CpG sites). Offspring adiposity was obtained via air displacement plethysmography. Linear regression modeled the association of DMRs potentially representing MEs with adiposity.Results &amp; conclusion: We identified two potential MEs, ZFP57, which associated with infant adiposity change and B4GALNT4, which associated with infancy and childhood adiposity change. Nine DMRs annotating to genes that annotated to MEs associated with change in offspring adiposity (false discovery rate &lt;0.05). Methylation of approximately 80% of DMRs identified associated with decreased change in adiposity

Joint MiRNA/mRNA expression profiling reveals rhanges consistent with development of dysfunctional corpus luteum after weight gain

<div><p>Obese women exhibit decreased fertility, high miscarriage rates and dysfunctional corpus luteum (CL), but molecular mechanisms are poorly defined. We hypothesized that weight gain induces alterations in CL gene expression. RNA sequencing was used to identify changes in the CL transcriptome in the vervet monkey (<i>Chlorocebus aethiops</i>) during weight gain. 10 months of high-fat, high-fructose diet (HFHF) resulted in a 20% weight gain for HFHF animals vs. 2% for controls (p = 0.03) and a 66% increase in percent fat mass for HFHF group. Ovulation was confirmed at baseline and after intervention in all animals. CL were collected on luteal day 7–9 based on follicular phase estradiol peak. 432 mRNAs and 9 miRNAs were differentially expressed in response to HFHF diet. Specifically, miR-28, miR-26, and let-7b previously shown to inhibit sex steroid production in human granulosa cells, were up-regulated. Using integrated miRNA and gene expression analysis, we demonstrated changes in 52 coordinately regulated mRNA targets corresponding to opposite changes in miRNA. Specifically, 2 targets of miR-28 and 10 targets of miR-26 were down-regulated, including genes linked to follicular development, steroidogenesis, granulosa cell proliferation and survival. To the best of our knowledge, this is the first report of dietary-induced responses of the ovulating ovary to developing adiposity. The observed HFHF diet-induced changes were consistent with development of a dysfunctional CL and provide new mechanistic insights for decreased sex steroid production characteristic of obese women. MiRNAs may represent novel biomarkers of obesity-related subfertility and potential new avenues for therapeutic intervention.</p></div

Accelerated epigenetic age at birth and child emotional and behavioural development in early childhood: a meta-analysis of four prospective cohort studies in ECHO

Background: ‘Epigenetic clocks’ have been developed to accurately predict chronologic gestational age and have been associated with child health outcomes in prior work. Methods: We meta-analysed results from four prospective U.S cohorts investigating the association between epigenetic age acceleration estimated using blood DNA methylation collected at birth and preschool age Childhood Behavior Checklist (CBCL) scores. Results: Epigenetic ageing was not significantly associated with CBCL total problem scores (β = 0.33, 95% CI: −0.95, 0.28) and DSM-oriented pervasive development problem scores (β = −0.23, 95% CI: −0.61, 0.15). No associations were observed for other DSM-oriented subscales. Conclusions: The meta-analysis results suggest that epigenetic gestational age acceleration is not associated with child emotional and behavioural functioning for preschool age group. These findings may relate to our study population, which includes two cohorts enriched for ASD and one preterm birth cohort.; future work should address the role of epigenetic age in child health in other study populations

The association of plasma biomarkers with computed tomography-assessed emphysema phenotypes

Rationale: Chronic obstructive pulmonary disease (COPD) is a phenotypically heterogeneous disease. In COPD, the presence of emphysema is associated with increased mortality and risk of lung cancer. High resolution computed tomography (HRCT) scans are useful in quantifying emphysema but are associated with radiation exposure and high incidence of false positive findings (i.e., nodules). Using a comprehensive biomarker panel, we sought to determine if there was a peripheral blood biomarker signature of emphysema. Methods: 114 plasma biomarkers were measured using a custom assay in 588 individuals enrolled in the COPDGene study. Quantitative emphysema measurements included percent low lung attenuation (%LAA) ≤ −950 HU, ≤ − 910 HU and mean lung attenuation at the 15th percentile on lung attenuation curve (LP15A). Multiple regression analysis was performed to determine plasma biomarkers associated with emphysema independent of covariates age, gender, smoking status, body mass index and FEV1. The findings were subsequently validated using baseline blood samples from a separate cohort of 388 subjects enrolled in the Treatment of Emphysema with a Selective Retinoid Agonist (TESRA) study. Results: Regression analysis identified multiple biomarkers associated with CT-assessed emphysema in COPDGene, including advanced glycosylation end-products receptor (AGER or RAGE, p < 0.001), intercellular adhesion molecule 1 (ICAM, p < 0.001), and chemokine ligand 20 (CCL20, p < 0.001). Validation in the TESRA cohort revealed significant associations with RAGE, ICAM1, and CCL20 with radiologic emphysema (p < 0.001 after meta-analysis). Other biomarkers that were associated with emphysema include CDH1, CDH 13 and SERPINA7, but were not available for validation in the TESRA study. Receiver operating characteristics analysis demonstrated a benefit of adding a biomarker panel to clinical covariates for detecting emphysema, especially in those without severe airflow limitation (AUC 0.85). Conclusions: Our findings, suggest that a panel of blood biomarkers including sRAGE, ICAM1 and CCL20 may serve as a useful surrogate measure of emphysema, and when combined with clinical covariates, may be useful clinically in predicting the presence of emphysema compared to just using covariates alone, especially in those with less severe COPD. Ultimately biomarkers may shed light on disease pathogenesis, providing targets for new treatments. Electronic supplementary material The online version of this article (doi:10.1186/s12931-014-0127-9) contains supplementary material, which is available to authorized users