1
|
Singh P, Kumar K, Jha AK, Yadava P, Pal M, Rakshit S, Singh I. Global gene expression profiling under nitrogen stress identifies key genes involved in nitrogen stress adaptation in maize (Zea mays L.). Sci Rep 2022; 12:4211. [PMID: 35273237 PMCID: PMC8913646 DOI: 10.1038/s41598-022-07709-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2021] [Accepted: 02/14/2022] [Indexed: 11/25/2022] Open
Abstract
Maize is a heavy consumer of fertilizer nitrogen (N) which not only results in the high cost of cultivation but may also lead to environmental pollution. Therefore, there is a need to develop N-use efficient genotypes, a prerequisite for which is a greater understanding of N-deficiency stress adaptation. In this study, comparative transcriptome analysis was performed using leaf and root tissues from contrasting inbred lines, viz., DMI 56 (tolerant to N stress) and DMI 81 (susceptible to N stress) to delineate the differentially expressed genes (DEGs) under low-N stress. The contrasting lines were grown hydroponically in modified Hoagland solution having either sufficient- or deficient-N, followed by high-throughput RNA-sequencing. A total of 8 sequencing libraries were prepared and 88–97% of the sequenced raw reads were mapped to the reference B73 maize genome. Genes with a p value ≤ 0.05 and fold change of ≥ 2.0 or ≤ − 2 were considered as DEGs in various combinations performed between susceptible and tolerant genotypes. DEGs were further classified into different functional categories and pathways according to their putative functions. Gene Ontology based annotation of these DEGs identified three different functional categories: biological processes, molecular function, and cellular component. The KEGG and Mapman based analysis revealed that most of the DEGs fall into various metabolic pathways, biosynthesis of secondary metabolites, signal transduction, amino acid metabolism, N-assimilation and metabolism, and starch metabolism. Some of the key genes involved in N uptake (high-affinity nitrate transporter 2.2 and 2.5), N assimilation and metabolism (glutamine synthetase, asparagine synthetase), redox homeostasis (SOD, POX), and transcription factors (MYB36, AP2-EREBP) were found to be highly expressed in the tolerant genotype compared to susceptible one. The candidate genes identified in the present study might be playing a pivotal role in low-N stress adaptation in maize and hence could be useful in augmenting further research on N metabolism and development of N-deficiency tolerant maize cultivars.
Collapse
Affiliation(s)
- Prabha Singh
- Indian Council of Agricultural Research-Indian Institute of Maize Research, Pusa Campus, New Delhi, 110012, India.,Indian Council of Agricultural Research-Indian Agricultural Research Institute, Pusa Campus, New Delhi, 110012, India.,Indian Council of Agricultural Research-Indian Grassland and Fodder Research Institute, Jhansi, 284003, India
| | - Krishan Kumar
- Indian Council of Agricultural Research-Indian Institute of Maize Research, Pusa Campus, New Delhi, 110012, India
| | - Abhishek Kumar Jha
- Indian Council of Agricultural Research-Indian Institute of Maize Research, Pusa Campus, New Delhi, 110012, India
| | - Pranjal Yadava
- Indian Council of Agricultural Research-Indian Agricultural Research Institute, Pusa Campus, New Delhi, 110012, India
| | - Madan Pal
- Indian Council of Agricultural Research-Indian Agricultural Research Institute, Pusa Campus, New Delhi, 110012, India
| | - Sujay Rakshit
- Indian Council of Agricultural Research-Indian Institute of Maize Research, Pusa Campus, New Delhi, 110012, India
| | - Ishwar Singh
- Indian Council of Agricultural Research-Indian Institute of Maize Research, Pusa Campus, New Delhi, 110012, India.
| |
Collapse
|
2
|
Ho EKH, Bellis ES, Calkins J, Adrion JR, Latta IV LC, Schaack S. Engines of change: Transposable element mutation rates are high and variable within Daphnia magna. PLoS Genet 2021; 17:e1009827. [PMID: 34723969 PMCID: PMC8594854 DOI: 10.1371/journal.pgen.1009827] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2021] [Revised: 11/16/2021] [Accepted: 09/16/2021] [Indexed: 12/22/2022] Open
Abstract
Transposable elements (TEs) represent a major portion of most eukaryotic genomes, yet little is known about their mutation rates or how their activity is shaped by other evolutionary forces. Here, we compare short- and long-term patterns of genome-wide mutation accumulation (MA) of TEs among 9 genotypes from three populations of Daphnia magna from across a latitudinal gradient. While the overall proportion of the genome comprised of TEs is highly similar among genotypes from Finland, Germany, and Israel, populations are distinguishable based on patterns of insertion site polymorphism. Our direct rate estimates indicate TE movement is highly variable (net rates ranging from -11.98 to 12.79 x 10-5 per copy per generation among genotypes), differing both among populations and TE families. Although gains outnumber losses when selection is minimized, both types of events appear to be highly deleterious based on their low frequency in control lines where propagation is not limited to random, single-progeny descent. With rate estimates 4 orders of magnitude higher than base substitutions, TEs clearly represent a highly mutagenic force in the genome. Quantifying patterns of intra- and interspecific variation in TE mobility with and without selection provides insight into a powerful mechanism generating genetic variation in the genome.
Collapse
Affiliation(s)
- Eddie K. H. Ho
- Department of Biology, Reed College, Portland, Oregon, United States of America
| | - Emily S. Bellis
- Department of Biology, Reed College, Portland, Oregon, United States of America
- Department of Computer Science, Arkansas State University, Jonesboro, Arkansas, United States of America
| | - Jaclyn Calkins
- Department of Biology, Reed College, Portland, Oregon, United States of America
- College of Human Medicine, Michigan State University, East Lansing, Michigan, United States of America
| | - Jeffrey R. Adrion
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, United States of America
| | - Leigh C. Latta IV
- Department of Biology, Reed College, Portland, Oregon, United States of America
- Lewis-Clark State College, Lewiston, Idaho, United States of America
| | - Sarah Schaack
- Department of Biology, Reed College, Portland, Oregon, United States of America
| |
Collapse
|
3
|
Jaron KS, Bast J, Nowell RW, Ranallo-Benavidez TR, Robinson-Rechavi M, Schwander T. Genomic Features of Parthenogenetic Animals. J Hered 2021; 112:19-33. [PMID: 32985658 PMCID: PMC7953838 DOI: 10.1093/jhered/esaa031] [Citation(s) in RCA: 35] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2020] [Accepted: 08/17/2020] [Indexed: 12/21/2022] Open
Abstract
Evolution without sex is predicted to impact genomes in numerous ways. Case studies of individual parthenogenetic animals have reported peculiar genomic features that were suggested to be caused by their mode of reproduction, including high heterozygosity, a high abundance of horizontally acquired genes, a low transposable element load, or the presence of palindromes. We systematically characterized these genomic features in published genomes of 26 parthenogenetic animals representing at least 18 independent transitions to asexuality. Surprisingly, not a single feature was systematically replicated across a majority of these transitions, suggesting that previously reported patterns were lineage-specific rather than illustrating the general consequences of parthenogenesis. We found that only parthenogens of hybrid origin were characterized by high heterozygosity levels. Parthenogens that were not of hybrid origin appeared to be largely homozygous, independent of the cellular mechanism underlying parthenogenesis. Overall, despite the importance of recombination rate variation for the evolution of sexual animal genomes, the genome-wide absence of recombination does not appear to have had the dramatic effects which are expected from classical theoretical models. The reasons for this are probably a combination of lineage-specific patterns, the impact of the origin of parthenogenesis, and a survivorship bias of parthenogenetic lineages.
Collapse
Affiliation(s)
- Kamil S Jaron
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Jens Bast
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| | - Reuben W Nowell
- Department of Life Sciences, Imperial College London, Ascot, Berkshire, UK
- Reuben W. Nowell is now at the Department of Zoology, University of Oxford, Oxford, UK
| | | | - Marc Robinson-Rechavi
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Tanja Schwander
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| |
Collapse
|
4
|
A Survey of Transposon Landscapes in the Putative Ancient Asexual Ostracod Darwinula stevensoni. Genes (Basel) 2021; 12:genes12030401. [PMID: 33799706 PMCID: PMC7998251 DOI: 10.3390/genes12030401] [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: 02/10/2021] [Revised: 03/02/2021] [Accepted: 03/06/2021] [Indexed: 11/17/2022] Open
Abstract
How asexual reproduction shapes transposable element (TE) content and diversity in eukaryotic genomes remains debated. We performed an initial survey of TE load and diversity in the putative ancient asexual ostracod Darwinula stevensoni. We examined long contiguous stretches of DNA in clones from a genomic fosmid library, totaling about 2.5 Mb, and supplemented these data with results on TE abundance and diversity from an Illumina draft genome. In contrast to other TE studies in putatively ancient asexuals, which revealed relatively low TE content, we found that at least 19% of the fosmid dataset and 26% of the genome assembly corresponded to known transposons. We observed a high diversity of transposon families, including LINE, gypsy, PLE, mariner/Tc, hAT, CMC, Sola2, Ginger, Merlin, Harbinger, MITEs and helitrons, with the prevalence of DNA transposons. The predominantly low levels of sequence diversity indicate that many TEs are or have recently been active. In the fosmid data, no correlation was found between telomeric repeats and non-LTR retrotransposons, which are present near telomeres in other taxa. Most TEs in the fosmid data were located outside of introns and almost none were found in exons. We also report an N-terminal Myb/SANT-like DNA-binding domain in site-specific R4/Dong non-LTR retrotransposons. Although initial results on transposable loads need to be verified with high quality draft genomes, this study provides important first insights into TE dynamics in putative ancient asexual ostracods.
Collapse
|
5
|
Bast J, Jaron KS, Schuseil D, Roze D, Schwander T. Asexual reproduction reduces transposable element load in experimental yeast populations. eLife 2019; 8:48548. [PMID: 31486772 PMCID: PMC6783261 DOI: 10.7554/elife.48548] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2019] [Accepted: 09/04/2019] [Indexed: 01/04/2023] Open
Abstract
Theory predicts that sexual reproduction can either facilitate or restrain transposable element (TE) accumulation by providing TEs with a means of spreading to all individuals in a population, versus facilitating TE load reduction via purifying selection. By quantifying genomic TE loads over time in experimental sexual and asexual Saccharomyces cerevisiae populations, we provide direct evidence that TE loads decrease rapidly under asexual reproduction. We show, using simulations, that this reduction may occur via evolution of TE activity, most likely via increased excision rates. Thus, sex is a major driver of genomic TE loads and at the root of the success of TEs. The genetic information of most living organisms contains parasitic invaders known as transposable elements. These genetic sequences multiply by copying and pasting themselves through the genome, but this process can disrupt the activity of important genes and put the organism at risk. How transposable elements proliferate in a population depends on the way organisms reproduce. If they simply clone themselves asexually, the selfish elements cannot spread between the different clones. If the organisms mate together their respective transposable elements get mixed, which helps the sequences to spread more easily and to potentially become more virulent. However, sexual reproduction also comes with mechanisms that keep transposable elements in check. Bast, Jaron et al. took advantage of the fact that yeasts can reproduce with or without mating to explore whether sexual or asexual organisms are better at controlling the spread of transposable elements. The number of copies of transposable elements in the genomes of yeast grown sexually or asexually was assessed. The results showed that sexual populations kept constant numbers of selfish elements, while asexual organisms lost these genomic parasites over time. Simulations then revealed that this difference emerged because a defense gene that helps to delete transposable elements was spreading more quickly in the asexual group. The work by Bast, Jaron et al. therefore suggests that sex is responsible for the evolutionary success of transposable elements, while asexual populations can discard these sequences over time. Sex therefore helps genetic parasites, somewhat similar to sexually transmitted diseases, to spread between individuals and remain virulent.
Collapse
Affiliation(s)
- Jens Bast
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| | - Kamil S Jaron
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Donovan Schuseil
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| | - Denis Roze
- UMI3614 Evolutionary Biology and Ecology of Alga, Centre National de la Recherche Scientifique (CNRS), Roscoff, France.,Station Biologique de Roscoff, Sorbonne Université, Roscoff, France
| | - Tanja Schwander
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| |
Collapse
|
6
|
Abstract
Daphnia normally reproduce by cyclical parthenogenesis, with offspring sex being determined by environmental cues. However, some females have lost the ability to produce males. Our results demonstrate that this loss of male-producing ability is controlled by a dominant allele at a single locus. We identified the locus by comparing whole-genome sequences of 67 nonmale-producing (NMP) and 100 male-producing (MP) clones from 5 Daphnia pulex populations, revealing 132 NMP-linked SNPs and 59 NMP-linked indels within a single 1.1-Mb nonrecombining region on chromosome I. These markers include 7 nonsynonymous mutations, all of which are located within one unannotated protein-coding gene (gene 8960). Within this single gene, all of the marker-linked NMP haplotypes from different populations form a monophyletic clade, suggesting a single origin of the NMP phenotype, with the NMP haplotype originating by introgression from a sister species, Daphnia pulicaria Methyl farnesoate (MF) is the innate juvenile hormone in daphnids, which induces the production of males and whose inhibition results in female-only production. Gene 8960 is sensitive to treatment by MF in MP clones, but such responsiveness is greatly reduced in NMP clones. Thus, we hypothesize that gene 8960 is located downstream of the MF-signaling pathway in D. pulex, with the NMP phenotype being caused by expression change of gene 8960.
Collapse
|
7
|
Jiang X, Tang H, Mohammed Ismail W, Lynch M. A Maximum-Likelihood Approach to Estimating the Insertion Frequencies of Transposable Elements from Population Sequencing Data. Mol Biol Evol 2018; 35:2560-2571. [PMID: 30099533 PMCID: PMC6188571 DOI: 10.1093/molbev/msy152] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
Transposable elements (TEs) contribute to a large fraction of the expansion of many eukaryotic genomes due to the capability of TEs duplicating themselves through transposition. A first step to understanding the roles of TEs in a eukaryotic genome is to characterize the population-wide variation of TE insertions in the species. Here, we present a maximum-likelihood (ML) method for estimating allele frequencies and detecting selection on TE insertions in a diploid population, based on the genotypes at TE insertion sites detected in multiple individuals sampled from the population using paired-end (PE) sequencing reads. Tests of the method on simulated data show that it can accurately estimate the allele frequencies of TE insertions even when the PE sequencing is conducted at a relatively low coverage (=5X). The method can also detect TE insertions under strong selection, and the detection ability increases with sample size in a population, although a substantial fraction of actual TE insertions under selection may be undetected. Application of the ML method to genomic sequencing data collected from a natural Daphnia pulex population shows that, on the one hand, most (>90%) TE insertions present in the reference D. pulex genome are either fixed or nearly fixed (with allele frequencies >0.95); on the other hand, among the nonreference TE insertions (i.e., those detected in some individuals in the population but absent from the reference genome), the majority (>70%) are still at low frequencies (<0.1). Finally, we detected a substantial fraction (∼9%) of nonreference TE insertions under selection.
Collapse
Affiliation(s)
- Xiaoqian Jiang
- Department of Biology, Indiana University, Bloomington, IN
| | - Haixu Tang
- School of Informatics and Computing, Indiana University, Bloomington, IN
| | | | - Michael Lynch
- Center for Mechanisms of Evolution, Arizona State University, Temple, AZ
| |
Collapse
|
8
|
Abstract
Comparing genomes of closely related genotypes from populations with distinct demographic histories can help reveal the impact of effective population size on genome evolution. For this purpose, we present a high quality genome assembly of Daphnia pulex (PA42), and compare this with the first sequenced genome of this species (TCO), which was derived from an isolate from a population with >90% reduction in nucleotide diversity. PA42 has numerous similarities to TCO at the gene level, with an average amino acid sequence identity of 98.8 and >60% of orthologous proteins identical. Nonetheless, there is a highly elevated number of genes in the TCO genome annotation, with ∼7000 excess genes appearing to be false positives. This view is supported by the high GC content, lack of introns, and short length of these suspicious gene annotations. Consistent with the view that reduced effective population size can facilitate the accumulation of slightly deleterious genomic features, we observe more proliferation of transposable elements (TEs) and a higher frequency of gained introns in the TCO genome.
Collapse
|