Cytological characterization of the tandem repetitive sequences and their methylation status in the Antirrhinum majus genome

Tandem repetitive sequences are DNA motifs common in the genomes of eukaryotic species and are often embedded in heterochromatic regions. In most eukaryotes, ribosomal genes, as well as centromeres and telomeres or subtelomeres, are associated with abundant tandem arrays of repetitive sequences and typically represent the final barriers to completion of whole-genome sequencing. The nature of these repeats makes it difficult to estimate their actual sizes. In this study, combining the two cytological techniques DNA fiber-FISH and pachytene chromosome FISH allowed us to characterize the tandem repeats distributed genome wide in Antirrhinum majus and identify four types of tandem repeats, 45S rDNA, 5S rDNA, CentA1, and CentA2, representing the major tandem repetitive components, which were estimated to have a total length of 18.50 Mb and account for 3.59% of the A. majus genome. FISH examination revealed that all the tandem repeats correspond to heterochromatic knobs along the pachytene chromosomes. Moreover, the methylation status of the tandem repeats was investigated in both somatic cells and pollen mother cells from anther tissues using an antibody against 5-methylcytosine combined with sequential FISH analyses. Our results showed that these repeats were hypomethylated in anther tissues, especially in the pollen mother cells at pachytene stage.     

DNA methylation of retrotransposons,DNA transposons and genes in sugar beet (Beta vulgaris L.)

The methylation of cytosines shapes the epigenetic landscape of plant genomes, coordinates transgenerational epigenetic inheritance, represses the activity of transposable elements (TEs), affects gene expression and, hence, can influence the phenotype. Sugar beet (Beta vulgaris ssp. vulgaris), an important crop that accounts for 30% of worldwide sugar needs, has a relatively small genome size (758 Mbp) consisting of approximately 485 Mbp repetitive DNA (64%), in particular satellite DNA, retrotransposons and DNA transposons. Genome‐wide cytosine methylation in the sugar beet genome was studied in leaves and leaf‐derived callus with a focus on repetitive sequences, including retrotransposons and DNA transposons, the major groups of repetitive DNA sequences, and compared with gene methylation. Genes showed a specific methylation pattern for CG, CHG (H = A, C, and T) and CHH sites, whereas the TE pattern differed, depending on the TE class (class 1, retrotransposons and class 2, DNA transposons). Along genes and TEs, CG and CHG methylation was higher than that of adjacent genomic regions. In contrast to the relatively low CHH methylation in retrotransposons and genes, the level of CHH methylation in DNA transposons was strongly increased, pointing to a functional role of asymmetric methylation in DNA transposon silencing. Comparison of genome‐wide DNA methylation between sugar beet leaves and callus revealed a differential methylation upon tissue culture. Potential epialleles were hypomethylated (lower methylation) at CG and CHG sites in retrotransposons and genes and hypermethylated (higher methylation) at CHH sites in DNA transposons of callus when compared with leaves.     

The holocentric species Luzula elegans shows interplay between centromere and large‐scale genome organization

In higher plants, the large‐scale structure of monocentric chromosomes consists of distinguishable eu‐ and heterochromatic regions, the proportions and organization of which depend on a species' genome size. To determine whether the same interplay is maintained for holocentric chromosomes, we investigated the distribution of repetitive sequences and epigenetic marks in the woodrush Luzula elegans (3.81 Gbp/1C). Sixty‐one per cent of the L. elegans genome is characterized by highly repetitive DNA, with over 30 distinct sequence families encoding an exceptionally high diversity of satellite repeats. Over 33% of the genome is composed of the Angela clade of Ty1/copia LTR retrotransposons, which are uniformly dispersed along the chromosomes, while the satellite repeats occur as bands whose distribution appears to be biased towards the chromosome termini. No satellite showed an almost chromosome‐wide distribution pattern as expected for a holocentric chromosome and no typical centromere‐associated LTR retrotransposons were found either. No distinguishable large‐scale patterns of eu‐ and heterochromatin‐typical epigenetic marks or early/late DNA replicating domains were found along mitotic chromosomes, although super‐high‐resolution light microscopy revealed distinguishable interspersed units of various chromatin types. Our data suggest a correlation between the centromere and overall genome organization in species with holocentric chromosomes.     

Rapid proliferation and nucleolar organizer targeting centromeric retrotransposons in cotton

Centromeric chromatin in most eukaryotes is composed of highly repetitive centromeric retrotransposons and satellite repeats that are highly variable even among closely related species. The evolutionary mechanisms that underlie the rapid evolution of centromeric repeats remain unknown. To obtain insight into the evolution of centromeric repeats following polyploidy, we studied a model diploid progenitor (Gossypium raimondii, D‐genome) of the allopolyploid (AD‐genome) cottons, G. hirsutum and G. barbadense. Sequence analysis of chromatin‐immunoprecipitated DNA showed that the G. raimondii centromeric repeats originated from retrotransposon‐related sequences. Comparative analysis showed that nine of the 10 analyzed centromeric repeats were absent from the centromeres in the A‐genome and related diploid species (B‐, F‐ and G‐genomes), indicating that they colonized the centromeres of D‐genome lineage after the divergence of the A‐ and D‐ ancestral species or that they were ancestrally retained prior to the origin of Gossypium. Notably, six of the nine repeats were present in both the A‐ and D‐subgenomes in tetraploid G. hirsutum, and increased in abundance in both subgenomes. This finding suggests that centromeric repeats may spread and proliferate between genomes subsequent to polyploidization. Two repeats, Gr334 and Gr359 occurred in both the centromeres and nucleolar organizer regions (NORs) in D‐ and AD‐genome species, yet localized to just the NORs in A‐, B‐, F‐, and G‐genome species. Contained within is a story of an established centromeric repeat that is eliminated and allopolyploidization provides an opportunity for reinvasion and reestablishment, which broadens our evolutionary understanding behind the cycles of centromeric repeat establishment and targeting.     

Diploidization and genome size change in allopolyploids is associated with differential dynamics of low‐ and high‐copy sequences

Recent advances have highlighted the ubiquity of whole‐genome duplication (polyploidy) in angiosperms, although subsequent genome size change and diploidization (returning to a diploid‐like condition) are poorly understood. An excellent system to assess these processes is provided by Nicotiana section Repandae, which arose via allopolyploidy (approximately 5 million years ago) involving relatives of Nicotiana sylvestris and Nicotiana obtusifolia. Subsequent speciation in Repandae has resulted in allotetraploids with divergent genome sizes, including Nicotiana repanda and Nicotiana nudicaulis studied here, which have an estimated 23.6% genome expansion and 19.2% genome contraction from the early polyploid, respectively. Graph‐based clustering of next‐generation sequence data enabled assessment of the global genome composition of these allotetraploids and their diploid progenitors. Unexpectedly, in both allotetraploids, over 85% of sequence clusters (repetitive DNA families) had a lower abundance than predicted from their diploid relatives; a trend seen particularly in low‐copy repeats. The loss of high‐copy sequences predominantly accounts for the genome downsizing in N. nudicaulis. In contrast, N. repanda shows expansion of clusters already inherited in high copy number (mostly chromovirus‐like Ty3/Gypsy retroelements and some low‐complexity sequences), leading to much of the genome upsizing predicted. We suggest that the differential dynamics of low‐ and high‐copy sequences reveal two genomic processes that occur subsequent to allopolyploidy. The loss of low‐copy sequences, common to both allopolyploids, may reflect genome diploidization, a process that also involves loss of duplicate copies of genes and upstream regulators. In contrast, genome size divergence between allopolyploids is manifested through differential accumulation and/or deletion of high‐copy‐number sequences.     

Maintenance of the induced hypomethylated state of tobacco nuclear repetitive DNA sequences in the course of protoplast and plant regeneration

Previously, we established experimental conditions allowing us to induce hypomethylation of tandem arrays of highly repetitive DNA sequences in the Nicotiana tabacum L. nuclear genome (M. Bezdk et al. 1991, Planta 184, 487–490). In this paper, we demonstrate that loci containing the highly repetitive sequences of the HRS60 family can maintain the induced hypomethylated state through protoplast regeneration, non-differentiated callus growth, and plant regeneration. The hypomethylation, induced with 5-azacytidine and monitored on these sequences, did not substantially alter the capacity of calli to produce plants. It can be concluded that, in contradistinction of multiple copies of transgenes or ectopic genes which are usually recognized as methylation targets, endogenous tandem repeats, such as the HRS60, present at 10⁵ copies in the genome, can escape de-novo methylation.Abbreviation AzaC 5-azacytidine This project was supported by the Grant Agency of the Academy of Sciences of the Czech Republic. We thank Ms. Libue Jedliková, Ms. Hana Suchánková, and Ms. Emilie Koudelková for technical assistance.     

Centromere and telomere sequence alterations reflect the rapid genome evolution within the carnivorous plant genus Genlisea

Linear chromosomes of eukaryotic organisms invariably possess centromeres and telomeres to ensure proper chromosome segregation during nuclear divisions and to protect the chromosome ends from deterioration and fusion, respectively. While centromeric sequences may differ between species, with arrays of tandemly repeated sequences and retrotransposons being the most abundant sequence types in plant centromeres, telomeric sequences are usually highly conserved among plants and other organisms. The genome size of the carnivorous genus Genlisea (Lentibulariaceae) is highly variable. Here we study evolutionary sequence plasticity of these chromosomal domains at an intrageneric level. We show that Genlisea nigrocaulis (1C = 86 Mbp; 2n = 40) and G. hispidula (1C = 1550 Mbp; 2n = 40) differ as to their DNA composition at centromeres and telomeres. G. nigrocaulis and its close relative G. pygmaea revealed mainly 161 bp tandem repeats, while G. hispidula and its close relative G. subglabra displayed a combination of four retroelements at centromeric positions. G. nigrocaulis and G. pygmaea chromosome ends are characterized by the Arabidopsis‐type telomeric repeats (TTTAGGG); G. hispidula and G. subglabra instead revealed two intermingled sequence variants (TTCAGG and TTTCAGG). These differences in centromeric and, surprisingly, also in telomeric DNA sequences, uncovered between groups with on average a > 9‐fold genome size difference, emphasize the fast genome evolution within this genus. Such intrageneric evolutionary alteration of telomeric repeats with cytosine in the guanine‐rich strand, not yet known for plants, might impact the epigenetic telomere chromatin modification.     

Hypomethylation of classical satellite DNA and chromosome instability in lymphoblastoid cell lines

To determine possible relationships between DNA hypomethylation and chromosome instability, human lymphoblastoid cell lines from different genetic constitutions were studied with regard to 1) uncoiling and rearrangements, which preferentially affect the heterochromatic segments of chromosomes 1 and 16; 2) the methylation status of the tandemly repetitive sequences (classical satellite and alphoid DNAs) from chromosomes 1 and 16, and of the L1Hs interspersed repetitive sequences. The methylation status largely varied from cell line to cell line, but for a given cell line, the degree of methylation was similar for all the repetitive DNAs studied. Two cell lines, one obtained from a Fanconi anemia patient and the other from an ataxia telangiectasia patient were found to be heavily hypomethylated. The heterochromatic segments of their chromosomes 1 and 16 were more frequently elongated and rearranged than those from other cell lines, which were found to be less hypomethylated. Thus, in these lymphoblastoid cell lines, alterations characterized by uncoiling and rearrangements of heterochromatic segments from chromosomes 1 and 16 seem to correlate with the hypomethylation of their repetitive DNAs. Two-color in situ hybridizations demonstrated that these elongations and rearrangements involved only classical satellite-DNA-containing heterochromatin. This specificity may be related to the excess of breakages affecting the chromosomes carrying these structures in a variety of pathological conditions.     

Sugarcane genome sequencing by methylation filtration provides tools for genomic research in the genus Saccharum

Many economically important crops have large and complex genomes that hamper their sequencing by standard methods such as whole genome shotgun (WGS). Large tracts of methylated repeats occur in plant genomes that are interspersed by hypomethylated gene‐rich regions. Gene‐enrichment strategies based on methylation profiles offer an alternative to sequencing repetitive genomes. Here, we have applied methyl filtration with McrBC endonuclease digestion to enrich for euchromatic regions in the sugarcane genome. To verify the efficiency of methylation filtration and the assembly quality of sequences submitted to gene‐enrichment strategy, we have compared assemblies using methyl‐filtered (MF) and unfiltered (UF) libraries. The use of methy filtration allowed a better assembly by filtering out 35% of the sugarcane genome and by producing 1.5× more scaffolds and 1.7× more assembled Mb in length compared with unfiltered dataset. The coverage of sorghum coding sequences (CDS) by MF scaffolds was at least 36% higher than by the use of UF scaffolds. Using MF technology, we increased by 134× the coverage of gene regions of the monoploid sugarcane genome. The MF reads assembled into scaffolds that covered all genes of the sugarcane bacterial artificial chromosomes (BACs), 97.2% of sugarcane expressed sequence tags (ESTs), 92.7% of sugarcane RNA‐seq reads and 98.4% of sorghum protein sequences. Analysis of MF scaffolds from encoded enzymes of the sucrose/starch pathway discovered 291 single‐nucleotide polymorphisms (SNPs) in the wild sugarcane species, S. spontaneum and S. officinarum. A large number of microRNA genes was also identified in the MF scaffolds. The information achieved by the MF dataset provides a valuable tool for genomic research in the genus Saccharum and for improvement of sugarcane as a biofuel crop.     

DNA regions flanking the major Arabidopsis thaliana satellite are principally enriched in Athila retroelement sequences

An analysis of Arabidopsis thaliana heterochromatic regions revealed that genomic sequences immediately flanking the major 180 bp satellite are essentially made of middle repetitive sequences and that most of these sequences correspond to defective Athila retroelements. Using YAC and clones, we evaluated the distribution of Athila elements in the Arabidopsis genome and showed that, despite the presence of numerous euchromatic copies, these elements are especially concentrated in or near heterochromatic regions. Sequencing of the various DNA transitions between satellite and Athila repeats provides strong evidence that most of the heterochromatic elements retrotransposed directly into 180 bp satellite clusters.     

Hypomethylated sequences: Characterization of the duplicate soybean genome

Soybean is believed to be a diploidized tetraploid generated from an allotetraploid ancestor. In this study, we used hypomethylated genomic DNA as a source of probes to investigate the genomic structure and methylation patterns of duplicated sequences. Forty-five genomic clones from Phaseolus vulgaris and 664 genomic clones from Glycine max were used to examine the duplicated regions in the soybean genome. Southern analysis of genomic DNA using probes from both sources revealed that greater than 15% of the hypomethylated genomic regions were only present once in the soybean genome. The remaining ca. 85% of the hypomethylated regions comprise duplicated or middle repetitive DNA sequences. If only the ratio of single to duplicate probe patterns is considered, it appears that 25% of the single-copy sequences have been lost. By using a subset of probes that only detected duplicated sequences, we examined the methylation status of the homeologous genomes with the restriction enzymes MspI and HpaII. We found that in all cases both copies of these regions were hypomethylated, although there were examples of low-level methylation. It appears that duplicate sequences are being eliminated in the diploidization process. Our data reveal no evidence that duplicated sequences are being silenced by inactivation correlated with methylation patterns.     

Analysis of transposable elements and organellar DNA in male and female genomes of a species with a huge Y chromosome reveals distinct Y centromeres

Few angiosperms have distinct Y chromosomes. Among those that do are Silene latifolia (Caryophyllaceae), Rumex acetosa (Polygonaceae) and Coccinia grandis (Cucurbitaceae), the latter having a male/female difference of 10% of the total genome (female individuals have a 0.85 pg genome, male individuals 0.94 pg), due to a Y chromosome that arose about 3 million years ago. We compared the sequence composition of male and female C. grandis plants and determined the chromosomal distribution of repetitive and organellar DNA with probes developed from 21 types of repetitive DNA, including 16 mobile elements. The size of the Y chromosome is largely due to the accumulation of certain repeats, such as members of the Ty1/copia and Ty3/gypsy superfamilies, an unclassified element and a satellite, but also plastome‐ and chondriome‐derived sequences. An abundant tandem repeat with a unit size of 144 bp stains the centromeres of the X chromosome and the autosomes, but is absent from the Y centromere. Immunostaining with pericentromere‐specific markers for anti‐histone H3Ser10ph and H2AThr120ph revealed a Y‐specific extension of these histone marks. That the Y centromere has a different make‐up from all the remaining centromeres raises questions about its spindle attachment, and suggests that centromeric or pericentromeric chromatin might be involved in the suppression of recombination.     

Heterochromatic regions in different Drosophila melanogaster stocks contain similar arrangements of moderate repeats with inserted copia-like elements (MDG1)

Seven out of twenty 30–50 kb genome fragments with an MDG1 copia-like element cloned in cosmids were found to carry homologous sequences which belong to a new family of non-mobile heterochromatic moderate repeats (the HMR family). These repeats along with the MDG1 copies inserted in them are under-replicated in polytene chromosomes. Such repeats may also be located in the intercalary heterochromatin site 12E of the X chromosome. Chromosomal heterochromatic regions are enriched with one of the two main genomic variants of MDG1, MDG1^het, identifiable by EcoRI restriction. From Southern DNA blot analysis the number of MDG1^het copies and their sites within the heterochromatin are invariant in all the stocks examined, while there is not a single MDG1 site along the polytene chromosomes shared by all the stocks in question.     

Evolution and selection of Rhg1, a copy‐number variant nematode‐resistance locus

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2‐kb unit contains four genes. One allele of Rhg1, Rhg1‐b, is responsible for protecting most US soybean production from SCN. Whole‐genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2‐kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high‐density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non‐neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.     

Chromosome‐level genome assembly of the razor clam Sinonovacula constricta (Lamarck, 1818)

Bivalves, a highly diverse and the most evolutionarily successful class of invertebrates native to aquatic habitats, provide valuable molecular resources for understanding the evolutionary adaptation and aquatic ecology. Here, we reported a high‐quality chromosome‐level genome assembly of the razor clam Sinonovacula constricta using Pacific Bioscience single‐molecule real‐time sequencing, Illumina paired‐end sequencing, 10X Genomics linked‐reads and Hi‐C reads. The genome size was 1,220.85 Mb, containing scaffold N50 of 65.93 Mb and contig N50 of 976.94 Kb. A total of 899 complete (91.92%) and seven partial (0.72%) matches of the 978 metazoa Benchmarking Universal Single‐Copy Orthologs were determined in this genome assembly. And Hi‐C scaffolding of the genome resulted in 19 pseudochromosomes. A total of 28,594 protein‐coding genes were predicted in the S. constricta genome, of which 25,413 genes (88.88%) were functionally annotated. In addition, 39.79% of the assembled genome was composed of repetitive sequences, and 4,372 noncoding RNAs were identified. The enrichment analyses of the significantly expanded and contracted genes suggested an evolutionary adaptation of S. constricta to highly stressful living environments. In summary, the genomic resources generated in this work not only provide a valuable reference genome for investigating the molecular mechanisms of S. constricta biological functions and evolutionary adaptation, but also facilitate its genetic improvement and disease treatment. Meanwhile, the obtained genome greatly improves our understanding of the genetics of molluscs and their comparative evolution.     

Chromosome‐level genome assembly of Triplophysa tibetana,a fish adapted to the harsh high‐altitude environment of the Tibetan Plateau

Triplophysa is an endemic fish genus of the Tibetan Plateau in China. Triplophysa tibetana, which lives at a recorded altitude of ~4,000 m and plays an important role in the highland aquatic ecosystem, serves as an excellent model for investigating high‐altitude environmental adaptation. However, evolutionary and conservation studies of T. tibetana have been limited by scarce genomic resources for the genus Triplophysa. In the present study, we applied PacBio sequencing and the Hi‐C technique to assemble the T. tibetana genome. A 652‐Mb genome with 1,325 contigs with an N50 length of 3.1 Mb was obtained. The 1,137 contigs were further assembled into 25 chromosomes, representing 98.7% and 80.47% of all contigs at the base and sequence number level, respectively. Approximately 260 Mb of sequence, accounting for ~39.8% of the genome, was identified as repetitive elements. DNA transposons (16.3%), long interspersed nuclear elements (12.4%) and long terminal repeats (11.0%) were the most repetitive types. In total, 24,372 protein‐coding genes were predicted in the genome, and ~95% of the genes were functionally annotated via a search in public databases. Using whole genome sequence information, we found that T. tibetana diverged from its common ancestor with Danio rerio ~121.4 million years ago. The high‐quality genome assembled in this work not only provides a valuable genomic resource for future population and conservation studies of T. tibetana, but it also lays a solid foundation for further investigation into the mechanisms of environmental adaptation of endemic fishes in the Tibetan Plateau.