201
|
Bankier S, Michoel T. eQTLs as causal instruments for the reconstruction of hormone linked gene networks. Front Endocrinol (Lausanne) 2022; 13:949061. [PMID: 36060942 PMCID: PMC9428692 DOI: 10.3389/fendo.2022.949061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/20/2022] [Accepted: 07/25/2022] [Indexed: 11/17/2022] Open
Abstract
Hormones act within in highly dynamic systems and much of the phenotypic response to variation in hormone levels is mediated by changes in gene expression. The increase in the number and power of large genetic association studies has led to the identification of hormone linked genetic variants. However, the biological mechanisms underpinning the majority of these loci are poorly understood. The advent of affordable, high throughput next generation sequencing and readily available transcriptomic databases has shown that many of these genetic variants also associate with variation in gene expression levels as expression Quantitative Trait Loci (eQTLs). In addition to further dissecting complex genetic variation, eQTLs have been applied as tools for causal inference. Many hormone networks are driven by transcription factors, and many of these genes can be linked to eQTLs. In this mini-review, we demonstrate how causal inference and gene networks can be used to describe the impact of hormone linked genetic variation upon the transcriptome within an endocrinology context.
Collapse
|
202
|
Ha D, Kim D, Kim I, Oh Y, Kong J, Han S, Kim S. OUP accepted manuscript. Nucleic Acids Res 2022; 50:1849-1863. [PMID: 35137181 PMCID: PMC8887464 DOI: 10.1093/nar/gkac050] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2021] [Revised: 01/14/2022] [Accepted: 01/25/2022] [Indexed: 11/14/2022] Open
Abstract
Mouse models have been engineered to reveal the biological mechanisms of human diseases based on an assumption. The assumption is that orthologous genes underlie conserved phenotypes across species. However, genetically modified mouse orthologs of human genes do not often recapitulate human disease phenotypes which might be due to the molecular evolution of phenotypic differences across species from the time of the last common ancestor. Here, we systematically investigated the evolutionary divergence of regulatory relationships between transcription factors (TFs) and target genes in functional modules, and found that the rewiring of gene regulatory networks (GRNs) contributes to the phenotypic discrepancies that occur between humans and mice. We confirmed that the rewired regulatory networks of orthologous genes contain a higher proportion of species-specific regulatory elements. Additionally, we verified that the divergence of target gene expression levels, which was triggered by network rewiring, could lead to phenotypic differences. Taken together, a careful consideration of evolutionary divergence in regulatory networks could be a novel strategy to understand the failure or success of mouse models to mimic human diseases. To help interpret mouse phenotypes in human disease studies, we provide quantitative comparisons of gene expression profiles on our website (http://sbi.postech.ac.kr/w/RN).
Collapse
Affiliation(s)
- Doyeon Ha
- Department of Life Sciences, Pohang University of Science and Technology, Pohang, Korea
| | - Donghyo Kim
- Department of Life Sciences, Pohang University of Science and Technology, Pohang, Korea
| | | | - Youngchul Oh
- Department of Life Sciences, Pohang University of Science and Technology, Pohang, Korea
| | - JungHo Kong
- Department of Life Sciences, Pohang University of Science and Technology, Pohang, Korea
| | - Seong Kyu Han
- Department of Life Sciences, Pohang University of Science and Technology, Pohang, Korea
| | - Sanguk Kim
- To whom correspondence should be addressed. Tel: +82 54 279 2348; Fax: +82 54 279 2199;
| |
Collapse
|
203
|
Roscito JG, Sameith K, Kirilenko BM, Hecker N, Winkler S, Dahl A, Rodrigues MT, Hiller M. Convergent and lineage-specific genomic differences in limb regulatory elements in limbless reptile lineages. Cell Rep 2022; 38:110280. [DOI: 10.1016/j.celrep.2021.110280] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2021] [Revised: 11/24/2021] [Accepted: 12/27/2021] [Indexed: 01/02/2023] Open
|
204
|
Yang X, Bu Y, Niu F, Cun Y, Zhang L, Song X. Comprehensive analysis of LIM gene family in wheat reveals the involvement of TaLIM2 in pollen development. PLANT SCIENCE : AN INTERNATIONAL JOURNAL OF EXPERIMENTAL PLANT BIOLOGY 2022; 314:111101. [PMID: 34895538 DOI: 10.1016/j.plantsci.2021.111101] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/03/2021] [Revised: 10/19/2021] [Accepted: 10/21/2021] [Indexed: 06/14/2023]
Abstract
LIM domain proteins were involved in organizing the cytoskeleton, adjusting the metabolism and gene expression, some of them were specific express in pollen. LIM gene family in plants were studied in sunflower, tobacco, foxtail millet, rape, rice and Arabidopsis thaliana, however, it has not been investigated in wheat to date. In the present study, we totally characterized 29 TaLIM genes through genome-wide analysis, which were divided into two categories and five subclasses according to phylogenetic analysis. RNA-Seq analysis indicated the expression patterns of TaLIM genes have specific temporal and spatial characteristics, especially TaLIM2 was highly expressed in fertility anthers. Phenotypic and cytological of BSMV: TaLIM2 showed that it had defects in the later stage of pollen development and germination, which further testified that TaLIM2 was closely related to fertility conversion. These findings will be useful for functional analysis of LIM genes in wheat fertility and contribute to hybrid wheat breeding.
Collapse
Affiliation(s)
- Xuetong Yang
- College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China.
| | - Yaning Bu
- College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China.
| | - Fuqiang Niu
- College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China.
| | - Yujie Cun
- College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China.
| | - Lingli Zhang
- College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China.
| | - Xiyue Song
- College of Agronomy, Northwest A&F University, Yangling, 712100, Shaanxi, China.
| |
Collapse
|
205
|
El Ghamrasni S, Quevedo R, Hawley J, Mazrooei P, Hanna Y, Cirlan I, Zhu H, Bruce JP, Oldfield LE, Yang SYC, Guilhamon P, Reimand J, Cescon DW, Done SJ, Lupien M, Pugh TJ. Mutations in Noncoding Cis-Regulatory Elements Reveal Cancer Driver Cistromes in Luminal Breast Cancer. Mol Cancer Res 2022; 20:102-113. [PMID: 34556523 PMCID: PMC9398156 DOI: 10.1158/1541-7786.mcr-21-0471] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Revised: 07/31/2021] [Accepted: 09/17/2021] [Indexed: 01/07/2023]
Abstract
Whole-genome sequencing of primary breast tumors enabled the identification of cancer driver genes and noncoding cancer driver plexuses from somatic mutations. However, differentiating driver from passenger events among noncoding genetic variants remains a challenge. Herein, we reveal cancer-driver cis-regulatory elements linked to transcription factors previously shown to be involved in development of luminal breast cancers by defining a tumor-enriched catalogue of approximately 100,000 unique cis-regulatory elements from 26 primary luminal estrogen receptor (ER)+ progesterone receptor (PR)+ breast tumors. Integrating this catalog with somatic mutations from 350 publicly available breast tumor whole genomes, we uncovered cancer driver cistromes, defined as the sum of binding sites for a transcription factor, for ten transcription factors in luminal breast cancer such as FOXA1 and ER, nine of which are essential for growth in breast cancer with four exclusive to the luminal subtype. Collectively, we present a strategy to find cancer driver cistromes relying on quantifying the enrichment of noncoding mutations over cis-regulatory elements concatenated into a functional unit. IMPLICATIONS: Mapping the accessible chromatin of luminal breast cancer led to discovery of an accumulation of mutations within cistromes of transcription factors essential to luminal breast cancer. This demonstrates coopting of regulatory networks to drive cancer and provides a framework to derive insight into the noncoding space of cancer.
Collapse
Affiliation(s)
- Samah El Ghamrasni
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
| | - Rene Quevedo
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
| | - James Hawley
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
| | - Parisa Mazrooei
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
- Genentech, South San Francisco, California
| | - Youstina Hanna
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
| | - Iulia Cirlan
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
| | - Helen Zhu
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
- Vector Institute, Toronto, Ontario, Canada
- Ontario Institute for Cancer Research, Toronto, Ontario, Canada
| | - Jeff P Bruce
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
| | - Leslie E Oldfield
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
| | - S Y Cindy Yang
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
| | - Paul Guilhamon
- Developmental and Stem Cell Biology Program, The Hospital for Sick Children, Toronto, Ontario, Canada
- Arthur and Sonia Labatt Brain Tumor Research Centre, The Hospital for Sick Children, Toronto, Ontario, Canada
| | - Jüri Reimand
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
- Ontario Institute for Cancer Research, Toronto, Ontario, Canada
- Department of Molecular Genetics, University of Toronto, Toronto, Ontario, Canada
| | - Dave W Cescon
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
| | - Susan J Done
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
- Department of Laboratory Medicine & Pathobiology, University of Toronto, Toronto, Ontario, Canada
| | - Mathieu Lupien
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada.
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
- Ontario Institute for Cancer Research, Toronto, Ontario, Canada
| | - Trevor J Pugh
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada.
- Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
- Ontario Institute for Cancer Research, Toronto, Ontario, Canada
| |
Collapse
|
206
|
Zhang Y, Chen W, Sang X, Wang T, Gong H, Zhao Y, Zhao P, Wang H. Genome-Wide Identification of the Thaumatin-like Protein Family Genes in Gossypium barbadense and Analysis of Their Responses to Verticillium dahliae Infection. PLANTS (BASEL, SWITZERLAND) 2021; 10:plants10122647. [PMID: 34961118 PMCID: PMC8708996 DOI: 10.3390/plants10122647] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/14/2021] [Revised: 11/05/2021] [Accepted: 11/05/2021] [Indexed: 06/14/2023]
Abstract
(1) Background: Plants respond to pathogen challenge by activating a defense system involving pathogenesis-related (PR) proteins. The PR-5 family includes thaumatin, thaumatin-like proteins (TLPs), and other related proteins. TLPs play an important role in response to biotic and abiotic stresses. Many TLP-encoding genes have been identified and functionally characterized in the model plant species. (2) Results: We identified a total of 90 TLP genes in the G. barbadense genome. They were phylogenetically classified into 10 subfamilies and distributed across 19 chromosomes and nine scaffolds. The genes were characterized by examining their exon-intron structures, promoter cis-elements, conserved domains, synteny and collinearity, gene family evolution, and gene duplications. Several TLP genes were predicted to be targets of miRNAs. Investigation of expression changes of 21 GbTLPs in a G. barbadense cultivar (Hai7124) resistance to Verticillium dahliae revealed 13 GbTLPs being upregulated in response to V. dahliae infection, suggesting a potential role of these GbTLP genes in disease response. (3) Conclusions: The results of this study allow insight into the GbTLP gene family, identify GbTLP genes responsive to V. dahliae infection, and provide candidate genes for future studies of their roles in disease resistance.
Collapse
Affiliation(s)
- Yilin Zhang
- Zhengzhou Research Base, State Key Laboratory of Cotton Biology, School of Agricultural Sciences, Zhengzhou University, Zhengzhou 450001, China; (Y.Z.); (T.W.)
| | - Wei Chen
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Anyang 455000, China; (W.C.); (X.S.); (H.G.)
| | - Xiaohui Sang
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Anyang 455000, China; (W.C.); (X.S.); (H.G.)
| | - Ting Wang
- Zhengzhou Research Base, State Key Laboratory of Cotton Biology, School of Agricultural Sciences, Zhengzhou University, Zhengzhou 450001, China; (Y.Z.); (T.W.)
| | - Haiyan Gong
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Anyang 455000, China; (W.C.); (X.S.); (H.G.)
| | - Yunlei Zhao
- Zhengzhou Research Base, State Key Laboratory of Cotton Biology, School of Agricultural Sciences, Zhengzhou University, Zhengzhou 450001, China; (Y.Z.); (T.W.)
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Anyang 455000, China; (W.C.); (X.S.); (H.G.)
| | - Pei Zhao
- Zhengzhou Research Base, State Key Laboratory of Cotton Biology, School of Agricultural Sciences, Zhengzhou University, Zhengzhou 450001, China; (Y.Z.); (T.W.)
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Anyang 455000, China; (W.C.); (X.S.); (H.G.)
| | - Hongmei Wang
- Zhengzhou Research Base, State Key Laboratory of Cotton Biology, School of Agricultural Sciences, Zhengzhou University, Zhengzhou 450001, China; (Y.Z.); (T.W.)
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Anyang 455000, China; (W.C.); (X.S.); (H.G.)
| |
Collapse
|
207
|
Schweizer G, Wagner A. Both Binding Strength and Evolutionary Accessibility Affect the Population Frequency of Transcription Factor Binding Sequences in Arabidopsis thaliana. Genome Biol Evol 2021; 13:6459646. [PMID: 34894231 PMCID: PMC8712246 DOI: 10.1093/gbe/evab273] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/06/2021] [Indexed: 11/22/2022] Open
Abstract
Mutations in DNA sequences that bind transcription factors and thus modulate gene expression are a source of adaptive variation in gene expression. To understand how transcription factor binding sequences evolve in natural populations of the thale cress Arabidopsis thaliana, we integrated genomic polymorphism data for loci bound by transcription factors with in vitro data on binding affinity for these transcription factors. Specifically, we studied 19 different transcription factors, and the allele frequencies of 8,333 genomic loci bound in vivo by these transcription factors in 1,135 A. thaliana accessions. We find that transcription factor binding sequences show very low genetic diversity, suggesting that they are subject to purifying selection. High frequency alleles of such binding sequences tend to bind transcription factors strongly. Conversely, alleles that are absent from the population tend to bind them weakly. In addition, alleles with high frequencies also tend to be the endpoints of many accessible evolutionary paths leading to these alleles. We show that both high affinity and high evolutionary accessibility contribute to high allele frequency for at least some transcription factors. Although binding sequences with stronger affinity are more frequent, we did not find them to be associated with higher gene expression levels. Epistatic interactions among individual mutations that alter binding affinity are pervasive and can help explain variation in accessibility among binding sequences. In summary, combining in vitro binding affinity data with in vivo binding sequence data can help understand the forces that affect the evolution of transcription factor binding sequences in natural populations.
Collapse
Affiliation(s)
- Gabriel Schweizer
- Department of Evolutionary Biology and Environmental Studies, University of Zürich, Switzerland.,Swiss Institute of Bioinformatics, Quartier Sorge-Batiment Genopode, Lausanne, Switzerland
| | - Andreas Wagner
- Department of Evolutionary Biology and Environmental Studies, University of Zürich, Switzerland.,Swiss Institute of Bioinformatics, Quartier Sorge-Batiment Genopode, Lausanne, Switzerland.,Santa Fe Institute, Santa Fe, New Mexico, USA.,Stellenbosch Institute for Advanced Study (STIAS), Wallenberg Research Centre at Stellenbosch University, South Africa
| |
Collapse
|
208
|
Yang W, Zhao P, Liu Y, Cao P, Ji X, Gao Y, Li P, Cheng J. Transcriptome analysis of lncRNA expression patterns in human congenital lung malformations. BMC Genomics 2021; 22:861. [PMID: 34844556 PMCID: PMC8630864 DOI: 10.1186/s12864-021-08204-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2021] [Accepted: 11/18/2021] [Indexed: 11/17/2022] Open
Abstract
Objectives To explore the long non-coding RNA (lncRNA) expression pattern of congenital lung malformations on a genome-wide scale and investigate their potential biological function in four subtypes of congenital lung malformations. Methods We obtained both lesions and normal lung control tissues from the patients diagnosed with CPAM-I, CPAM-II, ILS, and ILS-CPAM, and underwent lobectomy (i.e., surgical removal of the whole lobe which contains the localized lesion as well as normal lung tissue). Then, we performed lncRNA transcriptome profiling in these tissues by RNA sequencing (RNA-seq). A comprehensive bioinformatics analysis was conducted to characterize the expression profiles and relevant biological functions and for multiple comparisons of lncRNA expression in the different subtypes of congenital lung malformation tissues. Furthermore, the lncRNA-mRNA co-expression network was constructed, and dysregulated mRNAs were functionally analyzed. Finally, gene set enrichment analysis (GSEA) was used to predict the potential molecular mechanism of the identified lncRNAs. Results A total of 5921 lncRNA transcripts were identified between congenital lung malformations tissues and normal lung control tissues. Compared with normal lung control, 481of these expressed lncRNAs were upregulated and 142 were downregulated in CPAM-I, 91 were upregulated and 14 were downregulated in CPAM-II, 39 were upregulated and 38 were downregulated in ILS, and 201 were upregulated and 38 were downregulated in ILS-CPAM. Unsupervised clustering and principal component analysis of the expressed lncRNAs visualized the differences between normal lung control and different subtypes of congenital lung malformations samples. We also confirmed significant differences in the composition of differentially expressed genes (DEGs) and the differentially expressed lncRNAs (DE lncRNAs) between CPAM-I and other subtypes of congenital lung malformations, as well as in normal lung control tissues, and observed enrichment of DEGs in the regulation of the immune system, cell projection organization, and inflammatory pathways. Finally, we identified the lncRNA FLJ26850 might be related to congenital lung malformations via ZNF473. Conclusions Significant differences in lncRNAs expression patterns were observed between different subtypes of congenital lung malformations and normal control. The lncRNA FLJ26850 might be related to congenital lung malformations via ZNF473. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-08204-x.
Collapse
Affiliation(s)
- Weili Yang
- Department of Pediatric Surgery, the Second Affiliated Hospital of Xi'an Jiaotong University, No. 157, Xiwu Road, Xi'an, 710004, Shaanxi, China
| | - Pu Zhao
- Department of Neonatology, the third Affiliated Hospital, Xi'an Jiaotong University, Xi'an, 710068, Shaanxi Province, China
| | - Yun Liu
- Department of Respiratory Medicine, the Second Affiliated Hospital of Xi'an Jiaotong University, Xi'an, 710004, Shaanxi Province, China
| | - Ping Cao
- Department of Pediatric Surgery, the Second Affiliated Hospital of Xi'an Jiaotong University, No. 157, Xiwu Road, Xi'an, 710004, Shaanxi, China
| | - Xiang Ji
- Department of Pediatric Surgery, the Second Affiliated Hospital of Xi'an Jiaotong University, No. 157, Xiwu Road, Xi'an, 710004, Shaanxi, China
| | - Ya Gao
- Department of Pediatric Surgery, the Second Affiliated Hospital of Xi'an Jiaotong University, No. 157, Xiwu Road, Xi'an, 710004, Shaanxi, China
| | - Peng Li
- Department of Pediatric Surgery, the Second Affiliated Hospital of Xi'an Jiaotong University, No. 157, Xiwu Road, Xi'an, 710004, Shaanxi, China.
| | - Jiwen Cheng
- Department of Pediatric Surgery, the Second Affiliated Hospital of Xi'an Jiaotong University, No. 157, Xiwu Road, Xi'an, 710004, Shaanxi, China.
| |
Collapse
|
209
|
Ramírez-Ramírez JA, Madrigal Y, Alzate JF, Pabón-Mora N. Evolution and expression of the MADS-box flowering transition genes AGAMOUS-like 24/SHORT VEGETATIVE PHASE with emphasis in selected Neotropical orchids. Cells Dev 2021; 168:203755. [PMID: 34758403 DOI: 10.1016/j.cdev.2021.203755] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2021] [Revised: 10/24/2021] [Accepted: 10/31/2021] [Indexed: 11/18/2022]
Abstract
In angiosperms the reproductive transition results in the transformation of a vegetative apical meristem (SAM) into an inflorescence meristem (IM), capable of forming floral meristems (FM). Two key players in the flowering transition are AGAMOUS-like 24 (AGL24) and SHORT VEGETATIVE PHASE (SVP). They are eudicot MADS-box paralogs performing opposite roles, as AGL24 positively regulates flowering while SVP represses the reproductive transition in Arabidopsis. We confirm that the Arabidopsis functional reference cannot be readily extrapolated to all eudicots as there are additional duplications of AGL24 in early divergent eudicots and core eudicots with significant sequence variation. In addition, we found that in monocots, two additional independent duplication events have resulted in at least three clades of AGL24/SVP homologs, some only found in Orchidaceae. Protein sequence analyses and comparative evolutionary rates point to higher rates of relaxed negative selection in the Core Eudicot AGL24 B and the Orch SVP-like B clades, in eudicots and monocots respectively. On the other hand, expression data points to plesiomorphic pleiotropic roles of AGL24/SVP genes likely similar to SVP core eudicot genes, and the acquisition of new roles as flowering positive regulators in Core Eudicot AGL24 A genes. Our research presents evidence on the diversification and recruitment of AGL24/SVP homologs in flowering transition in orchids. Although, broad expression of most copies does not allow to determine if they act as flowering repressors or promoters, the restricted expression of some homologs in the SAM suggests putative roles in maintaining the vegetative phase. If so studying in detail the function of AGL24/SVP homologs in orchids is critical to identify putative flowering repressors in a lineage where other canonical repressors remain elusive.
Collapse
Affiliation(s)
- Jessica A Ramírez-Ramírez
- Facultad de Ciencias Exactas y Naturales, Instituto de Biología, Universidad de Antioquia, Medellín, Colombia.
| | - Yesenia Madrigal
- Facultad de Ciencias Exactas y Naturales, Instituto de Biología, Universidad de Antioquia, Medellín, Colombia.
| | - Juan F Alzate
- Centro Nacional de Secuenciación Genómica, Sede de Investigación Universitaria, Departamento de Microbiología y Parasitología, Facultad de Medicina, Universidad de Antioquia, Medellín, Colombia.
| | - Natalia Pabón-Mora
- Facultad de Ciencias Exactas y Naturales, Instituto de Biología, Universidad de Antioquia, Medellín, Colombia.
| |
Collapse
|
210
|
Contribution of Dysregulated DNA Methylation to Autoimmunity. Int J Mol Sci 2021; 22:ijms222111892. [PMID: 34769338 PMCID: PMC8584328 DOI: 10.3390/ijms222111892] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2021] [Revised: 10/22/2021] [Accepted: 10/29/2021] [Indexed: 12/24/2022] Open
Abstract
Epigenetic mechanisms, such as DNA methylation, histone modifications, and non-coding RNAs are known regulators of gene expression and genomic stability in cell growth, development, and differentiation. Because epigenetic mechanisms can regulate several immune system elements, epigenetic alterations have been found in several autoimmune diseases. The purpose of this review is to discuss the epigenetic modifications, mainly DNA methylation, involved in autoimmune diseases in which T cells play a significant role. For example, Rheumatoid Arthritis and Systemic Lupus Erythematosus display differential gene methylation, mostly hypomethylated 5′-C-phosphate-G-3′ (CpG) sites that may associate with disease activity. However, a clear association between DNA methylation, gene expression, and disease pathogenesis must be demonstrated. A better understanding of the impact of epigenetic modifications on the onset of autoimmunity will contribute to the design of novel therapeutic approaches for these diseases.
Collapse
|
211
|
Massey JH, Li J, Stern DL, Wittkopp PJ. Distinct genetic architectures underlie divergent thorax, leg, and wing pigmentation between Drosophila elegans and D. gunungcola. Heredity (Edinb) 2021; 127:467-474. [PMID: 34537820 PMCID: PMC8551284 DOI: 10.1038/s41437-021-00467-0] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Accepted: 08/09/2021] [Indexed: 02/07/2023] Open
Abstract
Pigmentation divergence between Drosophila species has emerged as a model trait for studying the genetic basis of phenotypic evolution, with genetic changes contributing to pigmentation differences often mapping to genes in the pigment synthesis pathway and their regulators. These studies of Drosophila pigmentation have tended to focus on pigmentation changes in one body part for a particular pair of species, but changes in pigmentation are often observed in multiple body parts between the same pair of species. The similarities and differences of genetic changes responsible for divergent pigmentation in different body parts of the same species thus remain largely unknown. Here we compare the genetic basis of pigmentation divergence between Drosophila elegans and D. gunungcola in the wing, legs, and thorax. Prior work has shown that regions of the genome containing the pigmentation genes yellow and ebony influence the size of divergent male-specific wing spots between these two species. We find that these same two regions of the genome underlie differences in leg and thorax pigmentation; however, divergent alleles in these regions show differences in allelic dominance and epistasis among the three body parts. These complex patterns of inheritance can be explained by a model of evolution involving tissue-specific changes in the expression of Yellow and Ebony between D. elegans and D. gunungcola.
Collapse
Affiliation(s)
- Jonathan H Massey
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, USA
- Janelia Research Campus of the Howard Hughes Medical Institute, Ashburn, VA, USA
| | - Jun Li
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, USA
- Institute of Evolution and Ecology, School of Life Sciences, Central China Normal University, Wuhan, China
| | - David L Stern
- Janelia Research Campus of the Howard Hughes Medical Institute, Ashburn, VA, USA
| | - Patricia J Wittkopp
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, USA.
- Department of Molecular, Cellular, and Developmental Biology, University of Michigan, Ann Arbor, MI, USA.
| |
Collapse
|
212
|
Pu J, Wang Z, Cong H, Chin JSR, Justen J, Finet C, Yew JY, Chung H. Repression precedes independent evolutionary gains of a highly specific gene expression pattern. Cell Rep 2021; 37:109896. [PMID: 34706247 PMCID: PMC8578697 DOI: 10.1016/j.celrep.2021.109896] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2020] [Revised: 08/24/2021] [Accepted: 10/06/2021] [Indexed: 12/12/2022] Open
Abstract
Highly specific expression patterns can be caused by the overlapping activities of activator and repressor sequences in enhancers. However, few studies illuminate how these sequences evolve in the origin of new enhancers. Here, we show that expression of the bond gene in the semicircular wall epithelium (swe) of the Drosophila melanogaster male ejaculatory bulb (EB) is controlled by an enhancer consisting of an activator region that requires Abdominal-B driving expression in the entire EB and a repressor region that restricts this expression to the EB swe. Although this expression pattern is independently gained in the distantly related Scaptodrosophila lebanonensis and does not require Abdominal-B, we show that functionally similar repressor sequences are present in Scaptodrosophila and also in species that do not express bond in the EB. We suggest that during enhancer evolution, repressor sequences can precede the evolution of activator sequences and may lead to similar but independently evolved expression patterns. Pu et al. show that the independent gain of a highly specific expression pattern across distantly related species may be because of the preexistence of repressor sequences that precedes the diversification of these species. This may reflect a general mechanism underlying the evolution of highly specific enhancers.
Collapse
Affiliation(s)
- Jian Pu
- Department of Entomology, Michigan State University, East Lansing, MI 48824, USA.
| | - Zinan Wang
- Department of Entomology, Michigan State University, East Lansing, MI 48824, USA; Ecology, Evolution, and Behavior Program, Michigan State University, East Lansing, MI 48824, USA
| | - Haosu Cong
- Department of Entomology, Michigan State University, East Lansing, MI 48824, USA
| | - Jacqueline S R Chin
- Singapore Institute for Clinical Sciences, Agency for Science, Technology and Research (A(∗)STAR), Brenner Centre for Molecular Medicine, Singapore 117609, Singapore
| | - Jessa Justen
- Laboratory of Cellular and Molecular Biology, University of Wisconsin, Madison, WI 53706, USA
| | - Cédric Finet
- Yale-NUS College, 16 College Avenue West, Singapore 138527, Singapore
| | - Joanne Y Yew
- Pacific Biosciences Research Center, University of Hawai'i at Mānoa, Honolulu, HI 96822, USA
| | - Henry Chung
- Department of Entomology, Michigan State University, East Lansing, MI 48824, USA; Ecology, Evolution, and Behavior Program, Michigan State University, East Lansing, MI 48824, USA.
| |
Collapse
|
213
|
Liu Y, Hua YP, Chen H, Zhou T, Yue CP, Huang JY. Genome-scale identification of plant defensin ( PDF) family genes and molecular characterization of their responses to diverse nutrient stresses in allotetraploid rapeseed. PeerJ 2021; 9:e12007. [PMID: 34603847 PMCID: PMC8445089 DOI: 10.7717/peerj.12007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2021] [Accepted: 07/27/2021] [Indexed: 11/22/2022] Open
Abstract
Plant defensins (PDFs), short peptides with strong antibacterial activity, play important roles in plant growth, development, and stress resistance. However, there are few systematic analyses on PDFs in Brassica napus. Here, bioinformatics methods were used to identify genome-wide PDFs in Brassica napus, and systematically analyze physicochemical properties, expansion pattern, phylogeny, and expression profiling of BnaPDFs under diverse nutrient stresses. A total of 37 full-length PDF homologs, divided into two subgroups (PDF1s and PDF2s), were identified in the rapeseed genome. A total of two distinct clades were identified in the BnaPDF phylogeny. Clade specific conserved motifs were identified within each clade respectively. Most BnaPDFs were proved to undergo powerful purified selection. The PDF members had enriched cis-elements related to growth and development, hormone response, environmental stress response in their promoter regions. GO annotations indicate that the functional pathways of BnaPDFs are mainly involved in cells killing and plant defense responses. In addition, bna-miRNA164 and bna-miRNA172 respectively regulate the expression of their targets BnaA2.PDF2.5 and BnaC7.PDF2.6. The expression patterns of BnaPDFs were analyzed in different tissues. BnaPDF1.2bs was mainly expressed in the roots, whereas BnaPDF2.2s and BnaPDF2.3s were both expressed in stamen, pericarp, silique, and stem. However, the other BnaPDF members showed low expression levels in various tissues. Differential expression of BnaPDFs under nitrate limitation, ammonium excess, phosphorus starvation, potassium deficiency, cadmium toxicity, and salt stress indicated that they might participate in different nutrient stress resistance. The genome-wide identification and characterization of BnaPDFs will enrich understanding of their molecular characteristics and provide elite gene resources for genetic improvement of rapeseed resistance to nutrient stresses.
Collapse
Affiliation(s)
- Ying Liu
- School of Agricultural Sciences, Zhengzhou University, Zhengzhou, China
| | - Ying-Peng Hua
- School of Agricultural Sciences, Zhengzhou University, Zhengzhou, China
| | - Huan Chen
- National Tobacco Quality Supervision and Inspection Center, Zhengzhou, China
| | - Ting Zhou
- School of Agricultural Sciences, Zhengzhou University, Zhengzhou, China
| | - Cai-Peng Yue
- School of Agricultural Sciences, Zhengzhou University, Zhengzhou, China
| | - Jin-Yong Huang
- School of Agricultural Sciences, Zhengzhou University, Zhengzhou, China
| |
Collapse
|
214
|
El-Sappah AH, Elbaiomy RG, Elrys AS, Wang Y, Zhu Y, Huang Q, Yan K, Xianming Z, Abbas M, El-Tarabily KA, Li J. Genome-Wide Identification and Expression Analysis of Metal Tolerance Protein Gene Family in Medicago truncatula Under a Broad Range of Heavy Metal Stress. Front Genet 2021; 12:713224. [PMID: 34603378 PMCID: PMC8482800 DOI: 10.3389/fgene.2021.713224] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2021] [Accepted: 06/30/2021] [Indexed: 01/02/2023] Open
Abstract
Metal tolerance proteins (MTPs) encompass plant membrane divalent cation transporters to specifically participate in heavy metal stress resistance and mineral acquisition. However, the molecular behaviors and biological functions of this family in Medicago truncatula are scarcely known. A total of 12 potential MTP candidate genes in the M. truncatula genome were successfully identified and analyzed for a phylogenetic relationship, chromosomal distributions, gene structures, docking analysis, gene ontology, and previous gene expression. M. truncatula MTPs (MtMTPs) were further classified into three major cation diffusion facilitator (CDFs) groups: Mn-CDFs, Zn-CDFs, and Fe/Zn-CDFs. The structural analysis of MtMTPs displayed high gene similarity within the same group where all of them have cation_efflux domain or ZT_dimer. Cis-acting element analysis suggested that various abiotic stresses and phytohormones could induce the most MtMTP gene transcripts. Among all MTPs, PF16916 is the specific domain, whereas GLY, ILE, LEU, MET, ALA, SER, THR, VAL, ASN, and PHE amino acids were predicted to be the binding residues in the ligand-binding site of all these proteins. RNA-seq and gene ontology analysis revealed the significant role of MTP genes in the growth and development of M. truncatula. MtMTP genes displayed differential responses in plant leaves, stems, and roots under five divalent heavy metals (Cd2+, Co2+, Mn2+, Zn2+, and Fe2+). Ten, seven, and nine MtMTPs responded to at least one metal ion treatment in the leaves, stems, and roots, respectively. Additionally, MtMTP1.1, MtMTP1.2, and MtMTP4 exhibited the highest expression responses in most heavy metal treatments. Our results presented a standpoint on the evolution of MTPs in M. truncatula. Overall, our study provides a novel insight into the evolution of the MTP gene family in M. truncatula and paves the way for additional functional characterization of this gene family.
Collapse
Affiliation(s)
- Ahmed H El-Sappah
- School of Agriculture, Forestry and Food Engineering, Yibin University, Yibin, China.,Genetics Department, Faculty of Agriculture, Zagazig University, Zagazig, Egypt
| | | | - Ahmed S Elrys
- Soil Science Department, Faculty of Agriculture, Zagazig University, Zagazig, Egypt
| | - Yu Wang
- School of Agriculture, Forestry and Food Engineering, Yibin University, Yibin, China
| | - Yumin Zhu
- School of Agriculture, Forestry and Food Engineering, Yibin University, Yibin, China
| | - Qiulan Huang
- College of Tea Science, Yibin University, Yibin, China
| | - Kuan Yan
- School of Agriculture, Forestry and Food Engineering, Yibin University, Yibin, China
| | - Zhao Xianming
- School of Agriculture, Forestry and Food Engineering, Yibin University, Yibin, China
| | - Manzar Abbas
- School of Agriculture, Forestry and Food Engineering, Yibin University, Yibin, China
| | - Khaled A El-Tarabily
- Department of Biology, College of Science, United Arab Emirates University, Al Ain, United Arab Emirates.,Harry Butler Institute, Murdoch University, Murdoch, WA, Australia
| | - Jia Li
- School of Agriculture, Forestry and Food Engineering, Yibin University, Yibin, China
| |
Collapse
|
215
|
Niu H, Xia P, Hu Y, Zhan C, Li Y, Gong S, Li Y, Ma D. Genome-wide identification of ZF-HD gene family in Triticum aestivum: Molecular evolution mechanism and function analysis. PLoS One 2021; 16:e0256579. [PMID: 34559835 PMCID: PMC8462724 DOI: 10.1371/journal.pone.0256579] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2021] [Accepted: 08/11/2021] [Indexed: 12/04/2022] Open
Abstract
ZF-HD family genes play important roles in plant growth and development. Studies about the whole genome analysis of ZF-HD gene family have been reported in some plant species. In this study, the whole genome identification and expression profile of the ZF-HD gene family were analyzed for the first time in wheat. A total of 37 TaZF-HD genes were identified and divided into TaMIF and TaZHD subfamilies according to the conserved domain. The phylogeny tree of the TaZF-HD proteins was further divided into six groups based on the phylogenetic relationship. The 37 TaZF-HDs were distributed on 18 of 21 chromosomes, and almost all the genes had no introns. Gene duplication and Ka/Ks analysis showed that the gene family may have experienced powerful purification selection pressure during wheat evolution. The qRT-PCR analysis showed that TaZF-HD genes had significant expression patterns in different biotic stress and abiotic stress. Through subcellular localization experiments, we found that TaZHD6-3B was located in the nucleus, while TaMIF4-5D was located in the cell membrane and nucleus. Our research contributes to a comprehensive understanding of the TaZF-HD family, provides a new perspective for further research on the biological functions of TaZF-HD genes in wheat.
Collapse
Affiliation(s)
- Hongli Niu
- Hubei Collaborative Innovation Center for Grain Industry/Engineering Research Center of Ecology and Agricultural Use of Wetland, Ministry of Education/College of Agriculture, Yangtze University, Jingzhou, China
| | - Pengliang Xia
- Enshi Tobacco Company of Hubei Province, Enshi, China
| | - Yifeng Hu
- Hubei Collaborative Innovation Center for Grain Industry/Engineering Research Center of Ecology and Agricultural Use of Wetland, Ministry of Education/College of Agriculture, Yangtze University, Jingzhou, China
| | - Chuang Zhan
- Hubei Collaborative Innovation Center for Grain Industry/Engineering Research Center of Ecology and Agricultural Use of Wetland, Ministry of Education/College of Agriculture, Yangtze University, Jingzhou, China
| | - Yiting Li
- Hubei Collaborative Innovation Center for Grain Industry/Engineering Research Center of Ecology and Agricultural Use of Wetland, Ministry of Education/College of Agriculture, Yangtze University, Jingzhou, China
| | - Shuangjun Gong
- Hubei Collaborative Innovation Center for Grain Industry/Engineering Research Center of Ecology and Agricultural Use of Wetland, Ministry of Education/College of Agriculture, Yangtze University, Jingzhou, China
| | - Yan Li
- Hubei Collaborative Innovation Center for Grain Industry/Engineering Research Center of Ecology and Agricultural Use of Wetland, Ministry of Education/College of Agriculture, Yangtze University, Jingzhou, China
- * E-mail: (YL); (DM)
| | - Dongfang Ma
- Hubei Collaborative Innovation Center for Grain Industry/Engineering Research Center of Ecology and Agricultural Use of Wetland, Ministry of Education/College of Agriculture, Yangtze University, Jingzhou, China
- Key Laboratory of Integrated Pest Management on Crop in Central China, Ministry of Agriculture/Hubei Province Key Laboratory for Control of Crop Diseases, Pest and Weeds/Institute of Plant Protection and Soil Science, Hubei Academy of Agricultural Sciences, Wuhan, China
- * E-mail: (YL); (DM)
| |
Collapse
|
216
|
Ullastres A, Merenciano M, González J. Regulatory regions in natural transposable element insertions drive interindividual differences in response to immune challenges in Drosophila. Genome Biol 2021; 22:265. [PMID: 34521452 PMCID: PMC8439047 DOI: 10.1186/s13059-021-02471-3] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Accepted: 08/19/2021] [Indexed: 02/08/2023] Open
Abstract
Background Variation in gene expression underlies interindividual variability in relevant traits including immune response. However, the genetic variation responsible for these gene expression changes remains largely unknown. Among the non-coding variants that could be relevant, transposable element insertions are promising candidates as they have been shown to be a rich and diverse source of cis-regulatory elements. Results In this work, we use a population genetics approach to identify transposable element insertions likely to increase the tolerance of Drosophila melanogaster to bacterial infection by affecting the expression of immune-related genes. We identify 12 insertions associated with allele-specific expression changes in immune-related genes. We experimentally validate three of these insertions including one likely to be acting as a silencer, one as an enhancer, and one with a dual role as enhancer and promoter. The direction in the change of gene expression associated with the presence of several of these insertions is consistent with an increased survival to infection. Indeed, for one of the insertions, we show that this is the case by analyzing both natural populations and CRISPR/Cas9 mutants in which the insertion is deleted from its native genomic context. Conclusions We show that transposable elements contribute to gene expression variation in response to infection in D. melanogaster and that this variation is likely to affect their survival capacity. Because the role of transposable elements as regulatory elements is not restricted to Drosophila, transposable elements are likely to play a role in immune response in other organisms as well. Supplementary Information The online version contains supplementary material available at 10.1186/s13059-021-02471-3.
Collapse
Affiliation(s)
- Anna Ullastres
- Institute of Evolutionary Biology (CSIC-Universitat Pompeu Fabra), Passeig Marítim de la Barceloneta 37-49, 08003, Barcelona, Spain
| | - Miriam Merenciano
- Institute of Evolutionary Biology (CSIC-Universitat Pompeu Fabra), Passeig Marítim de la Barceloneta 37-49, 08003, Barcelona, Spain
| | - Josefa González
- Institute of Evolutionary Biology (CSIC-Universitat Pompeu Fabra), Passeig Marítim de la Barceloneta 37-49, 08003, Barcelona, Spain.
| |
Collapse
|
217
|
Ahmed FF, Hossen MI, Sarkar MAR, Konak JN, Zohra FT, Shoyeb M, Mondal S. Genome-wide identification of DCL, AGO and RDR gene families and their associated functional regulatory elements analyses in banana (Musa acuminata). PLoS One 2021; 16:e0256873. [PMID: 34473743 PMCID: PMC8412350 DOI: 10.1371/journal.pone.0256873] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2021] [Accepted: 08/17/2021] [Indexed: 12/15/2022] Open
Abstract
RNA silencing is mediated through RNA interference (RNAi) pathway gene families, i.e., Dicer-Like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) and their cis-acting regulatory elements. The RNAi pathway is also directly connected with the post-transcriptional gene silencing (PTGS) mechanism, and the pathway controls eukaryotic gene regulation during growth, development, and stress response. Nevertheless, genome-wide identification of RNAi pathway gene families such as DCL, AGO, and RDR and their regulatory network analyses related to transcription factors have not been studied in many fruit crop species, including banana (Musa acuminata). In this study, we studied in silico genome-wide identification and characterization of DCL, AGO, and RDR genes in bananas thoroughly via integrated bioinformatics approaches. A genome-wide analysis identified 3 MaDCL, 13 MaAGO, and 5 MaRDR candidate genes based on multiple sequence alignment and phylogenetic tree related to the RNAi pathway in banana genomes. These genes correspond to the Arabidopsis thaliana RNAi silencing genes. The analysis of the conserved domain, motif, and gene structure (exon-intron numbers) for MaDCL, MaAGO, and MaRDR genes showed higher homogeneity within the same gene family. The Gene Ontology (GO) enrichment analysis exhibited that the identified RNAi genes could be involved in RNA silencing and associated metabolic pathways. A number of important transcription factors (TFs), e.g., ERF, Dof, C2H2, TCP, GATA and MIKC_MADS families, were identified by network and sub-network analyses between TFs and candidate RNAi gene families. Furthermore, the cis-acting regulatory elements related to light-responsive (LR), stress-responsive (SR), hormone-responsive (HR), and other activities (OT) functions were identified in candidate MaDCL, MaAGO, and MaRDR genes. These genome-wide analyses of these RNAi gene families provide valuable information related to RNA silencing, which would shed light on further characterization of RNAi genes, their regulatory elements, and functional roles, which might be helpful for banana improvement in the breeding program.
Collapse
Affiliation(s)
- Fee Faysal Ahmed
- Faculty of Science, Department of Mathematics, Jashore University of Science and Technology, Jashore, Bangladesh
- * E-mail:
| | - Md. Imran Hossen
- Faculty of Science, Department of Mathematics, Jashore University of Science and Technology, Jashore, Bangladesh
| | - Md. Abdur Rauf Sarkar
- Faculty of Biological Science and Technology, Department of Genetic Engineering and Biotechnology, Jashore University of Science and Technology, Jashore, Bangladesh
| | - Jesmin Naher Konak
- Faculty of Life Science, Department of Biochemistry and Molecular Biology, Mawlana Bhashani Science and Technology University, Tangail, Bangladesh
| | - Fatema Tuz Zohra
- Faculty of Agriculture, Laboratory of Fruit Science, Saga University, Honjo-machi, Saga, Japan
| | - Md. Shoyeb
- Faculty of Biological Science and Technology, Department of Genetic Engineering and Biotechnology, Jashore University of Science and Technology, Jashore, Bangladesh
| | - Samiran Mondal
- Faculty of Science, Department of Mathematics, Jashore University of Science and Technology, Jashore, Bangladesh
| |
Collapse
|
218
|
Van Dyke K, Lutz S, Mekonnen G, Myers CL, Albert FW. Trans-acting genetic variation affects the expression of adjacent genes. Genetics 2021; 217:6126816. [PMID: 33789351 DOI: 10.1093/genetics/iyaa051] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2020] [Accepted: 12/16/2020] [Indexed: 11/13/2022] Open
Abstract
Gene expression differences among individuals are shaped by trans-acting expression quantitative trait loci (eQTLs). Most trans-eQTLs map to hotspot locations that influence many genes. The molecular mechanisms perturbed by hotspots are often assumed to involve "vertical" cascades of effects in pathways that can ultimately affect the expression of thousands of genes. Here, we report that trans-eQTLs can affect the expression of adjacent genes via "horizontal" mechanisms that extend along a chromosome. Genes affected by trans-eQTL hotspots in the yeast Saccharomyces cerevisiae were more likely to be located next to each other than expected by chance. These paired hotspot effects tended to occur at adjacent genes that also show coexpression in response to genetic and environmental perturbations, suggesting shared mechanisms. Physical proximity and shared chromatin state, in addition to regulation of adjacent genes by similar transcription factors, were independently associated with paired hotspot effects among adjacent genes. Paired effects of trans-eQTLs can occur at neighboring genes even when these genes do not share a common function. This phenomenon could result in unexpected connections between regulatory genetic variation and phenotypes.
Collapse
Affiliation(s)
- Krisna Van Dyke
- Department of Genetics, Cell Biology, and Development, University of Minnesota, Minneapolis, MN 55455, USA
| | - Sheila Lutz
- Department of Genetics, Cell Biology, and Development, University of Minnesota, Minneapolis, MN 55455, USA
| | - Gemechu Mekonnen
- Department of Genetics, Cell Biology, and Development, University of Minnesota, Minneapolis, MN 55455, USA
| | - Chad L Myers
- Department of Computer Science and Engineering, University of Minnesota, Minneapolis, MN 55455, USA
| | - Frank W Albert
- Department of Genetics, Cell Biology, and Development, University of Minnesota, Minneapolis, MN 55455, USA
| |
Collapse
|
219
|
Abstract
Saccharomyces cerevisiae rewires its transcriptional output to survive stressful environments, such as nitrogen scarcity under fermentative conditions. Although divergence in nitrogen metabolism among natural yeast populations has been reported, the impact of regulatory genetic variants modulating gene expression and nitrogen consumption remains to be investigated. Here, we employed an F1 hybrid from two contrasting S. cerevisiae strains, providing a controlled genetic environment to map cis factors involved in the divergence of gene expression regulation in response to nitrogen scarcity. We used a dual approach to obtain genome-wide allele-specific profiles of chromatin accessibility, transcription factor binding, and gene expression through ATAC-seq (assay for transposase accessible chromatin) and RNA-seq (transcriptome sequencing). We observed large variability in allele-specific expression and accessibility between the two genetic backgrounds, with a third of these differences specific to a deficient nitrogen environment. Furthermore, we discovered events of allelic bias in gene expression correlating with allelic bias in transcription factor binding solely under nitrogen scarcity, where the majority of these transcription factors orchestrates the nitrogen catabolite repression regulatory pathway and demonstrates a cis × environment-specific response. Our approach allowed us to find cis variants modulating gene expression, chromatin accessibility, and allelic differences in transcription factor binding in response to low nitrogen culture conditions. IMPORTANCE Historically, coding variants were prioritized when searching for causal mechanisms driving adaptation of natural populations to stressful environments. However, the recent focus on noncoding variants demonstrated their ubiquitous role in adaptation. Here, we performed genome-wide regulatory variation profiles between two divergent yeast strains when facing nitrogen nutritional stress. The open chromatin availability of several regulatory regions changes in response to nitrogen scarcity. Importantly, we describe regulatory events that deviate between strains. Our results demonstrate a widespread variation in gene expression regulation between naturally occurring populations in response to stressful environments.
Collapse
|
220
|
Zhang W, Liao L, Xu J, Han Y, Li L. Genome-wide identification, characterization and expression analysis of MATE family genes in apple (Malus × domestica Borkh). BMC Genomics 2021; 22:632. [PMID: 34461821 PMCID: PMC8406601 DOI: 10.1186/s12864-021-07943-1] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Accepted: 08/12/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND As an important group of the multidrug efflux transporter family, the multidrug and toxic compound extrusion (MATE) family has a wide range of functions and is distributed in all kingdoms of living organisms. However, only two MATE genes in apple have been analyzed and genome-wide comprehensive analysis of MATE family is needed. RESULTS In this study, a total of 66 MATE (MdMATE) candidates encoding putative MATE transporters were identified in the apple genome. These MdMATE genes were classified into four groups by phylogenetic analysis with MATE genes in Arabidopsis. Synteny analysis reveals that whole genome duplication (WGD) and segmental duplication events played a major role in the expansion of MATE gene family in apple. MdMATE genes show diverse expression patterns in different tissues/organs and developmental stages. Analysis of cis-regulatory elements in MdMATE promoter regions indicates that the function of MdMATE genes is mainly related to stress response. Besides, the changes of gene expression levels upon different pathogen infections reveal that MdMATE genes are involved in biotic stress response. CONCLUSIONS In this work, we systematically identified MdMATE genes in apple genome using a set of bioinformatics approaches. Our comprehensive analysis provided valuable resources for improving disease resistance in apple and further functional characterization of MATE genes in other species.
Collapse
Affiliation(s)
- Weihan Zhang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, 430070, People's Republic of China
| | - Liao Liao
- CAS Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture, Wuhan Botanical Garden, The Innovative Academy of Seed Design, Chinese Academy of Sciences, Wuhan, 430074, People's Republic of China
| | - Jinsheng Xu
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, 430070, People's Republic of China
| | - Yuepeng Han
- CAS Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture, Wuhan Botanical Garden, The Innovative Academy of Seed Design, Chinese Academy of Sciences, Wuhan, 430074, People's Republic of China
| | - Li Li
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, 430070, People's Republic of China. .,Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan, 430070, People's Republic of China.
| |
Collapse
|
221
|
Wang L, Sun F, Wan ZY, Ye B, Wen Y, Liu H, Yang Z, Pang H, Meng Z, Fan B, Alfiko Y, Shen Y, Bai B, Lee MSQ, Piferrer F, Schartl M, Meyer A, Yue GH. Genomic Basis of Striking Fin Shapes and Colors in the Fighting Fish. Mol Biol Evol 2021; 38:3383-3396. [PMID: 33871625 PMCID: PMC8321530 DOI: 10.1093/molbev/msab110] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Resolving the genomic basis underlying phenotypic variations is a question of great importance in evolutionary biology. However, understanding how genotypes determine the phenotypes is still challenging. Centuries of artificial selective breeding for beauty and aggression resulted in a plethora of colors, long-fin varieties, and hyper-aggressive behavior in the air-breathing Siamese fighting fish (Betta splendens), supplying an excellent system for studying the genomic basis of phenotypic variations. Combining whole-genome sequencing, quantitative trait loci mapping, genome-wide association studies, and genome editing, we investigated the genomic basis of huge morphological variation in fins and striking differences in coloration in the fighting fish. Results revealed that the double tail, elephant ear, albino, and fin spot mutants each were determined by single major-effect loci. The elephant ear phenotype was likely related to differential expression of a potassium ion channel gene, kcnh8. The albinotic phenotype was likely linked to a cis-regulatory element acting on the mitfa gene and the double-tail mutant was suggested to be caused by a deletion in a zic1/zic4 coenhancer. Our data highlight that major loci and cis-regulatory elements play important roles in bringing about phenotypic innovations and establish Bettas as new powerful model to study the genomic basis of evolved changes.
Collapse
Affiliation(s)
- Le Wang
- Molecular Population Genetics & Breeding Group, Temasek Life Sciences Laboratory, Singapore, Singapore
| | - Fei Sun
- Molecular Population Genetics & Breeding Group, Temasek Life Sciences Laboratory, Singapore, Singapore
| | - Zi Yi Wan
- Molecular Population Genetics & Breeding Group, Temasek Life Sciences Laboratory, Singapore, Singapore
| | - Baoqing Ye
- Molecular Population Genetics & Breeding Group, Temasek Life Sciences Laboratory, Singapore, Singapore
| | - Yanfei Wen
- Molecular Population Genetics & Breeding Group, Temasek Life Sciences Laboratory, Singapore, Singapore
| | - Huiming Liu
- Molecular Population Genetics & Breeding Group, Temasek Life Sciences Laboratory, Singapore, Singapore
| | - Zituo Yang
- Molecular Population Genetics & Breeding Group, Temasek Life Sciences Laboratory, Singapore, Singapore
| | - Hongyan Pang
- Molecular Population Genetics & Breeding Group, Temasek Life Sciences Laboratory, Singapore, Singapore
| | - Zining Meng
- School of Life Sciences, Sun Yat-sen University, Guangzhou, China
| | - Bin Fan
- Department of Food and Environmental Engineering, Yangjiang Polytechnic, Yangjiang, China
| | - Yuzer Alfiko
- Biotech Lab, Wilmar International, Jakarta, Indonesia
| | - Yubang Shen
- Key Laboratory of Exploration and Utilization of Aquatic Genetic Resources, Shanghai Ocean University, Shanghai, China
| | - Bin Bai
- Molecular Population Genetics & Breeding Group, Temasek Life Sciences Laboratory, Singapore, Singapore
| | - May Shu Qing Lee
- Molecular Population Genetics & Breeding Group, Temasek Life Sciences Laboratory, Singapore, Singapore
| | - Francesc Piferrer
- Institute of Marine Sciences (ICM), Spanish National Research Council (CSIC), Barcelona, Spain
| | - Manfred Schartl
- Developmental Biochemistry, Biocenter, University of Wuerzburg, Wuerzburg, Germany
- The Xiphophorus Genetic Stock Center, Department of Chemistry and Biochemistry, Texas State University, San Marcos, TX, USA
| | - Axel Meyer
- Department of Biology, University of Konstanz, Konstanz, Germany
| | - Gen Hua Yue
- Molecular Population Genetics & Breeding Group, Temasek Life Sciences Laboratory, Singapore, Singapore
- Department of Biological Sciences, National University of Singapore, Singapore, Singapore
- School of Biological Sciences, Nanyang Technological University, Singapore, Singapore
| |
Collapse
|
222
|
Wang Y, Hua YP, Zhou T, Huang JY, Yue CP. Genomic identification of nitrogen assimilation-related genes and transcriptional characterization of their responses to nitrogen in allotetraploid rapeseed. Mol Biol Rep 2021; 48:5977-5992. [PMID: 34327662 DOI: 10.1007/s11033-021-06599-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2021] [Accepted: 07/25/2021] [Indexed: 11/25/2022]
Abstract
BACKGROUND Nitrogen (N) is an essential macronutrient to maintain plant growth and development. Plants absorb nitrate-N or ammonium-N in the environment and undergo reduction reactions catalyzed by nitrate reductase (NR), nitrite reductase (NIR), glutamine synthetase (GS), and glutamine oxoglutarate aminotransferase (GOGAT) within plants. METHODS AND RESULTS A total of 42 N assimilation-related genes (NAG) members were identified in rapeseed. Darwin's evolutionary pressure analysis showed that rapeseed NAGs underwent purification selection. Cis-element analysis revealed differences in the transcriptional regulation of NAGs between Arabidopsis and rapeseed. Expression analyses revealed that NRs were expressed mainly in old leaves, NIRs were expressed mainly in old leaves and lower stem peels, while the expression situation between different subfamilies of GSs and GOGATs was more complicated. CONCLUSIONS Differential expression of NAGs suggested that they might be involved in abiotic stresses. The above results greatly enriched our understanding of NAGs' molecular characteristics and provided central gene resources for NAGs-mediated NUE improvement in rapeseed.
Collapse
Affiliation(s)
- Yue Wang
- School of Agricultural Sciences, Zhengzhou University, Zhengzhou, 450001, China
| | - Ying-Peng Hua
- School of Agricultural Sciences, Zhengzhou University, Zhengzhou, 450001, China
| | - Ting Zhou
- School of Agricultural Sciences, Zhengzhou University, Zhengzhou, 450001, China
| | - Jin-Yong Huang
- School of Agricultural Sciences, Zhengzhou University, Zhengzhou, 450001, China
| | - Cai-Peng Yue
- School of Agricultural Sciences, Zhengzhou University, Zhengzhou, 450001, China.
| |
Collapse
|
223
|
Duan WJ, Liu ZH, Bai JF, Yuan SH, Li YM, Lu FK, Zhang TB, Sun JH, Zhang FT, Zhao CP, Zhang LP. Comprehensive analysis of formin gene family highlights candidate genes related to pollen cytoskeleton and male fertility in wheat (Triticum aestivum L.). BMC Genomics 2021; 22:570. [PMID: 34303338 PMCID: PMC8305537 DOI: 10.1186/s12864-021-07878-7] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Accepted: 07/01/2021] [Indexed: 11/19/2022] Open
Abstract
Background Formin, a highly conserved multi-domain protein, interacts with microfilaments and microtubules. Although specifically expressed formin genes in anthers are potentially significant in research on male sterility and hybrid wheat breeding, similar reports in wheat, especially in thermo-sensitive genic male sterile (TGMS) wheat, remain elusive. Results Herein, we systematically characterized the formin genes in TGMS wheat line BS366 named TaFormins (TaFHs) and predicted their functions in inducing stress response. In total, 25 TaFH genes were uncovered, majorly localized in 2A, 2B, and 2D chromosomes. According to the neighbor-joining (NJ) method, all TaFH proteins from wheat and other plants clustered in 6 sub-groups (A-F). The modeled 3D structures of TaFH1-A/B, TaFH2-A/B, TaFH3-A/B and TaFH3-B/D were validated. And different numbers of stress and hormone-responsive regulatory elements in their 1500 base pair promoter regions were contained in the TaFH genes copies. TaFHs had specific temporal and spatial expression characteristics, whereby TaFH1, TaFH4, and TaFH5 were expressed highly in the stamen of BS366. Besides, the accumulation of TaFHs was remarkably lower in a low-temperature sterile condition (Nanyang) than fertile condition (Beijing), particularly at the early stamen development stage. The pollen cytoskeleton of BS366 was abnormal in the three stages under sterile and fertile environments. Furthermore, under different stress levels, TaFHs expression could be induced by drought, salt, abscisic acid (ABA), salicylic acid (SA), methyl jasmonate (MeJA), indole-3-acetic acid (IAA), polyethylene glycol (PEG), and low temperature. Some miRNAs, including miR167, miR1120, and miR172, interacts with TaFH genes; thus, we constructed an interaction network between microRNAs, TaFHs, phytohormone responses, and distribution of cytoskeleton to reveal the regulatory association between upstream genes of TaFH family members and sterile. Conclusions Collectively, this comprehensive analysis provides novel insights into TaFHs and miRNA resources for wheat breeding. These findings are, therefore, valuable in understanding the mechanism of TGMS fertility conversion in wheat. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07878-7.
Collapse
Affiliation(s)
- Wen-Jing Duan
- Beijing Engineering and Technique Research Center for Hybrid Wheat, Beijing Academy of Agriculture and Forestry Sciences, Beijing, 100097, China.,The Municipal Key Laboratory of Molecular Genetic of Hybrid Wheat, Beijing, 100097, China.,College of Life Science, Capital Normal University, Beijing, 100048, China
| | - Zi-Han Liu
- Beijing Engineering and Technique Research Center for Hybrid Wheat, Beijing Academy of Agriculture and Forestry Sciences, Beijing, 100097, China.,The Municipal Key Laboratory of Molecular Genetic of Hybrid Wheat, Beijing, 100097, China
| | - Jian-Fang Bai
- Beijing Engineering and Technique Research Center for Hybrid Wheat, Beijing Academy of Agriculture and Forestry Sciences, Beijing, 100097, China.,The Municipal Key Laboratory of Molecular Genetic of Hybrid Wheat, Beijing, 100097, China
| | - Shao-Hua Yuan
- Beijing Engineering and Technique Research Center for Hybrid Wheat, Beijing Academy of Agriculture and Forestry Sciences, Beijing, 100097, China.,The Municipal Key Laboratory of Molecular Genetic of Hybrid Wheat, Beijing, 100097, China
| | - Yan-Mei Li
- Beijing Engineering and Technique Research Center for Hybrid Wheat, Beijing Academy of Agriculture and Forestry Sciences, Beijing, 100097, China.,The Municipal Key Laboratory of Molecular Genetic of Hybrid Wheat, Beijing, 100097, China
| | - Feng-Kun Lu
- Beijing Engineering and Technique Research Center for Hybrid Wheat, Beijing Academy of Agriculture and Forestry Sciences, Beijing, 100097, China.,The Municipal Key Laboratory of Molecular Genetic of Hybrid Wheat, Beijing, 100097, China.,College of Life Science, Capital Normal University, Beijing, 100048, China
| | - Tian-Bao Zhang
- Beijing Engineering and Technique Research Center for Hybrid Wheat, Beijing Academy of Agriculture and Forestry Sciences, Beijing, 100097, China.,The Municipal Key Laboratory of Molecular Genetic of Hybrid Wheat, Beijing, 100097, China
| | - Jia-Hui Sun
- Beijing Engineering and Technique Research Center for Hybrid Wheat, Beijing Academy of Agriculture and Forestry Sciences, Beijing, 100097, China.,The Municipal Key Laboratory of Molecular Genetic of Hybrid Wheat, Beijing, 100097, China
| | - Feng-Ting Zhang
- Beijing Engineering and Technique Research Center for Hybrid Wheat, Beijing Academy of Agriculture and Forestry Sciences, Beijing, 100097, China. .,The Municipal Key Laboratory of Molecular Genetic of Hybrid Wheat, Beijing, 100097, China.
| | - Chang-Ping Zhao
- Beijing Engineering and Technique Research Center for Hybrid Wheat, Beijing Academy of Agriculture and Forestry Sciences, Beijing, 100097, China. .,The Municipal Key Laboratory of Molecular Genetic of Hybrid Wheat, Beijing, 100097, China.
| | - Li-Ping Zhang
- Beijing Engineering and Technique Research Center for Hybrid Wheat, Beijing Academy of Agriculture and Forestry Sciences, Beijing, 100097, China. .,The Municipal Key Laboratory of Molecular Genetic of Hybrid Wheat, Beijing, 100097, China.
| |
Collapse
|
224
|
Johnson KM, Jones HR, Casas SM, La Peyre JF, Kelly MW. Transcriptomic signatures of temperature adaptation in the eastern oyster Crassostrea virginica. J Evol Biol 2021; 34:1212-1224. [PMID: 33837581 DOI: 10.1111/jeb.13789] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2020] [Revised: 02/09/2021] [Accepted: 03/18/2021] [Indexed: 12/20/2022]
Abstract
The large geographic distribution of the eastern oyster, Crassostrea virginica, makes it an ideal species to test how populations have adapted to latitudinal gradients in temperature. Despite inhabiting distinct thermal regimes, populations of C. virginica near the species' southern and northern geographic range show no population differences in their physiological response to temperature. In this study, we used comparative transcriptomics to understand how oysters from either end of the species' range maintain enantiostasis across three acclimation temperatures (10, 20, and 30°C). With this approach, we identified genes that were differentially expressed in response to temperature between individuals of C. virginica collected from New Brunswick, Canada and Louisiana, USA. We observed a core set of genes whose expression responded to temperature in both populations, but also an even larger set of genes with expression patterns that were unique to each population. Intriguingly, the genes with population-specific responses to temperature had elevated FST and Ka/Ks ratios compared to the genome-wide average. In contrast, genes showing only a response to temperature were found to only have elevated FST values suggesting that divergent FST may be due to selection on linked regulatory regions rather than positive selection on protein coding regions. Taken together, our results suggest that, despite coarse-scale physiological similarities, natural selection has shaped divergent gene expression responses to temperature in geographically separated populations of this broadly eurythermal marine invertebrate.
Collapse
Affiliation(s)
- Kevin M Johnson
- Department of Biological Sciences, Louisiana State University, Baton Rouge, LA, USA.,Center for Coastal Marine Sciences, California Polytechnic State University, San Luis Obispo, CA, USA.,California Sea Grant, University of California San Diego, La Jolla, CA, USA
| | - Hollis R Jones
- Department of Biological Sciences, Louisiana State University, Baton Rouge, LA, USA.,Department of Animal Science, University of California, Davis, CA, USA
| | - Sandra M Casas
- School of Animal Sciences, Louisiana State University Ag Center, Baton Rouge, LA, USA
| | - Jerome F La Peyre
- School of Animal Sciences, Louisiana State University Ag Center, Baton Rouge, LA, USA
| | - Morgan W Kelly
- Department of Biological Sciences, Louisiana State University, Baton Rouge, LA, USA
| |
Collapse
|
225
|
Wang Y, Jin R, Liu C, Gao Y, Deng X, Wan H, Li J. Functional characterization of the transcription factors AhR and ARNT in Nilaparvata lugens. PESTICIDE BIOCHEMISTRY AND PHYSIOLOGY 2021; 176:104875. [PMID: 34119220 DOI: 10.1016/j.pestbp.2021.104875] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/25/2021] [Revised: 05/09/2021] [Accepted: 05/10/2021] [Indexed: 06/12/2023]
Abstract
In the present study, the aryl hydrocarbon receptor (AhR) and aryl hydrocarbon receptor nuclear translocator (ARNT) of Nilaparvata lugens were cloned and identified. The NlAhR and NlARNT expression levels significantly increased after imidacloprid, etofenprox and isoprocarb treatments. Knockdowns of NlAhR and NlARNT increased the susceptibility of N. lugens to imidacloprid, etofenprox and isoprocarb, and the detoxification enzyme activities were also significantly decreased. In addition, NlCYP301A1, NlGSTt1 and NlCarE7 were significantly down-regulated after injections of dsNlAhR and dsNlARNT, with the NlCarE7 expression decreasing by greater than 80%. Moreover, after knocking down NlCarE7, the susceptibility of N. lugens to etofenprox and isoprocarb significantly increased. Both NlAhR and NlARNT bound the NlCarE7 promoter and significantly enhanced the transcriptional activity. Our research revealed the functional roles of transcription factors NlAhR and NlARNT in the detoxification metabolism of N. lugens. The results provide a theoretical basis for the pest management and comprehensive control of N. lugens and increase our knowledge of insect toxicology.
Collapse
Affiliation(s)
- Yue Wang
- Hubei Insect Resources Utilization and Sustainable Pest Management Key Laboratory, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, PR China
| | - Ruoheng Jin
- Hubei Insect Resources Utilization and Sustainable Pest Management Key Laboratory, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, PR China
| | - Chaoya Liu
- Hubei Insect Resources Utilization and Sustainable Pest Management Key Laboratory, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, PR China
| | - Yuanyuan Gao
- Hubei Insect Resources Utilization and Sustainable Pest Management Key Laboratory, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, PR China
| | - Xiaoqian Deng
- Hubei Insect Resources Utilization and Sustainable Pest Management Key Laboratory, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, PR China
| | - Hu Wan
- Hubei Insect Resources Utilization and Sustainable Pest Management Key Laboratory, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, PR China
| | - Jianhong Li
- Hubei Insect Resources Utilization and Sustainable Pest Management Key Laboratory, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, PR China.
| |
Collapse
|
226
|
Genome-Wide Identification, Structure Characterization, Expression Pattern Profiling, and Substrate Specificity of the Metal Tolerance Protein Family in Canavalia rosea (Sw.) DC. PLANTS 2021; 10:plants10071340. [PMID: 34209283 PMCID: PMC8309081 DOI: 10.3390/plants10071340] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/31/2021] [Revised: 06/22/2021] [Accepted: 06/29/2021] [Indexed: 11/16/2022]
Abstract
Plant metal tolerance proteins (MTPs) play key roles in heavy metal absorption and homeostasis in plants. By using genome-wide and phylogenetic approaches, the origin and diversification of MTPs from Canavalia rosea (Sw.) DC. was explored. Canavalia rosea (bay bean) is an extremophile halophyte with strong adaptability to seawater and drought and thereby shows specific metal tolerance with a potential phytoremediation ability. However, MTP genes in leguminous plants remain poorly understood. In our study, a total of 12 MTP genes were identified in C. rosea. Multiple sequence alignments showed that all CrMTP proteins possessed the conserved transmembrane domains (TM1 to TM6) and could be classified into three subfamilies: Zn-CDFs (five members), Fe/Zn-CDFs (five members), and Mn-CDFs (two members). Promoter cis-acting element analyses revealed that a distinct number and composition of heavy metal regulated elements and other stress-responsive elements existed in different promoter regions of CrMTPs. Analysis of transcriptome data revealed organ-specific expression of CrMTP genes and the involvement of this family in heavy metal stress responses and adaptation of C. rosea to extreme coral reef environments. Furthermore, the metal-specific activity of several functionally unknown CrMTPs was investigated in yeast. These results will contribute to uncovering the potential functions and molecular mechanisms of heavy metal absorption, translocation, and accumulation in C. rosea plants.
Collapse
|
227
|
Fotoukkiaii SM, Wybouw N, Kurlovs AH, Tsakireli D, Pergantis SA, Clark RM, Vontas J, Van Leeuwen T. High-resolution genetic mapping reveals cis-regulatory and copy number variation in loci associated with cytochrome P450-mediated detoxification in a generalist arthropod pest. PLoS Genet 2021; 17:e1009422. [PMID: 34153029 PMCID: PMC8248744 DOI: 10.1371/journal.pgen.1009422] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2021] [Revised: 07/01/2021] [Accepted: 05/28/2021] [Indexed: 12/11/2022] Open
Abstract
Chemical control strategies are driving the evolution of pesticide resistance in pest populations. Understanding the genetic mechanisms of these evolutionary processes is of crucial importance to develop sustainable resistance management strategies. The acaricide pyflubumide is one of the most recently developed mitochondrial complex II inhibitors with a new mode of action that specifically targets spider mite pests. In this study, we characterize the molecular basis of pyflubumide resistance in a highly resistant population of the spider mite Tetranychus urticae. Classical genetic crosses indicated that pyflubumide resistance was incompletely recessive and controlled by more than one gene. To identify resistance loci, we crossed the resistant population to a highly susceptible T. urticae inbred strain and propagated resulting populations with and without pyflubumide exposure for multiple generations in an experimental evolution set-up. High-resolution genetic mapping by a bulked segregant analysis approach led to the identification of three quantitative trait loci (QTL) linked to pyflubumide resistance. Two QTLs were found on the first chromosome and centered on the cytochrome P450 CYP392A16 and a cluster of CYP392E6-8 genes. Comparative transcriptomics revealed a consistent overexpression of CYP392A16 and CYP392E8 in the experimental populations that were selected for pyflubumide resistance. We further corroborated the involvement of CYP392A16 in resistance by in vitro functional expression and metabolism studies. Collectively, these experiments uncovered that CYP392A16 N-demethylates the toxic carboxamide form of pyflubumide to a non-toxic compound. A third QTL coincided with cytochrome P450 reductase (CPR), a vital component of cytochrome P450 metabolism. We show here that the resistant population harbors three gene copies of CPR and that this copy number variation is associated with higher mRNA abundance. Together, we provide evidence for detoxification of pyflubumide by cytochrome P450s that is likely synergized by gene amplification of CPR. Our understanding of the causal genetic variants that drive the evolution of quantitative traits, such as polygenic pesticide resistance, remains very limited. Here, we followed a high-resolution genetic mapping approach to localize the genetic variants that cause pyflubumide resistance in the two-spotted spider mite Tetranychus urticae. Three well-supported QTL were uncovered and pointed towards a major role for cytochrome P450-mediated detoxification. Cis-regulatory variation for cytochrome P450s was observed, and in vitro cytochrome P450 experiments showed that pyflubumide was metabolized into a non-toxic derivate. A third QTL centered on cytochrome P450 reductase (CPR), which is required for cytochrome P450 activity, and is amplified in pyflubumide resistant populations. Our results indicate that pyflubumide resistance is mediated by cytochrome P450 detoxification that is enhanced by gene amplification at the CPR locus.
Collapse
Affiliation(s)
- Seyedeh Masoumeh Fotoukkiaii
- Department of Evolutionary and Population Biology, Institute for Biodiversity and Ecosystem Dynamics, University of Amsterdam, Amsterdam, the Netherlands
| | - Nicky Wybouw
- Laboratory of Agrozoology, Department of Plants and Crops, Faculty of Bioscience Engineering, Ghent University, Ghent, Belgium
- Terrestrial Ecology Unit, Department of Biology, Faculty of Sciences, Ghent University, Ghent, Belgium
| | - Andre H. Kurlovs
- Laboratory of Agrozoology, Department of Plants and Crops, Faculty of Bioscience Engineering, Ghent University, Ghent, Belgium
- School of Biological Sciences, University of Utah, Salt Lake City, Utah, United States of America
| | - Dimitra Tsakireli
- Institute of Molecular Biology & Biotechnology, Foundation for Research & Technology Hellas, Heraklion, Crete, Greece
- Laboratory of Pesticide Science, Department of Crop Science, Agricultural University of Athens, Athens, Greece
| | | | - Richard M. Clark
- School of Biological Sciences, University of Utah, Salt Lake City, Utah, United States of America
- Henry Eyring Center for Cell and Genome Science, University of Utah, Salt Lake City, Utah, United States of America
| | - John Vontas
- Institute of Molecular Biology & Biotechnology, Foundation for Research & Technology Hellas, Heraklion, Crete, Greece
- Laboratory of Pesticide Science, Department of Crop Science, Agricultural University of Athens, Athens, Greece
| | - Thomas Van Leeuwen
- Laboratory of Agrozoology, Department of Plants and Crops, Faculty of Bioscience Engineering, Ghent University, Ghent, Belgium
- * E-mail:
| |
Collapse
|
228
|
Fruciano C, Franchini P, Jones JC. Capturing the rapidly evolving study of adaptation. J Evol Biol 2021; 34:856-865. [PMID: 34145685 DOI: 10.1111/jeb.13871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2021] [Revised: 05/12/2021] [Accepted: 05/12/2021] [Indexed: 11/30/2022]
Abstract
Research on the genomics of adaptation is rapidly changing. In the last few decades, progress in this area has been driven by methodological advances, not only in the way increasingly large amounts of molecular data are generated (e.g. with high-throughput sequencing), but also in the way these data are analysed. This includes a growing appreciation and quantitative treatment of covariation among units within the same data type (e.g. genes) or across data types (e.g. genes and phenotypes). The development and adoption of more and more integrative tools have resulted in richer and more interesting empirical work. This special issue - comprising methodological, empirical, and review papers - aims to capture a 'snapshot' of this rapidly evolving field. We discuss in particular three important themes in the study of adaptation: the genetic architecture of adaptive variation, protein-coding and regulatory changes, and parallel evolution. We highlight how more traditional key themes in the study of genetic architecture (e.g. the number of loci underlying adaptive traits and the distribution of their effects) are now being complemented by other factors (e.g. how patterns of linkage and number of loci interact to affect the ability to adapt). Similarly, apart from addressing the relative importance of protein-coding and regulatory changes, we now have the tools to look in-depth at specific types of regulatory variation to gain a clearer picture of regulatory networks. Finally, parallel evolution has always been central to the study of adaptation, but now we are often able to address the question of whether - and to what extent - parallelism at the organismal or phenotypic level is matched by parallelism at the genetic level. Perhaps most importantly, we can now determine what mechanisms are driving parallelism (or lack thereof) across levels of biological organization. All these recent methodological developments open up new directions for future studies of adaptive changes across traits, levels of biological organization, demographic contexts and time scales.
Collapse
Affiliation(s)
- Carmelo Fruciano
- National Research Council - Institute of Marine Biological Resources and Biotechnologies, Messina, Italy.,Institut de biologie de l'Ecole normale supérieure (IBENS), Ecole normale supérieure, CNRS, PSL Université Paris, Paris, France.,School of Biological Sciences, University of Portsmouth, Portsmouth, UK
| | - Paolo Franchini
- Department of Biology, University of Konstanz, Konstanz, Germany
| | - Julia C Jones
- School of Biology and Environmental Science, University College Dublin, Dublin, Ireland
| |
Collapse
|
229
|
Zrimec J, Buric F, Kokina M, Garcia V, Zelezniak A. Learning the Regulatory Code of Gene Expression. Front Mol Biosci 2021; 8:673363. [PMID: 34179082 PMCID: PMC8223075 DOI: 10.3389/fmolb.2021.673363] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2021] [Accepted: 05/24/2021] [Indexed: 11/13/2022] Open
Abstract
Data-driven machine learning is the method of choice for predicting molecular phenotypes from nucleotide sequence, modeling gene expression events including protein-DNA binding, chromatin states as well as mRNA and protein levels. Deep neural networks automatically learn informative sequence representations and interpreting them enables us to improve our understanding of the regulatory code governing gene expression. Here, we review the latest developments that apply shallow or deep learning to quantify molecular phenotypes and decode the cis-regulatory grammar from prokaryotic and eukaryotic sequencing data. Our approach is to build from the ground up, first focusing on the initiating protein-DNA interactions, then specific coding and non-coding regions, and finally on advances that combine multiple parts of the gene and mRNA regulatory structures, achieving unprecedented performance. We thus provide a quantitative view of gene expression regulation from nucleotide sequence, concluding with an information-centric overview of the central dogma of molecular biology.
Collapse
Affiliation(s)
- Jan Zrimec
- Department of Biology and Biological Engineering, Chalmers University of Technology, Gothenburg, Sweden
| | - Filip Buric
- Department of Biology and Biological Engineering, Chalmers University of Technology, Gothenburg, Sweden
| | - Mariia Kokina
- Department of Biology and Biological Engineering, Chalmers University of Technology, Gothenburg, Sweden
- Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Victor Garcia
- School of Life Sciences and Facility Management, Zurich University of Applied Sciences, Wädenswil, Switzerland
| | - Aleksej Zelezniak
- Department of Biology and Biological Engineering, Chalmers University of Technology, Gothenburg, Sweden
- Science for Life Laboratory, Stockholm, Sweden
| |
Collapse
|
230
|
Mendelevich A, Vinogradova S, Gupta S, Mironov AA, Sunyaev SR, Gimelbrant AA. Replicate sequencing libraries are important for quantification of allelic imbalance. Nat Commun 2021; 12:3370. [PMID: 34099647 PMCID: PMC8184992 DOI: 10.1038/s41467-021-23544-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2020] [Accepted: 04/30/2021] [Indexed: 12/13/2022] Open
Abstract
A sensitive approach to quantitative analysis of transcriptional regulation in diploid organisms is analysis of allelic imbalance (AI) in RNA sequencing (RNA-seq) data. A near-universal practice in such studies is to prepare and sequence only one library per RNA sample. We present theoretical and experimental evidence that data from a single RNA-seq library is insufficient for reliable quantification of the contribution of technical noise to the observed AI signal; consequently, reliance on one-replicate experimental design can lead to unaccounted-for variation in error rates in allele-specific analysis. We develop a computational approach, Qllelic, that accurately accounts for technical noise by making use of replicate RNA-seq libraries. Testing on new and existing datasets shows that application of Qllelic greatly decreases false positive rate in allele-specific analysis while conserving appropriate signal, and thus greatly improves reproducibility of AI estimates. We explore sources of technical overdispersion in observed AI signal and conclude by discussing design of RNA-seq studies addressing two biologically important questions: quantification of transcriptome-wide AI in one sample, and differential analysis of allele-specific expression between samples.
Collapse
Affiliation(s)
- Asia Mendelevich
- Skolkovo Institute of Science and Technology, Moscow, Russia.
- Center for Cancer Systems Biology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, USA.
| | - Svetlana Vinogradova
- Center for Cancer Systems Biology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, USA
| | - Saumya Gupta
- Center for Cancer Systems Biology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, USA
- Broad Institute of Harvard and MIT, Cambridge, USA
| | - Andrey A Mironov
- Lomonosov Moscow State University, Faculty of Bioengineering and Bioinformatics, Moscow, Russia
- Institute of Information Transmission Problems, Russian Academy of Sciences, Moscow, Russia
| | - Shamil R Sunyaev
- Department of Biomedical Informatics, Harvard Medical School, Boston, USA
- Division of Genetics, Brigham and Women's Hospital, Boston, USA
| | - Alexander A Gimelbrant
- Center for Cancer Systems Biology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, USA.
- Broad Institute of Harvard and MIT, Cambridge, USA.
| |
Collapse
|
231
|
Qin J, Ye F, Xu L, Zhou X, Crickmore N, Zhou X, Zhang Y, Guo Z. A cis-Acting Mutation in the PxABCG1 Promoter Is Associated with Cry1Ac Resistance in Plutella xylostella (L.). Int J Mol Sci 2021; 22:6106. [PMID: 34198929 PMCID: PMC8201282 DOI: 10.3390/ijms22116106] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2021] [Revised: 06/02/2021] [Accepted: 06/03/2021] [Indexed: 12/13/2022] Open
Abstract
The molecular mechanisms of insect resistance to Cry toxins generated from the bacterium Bacillus thuringiensis (Bt) urgently need to be elucidated to enable the improvement and sustainability of Bt-based products. Although downregulation of the expression of midgut receptor genes is a pivotal mechanism of insect resistance to Bt Cry toxins, the underlying transcriptional regulation of these genes remains elusive. Herein, we unraveled the regulatory mechanism of the downregulation of the ABC transporter gene PxABCG1 (also called Pxwhite), a functional midgut receptor of the Bt Cry1Ac toxin in Plutella xylostella. The PxABCG1 promoters of Cry1Ac-susceptible and Cry1Ac-resistant strains were cloned and analyzed, and they showed clear differences in activity. Subsequently, a dual-luciferase reporter assay, a yeast one-hybrid (Y1H) assay, and RNA interference (RNAi) experiments demonstrated that a cis-mutation in a binding site of the Hox transcription factor Antennapedia (Antp) decreased the promoter activity of the resistant strain and eliminated the binding and regulation of Antp, thereby enhancing the resistance of P. xylostella to the Cry1Ac toxin. These results advance our knowledge of the roles of cis- and trans-regulatory variations in the regulation of midgut Cry receptor genes and the evolution of Bt resistance, contributing to a more complete understanding of the Bt resistance mechanism.
Collapse
Affiliation(s)
- Jianying Qin
- Longping Branch, Graduate School of Hunan University, Changsha 410125, China; (J.Q.); (X.Z.)
- Department of Plant Protection, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China; (F.Y.); (L.X.); (Y.Z.)
| | - Fan Ye
- Department of Plant Protection, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China; (F.Y.); (L.X.); (Y.Z.)
| | - Linzheng Xu
- Department of Plant Protection, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China; (F.Y.); (L.X.); (Y.Z.)
| | - Xuguo Zhou
- Department of Entomology, University of Kentucky, Lexington, KY 40546-0091, USA;
| | - Neil Crickmore
- School of Life Sciences, University of Sussex, Brighton BN1 9QG, UK;
| | - Xiaomao Zhou
- Longping Branch, Graduate School of Hunan University, Changsha 410125, China; (J.Q.); (X.Z.)
| | - Youjun Zhang
- Department of Plant Protection, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China; (F.Y.); (L.X.); (Y.Z.)
| | - Zhaojiang Guo
- Department of Plant Protection, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China; (F.Y.); (L.X.); (Y.Z.)
| |
Collapse
|
232
|
Manna M, Thakur T, Chirom O, Mandlik R, Deshmukh R, Salvi P. Transcription factors as key molecular target to strengthen the drought stress tolerance in plants. PHYSIOLOGIA PLANTARUM 2021; 172:847-868. [PMID: 33180329 DOI: 10.1111/ppl.13268] [Citation(s) in RCA: 92] [Impact Index Per Article: 30.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/20/2020] [Revised: 10/23/2020] [Accepted: 11/07/2020] [Indexed: 05/03/2023]
Abstract
Amid apprehension of global climate change, crop plants are inevitably confronted with a myriad of abiotic stress factors during their growth that inflicts a serious threat to their development and overall productivity. These abiotic stresses comprise extreme temperature, pH, high saline soil, and drought stress. Among different abiotic stresses, drought is considered the most calamitous stressor with its serious impact on the crops' yield stability. The development of climate-resilient crops that withstands reduced water availability is a major focus of the scientific fraternity to ensure the food security of the sharply increasing population. Numerous studies aim to recognize the key regulators of molecular and biochemical processes associated with drought stress tolerance response. A few potential candidates are now considered as promising targets for crop improvement. Transcription factors act as a key regulatory switch controlling the gene expression of diverse biological processes and, eventually, the metabolic processes. Understanding the role and regulation of the transcription factors will facilitate the crop improvement strategies intending to develop and deliver agronomically-superior crops. Therefore, in this review, we have emphasized the molecular avenues of the transcription factors that can be exploited to engineer drought tolerance potential in crop plants. We have discussed the molecular role of several transcription factors, such as basic leucine zipper (bZIP), dehydration responsive element binding (DREB), DNA binding with one finger (DOF), heat shock factor (HSF), MYB, NAC, TEOSINTE BRANCHED1/CYCLOIDEA/PCF (TCP), and WRKY. We have also highlighted candidate transcription factors that can be used for the development of drought-tolerant crops.
Collapse
Affiliation(s)
- Mrinalini Manna
- National Institute of Plant Genome Research, New Delhi, India
| | - Tanika Thakur
- Agriculture Biotechnology Department, National Agri-Food Biotechnology Institute, Mohali, Punjab, India
| | - Oceania Chirom
- National Institute of Plant Genome Research, New Delhi, India
| | - Rushil Mandlik
- Agriculture Biotechnology Department, National Agri-Food Biotechnology Institute, Mohali, Punjab, India
| | - Rupesh Deshmukh
- Agriculture Biotechnology Department, National Agri-Food Biotechnology Institute, Mohali, Punjab, India
| | - Prafull Salvi
- Agriculture Biotechnology Department, National Agri-Food Biotechnology Institute, Mohali, Punjab, India
| |
Collapse
|
233
|
Wu TY, Goh H, Azodi CB, Krishnamoorthi S, Liu MJ, Urano D. Evolutionarily conserved hierarchical gene regulatory networks for plant salt stress response. NATURE PLANTS 2021; 7:787-799. [PMID: 34045707 DOI: 10.1038/s41477-021-00929-7] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/19/2020] [Accepted: 04/23/2021] [Indexed: 06/12/2023]
Abstract
Plant cells constantly alter their gene expression profiles to respond to environmental fluctuations. These continuous adjustments are regulated by multi-hierarchical networks of transcription factors. To understand how such gene regulatory networks (GRNs) have stabilized evolutionarily while allowing for species-specific responses, we compare the GRNs underlying salt response in the early-diverging and late-diverging plants Marchantia polymorpha and Arabidopsis thaliana. Salt-responsive GRNs, constructed on the basis of the temporal transcriptional patterns in the two species, share common trans-regulators but exhibit an evolutionary divergence in cis-regulatory sequences and in the overall network sizes. In both species, WRKY-family transcription factors and their feedback loops serve as central nodes in salt-responsive GRNs. The divergent cis-regulatory sequences of WRKY-target genes are probably associated with the expansion in network size, linking salt stress to tissue-specific developmental and physiological responses. The WRKY modules and highly linked WRKY feedback loops have been preserved widely in other plants, including rice, while keeping their binding-motif sequences mutable. Together, the conserved trans-regulators and the quickly evolving cis-regulatory sequences allow salt-responsive GRNs to adapt over a long evolutionary timescale while maintaining some consistent regulatory structure. This strategy may benefit plants as they adapt to changing environments.
Collapse
Affiliation(s)
- Ting-Ying Wu
- Temasek Life Sciences Laboratory, National University of Singapore, Singapore, Singapore.
| | - HonZhen Goh
- Temasek Life Sciences Laboratory, National University of Singapore, Singapore, Singapore
| | - Christina B Azodi
- Department of Plant Biology, Michigan State University, East Lansing, MI, USA
| | - Shalini Krishnamoorthi
- Temasek Life Sciences Laboratory, National University of Singapore, Singapore, Singapore
| | - Ming-Jung Liu
- Biotechnology Center in Southern Taiwan, Academia Sinica, Tainan, Taiwan
- Agricultural Biotechnology Research Center, Academia Sinica, Taipei, Taiwan
| | - Daisuke Urano
- Temasek Life Sciences Laboratory, National University of Singapore, Singapore, Singapore.
- Department of Biological Sciences, National University of Singapore, Singapore, Singapore.
- Singapore-MIT Alliance for Research and Technology, Singapore, Singapore.
| |
Collapse
|
234
|
Massel K, Lam Y, Wong ACS, Hickey LT, Borrell AK, Godwin ID. Hotter, drier, CRISPR: the latest edit on climate change. TAG. THEORETICAL AND APPLIED GENETICS. THEORETISCHE UND ANGEWANDTE GENETIK 2021; 134:1691-1709. [PMID: 33420514 DOI: 10.1007/s00122-020-03764-0] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/16/2020] [Accepted: 12/30/2020] [Indexed: 05/23/2023]
Abstract
Integrating CRISPR/Cas9 genome editing into modern breeding programs for crop improvement in cereals. Global climate trends in many agricultural regions have been rapidly changing over the past decades, and major advances in global food systems are required to ensure food security in the face of these emerging challenges. With increasing climate instability due to warmer temperatures and rising CO2 levels, the productivity of global agriculture will continue to be negatively impacted. To combat these growing concerns, creative approaches will be required, utilising all the tools available to produce more robust and tolerant crops with increased quality and yields under more extreme conditions. The integration of genome editing and transgenics into current breeding strategies is one promising solution to accelerate genetic gains through targeted genetic modifications, producing crops that can overcome the shifting climate realities. This review focuses on how revolutionary genome editing tools can be directly implemented into breeding programs for cereal crop improvement to rapidly counteract many of the issues affecting agriculture production in the years to come.
Collapse
Affiliation(s)
- Karen Massel
- Centre for Crop Science, Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, Brisbane, QLD, 4072, Australia.
| | - Yasmine Lam
- Centre for Crop Science, Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, Brisbane, QLD, 4072, Australia
| | - Albert C S Wong
- Centre for Crop Science, Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, Brisbane, QLD, 4072, Australia
| | - Lee T Hickey
- Centre for Crop Science, Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, Brisbane, QLD, 4072, Australia
| | - Andrew K Borrell
- Centre for Crop Science, Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, Brisbane, QLD, 4072, Australia
| | - Ian D Godwin
- Centre for Crop Science, Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, Brisbane, QLD, 4072, Australia
| |
Collapse
|
235
|
Marand AP, Chen Z, Gallavotti A, Schmitz RJ. A cis-regulatory atlas in maize at single-cell resolution. Cell 2021; 184:3041-3055.e21. [PMID: 33964211 DOI: 10.1101/2020.09.27.315499] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2020] [Revised: 03/04/2021] [Accepted: 04/07/2021] [Indexed: 05/22/2023]
Abstract
cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.
Collapse
Affiliation(s)
| | - Zongliang Chen
- Waksman Institute, Rutgers University, Piscataway, NJ 08854, USA
| | - Andrea Gallavotti
- Waksman Institute, Rutgers University, Piscataway, NJ 08854, USA; Department of Plant Biology, Rutgers University, New Brunswick, NJ 08901, USA
| | - Robert J Schmitz
- Department of Genetics, University of Georgia, Athens, GA 30602, USA.
| |
Collapse
|
236
|
Wang M, Chen B, Zhou W, Xie L, Wang L, Zhang Y, Zhang Q. Genome-wide identification and expression analysis of the AT-hook Motif Nuclear Localized gene family in soybean. BMC Genomics 2021; 22:361. [PMID: 34006214 PMCID: PMC8132359 DOI: 10.1186/s12864-021-07687-y] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2021] [Accepted: 05/04/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Soybean is an important legume crop and has significant agricultural and economic value. Previous research has shown that the AT-Hook Motif Nuclear Localized (AHL) gene family is highly conserved in land plants, playing crucial roles in plant growth and development. To date, however, the AHL gene family has not been studied in soybean. RESULTS To investigate the roles played by the AHL gene family in soybean, genome-wide identification, expression patterns and gene structures were performed to analyze. We identified a total of 63 AT-hook motif genes, which were characterized by the presence of the AT-hook motif and PPC domain in soybean. The AT-hook motif genes were distributed on 18 chromosomes and formed two distinct clades (A and B), as shown by phylogenetic analysis. All the AHL proteins were further classified into three types (I, II and III) based on the AT-hook motif. Type-I was belonged to Clade-A, while Type-II and Type-III were belonged to Clade-B. Our results also showed that the main type of duplication in the soybean AHL gene family was segmented duplication event. To discern whether the AHL gene family was involved in stress response in soybean, we performed cis-acting elements analysis and found that AHL genes were associated with light responsiveness, anaerobic induction, MYB and gibberellin-responsiveness elements. This suggest that AHL genes may participate in plant development and mediate stress response. Moreover, a co-expression network analysis showed that the AHL genes were also involved in energy transduction, and the associated with the gibberellin pathway and nuclear entry signal pathways in soybean. Transcription analysis revealed that AHL genes in Jack and Williams82 have a common expression pattern and are mostly expressed in roots, showing greater sensitivity under drought and submergence stress. Hence, the AHL gene family mainly reacts on mediating stress responses in the roots and provide comprehensive information for further understanding of the AT-hook motif gene family-mediated stress response in soybean. CONCLUSION Sixty-three AT-hook motif genes were identified in the soybean genome. These genes formed into two distinct phylogenetic clades and belonged to three different types. Cis-acting elements and co-expression network analyses suggested that AHL genes participated in significant biological processes. This work provides important theoretical basis for the understanding of AHLs biological functions in soybean.
Collapse
Affiliation(s)
- Min Wang
- College of Life Sciences, Northeast Forestry University, Harbin, 150040, People's Republic of China
- Key Laboratory of Saline-Alkali Vegetative Ecology Restoration, Ministry of Education, College of Life Science, Northeast Forestry University, Harbin, 150040, China
| | - Bowei Chen
- College of Life Sciences, Northeast Forestry University, Harbin, 150040, People's Republic of China
- Key Laboratory of Saline-Alkali Vegetative Ecology Restoration, Ministry of Education, College of Life Science, Northeast Forestry University, Harbin, 150040, China
| | - Wei Zhou
- College of Life Sciences, Northeast Forestry University, Harbin, 150040, People's Republic of China
- Key Laboratory of Saline-Alkali Vegetative Ecology Restoration, Ministry of Education, College of Life Science, Northeast Forestry University, Harbin, 150040, China
| | - Linan Xie
- College of Life Sciences, Northeast Forestry University, Harbin, 150040, People's Republic of China
- Key Laboratory of Saline-Alkali Vegetative Ecology Restoration, Ministry of Education, College of Life Science, Northeast Forestry University, Harbin, 150040, China
| | - Lishan Wang
- College of Life Sciences, Northeast Forestry University, Harbin, 150040, People's Republic of China
- Key Laboratory of Saline-Alkali Vegetative Ecology Restoration, Ministry of Education, College of Life Science, Northeast Forestry University, Harbin, 150040, China
| | - Yonglan Zhang
- College of Life Sciences, Northeast Forestry University, Harbin, 150040, People's Republic of China
- Key Laboratory of Saline-Alkali Vegetative Ecology Restoration, Ministry of Education, College of Life Science, Northeast Forestry University, Harbin, 150040, China
| | - Qingzhu Zhang
- College of Life Sciences, Northeast Forestry University, Harbin, 150040, People's Republic of China.
- State Key Laboratory of Tree Genetics and Breeding, Northeast Forestry University, Harbin, 150040, People's Republic of China.
| |
Collapse
|
237
|
Jovanovic VM, Sarfert M, Reyna-Blanco CS, Indrischek H, Valdivia DI, Shelest E, Nowick K. Positive Selection in Gene Regulatory Factors Suggests Adaptive Pleiotropic Changes During Human Evolution. Front Genet 2021; 12:662239. [PMID: 34079582 PMCID: PMC8166252 DOI: 10.3389/fgene.2021.662239] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2021] [Accepted: 04/19/2021] [Indexed: 01/09/2023] Open
Abstract
Gene regulatory factors (GRFs), such as transcription factors, co-factors and histone-modifying enzymes, play many important roles in modifying gene expression in biological processes. They have also been proposed to underlie speciation and adaptation. To investigate potential contributions of GRFs to primate evolution, we analyzed GRF genes in 27 publicly available primate genomes. Genes coding for zinc finger (ZNF) proteins, especially ZNFs with a Krüppel-associated box (KRAB) domain were the most abundant TFs in all genomes. Gene numbers per TF family differed between all species. To detect signs of positive selection in GRF genes we investigated more than 3,000 human GRFs with their more than 70,000 orthologs in 26 non-human primates. We implemented two independent tests for positive selection, the branch-site-model of the PAML suite and aBSREL of the HyPhy suite, focusing on the human and great ape branch. Our workflow included rigorous procedures to reduce the number of false positives: excluding distantly similar orthologs, manual corrections of alignments, and considering only genes and sites detected by both tests for positive selection. Furthermore, we verified the candidate sites for selection by investigating their variation within human and non-human great ape population data. In order to approximately assign a date to positively selected sites in the human lineage, we analyzed archaic human genomes. Our work revealed with high confidence five GRFs that have been positively selected on the human lineage and one GRF that has been positively selected on the great ape lineage. These GRFs are scattered on different chromosomes and have been previously linked to diverse functions. For some of them a role in speciation and/or adaptation can be proposed based on the expression pattern or association with human diseases, but it seems that they all contributed independently to human evolution. Four of the positively selected GRFs are KRAB-ZNF proteins, that induce changes in target genes co-expression and/or through arms race with transposable elements. Since each positively selected GRF contains several sites with evidence for positive selection, we suggest that these GRFs participated pleiotropically to phenotypic adaptations in humans.
Collapse
Affiliation(s)
- Vladimir M Jovanovic
- Human Biology and Primate Evolution, Freie Universität Berlin, Berlin, Germany.,Bioinformatics Solution Center, Freie Universität Berlin, Berlin, Germany
| | - Melanie Sarfert
- Human Biology and Primate Evolution, Freie Universität Berlin, Berlin, Germany
| | - Carlos S Reyna-Blanco
- Department of Biology, University of Fribourg, Fribourg, Switzerland.,Swiss Institute of Bioinformatics, Fribourg, Switzerland
| | - Henrike Indrischek
- Max Planck Institute of Molecular Cell Biology and Genetics, Dresden, Germany.,Max Planck Institute for the Physics of Complex Systems, Dresden, Germany.,Center for Systems Biology Dresden, Dresden, Germany
| | - Dulce I Valdivia
- Evolutionary Genomics Laboratory and Genome Topology and Regulation Laboratory, Genetic Engineering Department, Center for Research and Advanced Studies of the National Polytechnic Institute (CINVESTAV-Irapuato), Irapuato, Mexico
| | - Ekaterina Shelest
- Centre for Enzyme Innovation, University of Portsmouth, Portsmouth, United Kingdom
| | - Katja Nowick
- Human Biology and Primate Evolution, Freie Universität Berlin, Berlin, Germany
| |
Collapse
|
238
|
Huang Y, Lack JB, Hoppel GT, Pool JE. Parallel and Population-specific Gene Regulatory Evolution in Cold-Adapted Fly Populations. Genetics 2021; 218:6275754. [PMID: 33989401 PMCID: PMC8864734 DOI: 10.1093/genetics/iyab077] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2020] [Accepted: 05/10/2021] [Indexed: 11/15/2022] Open
Abstract
Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find evidence for parallel expression evolution between them, with stronger parallelism at larval and adult stages than for pupae. We also implement a flexible method to estimate cis- vs trans-encoded contributions to expression or splicing differences at the adult stage. The apparent contributions of cis- vs trans-regulation to adaptive evolution vary substantially among population pairs. While two of three population pairs show a greater enrichment of cis-regulatory differences among adaptation candidates, trans-regulatory differences are more likely to be implicated in parallel expression changes between population pairs. Genes with significant cis-effects are enriched for signals of elevated genetic differentiation between cold- and warm-adapted populations, suggesting that they are potential targets of local adaptation. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.
Collapse
Affiliation(s)
- Yuheng Huang
- Laboratory of Genetics, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA.,Department of Ecology and Evolutionary Biology, University of California, Irvine, Irvine, CA 92697, USA
| | - Justin B Lack
- Laboratory of Genetics, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA.,Advanced Biomedical Computational Science, Frederick National Laboratory for Cancer Research, Frederick, MD 21701, USA
| | - Grant T Hoppel
- Laboratory of Genetics, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - John E Pool
- Laboratory of Genetics, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
239
|
Campana MG, Corvelo A, Shelton J, Callicrate TE, Bunting KL, Riley-Gillis B, Wos F, DeGrazia J, Jarvis ED, Fleischer RC. Adaptive Radiation Genomics of Two Ecologically Divergent Hawai'ian Honeycreepers: The 'akiapōlā'au and the Hawai'i 'amakihi. J Hered 2021; 111:21-32. [PMID: 31723957 DOI: 10.1093/jhered/esz057] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2019] [Accepted: 09/30/2019] [Indexed: 12/14/2022] Open
Abstract
The Hawai'ian honeycreepers (drepanids) are a classic example of adaptive radiation: they adapted to a variety of novel dietary niches, evolving a wide range of bill morphologies. Here we investigated genomic diversity, demographic history, and genes involved in bill morphology phenotypes in 2 honeycreepers: the 'akiapōlā'au (Hemignathus wilsoni) and the Hawai'i 'amakihi (Chlorodrepanis virens). The 'akiapōlā'au is an endangered island endemic, filling the "woodpecker" niche by using a unique bill morphology, while the Hawai'i 'amakihi is a dietary generalist common on the islands of Hawai'i and Maui. We de novo sequenced the 'akiapōlā'au genome and compared it to the previously sequenced 'amakihi genome. The 'akiapōlā'au is far less heterozygous and has a smaller effective population size than the 'amakihi, which matches expectations due to its smaller census population and restricted ecological niche. Our investigation revealed genomic islands of divergence, which may be involved in the honeycreeper radiation. Within these islands of divergence, we identified candidate genes (including DLK1, FOXB1, KIF6, MAML3, PHF20, RBP1, and TIMM17A) that may play a role in honeycreeper adaptations. The gene DLK1, previously shown to influence Darwin's finch bill size, may be related to honeycreeper bill morphology evolution, while the functions of the other candidates remain unknown.
Collapse
Affiliation(s)
- Michael G Campana
- Center for Conservation Genomics, Smithsonian National Zoo and Conservation Biology Institute, Washington, DC
| | | | | | - Taylor E Callicrate
- Center for Conservation Genomics, Smithsonian National Zoo and Conservation Biology Institute, Washington, DC.,Species Conservation Toolkit Initiative, Chicago Zoological Society, Brookfield, IL
| | | | | | - Frank Wos
- New York Genome Center, New York, NY
| | | | - Erich D Jarvis
- The Rockefeller University, New York, NY.,Howard Hughes Medical Institute, Chevy Chase, MD
| | - Robert C Fleischer
- Center for Conservation Genomics, Smithsonian National Zoo and Conservation Biology Institute, Washington, DC
| |
Collapse
|
240
|
Fuqua T, Jordan J, Halavatyi A, Tischer C, Richter K, Crocker J. An open-source semi-automated robotics pipeline for embryo immunohistochemistry. Sci Rep 2021; 11:10314. [PMID: 33986394 PMCID: PMC8119710 DOI: 10.1038/s41598-021-89676-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2020] [Accepted: 04/19/2021] [Indexed: 11/09/2022] Open
Abstract
A significant challenge for developmental systems biology is balancing throughput with controlled conditions that minimize experimental artifacts. Large-scale developmental screens such as unbiased mutagenesis surveys have been limited in their applicability to embryonic systems, as the technologies for quantifying precise expression patterns in whole animals has not kept pace with other sequencing-based technologies. Here, we outline an open-source semi-automated pipeline to chemically fixate, stain, and 3D-image Drosophila embryos. Central to this pipeline is a liquid handling robot, Flyspresso, which automates the steps of classical embryo fixation and staining. We provide the schematics and an overview of the technology for an engineer or someone equivalently trained to reproduce and further improve upon Flyspresso, and highlight the Drosophila embryo fixation and colorimetric or antibody staining protocols. Additionally, we provide a detailed overview and stepwise protocol for our adaptive-feedback pipeline for automated embryo imaging on confocal microscopes. We demonstrate the efficiency of this pipeline compared to classical techniques, and how it can be repurposed or scaled to other protocols and biological systems. We hope our pipeline will serve as a platform for future research, allowing a broader community of users to build, execute, and share similar experiments.
Collapse
Affiliation(s)
- Timothy Fuqua
- European Molecular Biology Laboratory, Heidelberg, Germany.,Collaboration for Joint PhD Degree Between EMBL and Heidelberg University, Faculty of Biosciences, Heidelberg, Germany
| | - Jeff Jordan
- Janelia Research Campus, 19700 Helix Dr, Ashburn, VA, 20147, USA
| | | | | | | | - Justin Crocker
- European Molecular Biology Laboratory, Heidelberg, Germany.
| |
Collapse
|
241
|
Okubo S, Terauchi K, Okada S, Saito Y, Yamaura T, Misaka T, Nakajima KI, Abe K, Asakura T. De novo transcriptome analysis and comparative expression profiling of genes associated with the taste-modifying protein neoculin in Curculigo latifolia and Curculigo capitulata fruits. BMC Genomics 2021; 22:347. [PMID: 33985426 PMCID: PMC8120819 DOI: 10.1186/s12864-021-07674-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2020] [Accepted: 05/05/2021] [Indexed: 11/29/2022] Open
Abstract
Background Curculigo latifolia is a perennial plant endogenous to Southeast Asia whose fruits contain the taste-modifying protein neoculin, which binds to sweet receptors and makes sour fruits taste sweet. Although similar to snowdrop (Galanthus nivalis) agglutinin (GNA), which contains mannose-binding sites in its sequence and 3D structure, neoculin lacks such sites and has no lectin activity. Whether the fruits of C. latifolia and other Curculigo plants contain neoculin and/or GNA family members was unclear. Results Through de novo RNA-seq assembly of the fruits of C. latifolia and the related C. capitulata and detailed analysis of the expression patterns of neoculin and neoculin-like genes in both species, we assembled 85,697 transcripts from C. latifolia and 76,775 from C. capitulata using Trinity and annotated them using public databases. We identified 70,371 unigenes in C. latifolia and 63,704 in C. capitulata. In total, 38.6% of unigenes from C. latifolia and 42.6% from C. capitulata shared high similarity between the two species. We identified ten neoculin-related transcripts in C. latifolia and 15 in C. capitulata, encoding both the basic and acidic subunits of neoculin in both plants. We aligned these 25 transcripts and generated a phylogenetic tree. Many orthologs in the two species shared high similarity, despite the low number of common genes, suggesting that these genes likely existed before the two species diverged. The relative expression levels of these genes differed considerably between the two species: the transcripts per million (TPM) values of neoculin genes were 60 times higher in C. latifolia than in C. capitulata, whereas those of GNA family members were 15,000 times lower in C. latifolia than in C. capitulata. Conclusions The genetic diversity of neoculin-related genes strongly suggests that neoculin genes underwent duplication during evolution. The marked differences in their expression profiles between C. latifolia and C. capitulata may be due to mutations in regions involved in transcriptional regulation. Comprehensive analysis of the genes expressed in the fruits of these two Curculigo species helped elucidate the origin of neoculin at the molecular level. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07674-3.
Collapse
Affiliation(s)
- Satoshi Okubo
- The Yamashina Botanical Research Institute, Nippon Shinyaku Co., Ltd., Oyake Sakanotsuji-cho 39, Yamashina-ku, Kyoto, 607-8182, Japan
| | - Kaede Terauchi
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1, Yayoi, Bunkyo-ku, Tokyo, 113-8657, Japan
| | - Shinji Okada
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1, Yayoi, Bunkyo-ku, Tokyo, 113-8657, Japan
| | - Yoshikazu Saito
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1, Yayoi, Bunkyo-ku, Tokyo, 113-8657, Japan
| | - Takao Yamaura
- The Yamashina Botanical Research Institute, Nippon Shinyaku Co., Ltd., Oyake Sakanotsuji-cho 39, Yamashina-ku, Kyoto, 607-8182, Japan
| | - Takumi Misaka
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1, Yayoi, Bunkyo-ku, Tokyo, 113-8657, Japan
| | - Ken-Ichiro Nakajima
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1, Yayoi, Bunkyo-ku, Tokyo, 113-8657, Japan.,Present address: Division of Endocrinology and Metabolism, National Institute for Physiological Sciences, 38 Nishigonaka, Myodaiji, Okazaki, Aichi, 444-8585, Japan
| | - Keiko Abe
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1, Yayoi, Bunkyo-ku, Tokyo, 113-8657, Japan.,Kanagawa Institute of Industrial Science and Technology (KISTEC), 3-25-13 Tonomachi, Kawasaki-ku, Kawasaki, Kanagawa, 210-0821, Japan
| | - Tomiko Asakura
- Graduate School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1, Yayoi, Bunkyo-ku, Tokyo, 113-8657, Japan.
| |
Collapse
|
242
|
Fong SL, Capra JA. Modeling the evolutionary architectures of transcribed human enhancer sequences reveals distinct origins, functions, and associations with human-trait variation. Mol Biol Evol 2021; 38:3681-3696. [PMID: 33973014 PMCID: PMC8382917 DOI: 10.1093/molbev/msab138] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Despite the importance of gene regulatory enhancers in human biology and evolution, we lack a comprehensive model of enhancer evolution and function. This substantially limits our understanding of the genetic basis of species divergence and our ability to interpret the effects of noncoding variants on human traits. To explore enhancer sequence evolution and its relationship to regulatory function, we traced the evolutionary origins of transcribed human enhancer sequences with activity across diverse tissues and cellular contexts from the FANTOM5 consortium. The transcribed enhancers are enriched for sequences of a single evolutionary age (“simple” evolutionary architectures) compared with enhancers that are composites of sequences of multiple evolutionary ages (“complex” evolutionary architectures), likely indicating constraint against genomic rearrangements. Complex enhancers are older, more pleiotropic, and more active across species than simple enhancers. Genetic variants within complex enhancers are also less likely to associate with human traits and biochemical activity. Transposable-element-derived sequences (TEDS) have made diverse contributions to enhancers of both architectures; the majority of TEDS are found in enhancers with simple architectures, while a minority have remodeled older sequences to create complex architectures. Finally, we compare the evolutionary architectures of transcribed enhancers with histone-mark-defined enhancers. Our results reveal that most human transcribed enhancers are ancient sequences of a single age, and thus the evolution of most human enhancers was not driven by increases in evolutionary complexity over time. Our analyses further suggest that considering enhancer evolutionary histories provides context that can aid interpretation of the effects of variants on enhancer function. Based on these results, we propose a framework for analyzing enhancer evolutionary architecture.
Collapse
Affiliation(s)
- Sarah L Fong
- Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN, USA
| | - John A Capra
- Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN, USA.,Department of Biological Sciences, Vanderbilt University, Nashville, TN, USA.,Bakar Computational Health Sciences Institute and Department of Epidemiology and Biostatistics, University of California, San Francisco, USA
| |
Collapse
|
243
|
The Expression Characteristics of NPF Genes and Their Response to Vernalization and Nitrogen Deficiency in Rapeseed. Int J Mol Sci 2021; 22:ijms22094944. [PMID: 34066572 PMCID: PMC8125141 DOI: 10.3390/ijms22094944] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Revised: 05/01/2021] [Accepted: 05/03/2021] [Indexed: 11/24/2022] Open
Abstract
The NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER FAMILY (NPF) genes, initially characterized as nitrate or peptide transporters in plants, are involved in the transport of a large variety of substrates, including amino acids, nitrate, auxin (IAA), jasmonates (JAs), abscisic acid (ABA) and gibberellins (GAs) and glucosinolates. A total of 169 potential functional NPF genes were excavated in Brassica napus, and they showed diversified expression patterns in 90 different organs or tissues based on transcriptome profile data. The complex time-serial expression changes were found for most functional NPF genes in the development process of leaves, silique walls and seeds, which indicated that the expression of Brassica napus NPF (BnaNPF) genes may respond to altered phytohormone and secondary metabolite content through combining with promoter element enrichment analysis. Furthermore, many BnaNPF genes were detected to respond to vernalization with two different patterns, and 20 BnaNPF genes responded to nitrate deficiency. These results will provide useful information for further investigation of the biological function of BnaNPF genes for growth and development in rapeseed.
Collapse
|
244
|
Gao Z, Li H, Yang X, Yang P, Chen J, Shi T. Biased allelic expression in tissues of F1 hybrids between tropical and temperate lotus (Nelumbo nuicfera). PLANT MOLECULAR BIOLOGY 2021; 106:207-220. [PMID: 33738679 DOI: 10.1007/s11103-021-01138-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/03/2020] [Accepted: 03/05/2021] [Indexed: 06/12/2023]
Abstract
The genome-wide allele-specific expression in F1 hybrids from the cross of tropical and temperate lotus unveils how cis-regulatory divergences affect genes in key pathways related to ecotypic divergence. Genetic variation, particularly cis-regulatory variation, plays a crucial role in phenotypic variation and adaptive evolution in plants. Temperate and tropical lotus, the two ecotypes of Nelumbo nucifera, show distinction in the degree of rhizome enlargement, which is associated with winter dormancy. To understand the roles of genome-wide cis-regulatory divergences on adaptive evolution of temperate and tropical lotus (Nelumbo nucifera), here we performed allele-specific expression (ASE) analyses on the tissues including flowers, leaves and rhizome from F1 hybrids of tropical and temperate lotus. For all investigated tissues in F1s, about 36% of genes showed ASE and about 3% of genes showed strong consistent ASE. Most of ASEs were biased towards the tropical parent in all surveyed samples, indicating that the tropical genome might be dominant over the temperate genome in gene expression of tissues from their F1 hybrids. We found that promoter sequences with similar allelic expression are more conserved than genes with significant or conditional ASE, suggesting the cis-regulatory sequence divergence underlie the allelic expression bias. We further uncovered biased genes being related to phenotypic differentiation between two lotus ecotypes, especially metabolic and phytohormone-related pathways in the rhizome. Overall, our study provides a global landscape of cis-regulatory variations between two lotus ecotypes and highlights their roles in rhizome growth variation for the climatic adaptation.
Collapse
Affiliation(s)
- Zhiyan Gao
- Key Laboratory of Aquatic Botany and Watershed Ecology, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, 430074, China
- Center of Conservation Biology, Core Botanical Gardens, Chinese Academy of Sciences, Wuhan, 430074, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Hui Li
- Key Laboratory of Aquatic Botany and Watershed Ecology, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, 430074, China
- Center of Conservation Biology, Core Botanical Gardens, Chinese Academy of Sciences, Wuhan, 430074, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Xingyu Yang
- Wuhan Institute of Landscape Architecture, Wuhan, 430081, China
| | - Pingfang Yang
- State Key Laboratory of Biocatalysis and Enzyme Engineering, School of Life Sciences, Hubei University, Wuhan, 430062, China
| | - Jinming Chen
- Key Laboratory of Aquatic Botany and Watershed Ecology, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, 430074, China.
- Center of Conservation Biology, Core Botanical Gardens, Chinese Academy of Sciences, Wuhan, 430074, China.
| | - Tao Shi
- Key Laboratory of Aquatic Botany and Watershed Ecology, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, 430074, China.
- Center of Conservation Biology, Core Botanical Gardens, Chinese Academy of Sciences, Wuhan, 430074, China.
| |
Collapse
|
245
|
Froehlich JJ, Uyar B, Herzog M, Theil K, Glažar P, Akalin A, Rajewsky N. Parallel genetics of regulatory sequences using scalable genome editing in vivo. Cell Rep 2021; 35:108988. [PMID: 33852857 DOI: 10.1016/j.celrep.2021.108988] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2020] [Revised: 01/13/2021] [Accepted: 03/23/2021] [Indexed: 12/27/2022] Open
Abstract
How regulatory sequences control gene expression is fundamental for explaining phenotypes in health and disease. Regulatory elements must ultimately be understood within their genomic environment and development- or tissue-specific contexts. Because this is technically challenging, few regulatory elements have been characterized in vivo. Here, we use inducible Cas9 and multiplexed guide RNAs to create hundreds of mutations in enhancers/promoters and 3' UTRs of 16 genes in C. elegans. Our software crispr-DART analyzes indel mutations in targeted DNA sequencing. We quantify the impact of mutations on expression and fitness by targeted RNA sequencing and DNA sampling. When applying our approach to the lin-41 3' UTR, generating hundreds of mutants, we find that the two adjacent binding sites for the miRNA let-7 can regulate lin-41 expression independently of each other. Finally, we map regulatory genotypes to phenotypic traits for several genes. Our approach enables parallel analysis of regulatory sequences directly in animals.
Collapse
Affiliation(s)
- Jonathan J Froehlich
- Systems Biology of Gene Regulatory Elements, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Hannoversche Str. 28, 10115 Berlin, Germany
| | - Bora Uyar
- Bioinformatics and Omics Data Science Platform, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Hannoversche Str. 28, 10115 Berlin, Germany
| | - Margareta Herzog
- Systems Biology of Gene Regulatory Elements, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Hannoversche Str. 28, 10115 Berlin, Germany
| | - Kathrin Theil
- Systems Biology of Gene Regulatory Elements, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Hannoversche Str. 28, 10115 Berlin, Germany
| | - Petar Glažar
- Systems Biology of Gene Regulatory Elements, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Hannoversche Str. 28, 10115 Berlin, Germany
| | - Altuna Akalin
- Bioinformatics and Omics Data Science Platform, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Hannoversche Str. 28, 10115 Berlin, Germany
| | - Nikolaus Rajewsky
- Systems Biology of Gene Regulatory Elements, Berlin Institute for Medical Systems Biology, Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Hannoversche Str. 28, 10115 Berlin, Germany.
| |
Collapse
|
246
|
Gu H, Zhang P, Xu M, Liang D. Amplicon genome fishing (AGF): a rapid and efficient method for sequencing target cis-regulatory regions in nonmodel organisms. Mol Genet Genomics 2021; 296:527-539. [PMID: 33797587 DOI: 10.1007/s00438-021-01775-0] [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: 01/27/2021] [Accepted: 03/15/2021] [Indexed: 10/21/2022]
Abstract
Cis-regulatory sequences play a crucial role in regulating gene expression and are evolutionary hot spots that drive phenotypic divergence among organisms. Sequencing some cis-regulatory regions of interest in many different species is common in comparative genetic studies. For nonmodel organisms lacking genomic data, genome walking is often the preferred method for this type of application. However, applying genome walking will be laborious and time-consuming when the number of cis-regulatory regions and species to be analyzed is large. In this study, we propose a novel method called amplicon genome fishing (AGF), which can isolate and sequence cis-regulatory regions of interest for any organism. The main idea of the AGF method is to use fragments amplified from the target cis-regulatory regions as enrichment baits to capture and sequence the whole target cis-regulatory regions from genomic library pools. Unlike genome walking, the AGF method is based on hybridization capture and high-throughput sequencing, which makes this method rapid and efficient for projects where some cis-regulatory regions have to be sequenced for many species. We used human amplicons as capture baits and successfully sequenced five target enhancer regions of Homo sapiens, Mus musculus, Gallus gallus, and Xenopus tropicalis, proving the feasibility and repeatability of AGF. To show the utility of the AGF method in real studies, we used it to sequence the ZRS enhancer, a cis-regulatory region associated with the limb loss of snakes, for twenty-three vertebrate species (includes many limbless species never sequenced before). The newly obtained ZRS sequences provide new perspectives into the relationship between the ZRS enhancer's evolution and limb loss in major tetrapod lineages.
Collapse
Affiliation(s)
- HanMei Gu
- State Key Laboratory of Biocontrol, Higher Education Mega Center, School of Life Sciences, College of Ecology and Evolution, Sun Yat-Sen University, #434, Guangzhou, 510006, China
| | - Peng Zhang
- State Key Laboratory of Biocontrol, Higher Education Mega Center, School of Life Sciences, College of Ecology and Evolution, Sun Yat-Sen University, #434, Guangzhou, 510006, China
| | - ManHao Xu
- State Key Laboratory of Biocontrol, Higher Education Mega Center, School of Life Sciences, College of Ecology and Evolution, Sun Yat-Sen University, #434, Guangzhou, 510006, China
| | - Dan Liang
- State Key Laboratory of Biocontrol, Higher Education Mega Center, School of Life Sciences, College of Ecology and Evolution, Sun Yat-Sen University, #434, Guangzhou, 510006, China.
| |
Collapse
|
247
|
Díaz-Valenzuela E, Sawers RH, Cibrián-Jaramillo A. Cis- and Trans-Regulatory Variations in the Domestication of the Chili Pepper Fruit. Mol Biol Evol 2021; 37:1593-1603. [PMID: 32031611 PMCID: PMC7253206 DOI: 10.1093/molbev/msaa027] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
The process of domestication requires the rapid transformation of the wild morphology into the cultivated forms that humans select for. This process often takes place through changes in the regulation of genes, yet, there is no definite pattern on the role of cis- and trans-acting regulatory variations in the domestication of the fruit among crops. Using allele-specific expression and network analyses, we characterized the regulatory patterns and the inheritance of gene expression in wild and cultivated accessions of chili pepper, a crop with remarkable fruit morphological variation. We propose that gene expression differences associated to the cultivated form are best explained by cis-regulatory hubs acting through trans-regulatory cascades. We show that in cultivated chili, the expression of genes associated with fruit morphology is partially recessive with respect to those in the wild relative, consistent with the hybrid fruit phenotype. Decreased expression of fruit maturation and growth genes in cultivated chili suggest that selection for loss-of-function took place in its domestication. Trans-regulatory changes underlie the majority of the genes showing regulatory divergence and had larger effect sizes on gene expression than cis-regulatory variants. Network analysis of selected cis-regulated genes, including ARP9 and MED25, indicated their interaction with many transcription factors involved in organ growth and fruit ripening. Differentially expressed genes linked to cis-regulatory variants and their interactions with downstream trans-acting genes have the potential to drive the morphological differences observed between wild and cultivated fruits and provide an attractive mechanism of morphological transformation during the domestication of the chili pepper.
Collapse
Affiliation(s)
- Erik Díaz-Valenzuela
- Ecological and Evolutionary Genomics Laboratory, Unidad de Genomica Avanzada (Langebio), Irapuato, Guanajuato, México
| | - Ruairidh H Sawers
- Department of Plant Science, The Pennsylvania State University, University Park State College, University Park, PA
| | - Angélica Cibrián-Jaramillo
- Ecological and Evolutionary Genomics Laboratory, Unidad de Genomica Avanzada (Langebio), Irapuato, Guanajuato, México
| |
Collapse
|
248
|
Genome-wide identification and expression pattern analysis of the ribonuclease T2 family in Eucommia ulmoides. Sci Rep 2021; 11:6900. [PMID: 33767357 PMCID: PMC7994793 DOI: 10.1038/s41598-021-86337-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Accepted: 03/11/2021] [Indexed: 11/13/2022] Open
Abstract
The 2′,3′-cycling ribonuclease (RNase) genes are catalysts of RNA cleavage and include the RNase T2 gene family. RNase T2 genes perform important roles in plants and have been conserved in the genome of eukaryotic organisms. In this study we identified 21 EURNS genes in Eucommia ulmoides Oliver (E. ulmoides) and analyzed their structure, chromosomal location, phylogenetic tree, gene duplication, stress-related cis-elements, and expression patterns in different tissues. The length of 21 predicted EURNS proteins ranged from 143 to 374 amino acids (aa), their molecular weight (MW) ranged from 16.21 to 42.38 kDa, and their isoelectric point (PI) value ranged from 5.08 to 9.09. Two classifications (class I and class III) were obtained from the conserved domains analysis and phylogenetic tree. EURNS proteins contained a total of 15 motifs. Motif 1, motif 2, motif 3, and motif 7 were distributed in multiple sequences and were similar to the conserved domain of RNase T2. EURNS genes with similar structure and the predicted EURNS proteins with conserved motif compositions are in the same group in the phylogenetic tree. The results of RT-PCR and transcription data showed that EURNS genes have tissue-specific expression and exhibited obvious trends in different developmental stages. Gene duplication analysis results indicated that segment duplication may be the dominant duplication mode in this gene family. This study provides a theoretical basis for research on the RNase T2 gene family and lays a foundation for the further study of EURNS genes.
Collapse
|
249
|
Human-chimpanzee fused cells reveal cis-regulatory divergence underlying skeletal evolution. Nat Genet 2021; 53:467-476. [PMID: 33731941 PMCID: PMC8038968 DOI: 10.1038/s41588-021-00804-3] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Accepted: 01/26/2021] [Indexed: 01/06/2023]
Abstract
Gene regulatory divergence is thought to play a central role in determining human-specific traits. However, our ability to link divergent regulation to divergent phenotypes is limited. Here, we utilized human-chimpanzee hybrid induced pluripotent stem cells to study gene expression separating these species. The tetraploid hybrid cells allowed us to separate cis- from trans-regulatory effects, and to control for non-genetic confounding factors. We differentiated these cells into cranial neural crest cells (CNCCs), the primary cell type giving rise to the face. We discovered evidence of lineage-specific selection on the hedgehog signaling pathway, including a human-specific 6-fold down-regulation of EVC2 (LIMBIN), a key hedgehog gene. Inducing a similar down-regulation of EVC2 substantially reduced hedgehog signaling output. Mice and humans lacking functional EVC2 show striking phenotypic parallels to human-chimpanzee craniofacial differences, suggesting that the regulatory divergence of hedgehog signaling may have contributed to the unique craniofacial morphology of humans.
Collapse
|
250
|
Dao FY, Lv H, Su W, Sun ZJ, Huang QL, Lin H. iDHS-Deep: an integrated tool for predicting DNase I hypersensitive sites by deep neural network. Brief Bioinform 2021; 22:6158360. [PMID: 33751027 DOI: 10.1093/bib/bbab047] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2020] [Revised: 01/28/2021] [Accepted: 01/29/2021] [Indexed: 01/09/2023] Open
Abstract
DNase I hypersensitive site (DHS) refers to the hypersensitive region of chromatin for the DNase I enzyme. It is an important part of the noncoding region and contains a variety of regulatory elements, such as promoter, enhancer, and transcription factor-binding site, etc. Moreover, the related locus of disease (or trait) are usually enriched in the DHS regions. Therefore, the detection of DHS region is of great significance. In this study, we develop a deep learning-based algorithm to identify whether an unknown sequence region would be potential DHS. The proposed method showed high prediction performance on both training datasets and independent datasets in different cell types and developmental stages, demonstrating that the method has excellent superiority in the identification of DHSs. Furthermore, for the convenience of related wet-experimental researchers, the user-friendly web-server iDHS-Deep was established at http://lin-group.cn/server/iDHS-Deep/, by which users can easily distinguish DHS and non-DHS and obtain the corresponding developmental stage ofDHS.
Collapse
Affiliation(s)
- Fu-Ying Dao
- Informational Biology at University of Electronic Science and Technology of China, China
| | - Hao Lv
- Informational Biology at University of Electronic Science and Technology of China, China
| | - Wei Su
- Informational Biology at University of Electronic Science and Technology of China, China
| | - Zi-Jie Sun
- Informational Biology at University of Electronic Science and Technology of China, China
| | - Qin-Lai Huang
- Informational Biology at University of Electronic Science and Technology of China, China
| | - Hao Lin
- Informational Biology at University of Electronic Science and Technology of China, China
| |
Collapse
|