Telomeric retrotransposons show propensity to form G-quadruplexes in various eukaryotic species

BACKGROUND

Canonical telomeres (telomerase-synthetised) are readily forming G-quadruplexes (G4) on the G-rich strand. However, there are examples of non-canonical telomeres among eukaryotes where telomeric tandem repeats are invaded by specific retrotransposons. Drosophila melanogaster represents an extreme example with telomeres composed solely by three retrotransposons-Het-A, TAHRE and TART (HTT). Even though non-canonical telomeres often show strand biased G-distribution, the evidence for the G4-forming potential is limited.

RESULTS

Using circular dichroism spectroscopy and UV absorption melting assay we have verified in vitro G4-formation in the HTT elements of D. melanogaster. Namely 3 in Het-A, 8 in TART and 2 in TAHRE. All the G4s are asymmetrically distributed as in canonical telomeres. Bioinformatic analysis showed that asymmetric distribution of potential quadruplex sequences (PQS) is common in telomeric retrotransposons in other Drosophila species. Most of the PQS are located in the gag gene where PQS density correlates with higher DNA sequence conservation and codon selection favoring G4-forming potential. The importance of G4s in non-canonical telomeres is further supported by analysis of telomere-associated retrotransposons from various eukaryotic species including green algae, Diplomonadida, fungi, insects and vertebrates. Virtually all analyzed telomere-associated retrotransposons contained PQS, frequently with asymmetric strand distribution. Comparison with non-telomeric elements showed independent selection of PQS-rich elements from four distinct LINE clades.

CONCLUSION

Our findings of strand-biased G4-forming motifs in telomere-associated retrotransposons from various eukaryotic species support the G4-formation as one of the prerequisites for the recruitment of specific retrotransposons to chromosome ends and call for further experimental studies.

Variation in G-quadruplex sequence and topology differentially impacts human DNA polymerase fidelity

G-quadruplexes (G4s), a type of non-B DNA, play important roles in a wide range of molecular processes, including replication, transcription, and translation. Genome integrity relies on efficient and accurate DNA synthesis, and is compromised by various stressors, to which non-B DNA structures such as G4s can be particularly vulnerable. However, the impact of G4 structures on DNA polymerase fidelity is largely unknown. Using an in vitro forward mutation assay, we investigated the fidelity of human DNA polymerases delta (δ4, four-subunit), eta (η), and kappa (κ) during synthesis of G4 motifs representing those in the human genome. The motifs differ in sequence, topology, and stability, features that may affect DNA polymerase errors. Polymerase error rate hierarchy (δ4 < κ < η) is largely maintained during G4 synthesis. Importantly, we observed unique polymerase error signatures during synthesis of VEGF G4 motifs, stable G4s which form parallel topologies. These statistically significant errors occurred within, immediately flanking, and encompassing the G4 motif. For pol δ4, the errors were deletions, insertions and complex errors within the G4 or encompassing the G4 motif and surrounding sequence. For pol η, the errors occurred in 3' sequences flanking the G4 motif. For pol κ, the errors were frameshift mutations within G-tracts of the G4. Because these error signatures were not observed during synthesis of an antiparallel G4 and, to a lesser extent, a hybrid G4, we suggest that G4 topology and/or stability could influence polymerase fidelity. Using in silico analyses, we show that most polymerase errors are predicted to have minimal effects on predicted G4 stability. Our results provide a unique view of G4s not previously elucidated, showing that G4 motif heterogeneity differentially influences polymerase fidelity within the motif and flanking sequences. Thus, our study advances the understanding of how DNA polymerase errors contribute to G4 mutagenesis.

Overexpression of a flower-specific aerolysin-like protein from the dioecious plant Rumex acetosa alters flower development and induces male sterility in transgenic tobacco

Sex determination in Rumex acetosa, a dioecious plant with a complex XY1 Y2 sex chromosome system (females are XX and males are XY1 Y2 ), is not controlled by an active Y chromosome but depends on the ratio between the number of X chromosomes and autosomes. To gain insight into the molecular mechanisms of sex determination, we generated a subtracted cDNA library enriched in genes specifically or predominantly expressed in female floral buds in early stages of development, when sex determination mechanisms come into play. In the present paper, we report the molecular and functional characterization of FEM32, a gene encoding a protein that shares a common architecture with proteins in different plants, animals, bacteria and fungi of the aerolysin superfamily; many of these function as β pore-forming toxins. The expression analysis, assessed by northern blot, RT-PCR and in situ hybridization, demonstrates that this gene is specifically expressed in flowers in both early and late stages of development, although its transcripts accumulate much more in female flowers than in male flowers. The ectopic expression of FEM32 under both the constitutive promoter 35S and the flower-specific promoter AP3 in transgenic tobacco showed no obvious alteration in vegetative development but was able to alter floral organ growth and pollen fertility. The 35S::FEM32 and AP3::FEM32 transgenic lines showed a reduction in stamen development and pollen viability, as well as a diminution in fruit set, fruit development and seed production. Compared with other floral organs, pistil development was, however, enhanced in plants overexpressing FEM32. According to these effects, it is likely that FEM32 functions in Rumex by arresting stamen and pollen development during female flower development. The aerolysin-like pore-forming proteins of eukaryotes are mainly involved in defence mechanisms against bacteria, fungi and insects and are also involved in apoptosis and programmed cell death (PCD), a mechanism that could explain the role of FEM32 in Rumex sex determination.

From Mendel's discovery on pea to today's plant genetics and breeding : Commemorating the 150th anniversary of the reading of Mendel's discovery

KEY MESSAGE

This work discusses several selected topics of plant genetics and breeding in relation to the 150th anniversary of the seminal work of Gregor Johann Mendel. In 2015, we celebrated the 150th anniversary of the presentation of the seminal work of Gregor Johann Mendel. While Darwin's theory of evolution was based on differential survival and differential reproductive success, Mendel's theory of heredity relies on equality and stability throughout all stages of the life cycle. Darwin's concepts were continuous variation and "soft" heredity; Mendel espoused discontinuous variation and "hard" heredity. Thus, the combination of Mendelian genetics with Darwin's theory of natural selection was the process that resulted in the modern synthesis of evolutionary biology. Although biology, genetics, and genomics have been revolutionized in recent years, modern genetics will forever rely on simple principles founded on pea breeding using seven single gene characters. Purposeful use of mutants to study gene function is one of the essential tools of modern genetics. Today, over 100 plant species genomes have been sequenced. Mapping populations and their use in segregation of molecular markers and marker-trait association to map and isolate genes, were developed on the basis of Mendel's work. Genome-wide or genomic selection is a recent approach for the development of improved breeding lines. The analysis of complex traits has been enhanced by high-throughput phenotyping and developments in statistical and modeling methods for the analysis of phenotypic data. Introgression of novel alleles from landraces and wild relatives widens genetic diversity and improves traits; transgenic methodologies allow for the introduction of novel genes from diverse sources, and gene editing approaches offer possibilities to manipulate gene in a precise manner.

Quadruplex-forming sequences occupy discrete regions inside plant LTR retrotransposons

Retrotransposons with long terminal repeats (LTR) form a significant proportion of eukaryotic genomes, especially in plants. They have gag and pol genes and several regulatory regions necessary for transcription and reverse transcription. We searched for potential quadruplex-forming sequences (PQSs) and potential triplex-forming sequences (PTSs) in 18 377 full-length LTR retrotransposons collected from 21 plant species. We found that PQSs were often located in LTRs, both upstream and downstream of promoters from which the whole retrotransposon is transcribed. Upstream-located guanine PQSs were dominant in the minus DNA strand, whereas downstream-located guanine PQSs prevailed in the plus strand, indicating their role both at transcriptional and post-transcriptional levels. Our circular dichroism spectroscopy measurements confirmed that these PQSs readily adopted guanine quadruplex structures-some of them were paralell-stranded, while others were anti-parallel-stranded. The PQS often formed doublets at a mutual distance of up to 400 bp. PTSs were most abundant in 3'UTR (but were also present in 5'UTR). We discuss the potential role of quadruplexes and triplexes as the regulators of various processes participating in LTR retrotransposon life cycle and as potential recombination sites during post-insertional retrotransposon-based genome rearrangements.

Expansion of microsatellites on evolutionary young Y chromosome

Sex chromosomes are an ideal system to study processes connected with suppressed recombination. We found evidence of microsatellite expansion, on the relatively young Y chromosome of the dioecious plant sorrel (Rumex acetosa, XY1Y2 system), but no such expansion on the more ancient Y chromosomes of liverwort (Marchantia polymorpha) and human. The most expanding motifs were AC and AAC, which also showed periodicity of array length, indicating the importance of beginnings and ends of arrays. Our data indicate that abundance of microsatellites in genomes depends on the inherent expansion potential of specific motifs, which could be related to their stability and ability to adopt unusual DNA conformations. We also found that the abundance of microsatellites is higher in the neighborhood of transposable elements (TEs) suggesting that microsatellites are probably targets for TE insertions. This evidence suggests that microsatellite expansion is an early event shaping the Y chromosome where this process is not opposed by recombination, while accumulation of TEs and chromosome shrinkage predominate later.

The key role of repeated DNAs in sex chromosome evolution in two fish species with ZW sex chromosome system

Despite substantial progress, there are still several gaps in our knowledge about the process of sex chromosome differentiation. The degeneration of sex-specific chromosome in some species is well documented, but it is not clear if all species follow the same evolutionary pathway. The accumulation of repetitive DNA sequences, however, is a common feature. To better understand this involvement, fish species emerge as excellent models because they exhibit a wide variety of sex chromosome and sex determining systems. Besides, they have much younger sex chromosomes compared to higher vertebrates, making it possible to follow early steps of differentiation. Here, we analyzed the arrangement of 9 repetitive DNA sequences in the W chromosomes of 2 fish species, namely Leporinus reinhardti and Triportheus auritus, which present well-differentiated ZZ/ZW sex system, but differ in respect to the size of the sex-specific chromosome. Both W chromosomes are almost fully heterochromatic, with accumulation of repeated DNAs in their heterochromatic regions. We found that microsatellites have strongly accumulated on the large W chromosome of L. reinhardti but not on the reduced-size W chromosome of T. auritus and are therefore important players of the W chromosome expansion. The present data highlight that the evolution of the sex chromosomes can diverge even in the same type of sex system, with and without the degeneration of the specific-sex chromosome, being more dynamic than traditionally appreciated.

Next generation sequencing-based analysis of repetitive DNA in the model dioecious [corrected] plant Silene latifolia

Background

Silene latifolia is a dioceous plant with well distinguished X and Y chromosomes that is used as a model to study sex determination and sex chromosome evolution in plants. However, efficient utilization of this species has been hampered by the lack of large-scale sequencing resources and detailed analysis of its genome composition, especially with respect to repetitive DNA, which makes up the majority of the genome.

Methodology/Principal Findings

We performed low-pass 454 sequencing followed by similarity-based clustering of 454 reads in order to identify and characterize sequences of all major groups of S. latifolia repeats. Illumina sequencing data from male and female genomes were also generated and employed to quantify the genomic proportions of individual repeat families. The majority of identified repeats belonged to LTR-retrotransposons, constituting about 50% of genomic DNA, with Ty3/gypsy elements being more frequent than Ty1/copia. While there were differences between the male and female genome in the abundance of several repeat families, their overall repeat composition was highly similar. Specific localization patterns on sex chromosomes were found for several satellite repeats using in situ hybridization with probes based on k-mer frequency analysis of Illumina sequencing data.

Conclusions/Significance

This study provides comprehensive information about the sequence composition and abundance of repeats representing over 60% of the S. latifolia genome. The results revealed generally low divergence in repeat composition between the sex chromosomes, which is consistent with their relatively recent origin. In addition, the study generated various data resources that are available for future exploration of the S. latifolia genome.

Microsatellite distribution on sex chromosomes at different stages of heteromorphism and heterochromatinization in two lizard species (Squamata: Eublepharidae: Coleonyx elegans and lacertidae: Eremias velox)

BACKGROUND

The accumulation of repetitive sequences such as microsatellites during the differentiation of sex chromosomes has not been studied in most squamate reptiles (lizards, amphisbaenians and snakes), a group which has a large diversity of sex determining systems. It is known that the Bkm repeats containing tandem arrays of GATA tetranucleotides are highly accumulated on the degenerated W chromosomes in advanced snakes. Similar, potentially homologous, repetitive sequences were found on sex chromosomes in other vertebrates. Using FISH with probes containing all possible mono-, di-, and tri-nucleotide sequences and GATA, we studied the genome distribution of microsatellite repeats on sex chromosomes in two lizard species (the gecko Coleonyx elegans and the lacertid Eremias velox) with independently evolved sex chromosomes. The gecko possesses heteromorphic euchromatic sex chromosomes, while sex chromosomes in the lacertid are homomorphic and the W chromosome is highly heterochromatic. Our aim was to test whether microsatellite distribution on sex chromosomes corresponds to the stage of their heteromorphism or heterochromatinization. Moreover, because the lizards lie phylogenetically between snakes and other vertebrates with the Bkm-related repeats on sex chromosomes, the knowledge of their repetitive sequence is informative for the determination of conserved versus convergently evolved repetitive sequences across vertebrate lineages.

RESULTS

Heteromorphic sex chromosomes of C. elegans do not show any sign of microsatellite accumulation. On the other hand, in E. velox, certain microsatellite sequences are extensively accumulated over the whole length or parts of the W chromosome, while others, including GATA, are absent on this heterochromatinized sex chromosome.

CONCLUSION

The accumulation of microsatellite repeats corresponds to the stage of heterochromatinization of sex chromosomes rather than to their heteromorphism. The lack of GATA repeats on the sex chromosomes of both lizards suggests that the Bkm-related repeats on sex chromosomes in snakes and other vertebrates evolved convergently. The comparison of microsatellite sequences accumulated on sex chromosomes in E. velox and in other eukaryotic organisms suggests that historical contingency, not characteristics of particular sequences, plays a major role in the determination of which microsatellite sequence is accumulated on the sex chromosomes in a particular lineage.

Factors influencing resistance of UV-irradiated DNA to the restriction endonuclease cleavage

DNA molecules of pUC19, pBR322 and PhiX174 were irradiated by various doses of UV light and the irradiated molecules were cleaved by about two dozen type II restrictases. The irradiation generally blocked the cleavage in a dose-dependent way. In accordance with previous studies, the (A + T)-richness and the (PyPy) dimer content of the restriction site belongs among the factors that on average, cause an increase in the resistance of UV damaged DNA to the restrictase cleavage. However, we observed strong effects of UV irradiation even with (G + C)-rich and (PyPy)-poor sites. In addition, sequences flanking the restriction site influenced the protection in some cases (e.g. HindIII), but not in others (e.g. SalI), whereas neoschizomer couples SmaI and AvaI, or SacI and Ecl136II, cleaved the UV-irradiated DNA similarly. Hence the intrastrand thymine dimers located in the recognition site are not the only photoproduct blocking the restrictases. UV irradiation of the A-form generally made the irradiated DNA less resistant to restrictase cleavage than irradiation in the B-form and in some cases, the A-form completely protected the UV-irradiated DNA against the damage recognized by the restrictases. The present results also demonstrate that the UV irradiation approach used to generate partial digests in genomic DNA studies, can be extended to the (G + C)-rich and (PyPy)-poor restriction sites. The present extensive and quantitative data can be used in genomic applications of UV damage probing by restrictases.

Heterogeneity of rDNA distribution and genome size in Silene spp

Genus Silene L. (Caryophyllaceae) contains about 700 species divided into 44 sections. According to recent taxonomic classification this genus also includes taxa previously classified in genera Lychnis and Melandrium. In this work, four Silene species belonging to different sections were studied: S. latifolia (syn. Melandrium album, Section Elisanthe), S. vulgaris (Inflatae), S. pendula (Erectorefractae), and S. chalcedonica (syn. Lychnis chalcedonica, Lychnidiformes). Flow cytometric analysis revealed a genome size of 2.25 and 2.35 pg/2C for S. vulgaris and S. pendula and of 5.73 and 6.59 pg/2C for S. latifolia and S. chalcedonica. All four species have the same chromosome number including the pair of sex chromosomes of the dioecious S. latifolia (2n = 2x = 24). Double target fluorescence in-situ hybridization revealed the chromosomal locations of 25S rDNA and 5S rDNA. A marked variation in number and localization of rDNA loci but no correlation between the numbers of rDNA clusters and genome size was found. FISH and genome size data indicate that nuclear genomes of Silene species are highly diversified as a result of numerous DNA amplifications and translocations.

Localization of male-specifically expressed MROS genes of Silene latifolia by PCR on flow-sorted sex chromosomes and autosomes

The dioecious white campion Silene latifolia (syn. Melandrium album) has heteromorphic sex chromosomes, XX in females and XY in males, that are larger than the autosomes and enable their separation by flow sorting. The group of MROS genes, the first male-specifically expressed genes in dioecious plants, was recently identified in S. latifolia. To localize the MROS genes, we used the flow-sorted X chromosomes and autosomes as a template for PCR with internal primers. Our results indicate that the MROS3 gene is located in at least two copies tandemly arranged on the X chromosome with additional copy(ies) on the autosome(s), while MROS1, MROS2, and MROS4 are exclusively autosomal. The specificity of PCR products was checked by digestion with a restriction enzyme or reamplification using nested primers. Homology search of databases has shown the presence of five MROS3 homologues in A. thaliana, four of them arranged in two tandems, each consisting of two copies. We conclude that MROS3 is a low-copy gene family, connected with the proper pollen development, which is present not only in dioecious but also in other dicot plant species.

DNA extraction by zinc

A fast, very simple and efficient method of DNA extraction is described which takes advantage of DNA sedimentation induced by millimolar concentrations of ZnCl2. The zinc-induced sedimentation is furthermore strongly promoted by submillimolar phosphate anion concentrations. Within <30 min, the method recovers >90% of DNA irrespective of whether a plasmid DNA or short oligonucleotides are the extracted material. The method works with plasmid DNA and oligonucleotide concentrations as low as 100 ng/ml and 10 microg/ml, respectively, without using any expensive facilities or toxic chemicals.

UV light-induced crosslinking of short DNA duplex strands: nucleotide sequence preferences and a prominent role of the duplex ends

By various doses of UV light, we irradiated 32 short DNA duplexes having between 12 and 40 nucleotide residues in length, and monitored the induced crosslinks between the complementary DNA strands by denaturing polyacrylamide gel electrophoresis. The experiments revealed that the crosslinking was strongest with the alternating sequence of T and A and weaker with the alternating sequence of T and G (C and A in the complementary strand). On the other hand, GC blocks of any sequence provided undetectable amounts of interstrand crosslinks even at the highest doses of UV irradiation. The amount of crosslinked strands logarithmically increased with the UV dose but it did not depend on the oligonucleotide concentration, ionic strength, divalent magnesium or manganese cations and pH at least within the examined regions of the experimental conditions, unless the oligonucleotide denatured or isomerized into a unimolecular foldback. The extent of crosslinking also did not depend on the (dT-dA)n duplex length to indicate that the crosslink was predominantly localized at a specific duplex locus. Experiments with (dT-dA)8 "mutants" in which AT pairs were systematically replaced by GC pairs at various molecule positions, revealed that the crosslinking predominantly occurred at the oligo(dT-dA) duplex ends. The crosslinking is a direct method to detect duplexes of DNA, which is here, for example, demonstrated with the heteroduplex of (dT-dA)12 and (dT-dA)16.

Tris buffer protects DNA backbone against breakage upon irradiation with ultraviolet light

We show that Tris molecules protect DNA against nicking upon irradiation with ultraviolet light. However, the protective effect only concerns DNA backbone but not bases and it is observed in aqueous solution but not in formamide. Changes of pH or ionic strength due to Tris have no effect on the protection. The present observation has a practical importance for photofootprinting studies of DNA and its complexes with proteins but it can also serve as a basis for a development of a novel method reflecting DNA hydration and conformation.