Diversity of PBI-DdeI satellite DNA in snakes correlates with rapid independent evolution and different functional roles.

To better understand PBI-DdeI satellite DNA located in the centromeric region of python, molecular evolution analysis was conducted on 40 snake species. A ladder-like pattern of DNA bands with repetition of the 194-210 bp monomer was observed in 15 species using PCR. Molecular cloning was performed to obtain 97 AT-rich monomer sequences. Phylogenetic and network analyses showed three PBI-DdeI subfamilies with sequences grouped in species-specific clusters, suggesting rapid evolution. Slow evolution was found in eight species with shared PBI-DdeI sequences, suggesting recent species diversification, allowing PBI-DdeI no time to diverge, with limited homogenization and fixation processes. Quantitative real-time PCR showed large differences in copy number between Python bivittatus and other snakes, consistent with repeat scanning of whole genome sequences. Copy numbers were significantly higher in female Naja kaouthia than in males, concurring with chromosomal distribution of PBI-DdeI specifically localized to female W chromosomes. PBI-DdeI might act as an evolutionary driver with several repeats to promote W chromosome differentiation and heterochromatinization in N. kaouthia. Analysis revealed PBI-DdeI with a reduced copy number, compared to P. bivittatus, in most snakes studied, and it is possible that it subsequently dispersed and amplified on W chromosomes with different functional roles in N. kaouthia.

Low coverage sequencing for repetitive DNA analysis in Passiflora edulis Sims: citogenomic characterization of transposable elements and satellite DNA.

BACKGROUND: The cytogenomic study of repetitive regions is fundamental for the understanding of morphofunctional mechanisms and genome evolution. Passiflora edulis a species of relevant agronomic value, this work had its genome sequenced by next generation sequencing and bioinformatics analysis performed by RepeatExplorer pipeline. The clusters allowed the identification and characterization of repetitive elements (predominant contributors to most plant genomes). The aim of this study was to identify, characterize and map the repetitive DNA of P. edulis, providing important cytogenomic markers, especially sequences associated with the centromere. RESULTS: Three clusters of satellite DNAs (69, 118 and 207) and seven clusters of Long Terminal Repeat (LTR) retrotransposons of the superfamilies Ty1/Copy and Ty3/Gypsy and families Angela, Athila, Chromovirus and Maximus-Sire (6, 11, 36, 43, 86, 94 and 135) were characterized and analyzed. The chromosome mapping of satellite DNAs showed two hybridization sites co-located in the 5S rDNA region (PeSat_1), subterminal hybridizations (PeSat_3) and hybridization in four sites, co-located in the 45S rDNA region (PeSat_2). Most of the retroelements hybridizations showed signals scattered in the chromosomes, diverging in abundance, and only the cluster 6 presented pericentromeric regions marking. No satellite DNAs and retroelement associated with centromere was observed. CONCLUSION: P. edulis has a highly repetitive genome, with the predominance of Ty3/Gypsy LTR retrotransposon. The satellite DNAs and LTR retrotransposon characterized are promising markers for investigation of the evolutionary patterns and genetic distinction of species and hybrids of Passiflora.

Distribution of DTHS3 satellite DNA across 12 bivalve species.

In this study, characterization of DTHS3 satellite DNA (satDNA) was further expanded within the class Bivalvia. Monomer variants of DTHS3 satDNA were compared in 12 bivalve species belonging to two different subclasses, Heterodonta and Pteriomorphia. This satDNA, whose age is estimated to a minimum of 516 Ma, is contributing to the concept of the dual character of satDNA sequences: their sequence preservation throughout long evolutionary periods and generation of species-specific variants of the same satDNA family.

Fangorn Forest (F2): a machine learning approach to classify genes and genera in the family Geminiviridae.

BACKGROUND: Geminiviruses infect a broad range of cultivated and non-cultivated plants, causing significant economic losses worldwide. The studies of the diversity of species, taxonomy, mechanisms of evolution, geographic distribution, and mechanisms of interaction of these pathogens with the host have greatly increased in recent years. Furthermore, the use of rolling circle amplification (RCA) and advanced metagenomics approaches have enabled the elucidation of viromes and the identification of many viral agents in a large number of plant species. As a result, determining the nomenclature and taxonomically classifying geminiviruses turned into complex tasks. In addition, the gene responsible for viral replication (particularly, the viruses belonging to the genus Mastrevirus) may be spliced due to the use of the transcriptional/splicing machinery in the host cells. However, the current tools have limitations concerning the identification of introns. RESULTS: This study proposes a new method, designated Fangorn Forest (F2), based on machine learning approaches to classify genera using an ab initio approach, i.e., using only the genomic sequence, as well as to predict and classify genes in the family Geminiviridae. In this investigation, nine genera of the family Geminiviridae and their related satellite DNAs were selected. We obtained two training sets, one for genus classification, containing attributes extracted from the complete genome of geminiviruses, while the other was made up to classify geminivirus genes, containing attributes extracted from ORFs taken from the complete genomes cited above. Three ML algorithms were applied on those datasets to build the predictive models: support vector machines, using the sequential minimal optimization training approach, random forest (RF), and multilayer perceptron. RF demonstrated a very high predictive power, achieving 0.966, 0.964, and 0.995 of precision, recall, and area under the curve (AUC), respectively, for genus classification. For gene classification, RF could reach 0.983, 0.983, and 0.998 of precision, recall, and AUC, respectively. CONCLUSIONS: Therefore, Fangorn Forest is proven to be an efficient method for classifying genera of the family Geminiviridae with high precision and effective gene prediction and classification. The method is freely accessible at www.geminivirus.org:8080/geminivirusdw/discoveryGeminivirus.jsp .

TAREAN: a computational tool for identification and characterization of satellite DNA from unassembled short reads.

Satellite DNA is one of the major classes of repetitive DNA, characterized by tandemly arranged repeat copies that form contiguous arrays up to megabases in length. This type of genomic organization makes satellite DNA difficult to assemble, which hampers characterization of satellite sequences by computational analysis of genomic contigs. Here, we present tandem repeat analyzer (TAREAN), a novel computational pipeline that circumvents this problem by detecting satellite repeats directly from unassembled short reads. The pipeline first employs graph-based sequence clustering to identify groups of reads that represent repetitive elements. Putative satellite repeats are subsequently detected by the presence of circular structures in their cluster graphs. Consensus sequences of repeat monomers are then reconstructed from the most frequent k-mers obtained by decomposing read sequences from corresponding clusters. The pipeline performance was successfully validated by analyzing low-pass genome sequencing data from five plant species where satellite DNA was previously experimentally characterized. Moreover, novel satellite repeats were predicted for the genome of Vicia faba and three of these repeats were verified by detecting their sequences on metaphase chromosomes using fluorescence in situ hybridization.

High genetic homogeneity points to a single introduction event responsible for invasion of Cotton leaf curl Multan virus and its associated betasatellite into China.

BACKGROUND: Cotton leaf curl Multan virus (CLCuMuV) is a Whitefly Transmitted Geminivirus (WTG) endemic to the India subcontinent and is notorious as a causal agent of cotton leaf curl disease (CLCuD), a major constraint to cotton production in south Asia. We found CLCuMuV infecting Hibiscus rosa-sinensis in Guangzhou, China in 2006. The spread and evolution of the invading CLCuMuV were monitored in the following nine years. FINDINGS: CLCuMuV spread rapidly in the last nine years and became established in Southern China. It infects at least five malvaceous plant species, H. rosa-sinensis, H. esculentus, Malvaiscus arboreus, Gossypium hirsutum and H. cannabinus. Complete nucleotide sequences of 34 geographically and/or temporally distinct CLCuMuV isolates were determined and analyzed together with six other publicly available genomes of CLCuMuV occurring in China. The 40 CLCuMuV isolates were found to share > 99 % nucleotide sequence identity with each other. In all cases tested, the CLCuMuVs were associated with a CLCuMuB. The 36 CLCuMuBs (30 sequenced by us) shared > 98 % nucleotide sequence identity. CONCLUSION: The high genetic homogeneity of CLCuMuV and CLCuMuB in China suggests the establishment of them from a single founder event.

Cotton leaf curl Burewala virus with intact or mutant transcriptional activator proteins: complexity of cotton leaf curl disease.

Cotton leaf curl disease (CLCuD) is a serious disease of cotton on the Indian subcontinent. In the present study, three cotton leaf curl viruses, cotton leaf curl Burewala virus (CLCuBuV), cotton leaf curl Kokhran virus (CLCuKoV) and cotton leaf curl Multan virus (CLCuMV), and their associated satellites, cotton leaf curl Multan betasatellite (CLCuMB) and cotton leaf curl Multan alphasatellite (CLCuMA), were detected. CLCuBuV with either intact (CLCuBuV-1) or mutant (CLCuBuV-2) transcriptional activator protein (TrAP) were detected in different plants. Agroinoculation with CLCuBuV-1 or CLCuBuV-2 together with CLCuMB and CLCuMA, resulted in typical leaf curling and stunting of tobacco plants. Inoculation with CLCuKoV or an isolate of CLCuMV (CLCuMV-2), together with CLCuMB and CLCuMA, induced severe leaf curling, while the other isolate of CLCuMV (CLCuMV-1), which was recombinant in origin, showed mild leaf curling in tobacco. To investigate the effect of intact or mutant TrAP and also the recombination events, CLCuBuV-1, CLCuBuV-2, CLCuMV-1 or CLCuMV-2 together with the satellites (CLCuMA and CLCuMB) were transferred to cotton via whitefly-mediated transmission. Cotton plants containing CLCuBuV-1, CLCuBuV-2 or CLCuMV-2 together with satellites showed curling and stunting, whereas the plants having CLCuMV-1 and the satellites showed only mild and indistinguishable symptoms. CLCuBuV-1 (intact TrAP) showed severe symptoms in comparison to CLCuBuV-2 (mutant TrAP). The present study reveals that two types of CLCuBuV, one with an intact TrAP and the other with a mutant TrAP, exist in natural infection of cotton in India. Additionally, CLCuMuV-1, which has a recombinant origin, induces mild symptoms in comparison to the other CLCuMV isolates.

Circomics of Cuban geminiviruses reveals the first alpha-satellite DNA in the Caribbean.

Circomics (circular DNA genomics), the combination of rolling circle amplification (RCA), restriction fragment length polymorphism (RFLP) analysis and pyro-sequencing, has been used recently to identify geminiviruses with high efficiency and low costs. Circular DNAs associated with Cuban geminiviruses were characterised by RCA/RFLP analysis and 454 sequencing of two batches of DNA amplified from selected plant samples as well as individual cloning and Sanger sequencing of DNA components and compared to other geminiviral DNAs by phylogenetic analysis. Cuban geminiviruses that were closely related to each other challenged the circomics approach. Ten geminiviral components and one alpha-satellite DNA were determined and compared to three geminiviral components obtained by conventional cloning. New strains of Sida yellow mottle virus (SiYMoV), tomato yellow distortion leaf virus (ToYDLV), Sida golden mosaic Florida virus (SiGMFV) and Sida golden mosaic Liguanea virus (SiGMLV) are described with host plant species being classified by molecular PCR-based bar coding. A new virus species is named Peristrophe mosaic virus. The first alpha-satellite found in Middle America establishes the New World branch of these elements which are related to nanoviruses and were previously thought to be restricted to the Old World. In conclusion, circomics is efficient for complex infections and closely related viruses to detected unexpected viral DNAs, but may need some scrutinisation by direct sequencing and cloning of individual components for certain cases.

Viral metagenomics: analysis of begomoviruses by illumina high-throughput sequencing.

Traditional DNA sequencing methods are inefficient, lack the ability to discern the least abundant viral sequences, and ineffective for determining the extent of variability in viral populations. Here, populations of single-stranded DNA plant begomoviral genomes and their associated beta- and alpha-satellite molecules (virus-satellite complexes) (genus, Begomovirus; family, Geminiviridae) were enriched from total nucleic acids isolated from symptomatic, field-infected plants, using rolling circle amplification (RCA). Enriched virus-satellite complexes were subjected to Illumina-Next Generation Sequencing (NGS). CASAVA and SeqMan NGen programs were implemented, respectively, for quality control and for de novo and reference-guided contig assembly of viral-satellite sequences. The authenticity of the begomoviral sequences, and the reproducibility of the Illumina-NGS approach for begomoviral deep sequencing projects, were validated by comparing NGS results with those obtained using traditional molecular cloning and Sanger sequencing of viral components and satellite DNAs, also enriched by RCA or amplified by polymerase chain reaction. As the use of NGS approaches, together with advances in software development, make possible deep sequence coverage at a lower cost; the approach described herein will streamline the exploration of begomovirus diversity and population structure from naturally infected plants, irrespective of viral abundance. This is the first report of the implementation of Illumina-NGS to explore the diversity and identify begomoviral-satellite SNPs directly from plants naturally-infected with begomoviruses under field conditions.

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.

A survey of geminiviruses and associated satellite DNAs in the cotton-growing areas of northwestern India.

Severe symptoms of cotton leaf curl disease (CLCuD) are caused by the association of a single-stranded circular DNA satellite (betasatellite) with a helper begomovirus. In this study, we analyzed 40 leaf samples (primarily cotton with CLCuD symptoms and other plants growing close by) from four sites between New Delhi and the Pakistan/India border, using rolling-circle amplification (RCA) and PCR. In total, the complete sequences of 12 different helper viruses, eight alphasatellites, and one betasatellite from five different plant species were obtained. A recombinant helper virus molecule found in okra and a novel alphasatellite-related DNA from croton are also described. This is the first report of the presence of both DNA components (helper virus and betasatellite) associated with resistance-breaking CLCuD in India, and it highlights the need for further work to combat its damage and spread.

A satellite DNA evolutionary analysis in the North American endemic dioecious plant Rumex hastatulus (Polygonaceae).

We studied the evolution of RAE180 satellite DNA family in the North American endemic dioecious plant Rumex hastatulus. In this species, the Texas race is characterized by a single XX/XY sex chromosome system, whereas the North Carolina race has evolved a derived complex XX/XY(1)Y(2) sex chromosome system. RAE180 repeats were autosomic and poorly represented (2 × 10(-4)% of the genome) with no differences between individuals of different genders or different races of R. hastatulus. In fact, the sex chromosomes of the North Carolina race are still euchromatic, and they have not accumulated satellite DNA sequences, which contrasts with that occurring in the rest of dioecious XX/XY(1)Y(2) Rumex species. In R. hastatulus, we detected the existence of three RAE180 subfamilies. Notwithstanding, while in the Texas race the TX1/NC1 subfamily is the most frequent, the TX2/NC2 subfamily is the most abundant in the North Carolina race. Additionally, the third, less represented subfamily (TX3/NC3) appears currently as relict sequences in both genomes. A common feature of RAE180 satellite is the sudden replacement of one sequence variant by another in different species (or populations as in R. hastatulus races). Thus, the phylogenetic analysis of RAE180 repeats from six dioecious Rumex species supports the "library" hypothesis. According to this hypothesis, we assume that a set of divergent RAE180 variants were present in the ancestral genome of dioecious Rumex species, from which novel tandem arrays originated by the amplification of different variants in different lineages. Differential levels of RAE180 satellite DNA amplification in each lineage, at different evolutionary times, and in different chromosomal positions gave rise to differential patterns of sequence evolution.

Megasatellites: a new class of large tandem repeats discovered in the pathogenic yeast Candida glabrata.

Megasatellites are DNA tandem arrays made of large motifs; they were discovered in the yeast Candida glabrata. They are widespread in this species (40 copies) but are not found in any other hemiascomycete so far, raising the intriguing question of their origin. They are found mainly in genes encoding cell wall products, suggesting that megasatellites were selected for a function linked to cell-cell adhesion or to pathogenicity. Their putative role in promoting genome rearrangements by interfering with DNA replication will also be discussed.

Low rates of homogenization of the DBC-150 satellite DNA family restricted to a single pair of microchromosomes in species from the Drosophila buzzatii cluster.

A satellite DNA family, termed DBC-150, comprises slightly GC-rich repeat units of approximately 150 bp that were isolated (by DNA digestions or PCR) from the genome of all seven Drosophila species from the buzzatii cluster (repleta group). The presence of subrepeats suggests that part of the extant DBC-150 monomer originated by the duplication of small sequence motifs. The DBC-150 family is compared to the previously described pBuM satDNA family, an abundant component of the genome of five species of the cluster. The two families are different in several aspects, including primary structure, A + T content, intraspecific and interspecific variability and rates of homogenization (or nucleotide spread). The data indicate a lower rate of homogenization (and absence of complete concerted evolution) of the DBC-150 compared to the pBuM family. FISH on metaphase chromosomes revealed that the DBC-150 family is located exclusively in the microchromosomes. To our knowledge this is the first record of a complex Drosophila satDNA restricted to a single pair of microchromosomes. The observed low rates of homogenization of the DBC-150 family might be related to a presumed reduction or suppression of meiotic recombination in the microchromosomes.

The Evolution of satellite III DNA subfamilies among primates.

We demonstrate that satellite III (SatIII) DNA subfamilies cloned from human acrocentric chromosomes arose in the Hominoidea superfamily. Two groups, distinguished by sequence composition, evolved nonconcurrently, with group 2 evolving 16-23 million years ago (MYA) and the more recent group 1 sequences emerging approximately 4.5 MYA. We also show the relative order of emergence of each group 2 subfamily in the various primate species. Our results show that each SatIII subfamily is an independent evolutionary unit, that the rate of evolution is not uniform between species, and that the evolution within a species is not uniform between chromosomes.

Geminivirus disease complexes: the threat is spreading.

Symptom-modulating DNA satellites associated with geminiviruses have come to our attention only recently but have proven to be widespread, associated with many diseases throughout the Old World, and economically significant, particularly in developing countries. Recent developments are elucidating the role played by these novel molecules in pathogenicity and in overcoming host plant defense. Further investigation into the promiscuous nature of these satellites and their ability to recruit further begomoviruses indicates that regions not yet affected by such begomovirus-satellite complexes are at great risk.

Sequence divergence and conservation in organizationally distinct subfamilies of Donax trunculus satellite DNA.

Characterization of a low-copy number DTF1 satellite DNA detected in the bivalve mollusk Donax trunculus revealed extensive grouping of monomer sequence variants into subfamilies identified by distinctive combinations of diagnostic nucleotides. It can be anticipated that a large number of subfamilies exists in the genome. In addition to the tandem organization of 169 bp long monomers, at least one subfamily was created through amplification of adjacent repeats in a higher order register. This complex satellite unit consists of two distinctive monomer variants that differ both in specific nucleotide changes and in a deleted segment partially substituted with a short unrelated sequence element. Most of the nucleotide substitutions differing between subfamilies are highly homogenized within a corresponding group of monomer variants, and intra-subfamily variability in general is low. Nucleotide diversity analysis of all sequenced variants of DTF1 satellite revealed the presence of two conserved segments, while the rest of the monomer sequence shows uniform and considerably higher level of variability. The persistence of conserved segments stands in contrast to the sequence and organizational divergence of monomer variant groups, and may indicate constraints in the evolution of DTF1 satellite repeats.

Comparative analysis of genomic sequences suggests that Trypanosoma cruzi CL Brener contains two sets of non-intercalated repeats of satellite DNA that correspond to T. cruzi I and T. cruzi II types.

Approximately 10% of the Trypanosoma cruzi genome is formed by a satellite DNA, composed by 195-bp repeats organized in 30+/-10 kb clusters in some, but not all chromosomes. Here, the satellite DNA of six representative T. cruzi strains was sequenced and used for phylogenetic inference. The results show that CL Brener contains satellite repeats from T. cruzi I and T. cruzi II strains, although type II sequences are more abundant. The presence of types I and II sequences extends previous propositions that genetic exchange between the two major T. cruzi lineages have occurred in CL Brener, although our data accommodate alternative scenarios of hybridization within T. cruzi II, as proposed by others. Altogether, present data suggest a complex origin for CL Brener. Sequence analysis of satellites isolated from chromosomal bands indicates that satellite DNA sequences are not chromosome specific. Neighbor analysis of in tandem satellite DNAs containing up to five repeats shows that each cluster contains only one type of sequence. Consequently, clusters with intercalated types I and II repeats were not found. We propose that the CL Brener genome contains large pieces of satellite DNA originated mainly from chromosomes of T. cruzi II with introgression of T. cruzi I lineage.

Organization of the KpnI family of chromosomal distal-end satellite DNAs in Silene latifolia.

The major satellite DNAs of the dioecious plant Silene latifolia are represented by the repetitive sequences X43.1, RMY1 and members of the SacI family, which are located at the distal ends of chromosomes. To characterize the satellite DNAs at the distal ends of the chromosomes in S. latifolia ( Sl-distal-satDNA), we isolated a bacterial artificial chromosome clone (number 15B12) that contained multiple repeat sequences with KpnI restriction sites, and subcloned a portion of this sequence into a plasmid vector. Sequencing analysis confirmed that recognition or degenerate sites for KpnI were repeated 26 times at intervals of 310-324 bp in the inserted DNA. The phylogenetic tree that was constructed with the 26 KpnI repeat units contained clustered branches that were independent of the SacI family. It is clear that the KpnI repeat belongs to an Sl-distal-satDNA family that is distinct from the SacI family. We designated this family as " KpnI" after the restriction enzyme that does not have a site in the SacI family. Multi-colored fluorescent in situ hybridization was performed with the KpnI family and RMY1 probes under high stringency conditions. The results suggest that chromosome 7 is unique and that it carries the KpnI family at only one end.

Extrachromosomal circular DNA of tandemly repeated genomic sequences in Drosophila.

One characteristic of genomic plasticity is the presence of extrachromosomal circular DNA (eccDNA). This DNA is found in various eukaryotes from yeast to humans, and its levels are elevated by exposure to carcinogens. eccDNA is heterogeneous in size and composed of chromosomal sequences. In this study we used two-dimensional gel electrophoresis to detect and characterize eccDNA in Drosophila. We found eccDNA throughout the fly's life cycle. These molecules comprise up to 10% of the total repetitive DNA content, and their size ranges from <1 kb to >20 kb. The eccDNA population contains circular multimers of tandemly repeated genes such as histones, rDNA, Stellate, and the Suppressor of Stellate. Multimers of centromeric heterochromatin sequences are included in eccDNA as well. Our findings are consistent with the hypothesis that intramolecular homologous recombination between direct tandem repeats is a favorite mechanism for eccDNA formation. The level of eccDNA increased following MMS treatment of wild-type larvae, consistent with phenomena observed in cultured mammalian cells. This shows mutagen-induced eccDNA formation in the context of the whole organism for the first time. Mutations in the genes okra, mus309, and mei41 did not affect eccDNA under normal conditions or following mutagen treatment, implying that eccDNA formation is different from known pathways of DNA repair.