Discovery and Expansion of Gene Modules by Seeking Isolated Groups in a Random Graph Process

BACKGROUND: A central problem in systems biology research is the identification and extension of biological modules-groups of genes or proteins participating in a common cellular process or physical complex. As a result, there is a persistent need for practical, principled methods to infer the modular organization of genes from genome-scale data. RESULTS: We introduce a novel approach for the identification of modules based on the persistence of isolated gene groups within an evolving graph process. First, the underlying genomic data is summarized in the form of ranked gene-gene relationships, thereby accommodating studies that quantify the relevant biological relationship directly or indirectly. Then, the observed gene-gene relationship ranks are viewed as the outcome of a random graph process and candidate modules are given by the identifiable subgraphs that arise during this process. An isolation index is computed for each module, which quantifies the statistical significance of its survival time. CONCLUSIONS: The Miso (module isolation) method predicts gene modules from genomic data and the associated isolation index provides a module-specific measure of confidence. Improving on existing alternative, such as graph clustering and the global pruning of dendrograms, this index offers two intuitively appealing features: (1) the score is module-specific; and (2) different choices of threshold correlate logically with the resulting performance, i.e. a stringent cutoff yields high quality predictions, but low sensitivity. Through the analysis of yeast phenotype data, the Miso method is shown to outperform existing alternatives, in terms of the specificity and sensitivity of its predictions

Transcranial Low-Level Laser Therapy Improves Neurological Performance in Traumatic Brain Injury in Mice: Effect of Treatment Repetition Regimen

Low-level laser (light) therapy (LLLT) has been clinically applied around the world for a spectrum of disorders requiring healing, regeneration and prevention of tissue death. One area that is attracting growing interest in this scope is the use of transcranial LLLT to treat stroke and traumatic brain injury (TBI). We developed a mouse model of severe TBI induced by controlled cortical impact and explored the effect of different treatment schedules. Adult male BALB/c mice were divided into 3 broad groups (a) sham-TBI sham-treatment, (b) real-TBI sham-treatment, and (c) real-TBI active-treatment. Mice received active-treatment (transcranial LLLT by continuous wave 810 nm laser, 25 mW/cm[superscript 2], 18 J/cm[superscript 2], spot diameter 1 cm) while sham-treatment was immobilization only, delivered either as a single treatment at 4 hours post TBI, as 3 daily treatments commencing at 4 hours post TBI or as 14 daily treatments. Mice were sacrificed at 0, 4, 7, 14 and 28 days post-TBI for histology or histomorphometry, and injected with bromodeoxyuridine (BrdU) at days 21–27 to allow identification of proliferating cells. Mice with severe TBI treated with 1-laser Tx (and to a greater extent 3-laser Tx) had significant improvements in neurological severity score (NSS), and wire-grip and motion test (WGMT). However 14-laser Tx provided no benefit over TBI-sham control. Mice receiving 1- and 3-laser Tx had smaller lesion size at 28-days (although the size increased over 4 weeks in all TBI-groups) and less Fluoro-Jade staining for degenerating neurons (at 14 days) than in TBI control and 14-laser Tx groups. There were more BrdU-positive cells in the lesion in 1- and 3-laser groups suggesting LLLT may increase neurogenesis. Transcranial NIR laser may provide benefit in cases of acute TBI provided the optimum treatment regimen is employed.National Institutes of Health (U.S.) (Grant R01AI050875)Center for Integration of Medicine and Innovative Technology (DAMD17-02-2-0006)United States. Dept. of Defense. Congressionally Directed Medical Research Programs (W81XWH-09-1-0514)United States. Air Force Office of Scientific Research. Military Photomedicine Program (FA9550-11-1-0331

Dual-Level Regulation of ACC Synthase Activity by MPK3/MPK6 Cascade and Its Downstream WRKY Transcription Factor during Ethylene Induction in Arabidopsis

Plants under pathogen attack produce high levels of ethylene, which plays important roles in plant immunity. Previously, we reported the involvement of ACS2 and ACS6, two Type I ACS isoforms, in Botrytis cinerea–induced ethylene biosynthesis and their regulation at the protein stability level by MPK3 and MPK6, two Arabidopsis pathogen-responsive mitogen-activated protein kinases (MAPKs). The residual ethylene induction in the acs2/acs6 double mutant suggests the involvement of additional ACS isoforms. It is also known that a subset of ACS genes, including ACS6, is transcriptionally induced in plants under stress or pathogen attack. However, the importance of ACS gene activation and the regulatory mechanism(s) are not clear. In this report, we demonstrate using genetic analysis that ACS7 and ACS11, two Type III ACS isoforms, and ACS8, a Type II ACS isoform, also contribute to the B. cinerea–induced ethylene production. In addition to post-translational regulation, transcriptional activation of the ACS genes also plays a critical role in sustaining high levels of ethylene induction. Interestingly, MPK3 and MPK6 not only control the stability of ACS2 and ACS6 proteins via direct protein phosphorylation but also regulate the expression of ACS2 and ACS6 genes. WRKY33, another MPK3/MPK6 substrate, is involved in the MPK3/MPK6-induced ACS2/ACS6 gene expression based on genetic analyses. Furthermore, chromatin-immunoprecipitation assay reveals the direct binding of WRKY33 to the W-boxes in the promoters of ACS2 and ACS6 genes in vivo, suggesting that WRKY33 is directly involved in the activation of ACS2 and ACS6 expression downstream of MPK3/MPK6 cascade in response to pathogen invasion. Regulation of ACS activity by MPK3/MPK6 at both transcriptional and protein stability levels plays a key role in determining the kinetics and magnitude of ethylene induction

EMT and Stem Cell-Like Properties Associated with HIF-2α Are Involved in Arsenite-Induced Transformation of Human Bronchial Epithelial Cells

Arsenic is well-established as a human carcinogen, but the molecular mechanisms leading to arsenic-induced carcinogenesis are complex and elusive. It is not been determined if the epithelial-mesenchymal transition (EMT) and stem cell-like properties contribute in causing to carcinogen-induced malignant transformation and subsequent tumor formation.To investigate the molecular mechanisms underlying EMT and the emergence of cancer stem cell-like properties during neoplastic transformation of human bronchial epithelial (HBE) cells induced by chronic exposure to arsenite. HBE cells were continuously exposed to arsenite. Spheroid formation assays and analyses of side populations (SPs) were performed to confirm that arsenite induces the acquired EMT and cancer stem cell-like phenotype. Treated HBE cells were molecularly characterized by RT-PCR, Western blots, immunofluorescence, Southwestern assays, reporter assays, and chromatin immunoprecipitation.With chronic exposure to arsenite, HBE cells undergo an EMT and then acquire a malignant cancer stem cell-like phenotype. Twist1 and Bmi1 are involved in arsenite-induced EMT. The process is directly regulated by HIF-2α. The self-renewal genes, Oct4, Bmi1, and ALDH1, are necessary for arsenite-mediated maintenance of stem cells.EMT, regulated by HIF-2α, and the development of a cancer stem cell-like phenotype are associated with arsenite-induced transformation of HBE cells

Gene discovery in EST sequences from the wheat leaf rust fungus Puccinia triticina sexual spores, asexual spores and haustoria, compared to other rust and corn smut fungi

© 2011 Xu et al.; licensee BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.DOI: 10.1186/1471-2164-12-161Background.Rust fungi are biotrophic basidiomycete plant pathogens that cause major diseases on plants and trees world-wide, affecting agriculture and forestry. Their biotrophic nature precludes many established molecular genetic manipulations and lines of research. The generation of genomic resources for these microbes is leading to novel insights into biology such as interactions with the hosts and guiding directions for breakthrough research in plant pathology. Results. To support gene discovery and gene model verification in the genome of the wheat leaf rust fungus, Puccinia triticina (Pt), we have generated Expressed Sequence Tags (ESTs) by sampling several life cycle stages. We focused on several spore stages and isolated haustorial structures from infected wheat, generating 17,684 ESTs. We produced sequences from both the sexual (pycniospores, aeciospores and teliospores) and asexual (germinated urediniospores) stages of the life cycle. From pycniospores and aeciospores, produced by infecting the alternate host, meadow rue (Thalictrum speciosissimum), 4,869 and 1,292 reads were generated, respectively. We generated 3,703 ESTs from teliospores produced on the senescent primary wheat host. Finally, we generated 6,817 reads from haustoria isolated from infected wheat as well as 1,003 sequences from germinated urediniospores. Along with 25,558 previously generated ESTs, we compiled a database of 13,328 non-redundant sequences (4,506 singlets and 8,822 contigs). Fungal genes were predicted using the EST version of the self-training GeneMarkS algorithm. To refine the EST database, we compared EST sequences by BLASTN to a set of 454 pyrosequencing-generated contigs and Sanger BAC-end sequences derived both from the Pt genome, and to ESTs and genome reads from wheat. A collection of 6,308 fungal genes was identified and compared to sequences of the cereal rusts, Puccinia graminis f. sp. tritici (Pgt) and stripe rust, P. striiformis f. sp. tritici (Pst), and poplar leaf rust Melampsora species, and the corn smut fungus, Ustilago maydis (Um). While extensive homologies were found, many genes appeared novel and species-specific; over 40% of genes did not match any known sequence in existing databases. Focusing on spore stages, direct comparison to Um identified potential functional homologs, possibly allowing heterologous functional analysis in that model fungus. Many potentially secreted protein genes were identified by similarity searches against genes and proteins of Pgt and Melampsora spp., revealing apparent orthologs. Conclusions. The current set of Pt unigenes contributes to gene discovery in this major cereal pathogen and will be invaluable for gene model verification in the genome sequence

Context Differences Reveal Insulator and Activator Functions of a Su(Hw) Binding Region

Insulators are DNA elements that divide chromosomes into independent transcriptional domains. The Drosophila genome contains hundreds of binding sites for the Suppressor of Hairy-wing [Su(Hw)] insulator protein, corresponding to locations of the retroviral gypsy insulator and non-gypsy binding regions (BRs). The first non-gypsy BR identified, 1A-2, resides in cytological region 1A. Using a quantitative transgene system, we show that 1A-2 is a composite insulator containing enhancer blocking and facilitator elements. We discovered that 1A-2 separates the yellow (y) gene from a previously unannotated, non-coding RNA gene, named yar for y-achaete (ac) intergenic RNA. The role of 1A-2 was elucidated using homologous recombination to excise these sequences from the natural location, representing the first deletion of any Su(Hw) BR in the genome. Loss of 1A-2 reduced yar RNA accumulation, without affecting mRNA levels from the neighboring y and ac genes. These data indicate that within the 1A region, 1A-2 acts an activator of yar transcription. Taken together, these studies reveal that the properties of 1A-2 are context-dependent, as this element has both insulator and enhancer activities. These findings imply that the function of non-gypsy Su(Hw) BRs depends on the genomic environment, predicting that Su(Hw) BRs represent a diverse collection of genomic regulatory elements

Assessment of gene-by-sex interaction effect on bone mineral density

To access publisher's full text version of this article. Please click on the hyperlink in Additional Links field.Sexual dimorphism in various bone phenotypes, including bone mineral density (BMD), is widely observed; however, the extent to which genes explain these sex differences is unclear. To identify variants with different effects by sex, we examined gene-by-sex autosomal interactions genome-wide, and performed expression quantitative trait loci (eQTL) analysis and bioinformatics network analysis. We conducted an autosomal genome-wide meta-analysis of gene-by-sex interaction on lumbar spine (LS) and femoral neck (FN) BMD in 25,353 individuals from 8 cohorts. In a second stage, we followed up the 12 top single-nucleotide polymorphisms (SNPs; p < 1 × 10(-5) ) in an additional set of 24,763 individuals. Gene-by-sex interaction and sex-specific effects were examined in these 12 SNPs. We detected one novel genome-wide significant interaction associated with LS-BMD at the Chr3p26.1-p25.1 locus, near the GRM7 gene (male effect = 0.02 and p = 3.0 × 10(-5) ; female effect = -0.007 and p = 3.3 × 10(-2) ), and 11 suggestive loci associated with either FN- or LS-BMD in discovery cohorts. However, there was no evidence for genome-wide significant (p < 5 × 10(-8) ) gene-by-sex interaction in the joint analysis of discovery and replication cohorts. Despite the large collaborative effort, no genome-wide significant evidence for gene-by-sex interaction was found to influence BMD variation in this screen of autosomal markers. If they exist, gene-by-sex interactions for BMD probably have weak effects, accounting for less than 0.08% of the variation in these traits per implicated SNP. © 2012 American Society for Bone and Mineral Research.Medtronic NIH R01 AG18728 R01HL088119 R01AR046838 U01 HL084756 R01 AR43351 P01-HL45522 R01-MH-078111 R01-MH-083824 Nutrition and Obesity Research Center of Maryland P30DK072488 NIAMS/NIH F32AR059469 Instituto de Salud Carlos III-FIS (Spanish Health Ministry) PI 06/0034 PI08/0183 Canadian Institutes of Health Research (CIHR) NHLBI HHSN268201200036C N01-HC-85239 N01-HC-85079 N01-HC-85086 N01-HC-35129 N01 HC15103 N01 HC-55222 N01-HC-75150 N01-HC-45133 HL080295 HL087652 HL105756 NIA AG-023629 AG-15928 AG-20098 AG-027058 N01AG62101 N01AG62103 N01AG62106 1R01AG032098-01A1 National Center of Advancing Translational Technologies CTSI UL1TR000124 National Institute of Diabetes and Digestive and Kidney Diseases DK063491 EUROSPAN (European Special Populations Research Network) European Commission FP6 STRP grant 018947 LSHG-CT-2006-01947 Netherlands Organisation for Scientific Research Erasmus MC Centre for Medical Systems Biology (CMSB) Netherlands Brain Foundation (HersenStichting Nederland) US National Institute for Arthritis, Musculoskeletal and Skin Diseases National Institute on Aging R01 AR/AG41398 R01 AR050066 R21 AR056405 National Heart, Lung, and Blood Institute's Framingham Heart Study N01-HC-25195 Affymetrix, Inc. N02-HL-6-4278 Canadian Institutes of Health Research from Institute of Aging 165446 Institute of Genetics 179433 Institute of Musculoskeletal health 221765 Intramural Research Program of the NIH, National Institute on Aging National Institutes of Health HHSN268200782096C Hong Kong Research Grant Council HKU 768610M Bone Health Fund of HKU Foundation KC Wong Education Foundation Small Project Funding 201007176237 Matching Grant CRCG Grant Osteoporosis and Endocrine Research Fund Genomics Strategic Research Theme of The University of Hong Kong Netherlands Organisation of Scientific Research NWO Investments 175.010.2005.011 911-03-012 Research Institute for Diseases in the Elderly 014-93-015 Netherlands Genomics Initiative (NGI)/Netherlands Consortium for Healthy Aging (NCHA) 050-060-810 Erasmus Medical Center and Erasmus University, Rotterdam Netherlands Organization for the Health Research and Development (ZonMw) Research Institute for Diseases in the Elderly (RIDE) Ministry of Education, Culture and Science Ministry for Health, Welfare and Sports European Commission (DG XII) Municipality of Rotterdam German Bundesministerium fur Forschung und Technology 01 AK 803 A-H 01 IG 07015

Genetic determinants of heel bone properties: genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n = 14 260), velocity of sound (VOS; n = 15 514) and BMD (n = 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n = 11 452) and new genotyping in 15 cohorts (de novo n = 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 × 10(-8)) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 × 10(-14)). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 × 10(-6) also had the expected direction of association with any fracture (P < 0.05), including three SNPs with P < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, this GWA study reveals the effect of several genes common to central DXA-derived BMD and heel ultrasound/DXA measures and points to a new genetic locus with potential implications for better understanding of osteoporosis pathophysiology