Satellite DNA in insects: a review

High-throughput analysis of the satellitome revealed enormous diversity of satellite DNAs in the neo-Y chromosome of the cricket Eneoptera surinamensis

Satellite DNAs (satDNAs) constitute large portion of eukaryote genomes, comprising non-protein-coding sequences tandemly repeated. They are mostly found in heterochromatic regions of chromosomes such as around centromere or near telomeres, in intercalary heterochromatin, and often in non-recombining segments of sex chromosomes. We examined the satellitome in the cricket Eneoptera surinamensis (2n = 9, neo-X1X2Y, males) to characterize the molecular evolution of its neo-sex chromosomes. To achieve this, we analyzed illumina reads using graph-based clustering and complementary analyses. We found an unusually high number of 45 families of satDNAs, ranging from 4 bp to 517 bp, accounting for about 14% of the genome and showing different modular structures and high diversity of arrays. FISH mapping revealed that satDNAs are located mostly in C-positive pericentromeric regions of the chromosomes. SatDNAs enrichment was also observed in the neo-sex chromosomes in comparison to autosomes. Especially astonishing accumulation of satDNAs loci was found in the highly differentiated neo-Y, including 39 satDNAs over-represented in this chromosome, which is the greatest satDNAs diversity yet reported for sex chromosomes. Our results suggest possible involvement of satDNAs in genome increasing and in molecular differentiation of the neo-sex chromosomes in this species, contributing to the understanding of sex chromosome composition and evolution in Orthoptera.

W-enriched satellite sequence in the Indian meal moth, Plodia interpunctella (Lepidoptera, Pyralidae)

The W chromosome of most lepidopteran species represents the largest heterochromatin entity in the female genome. Although satellite DNA is a typical component of constitutive heterochromatin, there are only a few known satellite DNAs (satDNAs) located on the W chromosome in moths and butterflies. In this study, we isolated and characterized new satDNA (PiSAT1) from microdissected W chromosomes of the Indian meal moth, Plodia interpunctella. Even though the PiSAT1 is mainly localized near the female-specific segment of the W chromosome, short arrays of this satDNA also occur on autosomes and/or the Z chromosome. Probably due to the predominant location in the non-recombining part of the genome, PiSAT1 exhibits a relatively large nucleotide variability in its monomers. However, at least a part of all predicted functional motifs is located in conserved regions. Moreover, we detected polyadenylated transcripts of PiSAT1 in all developmental stages and in both sexes (female and male larvae, pupae and adults). Our results suggest a potential structural and functional role of PiSAT1 in the P. interpunctella genome, which is consistent with accumulating evidence for the important role of satDNAs in eukaryotic genomes.

Transcription of highly repetitive tandemly organized DNA in amphibians and birds: A historical overview and modern concepts

Tandemly organized highly repetitive DNA sequences are crucial structural and functional elements of eukaryotic genomes. Despite extensive evidence, satellite DNA remains an enigmatic part of the eukaryotic genome, with biological role and significance of tandem repeat transcripts remaining rather obscure. Data on tandem repeats transcription in amphibian and avian model organisms is fragmentary despite their genomes being thoroughly characterized. Review systematically covers historical and modern data on transcription of amphibian and avian satellite DNA in somatic cells and during meiosis when chromosomes acquire special lampbrush form. We highlight how transcription of tandemly repetitive DNA sequences is organized in interphase nucleus and on lampbrush chromosomes. We offer LTR-activation hypotheses of widespread satellite DNA transcription initiation during oogenesis. Recent explanations are provided for the significance of high-yield production of non-coding RNA derived from tandemly organized highly repetitive DNA. In many cases the data on the transcription of satellite DNA can be extrapolated from lampbrush chromosomes to interphase chromosomes. Lampbrush chromosomes with applied novel technical approaches such as superresolution imaging, chromosome microdissection followed by high-throughput sequencing, dynamic observation in life-like conditions provide amazing opportunities for investigation mechanisms of the satellite DNA transcription.

Karyotypes and Repetitive DNA Evolution in Six Species of the Genus Mahanarva (Auchenorrhyncha: Cercopidae)

Insects of the Cercopidae family are widely distributed and comprise 59 genera and 431 species in the New World. They are xylemophagous, causing losses in agricultural and pasture grasses, and are considered as emerging pests. Chromosomally, these insects have been studied by standard techniques, revealing variable diploid numbers and primarily X0 sex chromosome systems (males). We performed chromosome studies in 6 Mahanarva (Cercopidae) species using standard and differential chromosome staining as well as mapping of repetitive DNAs. Moreover, the relationship between the repetitive DNAs was analyzed at the interspecific level. A diploid chromosome number of 2n = 19,X0 was documented, with chromosomes gradually decreasing in size. Neutral or GC-rich regions were detected which varied depending on the species. Fluorescence in situ hybridization with a (TTAGG)n telomeric motif probe revealed terminal signals, matching those of the Cot DNAs obtained from each species, that were also restricted to the terminal regions of all chromosomes. Dot blot analysis with the Cot fraction from M. quadripunctata showed that at least part of the repetitive genome is shared among the 6 species. Our data highlight the conservation of chromosomal features and organization of repetitive DNAs in the genus Mahanarva, suggesting a low differentiation for chromosomes and repetitive DNAs in most of the 6 species studied.

Evolutionary dynamics of two satellite DNA families in rock lizards of the genus Iberolacerta (Squamata, Lacertidae): different histories but common traits

Satellite DNAs compose a large portion of all higher eukaryotic genomes. The turnover of these highly repetitive sequences is an important element in genome organization and evolution. However, information about the structure and dynamics of reptilian satellite DNA is still scarce. Two satellite DNA families, HindIII and TaqI, have been previously characterized in four species of the genus Iberolacerta. These families showed different chromosomal locations, abundances, and evolutionary rates. Here, we extend the study of both satellite DNAs (satDNAs) to the remaining Iberolacerta species, with the aim to investigate the patterns of variability and factors influencing the evolution of these repetitive sequences. Our results revealed disparate patterns but also common traits in the evolutionary histories of these satellite families: (i) each satellite DNA is made up of a library of monomer variants or subfamilies shared by related species; (ii) species-specific profiles of satellite repeats are shaped by expansions and/or contractions of different variants from the library; (iii) different turnover rates, even among closely related species, result in great differences in overall sequence homogeneity and in concerted or non-concerted evolution patterns, which may not reflect the phylogenetic relationships among taxa. Contrasting turnover rates are possibly related to genomic constraints such as karyotype architecture and the interspersed organization of diverging repeat variants in satellite arrays. Moreover, rapid changes in copy number, especially in the centromeric HindIII satDNA, may have been associated with chromosomal rearrangements and even contributed to speciation within Iberolacerta.

A refined Panax ginseng karyotype based on an ultra-high copy 167-bp tandem repeat and ribosomal DNAs

Background

Panax ginseng Meyer (Asian ginseng) has a large nuclear genome size of > 3.5 Gbp in haploid genome equivalent of 24 chromosomes. Tandem repeats (TRs) occupy significant portions of the genome in many plants and are often found in specific genomic loci, making them a valuable molecular cytogenetic tool in discriminating chromosomes. In an effort to understand the P. ginseng genome structure, we characterized an ultrahigh copy 167-bp TR (Pg167TR) and explored its chromosomal distribution as well as its utility for chromosome identification.

Methods

Polymerase chain reaction amplicons of Pg167TR were labeled, along with 5S and 45S rDNA amplicons, using a direct nick-translation method. Direct fluorescence in situ hybridization (FISH) was used to analyze the chromosomal distribution of Pg167TR.

Results

Recently, we reported a method of karyotyping the 24 chromosome pairs of P. ginseng using rDNA and DAPI (4′,6-diamidino-2-phenylindole) bands. Here, a unique distribution of Pg167TR in all 24 P. ginseng chromosomes was observed, allowing easy identification of individual homologous chromosomes. Additionally, direct labeling of 5S and 45S rDNA probes allowed the identification of two additional 5S rDNA loci not previously reported, enabling the refinement of the P. ginseng karyotype.

Conclusion

Identification of individual P. ginseng chromosomes was achieved using Pg167TR-FISH. Chromosome identification is important in understanding the P. ginseng genome structure, and our method will be useful for future integration of genetic linkage maps and genome scaffold anchoring. Additionally, it is a good tool for comparative studies with related species in efforts to understand the evolution of P. ginseng.

CpSAT-1, a transcribed satellite sequence from the codling moth, Cydia pomonella

Satellite DNA (satDNA) is a non-coding component of eukaryotic genomes, located mainly in heterochromatic regions. Relevance of satDNA began to emerge with accumulating evidence of its potential yet hardly comprehensible role that it can play in the genome of many organisms. We isolated the first satDNA of the codling moth (Cydia pomonella, Tortricidae, Lepidoptera), a species with holokinetic chromosomes and a single large heterochromatic element, the W chromosome in females. The satDNA, called CpSAT-1, is located on all chromosomes of the complement, although in different amounts. Surprisingly, the satellite is almost missing in the heterochromatic W chromosome. Additionally, we isolated mRNA from all developmental stages (1st-5th instar larva, pupa, adult), both sexes (adult male and female) and several tissues (Malpighian tubules, gut, heart, testes, and ovaries) of the codling moth and showed the CpSAT-1 sequence was transcribed in all tested samples. Using CpSAT-1 specific primers we amplified, cloned and sequenced 40 monomers from cDNA and gDNA, respectively. The sequence analysis revealed a high mutation rate and the presence of potentially functional motifs, mainly in non-conserved regions of the monomers. Both the chromosomal distribution and the sequence analysis suggest that CPSAT-1 has no function in the C. pomonella genome.

Is radon emission in caves causing deletions in satellite DNA sequences of cave-dwelling crickets?

The most stable isotope of radon, ²²²Rn, represents the major source of natural radioactivity in confined environments such as mines, caves and houses. In this study, we explored the possible radon-related effects on the genome of Dolichopoda cave crickets (Orthoptera, Rhaphidophoridae) sampled in caves with different concentrations of radon. We analyzed specimens from ten populations belonging to two genetically closely related species, D. geniculata and D. laetitiae, and explored the possible association between the radioactivity dose and the level of genetic polymorphism in a specific family of satellite DNA (pDo500 satDNA). Radon concentration in the analyzed caves ranged from 221 to 26000 Bq/m³. Specimens coming from caves with the highest radon concentration showed also the highest variability estimates in both species, and the increased sequence heterogeneity at pDo500 satDNA level can be explained as an effect of the mutation pressure induced by radon in cave. We discovered a specific category of nuclear DNA, the highly repetitive satellite DNA, where the effects of the exposure at high levels of radon-related ionizing radiation are detectable, suggesting that the satDNA sequences might be a valuable tool to disclose harmful effects also in other organisms exposed to high levels of radon concentration.

The draft genome of the pest tephritid fruit fly Bactrocera tryoni: resources for the genomic analysis of hybridising species

Background

The tephritid fruit flies include a number of economically important pests of horticulture, with a large accumulated body of research on their biology and control. Amongst the Tephritidae, the genus Bactrocera, containing over 400 species, presents various species groups of potential utility for genetic studies of speciation, behaviour or pest control. In Australia, there exists a triad of closely-related, sympatric Bactrocera species which do not mate in the wild but which, despite distinct morphologies and behaviours, can be force-mated in the laboratory to produce fertile hybrid offspring. To exploit the opportunities offered by genomics, such as the efficient identification of genetic loci central to pest behaviour and to the earliest stages of speciation, investigators require genomic resources for future investigations.

Results

We produced a draft de novo genome assembly of Australia’s major tephritid pest species, Bactrocera tryoni. The male genome (650 -700 Mbp) includes approximately 150Mb of interspersed repetitive DNA sequences and 60Mb of satellite DNA. Assessment using conserved core eukaryotic sequences indicated 98% completeness. Over 16,000 MAKER-derived gene models showed a large degree of overlap with other Dipteran reference genomes. The sequence of the ribosomal RNA transcribed unit was also determined. Unscaffolded assemblies of B. neohumeralis and B. jarvisi were then produced; comparison with B. tryoni showed that the species are more closely related than any Drosophila species pair. The similarity of the genomes was exploited to identify 4924 potentially diagnostic indels between the species, all of which occur in non-coding regions.

Conclusions

This first draft B. tryoni genome resembles other dipteran genomes in terms of size and putative coding sequences. For all three species included in this study, we have identified a comprehensive set of non-redundant repetitive sequences, including the ribosomal RNA unit, and have quantified the major satellite DNA families. These genetic resources will facilitate the further investigations of genetic mechanisms responsible for the behavioural and morphological differences between these three species and other tephritids. We have also shown how whole genome sequence data can be used to generate simple diagnostic tests between very closely-related species where only one of the species is scaffolded.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1153) contains supplementary material, which is available to authorized users.

A step to the gigantic genome of the desert locust: chromosome sizes and repeated DNAs

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.

Distribution and evolution of repeated sequences in genomes of Triatominae (Hemiptera-Reduviidae) inferred from genomic in situ hybridization

The subfamily Triatominae, vectors of Chagas disease, comprises 140 species characterized by a highly homogeneous chromosome number. We analyzed the chromosomal distribution and evolution of repeated sequences in Triatominae genomes by Genomic in situ Hybridization using Triatoma delpontei and Triatoma infestans genomic DNAs as probes. Hybridizations were performed on their own chromosomes and on nine species included in six genera from the two main tribes: Triatomini and Rhodniini. Genomic probes clearly generate two different hybridization patterns, dispersed or accumulated in specific regions or chromosomes. The three used probes generate the same hybridization pattern in each species. However, these patterns are species-specific. In closely related species, the probes strongly hybridized in the autosomal heterochromatic regions, resembling C-banding and DAPI patterns. However, in more distant species these co-localizations are not observed. The heterochromatic Y chromosome is constituted by highly repeated sequences, which is conserved among 10 species of Triatomini tribe suggesting be an ancestral character for this group. However, the Y chromosome in Rhodniini tribe is markedly different, supporting the early evolutionary dichotomy between both tribes. In some species, sex chromosomes and autosomes shared repeated sequences, suggesting meiotic chromatin exchanges among these heterologous chromosomes. Our GISH analyses enabled us to acquire not only reliable information about autosomal repeated sequences distribution but also an insight into sex chromosome evolution in Triatominae. Furthermore, the differentiation obtained by GISH might be a valuable marker to establish phylogenetic relationships and to test the controversial origin of the Triatominae subfamily.

Molecular and physical characterization of the complex pericentromeric heterochromatin of the vole species Microtus thomasi

A new repeated DNA from Microtus thomasi, Mth-Alu2.2, was cloned and characterized and is presented here for the first time. Digestion of genomic DNA from M. thomasi with AluI restriction enzyme revealed a 2.2-kb repetitive DNA sequence with a high AT content (69%). This sequence consists of a tandemly repeated nonanucleotide of the consensus sequence CACAATGTA, which constitutes approximately 93-95% of the total unit length. The location of the Mth-Alu2.2 sequence in the karyotype was determined by FISH, demonstrating strong hybridization signals in the pericentromeric regions of all chromosomes and in the heterochromatin blocks of several X chromosome variants. In addition, the distribution of the 4 pericentromeric repeat sequences Msat-160, Mth-Alu900, Mth-Alu2.2, and interstitial telomeric repeats was analyzed by in situ hybridization in M. thomasi, in order to shed light on the complex composition of the chromosomal pericentromeric regions in this species. The order and organization of these sequences in the pericentromeric regions are conserved, with slight variations in both the degree of overlapping and the amount of each repeated DNA in the chromosomes. Specifically, Mth-Alu2.2 is localized in the terminal regions of the chromosomes, with Msat-160 occupying the immediately inner region, partially intermixed with Mth-Alu2.2. The sequence Mth-Alu900 is found in internal positions below Msat-160, and the interstitial telomeric repeats are located close to the long-arm euchromatin of the chromosomes.

Evolutionary story of a satellite DNA from Phodopus sungorus (Rodentia, Cricetidae)

With the goal to contribute for the understanding of satellite DNA evolution and its genomic involvement, in this work it was isolated and characterized the first satellite DNA (PSUcentSat) from Phodopus sungorus (Cricetidae). Physical mapping of this sequence in P. sungorus showed large PSUcentSat arrays located at the heterochromatic (peri)centromeric region of five autosomal pairs and Y-chromosome. The presence of orthologous PSUcentSat sequences in the genomes of other Cricetidae and Muridae rodents was also verified, presenting however, an interspersed chromosomal distribution. This distribution pattern suggests a PSUcentSat-scattered location in an ancestor of Muridae/Cricetidae families, that assumed afterwards, in the descendant genome of P. sungorus a restricted localization to few chromosomes in the (peri)centromeric region. We believe that after the divergence of the studied species, PSUcentSat was most probably highly amplified in the (peri)centromeric region of some chromosome pairs of this hamster by recombinational mechanisms. The bouquet chromosome configuration (prophase I) possibly displays an important role in this selective amplification, providing physical proximity of centromeric regions between chromosomes with similar size and/or morphology. This seems particularly evident for the acrocentric chromosomes of P. sungorus (including the Y-chromosome), all presenting large PSUcentSat arrays at the (peri)centromeric region. The conservation of this sequence in the studied genomes and its (peri)centromeric amplification in P. sungorus strongly suggests functional significance, possibly displaying this satellite family different functions in the different genomes. The verification of PSUcentSat transcriptional activity in normal proliferative cells suggests that its transcription is not stage-limited, as described for some other satellites.

Effects of histone deacetylase inhibitor oxamflatin on in vitro porcine somatic cell nuclear transfer embryos

Low cloning efficiency is considered to be caused by the incomplete or aberrant epigenetic reprogramming of differentiated donor cells in somatic cell nuclear transfer (SCNT) embryos. Oxamflatin, a novel class of histone deacetylase inhibitor (HDACi), has been found to improve the in vitro and full-term developmental potential of SCNT embryos. In the present study, we studied the effects of oxamflatin treatment on in vitro porcine SCNT embryos. Our results indicated that the rate of in vitro blastocyst formation of SCNT embryos treated with 1 μM oxamflatin for 15 h postactivation was significantly higher than all other treatments. Treatment of oxamflatin decreased the relative histone deacetylase (HDAC) activity in cloned embryos and resulted in hyperacetylation levels of histone H3 at lysine 9 (AcH3K9) and histone H4 at lysine 5 (AcH4K5) at pronuclear, two-cell, and four-cell stages partly through downregulating HDAC1. The suppression of HDAC6 through oxamflatin increased the nonhistone acetylation level of α-tubulin during the mitotic cell cycle of early SCNT embryos. In addition, we demonstrated that oxamflatin downregulated DNA methyltransferase 1 (DNMT1) expression and global DNA methylation level (5-methylcytosine) in two-cell-stage porcine SCNT embryos. The pluripotency-related gene POU5F1 was found to be upregulated in the oxamflatin-treated group with a decreased DNA methylation tendency in its promoter regions. Treatment of oxamflatin did not change the locus-specific DNA methylation levels of Sus scrofa heterochromatic satellite DNA sequences at the blastocyst stage. Meanwhile, our findings suggest that treatment with HDACi may contribute to maintaining the stable status of cytoskeleton-associated elements, such as acetylated α-tubulin, which may be the crucial determinants of donor nuclear reprogramming in early SCNT embryos. In summary, oxamflatin treatment improves the developmental potential of porcine SCNT embryos in vitro.

Cytosine methylation of an ancient satellite family in the wild beet Beta procumbens

DNA methylation is an essential epigenetic feature for the regulation and maintenance of heterochromatin. Satellite DNA is a repetitive sequence component that often occurs in large arrays in heterochromatin of subtelomeric, intercalary and centromeric regions. Knowledge about the methylation status of satellite DNA is important for understanding the role of repetitive DNA in heterochromatization. In this study, we investigated the cytosine methylation of the ancient satellite family pEV in the wild beet Beta procumbens. The pEV satellite is widespread in species-specific pEV subfamilies in the genus Beta and most likely originated before the radiation of the Betoideae and Chenopodioideae. In B. procumbens, the pEV subfamily occurs abundantly and spans intercalary and centromeric regions. To uncover its cytosine methylation, we performed chromosome-wide immunostaining and bisulfite sequencing of pEV satellite repeats. We found that CG and CHG sites are highly methylated while CHH sites show only low levels of methylation. As a consequence of the low frequency of CG and CHG sites and the preferential occurrence of most cytosines in the CHH motif in pEV monomers, this satellite family displays only low levels of total cytosine methylation.

Centromere identity from the DNA point of view

The centromere is a chromosomal locus responsible for the faithful segregation of genetic material during cell division. It has become evident that centromeres can be established literally on any DNA sequence, and the possible synergy between DNA sequences and the most prominent centromere identifiers, protein components, and epigenetic marks remains uncertain. However, some evolutionary preferences seem to exist, and long-term established centromeres are frequently formed on long arrays of satellite DNAs and/or transposable elements. Recent progress in understanding functional centromere sequences is based largely on the high-resolution DNA mapping of sequences that interact with the centromere-specific histone H3 variant, the most reliable marker of active centromeres. In addition, sequence assembly and mapping of large repetitive centromeric regions, as well as comparative genome analyses offer insight into their complex organization and evolution. The rapidly advancing field of transcription in centromere regions highlights the functional importance of centromeric transcripts. Here, we comprehensively review the current state of knowledge on the composition and functionality of DNA sequences underlying active centromeres and discuss their contribution to the functioning of different centromere types in higher eukaryotes.

An ancient satellite DNA has maintained repetitive units of the original structure in most species of the living fossil plant genus Zamia

ZpS1 satellite DNA is specific to the genus Zamia and presents repetitive units organized as long arrays and also as very short arrays dispersed in the genome. We have characterized the structure of the ZpS1 repeats in 12 species representative of the whole geographic distribution of the genus. In most species, the clone most common sequences (cMCS) were so similar that a general most common sequence (GMCS) of the ZpS1 repetitive unit in the genus could be obtained. The few partial variations from the GMCS found in cMCS of some species correspond to variable positions present in most other species, as indicated by the clone consensus sequences (cCS). Two species have an additional species-specific variety of ZpS1 satellite. The dispersed repeats were found to contain more mutations than repeats from long arrays. Our results indicate that all or most species of Zamia inherited the ZpS1 satellite from a common ancestor in Miocene and have maintained repetitive units of the original structure till present. The features of ZpS1 satellite in the genus Zamia are poorly compatible with the model of concerted evolution, but they are perfectly consistent with a new model of satellite evolution based on experimental evidences indicating that a specific amplification-substitution repair mechanism maintains the homogeneity and stability of the repeats structure in each satellite DNA originally present in a species as long as the species exists.

The cytogenetic architecture of the aphid genome

In recent years aphids, with their well-defined polyphenism, have become favoured as model organisms for the study of epigenetic processes. The availability of the pea aphid (Acyrthosiphon pisum) genome sequence has engendered much research aimed at elucidating the mechanisms by which the phenotypic plasticity of aphids is inherited and controlled. Yet so far this research effort has paid little attention to the cytogenetic processes that play a vital part in the organisation, expression and inheritance of the aphid genome. Aphids have holocentric chromosomes, which have very different properties from the chromosomes with localised centromeres that are found in most other organisms. Here we review the diverse forms of aphid chromosome behaviour that occur during sex determination and male and female meiosis, often in response to environmental changes and mediated by endocrine factors. Remarkable differences occur, even between related species, that could have significant effects on the inheritance of all or parts of the genome. In relation to this, we review the particular features of the distribution of heterochromatin, rDNA genes and other repetitive DNA in aphid chromosomes, and discuss the part that these may play in the epigenetic modification of chromatin structure and function.

Origin and distribution of AT-rich repetitive DNA families in Triatoma infestans (Heteroptera)

Triatoma infestans, one of the most important vectors of Trypanosoma cruzi, is very interesting model, because it shows large interpopulation variation in the amount and distribution of heterochromatin. This polymorphism involved the three large pairs up to almost all autosomal pairs, including the sex chromosomes. To understand the dynamics of heterochromatin variation in T. infestans, we isolated the AT-rich satDNA portion of this insect using reassociation kinetics (C0t), followed by cloning, sequencing and FISH. After chromosome localization, immunolabeling with anti-5-methylcytosine, anti-H4K5ac and anti-H3K9me2 antibodies was performed to determine the functional characteristics of heterochromatin. The results allowed us to reorganize the karyotype of T. infestans in accordance with the distribution of the families of repetitive DNA using seven different markers. We found that two arrays with lengths of 79 and 33bp have a strong relationship with transposable element sequences, suggesting that these two families of satDNA probably originated from Polintons. The results also allowed us to identify at least four chromosome rearrangements involved in the amplification/dispersion of AT-rich satDNA of T. infestans. These data should be very useful in new studies including those examining the cytogenomic and population aspects of this very important species of insect.

Disparate molecular evolution of two types of repetitive DNAs in the genome of the grasshopper Eyprepocnemis plorans

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.

Molecular characterization and chromosomal distribution of a species-specific transcribed centromeric satellite repeat from the olive fruit fly, Bactrocera oleae

Satellite repetitive sequences that accumulate in the heterochromatin consist a large fraction of a genome and due to their properties are suggested to be implicated in centromere function. Current knowledge of heterochromatic regions of Bactrocera oleae genome, the major pest of the olive tree, is practically nonexistent. In our effort to explore the repetitive DNA portion of B. oleae genome, a novel satellite sequence designated BoR300 was isolated and cloned. The present study describes the genomic organization, abundance and chromosomal distribution of BoR300 which is organized in tandem, forming arrays of 298 bp-long monomers. Sequence analysis showed an AT content of 60.4%, a CENP-B like-motif and a high curvature value based on predictive models. Comparative analysis among randomly selected monomers demonstrated a high degree of sequence homogeneity (88%-97%) of BoR300 repeats, which are present at approximately 3,000 copies per haploid genome accounting for about 0.28% of the total genomic DNA, based on two independent qPCR approaches. In addition, expression of the repeat was also confirmed through RT-PCR, by which BoR300 transcripts were detected in both sexes. Fluorescence in situ hybridization (FISH) of BoR300 on mitotic metaphases and polytene chromosomes revealed signals to the centromeres of two out of the six chromosomes which indicated a chromosome-specific centromeric localization. Moreover, BoR300 is not conserved in the closely related Bactrocera species tested and it is also absent in other dipterans, but it's rather restricted to the B. oleae genome. This feature of species-specificity attributed to BoR300 satellite makes it a good candidate as an identification probe of the insect among its relatives at early development stages.

Medium-sized tandem repeats represent an abundant component of the Drosophila virilis genome

Background

Previously, we developed a simple method for carrying out a restriction enzyme analysis of eukaryotic DNA in silico, based on the known DNA sequences of the genomes. This method allows the user to calculate lengths of all DNA fragments that are formed after a whole genome is digested at the theoretical recognition sites of a given restriction enzyme. A comparison of the observed peaks in distribution diagrams with the results from DNA cleavage using several restriction enzymes performed in vitro have shown good correspondence between the theoretical and experimental data in several cases. Here, we applied this approach to the annotated genome of Drosophila virilis which is extremely rich in various repeats.

Results

Here we explored the combined approach to perform the restriction analysis of D. virilis DNA. This approach enabled to reveal three abundant medium-sized tandem repeats within the D. virilis genome. While the 225 bp repeats were revealed previously in intergenic non-transcribed spacers between ribosomal genes of D. virilis, two other families comprised of 154 bp and 172 bp repeats were not described. Tandem Repeats Finder search demonstrated that 154 bp and 172 bp units are organized in multiple clusters in the genome of D. virilis. Characteristically, only 154 bp repeats derived from Helitron transposon are transcribed.

Conclusion

Using in silico digestion in combination with conventional restriction analysis and sequencing of repeated DNA fragments enabled us to isolate and characterize three highly abundant families of medium-sized repeats present in the D. virilis genome. These repeats comprise a significant portion of the genome and may have important roles in genome function and structural integrity. Therefore, we demonstrated an approach which makes possible to investigate in detail the gross arrangement and expression of medium-sized repeats basing on sequencing data even in the case of incompletely assembled and/or annotated genomes.

Chromosomal diversification of diploid number, heterochromatin and rDNAs in two species of Phanaeus beetles (Scarabaeidae, Scarabaeinae)

The genus Phanaeus is included in the tribe Phanaeini, one of the most diverse tribes within the subfamily Scarabaeinae in terms of chromosomal characteristics. However, so far the species of this genus were not studied with differential cytogenetic techniques, limiting any inference of the probable mechanisms responsible for this diversity. In this work, several techniques were applied with the aim of cytogenetically characterizing two Phanaeus species. The karyotype found for Phanaeus (Notiophanaeus) chalcomelas was 2n = 12, neo-XY, and that of P. (N.) splendidulus was 2n = 20, Xy_p, considered primitive for the family Scarabaeidae. The chromosomes of both species showed a high amount of constitutive heterochromatin (CH), with blocks rich in base pairs GC (CMA₃⁺). Moreover, in P. (N.) chalcomelas the marks revealed by C-banding and fluorochrome staining were different in size, showing CH variability. Sites of 18S ribosomal DNA (rDNA) were identified in one autosomal pair of P. (N.) chalcomelas and in five autosomal pairs of P. (N.) splendidulus. On the other hand, only one autosomal pair exhibited 5S rDNA sequences in these species. The results suggest that the karyotype differentiation of the Phanaeus species studied here involved pericentric inversions and centric fusions, as well as mechanisms related to amplification and dispersion of CH and rDNA sequences.

Characterization of two unrelated satellite DNA families in the Colorado potato beetle Leptinotarsa decemlineata (Coleoptera, Chrysomelidae)

The Colorado potato beetle (Leptinotarsa decemlineata, family Chrysomelidae),a phytophagous insect, which feeds preferably on potatoes, constitutes a serious pest of this crop and causes extensive damage to tomatoes and egg plants. It has a remarkable ability to develop resistance quickly against insecticides and shows a diversified and flexible life history. Consequently, the control of this pest has become difficult, requiring the development of new alternative biotechnology-based strategies. Such strategies require a thorough knowledge of the beetle’s genome,including the repetitive DNA. Satellite DNA (stDNA), composed of long arrays of tandemly arranged repeat units, constitutes the major component of heterochromatin and is located mainly in centromeric and telomeric chromosomal regions. We have studied two different unrelated satellite-DNA families of which the consensus sequences were 295 and 109bp in length, named LEDE-I and LEDE-II, respectively.Both were AT-rich (70.8% and 71.6%, respectively). Predictive models of sequence-dependent DNA bending and the study of electrophoretic mobility on non-denaturing polyacrylamide gels have shown that the DNA was curved in both satellite-DNA families. Among other features, the chromosome localization of both stDNAs has been studied. In situ hybridization performed on meiotic and mitoticnuclei showed chromosomes, including the X chromosome, with zero, one, or two stDNAs. In recent years, it has been proposed that the repetitive DNA may play a key role in biological diversification processes. This is the first molecular and cytogenetic study conducted on L. decemlineata repetitive DNA and specifically on stDNA, which is one of the important constituents of eukaryotic genomes.

Comparative analysis of tandem repeats from hundreds of species reveals unique insights into centromere evolution

Background

Centromeres are essential for chromosome segregation, yet their DNA sequences evolve rapidly. In most animals and plants that have been studied, centromeres contain megabase-scale arrays of tandem repeats. Despite their importance, very little is known about the degree to which centromere tandem repeats share common properties between different species across different phyla. We used bioinformatic methods to identify high-copy tandem repeats from 282 species using publicly available genomic sequence and our own data.

Results

Our methods are compatible with all current sequencing technologies. Long Pacific Biosciences sequence reads allowed us to find tandem repeat monomers up to 1,419 bp. We assumed that the most abundant tandem repeat is the centromere DNA, which was true for most species whose centromeres have been previously characterized, suggesting this is a general property of genomes. High-copy centromere tandem repeats were found in almost all animal and plant genomes, but repeat monomers were highly variable in sequence composition and length. Furthermore, phylogenetic analysis of sequence homology showed little evidence of sequence conservation beyond approximately 50 million years of divergence. We find that despite an overall lack of sequence conservation, centromere tandem repeats from diverse species showed similar modes of evolution.

Conclusions

While centromere position in most eukaryotes is epigenetically determined, our results indicate that tandem repeats are highly prevalent at centromeres of both animal and plant genomes. This suggests a functional role for such repeats, perhaps in promoting concerted evolution of centromere DNA across chromosomes.

Large tandem repeats make up the chromosome bar code: a hypothesis

Much of tandem repeats' functional nature in any genome remains enigmatic because there are only few tools available for dissecting and elucidating the functions of repeated DNA. The large tandem repeat arrays (satellite DNA) found in two mouse whole-genome shotgun assemblies were classified into 4 superfamilies, 8 families, and 62 subfamilies. With the simplified variant of chromosome positioning of different tandem repeats, we noticed the nonuniform distribution instead of the positions reported for mouse major and minor satellites. It is visible that each chromosome possesses a kind of unique code made up of different large tandem repeats. The reference genomes allow marking only internal tandem repeats, and even with such a limited data, the colored "bar code" made up of tandem repeats is visible. We suppose that tandem repeats bare the mechanism for chromosomes to recognize the regions to be associated. The associations, initially established via RNA, become fixed by histone modifications (the histone or chromatin code) and specific proteins. In such a way, associations, being at the beginning flexible and regulated, that is, adjustable, appear as irreversible and inheritable in cell generations. Tandem repeat multiformity tunes the developed nuclei 3D pattern by sequential steps of associations. Tandem repeats-based chromosome bar code could be the carrier of the genome structural information; that is, the order of precise tandem repeat association is the DNA morphogenetic program. Tandem repeats are the cores of the distinct 3D structures postulated in "gene gating" hypothesis.

Analysis of two abundant, highly related satellites in the allotetraploid Nicotiana arentsii using double-strand conformation polymorphism analysis and sequencing

• Allopolyploidy, a driving force in plant evolution, can induce rapid structural changes in parental subgenomes. Here, we examined the fate of homologous subtelomeric satellites in intrasection allotetraploid Nicotiana arentsii formed from N. undulata and N. wigandioides progenitors < 200,000 yr ago. • We cloned and sequenced a number of monomers from progenitors and the allotetraploid. Structural features of both cloned and genomic monomers were studied using double-strand conformation polymorphism analysis. • Two homologous satellites were isolated from N. undulata (called NUNSSP) and N. wigandioides (NWISSP). While the NUNSSP monomers were highly homogeneous in nucleotide sequences, the NWISSP monomers formed two separate clades. Likewise, the genomic NUNSSP monomers showed less DNA conformation heterogeneity than NWISSP monomers, with distinct conformations. While both satellites predominantly occupy subtelomeric positions, a fraction of the NWISSP repeats was found in an intercalary location, supporting the hypothesis that dispersion prevents the repeats becoming homogeneous. Sequence, structural and chromosomal features of the parental satellites were faithfully inherited by N. arentsii. • Our study revealed that intergenomic homogenization of subtelomeric satellite repeats does not occur in N. arentsii allotetraploid. We propose that the sequence and structural divergence of subtelomeric satellites may render allopolyploid chromosomes less vulnerable to intergenomic exchanges.

Tandemly repeated DNA families in the mouse genome

Background

Functional and morphological studies of tandem DNA repeats, that combine high portion of most genomes, are mostly limited due to the incomplete characterization of these genome elements. We report here a genome wide analysis of the large tandem repeats (TR) found in the mouse genome assemblies.

Results

Using a bioinformatics approach, we identified large TR with array size more than 3 kb in two mouse whole genome shotgun (WGS) assemblies. Large TR were classified based on sequence similarity, chromosome position, monomer length, array variability, and GC content; we identified four superfamilies, eight families, and 62 subfamilies - including 60 not previously described. 1) The superfamily of centromeric minor satellite is only found in the unassembled part of the reference genome. 2) The pericentromeric major satellite is the most abundant superfamily and reveals high order repeat structure. 3) Transposable elements related superfamily contains two families. 4) The superfamily of heterogeneous tandem repeats includes four families. One family is found only in the WGS, while two families represent tandem repeats with either single or multi locus location. Despite multi locus location, TRPC-21A-MM is placed into a separated family due to its abundance, strictly pericentromeric location, and resemblance to big human satellites.

To confirm our data, we next performed in situ hybridization with three repeats from distinct families. TRPC-21A-MM probe hybridized to chromosomes 3 and 17, multi locus TR-22A-MM probe hybridized to ten chromosomes, and single locus TR-54B-MM probe hybridized with the long loops that emerge from chromosome ends. In addition to in silico predicted several extra-chromosomes were positive for TR by in situ analysis, potentially indicating inaccurate genome assembly of the heterochromatic genome regions.

Conclusions

Chromosome-specific TR had been predicted for mouse but no reliable cytogenetic probes were available before. We report new analysis that identified in silico and confirmed in situ 3/17 chromosome-specific probe TRPC-21-MM. Thus, the new classification had proven to be useful tool for continuation of genome study, while annotated TR can be the valuable source of cytogenetic probes for chromosome recognition.

Structure and population dynamics of the major satellite DNA in the red flour beetle Tribolium castaneum

In the beetle genus Tribolium, satellite DNAs comprise a significant amount of pericentromeric heterochromatin and are characterized by rapid turnover resulting in species specific profiles. In the present work we characterize the major pericentromeric satellite DNA TCAST of the beetle T. castaneum and analyse its population dynamics. Using direct sequencing of genomic PCR products we show that the TCAST satellite exists in the form of two related subfamilies: Tcast1a and Tcast1b that make up 20 and 15% of the genome, respectively. Tcast1a and Tcast1b have consensus sequences of 377 and 362 bp respectively, share an average similarity of 79% and are characterized by a divergent, subfamily specific region of approximately 100 bp. The two subfamilies are prevalently organized in the interspersed form, although a portion exists in the form of homogenous tandem arrays composed of only Tcast1a or Tcast1b. The pattern of restriction enzyme digestion indicates that Tcast1a and Tcast1b are organized in composite higher order repeats. Comparison of sequence variability of Tcast1a and Tcast1b among ten strains reveals a difference in the frequency of particular mutations present at some positions. However, no difference in the organization and in the amount of subfamilies was detected among strains. The results show that direct genomic sequencing can be a useful method for the detection of population specific features of satellite DNA. In the case of TCAST satellite DNA, changes in the mutational profiles seem to represent the first step in the genesis of a population specific satellite profile.

Identification of two new repetitive elements and chromosomal mapping of repetitive DNA sequences in the fish Gymnothorax unicolor (Anguilliformes: Muraenidae)

Muraenidae is a species-rich family, with relationships among genera and species and taxonomy that have not been completely clarified. Few cytogenetic studies have been conducted on this family, and all of them showed the same diploid chromosome number (2n=42) but with conspicuous karyotypic variation among species. The Mediterranean moray eel Gymnothorax unicolor was previously cytogenetically studied using classical techniques that allowed the characterization of its karyotype structure and the constitutive heterochromatin and argyrophilic nucleolar organizer regions (Ag-NORs) distribution pattern. In the present study, we describe two new repetitive elements (called GuMboI and GuDdeI) obtained from restricted genomic DNA of G. unicolor that were characterized by Southern blot and physically localized by in situ hybridization on metaphase chromosomes. As they are highly repetitive DNA sequences, they map in heterochromatic regions. However, while GuDdeI was localized in the centromeric regions, the GuMboI fraction was distributed on some centromeres and was co-localized with the nucleolus organizer region (NOR). Comparative analysis with other Mediterranean species such as Muraena helena pointed out that these DNA fractions are species-specific and could potentially be used for species discrimination. As a new contribution to the karyotype of this species, we found that the major ribosomal genes are localized on acrocentric chromosome 9 and that the telomeres of each chromosome are composed of a tandem repeat derived from a poly-TTAGGG DNA sequence, as it occurs in most vertebrate species. The results obtained add new information useful in comparative genomics at the chromosomal level and contribute to the cytogenetic knowledge regarding this fish family, which has not been extensively studied.

Fosmid library end sequencing reveals a rarely known genome structure of marine shrimp Penaeus monodon

Background

The black tiger shrimp (Penaeus monodon) is one of the most important aquaculture species in the world, representing the crustacean lineage which possesses the greatest species diversity among marine invertebrates. Yet, we barely know anything about their genomic structure. To understand the organization and evolution of the P. monodon genome, a fosmid library consisting of 288,000 colonies and was constructed, equivalent to 5.3-fold coverage of the 2.17 Gb genome. Approximately 11.1 Mb of fosmid end sequences (FESs) from 20,926 non-redundant reads representing 0.45% of the P. monodon genome were obtained for repetitive and protein-coding sequence analyses.

Results

We found that microsatellite sequences were highly abundant in the P. monodon genome, comprising 8.3% of the total length. The density and the average length of microsatellites were evidently higher in comparison to those of other taxa. AT-rich microsatellite motifs, especially poly (AT) and poly (AAT), were the most abundant. High abundance of microsatellite sequences were also found in the transcribed regions. Furthermore, via self-BlastN analysis we identified 103 novel repetitive element families which were categorized into four groups, i.e., 33 WSSV-like repeats, 14 retrotransposons, 5 gene-like repeats, and 51 unannotated repeats. Overall, various types of repeats comprise 51.18% of the P. monodon genome in length. Approximately 7.4% of the FESs contained protein-coding sequences, and the Inhibitor of Apoptosis Protein (IAP) gene and the Innexin 3 gene homologues appear to be present in high abundance in the P. monodon genome.

Conclusions

The redundancy of various repeat types in the P. monodon genome illustrates its highly repetitive nature. In particular, long and dense microsatellite sequences as well as abundant WSSV-like sequences highlight the uniqueness of genome organization of penaeid shrimp from those of other taxa. These results provide substantial improvement to our current knowledge not only for shrimp but also for marine crustaceans of large genome size.

Composition and epigenetic markers of heterochromatin in the aphid Aphis nerii (Hemiptera: Aphididae)

A detailed karyotype analysis of the oleander aphid Aphis nerii focusing on the distribution, molecular composition and epigenetic modifications of heterochromatin was done in order to better understand the structure and evolution of holocentric/holokinetic chromosomes in aphids. The female karyotype (2n = 8) consisted of 3 pairs of autosomes and a pair of X chromosomes that were the longest elements in the karyotype and carried a single, terminally located nucleolar organizer region. Males showed 2n = 7 chromosomes due to the presence of a single X chromosome. Heterochromatin was located in the X chromosomes only and consisted of 4 satellite DNAs that have been identified. A. nerii constitutive heterochromatin was enriched in mono-, di- and tri-methylated H3 histones and HP1 proteins but, interestingly, it lacked DNA methylation that was widespread in euchromatic chromosomal regions. These results suggest that aphid heterochromatin is assembled and condensed without any involvement of DNA methylation.

Characterization of the satellite DNA Msat-160 from species of Terricola (Microtus) and Arvicola (Rodentia, Arvicolinae)

In the subfamily Arvicolinae (Cricetidae, Rodentia) the satellite DNA Msat-160 has been so far described in only some species from the genus Microtus and in one species from another genus, Chionomys nivalis. Here we cloned and characterized this satellite in two new arvicoline species, Microtus (Terricola) savii and Arvicola amphibius (terrestris). We have also demonstrated, by PCR and FISH, its existence in the genomes of several other species from both genera. These results suggest that Msat-160 already occurred in the common ancestor of the four genera/subgenera of Arvicolinae (Microtus, Chionomys, Arvicola, and Terricola). In Arvicola and Terricola, Msat-160 showed the basic monomer length of 160 bp, although a higher-order repeat (HORs) of 640 bp could have been probably replacing the original monomeric unit in A. a. terrestris. Msat-160 was localized by FISH mostly on the pericentromeric regions of the chromosomes, but the signal intensity and the number of carrier chromosomes varied extremely even between closely related species, resulting in a species-specific pattern of chromosomal distribution of this satellite. Such a variable pattern most likely is a consequence of a rapid amplification and contraction of particular repeats in the pericentromeric regions of chromosomes. In addition, we proposed that the rapid variation of pericentromeric repeats is strictly related to the prolific species radiation and diversification of karyotypes that characterize Arvicolinae lineage. Finally, we performed phylogenetic analysis in this group of related species based on Msat-160 that results to be in agreement with previously reported phylogenies, derived from other molecular markers.

Low satellite DNA variability in natural populations of Drosophila antonietae involved in different evolutionary events

Drosophila antonietae is a cactophilic species that is found in the mesophilic forest of the Paraná-Paraguay river basin and in the dunes of the South Atlantic coast of Brazil. Although the genetic structure of the Paraná-Paraguay river basin populations has already been established, the relationship between these populations and those on the Atlantic coast is controversial. In this study, we compared 33 repetitive units of pBuM-2 satellite DNA isolated from individuals from 8 populations of D. antonietae in these geographic regions, including some populations found within a contact zone with the closely related D. serido. The pBuM-2 sequences showed low interpopulational variability. This result was interpreted as a consequence of both gene flow among the populations and unequal crossing over promoting homogenization of the tandem arrays. The results presented here, together with those of previous studies, highlight the use of pBuM-2 for solving taxonomic conflicts within the D. buzzatii species cluster.

Epigenetic effects of polymorphic Y chromosomes modulate chromatin components, immune response, and sexual conflict

Genetic conflicts between sexes and generations provide a foundation for understanding the functional evolution of sex chromosomes and sexually dimorphic phenotypes. Y chromosomes of Drosophila contain multi-megabase stretches of satellite DNA repeats and a handful of protein-coding genes that are monomorphic within species. Nevertheless, polymorphic variation in heterochromatic Y chromosomes of Drosophila result in genome-wide gene expression variation. Here we show that such naturally occurring Y-linked regulatory variation (YRV) can be detected in somatic tissues and contributes to the epigenetic balance of heterochromatin/euchromatin at three distinct loci showing position-effect variegation (PEV). Moreover, polymorphic Y chromosomes differentially affect the expression of thousands of genes in XXY female genotypes in which Y-linked protein-coding genes are not transcribed. The data show a disproportionate influence of YRV on the variable expression of genes whose protein products localize to the nucleus, have nucleic-acid binding activity, and are involved in transcription, chromosome organization, and chromatin assembly. These include key components such as HP1, Trithorax-like (GAGA factor), Su(var)3-9, Brahma, MCM2, ORC2, and inner centromere protein. Furthermore, mitochondria-related genes, immune response genes, and transposable elements are also disproportionally affected by Y chromosome polymorphism. These functional clusterings may arise as a consequence of the involvement of Y-linked heterochromatin in the origin and resolution of genetic conflicts between males and females. Taken together, our results indicate that Y chromosome heterochromatin serves as a major source of epigenetic variation in natural populations that interacts with chromatin components to modulate the expression of biologically relevant phenotypic variation.

Cytogenetic analysis of two Coprophanaeus species (Scarabaeidae) revealing wide constitutive heterochromatin variability and the largest number of 45S rDNA sites among Coleoptera

The Coleopterans of Scarabaeinae clade presents Coprophanaeus (Megaphanaeus) ensifer and C. (Coprophanaeus) cyanescens (Scarabaeidae) when they are studied cytogenetically by different techniques. The species present symmetric karyotypes, diploid number of 2n=20, and meta-submetacentric chromosomes. C. (M.) ensifer present an XY sex-determining mechanism and C. (C.) cyanescens an XY(p) parachute mechanism. Analysis of constitutive heterochromatin (CH) in the two species revealed the presence of diphasic autosomes, with log arm heterochromatics. Moreover, an additional heterochromatic block in four autosomal bivalents were observed in C. (M.) ensifer. CMA(3)/DA/DAPI fluorochrome staining detected CMA(3) positive heterochromatic blocks restricted to the sex chromosomes in C. (C.) cyanescens, whereas in C. (M.) ensifer CMA(3) positive pericentromeric blocks were present in all autosomes, in the Y chromosome and in the four additional heterochromatic blocks. DAPI staining was neutral in both species. Silver nitrate (AgNO(3)) staining was inefficient for the detection of the nucleolar organizer region (NORs), but showed affinity for the heterochromatic regions. Fluorescence in situ hybridization (FISH) revealed the presence of 45S rDNA sites in the terminal region of the three autosomal bivalents of C. (C.) cyanescens and in seven bivalents and the Y chromosome of C. (M.) ensifer. These results contribute to a better understanding of chromosome evolution in the genus Coprophanaeus, and demonstrate a wide CH variability and the largest number of ribosomal sites among Coleoptera.

An expansion of the genome size dataset for the insect order Hymenoptera, with a first test of parasitism and eusociality as possible constraints

Although the Hymenoptera represent a remarkably diverse and socioeconomically important group that is of considerable interest in genome biology, they remain understudied in terms of genome size. This study reports new genome size estimates for 89 species of ants, bees and wasps, representing 17 families and four superfamilies. These are used in a test of the hypothesis that genome sizes are constrained by traits associated with parasitism or eusociality. Not all parasitoid wasps exhibit small genomes, though a relationship based on specific types of parasitism may still occur; by contrast, there was no convincing evidence of a constraint relating to eusociality. The data provided here can be used to guide future research aimed at understanding the evolution of large-scale genomic properties in this order.

Molecular evolution of the pDo500 satellite DNA family in Dolichopoda cave crickets (Rhaphidophoridae)

BACKGROUND

Non-coding satellite DNA (satDNA) usually has a high turn-over rate frequently leading to species specific patterns. However, some satDNA families evolve more slowly and can be found in several related species. Here, we analyzed the mode of evolution of the pDo500 satDNA family of Dolichopoda cave crickets. In addition, we discuss the potential of slowly evolving satDNAs as phylogenetic markers.

RESULTS

We sequenced 199 genomic or PCR amplified satDNA repeats of the pDo500 family from 12 Dolichopoda species. For the 38 populations under study, 39 pDo500 consensus sequences were deduced. Phylogenetic analyses using Bayesian, Maximum Parsimony, and Maximum Likelihood approaches yielded largely congruent tree topologies. The vast majority of pDo500 sequences grouped according to species designation. Scatter plots and statistical tests revealed a significant correlation between genetic distances for satDNA and mitochondrial DNA. Sliding window analyses showed species specific patterns of variable and conserved regions. The evolutionary rate of the pDo500 satDNA was estimated to be 1.63-1.78% per lineage per million years.

CONCLUSIONS

The pDo500 satDNA evolves gradually at a rate that is only slightly faster than previously published rates of insect mitochondrial COI sequences. The pDo500 phylogeny was basically congruent with the previously published mtDNA phylogenies. Accordingly, the slowly evolving pDo500 satDNA family is indeed informative as a phylogenetic marker.

Cytogenetic characterization of Eurysternus caribaeus (Coleoptera: Scarabaeidae): evidence of sex-autosome fusion and diploid number reduction prior to species dispersion

Mitotic and meiotic chromosomes of several populations of Eurysternus caribaeus (Coleoptera: Scarabaeidae) were analysed through conventional staining, C-banding, base-specific fluorochromes, silver nitrate staining and fluorescent in situ hybridization (FISH). All specimens showed 2n = 8 in their karyotypes, with a neo-XY sex system (Y is a submetacentric and X a metacentric) and three pairs of submetacentric autosomes. The analysis of constitutive heterochromatin (CH) revealed small blocks located in the centromeric region of all chromosomes which do not present positive staining under the fluorochromes CMA3 and DAPI. Silver nitrate staining revealed that the nucleolar organizer region (NORs) is associated with the sex chromosomes. The FISH technique revealed that rDNA sites in the X and Y are different in size. Data from different populations indicate that the diploid number reduction (2n = 8) observed in E. caribaeus is established and presumably has preceded the dispersion of this species. Moreover, this reduction occasioned the translocation of rDNA sites to the sex chromosomes, X and Y, an uncommon pattern in Scarabaeidae that was observed for the first time by the FISH in this work.

Evolutionary dynamics and sites of illegitimate recombination revealed in the interspersion and sequence junctions of two nonhomologous satellite DNAs in cactophilic Drosophila species

Satellite DNA (satDNA) is a major component of genomes but relatively little is known about the fine-scale organization of unrelated satDNAs residing at the same chromosome location, and the sequence structure and dynamics of satDNA junctions. We studied the organization and sequence junctions of two nonhomologous satDNAs, pBuM and DBC-150, in three species from the neotropical Drosophila buzzatii cluster (repleta group). In situ hybridization to microchromosomes, interphase nuclei and extended DNA fibers showed frequent interspersion of the two satellites in D. gouveai, D. antonietae and, to a lesser extent, D. seriema. We isolated by PCR six pBuM x DBC-150 junctions: four are exclusive to D. gouveai and two are exclusive to D. antonietae. The six junction breakpoints occur at different positions within monomers, suggesting independent origin. Four junctions showed abrupt transitions between the two satellites, whereas two junctions showed a distinct 10 bp tandem duplication before the junction. Unlike pBuM, DBC-150 junction repeats are more variable than randomly cloned monomers and showed diagnostic features in common to a 3-monomer higher-order repeat seen in the sister species D. serido. The high levels of interspersion between pBuM and DBC-150 repeats suggest extensive rearrangements between the two satellites, maybe favored by specific features of the microchromosomes. Our interpretation is that the junctions evolved by multiples events of illegitimate recombination between nonhomologous satDNA repeats, with subsequent rounds of unequal crossing-over expanding the copy number of some of the junctions.