1
|
Comparative analysis reveals within-population genome size variation in a rotifer is driven by large genomic elements with highly abundant satellite DNA repeat elements. BMC Biol 2021; 19:206. [PMID: 34530817 PMCID: PMC8447722 DOI: 10.1186/s12915-021-01134-w] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Accepted: 08/27/2021] [Indexed: 12/02/2022] Open
Abstract
BACKGROUND Eukaryotic genomes are known to display an enormous variation in size, but the evolutionary causes of this phenomenon are still poorly understood. To obtain mechanistic insights into such variation, previous studies have often employed comparative genomics approaches involving closely related species or geographically isolated populations within a species. Genome comparisons among individuals of the same population remained so far understudied-despite their great potential in providing a microevolutionary perspective to genome size evolution. The rotifer Brachionus asplanchnoidis represents one of the most extreme cases of within-population genome size variation among eukaryotes, displaying almost twofold variation within a geographic population. RESULTS Here, we used a whole-genome sequencing approach to identify the underlying DNA sequence differences by assembling a high-quality reference genome draft for one individual of the population and aligning short reads of 15 individuals from the same geographic population including the reference individual. We identified several large, contiguous copy number variable regions (CNVs), up to megabases in size, which exhibited striking coverage differences among individuals, and whose coverage overall scaled with genome size. CNVs were of remarkably low complexity, being mainly composed of tandemly repeated satellite DNA with only a few interspersed genes or other sequences, and were characterized by a significantly elevated GC-content. CNV patterns in offspring of two parents with divergent genome size and CNV patterns in several individuals from an inbred line differing in genome size demonstrated inheritance and accumulation of CNVs across generations. CONCLUSIONS By identifying the exact genomic elements that cause within-population genome size variation, our study paves the way for studying genome size evolution in contemporary populations rather than inferring patterns and processes a posteriori from species comparisons.
Collapse
|
2
|
Substantial intraspecific genome size variation in golden-brown algae and its phenotypic consequences. ANNALS OF BOTANY 2020; 126:1077-1087. [PMID: 32686820 PMCID: PMC7596369 DOI: 10.1093/aob/mcaa133] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/13/2020] [Accepted: 07/17/2020] [Indexed: 06/11/2023]
Abstract
BACKGROUND AND AIMS While nuclear DNA content variation and its phenotypic consequences have been well described for animals, vascular plants and macroalgae, much less about this topic is known regarding unicellular algae and protists in general. The dearth of data is especially pronounced when it comes to intraspecific genome size variation. This study attempts to investigate the extent of intraspecific variability in genome size and its adaptive consequences in a microalgal species. METHODS Propidium iodide flow cytometry was used to estimate the absolute genome size of 131 strains (isolates) of the golden-brown alga Synura petersenii (Chrysophyceae, Stramenopiles), identified by identical internal transcribed spacer (ITS) rDNA barcodes. Cell size, growth rate and genomic GC content were further assessed on a sub-set of strains. Geographic location of 67 sampling sites across the Northern hemisphere was used to extract climatic database data and to evaluate the ecogeographical distribution of genome size diversity. KEY RESULTS Genome size ranged continuously from 0.97 to 2.02 pg of DNA across the investigated strains. The genome size was positively associated with cell size and negatively associated with growth rate. Bioclim variables were not correlated with genome size variation. No clear trends in the geographical distribution of strains of a particular genome size were detected, and strains of different genome size occasionally coexisted at the same locality. Genomic GC content was significantly associated only with genome size via a quadratic relationship. CONCLUSIONS Genome size variability in S. petersenii was probably triggered by an evolutionary mechanism operating via gradual changes in genome size accompanied by changes in genomic GC content, such as, for example, proliferation of transposable elements. The variation was reflected in cell size and relative growth rate, possibly with adaptive consequences.
Collapse
|
3
|
Within-Population Genome Size Variation is Mediated by Multiple Genomic Elements That Segregate Independently during Meiosis. Genome Biol Evol 2019; 11:3424-3435. [PMID: 31742335 PMCID: PMC7145553 DOI: 10.1093/gbe/evz253] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/14/2019] [Indexed: 01/14/2023] Open
Abstract
Within-species variation in genome size has been documented in many animals and plants. Despite its importance for understanding eukaryotic genome diversity, there is only sparse knowledge about how individual-level processes mediate genome size variation in populations. Here, we study a natural population of the rotifer Brachionus asplanchnoidis whose members differ up to 1.9-fold in diploid genome size, but were still able to interbreed and produce viable offspring. We show that genome size is highly heritable and can be artificially selected up or down, but not below a certain basal diploid genome size for this species. Analyses of segregation patterns in haploid males reveal that large genomic elements (several megabases in size) provide the substrate of genome size variation. These elements, and their segregation patterns, explain the generation of new genome size variants, the short-term evolutionary potential of genome size change in populations, and some seemingly paradoxical patterns, like an increase in genome size variation among highly inbred lines. Our study suggests that a conceptual model involving only two variables, 1) a basal genome size of the population, and 2) a vector containing information on additional elements that may increase genome size in this population (size, number, and meiotic segregation behavior), can effectively address most scenarios of short-term evolutionary change of genome size in a population.
Collapse
|
4
|
Small, but surprisingly repetitive genomes: transposon expansion and not polyploidy has driven a doubling in genome size in a metazoan species complex. BMC Genomics 2019; 20:466. [PMID: 31174483 PMCID: PMC6555955 DOI: 10.1186/s12864-019-5859-y] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2019] [Accepted: 05/29/2019] [Indexed: 02/01/2023] Open
Abstract
BACKGROUND The causes and consequences of genome size variation across Eukaryotes, which spans five orders of magnitude, have been hotly debated since before the advent of genome sequencing. Previous studies have mostly examined variation among larger taxonomic units (e.g., orders, or genera), while comparisons among closely related species are rare. Rotifers of the Brachionus plicatilis species complex exhibit a seven-fold variation in genome size and thus represent a unique opportunity to study such changes on a relatively short evolutionary timescale. Here, we sequenced and analysed the genomes of four species of this complex with nuclear DNA contents spanning 110-422 Mbp. To establish the likely mechanisms of genome size change, we analysed both sequencing read libraries and assemblies for signatures of polyploidy and repetitive element content. We also compared these genomes to that of B. calyciflorus, the closest relative with a sequenced genome (293 Mbp nuclear DNA content). RESULTS Despite the very large differences in genome size, we saw no evidence of ploidy level changes across the B. plicatilis complex. However, repetitive element content explained a large portion of genome size variation (at least 54%). The species with the largest genome, B. asplanchnoidis, has a strikingly high 44% repetitive element content, while the smaller B. plicatilis genomes contain between 14 and 25% repetitive elements. According to our analyses, the B. calyciflorus genome contains 39% repetitive elements, which is substantially higher than previously reported (21%), and suggests that high repetitive element load could be widespread in monogonont rotifers. CONCLUSIONS Even though the genome sizes of these species are at the low end of the metazoan spectrum, their genomes contain substantial amounts of repetitive elements. Polyploidy does not appear to play a role in genome size variations in these species, and these variations can be mostly explained by changes in repetitive element content. This contradicts the naïve expectation that small genomes are streamlined, or less complex, and that large variations in nuclear DNA content between closely related species are due to polyploidy.
Collapse
|
5
|
Gene expression changes elicited by a parasitic B chromosome in the grasshopper Eyprepocnemis plorans are consistent with its phenotypic effects. Chromosoma 2019; 128:53-67. [PMID: 30617552 DOI: 10.1007/s00412-018-00689-y] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2018] [Revised: 12/20/2018] [Accepted: 12/21/2018] [Indexed: 01/16/2023]
Abstract
Parasitism evokes adaptive physiological changes in the host, many of which take place through gene expression changes. This response can be more or less local, depending on the organ or tissue affected by the parasite, or else systemic when the parasite affects the entire host body. The most extreme of the latter cases is intragenomic parasitism, where the parasite is present in all host nuclei as any other genomic element. Here, we show the molecular crosstalk between a parasitic chromosome (also named B chromosome) and the host genome, manifested through gene expression changes. The transcriptome analysis of 0B and 1B females of the grasshopper Eyprepocnemis plorans, validated by a microarray experiment performed on four B-lacking and five B-carrying females, revealed changes in gene expression for 188 unigenes being consistent in both experiments. Once discarded B-derived transcripts, there were 46 differentially expressed genes (30 up- and 16 downregulated) related with the adaptation of the host genome to the presence of the parasitic chromosome. Interestingly, the functions of these genes could explain some of the most important effects of B chromosomes, such as nucleotypic effects derived from the additional DNA they represent, chemical defense and detoxification, protein modification and response to stress, ovary function, and regulation of gene expression. Collectively, these changes uncover an intimate host-parasite interaction between A and B chromosomes during crucial steps of gene expression and protein function.
Collapse
|
6
|
Accumulation of Transposable Elements in Autosomes and Giant Sex Chromosomes of Omophoita (Chrysomelidae: Alticinae). Cytogenet Genome Res 2018; 156:215-222. [PMID: 30504708 DOI: 10.1159/000495199] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/07/2018] [Indexed: 12/28/2022] Open
Abstract
Coleoptera is the most diverse order among insects, and comparative molecular cytogenetic studies in this group are lacking. The species of Omophoita (Oedionychina) possess a karyotype of 2n = 22 = 10II+X+Y. They are interesting models for evolutionary cytogenetic studies due to giant sex chromosomes which are asynaptic during meiosis. Transposable elements (TEs) confer plasticity and mobility to genomes and are considered hotspots for DNA double-strand breaks and chromosomal rearrangements. The objective of the present study was to verify the role of TEs in the karyotype and in the size expansion of the giant sex chromosomes in Omophoita. Thus, different TEs were characterized in the Omophoita genome and localized in the chromosomes by fluorescence in situ hybridization (FISH). The DNA sequencing data revealed identity with TE families Tc1/Mariner and RTE/L1-56_XT. FISH showed signals of all TEs in the karyotypes and a high accumulation in the sex chromosomes of the 3 Omophoita species analyzed. These data suggest that the genome size expansion and the origin of the giant sex chromosomes of Omophoita are due to an intensive genomic invasion of TEs, as those characterized here as Tc1/Mariner-Ooc and RTE-Ooc. Differences in the chromosomal location of the TEs among the 3 species indicate that they have participated in the karyotype differentiation in Omophoita.
Collapse
|
7
|
High-throughput analysis of satellite DNA in the grasshopper Pyrgomorpha conica reveals abundance of homologous and heterologous higher-order repeats. Chromosoma 2018; 127:323-340. [PMID: 29549528 DOI: 10.1007/s00412-018-0666-9] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2017] [Revised: 02/13/2018] [Accepted: 03/06/2018] [Indexed: 12/25/2022]
Abstract
Satellite DNA (satDNA) constitutes an important fraction of repetitive DNA in eukaryotic genomes, but it is barely known in most species. The high-throughput analysis of satDNA in the grasshopper Pyrgomorpha conica revealed 87 satDNA variants grouped into 76 different families, representing 9.4% of the genome. Fluorescent in situ hybridization (FISH) analysis of the 38 most abundant satDNA families revealed four different patterns of chromosome distribution. Homology search between the 76 satDNA families showed the existence of 15 superfamilies, each including two or more families, with the most abundant superfamily representing more than 80% of all satDNA found in this species. This also revealed the presence of two types of higher-order repeats (HORs), one showing internal homologous subrepeats, as conventional HORs, and an additional type showing non-homologous internal subrepeats, the latter arising by the combination of a given satDNA family with a non-annotated sequence, or with telomeric DNA. Interestingly, the heterologous subrepeats included in these HORs showed higher divergence within the HOR than outside it, suggesting that heterologous HORs show poor homogenization, in high contrast with conventional (homologous) HORs. Finally, heterologous HORs can show high differences in divergence between their constituent subrepeats, suggesting the possibility of regional homogenization.
Collapse
|
8
|
Quantitative sequence characterization for repetitive DNA content in the supernumerary chromosome of the migratory locust. Chromosoma 2017; 127:45-57. [PMID: 28868580 DOI: 10.1007/s00412-017-0644-7] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2017] [Revised: 08/23/2017] [Accepted: 08/24/2017] [Indexed: 12/23/2022]
Abstract
Repetitive DNA is a major component in most eukaryotic genomes but is ignored in most genome sequencing projects. Here, we report the quantitative composition in repetitive DNA for a supernumerary (B) chromosome, in the migratory locust (Locusta migratoria), by Illumina sequencing of genomic DNA from B-carrying and B-lacking individuals and DNA obtained from a microdissected B chromosome, as well as the physical mapping of some elements. B chromosome DNA of 94.9% was repetitive, in high contrast with the 64.1% of standard (A) chromosomes. B chromosomes are enriched in satellite DNA (satDNA) (65.2% of B-DNA), with a single satellite (LmiSat02-176) comprising 55% of the B. Six satDNAs were visualized by FISH on the B chromosome, and the only A chromosome carrying all these satellites was autosome 9, pointing to this chromosome, along with autosome 8 (which shares histone genes with the B) as putative ancestors of the B chromosome. We found several transposable elements (TEs) showing nucleotidic variation specific to B-carrying individuals, which was also present in B-carrying transcriptomes. Remarkably, an interstitial region of the B chromosome included a 17 kb chimera composed of 29 different TEs, suggesting reiterative TE insertion in this B chromosome region.
Collapse
|
9
|
Protein-coding genes in B chromosomes of the grasshopper Eyprepocnemis plorans. Sci Rep 2017; 7:45200. [PMID: 28367986 PMCID: PMC5377258 DOI: 10.1038/srep45200] [Citation(s) in RCA: 50] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2016] [Accepted: 02/22/2017] [Indexed: 01/20/2023] Open
Abstract
For many years, parasitic B chromosomes have been considered genetically inert elements. Here we show the presence of ten protein-coding genes in the B chromosome of the grasshopper Eyprepocnemis plorans. Four of these genes (CIP2A, GTPB6, KIF20A, and MTG1) were complete in the B chromosome whereas the six remaining (CKAP2, CAP-G, HYI, MYCB2, SLIT and TOP2A) were truncated. Five of these genes (CIP2A, CKAP2, CAP-G, KIF20A, and MYCB2) were significantly up-regulated in B-carrying individuals, as expected if they were actively transcribed from the B chromosome. This conclusion is supported by three truncated genes (CKAP2, CAP-G and MYCB2) which showed up-regulation only in the regions being present in the B chromosome. Our results indicate that B chromosomes are not so silenced as was hitherto believed. Interestingly, the five active genes in the B chromosome code for functions related with cell division, which is the main arena where B chromosome destiny is played. This suggests that B chromosome evolutionary success can lie on its gene content.
Collapse
|
10
|
Satellite DNA content illuminates the ancestry of a supernumerary (B) chromosome. Chromosoma 2016; 126:487-500. [PMID: 27522227 DOI: 10.1007/s00412-016-0611-8] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2016] [Revised: 07/28/2016] [Accepted: 08/02/2016] [Indexed: 12/23/2022]
Abstract
B chromosomes are supernumerary genomic elements most likely derived from the standard (A) chromosomes, whose dispensability has freed their DNA sequences to evolve fast, thus making it difficult to uncover their ancestry. Here, we show the ancestry of a B chromosome in the grasshopper Eumigus monticola by means of the high-throughput analysis of the satellitome, i.e., the whole collection of satellite DNA (satDNA). The satellitome found in this species consists of 27 satDNA families, with monomer length between 5 and 325 nt and A + T content between 42.9 and 83.3 %. Two out of the 20 clustered satDNA families (EmoSat26-41 and EmoSat27-102) were observed only on the B chromosome. The A chromosome carrying the highest number of satDNA families was the megameric S8 (13 families), six of which were also present in the B chromosome, and three of these were exclusive of the S8 and B chromosomes. The absence in the B chromosome of the H3 histone gene cluster (located interstitially on S8) and three satDNA families (located distally on S8) allowed delimiting the possible origin of the B chromosome to the proximal third of the S8 autosome, through a breakpoint between EmoSat11-122 and the H3 cluster. Interestingly, bioinformatic analysis revealed the presence of seeds for the two B-specific satDNAs in the A chromosomes, suggesting their massive amplification in the B chromosome after its origin. Therefore, intraspecifically arisen B chromosomes can harbor DNA sequences apparently being B-specific.
Collapse
|
11
|
Comment on Schielzeth et al. (2014): "Genome size variation affects song attractiveness in grasshoppers: Evidence for sexual selection against large genomes". Evolution 2016; 70:1428-30. [PMID: 27327141 DOI: 10.1111/evo.12869] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2015] [Accepted: 01/13/2016] [Indexed: 11/29/2022]
Abstract
Schielzeth et al. (2014) concluded that attractive grasshopper singers have significantly smaller genomes thus suggesting a possible role for sexual selection on genome size. Whereas this conclusion could still be conceivably valid, it is not supported by the data presented due to some technical flaws. In addition, the interpretation of the results, speculating on the possible presence of B chromosomes, is not justified.
Collapse
|
12
|
Non-random expression of ribosomal DNA units in a grasshopper showing high intragenomic variation for the ITS2 region. INSECT MOLECULAR BIOLOGY 2015; 24:319-330. [PMID: 25565136 DOI: 10.1111/imb.12158] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We analyse intragenomic variation of the ITS2 internal transcribed spacer of ribosomal DNA (rDNA) in the grasshopper Eyprepocnemis plorans, by means of tagged PCR 454 amplicon sequencing performed on both genomic DNA (gDNA) and RNA-derived complementary DNA (cDNA), using part of the ITS2 flanking coding regions (5.8S and 28S rDNA) as an internal control for sequencing errors. Six different ITS2 haplotypes (i.e. variants for at least one nucleotide in the complete ITS2 sequence) were found in a single population, one of them (Hap4) being specific to a supernumerary (B) chromosome. The analysis of both gDNA and cDNA from the same individuals provided an estimate of the expression efficiency of the different haplotypes. We found random expression (i.e. about similar recovery in gDNA and cDNA) for three haplotypes (Hap1, Hap2 and Hap5), but significant underexpression for three others (Hap3, Hap4 and Hap6). Hap4 was the most extremely underexpressed and, remarkably, it showed the lowest sequence conservation for the flanking 5.8-28S coding regions in the gDNA reads but the highest conservation (100%) in the cDNA ones, suggesting the preferential expression of mutation-free rDNA units carrying this ITS2 haplotype. These results indicate that the ITS2 region of rDNA is far from complete homogenization in this species, and that the different rDNA units are not expressed at random, with some of them being severely downregulated.
Collapse
|
13
|
Intragenomic distribution of RTE retroelements suggests intrachromosomal movement. Chromosome Res 2015; 23:211-23. [PMID: 25605325 DOI: 10.1007/s10577-014-9461-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2014] [Revised: 11/25/2014] [Accepted: 12/18/2014] [Indexed: 11/25/2022]
Abstract
Much is known about the abundance of transposable elements (TEs) in eukaryotic genomes, but much is still unknown on their behaviour within cells. We employ here a combination of cytological, molecular and genomic approaches providing information on the intragenomic distribution and behaviour of non-long terminal repeat (LTR) retrotransposon-like elements (RTE). We microdissected every chromosome in a single first meiotic metaphase cell of the grasshopper Eyprepocnemis plorans and polymerase chain reaction (PCR) amplified a fragment of the RTE reverse transcriptase gene with specific primers. PCR products were cloned and 139 clones were sequenced. Analysis of molecular variance (AMOVA) showed significant intragenomic structure for these elements, with 4.6 % of molecular variance being found between chromosomes. A maximum likelihood tree built with the RTE sequences revealed the frequent presence of two or more elements showing very high similarity and being located on the same chromosome, thus suggesting intrachromosome movement. The 454 pyrosequencing of genomic DNA gave strong support to the microdissection results and provided evidence for the existence of 5' truncated elements. Our results thus indicate a tendency of RTE elements to reinsert into the same chromosome from where they were transcribed, which could be achieved if retrotranscription and insertion takes place immediately after transcription.
Collapse
|
14
|
Abstract
Cactophilic Drosophila species provide a valuable model to study gene–environment interactions and ecological adaptation. Drosophila buzzatii and Drosophila mojavensis are two cactophilic species that belong to the repleta group, but have very different geographical distributions and primary host plants. To investigate the genomic basis of ecological adaptation, we sequenced the genome and developmental transcriptome of D. buzzatii and compared its gene content with that of D. mojavensis and two other noncactophilic Drosophila species in the same subgenus. The newly sequenced D. buzzatii genome (161.5 Mb) comprises 826 scaffolds (>3 kb) and contains 13,657 annotated protein-coding genes. Using RNA sequencing data of five life-stages we found expression of 15,026 genes, 80% protein-coding genes, and 20% noncoding RNA genes. In total, we detected 1,294 genes putatively under positive selection. Interestingly, among genes under positive selection in the D. mojavensis lineage, there is an excess of genes involved in metabolism of heterocyclic compounds that are abundant in Stenocereus cacti and toxic to nonresident Drosophila species. We found 117 orphan genes in the shared D. buzzatii–D. mojavensis lineage. In addition, gene duplication analysis identified lineage-specific expanded families with functional annotations associated with proteolysis, zinc ion binding, chitin binding, sensory perception, ethanol tolerance, immunity, physiology, and reproduction. In summary, we identified genetic signatures of adaptation in the shared D. buzzatii–D. mojavensis lineage, and in the two separate D. buzzatii and D. mojavensis lineages. Many of the novel lineage-specific genomic features are promising candidates for explaining the adaptation of these species to their distinct ecological niches.
Collapse
|
15
|
A step to the gigantic genome of the desert locust: chromosome sizes and repeated DNAs. Chromosoma 2014; 124:263-75. [PMID: 25472934 DOI: 10.1007/s00412-014-0499-0] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2014] [Revised: 11/24/2014] [Accepted: 11/25/2014] [Indexed: 10/24/2022]
Abstract
The desert locust (Schistocerca gregaria) has been used as material for numerous cytogenetic studies. Its genome size is estimated to be 8.55 Gb of DNA comprised in 11 autosomes and the X chromosome. Its X0/XX sex chromosome determinism therefore results in females having 24 chromosomes whereas males have 23. Surprisingly, little is known about the DNA content of this locust's huge chromosomes. Here, we use the Feulgen Image Analysis Densitometry and C-banding techniques to respectively estimate the DNA quantity and heterochromatin content of each chromosome. We also identify three satellite DNAs using both restriction endonucleases and next-generation sequencing. We then use fluorescent in situ hybridization to determine the chromosomal location of these satellite DNAs as well as that of six tandem repeat DNA gene families. The combination of the results obtained in this work allows distinguishing between the different chromosomes not only by size, but also by the kind of repetitive DNAs that they contain. The recent publication of the draft genome of the migratory locust (Locusta migratoria), the largest animal genome hitherto sequenced, invites for sequencing even larger genomes. S. gregaria is a pest that causes high economic losses. It is thus among the primary candidates for genome sequencing. But this species genome is about 50 % larger than that of L. migratoria, and although next-generation sequencing currently allows sequencing large genomes, sequencing it would mean a greater challenge. The chromosome sizes and markers provided here should not only help planning the sequencing project and guide the assembly but would also facilitate assigning assembled linkage groups to actual chromosomes.
Collapse
|
16
|
Next generation sequencing and FISH reveal uneven and nonrandom microsatellite distribution in two grasshopper genomes. Chromosoma 2014; 124:221-34. [PMID: 25387401 DOI: 10.1007/s00412-014-0492-7] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2014] [Revised: 10/17/2014] [Accepted: 10/29/2014] [Indexed: 11/29/2022]
Abstract
Simple sequence repeats (SSRs), also known as microsatellites, are one of the prominent DNA sequences shaping the repeated fraction of eukaryotic genomes. In spite of their profuse use as molecular markers for a variety of genetic and evolutionary studies, their genomic location, distribution, and function are not yet well understood. Here we report the first thorough joint analysis of microsatellite motifs at both genomic and chromosomal levels in animal species, by a combination of 454 sequencing and fluorescent in situ hybridization (FISH) techniques performed on two grasshopper species. The in silico analysis of the 454 reads suggested that microsatellite expansion is not driving size increase of these genomes, as SSR abundance was higher in the species showing the smallest genome. However, the two species showed the same uneven and nonrandom location of SSRs, with clear predominance of dinucleotide motifs and association with several types of repetitive elements, mostly histone gene spacers, ribosomal DNA intergenic spacers (IGS), and transposable elements (TEs). The FISH analysis showed a dispersed chromosome distribution of microsatellite motifs in euchromatic regions, in coincidence with chromosome location patterns previously observed for many mobile elements in these species. However, some SSR motifs were clustered, especially those located in the histone gene cluster.
Collapse
|
17
|
Preferential occupancy of R2 retroelements on the B chromosomes of the grasshopper Eyprepocnemis plorans. PLoS One 2014; 9:e91820. [PMID: 24632855 PMCID: PMC3954772 DOI: 10.1371/journal.pone.0091820] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2014] [Accepted: 02/14/2014] [Indexed: 02/02/2023] Open
Abstract
R2 non-LTR retrotransposons exclusively insert into the 28S rRNA genes of their host, and are expressed by co-transcription with the rDNA unit. The grasshopper Eyprepocnemis plorans contains transcribed rDNA clusters on most of its A chromosomes, as well as non-transcribed rDNA clusters on the parasitic B chromosomes found in many populations. Here the structure of the E. plorans R2 element, its abundance relative to the number of rDNA units and its retrotransposition activity were determined. Animals screened from five populations contained on average over 12,000 rDNA units on their A chromosomes, but surprisingly only about 100 R2 elements. Monitoring the patterns of R2 insertions in individuals from these populations revealed only low levels of retrotransposition. The low rates of R2 insertion observed in E. plorans differ from the high levels of R2 insertion previously observed in insect species that have many fewer rDNA units. It is proposed that high levels of R2 are strongly selected against in E. plorans, because the rDNA transcription machinery in this species is unable to differentiate between R2-inserted and uninserted units. The B chromosomes of E. plorans contain an additional 7,000 to 15,000 rDNA units, but in contrast to the A chromosomes, from 150 to over 1,500 R2 elements. The higher concentration of R2 in the inactive B chromosomes rDNA clusters suggests these chromosomes can act as a sink for R2 insertions thus further reducing the level of insertions on the A chromosomes. These studies suggest an interesting evolutionary relationship between the parasitic B chromosomes and R2 elements.
Collapse
|
18
|
HP1 knockdown is associated with abnormal condensation of almost all chromatin types in a grasshopper (Eyprepocnemis plorans). Chromosome Res 2014; 22:253-66. [PMID: 24398928 DOI: 10.1007/s10577-013-9399-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2013] [Revised: 12/11/2013] [Accepted: 12/12/2013] [Indexed: 10/25/2022]
Abstract
Heterochromatin protein 1 (HP1) is a highly conserved family of eukaryotic proteins required for heterochromatic gene silencing and euchromatic gene transcription regulation. In addition, HP1 is involved in chromatin organization and protection of chromosome integrity during cell division. Here, we present a cytological and molecular analysis of the effects of HP1 knockdown in Eyprepocnemis plorans, a grasshopper species polymorphic for supernumerary heterochromatic chromosomes. Our results revealed contrasting effects of HP1 knockdown on gene activity. While the Bub1 gene decreased in expression level in HP1 knockdown animals, NOR activity, rRNA and, contrarily to previous reports in Drosophila, Hsp70 gene expression remained unchanged. Furthermore, HP1 knockdown resulted in abnormal chromatin condensation, chromosomal bridges, higher frequency of macrospermatids, loss of muscle mass and hemolymph amount as well as a low number of dividing cells and survival reduction. All these phenotypes are very likely due to the chromatin condensation disruption observed for almost all kinds of chromatin.
Collapse
|
19
|
Disparate molecular evolution of two types of repetitive DNAs in the genome of the grasshopper Eyprepocnemis plorans. Heredity (Edinb) 2013; 112:531-42. [PMID: 24346496 DOI: 10.1038/hdy.2013.135] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2013] [Accepted: 11/14/2013] [Indexed: 11/09/2022] Open
Abstract
Wide arrays of repetitive DNA sequences form an important part of eukaryotic genomes. These repeats appear to evolve as coherent families, where repeats within a family are more similar to each other than to other orthologous representatives in related species. The continuous homogenization of repeats, through selective and non-selective processes, is termed concerted evolution. Ascertaining the level of variation between repeats is crucial to determining which evolutionary model best explains the homogenization observed for these sequences. Here, for the grasshopper Eyprepocnemis plorans, we present the analysis of intragenomic diversity for two repetitive DNA sequences (a satellite DNA (satDNA) and the 45S rDNA) resulting from the independent microdissection of several chromosomes. Our results show different homogenization patterns for these two kinds of paralogous DNA sequences, with a high between-chromosome structure for rDNA but no structure at all for the satDNA. This difference is puzzling, considering the adjacent localization of the two repetitive DNAs on paracentromeric regions in most chromosomes. The disparate homogenization patterns detected for these two repetitive DNA sequences suggest that several processes participate in the concerted evolution in E. plorans, and that these mechanisms might not work as genome-wide processes but rather as sequence-specific ones.
Collapse
|
20
|
Ribosomal DNA is active in different B chromosome variants of the grasshopper Eyprepocnemis plorans. Genetica 2013; 141:337-45. [PMID: 24008810 DOI: 10.1007/s10709-013-9733-6] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2013] [Accepted: 08/31/2013] [Indexed: 11/26/2022]
Abstract
B chromosomes are considered to be genetically inert elements. However, some of them are able to show nucleolus organizer region (NOR) activity, as detected by both cytological and molecular means. The grasshopper Eyprepocnemis plorans shows a B chromosome polymorphism characterized by the existence of many B variants. One of them, B24, shows NOR activity in about half of B-carrying males in the Torrox population. Molecular data have suggested the recent origin for B chromosomes in this species, and on this basis it would be expected that NOR activity was widespread among the different B variants. Here we test this hypothesis in four different B chromosome variants (B1, B2, B5, and B24) from 11 natural populations of the grasshopper E. plorans covering the south and east of the Iberian Peninsula plus the Balearic Islands. We used two different approaches: (1) the cytological observation of nucleoli attached to the distal region of the B chromosome (where the rDNA is located), and (2) the molecular detection of the rDNA transcripts carrying an adenine insertion characteristic of B chromosome ITS2 sequences. The results showed NOR expression not only for B24 but also for the B1 and B2 variants. However, the level of B-NOR expression in these latter variants, measured by the proportion of cells showing nucleoli attached to the B chromosomes, was much lower than that previously reported for B24. This suggests the possibility that structural or genetic background conditions are enhancing the expressivity of the rDNA in the B24 variant.
Collapse
|
21
|
Chromosomal Mapping of Repetitive DNAs in the Grasshopper Abracris flavolineata Reveal Possible Ancestry of the B Chromosome and H3 Histone Spreading. PLoS One 2013; 8:e66532. [PMID: 23826099 PMCID: PMC3694960 DOI: 10.1371/journal.pone.0066532] [Citation(s) in RCA: 75] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2013] [Accepted: 05/06/2013] [Indexed: 11/18/2022] Open
Abstract
Supernumerary chromosomes (B chromosomes) occur in approximately 15% of eukaryote species. Although these chromosomes have been extensively studied, knowledge concerning their specific molecular composition is lacking in most cases. The accumulation of repetitive DNAs is one remarkable characteristic of B chromosomes, and the occurrence of distinct types of multigene families, satellite DNAs and some transposable elements have been reported. Here, we describe the organization of repetitive DNAs in the A complement and B chromosome system in the grasshopper species Abracris flavolineata using classical cytogenetic techniques and FISH analysis using probes for five multigene families, telomeric repeats and repetitive C0t-1 DNA fractions. The 18S rRNA and H3 histone multigene families are highly variable and well distributed in A. flavolineata chromosomes, which contrasts with the conservation of U snRNA genes and less variable distribution of 5S rDNA sequences. The H3 histone gene was an extensively distributed with clusters occurring in all chromosomes. Repetitive DNAs were concentrated in C-positive regions, including the pericentromeric region and small chromosomal arms, with some occurrence in C-negative regions, but abundance was low in the B chromosome. Finally, the first demonstration of the U2 snRNA gene in B chromosomes in A. flavolineata may shed light on its possible origin. These results provide new information regarding chromosomal variability for repetitive DNAs in grasshoppers and the specific molecular composition of B chromosomes.
Collapse
|
22
|
The Ku70 DNA-repair protein is involved in centromere function in a grasshopper species. Chromosome Res 2013; 21:393-406. [PMID: 23797468 DOI: 10.1007/s10577-013-9367-7] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2013] [Revised: 06/07/2013] [Accepted: 06/09/2013] [Indexed: 01/05/2023]
Abstract
The Ku70 protein is involved in numerous cell functions, the nonhomologous end joining (NHEJ) DNA repair pathway being the best known. Here, we report a novel function for this protein in the grasshopper Eyprepocnemis plorans. We observed the presence of large Ku70 foci on the centromeres of meiotic and mitotic chromosomes during the cell cycle stages showing the highest centromeric activity (i.e., metaphase and anaphase). The fact that colchicine treatment prevented centromeric location of Ku70, suggests a microtubule-dependent centromeric function for Ku70. Likewise, the absence of Ku70 at metaphase-anaphase centromeres from three males whose Ku70 gene had been knocked down using interference RNA, and the dramatic increase in the frequency of polyploid spermatids observed in these males, suggest that the centromeric presence of Ku70 is required for normal cytokinesis in this species. The centromeric function of Ku70 was not observed in 14 other grasshopper and locust species, or in the mouse, thus suggesting that it is an autapomorphy in E. plorans.
Collapse
|
23
|
Gypsy, RTE and Mariner transposable elements populate Eyprepocnemis plorans genome. Genetica 2012; 140:365-74. [PMID: 23073915 DOI: 10.1007/s10709-012-9686-1] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2012] [Accepted: 10/08/2012] [Indexed: 12/12/2022]
Abstract
We analyze here the presence and abundance of three types of transposable elements (TEs), i.e. Gypsy, RTE and Mariner, in the genome of the grasshopper Eyprepocnemis plorans. PCR experiments allowed amplification, cloning and sequencing of these elements (EploGypI, EploRTE5, EploMar20) from the E. plorans genome. Fluorescent in situ hybridization (FISH) showed that all three elements are restricted to euchromatic regions, thus being absent from the pericentromeric region of all A chromosomes, which contain a satellite DNA (satDNA) and ribosomal DNA (rDNA), and being very scarce in B chromosomes mostly made up of these two types of repetitive DNA. FISH suggested that EploGypI is the most abundant and EploMar20 is the least abundant, with EploRTE5 showing intermediate abundance. An estimation of copy number, by means of quantitative PCR, showed that EploGypI is, by far, the most abundant element, followed by EploRTE5 and EploMar20, in consistency with FISH results. RNA isolation and PCR experiments on complementary DNA (cDNA) showed the presence of transcripts for the three TE elements. The implications of the preferential location of these TE elements into euchromatin, the significance of TE abundance in the giant genome of this species, and a possible relationship between TEs and B chromosome mutability, are discussed.
Collapse
|
24
|
Isolation of a cancer-associated microchromosome in the sperm-dependent parthenogen Poecilia formosa. Cytogenet Genome Res 2011; 135:135-42. [PMID: 21952475 DOI: 10.1159/000331271] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/11/2011] [Indexed: 11/19/2022] Open
Abstract
In the asexual all-female fish species Poecilia formosa, the Amazon molly, supernumerary chromosomes have frequently been found in both laboratory-reared and wild-caught individuals. While wild-caught individuals with B chromosomes are phenotypically indifferent from conspecifics, individuals carrying B chromosomes from recent introgression events in the laboratory show phenotypic changes. Former analyses showed that the expression of a pigment cell locus is associated with the presence of these B chromosomes. In addition, they contain a so far unidentified locus that confers a higher susceptibility to tumor formation in the presence of pigmentation pattern. Isolation by microdissection and hybridization to metaphase chromosomes revealed that they contain one or several sequences with similarity to a highly repetitive pericentromeric and subtelomeric sequence in A chromosomes. Isolation of one particular sequence by AFLP showed that the B chromosomes contain at least 1 copy of an A-chromosomal region which is highly conserved in the whole genus Poecilia, i.e. more than 5 million years old. We propose it to be a single copy sequence.
Collapse
|