1
|
Cooper PO, Yang J, Wang HH, Broman MM, Awdalkreem GD, Cresswell GM, Wang L, Goossens E, Lanman NA, Doerge RW, Zheng F, Cheng L, Crist SA, Braun RE, Jerde TJ, Ratliff TL. Inflammation Impacts Androgen Receptor Signaling in Basal Prostate Stem Cells Through Interleukin 1 Receptor Antagonist. RESEARCH SQUARE 2023:rs.3.rs-3539806. [PMID: 38168414 PMCID: PMC10760215 DOI: 10.21203/rs.3.rs-3539806/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
The majority of patients with benign prostate hyperplasia (BPH) exhibit chronic prostate inflammation and the extent of inflammation correlates with the severity of symptoms. How inflammation contributes to prostate enlargement and/or BPH symptoms and the underlying mechanisms are not clearly understood. We established a unique mouse model Prostate Ovalbumin Expressing Transgenic 3 (POET3) that mimics chronic non-bacterial prostatitis in men to study the role of inflammation in prostate hyperplasia. After the injection of ovalbumin peptide-specific T cells, POET3 prostates exhibited an influx of inflammatory cells and an increase in pro-inflammatory cytokines that led to epithelial and stromal hyperplasia. We have previously demonstrated with the POET3 model that inflammation expands the basal prostate stem cell (bPSC) population and promotes bPSC differentiation in organoid cultures. In this study, we investigated the mechanisms underlying the impact of inflammation on bPSC. We found that AR activity was enhanced in inflamed bPSC and was essential for bPSC differentiation in organoid cultures. Most importantly, we identified, for the first time, interleukin 1 receptor antagonist (IL-1RA) as a key regulator of AR in basal stem cells. IL-1RA was one of the top genes upregulated by inflammation and inhibition of IL-1RA abrogated the enhanced AR nuclear accumulation and activity in organoids derived from inflamed bPSC. The mirroring effects of IL-1RA recombinant protein and IL-1α neutralizing antibody suggest that IL-1RA may function by antagonizing IL-1α inhibition of AR expression. Furthermore, we established a lineage tracing model to follow bPSC during inflammation and under castrate conditions. We found that inflammation induced bPSC proliferation and differentiation into luminal cells even under castrate conditions, indicating that AR activation driven by inflammation in bPSC is sufficient for their proliferation and differentiation under androgen-deprived conditions. However, proliferation of the differentiated bPSC in the luminal layer significantly diminished with castration, suggesting inflammation may not maintain AR activity in stromal cells, as stromal cells deprived of androgen after castration could no longer provide paracrine growth factors essential for luminal proliferation. Taken together, we have discovered novel mechanisms through which inflammation modulates AR signaling in bPSC and induces bPSC luminal differentiation that contributes to prostate hyperplasia.
Collapse
|
2
|
Wallin J, Bogdan M, Szulc PA, Doerge RW, Siegmund DO. Ghost QTL and hotspots in experimental crosses: novel approach for modeling polygenic effects. Genetics 2021; 217:6067404. [PMID: 33789342 DOI: 10.1093/genetics/iyaa041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2020] [Accepted: 12/10/2020] [Indexed: 11/14/2022] Open
Abstract
Ghost quantitative trait loci (QTL) are the false discoveries in QTL mapping, that arise due to the "accumulation" of the polygenic effects, uniformly distributed over the genome. The locations on the chromosome that are strongly correlated with the total of the polygenic effects depend on a specific sample correlation structure determined by the genotypes at all loci. The problem is particularly severe when the same genotypes are used to study multiple QTL, e.g. using recombinant inbred lines or studying the expression QTL. In this case, the ghost QTL phenomenon can lead to false hotspots, where multiple QTL show apparent linkage to the same locus. We illustrate the problem using the classic backcross design and suggest that it can be solved by the application of the extended mixed effect model, where the random effects are allowed to have a nonzero mean. We provide formulas for estimating the thresholds for the corresponding t-test statistics and use them in the stepwise selection strategy, which allows for a simultaneous detection of several QTL. Extensive simulation studies illustrate that our approach eliminates ghost QTL/false hotspots, while preserving a high power of true QTL detection.
Collapse
|
3
|
Chen X, Doerge RW. Comments on Dr. Aniket Biswas' Letter to the Editor. Biom J 2020; 62:2034-2035. [PMID: 32776371 DOI: 10.1002/bimj.202000224] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2020] [Revised: 07/13/2020] [Accepted: 07/26/2020] [Indexed: 11/11/2022]
|
4
|
Yao H, Srivastava S, Swyers N, Han F, Doerge RW, Birchler JA. Inbreeding Depression in Genotypically Matched Diploid and Tetraploid Maize. Front Genet 2020; 11:564928. [PMID: 33329701 PMCID: PMC7734256 DOI: 10.3389/fgene.2020.564928] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2020] [Accepted: 10/13/2020] [Indexed: 01/10/2023] Open
Abstract
The genetic and molecular basis of heterosis has long been studied but without a consensus about mechanism. The opposite effect, inbreeding depression, results from repeated self-pollination and leads to a reduction in vigor. A popular explanation for this reaction is the homozygosis of recessive, slightly deleterious alleles upon inbreeding. However, extensive studies in alfalfa indicated that inbreeding between diploids and autotetraploids was similar despite the fact that homozygosis of alleles would be dramatically different. The availability of tetraploid lines of maize generated directly from various inbred lines provided the opportunity to examine this issue in detail in perfectly matched diploid and tetraploid hybrids and their parallel inbreeding regimes. Identical hybrids at the diploid and tetraploid levels were inbred in triplicate for seven generations. At the conclusion of this regime, F1 hybrids and selected representative generations (S1, S3, S5, S7) were characterized phenotypically in randomized blocks during the same field conditions. Quantitative measures of the multiple generations of inbreeding provided little evidence for a distinction in the decline of vigor between the diploids and the tetraploids. The results suggest that the homozygosis of completely recessive, slightly deleterious alleles is an inadequate hypothesis to explain inbreeding depression in general.
Collapse
|
5
|
Chen X, Doerge RW, Sarkar SK. A weighted FDR procedure under discrete and heterogeneous null distributions. Biom J 2020; 62:1544-1563. [PMID: 32367597 DOI: 10.1002/bimj.201900216] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2019] [Revised: 01/21/2020] [Accepted: 01/22/2020] [Indexed: 11/11/2022]
Abstract
Multiple testing (MT) with false discovery rate (FDR) control has been widely conducted in the "discrete paradigm" where p-values have discrete and heterogeneous null distributions. However, in this scenario existing FDR procedures often lose some power and may yield unreliable inference, and for this scenario there does not seem to be an FDR procedure that partitions hypotheses into groups, employs data-adaptive weights and is nonasymptotically conservative. We propose a weighted p-value-based FDR procedure, "weighted FDR (wFDR) procedure" for short, for MT in the discrete paradigm that efficiently adapts to both heterogeneity and discreteness of p-value distributions. We theoretically justify the nonasymptotic conservativeness of the wFDR procedure under independence, and show via simulation studies that, for MT based on p-values of binomial test or Fisher's exact test, it is more powerful than six other procedures. The wFDR procedure is applied to two examples based on discrete data, a drug safety study, and a differential methylation study, where it makes more discoveries than two existing methods.
Collapse
|
6
|
Vickman RE, Yang J, Lanman NA, Cresswell GM, Zheng F, Zhang C, Doerge RW, Crist SA, Mesecar AD, Hu CD, Ratliff TL. Cholesterol Sulfotransferase SULT2B1b Modulates Sensitivity to Death Receptor Ligand TNFα in Castration-Resistant Prostate Cancer. Mol Cancer Res 2019; 17:1253-1263. [PMID: 30824526 PMCID: PMC6548593 DOI: 10.1158/1541-7786.mcr-18-1054] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2018] [Revised: 12/18/2018] [Accepted: 02/25/2019] [Indexed: 12/31/2022]
Abstract
Cholesterol sulfotransferase, SULT2B1b, has been demonstrated to modulate both androgen receptor activity and cell growth properties. However, the mechanism(s) by which SULT2B1b alters these properties within prostate cancer cells has not been described. Furthermore, specific advantages of SULT2B1b expression in prostate cancer cells are not understood. In these studies, single-cell mRNA sequencing was conducted to compare the transcriptomes of SULT2B1b knockdown (KD) versus Control KD LNCaP cells. Over 2,000 differentially expressed genes were identified along with alterations in numerous canonical pathways, including the death receptor signaling pathway. The studies herein demonstrate that SULT2B1b KD increases TNFα expression in prostate cancer cells and results in NF-κB activation in a TNF-dependent manner. More importantly, SULT2B1b KD significantly enhances TNF-mediated apoptosis in both TNF-sensitive LNCaP cells and TNF-resistant C4-2 cells. Overexpression of SULT2B1b in LNCaP cells also decreases sensitivity to TNF-mediated cell death, suggesting that SULT2B1b modulates pathways dictating the TNF sensitivity capacity of prostate cancer cells. Probing human prostate cancer patient datasets further supports this work by providing evidence that SULT2B1b expression is inversely correlated with TNF-related genes, including TNF, CD40LG, FADD, and NFKB1. Together, these data provide evidence that SULT2B1b expression in prostate cancer cells enhances resistance to TNF and may provide a growth advantage. In addition, targeting SULT2B1b may induce an enhanced therapeutic response to TNF treatment in advanced prostate cancer. IMPLICATIONS: These data suggest that SULT2B1b expression enhances resistance to TNF and may promote prostate cancer.
Collapse
|
7
|
Chen X, Doerge RW, Heyse JF. Multiple testing with discrete data: Proportion of true null hypotheses and two adaptive FDR procedures. Biom J 2018; 60:761-779. [PMID: 29748972 DOI: 10.1002/bimj.201700157] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2017] [Revised: 02/24/2018] [Accepted: 02/28/2018] [Indexed: 11/11/2022]
Abstract
We consider multiple testing with false discovery rate (FDR) control when p values have discrete and heterogeneous null distributions. We propose a new estimator of the proportion of true null hypotheses and demonstrate that it is less upwardly biased than Storey's estimator and two other estimators. The new estimator induces two adaptive procedures, that is, an adaptive Benjamini-Hochberg (BH) procedure and an adaptive Benjamini-Hochberg-Heyse (BHH) procedure. We prove that the adaptive BH (aBH) procedure is conservative nonasymptotically. Through simulation studies, we show that these procedures are usually more powerful than their nonadaptive counterparts and that the adaptive BHH procedure is usually more powerful than the aBH procedure and a procedure based on randomized p-value. The adaptive procedures are applied to a study of HIV vaccine efficacy, where they identify more differentially polymorphic positions than the BH procedure at the same FDR level.
Collapse
|
8
|
Cheng R, Doerge RW, Borevitz J. Novel Resampling Improves Statistical Power for Multiple-Trait QTL Mapping. G3 (BETHESDA, MD.) 2017; 7:813-822. [PMID: 28064191 PMCID: PMC5345711 DOI: 10.1534/g3.116.037531] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/16/2016] [Accepted: 12/29/2016] [Indexed: 01/13/2023]
Abstract
Multiple-trait analysis typically employs models that associate a quantitative trait locus (QTL) with all of the traits. As a result, statistical power for QTL detection may not be optimal if the QTL contributes to the phenotypic variation in only a small proportion of the traits. Excluding QTL effects that contribute little to the test statistic can improve statistical power. In this article, we show that an optimal power can be achieved when the number of QTL effects is best estimated, and that a stringent criterion for QTL effect selection may improve power when the number of QTL effects is small but can reduce power otherwise. We investigate strategies for excluding trivial QTL effects, and propose a method that improves statistical power when the number of QTL effects is relatively small, and fairly maintains the power when the number of QTL effects is large. The proposed method first uses resampling techniques to determine the number of nontrivial QTL effects, and then selects QTL effects by the backward elimination procedure for significance test. We also propose a method for testing QTL-trait associations that are desired for biological interpretation in applications. We validate our methods using simulations and Arabidopsis thaliana transcript data.
Collapse
|
9
|
Solhaug EM, Ihinger J, Jost M, Gamboa V, Marchant B, Bradford D, Doerge RW, Tyagi A, Replogle A, Madlung A. Environmental Regulation of Heterosis in the Allopolyploid Arabidopsis suecica. PLANT PHYSIOLOGY 2016; 170:2251-63. [PMID: 26896394 PMCID: PMC4825151 DOI: 10.1104/pp.16.00052] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2016] [Accepted: 02/17/2016] [Indexed: 05/27/2023]
Abstract
Allopolyploids are organisms possessing more than two complete sets of chromosomes from two or more species and are frequently more vigorous than their progenitors. To address the question why allopolyploids display hybrid vigor, we compared the natural allopolyploid Arabidopsis suecica to its progenitor species Arabidopsis thaliana and Arabidopsis arenosa. We measured chlorophyll content, CO2 assimilation, and carbohydrate production under varying light conditions and found that the allopolyploid assimilates more CO2 per unit chlorophyll than either of the two progenitor species in high intensity light. The increased carbon assimilation corresponds with greater starch accumulation, but only in strong light, suggesting that the strength of hybrid vigor is dependent on environmental conditions. In weaker light A. suecica tends to produce as much primary metabolites as the better progenitor. We found that gene expression of LIMIT DEXTRINASE1, a debranching enzyme that cleaves branch points within starch molecules, is at the same level in the allopolyploid as in the maternal progenitor A. thaliana and significantly more expressed than in the paternal progenitor A. arenosa. However, expression differences of β-amylases and GLUCAN-WATER DIKINASE1 were not statistically significantly elevated in the allopolyploid over progenitor expression levels. In contrast to allopolyploids, autopolyploid A. thaliana showed the same photosynthetic rate as diploids, indicating that polyploidization alone is likely not the reason for enhanced vigor in the allopolyploid. Taken together, our data suggest that the magnitude of heterosis in A. suecica is environmentally regulated, arises from more efficient photosynthesis, and, under specific conditions, leads to greater starch accumulation than in its progenitor species.
Collapse
|
10
|
Abstract
By incorporating annotation information into the analysis of next-generation sequencing DNA methylation data, we provide an improvement in performance over current testing procedures. Methylation analysis using genome information (MAGI) is applicable for both unreplicated and replicated data, and provides an effective analysis for studies with low sequencing depth. When compared with current tests, the annotation-informed tests provide an increase in statistical power and offer a significance-based interpretation of differential methylation.
Collapse
|
11
|
Pardo I, Lillemoe HA, Blosser RJ, Choi M, Sauder CAM, Doxey DK, Mathieson T, Hancock BA, Baptiste D, Atale R, Hickenbotham M, Zhu J, Glasscock J, Storniolo AMV, Zheng F, Doerge RW, Liu Y, Badve S, Radovich M, Clare SE. Next-generation transcriptome sequencing of the premenopausal breast epithelium using specimens from a normal human breast tissue bank. Breast Cancer Res 2014; 16:R26. [PMID: 24636070 PMCID: PMC4053088 DOI: 10.1186/bcr3627] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2013] [Accepted: 03/10/2014] [Indexed: 12/12/2022] Open
Abstract
Introduction Our efforts to prevent and treat breast cancer are significantly impeded by a lack of knowledge of the biology and developmental genetics of the normal mammary gland. In order to provide the specimens that will facilitate such an understanding, The Susan G. Komen for the Cure Tissue Bank at the IU Simon Cancer Center (KTB) was established. The KTB is, to our knowledge, the only biorepository in the world prospectively established to collect normal, healthy breast tissue from volunteer donors. As a first initiative toward a molecular understanding of the biology and developmental genetics of the normal mammary gland, the effect of the menstrual cycle and hormonal contraceptives on DNA expression in the normal breast epithelium was examined. Methods Using normal breast tissue from 20 premenopausal donors to KTB, the changes in the mRNA of the normal breast epithelium as a function of phase of the menstrual cycle and hormonal contraception were assayed using next-generation whole transcriptome sequencing (RNA-Seq). Results In total, 255 genes representing 1.4% of all genes were deemed to have statistically significant differential expression between the two phases of the menstrual cycle. The overwhelming majority (221; 87%) of the genes have higher expression during the luteal phase. These data provide important insights into the processes occurring during each phase of the menstrual cycle. There was only a single gene significantly differentially expressed when comparing the epithelium of women using hormonal contraception to those in the luteal phase. Conclusions We have taken advantage of a unique research resource, the KTB, to complete the first-ever next-generation transcriptome sequencing of the epithelial compartment of 20 normal human breast specimens. This work has produced a comprehensive catalog of the differences in the expression of protein-coding genes as a function of the phase of the menstrual cycle. These data constitute the beginning of a reference data set of the normal mammary gland, which can be consulted for comparison with data developed from malignant specimens, or to mine the effects of the hormonal flux that occurs during the menstrual cycle.
Collapse
|
12
|
Lin F, Zhao M, Baumann DD, Ping J, Sun L, Liu Y, Zhang B, Tang Z, Hughes E, Doerge RW, Hughes TJ, Ma J. Molecular response to the pathogen Phytophthora sojae among ten soybean near isogenic lines revealed by comparative transcriptomics. BMC Genomics 2014; 15:18. [PMID: 24410936 PMCID: PMC3893405 DOI: 10.1186/1471-2164-15-18] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2013] [Accepted: 01/07/2014] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Phytophthora root and stem rot (PRR) of soybean, caused by Phytophthora sojae, is controlled by Rps genes. However, little is known regarding the Rps-induced molecular responses to P. sojae and how they actually overlap. We thus sequenced, analyzed, and compared the transcriptomes of 10 near isogenic lines (NILs), each with a unique Rps gene/allele, and the susceptible parent Williams, pre- and post-inoculation with the pathogen. RESULTS A total of 4,330 differentially expressed genes (DEGs) were identified in Williams versus 2,014 to 5,499 DEGs in individual NILs upon inoculation with the pathogen. Comparisons of the DEGs between the NILs and Williams identified incompatible interaction genes (IIGs) and compatible interaction genes (CIGs). Hierarchical cluster and heatmap analyses consistently grouped the NILs into three clusters: Cluster I (Rps1-a), Cluster II (Rps1-b, 1-c and 1-k) and Cluster III (Rps3-a, 3-b, 3-c, 4, 5, and 6), suggesting an overlap in Rps-induced defense signaling among certain NILs. Gene ontology (GO) analysis revealed associations between members of the WRKY family and incompatible reactions and between a number of phytohormone signaling pathways and incompatible/compatible interactions. These associations appear to be distinguished according to the NIL clusters. CONCLUSIONS This study characterized genes and multiple branches of putative regulatory networks associated with resistance to P. sojae in ten soybean NILs, and depicted functional "fingerprints" of individual Rps-mediated resistance responses through comparative transcriptomic analysis. Of particular interest are dramatic variations of detected DEGs, putatively involved in ethylene (ET)-, jasmonic acid (JA)-, (reactive oxygen species) ROS-, and (MAP-kinase) MAPK- signaling, among these soybean NILs, implicating their important roles of these signaling in differentiating molecular defense responses. We hypothesize that different timing and robustness in defense signaling to the same pathogen may be largely responsible for such variations.
Collapse
|
13
|
Sardesai N, Lee LY, Chen H, Yi H, Olbricht GR, Stirnberg A, Jeffries J, Xiong K, Doerge RW, Gelvin SB. Cytokinins secreted by Agrobacterium promote transformation by repressing a plant myb transcription factor. Sci Signal 2013; 6:ra100. [PMID: 24255177 DOI: 10.1126/scisignal.2004518] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
Agrobacterium-mediated transformation is the most widely used technique for generating transgenic plants. However, many crops remain recalcitrant. We found that an Arabidopsis myb family transcription factor (MTF1) inhibited plant transformation susceptibility. Mutating MTF1 increased attachment of several Agrobacterium strains to roots and increased both stable and transient transformation in both susceptible and transformation-resistant Arabidopsis ecotypes. Cytokinins from Agrobacterium tumefaciens decreased the expression of MTF1 through activation of the cytokinin response regulator ARR3. Mutating AHK3 and AHK4, genes that encode cytokinin-responsive kinases, increased the expression of MTF1 and impaired plant transformation. Mutant mtf1 plants also had increased expression of AT14A, which encodes a putative transmembrane receptor for cell adhesion molecules. Plants overexpressing AT14A exhibited increased susceptibility to transformation, whereas at14a mutant plants exhibited decreased attachment of bacteria to roots and decreased transformation, suggesting that AT14A may serve as an anchor point for Agrobacteria. Thus, by promoting bacterial attachment and transformation of resistant plants and increasing such processes in susceptible plants, treating roots with cytokinins may help engineer crops with improved features or yield.
Collapse
|
14
|
Abstract
MOTIVATION The majority of next-generation sequencing technologies effectively sample small amounts of DNA or RNA that are amplified (i.e. copied) before sequencing. The amplification process is not perfect, leading to extreme bias in sequenced read counts. We present a novel procedure to account for amplification bias and demonstrate its effectiveness in mitigating gene length dependence when estimating true gene expression. RESULTS We tested the proposed method on simulated and real data. Simulations indicated that our method captures true gene expression more effectively than classic censoring-based approaches and leads to power gains in differential expression testing, particularly for shorter genes with high transcription rates. We applied our method to an unreplicated Arabidopsis RNA-seq dataset resulting in disparate gene ontologies arising from gene set enrichment analyses. AVAILABILITY AND IMPLEMENTATION R code to perform the RASTA procedures is freely available on the web at www.stat.purdue.edu/∼doerge/.
Collapse
|
15
|
Billingsley CN, Allen JR, Baumann DD, Deitz SL, Blazek JD, Newbauer A, Darrah A, Long BC, Young B, Clement M, Doerge RW, Roper RJ. Non-trisomic homeobox gene expression during craniofacial development in the Ts65Dn mouse model of Down syndrome. Am J Med Genet A 2013; 161A:1866-74. [PMID: 23843306 DOI: 10.1002/ajmg.a.36006] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2012] [Accepted: 04/08/2013] [Indexed: 01/25/2023]
Abstract
Trisomy 21 in humans causes cognitive impairment, craniofacial dysmorphology, and heart defects collectively referred to as Down syndrome. Yet, the pathophysiology of these phenotypes is not well understood. Craniofacial alterations may lead to complications in breathing, eating, and communication. Ts65Dn mice exhibit craniofacial alterations that model Down syndrome including a small mandible. We show that Ts65Dn embryos at 13.5 days gestation (E13.5) have a smaller mandibular precursor but a normal sized tongue as compared to euploid embryos, suggesting a relative instead of actual macroglossia originates during development. Neurological tissues were also altered in E13.5 trisomic embryos. Our array analysis found 155 differentially expressed non-trisomic genes in the trisomic E13.5 mandible, including 20 genes containing a homeobox DNA binding domain. Additionally, Sox9, important in skeletal formation and cell proliferation, was upregulated in Ts65Dn mandible precursors. Our results suggest trisomy causes altered expression of non-trisomic genes in development leading to structural changes associated with DS. Identification of genetic pathways disrupted by trisomy is an important step in proposing rational therapies at relevant time points to ameliorate craniofacial abnormalities in DS and other congenital disorders.
Collapse
|
16
|
Roulin A, Auer PL, Libault M, Schlueter J, Farmer A, May G, Stacey G, Doerge RW, Jackson SA. The fate of duplicated genes in a polyploid plant genome. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2013; 73:143-53. [PMID: 22974547 DOI: 10.1111/tpj.12026] [Citation(s) in RCA: 190] [Impact Index Per Article: 17.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/30/2012] [Revised: 08/09/2012] [Accepted: 09/10/2012] [Indexed: 05/18/2023]
Abstract
Polyploidy is generally not tolerated in animals, but is widespread in plant genomes and may result in extensive genetic redundancy. The fate of duplicated genes is poorly understood, both functionally and evolutionarily. Soybean (Glycine max L.) has undergone two separate polyploidy events (13 and 59 million years ago) that have resulted in 75% of its genes being present in multiple copies. It therefore constitutes a good model to study the impact of whole-genome duplication on gene expression. Using RNA-seq, we tested the functional fate of a set of approximately 18 000 duplicated genes. Across seven tissues tested, approximately 50% of paralogs were differentially expressed and thus had undergone expression sub-functionalization. Based on gene ontology and expression data, our analysis also revealed that only a small proportion of the duplicated genes have been neo-functionalized or non-functionalized. In addition, duplicated genes were often found in collinear blocks, and several blocks of duplicated genes were co-regulated, suggesting some type of epigenetic or positional regulation. We also found that transcription factors and ribosomal protein genes were differentially expressed in many tissues, suggesting that the main consequence of polyploidy in soybean may be at the regulatory level.
Collapse
|
17
|
Olbricht GR, Craig BA, Doerge RW. Incorporating Genomic Annotation into a Hidden Markov Model for DNA Methylation Tiling Array Data. Stat Appl Genet Mol Biol 2012. [DOI: 10.1515/1544-6115.1775] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
18
|
Auer PL, Srivastava S, Doerge RW. Differential expression--the next generation and beyond. Brief Funct Genomics 2011; 11:57-62. [PMID: 22210853 DOI: 10.1093/bfgp/elr041] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
RNA-sequencing (RNA-seq) technologies have not only pushed the boundaries of science, but also pushed the computational and analytic capacities of many laboratories. With respect to mapping and quantifying transcriptomes, RNA-seq has certainly established itself as the approach of choice. However, as the complexities of experiments continue to grow, there is still no standard practice that allows for design, processing, normalization, efficient dimension reduction and/or statistical analysis. With this in mind, we provide a brief review of some of the key challenges that are general to all RNA-seq experiments, namely experimental design, statistical analysis and dimensionality reduction.
Collapse
|
19
|
Long AA, Mahapatra CT, Woodruff EA, Rohrbough J, Leung HT, Shino S, An L, Doerge RW, Metzstein MM, Pak WL, Broadie K. The nonsense-mediated decay pathway maintains synapse architecture and synaptic vesicle cycle efficacy. J Cell Sci 2010; 123:3303-15. [PMID: 20826458 DOI: 10.1242/jcs.069468] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
A systematic Drosophila forward genetic screen for photoreceptor synaptic transmission mutants identified no-on-and-no-off transient C (nonC) based on loss of retinal synaptic responses to light stimulation. The cloned gene encodes phosphatidylinositol-3-kinase-like kinase (PIKK) Smg1, a regulatory kinase of the nonsense-mediated decay (NMD) pathway. The Smg proteins act in an mRNA quality control surveillance mechanism to selectively degrade transcripts containing premature stop codons, thereby preventing the translation of truncated proteins with dominant-negative or deleterious gain-of-function activities. At the neuromuscular junction (NMJ) synapse, an extended allelic series of Smg1 mutants show impaired structural architecture, with decreased terminal arbor size, branching and synaptic bouton number. Functionally, loss of Smg1 results in a ~50% reduction in basal neurotransmission strength, as well as progressive transmission fatigue and greatly impaired synaptic vesicle recycling during high-frequency stimulation. Mutation of other NMD pathways genes (Upf2 and Smg6) similarly impairs neurotransmission and synaptic vesicle cycling. These findings suggest that the NMD pathway acts to regulate proper mRNA translation to safeguard synapse morphology and maintain the efficacy of synaptic function.
Collapse
|
20
|
Stevens JR, Doerge RW. Meta-analysis combines affymetrix microarray results across laboratories. Comp Funct Genomics 2010; 6:116-22. [PMID: 18629222 PMCID: PMC2447518 DOI: 10.1002/cfg.460] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2005] [Accepted: 01/19/2005] [Indexed: 12/31/2022] Open
Abstract
With microarray technology becoming more prevalent in recent years, it is now common for several laboratories to employ the same microarray technology to identify differentially expressed genes that are related to the same phenomenon in the same species. Although experimental specifics may be similar, each laboratory will typically produce a slightly different list of statistically significant genes, which calls into question the validity of each gene list (i.e. which list is best). A statistically-based meta-analytic approach to microarray analysis systematically combines results from the different laboratories to provide a single estimate of the degree of differential expression for each gene. This approach provides a more precise view of genes that are of significant interest, while simultaneously allowing for differences between laboratories. The widely-used Affymetrix oligonucleotide array and its software are of particular interest because the results are naturally suited to a meta-analysis. A simulation model based on the Affymetrix platform is developed to examine the adaptive nature of the meta-analytic approach and to illustrate the utility of such an approach in combining microarray results across laboratories.
Collapse
|
21
|
Pignatta D, Dilkes BP, Yoo SY, Henry IM, Madlung A, Doerge RW, Jeffrey Chen Z, Comai L. Differential sensitivity of the Arabidopsis thaliana transcriptome and enhancers to the effects of genome doubling. THE NEW PHYTOLOGIST 2010; 186:194-206. [PMID: 20409178 DOI: 10.1111/j.1469-8137.2010.03198.x] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
Two fundamental types of polyploids are known: allopolyploids, in which different parental chromosome sets were combined by ancestral hybridization and duplication; and autopolyploids, which derive from multiplication of the same chromosome set. In autopolyploids, changes to the nuclear environment are not as profound as in allopolyploids, and therefore the effects of genome doubling on gene regulation remain unclear. To investigate the consequences of autopolyploidization per se, we performed a microarray analysis in three equivalent lineages of matched diploids and autotetraploids of Arabidopsis thaliana. Additionally, we compared the expression levels of GFP transgenes driven by endogenous enhancer elements (enhancer traps) in diploids and autotetraploid of 16 transgenic lines. We expected that true ploidy-dependent changes should occur in independently derived autopolyploid lineages. By this criterion, our microarray analysis detected few changes associated with polyploidization, while the enhancer-trap analysis revealed altered GFP expression at multiple plant life stages for 25% of the lines tested. Genes on individual traps were coordinately regulated while endogenous gene expression was not affected except for one line. The unique sensitivity of enhancer traps to ploidy, in contrast to the observed stability of genes, could derive from lower complexity of regulatory pathways acting on traps versus endogenous genes.
Collapse
|
22
|
Riddle NC, Jiang H, An L, Doerge RW, Birchler JA. Gene expression analysis at the intersection of ploidy and hybridity in maize. TAG. THEORETICAL AND APPLIED GENETICS. THEORETISCHE UND ANGEWANDTE GENETIK 2010; 120:341-53. [PMID: 19657617 DOI: 10.1007/s00122-009-1113-3] [Citation(s) in RCA: 62] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/05/2009] [Accepted: 07/13/2009] [Indexed: 05/24/2023]
Abstract
Heterosis and polyploidy are two important aspects of plant evolution. To examine these issues, we conducted a global gene expression study of a maize ploidy series as well as a set of tetraploid inbred and hybrid lines. This gene expression analysis complements an earlier phenotypic study of these same materials. We find that ploidy change affects a large fraction of the genome, albeit at low levels; gene expression changes rarely exceed 2-fold and are typically not statistically significant. The most common gene expression profile we detected is greater than linear increase from monoploid to diploid, and reductions from diploid to triploid and from triploid to tetraploid, a trend that mirrors plant stature. When examining heterosis in tetraploid maize lines, we found a large fraction of the genome impacted but the majority of changes were not statistically significant at 2-fold or less. Non-additive expression was common in the hybrids, and the extent of non-additivity increased both in number and magnitude from duplex to quadruplex hybrids. Overall, we find that gene expression trends mirror observations from the phenotypic studies; however, obvious mechanistic connections remain unknown.
Collapse
|
23
|
Tanurdzic M, Vaughn MW, Jiang H, Lee TJ, Slotkin RK, Sosinski B, Thompson WF, Doerge RW, Martienssen RA. Epigenomic consequences of immortalized plant cell suspension culture. PLoS Biol 2009; 6:2880-95. [PMID: 19071958 PMCID: PMC2596858 DOI: 10.1371/journal.pbio.0060302] [Citation(s) in RCA: 169] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2008] [Accepted: 10/23/2008] [Indexed: 11/19/2022] Open
Abstract
Plant cells grown in culture exhibit genetic and epigenetic instability. Using a combination of chromatin immunoprecipitation and DNA methylation profiling on tiling microarrays, we have mapped the location and abundance of histone and DNA modifications in a continuously proliferating, dedifferentiated cell suspension culture of Arabidopsis. We have found that euchromatin becomes hypermethylated in culture and that a small percentage of the hypermethylated genes become associated with heterochromatic marks. In contrast, the heterochromatin undergoes dramatic and very precise DNA hypomethylation with transcriptional activation of specific transposable elements (TEs) in culture. High throughput sequencing of small interfering RNA (siRNA) revealed that TEs activated in culture have increased levels of 21-nucleotide (nt) siRNA, sometimes at the expense of the 24-nt siRNA class. In contrast, TEs that remain silent, which match the predominant 24-nt siRNA class, do not change significantly in their siRNA profiles. These results implicate RNA interference and chromatin modification in epigenetic restructuring of the genome following the activation of TEs in immortalized cell culture.
Collapse
|
24
|
Leung HT, Tseng-Crank J, Kim E, Mahapatra C, Shino S, Zhou Y, An L, Doerge RW, Pak WL. DAG lipase activity is necessary for TRP channel regulation in Drosophila photoreceptors. Neuron 2008; 58:884-96. [PMID: 18579079 DOI: 10.1016/j.neuron.2008.05.001] [Citation(s) in RCA: 86] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2007] [Revised: 04/16/2008] [Accepted: 04/29/2008] [Indexed: 12/31/2022]
Abstract
In Drosophila, a phospholipase C-mediated signaling cascade links photoexcitation of rhodopsin to the opening of the TRP/TRPL channels. A lipid product of the cascade, diacylglycerol (DAG) and its metabolite(s), polyunsaturated fatty acids (PUFAs), have both been proposed as potential excitatory messengers. A crucial enzyme in the understanding of this process is likely to be DAG lipase (DAGL). However, DAGLs that might fulfill this role have not been previously identified in any organism. In this work, the Drosophila DAGL gene, inaE, has been identified from mutants that are defective in photoreceptor responses to light. The inaE-encoded protein isoforms show high sequence similarity to known mammalian DAG lipases, exhibit DAG lipase activity in vitro, and are highly expressed in photoreceptors. Analyses of norpA inaE double mutants and severe inaE mutants show that normal DAGL activity is required for the generation of physiologically meaningful photoreceptor responses.
Collapse
|
25
|
Bogdan M, Frommlet F, Biecek P, Cheng R, Ghosh JK, Doerge RW. Extending the modified bayesian information criterion (mBIC) to dense markers and multiple interval mapping. Biometrics 2008; 64:1162-9. [PMID: 18266892 DOI: 10.1111/j.1541-0420.2008.00989.x] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
SUMMARY The modified version of Bayesian Information Criterion (mBIC) is a relatively simple model selection procedure that can be used when locating multiple interacting quantitative trait loci (QTL). Our earlier work demonstrated the statistical properties of mBIC for situations where the average genetic map interval is at least 5 cM. In this work mBIC is adapted to genome searches based on a dense map and, more importantly, to the situation where consecutive QTL and interactions are located by multiple interval mapping. Easy to use formulas for the extended mBIC are given. A simulation study, as well as the analysis of real data, confirm the good properties of the extended mBIC.
Collapse
|