Correlated evolution of LTR retrotransposons and genome size in the genus Eleocharis

The pan-genome of Spodoptera frugiperda provides new insights into genome evolution and horizontal gene transfer

Spodoptera frugiperda is a common and severely damaging agricultural pest. In-depth analysis of its population genomics and transcriptomics is crucial for providing references for pest control efforts. This study, focused on the extensive variation in the genome size of S. frugiperda, constructed its pan-genome and identified 1.37 Gb of non-reference sequences, highlighting significant genetic variation within the population. Analysis of Long Terminal Repeat (LTR) Presence/Absence Variation (PAV) suggests that LTR alterations may be one of the driving factors for genome size variation. Additionally, population gene PAV analysis revealed that variable genes are enriched in functions like acetyltransferase activity, which might be associated with detoxification, implying diverse selection pressures related to detoxification in different S. frugiperda populations. Moreover, 19 horizontal gene transfer (HGT) acquired genes were identified in the reference genome used in this study, which responded to 16 different treatments. Notably, three HGT-acquired genes (SFR02618, SFR05248, and SFR05249) co-expressed with heat shock protein family and responded under treatments with Avermectin and Cypermethrin. This may indicate their involvement in a detoxification mechanism coordinated with heat shock proteins. These results offering new insights into its genomic evolution and the potential functions of HGT-acquired genes.

The genome of Eleocharis vivipara elucidates the genetics of C3-C4 photosynthetic plasticity and karyotype evolution in the Cyperaceae

Eleocharis vivipara, an amphibious sedge in the Cyperaceae family, has several remarkable properties, most notably its alternate use of C3 photosynthesis underwater and C4 photosynthesis on land. However, the absence of genomic data has hindered its utility for evolutionary and genetic research. Here, we present a high-quality genome for E. vivipara, representing the first chromosome-level genome for the Eleocharis genus, with an approximate size of 965.22 Mb mainly distributed across 10 chromosomes. Its Hi-C pattern, chromosome clustering results, and one-to-one genome synteny across two subgroups indicates a tetraploid structure with chromosome count 2n = 4x = 20. Phylogenetic analysis suggests that E. vivipara diverged from Cyperus esculentus approximately 32.96 million years ago (Mya), and underwent a whole-genome duplication (WGD) about 3.5 Mya. Numerous fusion and fission events were identified between the chromosomes of E. vivipara and its close relatives. We demonstrate that E. vivipara has holocentromeres, a chromosomal feature which can maintain the stability of such chromosomal rearrangements. Experimental transplantation and cross-section studies showed its terrestrial culms developed C4 Kranz anatomy with increased number of chloroplasts in the bundle sheath (BS) cells. Gene expression and weighted gene co-expression network analysis (WGCNA) showed overall elevated expression of core genes associated with the C4 pathway, and significant enrichment of genes related to modified culm anatomy and photosynthesis efficiency. We found evidence of mixed nicotinamide adenine dinucleotide - malic enzyme and phosphoenolpyruvate carboxykinase type C4 photosynthesis in E. vivipara, and hypothesize that the evolution of C4 photosynthesis predates the WGD event. The mixed type is dominated by subgenome A and supplemented by subgenome B. Collectively, our findings not only shed light on the evolution of E. vivipara and karyotype within the Cyperaceae family, but also provide valuable insights into the transition between C3 and C4 photosynthesis, offering promising avenues for crop improvement and breeding.

Look4LTRs: a Long terminal repeat retrotransposon detection tool capable of cross species studies and discovering recently nested repeats

Plant genomes include large numbers of transposable elements. One particular type of these elements is flanked by two Long Terminal Repeats (LTRs) and can translocate using RNA. Such elements are known as LTR-retrotransposons; they are the most abundant type of transposons in plant genomes. They have many important functions involving gene regulation and the rise of new genes and pseudo genes in response to severe stress. Additionally, LTR-retrotransposons have several applications in biotechnology. Due to the abundance and the importance of LTR-retrotransposons, multiple computational tools have been developed for their detection. However, none of these tools take advantages of the availability of related genomes; they process one chromosome at a time. Further, recently nested LTR-retrotransposons (multiple elements of the same family are inserted into each other) cannot be annotated accurately - or cannot be annotated at all - by the currently available tools. Motivated to overcome these two limitations, we built Look4LTRs, which can annotate LTR-retrotransposons in multiple related genomes simultaneously and discover recently nested elements. The methodology of Look4LTRs depends on techniques imported from the signal-processing field, graph algorithms, and machine learning with a minimal use of alignment algorithms. Four plant genomes were used in developing Look4LTRs and eight plant genomes for evaluating it in contrast to three related tools. Look4LTRs is the fastest while maintaining better or comparable F1 scores (the harmonic average of recall and precision) to those obtained by the other tools. Our results demonstrate the added benefit of annotating LTR-retrotransposons in multiple related genomes simultaneously and the ability to discover recently nested elements. Expert human manual examination of six elements - not included in the ground truth - revealed that three elements belong to known families and two elements are likely from new families. With respect to examining recently nested LTR-retrotransposons, three out of five were confirmed to be valid elements. Look4LTRs - with its speed, accuracy, and novel features - represents a true advancement in the annotation of LTR-retrotransposons, opening the door to many studies focused on understanding their functions in plants.

Repeat turnover meets stable chromosomes: repetitive DNA sequences mark speciation and gene pool boundaries in sugar beet and wild beets

Sugar beet and its wild relatives share a base chromosome number of nine and similar chromosome morphologies. Yet, interspecific breeding is impeded by chromosome and sequence divergence that is still not fully understood. Since repetitive DNAs are among the fastest evolving parts of the genome, we investigated, if repeatome innovations and losses are linked to chromosomal differentiation and speciation. We traced genome and chromosome-wide evolution across 13 beet species comprising all sections of the genera Beta and Patellifolia. For this, we combined short and long read sequencing, flow cytometry, and cytogenetics to build a comprehensive framework that spans the complete scale from DNA to chromosome to genome. Genome sizes and repeat profiles reflect the separation into three gene pools with contrasting evolutionary patterns. Among all repeats, satellite DNAs harbor most genomic variability, leading to fundamentally different centromere architectures, ranging from chromosomal uniformity in Beta and Patellifolia to the formation of patchwork chromosomes in Corollinae/Nanae. We show that repetitive DNAs are causal for the genome expansions and contractions across the beet genera, providing insights into the genomic underpinnings of beet speciation. Satellite DNAs in particular vary considerably between beet genomes, leading to the evolution of distinct chromosomal setups in the three gene pools, likely contributing to the barriers in beet breeding. Thus, with their isokaryotypic chromosome sets, beet genomes present an ideal system for studying the link between repeats, genomic variability, and chromosomal differentiation and provide a theoretical fundament for understanding barriers in any crop breeding effort.

The first homosporous lycophyte genome revealed the association between the recent dynamic accumulation of LTR-RTs and genome size variation

The contrasting genome size between homosporous and heterosporous plants is fascinating. Different from the heterosporous seed plants and mainly homosporous ferns, the lycophytes are either heterosporous (Isoetales and Selaginellales) or homosporous (Lycopodiales). Many lycophytes are the resource plants of Huperzine A (HupA) which is invaluable for treating Alzheimer's disease. For the seed-free vascular plants, several high-quality genomes of heterosporous Selaginella, homosporous ferns (maidenhair fern, monkey spider tree fern), and heterosporous ferns (Azolla) have been published and provided important insights into the origin and evolution of early land plants. However, the homosporous lycophyte genome has not been decoded. Here, we assembled the first homosporous lycophyte genome and conducted comparative genomic analyses by applying a reformed pipeline for filtering out non-plant sequences. The obtained genome size of Lycopodium clavatum is 2.30 Gb, distinguished in more than 85% repetitive elements of which 62% is long terminal repeat (LTR). This study disclosed a high birth rate and a low death rate of the LTR-RTs in homosporous lycophytes, but the opposite occurs in heterosporous lycophytes. we propose that the recent activity of LTR-RT is responsible for the immense genome size variation between homosporous and heterosporous lycophytes. By combing Ks analysis with a phylogenetic approach, we discovered two whole genome duplications (WGD). Morover, we identified all the five recognized key enzymes for the HupA biosynthetic pathway in the L. clavatum genome, but found this pathway incomplete in other major lineages of land plants. Overall, this study is of great importance for the medicinal utilization of lycophytes and the decoded genome data will be a key cornerstone to elucidate the evolution and biology of early vascular land plants.

Complex patterns of ploidy in a holocentric plant clade (Schoenus, Cyperaceae) in the Cape biodiversity hotspot

BACKGROUND AND AIMS

It is unclear how widespread polyploidy is throughout the largest holocentric plant family - the Cyperaceae. Because of the prevalence of chromosomal fusions and fissions, which affect chromosome number but not genome size, it can be impossible to distinguish if individual plants are polyploids in holocentric lineages based on chromosome count data alone. Furthermore, it is unclear how differences in genome size and ploidy levels relate to environmental correlates within holocentric lineages, such as the Cyperaceae.

METHODS

We focus our analyses on tribe Schoeneae, and more specifically the southern African clade of Schoenus. We examine broad-scale patterns of genome size evolution in tribe Schoeneae and focus more intensely on determining the prevalence of polyploidy across the southern African Schoenus by inferring ploidy level with the program ChromEvol, as well as interpreting chromosome number and genome size data. We further investigate whether there are relationships between genome size/ploidy level and environmental variables across the nutrient-poor and summer-arid Cape biodiversity hotspot.

KEY RESULTS

Our results show a large increase in genome size, but not chromosome number, within Schoenus compared to other species in tribe Schoeneae. Across Schoenus, there is a positive relationship between chromosome number and genome size, and our results suggest that polyploidy is a relatively common process throughout the southern African Schoenus. At the regional scale of the Cape, we show that polyploids are more often associated with drier locations that have more variation in precipitation between dry and wet months, but these results are sensitive to the classification of ploidy level.

CONCLUSIONS

Polyploidy is relatively common in the southern African Schoenus, where a positive relationship is observed between chromosome number and genome size. Thus, there may be a high incidence of polyploidy in holocentric plants, whose cell division properties differ from monocentrics.

Phylogenomics and genome size evolution in Amomum s. s. (Zingiberaceae): Comparison of traditional and modern sequencing methods

BACKGROUND AND AIMS

A targeted enrichment NGS approach was used to construct the phylogeny of Amomum Roxb. (Zingiberaceae). Phylogenies based on hundreds of nuclear genes, the whole plastome and the rDNA cistron were compared with an ITS-based phylogeny. Trends in genome size (GS) evolution were examined, chromosomes were counted and the geographical distribution of phylogenetic lineages was evaluated.

METHODS

In total, 92 accessions of 54 species were analysed. ITS was obtained for 79 accessions, 37 accessions were processed with Hyb-Seq and sequences from 449 nuclear genes, the whole cpDNA, and the rDNA cistron were analysed using concatenation, coalescence and supertree approaches. The evolution of absolute GS was analysed in a phylogenetic and geographical context. The chromosome numbers of 12 accessions were counted.

KEY RESULTS

Four groups were recognised in all datasets though their mutual relationships differ among datasets. While group A (A. subulatum and A. petaloideum) is basal to the remaining groups in the nuclear gene phylogeny, in the cpDNA topology it is sister to group B (A. repoeense and related species) and, in the ITS topology, it is sister to group D (the Elettariopsis lineage). The former Elettariopsis makes a monophyletic group. There is an increasing trend in GS during evolution. The largest GS values were found in group D in two tetraploid taxa, A. cinnamomeum and A. aff. biphyllum (both 2n = 96 chromosomes). The rest varied in GS (2C = 3.54-8.78 pg) with a constant chromosome number 2n = 48. There is a weak connection between phylogeny, GS and geography in Amomum.

CONCLUSIONS

Amomum consists of four groups, and the former Elettariopsis is monophyletic. Species in this group have the largest GS. Two polyploids were found and GS greatly varied in the rest of Amomum.

Chromosome size matters: genome evolution in the cyperid clade

BACKGROUND AND AIMS

While variation in genome size and chromosome numbers and their consequences are often investigated in plants, the biological relevance of variation in chromosome size remains poorly known. Here, we examine genome and mean chromosome size in the cyperid clade (families Cyperaceae, Juncaceae and Thurniaceae), which is the largest vascular plant lineage with predominantly holocentric chromosomes.

METHODS

We measured genome size in 436 species of cyperids using flow cytometry, and augment these data with previously published datasets. We then separately compared genome and mean chromosome sizes (2C/2n) amongst the major lineages of cyperids and analysed how these two genomic traits are associated with various environmental factors using phylogenetically informed methods.

KEY RESULTS

We show that cyperids have the smallest mean chromosome sizes recorded in seed plants, with a large divergence between the smallest and largest values. We found that cyperid species with smaller chromosomes have larger geographical distributions and that there is a strong inverse association between mean chromosome size and number across this lineage.

CONCLUSIONS

The distinct patterns in genome size and mean chromosome size across the cyperids might be explained by holokinetic drive. The numerous small chromosomes might function to increase genetic diversity in this lineage where crossovers are limited during meiosis.

Transposons and non-coding regions drive the intrafamily differences of genome size in insects

Genome size (GS) can vary considerably between phylogenetically close species, but the landscape of GS changes in insects remain largely unclear. To better understand the specific evolutionary factors that determine GS in insects, we examined flow cytometry-based published GS data from 1,326 insect species, spanning 700 genera, 155 families, and 21 orders. Model fitting showed that GS generally followed an Ornstein-Uhlenbeck adaptive evolutionary model in Insecta overall. Ancestral reconstruction indicated a likely GS of 1,069 Mb, suggesting that most insect clades appeared to undergo massive genome expansions or contractions. Quantification of genomic components in 56 species from nine families in four insect orders revealed that the proliferation of transposable elements contributed to high variation in GS between close species, such as within Coleoptera. This study sheds lights on the pattern of GS variation in insects and provides a better understanding of insect GS evolution.

Chromosome-level genome assembly of the aquatic plant Nymphoides indica reveals transposable element bursts and NBS-LRR gene family expansion shedding light on its invasiveness

Nymphoides indica, an aquatic plant, is an invasive species that causes both ecological and economic damage in North America and elsewhere. However, the lack of genomic data of N. indica limits the in-depth analysis of this invasive species. Here, we report a chromosome-level genome assembly of nine pseudochromosomes of N. indica with a total size of ∼ 520 Mb. More than half of the N. indica genome consists of transposable elements (TEs), and a higher density of TEs around genes may play a significant role in response to an ever-changing environment by regulating the nearby gene. Additionally, our analysis revealed that N. indica only experienced a gamma (γ) whole-genome triplication event. Functional enrichment of the N. indica-specific and expanded gene families highlighted genes involved in the responses to hypoxia and plant-pathogen interactions, which may strengthen the ability to adapt to external challenges and improve ecological fitness. Furthermore, we identified 160 members of the nucleotide-binding site and leucine-rich repeat gene family, which may be linked to the defence response. Collectively, the high-quality N. indica genome reported here opens a novel avenue to understand the evolution and rapid invasion of Nymphoides spp.

Hybrid Assembly and Annotation of the Genome of the Indian Punica granatum, a Superfood

The wonder fruit pomegranate (Punica granatum, family Lythraceae) is one of India’s economically important fruit crops that can grow in different agro-climatic conditions ranging from tropical to temperate regions. This study reports high-quality de novo draft hybrid genome assembly of diploid Punica cultivar “Bhagwa” and identifies its genomic features. This cultivar is most common among the farmers due to its high sustainability, glossy red color, soft seed, and nutraceutical properties with high market value. The draft genome assembly is about 361.76 Mb (N50 = 40 Mb), ∼9.0 Mb more than the genome size estimated by flow cytometry. The genome is 90.9% complete, and only 26.68% of the genome is occupied by transposable elements and has a relative abundance of 369.93 SSRs/Mb of the genome. A total of 30,803 proteins and their putative functions were predicted. Comparative whole-genome analysis revealed Eucalyptus grandis as the nearest neighbor. KEGG-KASS annotations indicated an abundance of genes involved in the biosynthesis of flavonoids, phenylpropanoids, and secondary metabolites, which are responsible for various medicinal properties of pomegranate, including anticancer, antihyperglycemic, antioxidant, and anti-inflammatory activities. The genome and gene annotations provide new insights into the pharmacological properties of the secondary metabolites synthesized in pomegranate. They will also serve as a valuable resource in mining biosynthetic pathways for key metabolites, novel genes, and variations associated with disease resistance, which can facilitate the breeding of new varieties with high yield and superior quality.

Horizontally Acquired nrDNAs Persist in Low Amounts in Host Hordeum Genomes and Evolve Independently of Native nrDNA

Nuclear ribosomal DNA (nrDNA) has displayed extraordinary dynamics during the evolution of plant species. However, the patterns and evolutionary significance of nrDNA array expansion or contraction are still relatively unknown. Moreover, only little is known of the fate of minority nrDNA copies acquired between species via horizontal transfer. The barley genus Hordeum (Poaceae) represents a good model for such a study, as species of section Stenostachys acquired nrDNA via horizontal transfer from at least five different panicoid genera, causing long-term co-existence of native (Hordeum-like) and non-native (panicoid) nrDNAs. Using quantitative PCR, we investigated copy number variation (CNV) of nrDNA in the diploid representatives of the genus Hordeum. We estimated the copy number of the foreign, as well as of the native ITS types (ribotypes), and followed the pattern of their CNV in relation to the genus' phylogeny, species' genomes size and the number of nrDNA loci. For the native ribotype, we encountered an almost 19-fold variation in the mean copy number among the taxa analysed, ranging from 1689 copies (per 2C content) in H. patagonicum subsp. mustersii to 31342 copies in H. murinum subsp. glaucum. The copy numbers did not correlate with any of the genus' phylogeny, the species' genome size or the number of nrDNA loci. The CNV was high within the recognised groups (up to 13.2 × in the American I-genome species) as well as between accessions of the same species (up to 4×). Foreign ribotypes represent only a small fraction of the total number of nrDNA copies. Their copy numbers ranged from single units to tens and rarely hundreds of copies. They amounted, on average, to between 0.1% (Setaria ribotype) and 1.9% (Euclasta ribotype) of total nrDNA. None of the foreign ribotypes showed significant differences with respect to phylogenetic groups recognised within the sect. Stenostachys. Overall, no correlation was found between copy numbers of native and foreign nrDNAs suggesting the sequestration and independent evolution of native and non-native nrDNA arrays. Therefore, foreign nrDNA in Hordeum likely poses a dead-end by-product of horizontal gene transfer events.

Holocentric Karyotype Evolution in Rhynchospora Is Marked by Intense Numerical, Structural, and Genome Size Changes

Cyperaceae is a family of Monocotyledons comprised of species with holocentric chromosomes that are associated with intense dysploidy and polyploidy events. Within this family the genus Rhynchospora has recently become the focus of several studies that characterize the organization of the holocentric karyotype and genome structures. To broaden our understanding of genome evolution in this genus, representatives of Rhynchospora were studied to contrast chromosome features, C-CMA/DAPI band distribution and genome sizes. Here, we carried out a comparative analysis for 35 taxa of Rhynchospora, and generated new genome size estimates for 20 taxa. The DNA 2C-values varied up to 22-fold, from 2C = 0.51 pg to 11.32 pg, and chromosome numbers ranged from 2n = 4 to 61. At least 37% of our sampling exhibited 2n different from the basic number x = 5, and chromosome rearrangements were also observed. A large variation in C-CMA/DAPI band accumulation and distribution was observed as well. We show that genome variation in Rhynchospora is much larger than previously reported. Phylogenetic analysis showed that most taxa were grouped in clades corresponding to previously described taxonomic sections. Basic chromosome numbers are the same within every section, however, changes appeared in all the clades. Ancestral chromosome number reconstruction revealed n = 5 as the most likely ancestral complements, but n = 10 appears as a new possibility. Chromosome evolution models point to polyploidy as the major driver of chromosome evolution in Rhynchospora, followed by dysploidy. A negative correlation between chromosome size and diploid number open the discussion for holokinetic drive-based genome evolution. This study explores relationships between karyotype differentiation and genome size variation in Rhynchospora, and contrasts it against the phylogeny of this holocentric group.

Comparative analysis of repetitive sequences among species from the potato and the tomato clades

BACKGROUND AND AIMS

The genus Solanum includes important vegetable crops and their wild relatives. Introgression of their useful traits into elite cultivars requires effective recombination between hom(e)ologues, which is partially determined by genome sequence differentiation. In this study we compared the repetitive genome fractions of wild and cultivated species of the potato and tomato clades in a phylogenetic context.

METHODS

Genome skimming followed by a clustering approach was used as implemented in the RepeatExplorer pipeline. Repeat classes were annotated and the sequences of their main domains were compared.

KEY RESULTS

Repeat abundance and genome size were correlated and the larger genomes of species in the tomato clade were found to contain a higher proportion of unclassified elements. Families and lineages of repetitive elements were largely conserved between the clades, but their relative proportions differed. The most abundant repeats were Ty3/Gypsy elements. Striking differences in abundance were found in the highly dynamic Ty3/Gypsy Chromoviruses and Ty1/Copia Tork elements. Within the potato clade, early branching Solanum cardiophyllum showed a divergent repeat profile. There were also contrasts between cultivated and wild potatoes, mostly due to satellite amplification in the cultivated species. Interspersed repeat profiles were very similar among potatoes. The repeat profile of Solanum etuberosum was more similar to that of the potato clade.

CONCLUSIONS

The repeat profiles in Solanum seem to be very similar despite genome differentiation at the level of collinearity. Removal of transposable elements by unequal recombination may have been responsible for structural rearrangements across the tomato clade. Sequence variability in the tomato clade is congruent with clade-specific amplification of repeats after its divergence from S. etuberosum and potatoes. The low differentiation among potato and its wild relatives at the level of interspersed repeats may explain the difficulty in discriminating their genomes by genomic in situ hybridization techniques.

Transcriptional activation of long terminal repeat retrotransposon sequences in the genome of pitaya under abiotic stress

Frequent somatic variations exist in pitaya (Hylocereus undatus) plants grown under abiotic stress conditions. Long terminal repeat (LTR) retrotransposons can be activated under stressful conditions and play key roles in plant genetic variation and evolution. However, whether LTR retrotransposons promotes pitaya somatic variations by regulating abiotic stress responses is still uncertain. In this study, transcriptionally active LTR retrotransposons were identified in pitaya after exposure to a number of stress factors, including in vitro culturing, osmotic changes, extreme temperatures and hormone treatments. In total, 26 LTR retrotransposon reverse transcriptase (RT) cDNA sequences were isolated and identified as belonging to 9 Ty1-copia and 4 Ty3-gypsy families. Several RT cDNA sequences had differing similarity levels with RTs from pitaya genomic DNA and other plant species, and were differentially expressed in pitaya under various stress conditions. LTR retrotransposons accounted for at least 13.07% of the pitaya genome. HuTy1P4 had a high copy number and low expression level in young stems of pitaya, and its expression level increased after exposure to hormones and abiotic stresses, including in vitro culturing, osmotic changes, cold and heat. HuTy1P4 may have been subjected to diverse transposon events in 13 pitaya plantlets successively subcultured for four cycles. Thus, the expression levels of these retrotransposons in pitaya were associated with stress responses and may be involved in the occurrence of the somaclonal variation in pitaya.

Analysis of retrotransposon abundance, diversity and distribution in holocentric Eleocharis (Cyperaceae) genomes

BACKGROUND AND AIMS

Long terminal repeat-retrotransposons (LTR-RTs) comprise a large portion of plant genomes, with massive repeat blocks distributed across the chromosomes. Eleocharis species have holocentric chromosomes, and show a positive correlation between chromosome numbers and the amount of nuclear DNA. To evaluate the role of LTR-RTs in karyotype diversity in members of Eleocharis (subgenus Eleocharis), the occurrence and location of different members of the Copia and Gypsy superfamilies were compared, covering interspecific variations in ploidy levels (considering chromosome numbers), DNA C-values and chromosomal arrangements.

METHODS

The DNA C-value was estimated by flow cytometry. Genomes of Eleocharis elegans and E. geniculata were partially sequenced using Illumina MiSeq assemblies, which were a source for searching for conserved proteins of LTR-RTs. POL domains were used for recognition, comparing families and for probe production, considering different families of Copia and Gypsy superfamilies. Probes were obtained by PCR and used in fluorescence in situ hybridization (FISH) against chromosomes of seven Eleocharis species.

KEY RESULTS

A positive correlation between ploidy levels and the amount of nuclear DNA was observed, but with significant variations between samples with the same ploidy levels, associated with repetitive DNA fractions. LTR-RTs were abundant in E. elegans and E. geniculata genomes, with a predominance of Copia Sirevirus and Gypsy Athila/Tat clades. FISH using LTR-RT probes exhibited scattered and clustered signals, but with differences in the chromosomal locations of Copia and Gypsy. The diversity in LTR-RT locations suggests that there is no typical chromosomal distribution pattern for retrotransposons in holocentric chromosomes, except the CRM family with signals distributed along chromatids.

CONCLUSIONS

These data indicate independent fates for each LTR-RT family, including accumulation between and within chromosomes and genomes. Differential activity and small changes in LTR-RTs suggest a secondary role in nuclear DNA variation, when compared with ploidy changes.

Molecular characterization of Fagaceae species using inter-primer binding site (iPBS) markers

Retrotransposons (RTNs) contribute for genome evolution, influencing its size and structure. We investigated the utility of the RTN-based markers inter-primer binding site (iPBS) for the molecular characterization of 25 Fagaceae species from genera Castanea, Fagus and Quercus. The assessment of genetic diversity, relationships and structure, as well as taxonomic classification of Fagaceae based on molecular data is important for definition of conservation, forestry management strategies and discrimination among natural hybrids and their parents since natural hybridization may increase with the climate changes. Here, iPBS primers designed by other authors were tested alone and combined. Some of them were discriminative, revealed polymorphism within and among taxa allowing the production of a total of 150 iPBS markers. In addition, several monomorphic iPBS markers were also amplified in each taxon. The UPGMA dendrogram based on the pooled iPBS data revealed 27% of genetic similarity among species. The individuals were clustered per genus and most of the oaks per infrageneric group corroborating the adopted taxonomy. Globally, the iPBS markers demonstrated suitability for DNA fingerprinting, determination of phylogenies and taxonomic discrimination in Fagaceae, and could constitute a useful and alternative tool for germplasm characterization, and for definition of conservation strategies and forestry management. Moreover, these markers would be useful for fingerprinting natural hybrids that share morphological similarities with their parents. Since iPBS markers could also enable insights about RTNs evolution, an eventual correlation among iPBS polymorphism, variability of RTN insertions and/or genome size in Fagaceae is discussed.

Are holocentrics doomed to change? Limited chromosome number variation in Rhynchospora Vahl (Cyperaceae)

Karyotype evolution in species with non-localised centromeres (holocentric chromosomes) is usually very dynamic and associated with recurrent fission and fusion (also termed agmatoploidy/symploidy) events. In Rhynchospora (Cyperaceae), one of the most species-rich sedge genera, all analysed species have holocentric chromosomes and their numbers range from 2n = 4 to 2n = 84. Agmatoploidy/symploidy and polyploidy were suggested as the main processes in the reshuffling of Rhynchospora karyotypes, although testing different scenarios of chromosome number evolution in a phylogenetic framework has not been attempted until now. Here, we used maximum likelihood and model-based analyses, in combination with genome size estimation and ribosomal DNA distribution, to understand chromosome evolution in Rhynchospora. Overall, chromosome number variation showed a significant phylogenetic signal and the majority of the lineages maintained a karyotype of 2n = 10 (~48% of the species), the most likely candidate for the ancestral number of the genus. Higher and lower chromosome numbers were restricted to specific clades, whilst polyploidy and/or fusion/fission events were present in specific branches. Variation in genome size and ribosomal DNA site number showed no correlation with ploidy level or chromosome number. Although different mechanisms of karyotype evolution (polyploidy, fusion and fission) seem to be acting in distinct lineages, the degree of chromosome variation and the main mechanisms involved are comparable to those found in some monocentric genera and lower than expected for a holocentric genus.

Cold-induced retrotransposition of fish LINEs

Classes of retrotransposons constitute a large portion of metazoan genome. There have been cases reported that genomic abundance of retrotransposons is correlated with the severity of low environmental temperatures. However, the molecular mechanisms underlying such correlation are unknown. We show here by cell transfection assays that retrotransposition (RTP) of a long interspersed nuclear element (LINE) from an Antarctic notothenioid fish Dissostichus mawsoni (dmL1) could be activated by low temperature exposure, causing increased dmL1 copies in the host cell genome. The cold-induced dmL1 propagation was demonstrated to be mediated by the mitogen-activated protein kinases (MAPK)/p38 signaling pathway, which is activated by accumulation of reactive oxygen species (ROS) in cold-stressed conditions. Surprisingly, dmL1 transfected cells showed an increase in the number of viable cells after prolonged cold exposures than non-transfected cells. Features of cold inducibility of dmL1 were recapitulated in LINEs of zebrafish origin both in cultured cell lines and tissues, suggesting existence of a common cold-induced LINE amplification in fishes. The findings reveal an important function of LINEs in temperature adaptation and provid insights into the MAPK/p38 stress responsive pathway that shapes LINE composition in fishes facing cold stresses.

Useful parasites: the evolutionary biology and biotechnology applications of transposable elements

Transposable elements usually comprise the most abundant nongenic fraction of eukaryotic genomes. Because of their capacity to selfreplicate and to induce a wide range of mutations, transposable elements have long been considered as 'parasitic' or 'selfish'. Today, we recognize that the findings about genomic changes affected by transposable elements have considerably altered our view of the ways in which genomes evolve and work. Numerous studies have provided evidences that mobile elements have the potential to act as agents of evolution by increasing, rearranging and diversifying the genetic repertoire of their hosts. With large-scale sequencing becoming increasingly available, more and more scientists come across transposable element sequences in their data. I will provide examples that transposable elements, although having signatures of 'selfish' DNA, play a significant biological role in the maintainance of genome integrity and providing novel regulatoty networks. These features, along with the transpositional and mutagenic capacity to produce a raw genetic diversity, make the genome mobile fraction, a key player in species adaptation and microevolution. The last but not least, transposable elements stand as informative DNA markers that may complement other conventional DNA markers. Altogether, transposable elements represent a promising, but still largely unexplored research niche and deserve to be included into the agenda of molecular ecologists, evolutionary geneticists, conservation biologists and plant breeders.

Evolution of genome size and genomic GC content in carnivorous holokinetics (Droseraceae)

BACKGROUND AND AIMS

Studies in the carnivorous family Lentibulariaceae in the last years resulted in the discovery of the smallest plant genomes and an unusual pattern of genomic GC content evolution. However, scarcity of genomic data in other carnivorous clades still prevents a generalization of the observed patterns. Here the aim was to fill this gap by mapping genome evolution in the second largest carnivorous family, Droseraceae, where this evolution may be affected by chromosomal holokinetism in Drosera METHODS: The genome size and genomic GC content of 71 Droseraceae species were measured by flow cytometry. A dated phylogeny was constructed, and the evolution of both genomic parameters and their relationship to species climatic niches were tested using phylogeny-based statistics.

KEY RESULTS

The 2C genome size of Droseraceae varied between 488 and 10 927 Mbp, and the GC content ranged between 37·1 and 44·7 %. The genome sizes and genomic GC content of carnivorous and holocentric species did not differ from those of their non-carnivorous and monocentric relatives. The genomic GC content positively correlated with genome size and annual temperature fluctuations. The genome size and chromosome numbers were inversely correlated in the Australian clade of Drosera CONCLUSIONS: Our results indicate that neither carnivory (nutrient scarcity) nor the holokinetism have a prominent effect on size and DNA base composition of Droseraceae genomes. However, the holokinetic drive seems to affect karyotype evolution in one of the major clades of Drosera Our survey confirmed that the evolution of GC content is tightly connected with the evolution of genome size and also with environmental conditions.

EARE-1, a Transcriptionally Active Ty1/Copia-Like Retrotransposon Has Colonized the Genome of Excoecaria agallocha through Horizontal Transfer

Long terminal repeat (LTR) retrotransposons constitute the majority of the content of angiosperm genomes, but their evolutionary dynamics remain poorly understood. Here, we report the isolation and characterization of a putative full-length (~9550 bp) Ty1/copia-like retrotransposon in Excoecaria agallocha and its evolution in Euphorbiaceae. The so-called EARE-1 is phylogenetically closely related to RIRE-1 from Oryza australiensis, and has proliferated recently (~7.19 Mya) in the E. agallocha genome. An RT-PCR analysis revealed substantial transcription of EARE-1 in all examined organs (leaves, staminate flowers, pistillate flowers, seeds, and roots) in unstressed E. agallocha plants and indications of elevated expression under stress. We conducted sequence analyses of 256 RT-RH fragments (~860 bp) of EARE-1 from 34 species representing four subfamilies of Euphorbiaceae that exist in China. EARE-1 copies from two Excoecaria species and Phyllanthus urinaria showed incongruent phylogeny with the host species and exhibited high sequence similarity to the host genes, suggesting a horizontal transfer from P. urinaria to the common ancestor of Excoecaria. However, SSAP analysis detected no new insertions of EARE-1 among full-sibling progeny plants of E. agallocha, despite considerable SSAP polymorphisms among half-siblings. EARE-1 is the first transcriptionally active Ty1/copia-like retrotransposon isolated from E. agallocha. Our results provide empirical evidence of the horizontal transfer of LTR retrotransposons in plants, and may suggest a significant role of post-transcriptional host control in the life cycles of transposable elements.

Agmatoploidy and symploidy: a critical review

Agmatoploidy is a type of chromosome rearrangement that involves the fragmentation of an entire chromosome complement, generating a diploid with double its original chromosome number. Agmatoploidy and other related karyotype changes, such as symploidy (the opposite change, promoted by chromosome fusion), partial agmatoploidy, polyagmatoploidy, etc., are restricted to species with holokinetic chromosomes and are assumed to play an important role in their karyotype evolution. However, a critical review of the literature shows that examples of chromosome number doubling by agmatoploidy are rare and not clearly demonstrated, while partial agmatoploidy and partial symploidy seem to be the same as ascending and descending disploidy, respectively. It is therefore proposed here that these terms should be avoided or even abolished.

Flow cytometry may allow microscope-independent detection of holocentric chromosomes in plants

Two chromosomal structures, known as monocentric and holocentric chromosomes, have evolved in eukaryotes. Acentric fragments of monocentric chromosomes are unequally distributed to daughter cells and/or lost, while holocentric fragments are inherited normally. In monocentric species, unequal distribution should generate chimeras of cells with different nuclear DNA content. We investigated whether such differences in monocentric species are detectable by flow cytometry (FCM) as (i) a decreased nuclear DNA content and (ii) an increased coefficient of variance (CV) of the G1 peak after gamma radiation-induced fragmentation. We compared 13 monocentric and 9 holocentric plant species. Unexpectedly, monocentrics and holocentrics did not differ with respect to parameters (i) and (ii) in their response to gamma irradiation. However, we found that the proportion of G2 nuclei was highly elevated in monocentrics after irradiation, while holocentrics were negligibly affected. Therefore, we hypothesize that DNA-damaging agents induce cell cycle arrest leading to endopolyploidy only in monocentric and not (or to much lesser extent) in holocentric plants. While current microscope-dependent methods for holocentrism detection are unreliable for small and numerous chromosomes, which are common in holocentrics, FCM can use somatic nuclei. Thus, FCM may be a rapid and reliable method of high-throughput screening for holocentric candidates across plant phylogeny.

The Small Nuclear Genomes of Selaginella Are Associated with a Low Rate of Genome Size Evolution

The haploid nuclear genome size (1C DNA) of vascular land plants varies over several orders of magnitude. Much of this observed diversity in genome size is due to the proliferation and deletion of transposable elements. To date, all vascular land plant lineages with extremely small nuclear genomes represent recently derived states, having ancestors with much larger genome sizes. The Selaginellaceae represent an ancient lineage with extremely small genomes. It is unclear how small nuclear genomes evolved in Selaginella We compared the rates of nuclear genome size evolution in Selaginella and major vascular plant clades in a comparative phylogenetic framework. For the analyses, we collected 29 new flow cytometry estimates of haploid genome size in Selaginella to augment publicly available data. Selaginella possess some of the smallest known haploid nuclear genome sizes, as well as the lowest rate of genome size evolution observed across all vascular land plants included in our analyses. Additionally, our analyses provide strong support for a history of haploid nuclear genome size stasis in Selaginella Our results indicate that Selaginella, similar to other early diverging lineages of vascular land plants, has relatively low rates of genome size evolution. Further, our analyses highlight that a rapid transition to a small genome size is only one route to an extremely small genome.

Holokinetic centromeres and efficient telomere healing enable rapid karyotype evolution

Species with holocentric chromosomes are often characterized by a rapid karyotype evolution. In contrast to species with monocentric chromosomes where acentric fragments are lost during cell division, breakage of holocentric chromosomes creates fragments with normal centromere activity. To decipher the mechanism that allows holocentric species an accelerated karyotype evolution via chromosome breakage, we analyzed the chromosome complements of irradiated Luzula elegans plants. The resulting chromosomal fragments and rearranged chromosomes revealed holocentromere-typical CENH3 and histone H2AThr120ph signals as well as the same mitotic mobility like unfragmented chromosomes. Newly synthesized telomeres at break points become detectable 3 weeks after irradiation. The presence of active telomerase suggests a telomerase-based mechanism of chromosome healing. A successful transmission of holocentric chromosome fragments across different generations was found for most offspring of irradiated plants. Hence, a combination of holokinetic centromere activity and the fast formation of new telomeres at break points enables holocentric species a rapid karyotype evolution involving chromosome fissions and rearrangements.

Franklin FC, Heckmann S. Atypical centromeres in plants-what they can tell us

The centromere, visible as the primary constriction of condensed metaphase chromosomes, is a defined chromosomal locus essential for genome stability. It mediates transient assembly of a multi-protein complex, the kinetochore, which enables interaction with spindle fibers and thus faithful segregation of the genetic information during nuclear divisions. Centromeric DNA varies in extent and sequence composition among organisms, but a common feature of almost all active eukaryotic centromeres is the presence of the centromeric histone H3 variant cenH3 (a.k.a. CENP-A). These typical centromere features apply to most studied species. However, a number of species display "atypical" centromeres, such as holocentromeres (centromere extension along almost the entire chromatid length) or neocentromeres (ectopic centromere activity). In this review, we provide an overview of different atypical centromere types found in plants including holocentromeres, de novo formed centromeres and terminal neocentromeres as well as di-, tri- and metapolycentromeres (more than one centromere per chromosomes). We discuss their specific and common features and compare them to centromere types found in other eukaryotic species. We also highlight new insights into centromere biology gained in plants with atypical centromeres such as distinct mechanisms to define a holocentromere, specific adaptations in species with holocentromeres during meiosis or various scenarios leading to neocentromere formation.

Transposable elements and genome size variations in plants

Although the number of protein-coding genes is not highly variable between plant taxa, the DNA content in their genomes is highly variable, by as much as 2,056-fold from a 1C amount of 0.0648 pg to 132.5 pg. The mean 1C-value in plants is 2.4 pg, and genome size expansion/contraction is lineage-specific in plant taxonomy. Transposable element fractions in plant genomes are also variable, as low as ~3% in small genomes and as high as ~85% in large genomes, indicating that genome size is a linear function of transposable element content. Of the 2 classes of transposable elements, the dynamics of class 1 long terminal repeat (LTR) retrotransposons is a major contributor to the 1C value differences among plants. The activity of LTR retrotransposons is under the control of epigenetic suppressing mechanisms. Also, genome-purging mechanisms have been adopted to counter-balance the genome size amplification. With a wealth of information on whole-genome sequences in plant genomes, it was revealed that several genome-purging mechanisms have been employed, depending on plant taxa. Two genera, Lilium and Fritillaria, are known to have large genomes in angiosperms. There were twice times of concerted genome size evolutions in the family Liliaceae during the divergence of the current genera in Liliaceae. In addition to the LTR retrotransposons, non-LTR retrotransposons and satellite DNAs contributed to the huge genomes in the two genera by possible failure of genome counter-balancing mechanisms.

Holokinetic drive: centromere drive in chromosomes without centromeres

Similar to how the model of centromere drive explains the size and complexity of centromeres in monocentrics (organisms with localized centromeres), our model of holokinetic drive is consistent with the divergent evolution of chromosomal size and number in holocentrics (organisms with nonlocalized centromeres) exhibiting holokinetic meiosis (holokinetics). Holokinetic drive is proposed to facilitate chromosomal fission and/or repetitive DNA removal (or any segmental deletion) when smaller homologous chromosomes are preferentially inherited or chromosomal fusion and/or repetitive DNA proliferation (or any segmental duplication) when larger homologs are preferred. The hypothesis of holokinetic drive is supported primarily by the negative correlation between chromosome number and genome size that is documented in holokinetic lineages. The supporting value of two older cross-experiments on holokinetic structural heterozygotes (the rush Luzula elegans and butterflies of the genus Antheraea) that indicate the presence of size-preferential homolog transmission via female meiosis for holokinetic drive is discussed, along with the further potential consequences of holokinetic drive in comparison with centromere drive.

Casposons: a new superfamily of self-synthesizing DNA transposons at the origin of prokaryotic CRISPR-Cas immunity

BACKGROUND

Diverse transposable elements are abundant in genomes of cellular organisms from all three domains of life. Although transposons are often regarded as junk DNA, a growing body of evidence indicates that they are behind some of the major evolutionary innovations. With the growth in the number and diversity of sequenced genomes, previously unnoticed mobile elements continue to be discovered.

RESULTS

We describe a new superfamily of archaeal and bacterial mobile elements which we denote casposons because they encode Cas1 endonuclease, a key enzyme of the CRISPR-Cas adaptive immunity systems of archaea and bacteria. The casposons share several features with self-synthesizing eukaryotic DNA transposons of the Polinton/Maverick class, including terminal inverted repeats and genes for B family DNA polymerases. However, unlike any other known mobile elements, the casposons are predicted to rely on Cas1 for integration and excision, via a mechanism similar to the integration of new spacers into CRISPR loci. We identify three distinct families of casposons that differ in their gene repertoires and evolutionary provenance of the DNA polymerases. Deep branching of the casposon-encoded endonuclease in the Cas1 phylogeny suggests that casposons played a pivotal role in the emergence of CRISPR-Cas immunity.

CONCLUSIONS

The casposons are a novel superfamily of mobile elements, the first family of putative self-synthesizing transposons discovered in prokaryotes. The likely contribution of capsosons to the evolution of CRISPR-Cas parallels the involvement of the RAG1 transposase in vertebrate immunoglobulin gene rearrangement, suggesting that recruitment of endonucleases from mobile elements as ready-made tools for genome manipulation is a general route of evolution of adaptive immunity.

The contributions of transposable elements to the structure, function, and evolution of plant genomes

Transposable elements (TEs) are the key players in generating genomic novelty by a combination of the chromosome rearrangements they cause and the genes that come under their regulatory sway. Genome size, gene content, gene order, centromere function, and numerous other aspects of nuclear biology are driven by TE activity. Although the origins and attitudes of TEs have the hallmarks of selfish DNA, there are numerous cases where TE components have been co-opted by the host to create new genes or modify gene regulation. In particular, epigenetic regulation has been transformed from a process to silence invading TEs and viruses into a key strategy for regulating plant genes. Most, perhaps all, of this epigenetic regulation is derived from TE insertions near genes or TE-encoded factors that act in trans. Enormous pools of genome data and new technologies for reverse genetics will lead to a powerful new era of TE analysis in plants.

Genome size expansion and the relationship between nuclear DNA content and spore size in the Asplenium monanthes fern complex (Aspleniaceae)

BACKGROUND

Homosporous ferns are distinctive amongst the land plant lineages for their high chromosome numbers and enigmatic genomes. Genome size measurements are an under exploited tool in homosporous ferns and show great potential to provide an overview of the mechanisms that define genome evolution in these ferns. The aim of this study is to investigate the evolution of genome size and the relationship between genome size and spore size within the apomictic Asplenium monanthes fern complex and related lineages.

RESULTS

Comparative analyses to test for a relationship between spore size and genome size show that they are not correlated. The data do however provide evidence for marked genome size variation between species in this group. These results indicate that Asplenium monanthes has undergone a two-fold expansion in genome size.

CONCLUSIONS

Our findings challenge the widely held assumption that spore size can be used to infer ploidy levels within apomictic fern complexes. We argue that the observed genome size variation is likely to have arisen via increases in both chromosome number due to polyploidy and chromosome size due to amplification of repetitive DNA (e.g. transposable elements, especially retrotransposons). However, to date the latter has not been considered to be an important process of genome evolution within homosporous ferns. We infer that genome evolution, at least in some homosporous fern lineages, is a more dynamic process than existing studies would suggest.

Comparative genomic paleontology across plant kingdom reveals the dynamics of TE-driven genome evolution

Long terminal repeat-retrotransposons (LTR-RTs) are the most abundant class of transposable elements (TEs) in plants. They strongly impact the structure, function, and evolution of their host genome, and, in particular, their role in genome size variation has been clearly established. However, the dynamics of the process through which LTR-RTs have differentially shaped plant genomes is still poorly understood because of a lack of comparative studies. Using a new robust and automated family classification procedure, we exhaustively characterized the LTR-RTs in eight plant genomes for which a high-quality sequence is available (i.e., Arabidopsis thaliana, A. lyrata, grapevine, soybean, rice, Brachypodium dystachion, sorghum, and maize). This allowed us to perform a comparative genome-wide study of the retrotranspositional landscape in these eight plant lineages from both monocots and dicots. We show that retrotransposition has recurrently occurred in all plant genomes investigated, regardless their size, and through bursts, rather than a continuous process. Moreover, in each genome, only one or few LTR-RT families have been active in the recent past, and the difference in genome size among the species studied could thus mostly be accounted for by the extent of the latest transpositional burst(s). Following these bursts, LTR-RTs are efficiently eliminated from their host genomes through recombination and deletion, but we show that the removal rate is not lineage specific. These new findings lead us to propose a new model of TE-driven genome evolution in plants.

Isolation and characterization of reverse transcriptase fragments of LTR retrotransposons from the genome of Chenopodium quinoa (Amaranthaceae)

High heterogeneity was observed among conserved domains of reverse transcriptase ( rt ) isolated from quinoa. Only one Ty1- copia rt was highly amplified. Reverse transcriptase sequences were located predominantly in pericentromeric region of quinoa chromosomes. The heterogeneity, genomic abundance, and chromosomal distribution of reverse transcriptase (rt)-coding fragments of Ty1-copia and Ty3-gypsy long terminal repeat retrotransposons were analyzed in the Chenopodium quinoa genome. Conserved domains of the rt gene were amplified and characterized using degenerate oligonucleotide primer pairs. Sequence analyses indicated that half of Ty1-copia rt (51 %) and 39 % of Ty3-gypsy rt fragments contained intact reading frames. High heterogeneity among rt sequences was observed for both Ty1-copia and Ty3-gypsy rt amplicons, with Ty1-copia more heterogeneous than Ty3-gypsy. Most of the isolated rt fragments were present in quinoa genome in low copy numbers, with only one highly amplified Ty1-copia rt sequence family. The gypsy-like RNase H fragments co-amplified with Ty1-copia-degenerate primers were shown to be highly amplified in the quinoa genome indicating either higher abundance of some gypsy families of which rt domains could not be amplified, or independent evolution of this gypsy-region in quinoa. Both Ty1-copia and Ty3-gypsy retrotransposons were preferentially located in pericentromeric heterochromatin of quinoa chromosomes. Phylogenetic analyses of newly amplified rt fragments together with well-characterized retrotransposon families from other organisms allowed identification of major lineages of retroelements in the genome of quinoa and provided preliminary insight into their evolutionary dynamics.

Isolation and characterization of novel Ty1-copia-like retrotransposons from lily

Species of the genus Lilium are well known for their large genomes. Although expansion of noncoding repeated DNA is believed to account for this genome size, retroelement del Ty3-gypsy is the only one described so far in the genus Lilium. We isolated Ty1-copia elements from Lilium longiflorum and named them LIREs (lily retrotransposons). The long terminal repeats, primer binding site, and polypurine tract sequences are highly similar among the LIRE elements, indicating that they are in the same lineage. Although the protein-coding regions were highly decayed, the sequence motifs of the integrase, reverse transcriptase, and RNase H domains were identifiable as belonging to the order of Ty1-copia elements. Phylogenetic analysis and primer binding site sequences revealed that these elements belonged to the Ale lineage among the six lineages of plant Ty1-copia elements. Base substitutions in the long terminal repeats estimated that the integration times of the LIRE Ty1-copia elements were between 0.7 and 5.5 mya. In situ hybridization showed that the LIRE elements were present in all the chromosomes of L. longiflorum and L. lancifolium, but absent in centromeres, telomeres, and 45S rRNA sites in both species. The LIRE elements were present very abundantly in species of the genus Lilium, but absent in other genera of the family Liliaceae, implying that the LIRE elements might have contributed to the expansion of the genome in the genus Lilium.

Evolution of genome size in Carex (Cyperaceae) in relation to chromosome number and genomic base composition

BACKGROUND AND AIMS

The genus Carex exhibits karyological peculiarities related to holocentrism, specifically extremely broad and almost continual variation in chromosome number. However, the effect of these peculiarities on the evolution of the genome (genome size, base composition) remains unknown. While in monocentrics, determining the arithmetic relationship between the chromosome numbers of related species is usually sufficient for the detection of particular modes of karyotype evolution (i.e. polyploidy and dysploidy), in holocentrics where chromosomal fission and fusion occur such detection requires knowledge of the DNA content.

METHODS

The genome size and GC content were estimated in 157 taxa using flow cytometry. The exact chromosome numbers were known for 96 measured samples and were taken from the available literature for other taxa. All relationships were tested in a phylogenetic framework using the ITS tree of 105 species.

KEY RESULTS

The 1C genome size varied between 0·24 and 1·64 pg in Carex secalina and C. cuspidata, respectively. The genomic GC content varied from 34·8 % to 40·6 % from C. secalina to C. firma. Both genomic parameters were positively correlated. Seven polyploid and two potentially polyploid taxa were detected in the core Carex clade. A strong negative correlation between genome size and chromosome number was documented in non-polyploid taxa. Non-polyploid taxa of the core Carex clade exhibited a higher rate of genome-size evolution compared with the Vignea clade. Three dioecious taxa exhibited larger genomes, larger chromosomes, and a higher GC content than their hermaphrodite relatives.

CONCLUSIONS

Genomes of Carex are relatively small and very GC-poor compared with other angiosperms. We conclude that the evolution of genome and karyotype in Carex is promoted by frequent chromosomal fissions/fusions, rare polyploidy and common repetitive DNA proliferation/removal.

Genomic organization and dynamics of repetitive DNA sequences in representatives of three Fagaceae genera

Oaks, chestnuts, and beeches are economically important species of the Fagaceae. To understand the relationship between these members of this family, a deep knowledge of their genome composition and organization is needed. In this work, we have isolated and characterized several AFLP fragments obtained from Quercus rotundifolia Lam. through homology searches in available databases. Genomic polymorphisms involving some of these sequences were evaluated in two species of Quercus , one of Castanea , and one of Fagus with specific primers. Comparative FISH analysis with generated sequences was performed in interphase nuclei of the four species, and the co-immunolocalization of 5-methylcytosine was also studied. Some of the sequences isolated proved to be genus-specific, while others were present in all the genera. Retroelements, either gypsy-like of the Tat/Athila clade or copia-like, are well represented, and most are dispersed in euchromatic regions of these species with no DNA methylation associated, pointing to an interspersed arrangement of these retroelements with potential gene-rich regions. A particular gypsy-sequence is dispersed in oaks and chestnut nuclei, but its confinement to chromocenters in beech evidences genome restructuring events during evolution of Fagaceae. Several sequences generated in this study proved to be good tools to comparatively study Fagaceae genome organization.

The application of LTR retrotransposons as molecular markers in plants

Retrotransposons are a major agent of genome evolution. Various molecular marker systems have been developed that exploit the ubiquitous nature of these genetic elements and their property of stable integration into dispersed chromosomal loci that are polymorphic within species. The key methods, SSAP, IRAP, REMAP, RBIP, and ISBP, all detect the sites at which the retrotransposon DNA, which is conserved between families of elements, is integrated into the genome. Marker systems exploiting these methods can be easily developed and inexpensively deployed in the absence of extensive genome sequence data. They offer access to the dynamic and polymorphic, nongenic portion of the genome and thereby complement methods, such as gene-derived SNPs, that target primarily the genic fraction.

LTR retrotransposons contribute to genomic gigantism in plethodontid salamanders

Among vertebrates, most of the largest genomes are found within the salamanders, a clade of amphibians that includes 613 species. Salamander genome sizes range from ~14 to ~120 Gb. Because genome size is correlated with nucleus and cell sizes, as well as other traits, morphological evolution in salamanders has been profoundly affected by genomic gigantism. However, the molecular mechanisms driving genomic expansion in this clade remain largely unknown. Here, we present the first comparative analysis of transposable element (TE) content in salamanders. Using high-throughput sequencing, we generated genomic shotgun data for six species from the Plethodontidae, the largest family of salamanders. We then developed a pipeline to mine TE sequences from shotgun data in taxa with limited genomic resources, such as salamanders. Our summaries of overall TE abundance and diversity for each species demonstrate that TEs make up a substantial portion of salamander genomes, and that all of the major known types of TEs are represented in salamanders. The most abundant TE superfamilies found in the genomes of our six focal species are similar, despite substantial variation in genome size. However, our results demonstrate a major difference between salamanders and other vertebrates: salamander genomes contain much larger amounts of long terminal repeat (LTR) retrotransposons, primarily Ty3/gypsy elements. Thus, the extreme increase in genome size that occurred in salamanders was likely accompanied by a shift in TE landscape. These results suggest that increased proliferation of LTR retrotransposons was a major molecular mechanism contributing to genomic expansion in salamanders.

On the Coevolution of Transposable Elements and Plant Genomes

Plant genomes are unique in an intriguing feature: the range of their size variation is unprecedented among living organisms. Although polyploidization contributes to this variability, transposable elements (TEs) seem to play the pivotal role. TEs, often considered intragenomic parasites, not only affect the genome size of the host, but also interact with other genes, disrupting and creating new functions and regulatory networks. Coevolution of plant genomes and TEs has led to tight regulation of TE activity, and growing evidence suggests their relationship became mutualistic. Although the expansions of TEs represent certain costs for the host genomes, they may also bring profits for populations, helping to overcome challenging environmental (biotic/abiotic stress) or genomic (hybridization and allopolyploidization) conditions. In this paper, we discuss the possibility that the possession of inducible TEs may provide a selective advantage for various plant populations.

Genome-wide analysis of eukaryote thaumatin-like proteins (TLPs) with an emphasis on poplar

BACKGROUND

Plant inducible immunity includes the accumulation of a set of defense proteins during infection called pathogenesis-related (PR) proteins, which are grouped into families termed PR-1 to PR-17. The PR-5 family is composed of thaumatin-like proteins (TLPs), which are responsive to biotic and abiotic stress and are widely studied in plants. TLPs were also recently discovered in fungi and animals. In the poplar genome, TLPs are over-represented compared with annual species and their transcripts strongly accumulate during stress conditions.

RESULTS

Our analysis of the poplar TLP family suggests that the expansion of this gene family was followed by diversification, as differences in expression patterns and predicted properties correlate with phylogeny. In particular, we identified a clade of poplar TLPs that cluster to a single 350 kb locus of chromosome I and that are up-regulated by poplar leaf rust infection. A wider phylogenetic analysis of eukaryote TLPs - including plant, animal and fungi sequences - shows that TLP gene content and diversity increased markedly during land plant evolution. Mapping the reported functions of characterized TLPs to the eukaryote phylogenetic tree showed that antifungal or glycan-lytic properties are widespread across eukaryote phylogeny, suggesting that these properties are shared by most TLPs and are likely associated with the presence of a conserved acidic cleft in their 3D structure. Also, we established an exhaustive catalog of TLPs with atypical architectures such as small-TLPs, TLP-kinases and small-TLP-kinases, which have potentially developed alternative functions (such as putative receptor kinases for pathogen sensing and signaling).

CONCLUSION

Our study, based on the most recent plant genome sequences, provides evidence for TLP gene family diversification during land plant evolution. We have shown that the diverse functions described for TLPs are not restricted to specific clades but seem to be universal among eukaryotes, with some exceptions likely attributable to atypical protein structures. In the perennial plant model Populus, we unravelled the TLPs likely involved in leaf rust resistance, which will provide the foundation for further functional investigations.