Origin, Diversity, and Evolution of Telomere Sequences in Plants

Analysis of Telomerase RNA Structure in Physcomitrium patens Indicates Functionally Relevant Transitions Between OPEN and CLOSED Conformations

Telomerase RNA (TR) conformation determines its function as a template for telomere synthesis and as a scaffold for the assembly of the telomerase nucleoprotein complex. Experimental analyses of TR secondary structure using DMS-Map Seq and SHAPE-Map Seq techniques show its CLOSED conformation as the consensus structure where the template region cannot perform its function. Our data show that the apparent discrepancy between experimental results and predicted TR functional conformation, mostly ignored in published studies, can be explained using data analysis based on single-molecule structure prediction from individual sequencing reads by the recently established DaVinci method. This method results in several clusters of secondary structures reflecting the structural dynamics of TR, possibly related to its multiple functional states. Interestingly, the presumed active (OPEN) conformation of TR corresponds to a minor fraction of TR under in vivo conditions. Therefore, structural polymorphism and dynamic TR transitions between CLOSED and OPEN conformations may be involved in telomerase activity regulation as a switch that functions independently of total TR transcript levels.

Leveraging whole-genome sequencing to estimate telomere length in plants

Changes in telomere length are increasingly used to indicate species' response to environmental stress across diverse taxa. Despite this broad use, few studies have explored telomere length in plants. Thus, evaluation of new approaches for measuring telomeres in plants is needed. Rapid advances in sequencing approaches and bioinformatic tools now allow estimation of telomere content from whole-genome sequencing (WGS) data, a proxy for telomere length. While telomere content has been quantified extensively using quantitative polymerase chain reaction (qPCR) and WGS in humans, no study to date has compared the effectiveness of WGS in estimating telomere length in plants relative to qPCR approaches. In this study, we use 100 Populus clones re-sequenced using short-read Illumina sequencing to quantify telomere length comparing three different bioinformatic approaches (Computel, K-seek and TRIP) in addition to qPCR. Overall, telomere length estimates varied across different bioinformatic approaches, but were highly correlated across methods for individual genotypes. A positive correlation was observed between WGS estimates and qPCR, however, Computel estimates exhibited the greatest correlation. Computel incorporates genome coverage into telomere length calculations, suggesting that genome coverage is likely important to telomere length quantification when using WGS data. Overall, telomere estimates from WGS provided greater precision and accuracy of telomere length estimates relative to qPCR. The findings suggest WGS is a promising approach for assessing telomere length and, as the field of telomere ecology evolves, may provide added value to assaying response to biotic and abiotic environments for plants needed to accelerate plant breeding and conservation management.

TeloBase: a community-curated database of telomere sequences across the tree of life

Discoveries over the recent decade have demonstrated the unexpected diversity of telomere DNA motifs in nature. However, currently available resources, 'Telomerase database' and 'Plant rDNA database', contain just fragments of all relevant literature published over decades of telomere research as they have a different primary focus and limited updates. To fill this gap, we gathered data about telomere DNA sequences from a thorough literature screen as well as by analysing publicly available NGS data, and we created TeloBase (http://cfb.ceitec.muni.cz/telobase/) as a comprehensive database of information about telomere motif diversity. TeloBase is supplemented by internal taxonomy utilizing popular on-line taxonomic resources that enables in-house data filtration and graphical visualisation of telomere DNA evolutionary dynamics in the form of heat tree plots. TeloBase avoids overreliance on administrators for future data updates by having a simple form and community-curation system for application and approval, respectively, of new telomere sequences by users, which should ensure timeliness of the database and topicality. To demonstrate TeloBase utility, we examined telomere motif diversity in species from the fungal genus Aspergillus, and discovered (TTTATTAGGG)n sequence as a putative telomere motif in the plant family Chrysobalanaceae. This was bioinformatically confirmed by analysing template regions of identified telomerase RNAs.

Telomere-to-telomere and gap-free genome assembly of a susceptible grapevine species (Thompson Seedless) to facilitate grape functional genomics

Grapes are globally recognized as economically significant fruit trees. Among grape varieties, Thompson Seedless holds paramount influence for fresh consumption and for extensive applications in winemaking, drying, and juicing. This variety is one of the most efficient genotypes for grape genetic modification. However, the lack of a high-quality genome has impeded effective breeding efforts. Here, we present the high-quality reference genome of Thompson Seedless with all 19 chromosomes represented as 19 contiguous sequences (N50 = 27.1 Mb) with zero gaps and prediction of all telomeres and centromeres. Compared with the previous assembly (TSv1 version), the new assembly incorporates an additional 31.5 Mb of high-quality sequenced data with annotation of a total of 30 397 protein-coding genes. We also performed a meticulous analysis to identify nucleotide-binding leucine-rich repeat genes (NLRs) in Thompson Seedless and two wild grape varieties renowned for their disease resistance. Our analysis revealed a significant reduction in the number of two types of NLRs, TIR-NB-LRR (TNL) and CC-NB-LRR (CNL), in Thompson Seedless, which may have led to its sensitivity to many fungal diseases, such as powdery mildew, and an increase in the number of a third type, RPW8 (resistance to powdery mildew 8)-NB-LRR (RNL). Subsequently, transcriptome analysis showed significant enrichment of NLRs during powdery mildew infection, emphasizing the pivotal role of these elements in grapevine's defense against powdery mildew. The successful assembly of a high-quality Thompson Seedless reference genome significantly contributes to grape genomics research, providing insight into the importance of seedlessness, disease resistance, and color traits, and these data can be used to facilitate grape molecular breeding efforts.

Transposition, duplication, and divergence of the telomerase RNA underlies the evolution of Mimulus telomeres

Telomeres are nucleoprotein complexes with a crucial role of protecting chromosome ends. It consists of simple repeat sequences and dedicated telomere-binding proteins. Because of its vital functions, components of the telomere, for example its sequence, should be under strong evolutionary constraint. But across all plants, telomere sequences display a range of variation and the evolutionary mechanism driving this diversification is largely unknown. Here, we discovered in Monkeyflower (Mimulus) the telomere sequence is even variable between species. We investigated the basis of Mimulus telomere sequence evolution by studying the long noncoding telomerase RNA (TR), which is a core component of the telomere maintenance complex and determines the telomere sequence. We conducted total RNA-based de novo transcriptomics from 16 Mimulus species and analyzed reference genomes from 6 species, and discovered Mimulus species have evolved at least three different telomere sequences: (AAACCCT)n, (AAACCCG)n, and (AAACCG)n. Unexpectedly, we discovered several species with TR duplications and the paralogs had functional consequences that could influence telomere evolution. For instance, M. lewisii had two sequence-divergent TR paralogs and synthesized a telomere with sequence heterogeneity, consisting of AAACCG and AAACCCG repeats. Evolutionary analysis of the M. lewisii TR paralogs indicated it had arisen from a transposition-mediate duplication process. Further analysis of the TR from multiple Mimulus species showed the gene had frequently transposed and inserted into new chromosomal positions during Mimulus evolution. From our results, we propose the TR transposition, duplication, and divergence model to explain the evolutionary sequence turnovers in Mimulus and potentially all plant telomeres.

Telomeric repeat evolution in the phylum Nematoda revealed by high-quality genome assemblies and subtelomere structures

Telomeres are composed of tandem arrays of telomeric-repeat motifs (TRMs) and telomere-binding proteins (TBPs), which are responsible for ensuring end-protection and end-replication of chromosomes. TRMs are highly conserved owing to the sequence specificity of TBPs, although significant alterations in TRM have been observed in several taxa, except Nematoda. We used public whole-genome sequencing data sets to analyze putative TRMs of 100 nematode species and determined that three distinct branches included specific novel TRMs, suggesting that evolutionary alterations in TRMs occurred in Nematoda. We focused on one of the three branches, the Panagrolaimidae family, and performed a de novo assembly of four high-quality draft genomes of the canonical (TTAGGC) and novel TRM (TTAGAC) isolates; the latter genomes revealed densely clustered arrays of the novel TRM. We then comprehensively analyzed the subtelomeric regions of the genomes to infer how the novel TRM evolved. We identified DNA damage-repair signatures in subtelomeric sequences that were representative of consequences of telomere maintenance mechanisms by alternative lengthening of telomeres. We propose a hypothetical scenario in which TTAGAC-containing units are clustered in subtelomeric regions and pre-existing TBPs capable of binding both canonical and novel TRMs aided the evolution of the novel TRM in the Panagrolaimidae family.

Extracellular oil production by Rhodotorula paludigena BS15 for biorefinery without complex downstream processes

To realize biomass refinery without complex downstream processes, we extensively screened for microbial strains that efficiently produce extracellular oil from sugars. Rhodotorula paludigena (formerly Rhodosporidium paludigenum) BS15 was found to efficiently produce polyol esters of fatty acids (PEFAs), which mainly comprised of 3-acetoxypalmitic acid and partially acetylated mannitol/arabinitol. To evaluate the performance of this strain, fed-batch fermentation was demonstrated on a flask scale, and 110 g/L PEFA and 103 g/L dry cells were produced in 12 days. To the best of our knowledge, the strain BS15 exhibited the highest PEFA titer (g/L) ever to be reported so far. Because the PEFA precipitated at the bottom of the culture broth, it could be easily recovered by simply discarding the upper phase. Various carbon sources can be utilized for cell growth and/or PEFA production, which signifies the potential for converting diverse biomass sources. Two different types of next-generation sequencers, Illumina HiSeq and Oxford Nanopore PromethION, were used to analyze the whole-genome sequence of the strain BS15. The integrative data analysis generated a high-quality and reliable reference genome for PEFA-producing R. paludigena. The 22.5-M base genome sequence and the estimated genes were registered in Genbank (accession numbers BQKY01000001-BQKY01000019). KEY POINTS: • R. paludigena BS15 was isolated after an extensive screening of extracellular oil producers from natural sources. • Fed-batch fermentation of R. paludigena BS15 yielded 110 g/L of PEFA, which is the highest titer ever reported to date. • Combined analysis using Illumina and Oxford Nanopore sequencers produced the near-complete genome sequence.

Single-haplotype comparative genomics provides insights into lineage-specific structural variation during cat evolution

The role of structurally dynamic genomic regions in speciation is poorly understood due to challenges inherent in diploid genome assembly. Here we reconstructed the evolutionary dynamics of structural variation in five cat species by phasing the genomes of three interspecies F1 hybrids to generate near-gapless single-haplotype assemblies. We discerned that cat genomes have a paucity of segmental duplications relative to great apes, explaining their remarkable karyotypic stability. X chromosomes were hotspots of structural variation, including enrichment with inversions in a large recombination desert with characteristics of a supergene. The X-linked macrosatellite DXZ4 evolves more rapidly than 99.5% of the genome clarifying its role in felid hybrid incompatibility. Resolved sensory gene repertoires revealed functional copy number changes associated with ecomorphological adaptations, sociality and domestication. This study highlights the value of gapless genomes to reveal structural mechanisms underpinning karyotypic evolution, reproductive isolation and ecological niche adaptation.

Chromosome-scale assembly of the wild wheat relative Aegilops umbellulata

Wild wheat relatives have been explored in plant breeding to increase the genetic diversity of bread wheat, one of the most important food crops. Aegilops umbellulata is a diploid U genome-containing grass species that serves as a genetic reservoir for wheat improvement. In this study, we report the construction of a chromosome-scale reference assembly of Ae. umbellulata accession TA1851 based on corrected PacBio HiFi reads and chromosome conformation capture. The total assembly size was 4.25 Gb with a contig N50 of 17.7 Mb. In total, 36,268 gene models were predicted. We benchmarked the performance of hifiasm and LJA, two of the most widely used assemblers using standard and corrected HiFi reads, revealing a positive effect of corrected input reads. Comparative genome analysis confirmed substantial chromosome rearrangements in Ae. umbellulata compared to bread wheat. In summary, the Ae. umbellulata assembly provides a resource for comparative genomics in Triticeae and for the discovery of agriculturally important genes.

Telomere sequence variability in genotypes from natural plant populations: unusual block-organized double-monomer terminal telomeric arrays

BACKGROUND

Telomeres are the nucleoprotein complexes that physically cap the ends of eukaryotic chromosomes. Most plants possess Arabidopsis-type telomere sequences (TSs). In addition to terminal TSs, more diverse interstitial TSs exists in plants. Although telomeres have been sufficiently studied, the actual diversity of TSs in land plants is underestimated.

RESULTS

We investigate genotypes from seven natural populations with contrasting environments of four Chenopodium species to reveal the variability in TSs by analyzing Oxford Nanopore reads. Fluorescent in situ hybridization was used to localize telomeric repeats on chromosomes. We identified a number of derivative monomers that arise in part of both terminal and interstitial telomeric arrays of a single genotype. The former presents a case of block-organized double-monomer telomers, where blocks of Arabidopsis-type TTTAGGG motifs were interspersed with blocks of derivative TTTAAAA motifs. The latter is an integral part of the satellitome with transformations specific to the inactive genome fraction.

CONCLUSIONS

We suggested two alternative models for the possible formation of derivative monomers from telomeric heptamer motifs of Arabidopsis-type. It was assumed that derivatization of TSs is a ubiquitous process in the plant genome but occurrence and frequencies of derivatives may be genotype-specific. We also propose that the formation of non-canonical arrays of TSs, especially at chromosomal termini, may be a source for genomic variability in nature.

Multi-omic analyses reveal the unique properties of chia (Salvia hispanica) seed metabolism

Chia (Salvia hispanica) is an emerging crop considered a functional food containing important substances with multiple potential applications. However, the molecular basis of some relevant chia traits, such as seed mucilage and polyphenol content, remains to be discovered. This study generates an improved chromosome-level reference of the chia genome, resolving some highly repetitive regions, describing methylation patterns, and refining genome annotation. Transcriptomic analysis shows that seeds exhibit a unique expression pattern compared to other organs and tissues. Thus, a metabolic and proteomic approach is implemented to study seed composition and seed-produced mucilage. The chia genome exhibits a significant expansion in mucilage synthesis genes (compared to Arabidopsis), and gene network analysis reveals potential regulators controlling seed mucilage production. Rosmarinic acid, a compound with enormous therapeutic potential, was classified as the most abundant polyphenol in seeds, and candidate genes for its complex pathway are described. Overall, this study provides important insights into the molecular basis for the unique characteristics of chia seeds.

Allium cytogenetics: a critical review on the Indian taxa

The genus Allium Linnaeus, 1753 (tribe Allieae) contains about 800 species worldwide of which almost 38 species are reported in India, including the globally important crops (onion, garlic, leek, shallot) and many wild species. A satisfactory chromosomal catalogue of Allium species is missing which has been considered in the review for the species occurring in India. The most prominent base number is x=8, with few records of x=7, 10, 11. The genome size has sufficient clues for divergence, ranging from 7.8 pg/1C to 30.0 pg/1C in diploid and 15.16 pg/1C to 41.78 pg/1C in polyploid species. Although the karyotypes are seemingly dominated by metacentrics, substantial variation in nucleolus organizing regions (NORs) is noteworthy. The chromosomal rearrangement between A.cepa Linnaeus, 1753 and its allied species has paved way to appreciate genomic evolution within Allium. The presence of a unique telomere sequence and its conservation in Allium sets this genus apart from all other Amaryllids and supports monophyletic origin. Any cytogenetic investigation regarding NOR variability, telomere sequence and genome size in the Indian species becomes the most promising field to decipher chromosome evolution against the background of species diversity and evolution, especially in the Indian subcontinent.

The telomere-to-telomere genome of Fragaria vesca reveals the genomic evolution of Fragaria and the origin of cultivated octoploid strawberry

Fragaria vesca, commonly known as wild or woodland strawberry, is the most widely distributed diploid Fragaria species and is native to Europe and Asia. Because of its small plant size, low heterozygosity, and relative ease of genetic transformation, F. vesca has been a model plant for fruit research since the publication of its Illumina-based genome in 2011. However, its genomic contribution to octoploid cultivated strawberry remains a long-standing question. Here, we de novo assembled and annotated a telomere-to-telomere, gap-free genome of F. vesca 'Hawaii 4', with all seven chromosomes assembled into single contigs, providing the highest completeness and assembly quality to date. The gap-free genome is 220 785 082 bp in length and encodes 36 173 protein-coding gene models, including 1153 newly annotated genes. All 14 telomeres and seven centromeres were annotated within the seven chromosomes. Among the three previously recognized wild diploid strawberry ancestors, F. vesca, F. iinumae, and F. viridis, phylogenomic analysis showed that F. vesca and F. viridis are the ancestors of the cultivated octoploid strawberry F. × ananassa, and F. vesca is its closest relative. Three subgenomes of F. × ananassa belong to the F. vesca group, and one is sister to F. viridis. We anticipate that this high-quality, telomere-to-telomere, gap-free F. vesca genome, combined with our phylogenomic inference of the origin of cultivated strawberry, will provide insight into the genomic evolution of Fragaria and facilitate strawberry genetics and molecular breeding.

Large-Scale Detection of Telomeric Motif Sequences in Genomic Data Using TelFinder

Telomeres are regions of tandem repeated sequences at the ends of linear chromosomes that protect against DNA damage and chromosome fusion. Telomeres are associated with senescence and cancers and have attracted the attention of an increasing number of researchers. However, few telomeric motif sequences are known. Given the mounting interest in telomeres, an efficient computational tool for the de novo detection of the telomeric motif sequence of new species is needed since experimental-based methods are costly in terms of time and effort. Here, we report the development of TelFinder, an easy-to-use and freely available tool for the de novo detection of telomeric motif sequences from genomic data. The vast quantity of readily available genomic data makes it possible to apply this tool to any species of interest, which will undoubtedly inspire studies requiring telomeric repeat information and improve the utilization of these genomic data sets. We have tested TelFinder on telomeric sequences available in the Telomerase Database, and the detection accuracy reaches 90%. In addition, variation analyses in telomere sequences can be performed by TelFinder for the first time. The telomere variation preference of different chromosomes and even at the ends of the chromosome can provide clues regarding the underlying mechanisms of telomeres. Overall, these results shed new light on the divergent evolution of telomeres. IMPORTANCE Telomeres are reported to be highly correlated with the cell cycle and aging. As a result, research on telomere composition and evolution has become more and more urgent. However, using experimental methods to detect telomeric motif sequences is slow and costly. To combat this challenge, we developed TelFinder, a computational tool for the de novo detection of the telomere composition only using genomic data. In this study, we showed that a lot of complicated telomeric motifs could be identified by TelFinder only using genomic data. In addition, TelFinder can be used to check variation analyses in telomere sequences, which could lead to a deeper understanding of telomere sequences.

A chromosome-level genome assembly of Plantago ovata

Plantago ovata is cultivated for production of its seed husk (psyllium). When wet, the husk transforms into a mucilage with properties suitable for pharmaceutical industries, utilised in supplements for controlling blood cholesterol levels, and food industries for making gluten-free products. There has been limited success in improving husk quantity and quality through breeding approaches, partly due to the lack of a reference genome. Here we constructed the first chromosome-scale reference assembly of P. ovata using a combination of 5.98 million PacBio and 636.5 million Hi-C reads. We also used corrected PacBio reads to estimate genome size and transcripts to generate gene models. The final assembly covers ~ 500 Mb with 99.3% gene set completeness. A total of 97% of the sequences are anchored to four chromosomes with an N50 of ~ 128.87 Mb. The P. ovata genome contains 61.90% repeats, where 40.04% are long terminal repeats. We identified 41,820 protein-coding genes, 411 non-coding RNAs, 108 ribosomal RNAs, and 1295 transfer RNAs. This genome will provide a resource for plant breeding programs to, for example, reduce agronomic constraints such as seed shattering, increase psyllium yield and quality, and overcome crop disease susceptibility.

Release 4.0 of the Plant rDNA Database: A Database on Plant Ribosomal DNA Loci Number, Their Position, and Organization: An Information Source for Comparative Cytogenetics

This paper presents the latest update to the Plant rDNA database (Release 4.0), a valuable resource for researchers in the field of plant cytogenetics. The database provides information on the number, position, and arrangement of ribosomal DNA loci in plants, including angiosperms, gymnosperms, bryophytes, and pteridophytes. The new release includes new data for 820 species coming from additional 173 papers. In the updated version of the Plant rDNA database, 4948 entries comprising 2760 organisms can be found. A brief guide on how to navigate the database and obtain the desired information is also provided. The regular updating of the database is important for ensuring the information it contains is accurate, up-to-date, and useful for the research community. The Plant rDNA database continues to be beneficial for phylogenetic and cytogenetic studies in a wide range of taxa including angiosperms, gymnosperms, and early diverging groups, such as bryophytes and lycophytes.

Identification of the Sequence and the Length of Telomere DNA

Telomeres are essential nucleoprotein structures at the very ends of linear eukaryote chromosomes. They shelter the terminal genome territories against degradation and prevent the natural chromosome ends from being recognized by repair mechanisms as double-strand DNA breaks.There are two basic characteristics of telomeric DNA, its sequence and its length. The telomere sequence is important as a "landing area" for specific telomere-binding proteins, which function as signals and moderate the interactions required for correct telomere function. While the sequence forms the proper "landing surface" of telomeric DNA, its length is similarly important. Too short or exceptionally long telomere DNA cannot perform its function properly. In this chapter, methods for the investigation of these two basic telomere DNA characteristics are described, namely, telomere motif identification and telomere length measurement.

Chromosome-Scale Genome Assembly and Annotation of Allotetraploid Annual Bluegrass (Poa annua L.)

Poa annua L. is a globally distributed grass with economic and horticultural significance as a weed and as a turfgrass. This dual significance, and its phenotypic plasticity and ecological adaptation, have made P. annua an intriguing plant for genetic and evolutionary studies. Because of the lack of genomic resources and its allotetraploid (2n = 4x = 28) nature, a reference genome sequence would be a valuable asset to better understand the significance and polyploid origin of P. annua. Here we report a genome assembly with scaffolds representing the 14 haploid chromosomes that are 1.78 Gb in length with an N50 of 112 Mb and 96.7% of BUSCO orthologs. Seventy percent of the genome was identified as repetitive elements, 91.0% of which were Copia- or Gypsy-like long-terminal repeats. The genome was annotated with 76,420 genes spanning 13.3% of the 14 chromosomes. The two subgenomes originating from Poa infirma (Knuth) and Poa supina (Schrad) were sufficiently divergent to be distinguishable but syntenic in sequence and annotation with repetitive elements contributing to the expansion of the P. infirma subgenome.

Telomerase RNA in Hymenoptera (Insecta) switched to plant/ciliate-like biogenesis

In contrast to the catalytic subunit of telomerase, its RNA subunit (TR) is highly divergent in size, sequence and biogenesis pathways across eukaryotes. Current views on TR evolution assume a common origin of TRs transcribed with RNA polymerase II in Opisthokonta (the supergroup including Animalia and Fungi) and Trypanosomida on one hand, and TRs transcribed with RNA polymerase III under the control of type 3 promoter, found in TSAR and Archaeplastida supergroups (including e.g. ciliates and Viridiplantae taxa, respectively). Here, we focus on unknown TRs in one of the largest Animalia order - Hymenoptera (Arthropoda) with more than 300 available representative genomes. Using a combination of bioinformatic and experimental approaches, we identify their TRs. In contrast to the presumed type of TRs (H/ACA box snoRNAs transcribed with RNA Polymerase II) corresponding to their phylogenetic position, we find here short TRs of the snRNA type, likely transcribed with RNA polymerase III under the control of the type 3 promoter. The newly described insect TRs thus question the hitherto assumed monophyletic origin of TRs across Animalia and point to an evolutionary switch in TR type and biogenesis that was associated with the divergence of Arthropods.

Mega-sized pericentromeric blocks of simple telomeric repeats and their variants reveal patterns of chromosome evolution in ancient Cycadales genomes

Simple telomeric repeats composed of six to seven iterating nucleotide units are important sequences typically found at the ends of chromosomes. Here we analyzed their abundance and homogeneity in 42 gymnosperm (29 newly sequenced), 29 angiosperm (one newly sequenced), and eight bryophytes using bioinformatics, conventional cytogenetic and molecular biology approaches to explore their diversity across land plants. We found more than 10 000-fold variation in the amounts of telomeric repeats among the investigated taxa. Repeat abundance was positively correlated with increasing intragenomic sequence heterogeneity and occurrence at non-telomeric positions, but there was no correlation with genome size. The highest abundance/heterogeneity was found in the gymnosperm genus Cycas (Cycadaceae), in which megabase-sized blocks of telomeric repeats (i.e., billions of copies) were identified. Fluorescent in situ hybridization experiments using variant-specific probes revealed canonical Arabidopsis-type telomeric TTTAGGG repeats at chromosome ends, while pericentromeric blocks comprised at least four major telomeric variants with decreasing abundance: TTTAGGG>TTCAGGG >TTTAAGG>TTCAAGG. Such a diversity of repeats was not found in the sister cycad family Zamiaceae or in any other species analyzed. Using immunocytochemistry, we showed that the pericentromeric blocks of telomeric repeats overlapped with histone H3 serine 10 phosphorylation signals. We show that species of Cycas have amplified their telomeric repeats in centromeric and telomeric positions on telocentric chromosomes to extraordinary high levels. The ancestral chromosome number reconstruction suggests their occurrence is unlikely to be the product of ancient Robertsonian chromosome fusions. We speculate as to how the observed chromosome dynamics may be associated with the diversification of cycads.

Telomeres and Their Neighbors

Telomeres are essential structures formed from satellite DNA repeats at the ends of chromosomes in most eukaryotes. Satellite DNA repeat sequences are useful markers for karyotyping, but have a more enigmatic role in the eukaryotic cell. Much work has been done to investigate the structure and arrangement of repetitive DNA elements in classical models with implications for species evolution. Still more is needed until there is a complete picture of the biological function of DNA satellite sequences, particularly when considering non-model organisms. Celebrating Gregor Mendel’s anniversary by going to the roots, this review is designed to inspire and aid new research into telomeres and satellites with a particular focus on non-model organisms and accessible experimental and in silico methods that do not require specialized equipment or expensive materials. We describe how to identify telomere (and satellite) repeats giving many examples of published (and some unpublished) data from these techniques to illustrate the principles behind the experiments. We also present advice on how to perform and analyse such experiments, including details of common pitfalls. Our examples are a selection of recent developments and underexplored areas of research from the past. As a nod to Mendel’s early work, we use many examples from plants and insects, especially as much recent work has expanded beyond the human and yeast models traditional in telomere research. We give a general introduction to the accepted knowledge of telomere and satellite systems and include references to specialized reviews for the interested reader.

Tieing together loose ends: telomere instability in cancer and aging

Telomere maintenance is essential for maintaining genome integrity in both normal and cancer cells. Without functional telomeres, chromosomes lose their protective structure and undergo fusion and breakage events that drive further genome instability, including cell arrest or death. One means by which this loss can be overcome in stem cells and cancer cells is via re‐addition of G‐rich telomeric repeats by the telomerase reverse transcriptase (TERT). During aging of somatic tissues, however, insufficient telomerase expression leads to a proliferative arrest called replicative senescence, which is triggered when telomeres reach a critically short threshold that induces a DNA damage response. Cancer cells express telomerase but do not entirely escape telomere instability as they often possess short telomeres; hence there is often selection for genetic alterations in the TERT promoter that result in increased telomerase expression. In this review, we discuss our current understanding of the consequences of telomere instability in cancer and aging, and outline the opportunities and challenges that lie ahead in exploiting the reliance of cells on telomere maintenance for preserving genome stability.

FISH Mapping of Telomeric and Non-Telomeric (AG3T3)3 Reveal the Chromosome Numbers and Chromosome Rearrangements of 41 Woody Plants

Data for the chromosomal FISH mapping localization of (AG₃T₃)₃ are compiled for 37 species belonging 27 families; for 24 species and 14 families, this is the first such report. The chromosome number and length ranged from 14–136 and 0.56–14.48 μm, respectively. A total of 23 woody plants presented chromosome length less than 3 μm, thus belonging to the small chromosome group. Telomeric signals were observed at each chromosome terminus in 38 plants (90.5%) and were absent at several chromosome termini in only four woody plants (9.5%). Non-telomeric signals were observed in the chromosomes of 23 plants (54.8%); in particular, abundant non-telomeric (AG₃T₃)₃ was obviously observed in Chimonanthus campanulatus. Telomeric signals outside of the chromosome were observed in 11 woody plants (26.2%). Overall, ten (AG₃T₃)₃ signal pattern types were determined, indicating the complex genome architecture of the 37 considered species. The variation in signal pattern was likely due to chromosome deletion, duplication, inversion, and translocation. In addition, large primary constriction was observed in some species, probably due to or leading to chromosome breakage and the formation of new chromosomes. The presented results will guide further research focused on determining the chromosome number and disclosing chromosome rearrangements of woody plants.

Plant telomere biology: The green solution to the end-replication problem

Telomere maintenance is a fundamental cellular process conserved across all eukaryotic lineages. Although plants and animals diverged over 1.5 billion years ago, lessons learned from plants continue to push the boundaries of science, revealing detailed molecular mechanisms in telomere biology with broad implications for human health, aging biology, and stress responses. Recent studies of plant telomeres have unveiled unexpected divergence in telomere sequence and architecture, and the proteins that engage telomeric DNA and telomerase. The discovery of telomerase RNA components in the plant kingdom and some algae groups revealed new insight into the divergent evolution and the universal core of telomerase across major eukaryotic kingdoms. In addition, resources cataloging the abundant natural variation in Arabidopsis thaliana, maize (Zea mays), and other plants are providing unparalleled opportunities to understand the genetic networks that govern telomere length polymorphism and, as a result, are uncovering unanticipated crosstalk between telomeres, environmental factors, organismal fitness, and plant physiology. Here we recap current advances in plant telomere biology and put this field in perspective relative to telomere and telomerase research in other eukaryotic lineages.

Five Fabaceae Karyotype and Phylogenetic Relationship Analysis Based on Oligo-FISH for 5S rDNA and (AG3T3)3

Most Fabaceae have nitrogen fixation abilities and are valuable forage and medicinal resources. However, cytogenetic data of many Fabaceae species are unclear. Karyotypes reveal cytological characteristics and are crucial to understanding the organization and evolution of chromosomes in species. Oligo-FISH can reveal genetic composition and karyotype variation patterns with rapid and efficient results. Karyotype analysis of five Fabaceae species by oligonucleotide probes showed that: Robinia pseudoacacia, karyotype formula 2n = 2x = 20m + 2sm, cytotype 2B, arm ratio 3.4821, eight chromosomes distributed 5S rDNA signal. The karyotype formula of Robinia pseudoacacia ‘idaho’ was 2n = 2x = 20m + 2sm, cytotype 1A, arm ratio 1.8997, and 5S rDNA signal was distributed on six chromosomes. Karyotype of Robinia pseudoacacia f. decaisneana 2n = 2x = 20m + 2sm, cytotype 1B, arm ratio 2.0787, the distribution of eight chromosomes with 5S rDNA signal. Karyotype formula of Styphnolobium japonicum 2n = 2x = 14m + 12sm + 2st, cytotype 2B, arm ratio 2.6847, two chromosomes have 5S rDNA signal. Amorpha fruticose karyotype 2n = 2x = 38m + 2sm, cytotype 1B, arm ratio 3.2058, four chromosomes possessed 5S rDNA signal. Both ends of all species’ chromosomes have (AG₃T₃)₃ signals. The results of this study provide chromosome numbers and a physical map, contributing to the construction of the Oligo-FISH barcode and providing molecular cytogenetics data for Fabaceae.

Profiles of telomeric repeats in Insecta reveal diverse forms of telomeric motifs in Hymenopterans

Telomeres consist of highly conserved simple tandem telomeric repeat motif (TRM): (TTAGG)n in arthropods, (TTAGGG)n in vertebrates, and (TTTAGGG)n in most plants. TRM can be detected from chromosome-level assembly, which typically requires long-read sequencing data. To take advantage of short-read data, we developed an ultra-fast Telomeric Repeats Identification Pipeline and evaluated its performance on 91 species. With proven accuracy, we applied Telomeric Repeats Identification Pipeline in 129 insect species, using 7 Tbp of short-read sequences. We confirmed (TTAGG)n as the TRM in 19 orders, suggesting it is the ancestral form in insects. Systematic profiling in Hymenopterans revealed a diverse range of TRMs, including the canonical 5-bp TTAGG (bees, ants, and basal sawflies), three independent losses of tandem repeat form TRM (Ichneumonoids, hunting wasps, and gall-forming wasps), and most interestingly, a common 8-bp (TTATTGGG)n in Chalcid wasps with two 9-bp variants in the miniature wasp (TTACTTGGG) and fig wasps (TTATTGGGG). Our results identified extraordinary evolutionary fluidity of Hymenopteran TRMs, and rapid evolution of TRM and repeat abundance at all evolutionary scales, providing novel insights into telomere evolution.

Born in the mitochondrion and raised in the nucleus: evolution of a novel tandem repeat family in Medicago polymorpha (Fabaceae)

Plant nuclear genomes harbor sequence elements derived from the organelles (mitochondrion and plastid) through intracellular gene transfer (IGT). Nuclear genomes also show a dramatic range of repeat content, suggesting that any sequence can be readily amplified. These two aspects of plant nuclear genomes are well recognized but have rarely been linked. Through investigation of 31 Medicago taxa we detected exceptionally high post-IGT amplification of mitochondrial (mt) DNA sequences containing rps10 in the nuclear genome of Medicago polymorpha and closely related species. The amplified sequences were characterized as tandem arrays of five distinct repeat motifs (2157, 1064, 987, 971, and 587 bp) that have diverged from the mt genome (mitogenome) in the M. polymorpha nuclear genome. The mt rps10-like arrays were identified in seven loci (six intergenic and one telomeric) of the nuclear chromosome assemblies and were the most abundant tandem repeat family, representing 1.6-3.0% of total genomic DNA, a value approximately three-fold greater than the entire mitogenome in M. polymorpha. Compared to a typical mt gene, the mt rps10-like sequence coverage level was 691.5-7198-fold higher in M. polymorpha and closely related species. In addition to the post-IGT amplification, our analysis identified the canonical telomeric repeat and the species-specific satellite arrays that are likely attributable to an ancestral chromosomal fusion in M. polymorpha. A possible relationship between chromosomal instability and the mt rps10-like tandem repeat family in the M. polymorpha clade is discussed.

Complete Genome Sequencing and Comparative Analysis of the Clinically-Derived Apiotrichum mycotoxinivorans Strain GMU1709

Over the past decade, Apiotrichum mycotoxinivorans has been recognized globally as a source of opportunistic infections. It is a yeast-like fungus, and its association as an uncommon pulmonary pathogen with cystic fibrosis patients has been previously reported. Immunocompromised patients are at the highest risk of A. mycotoxinivorans infections. Therefore, to investigate the genetic basis for the pathogenicity of A. mycotoxinivorans, we performed whole-genome sequencing and comparative genomic analysis of A. mycotoxinivorans GMU1709 that was isolated from sputum specimens of a pneumonia patient receiving cardiac repair surgery. The assembly of Oxford Nanopore reads from the GMU1709 strain and its subsequent correction using Illumina paired-end reads yielded a high-quality complete genome with a genome size of 30.5 Mb in length, which comprised six chromosomes and one mitochondrion. Subsequently, 8,066 protein-coding genes were predicted based on multiple pieces of evidence, including transcriptomes. Phylogenomic analysis indicated that A. mycotoxinivorans exhibited the closest evolutionary affinity to A. veenhuisii, and both the A. mycotoxinivorans strains and the formerly Trichosporon cutaneum ACCC 20271 strain occupied the same phylogenetic position. Further comparative analysis supported that the ACCC 20271 strain belonged to A. mycotoxinivorans. Comparisons of three A. mycotoxinivorans strains indicated that the differences between clinical and non-clinical strains in pathogenicity and drug resistance may be little or none. Based on the comparisons with strains of other species in the Trichosporonaceae family, we identified potential key genetic factors associated with A. mycotoxinivorans infection or pathogenicity. In addition, we also deduced that A. mycotoxinivorans had great potential to inactivate some antibiotics (e.g., tetracycline), which may affect the efficacy of these drugs in co-infection. In general, our analyses provide a better understanding of the classification and phylogeny of the Trichosporonaceae family, uncover the underlying genetic basis of A. mycotoxinivorans infections and associated drug resistance, and provide clues into potential targets for further research and the therapeutic intervention of infections.

Oligo-FISH Can Identify Chromosomes and Distinguish Hippophaë rhamnoides L. Taxa

Oligo-fluorescence in situ hybridization (FISH) facilitates precise chromosome identification and comparative cytogenetic analysis. Detection of autosomal chromosomes of Hippophaë rhamnoides has not been achieved using oligonucleotide sequences. Here, the chromosomes of five H. rhamnoides taxa in the mitotic metaphase and mitotic metaphase to anaphase were detected using the oligo-FISH probes (AG3T3)3, 5S rDNA, and (TTG)6. In total, 24 small chromosomes were clearly observed in the mitotic metaphase (0.89–3.03 μm), whereas 24–48 small chromosomes were observed in the mitotic metaphase to anaphase (0.94–3.10 μm). The signal number and intensity of (AG3T3)3, 5S rDNA, and (TTG)6 in the mitotic metaphase to anaphase chromosomes were nearly consistent with those in the mitotic metaphase chromosomes when the two split chromosomes were integrated as one unit. Of note, 14 chromosomes (there is a high chance that sex chromosomes are included) were exclusively identified by (AG3T3)3, 5S rDNA, and (TTG)6. The other 10 also showed a terminal signal with (AG3T3)3. Moreover, these oligo-probes were able to distinguish one wild H. rhamnoides taxon from four H. rhamnoides taxa. These chromosome identification and taxa differentiation data will help in elucidating visual and elaborate physical mapping and guide breeders’ utilization of wild resources of H. rhamnoides.

Interstitial Arabidopsis-Type Telomeric Repeats in Asteraceae

Tandem repeats of telomeric-like motifs at intra-chromosomal regions, known as interstitial telomeric repeats (ITR), have drawn attention as potential markers of structural changes, which might convey information about evolutionary relationships if preserved through time. Building on our previous work that reported outstanding ITR polymorphisms in the genus Anacyclus, we undertook a survey across 132 Asteraceae species, focusing on the six most speciose subfamilies and considering all the ITR data published to date. The goal was to assess whether the presence, site number, and chromosomal location of ITRs convey any phylogenetic signal. We conducted fluorescent in situ hybridization (FISH) using an Arabidopsis-type telomeric sequence as a probe on karyotypes obtained from mitotic chromosomes. FISH signals of ITR sites were detected in species of subfamilies Asteroideae, Carduoideae, Cichorioideae, Gymnarhenoideae, and Mutisioideae, but not in Barnadesioideae. Although six small subfamilies have not yet been sampled, altogether, our results suggest that the dynamics of ITR formation in Asteraceae cannot accurately trace the complex karyological evolution that occurred since the early diversification of this family. Thus, ITRs do not convey a reliable signal at deep or shallow phylogenetic levels and cannot help to delimitate taxonomic categories, a conclusion that might also hold for other important families such as Fabaceae.

Antibacterial kinetics and phylogenetic analysis of Aloe vera plants

Uncontrolled use of antibiotics has resulted in the emergence of resistant bacteria. It has necessitated the evaluation of antibacterial activities and phylo-diversity of Aloe vera (also called Aloe barbadensis) plants as antimicrobial agent in Nigeria. Biotyped enteric bacilli of 251 strains obtained from fecal samples of patients with various gastro-intestinal complications are profiled for antibiogram. Resistant biotypes were assayed for susceptibility to Aloe vera latex and further evaluated for time-kill kinetics and phylo-diversity. More than 30% of enteric bacilli, including Citrobacter freundii, Escherichia coli and Proteus mirabilis were resistant to cotrimoxazole, ciprofloxacin, and tetracycline respectively at MIC >16 µg/ml (p=0.004). Aloe vera latex significantly inhibited 39.5% resistant enteric biotypes with a significant average reduction of the viable count at 1xMIC and 2xMIC to less than 3.0 Log10CFU/mL after 24 hours. Flavonoids, alkaloids, terpenoids and anthraquinine in anti-enteric sap significantly correlated and regressed with antibacterial activity (p<0.05), while two of the antimicrobial Aloe vera plants showed phylogenetic relatedness with other homologous. Anti-bacteria efficacy of some Nigerian Aloe vera latex could provide alternative therapy, while its phylo-diversity and genomic profiling would offer a promising avenue for identification and development of antimicrobial agents as drug candidates for natural antibiotics.

Genome assembly of two nematode-resistant cotton lines (Gossypium hirsutum L.)

Upland cotton (Gossypium hirsutum L.) is susceptible to damage by the root-knot and the reniform nematodes, causing yield losses greater than 4% annually in the United States. In addition, these nematodes are synergistic with seeding disease and root rot pathogens that exacerbate diseases and subsequent yield losses. Production practices to minimize nematode damage include crop rotation and nematicides, but these techniques need to be repeated and are expensive. The use of resistant cultivars is deemed the most effective and economical approach for managing nematodes in cotton. Here, we describe the genomes of two nematode-resistant lines of cotton, BARBREN-713 and BAR 32-30. These genomes may expedite the development of DNA markers that can be used to efficiently introduce nematode resistance into commercially valuable Upland lines.

Telomere Roles in Fungal Genome Evolution and Adaptation

Telomeres form the ends of linear chromosomes and usually comprise protein complexes that bind to simple repeated sequence motifs that are added to the 3′ ends of DNA by the telomerase reverse transcriptase (TERT). One of the primary functions attributed to telomeres is to solve the “end-replication problem” which, if left unaddressed, would cause gradual, inexorable attrition of sequences from the chromosome ends and, eventually, loss of viability. Telomere-binding proteins also protect the chromosome from 5′ to 3′ exonuclease action, and disguise the chromosome ends from the double-strand break repair machinery whose illegitimate action potentially generates catastrophic chromosome aberrations. Telomeres are of special interest in the blast fungus, Pyricularia, because the adjacent regions are enriched in genes controlling interactions with host plants, and the chromosome ends show enhanced polymorphism and genetic instability. Previously, we showed that telomere instability in some P. oryzae strains is caused by novel retrotransposons (MoTeRs) that insert in telomere repeats, generating interstitial telomere sequences that drive frequent, break-induced rearrangements. Here, we sought to gain further insight on telomeric involvement in shaping Pyricularia genome architecture by characterizing sequence polymorphisms at chromosome ends, and surrounding internalized MoTeR loci (relics) and interstitial telomere repeats. This provided evidence that telomere dynamics have played historical, and likely ongoing, roles in shaping the Pyricularia genome. We further demonstrate that even telomeres lacking MoTeR insertions are poorly preserved, such that the telomere-adjacent sequences exhibit frequent presence/absence polymorphism, as well as exchanges with the genome interior. Using TERT knockout experiments, we characterized chromosomal responses to failed telomere maintenance which suggested that much of the MoTeR relic-/interstitial telomere-associated polymorphism could be driven by compromised telomere function. Finally, we describe three possible examples of a phenomenon known as “Adaptive Telomere Failure,” where spontaneous losses of telomere maintenance drive rapid accumulation of sequence polymorphism with possible adaptive advantages. Together, our data suggest that telomere maintenance is frequently compromised in Pyricularia but the chromosome alterations resulting from telomere failure are not as catastrophic as prior research would predict, and may, in fact, be potent drivers of adaptive polymorphism.

Evolution of plant telomerase RNAs: farther to the past, deeper to the roots

The enormous sequence heterogeneity of telomerase RNA (TR) subunits has thus far complicated their characterization in a wider phylogenetic range. Our recent finding that land plant TRs are, similarly to known ciliate TRs, transcribed by RNA polymerase III and under the control of the type-3 promoter, allowed us to design a novel strategy to characterize TRs in early diverging Viridiplantae taxa, as well as in ciliates and other Diaphoretickes lineages. Starting with the characterization of the upstream sequence element of the type 3 promoter that is conserved in a number of small nuclear RNAs, and the expected minimum TR template region as search features, we identified candidate TRs in selected Diaphoretickes genomes. Homologous TRs were then used to build covariance models to identify TRs in more distant species. Transcripts of the identified TRs were confirmed by transcriptomic data, RT-PCR and Northern hybridization. A templating role for one of our candidates was validated in Physcomitrium patens. Analysis of secondary structure demonstrated a deep conservation of motifs (pseudoknot and template boundary element) observed in all published TRs. These results elucidate the evolution of the earliest eukaryotic TRs, linking the common origin of TRs across Diaphoretickes, and underlying evolutionary transitions in telomere repeats.

Telomeric DNA sequences in beetle taxa vary with species richness

Telomeres are protective structures at the ends of eukaryotic chromosomes, and disruption of their nucleoprotein composition usually results in genome instability and cell death. Telomeric DNA sequences have generally been found to be exceptionally conserved in evolution, and the most common pattern of telomeric sequences across eukaryotes is (TxAyGz)n maintained by telomerase. However, telomerase-added DNA repeats in some insect taxa frequently vary, show unusual features, and can even be absent. It has been speculated about factors that might allow frequent changes in telomere composition in Insecta. Coleoptera (beetles) is the largest of all insect orders and based on previously available data, it seemed that the telomeric sequence of beetles varies to a great extent. We performed an extensive mapping of the (TTAGG)n sequence, the ancestral telomeric sequence in Insects, across the main branches of Coleoptera. Our study indicates that the (TTAGG)n sequence has been repeatedly or completely lost in more than half of the tested beetle superfamilies. Although the exact telomeric motif in most of the (TTAGG)n-negative beetles is unknown, we found that the (TTAGG)n sequence has been replaced by two alternative telomeric motifs, the (TCAGG)n and (TTAGGG)n, in at least three superfamilies of Coleoptera. The diversity of the telomeric motifs was positively related to the species richness of taxa, regardless of the age of the taxa. The presence/absence of the (TTAGG)n sequence highly varied within the Curculionoidea, Chrysomeloidea, and Staphylinoidea, which are the three most diverse superfamilies within Metazoa. Our data supports the hypothesis that telomere dysfunctions can initiate rapid genomic changes that lead to reproductive isolation and speciation.

Extraordinary diversity of telomeres, telomerase RNAs and their template regions in Saccharomycetaceae

Telomerase RNA (TR) carries the template for synthesis of telomere DNA and provides a scaffold for telomerase assembly. Fungal TRs are long and have been compared to higher eukaryotes, where they show considerable diversity within phylogenetically close groups. TRs of several Saccharomycetaceae were recently identified, however, many of these remained uncharacterised in the template region. Here we show that this is mainly due to high variability in telomere sequence. We predicted the telomere sequences using Tandem Repeats Finder and then we identified corresponding putative template regions in TR candidates. Remarkably long telomere units and the corresponding putative TRs were found in Tetrapisispora species. Notably, variable lengths of the annealing sequence of the template region (1–10 nt) were found. Consequently, species with the same telomere sequence may not harbour identical TR templates. Thus, TR sequence alone can be used to predict a template region and telomere sequence, but not to determine these exactly. A conserved feature of telomere sequences, tracts of adjacent Gs, led us to test the propensity of individual telomere sequences to form G4. The results show highly diverse values of G4-propensity, indicating the lack of ubiquitous conservation of this feature across Saccharomycetaceae.

Whole-genome sequencing and genome regions of special interest: Lessons from major histocompatibility complex, sex determination, and plant self-incompatibility

Whole‐genome sequencing of non‐model organisms is now widely accessible and has allowed a range of questions in the field of molecular ecology to be investigated with greater power. However, some genomic regions that are of high biological interest remain problematic for assembly and data‐handling. Three such regions are the major histocompatibility complex (MHC), sex‐determining regions (SDRs) and the plant self‐incompatibility locus (S‐locus). Using these as examples, we illustrate the challenges of both assembling and resequencing these highly polymorphic regions and how bioinformatic and technological developments are enabling new approaches to their study. Mapping short‐read sequences against multiple alternative references improves genotyping comprehensiveness at the S‐locus thereby contributing to more accurate assessments of allelic frequencies. Long‐read sequencing, producing reads of several tens to hundreds of kilobase pairs in length, facilitates the assembly of such regions as single sequences can span the multiple duplicated gene copies of the MHC region, and sequence through repetitive stretches and translocations in SDRs and S‐locus haplotypes. These advances are adding value to short‐read genome resequencing approaches by allowing, for example, more accurate haplotype phasing across longer regions. Finally, we assessed further technical improvements, such as nanopore adaptive sequencing and bioinformatic tools using pangenomes, which have the potential to further expand our knowledge of a number of genomic regions that remain challenging to study with classical resequencing approaches.

Comparative genomics of Chlamydomonas

Despite its role as a reference organism in the plant sciences, the green alga Chlamydomonas reinhardtii entirely lacks genomic resources from closely related species. We present highly contiguous and well-annotated genome assemblies for three unicellular C. reinhardtii relatives: Chlamydomonas incerta, Chlamydomonas schloesseri, and the more distantly related Edaphochlamys debaryana. The three Chlamydomonas genomes are highly syntenous with similar gene contents, although the 129.2 Mb C. incerta and 130.2 Mb C. schloesseri assemblies are more repeat-rich than the 111.1 Mb C. reinhardtii genome. We identify the major centromeric repeat in C. reinhardtii as a LINE transposable element homologous to Zepp (the centromeric repeat in Coccomyxa subellipsoidea) and infer that centromere locations and structure are likely conserved in C. incerta and C. schloesseri. We report extensive rearrangements, but limited gene turnover, between the minus mating type loci of these Chlamydomonas species. We produce an eight-species core-Reinhardtinia whole-genome alignment, which we use to identify several hundred false positive and missing genes in the C. reinhardtii annotation and >260,000 evolutionarily conserved elements in the C. reinhardtii genome. In summary, these resources will enable comparative genomics analyses for C. reinhardtii, significantly extending the analytical toolkit for this emerging model system.

Charting the genomic landscape of seed-free plants

During the past few years several high-quality genomes has been published from Charophyte algae, bryophytes, lycophytes and ferns. These genomes have not only elucidated the origin and evolution of early land plants, but have also provided important insights into the biology of the seed-free lineages. However, critical gaps across the phylogeny remain and many new questions have been raised through comparing seed-free and seed plant genomes. Here, we review the reference genomes available and identify those that are missing in the seed-free lineages. We compare patterns of various levels of genome and epigenomic organization found in seed-free plants to those of seed plants. Some genomic features appear to be fundamentally different. For instance, hornworts, Selaginella and most liverworts are devoid of whole-genome duplication, in stark contrast to other land plants. In addition, the distribution of genes and repeats appear to be less structured in seed-free genomes than in other plants, and the levels of gene body methylation appear to be much lower. Finally, we highlight the currently available (or needed) model systems, which are crucial to further our understanding about how changes in genes translate into evolutionary novelties.

Yarrowia lipolytica: a multitalented yeast species of ecological significance

Yarrowia lipolytica is characterized by GRAS (Generally regarded as safe) status, the versatile substrate utilization profile, rapid utilization rates, metabolic diversity and flexibility, the unique abilities to tolerate to extreme environments (acidic, alkaline, hypersaline, heavy metal-pollutions and others) and elevated biosynthesis and secreting capacities. These advantages of Y. lipolytica allow us to consider it as having great ecological significance. Unfortunately, there is still a paucity of relevant review data. This mini-review highlights ecological ubiquity of Y. lipolytica species, their ability to diversify and colonize specialized niches. Different Y. lipolytica strains, native and engineered, are beneficial in degrading many environmental pollutants causing serious ecological problems worldwide. In agriculture has a potential to be a bio-control agent by stimulating plant defense response, and an eco-friendly bio-fertilizer. Engineered strains of Y. lipolytica have become a very promising platform for eco-friendly production of biofuel, commodities, chemicals and secondary metabolites of plant origin, obtaining which by other method were limited or economically infeasible, or were accompanied by stringent environmental problems. Perspectives to use potential of Y. lipolytica's capacities for industrial scale production of valuable compounds in an eco-friendly manner are proposed.

Origin and Fates of TERT Gene Copies in Polyploid Plants

The gene coding for the telomerase reverse transcriptase (TERT) is essential for the maintenance of telomeres. Previously we described the presence of three TERT paralogs in the allotetraploid plant Nicotiana tabacum, while a single TERT copy was identified in the paleopolyploid model plant Arabidopsis thaliana. Here we examine the presence, origin and functional status of TERT variants in allotetraploid Nicotiana species of diverse evolutionary ages and their parental genome donors, as well as in other diploid and polyploid plant species. A combination of experimental and in silico bottom-up analyses of TERT gene copies in Nicotiana polyploids revealed various patterns of retention or loss of parental TERT variants and divergence in their functions. RT-qPCR results confirmed the expression of all the identified TERT variants. In representative plant and green algal genomes, our synteny analyses show that their TERT genes were located in a conserved locus that became advantageous after the divergence of eudicots, and the gene was later translocated in several plant groups. In various diploid and polyploid species, translocation of TERT became fixed in target loci that show ancient synapomorphy.

Telomeres and Subtelomeres Dynamics in the Context of Early Chromosome Interactions During Meiosis and Their Implications in Plant Breeding

Genomic architecture facilitates chromosome recognition, pairing, and recombination. Telomeres and subtelomeres play an important role at the beginning of meiosis in specific chromosome recognition and pairing, which are critical processes that allow chromosome recombination between homologs (equivalent chromosomes in the same genome) in later stages. In plant polyploids, these terminal regions are even more important in terms of homologous chromosome recognition, due to the presence of homoeologs (equivalent chromosomes from related genomes). Although telomeres interaction seems to assist homologous pairing and consequently, the progression of meiosis, other chromosome regions, such as subtelomeres, need to be considered, because the DNA sequence of telomeres is not chromosome-specific. In addition, recombination operates at subtelomeres and, as it happens in rye and wheat, homologous recognition and pairing is more often correlated with recombining regions than with crossover-poor regions. In a plant breeding context, the knowledge of how homologous chromosomes initiate pairing at the beginning of meiosis can contribute to chromosome manipulation in hybrids or interspecific genetic crosses. Thus, recombination in interspecific chromosome associations could be promoted with the aim of transferring desirable agronomic traits from related genetic donor species into crops. In this review, we summarize the importance of telomeres and subtelomeres on chromatin dynamics during early meiosis stages and their implications in recombination in a plant breeding framework.

Step-by-Step Evolution of Telomeres: Lessons from Yeasts

In virtually every eukaryotic species, the ends of nuclear chromosomes are protected by telomeres, nucleoprotein structures counteracting the end-replication problem and suppressing recombination and undue DNA repair. Although in most cases, the primary structure of telomeric DNA is conserved, there are several exceptions to this rule. One is represented by the telomeric repeats of ascomycetous yeasts, which encompass a great variety of sequences, whose evolutionary origin has been puzzling for several decades. At present, the key questions concerning the driving force behind their rapid evolution and the means of co-evolution of telomeric repeats and telomere-binding proteins remain largely unanswered. Previously published studies addressed mostly the general concepts of the evolutionary origin of telomeres, key properties of telomeric proteins as well as the molecular mechanisms of telomere maintenance; however, the evolutionary process itself has not been analyzed thoroughly. Here, we aimed to inspect the evolution of telomeres in ascomycetous yeasts from the subphyla Saccharomycotina and Taphrinomycotina, with special focus on the evolutionary origin of species-specific telomeric repeats. We analyzed the sequences of telomeric repeats from 204 yeast species classified into 20 families and as a result, we propose a step-by-step model, which integrates the diversity of telomeric repeats, telomerase RNAs, telomere-binding protein complexes and explains a propensity of certain species to generate the repeat heterogeneity within a single telomeric array.

A Reference Genome Sequence for Giant Sequoia

The giant sequoia (Sequoiadendron giganteum) of California are massive, long-lived trees that grow along the U.S. Sierra Nevada mountains. Genomic data are limited in giant sequoia and producing a reference genome sequence has been an important goal to allow marker development for restoration and management. Using deep-coverage Illumina and Oxford Nanopore sequencing, combined with Dovetail chromosome conformation capture libraries, the genome was assembled into eleven chromosome-scale scaffolds containing 8.125 Gbp of sequence. Iso-Seq transcripts, assembled from three distinct tissues, was used as evidence to annotate a total of 41,632 protein-coding genes. The genome was found to contain, distributed unevenly across all 11 chromosomes and in 63 orthogroups, over 900 complete or partial predicted NLR genes, of which 375 are supported by annotation derived from protein evidence and gene modeling. This giant sequoia reference genome sequence represents the first genome sequenced in the Cupressaceae family, and lays a foundation for using genomic tools to aid in giant sequoia conservation and management.

Composition and Function of Telomerase-A Polymerase Associated with the Origin of Eukaryotes

The canonical DNA polymerases involved in the replication of the genome are unable to fully replicate the physical ends of linear chromosomes, called telomeres. Chromosomal termini thus become shortened in each cell cycle. The maintenance of telomeres requires telomerase—a specific RNA-dependent DNA polymerase enzyme complex that carries its own RNA template and adds telomeric repeats to the ends of chromosomes using a reverse transcription mechanism. Both core subunits of telomerase—its catalytic telomerase reverse transcriptase (TERT) subunit and telomerase RNA (TR) component—were identified in quick succession in Tetrahymena more than 30 years ago. Since then, both telomerase subunits have been described in various organisms including yeasts, mammals, birds, reptiles and fish. Despite the fact that telomerase activity in plants was described 25 years ago and the TERT subunit four years later, a genuine plant TR has only recently been identified by our group. In this review, we focus on the structure, composition and function of telomerases. In addition, we discuss the origin and phylogenetic divergence of this unique RNA-dependent DNA polymerase as a witness of early eukaryotic evolution. Specifically, we discuss the latest information regarding the recently discovered TR component in plants, its conservation and its structural features.

Human-like telomeres in Zostera marina reveal a mode of transition from the plant to the human telomeric sequences

A previous study describing the genome of Zostera marina, the most widespread seagrass in the Northern hemisphere, revealed some genomic signatures of adaptation to the aquatic environment such as the loss of stomatal genes, while other functions such as an algal-like cell wall composition were acquired. Beyond these, the genome structure and organization were comparable with those of the majority of plant genomes sequenced, except for one striking feature that went unnoticed at that time: the presence of human-like instead of the expected plant-type telomeric sequences. By using different experimental approaches including fluorescence in situ hybridization (FISH), genome skimming by next-generation sequencing (NGS), and analysis of non-coding transcriptome, we have confirmed its telomeric location in the chromosomes of Z. marina. We have also identified its telomerase RNA (TR) subunit, confirming the presence of the human-type telomeric sequence in the template region. Remarkably, this region was found to be very variable even in clades with a highly conserved telomeric sequence across their species. Based on this observation, we propose that alternative annealing preferences in the template borders can explain the transition between the plant and human telomeric sequences. The further identification of paralogues of TR in several plant genomes led us to the hypothesis that plants may retain an increased ability to change their telomeric sequence. We discuss the implications of this occurrence in the evolution of telomeres while introducing a mechanistic model for the transition from the plant to the human telomeric sequences.

Telomerase as a Possible Candidate Targeting Therapy in Different Breast Cancer Cell Lines

Background:

Telomerase activity is up regulated in most breast cancer subtypes but not in the adjacent normal tissues. Thus, it is a promising target for anticancer therapy. The present work investigated the effects of telomerase inhibition by siRNA on breast cancer cell lines and studied the feasibility of whether the combined effect of doxorubicin with siRNA treatment on breast cancer cells potentiates a rapid cellular response to the cytotoxic effect of chemotherapy.

Methods:

This study was performed on Luminal A (MCF-7), triple negative (MDA-MB-468), and HER-2/neu (SKBR-3) human breast cancer cell lines, wherein telomerase activity inhibition by hTERT siRNA and doxorubicin was detected by measuring telomerase activity using Telomeric Repeat Amplification Protocol (TRAP assay), assessing cell viability through MTT assay, and evaluating apoptosis through scanning electron microscopy (SEM) and through estimating caspase-3 and -8 activities using enzyme-linked immunosorbent assay (ELISA).

Results:

In the present study, hTERT siRNA effectively reduced telomerase activity and cell viability to more than 90% and 60%, respectively, in most breast cancer cell lines within 72 hours after transfection. The combination of hTERT siRNA and doxorubicin showed a cumulative effect compared with either treatment alone (P < 0.05). Meanwhile, SEM demonstrated apoptotic morphologic cell changes.

Conclusion:

Telomerase inhibition is a promising strategy for the effective treatment of breast cancer. When used in combination with doxorubicin, it could potentiate the cytotoxic effect of the drug on breast cancer cells.