Genome evolution of wild barley (Hordeum spontaneum) by BARE-1 retrotransposon dynamics in response to sharp microclimatic divergence

Chromosome-level genome assembly of Hippophae tibetana provides insights into high-altitude adaptation and flavonoid biosynthesis

BACKGROUND

As an endemic shrub of the Qinghai-Tibetan Plateau (QTP), the distribution of Hippophae tibetana Schlecht. ranges between 2800 and 5200 m above sea level. As the most basal branch of the Hippophae genus, H. tibetana has an extensive evolutionary history. The H. tibetana is a valuable tree for studying the ecological evolution of species under extreme conditions.

RESULTS

Here, we generated a high-quality chromosome-level genome of H. tibetana. The total size of the assembly genome is 917 Mb. The phylogenomic analysis of 1064 single-copy genes showed a divergence between 3.4 and 12.8 Mya for H. tibetana. Multiple gene families associated with DNA repair and disease resistance were significantly expanded in H. tibetana. We also identified many genes related to DNA repair with signs of positive selection. These results showed expansion and positive selection likely play important roles in H. tibetana's adaptation to comprehensive extreme environments in the QTP. A comprehensive genomic and transcriptomic analysis identified 49 genes involved in the flavonoid biosynthesis pathway in H. tibetana. We generated transgenic sea buckthorn hairy root producing high levels of flavonoid.

CONCLUSIONS

Taken together, this H. tibetana high-quality genome provides insights into the plant adaptation mechanisms of plant under extreme environments and lay foundation for the functional genomic research and molecular breeding of H. tibetana.

Genome size variation and evolution during invasive range expansion in an introduced plant

Plants demonstrate exceptional variation in genome size across species, and their genome sizes can also vary dramatically across individuals and populations within species. This aspect of genetic variation can have consequences for traits and fitness, but few studies attributed genome size differentiation to ecological and evolutionary processes. Biological invasions present particularly useful natural laboratories to infer selective agents that might drive genome size shifts across environments and population histories. Here, we test hypotheses for the evolutionary causes of genome size variation across 14 invading populations of yellow starthistle, Centaurea solstitialis, in California, United States. We use a survey of genome sizes and trait variation to ask: (1) Is variation in genome size associated with developmental trait variation? (2) Are genome sizes smaller toward the leading edge of the expansion, consistent with selection for "colonizer" traits? Or alternatively, does genome size increase toward the leading edge of the expansion, consistent with predicted consequences of founder effects and drift? (3) Finally, are genome sizes smaller at higher elevations, consistent with selection for shorter development times? We found that 2C DNA content varied 1.21-fold among all samples, and was associated with flowering time variation, such that plants with larger genomes reproduced later, with lower lifetime capitula production. Genome sizes increased toward the leading edge of the invasion, but tended to decrease at higher elevations, consistent with genetic drift during range expansion but potentially strong selection for smaller genomes and faster development time at higher elevations. These results demonstrate how genome size variation can contribute to traits directly tied to reproductive success, and how selection and drift can shape that variation. We highlight the influence of genome size on dynamics underlying a rapid range expansion in a highly problematic invasive plant.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

The genome of the glasshouse plant noble rhubarb (Rheum nobile) provides a window into alpine adaptation

Glasshouse plants are species that trap warmth via specialized morphology and physiology, mimicking a human glasshouse. In the Himalayan alpine region, the highly specialized glasshouse morphology has independently evolved in distinct lineages to adapt to intensive UV radiation and low temperature. Here we demonstrate that the glasshouse structure - specialized cauline leaves - is highly effective in absorbing UV light but transmitting visible and infrared light, creating an optimal microclimate for the development of reproductive organs. We reveal that this glasshouse syndrome has evolved at least three times independently in the rhubarb genus Rheum. We report the genome sequence of the flagship glasshouse plant Rheum nobile and identify key genetic network modules in association with the morphological transition to specialized glasshouse leaves, including active secondary cell wall biogenesis, upregulated cuticular cutin biosynthesis, and suppression of photosynthesis and terpenoid biosynthesis. The distinct cell wall organization and cuticle development might be important for the specialized optical property of glasshouse leaves. We also find that the expansion of LTRs has likely played an important role in noble rhubarb adaptation to high elevation environments. Our study will enable additional comparative analyses to identify the genetic basis underlying the convergent occurrence of glasshouse syndrome.

Drought and recovery in barley: key gene networks and retrotransposon response

Introduction

During drought, plants close their stomata at a critical soil water content (SWC), together with making diverse physiological, developmental, and biochemical responses.

Methods

Using precision-phenotyping lysimeters, we imposed pre-flowering drought on four barley varieties (Arvo, Golden Promise, Hankkija 673, and Morex) and followed their physiological responses. For Golden Promise, we carried out RNA-seq on leaf transcripts before and during drought and during recovery, also examining retrotransposon BARE1expression. Transcriptional data were subjected to network analysis.

Results

The varieties differed by their critical SWC (ϴ_crit), Hankkija 673 responding at the highest and Golden Promise at the lowest. Pathways connected to drought and salinity response were strongly upregulated during drought; pathways connected to growth and development were strongly downregulated. During recovery, growth and development pathways were upregulated; altogether, 117 networked genes involved in ubiquitin-mediated autophagy were downregulated.

Discussion

The differential response to SWC suggests adaptation to distinct rainfall patterns. We identified several strongly differentially expressed genes not earlier associated with drought response in barley. BARE1 transcription is strongly transcriptionally upregulated by drought and downregulated during recovery unequally between the investigated cultivars. The downregulation of networked autophagy genes suggests a role for autophagy in drought response; its importance to resilience should be further investigated.

Epigenomics in stress tolerance of plants under the climate change

BACKGROUND

Climate change has had a tremendous impact on the environment in general as well as agricultural crops grown in these situations as time passed. Agricultural production of crops is less suited and of lower quality due to disturbances in plant metabolism brought on by sensitivity to environmental stresses, which are brought on by climate change. Abiotic stressors that are specific to climate change, including as drought, extremes in temperature, increasing CO₂, waterlogging from heavy rain, metal toxicity, and pH changes, are known to negatively affect an array of species. Plants adapt to these challenges by undergoing genome-wide epigenetic changes, which are frequently accompanied by differences in transcriptional gene expression. The sum of a cell's biochemical modifications to its nuclear DNA, post-translational modifications to histones, and variations in the synthesis of non-coding RNAs is called an epigenome. These modifications frequently lead to variations in gene expression that occur without any alteration in the underlying base sequence.

EPIGENETIC MECHANISMS AND MARKS

The methylation of homologous loci by three different modifications-genomic (DNA methylation), chromatin (histone modifications), and RNA-directed DNA methylation (RdDM)-could be regarded as epigenetic mechanisms that control the regulation of differential gene expression. Stresses from the environment cause chromatin remodelling, which enables plant cells to adjust their expression patterns temporarily or permanently. EPIGENOMICS' CONSEQUENCES FOR GENOME STABILITY AND GENE EXPRESSION: DNA methylation affects gene expression in response to abiotic stressors by blocking or suppressing transcription. Environmental stimuli cause changes in DNA methylation levels, either upward in the case of hypermethylation or downward in the case of hypomethylation. The type of stress response that occurs as a result also affects the degree of DNA methylation alterations. Stress is also influenced by DRM2 and CMT3 methylating CNN, CNG, and CG. Both plant development and stress reactions depend on histone changes. Gene up-regulation is associated with histone tail phosphorylation, ubiquitination, and acetylation, while gene down-regulation is associated with de-acetylation and biotinylation. Plants undergo a variety of dynamic changes to histone tails in response to abiotic stressors. The relevance of these transcripts against stress is highlighted by the accumulation of numerous additional antisense transcripts, a source of siRNAs, caused by abiotic stresses. The study highlights the finding that plants can be protected from a range of abiotic stresses by epigenetic mechanisms such DNA methylation, histone modification, and RNA-directed DNA methylation. TRANSGENERATIONAL INHERITANCE AND SOURCES OF EPIGENETIC VARIATION: Stress results in the formation of epialleles, which are either transient or enduring epigenetic stress memory in plants. After the stress is gone, the stable memory is kept for the duration of the plant's remaining developmental cycles or passed on to the next generations, leading to plant evolution and adaptability. The bulk of epigenetic changes brought on by stress are temporary and return to normal after the stress has passed. Some of the modifications, however, might be long-lasting and transmitted across mitotic or even meiotic cell divisions. Epialleles often have genetic or non-genetic causes. Epialleles can arise spontaneously due to improper methylation state maintenance, short RNA off-target effects, or other non-genetic causes. Developmental or environmental variables that influence the stability of epigenetic states or direct chromatin modifications may also be non-genetic drivers of epigenetic variation. Transposon insertions that change local chromatin and structural rearrangements, such copy number changes that are genetically related or unrelated, are two genetic sources of epialleles.

EPIGENOMICS IN CROP IMPROVEMENT

To include epigenetics into crop breeding, it is necessary to create epigenetic variation as well as to identify and evaluate epialleles. Epigenome editing or epi-genomic selection may be required for epiallele creation and identification. In order to combat the challenges given by changing environments, these epigenetic mechanisms have generated novel epialleles that can be exploited to develop new crop types that are more climate-resilient. Numerous techniques can be used to alter the epigenome generally or at specific target loci in order to induce the epigenetic alterations necessary for crop development. Technologies like CRISPR/Cas9 and dCas, which have recently advanced, have opened up new avenues for the study of epigenetics. Epialleles could be employed in epigenomics-assisted breeding in addition to sequence-based markers for crop breeding.

CONCLUSIONS AND FUTURE PROSPECTUS

A few of the exciting questions that still need to be resolved in the area of heritable epigenetic variation include a better understanding of the epigenetic foundation of characteristics, the stability and heritability of epialleles, and the sources of epigenetic variation in crops. Investigating long intergenic non-coding RNAs (lincRNAs) as an epigenetic process might open up a new path to understanding crop plant's ability to withstand abiotic stress. For many of these technologies and approaches to be more applicable and deployable at a lower cost, technological breakthroughs will also be necessary. Breeders will probably need to pay closer attention to crop epialleles and how they can affect future responses to climate changes. The development of epialleles suitable for particular environmental circumstances may be made possible by creating targeted epigenetic changes in pertinent genes and by comprehending the molecular underpinnings of trans generational epigenetic inheritance. More research on a wider variety of plant species is required in order to fully comprehend the mechanisms that produce and stabilise epigenetic variation in crops. In addition to a collaborative and multidisciplinary effort by researchers in many fields of plant science, this will require a greater integration of the epigenomic data gathered in many crops. Before it may be applied generally, more study is required.

The power of retrotransposons in high-throughput genotyping and sequencing

The use of molecular markers has become an essential part of molecular genetics through their application in numerous fields, which includes identification of genes associated with targeted traits, operation of backcrossing programs, modern plant breeding, genetic characterization, and marker-assisted selection. Transposable elements are a core component of all eukaryotic genomes, making them suitable as molecular markers. Most of the large plant genomes consist primarily of transposable elements; variations in their abundance contribute to most of the variation in genome size. Retrotransposons are widely present throughout plant genomes, and replicative transposition enables them to insert into the genome without removing the original elements. Various applications of molecular markers have been developed that exploit the fact that these genetic elements are present everywhere and their ability to stably integrate into dispersed chromosomal localities that are polymorphic within a species. The ongoing development of molecular marker technologies is directly related to the deployment of high-throughput genotype sequencing platforms, and this research is of considerable significance. In this review, the practical application to molecular markers, which is a use of technology of interspersed repeats in the plant genome were examined using genomic sources from the past to the present. Prospects and possibilities are also presented.

Genome size variation and polyploidy prevalence in the genus Eragrostis are associated with the global dispersal in arid area

BACKGROUND

Biologists have long debated the drivers of the genome size evolution and variation ever since Darwin. Assumptions for the adaptive or maladaptive consequences of the associations between genome sizes and environmental factors have been proposed, but the significance of these hypotheses remains controversial. Eragrostis is a large genus in the grass family and is often used as crop or forage during the dry seasons. The wide range and complex ploidy levels make Eragrostis an excellent model for investigating how the genome size variation and evolution is associated with environmental factors and how these changes can ben interpreted.

METHODS

We reconstructed the Eragrostis phylogeny and estimated genome sizes through flow cytometric analyses. Phylogenetic comparative analyses were performed to explore how genome size variation and evolution is related to their climatic niches and geographical ranges. The genome size evolution and environmental factors were examined using different models to study the phylogenetic signal, mode and tempo throughout evolutionary history.

RESULTS

Our results support the monophyly of Eragrostis. The genome sizes in Eragrostis ranged from ~0.66 pg to ~3.80 pg. We found that a moderate phylogenetic conservatism existed in terms of the genome sizes but was absent from environmental factors. In addition, phylogeny-based associations revealed close correlations between genome sizes and precipitation-related variables, indicating that the genome size variation mainly caused by polyploidization may have evolved as an adaptation to various environments in the genus Eragrostis.

CONCLUSION

This is the first study to take a global perspective on the genome size variation and evolution in the genus Eragrostis. Our results suggest that the adaptation and conservatism are manifested in the genome size variation, allowing the arid species of Eragrostis to spread the xeric area throughout the world.

The role of LTR retrotransposons in plant genetic engineering: how to control their transposition in the genome

We briefly discuss that the similarity of LTR retrotransposons to retroviruses is a great opportunity for the development of a genetic engineering tool that exploits intragenic elements in the plant genome for plant genetic improvement. Long terminal repeat (LTR) retrotransposons are very similar to retroviruses but do not have the property of being infectious. While spreading between its host cells, a retrovirus inserts a DNA copy of its genome into the cells. The ability of retroviruses to cause infection with genome integration allows genes to be delivered to cells and tissues. Retrovirus vectors are, however, only specific to animals and insects, and, thus, are not relevant to plant genetic engineering. However, the similarity of LTR retrotransposons to retroviruses is an opportunity to explore the former as a tool for genetic engineering. Although recent long-read sequencing technologies have advanced the knowledge about transposable elements (TEs), the integration of TEs is still unable either to control them or to direct them to specific genomic locations. The use of existing intragenic elements to achieve the desired genome composition is better than using artificial constructs like vectors, but it is not yet clear how to control the process. Moreover, most LTR retrotransposons are inactive and unable to produce complete proteins. They are also highly mutable. In addition, it is impossible to find a full active copy of a LTR retrotransposon out of thousands of its own copies. Theoretically, if these elements were directly controlled and turned on or off using certain epigenetic mechanisms (inducing by stress or infection), LTR retrotransposons could be a great opportunity to develop a genetic engineering tool using intragenic elements in the plant genome. In this review, the recent developments in uncovering the nature of LTR retrotransposons and the possibility of using these intragenic elements as a tool for plant genetic engineering are briefly discussed.

The ecology of palm genomes: repeat-associated genome size expansion is constrained by aridity

Genome size varies 2400-fold across plants, influencing their evolution through changes in cell size and cell division rates which impact plants' environmental stress tolerance. Repetitive element expansion explains much genome size diversity, and the processes structuring repeat 'communities' are analogous to those structuring ecological communities. However, which environmental stressors influence repeat community dynamics has not yet been examined from an ecological perspective. We measured genome size and leveraged climatic data for 91% of genera within the ecologically diverse palm family (Arecaceae). We then generated genomic repeat profiles for 141 palm species, and analysed repeats using phylogenetically informed linear models to explore relationships between repeat dynamics and environmental factors. We show that palm genome size and repeat 'community' composition are best explained by aridity. Specifically, Ty3-gypsy and TIR elements were more abundant in palm species from wetter environments, which generally had larger genomes, suggesting amplification. By contrast, Ty1-copia and LINE elements were more abundant in drier environments. Our results suggest that water stress inhibits repeat expansion through selection on upper genome size limits. However, elements that may associate with stress-response genes (e.g. Ty1-copia) have amplified in arid-adapted palm species. Overall, we provide novel evidence of climate influencing the assembly of repeat 'communities'.

Do water and soil nutrient scarcities differentially impact the performance of diploid and tetraploid Solidago gigantea (Giant Goldenrod, Asteraceae)?

Plants require water and nutrients for survival, although the effects of their availabilities on plant fitness differ amongst species. Genome size variation, within and across species, is suspected to influence plant water and nutrient requirements, but little is known about how variations in these resources concurrently affect plant fitness based on genome size. We examined how genome size variation between autopolyploid cytotypes influences plant morphological and physiological traits, and whether cytotype-specific trait responses differ based on water and/or nutrient availability. Diploid and autotetraploid Solidago gigantea (Giant Goldenrod) were grown in a greenhouse under four soil water:N+P treatments (L:L, L:H, H:L, H:H), and stomata characteristics (size, density), growth (above- and belowground biomass, R/S), and physiological (A_net , E, WUE) responses were measured. Resource availabilities and cytotype identity influenced some plant responses but their effects were independent of each other. Plants grown in high-water and nutrient treatments were larger, plants grown in low-water or high-nutrient treatments had higher WUE but lower E, and A_net and E rates decreased as plants aged. Autotetraploids also had larger and fewer stomata, higher biomass and larger A_net than diploids. Nutrient and water availability could influence intra- and interspecific competitive outcomes. Although S. gigantea cytotypes were not differentially affected by resource treatments, genome size may influence cytogeographic range patterning and population establishment likelihood. For instance, the larger size of autotetraploid S. gigantea might render them more competitive for resources and niche space than diploids.

Long terminal repeats (LTR) and transcription factors regulate PHRE1 and PHRE2 activity in Moso bamboo under heat stress

BACKGROUND

LTR retrotransposons play a significant role in plant growth, genome evolution, and environmental stress response, but their regulatory response to heat stress remains unclear. We have investigated the activities of two LTR retrotransposons, PHRE1 and PHRE2, of moso bamboo (Phyllostachys edulis) in response to heat stress.

RESULTS

The differential overexpression of PHRE1 and PHRE2 with or without CaMV35s promoter showed enhanced expression under heat stress in transgenic plants. The transcriptional activity studies showed an increase in transposition activity and copy number among moso bamboo wild type and Arabidopsis transgenic plants under heat stress. Comparison of promoter activity in transgenic plants indicated that 5'LTR promoter activity was higher than CaMV35s promoter. Additionally, yeast one-hybrid (Y1H) system and in planta biomolecular fluorescence complementation (BiFC) assay revealed interactions of heat-dependent transcription factors (TFs) with 5'LTR sequence and direct interactions of TFs with pol and gag.

CONCLUSIONS

Our results conclude that the 5'LTR acts as a promoter and could regulate the LTR retrotransposons in moso bamboo under heat stress.

Primer Binding Site (PBS) Profiling of Genetic Diversity of Natural Populations of Endemic Species Allium ledebourianum Schult

Endemic species are especially vulnerable to biodiversity loss caused by isolation or habitat specificity, small population size, and anthropogenic factors. Endemic species biodiversity analysis has a critically important global value for the development of conservation strategies. The rare onion Allium ledebourianum is a narrow-lined endemic species, with natural populations located in the extreme climatic conditions of the Kazakh Altai. A. ledebourianum populations are decreasing everywhere due to anthropogenic impact, and therefore, this species requires preservation and protection. Conservation of this rare species is associated with monitoring studies to investigate the genetic diversity of natural populations. Fundamental components of eukaryote genome include multiple classes of interspersed repeats. Various PCR-based DNA fingerprinting methods are used to detect chromosomal changes related to recombination processes of these interspersed elements. These methods are based on interspersed repeat sequences and are an effective approach for assessing the biological diversity of plants and their variability. We applied DNA profiling approaches based on conservative sequences of interspersed repeats to assess the genetic diversity of natural A. ledebourianum populations located in the territory of Kazakhstan Altai. The analysis of natural A. ledebourianum populations, carried out using the DNA profiling approach, allowed the effective differentiation of the populations and assessment of their genetic diversity. We used conservative sequences of tRNA primer binding sites (PBS) of the long-terminal repeat (LTR) retrotransposons as PCR primers. Amplification using the three most effective PBS primers generated 628 PCR amplicons, with an average of 209 amplicons. The average polymorphism level varied from 34% to 40% for all studied samples. Resolution analysis of the PBS primers showed all of them to have high or medium polymorphism levels, which varied from 0.763 to 0.965. Results of the molecular analysis of variance showed that the general biodiversity of A. ledebourianum populations is due to interpopulation (67%) and intrapopulation (33%) differences. The revealed genetic diversity was higher in the most distant population of A. ledebourianum LD64, located on the Sarymsakty ridge of Southern Altai. This is the first genetic diversity study of the endemic species A. ledebourianum using DNA profiling approaches. This work allowed us to collect new genetic data on the structure of A. ledebourianum populations in the Altai for subsequent development of preservation strategies to enhance the reproduction of this relict species. The results will be useful for the conservation and exploitation of this species, serving as the basis for further studies of its evolution and ecology.

Transposable element landscapes illuminate past evolutionary events in the endangered fern Vandenboschia speciosa

Vandenboschia speciosa is an endangered tetraploid fern species with a large genome (10.5 Gb). Its geographical distribution is characterized by disjoined tertiary flora refuges, with relict populations that survived past climate crises. Here, we analyzed the transposable elements (TEs) and found that they comprise approximately 76% of the V. speciosa genome, thus being the most abundant type of DNA sequence in this gigantic genome. The V. speciosa genome is composed of 51% and 5.6% of Class I and Class II elements, respectively. LTR retrotransposons were the most abundant TEs in this species (at least 42% of the genome), followed by non-LTR retrotransposons, which constituted at least 8.7% of the genome of this species. We introduce an additional analysis to identify the nature of non-annotated elements (19% of the genome). A BLAST search of the non-annotated contigs against the V. speciosa TE database allowed for the identification of almost half of them, which were most likely diverged sequence variants of the annotated TEs. In general, the TE composition in V. speciosa resembles the TE composition in seed plants. In addition, repeat landscapes revealed three episodes of amplification for all TEs, most likely due to demographic changes associated with past climate crises.

Genome Size Dynamics in Marine Ribbon Worms (Nemertea, Spiralia)

Nemertea is a phylum consisting of 1300 mostly marine species. Nemertea is distinguished by an eversible muscular proboscis, and most of the species are venomous. Genomic resources for this phylum are scarce despite their value in understanding biodiversity. Here, we present genome size estimates of Nemertea based on flow cytometry and their relationship to different morphological and developmental traits. Ancestral genome size estimations were done across the nemertean phylogeny. The results increase the available genome size estimates for Nemertea three-fold. Our analyses show that Nemertea has a narrow genome size range (0.43-3.89 pg) compared to other phyla in Lophotrochozoa. A relationship between genome size and evolutionary rate, developmental modes, and habitat was found. Trait analyses show that the highest evolutionary rate of genome size is found in upper intertidal, viviparous species with direct development. Despite previous findings, body size in nemerteans was not correlated with genome size. A relatively small genome (1.18 pg) is assumed for the most recent common ancestor of all extant nemerteans. The results provide an important basis for future studies in nemertean genomics, which will be instrumental to understanding the evolution of this enigmatic and often neglected phylum.

Genomic basis of high-altitude adaptation in Tibetan Prunus fruit trees

The Great Himalayan Mountains and their foothills are believed to be the place of origin and development of many plant species. The genetic basis of adaptation to high plateaus is a fascinating topic that is poorly understood at the population level. We comprehensively collected and sequenced 377 accessions of Prunus germplasm along altitude gradients ranging from 2,067 to 4,492 m in the Himalayas. We de novo assembled three high-quality genomes of Tibetan Prunus species. A comparative analysis of Prunus genomes indicated a remarkable expansion of the SINE retrotransposons occurred in the genomes of Tibetan species. We observed genetic differentiation between Tibetan peaches from high and low altitudes and that genes associated with light stress signaling, especially UV stress signaling, were enriched in the differentiated regions. By profiling the metabolomes of Tibetan peach fruit, we determined 379 metabolites had significant genetic correlations with altitudes and that in particular phenylpropanoids were positively correlated with altitudes. We identified 62 Tibetan peach-specific SINEs that colocalized with metabolites differentially accumualted in Tibetan relative to cultivated peach. We demonstrated that two SINEs were inserted in a locus controlling the accumulation of 3-O-feruloyl quinic acid. SINE1 was specific to Tibetan peach. SINE2 was predominant in high altitudes and associated with the accumulation of 3-O-feruloyl quinic acid. These genomic and metabolic data for Prunus populations native to the Himalayan region indicate that the expansion of SINE retrotransposons helped Tibetan Prunus species adapt to the harsh environment of the Himalayan plateau by promoting the accumulation of beneficial metabolites.

Genome size evolution is associated with climate seasonality and glucosinolates, but not life history, soil nutrients or range size, across a clade of mustards

BACKGROUND AND AIMS

We investigate patterns of evolution of genome size across a morphologically and ecologically diverse clade of Brassicaceae, in relation to ecological and life history traits. While numerous hypotheses have been put forward regarding autecological and environmental factors that could favour small vs. large genomes, a challenge in understanding genome size evolution in plants is that many hypothesized selective agents are intercorrelated.

METHODS

We contribute genome size estimates for 47 species of Streptanthus Nutt. and close relatives, and take advantage of many data collections for this group to assemble data on climate, life history, soil affinity and composition, geographic range and plant secondary chemistry to identify simultaneous correlates of variation in genome size in an evolutionary framework. We assess models of evolution across clades and use phylogenetically informed analyses as well as model selection and information criteria approaches to identify variables that can best explain genome size variation in this clade.

KEY RESULTS

We find differences in genome size and heterogeneity in its rate of evolution across subclades of Streptanthus and close relatives. We show that clade-wide genome size is positively associated with climate seasonality and glucosinolate compounds. Model selection and information criteria approaches identify a best model that includes temperature seasonality and fraction of aliphatic glucosinolates, suggesting a possible role for genome size in climatic adaptation or a role for biotic interactions in shaping the evolution of genome size. We find no evidence supporting hypotheses of life history, range size or soil nutrients as forces shaping genome size in this system.

CONCLUSIONS

Our findings suggest climate seasonality and biotic interactions as potential forces shaping the evolution of genome size and highlight the importance of evaluating multiple factors in the context of phylogeny to understand the effect of possible selective agents on genome size.

Androgenic-Induced Transposable Elements Dependent Sequence Variation in Barley

A plant genome usually encompasses different families of transposable elements (TEs) that may constitute up to 85% of nuclear DNA. Under stressful conditions, some of them may activate, leading to sequence variation. In vitro plant regeneration may induce either phenotypic or genetic and epigenetic changes. While DNA methylation alternations might be related, i.e., to the Yang cycle problems, DNA pattern changes, especially DNA demethylation, may activate TEs that could result in point mutations in DNA sequence changes. Thus, TEs have the highest input into sequence variation (SV). A set of barley regenerants were derived via in vitro anther culture. High Performance Liquid Chromatography (RP-HPLC), used to study the global DNA methylation of donor plants and their regenerants, showed that the level of DNA methylation increased in regenerants by 1.45% compared to the donors. The Methyl-Sensitive Transposon Display (MSTD) based on methylation-sensitive Amplified Fragment Length Polymorphism (metAFLP) approach demonstrated that, depending on the selected elements belonging to the TEs family analyzed, varying levels of sequence variation were evaluated. DNA sequence contexts may have a different impact on SV generated by distinct mobile elements belonged to various TE families. Based on the presented study, some of the selected mobile elements contribute differently to TE-related SV. The surrounding context of the TEs DNA sequence is possibly important here, and the study explained some part of SV related to those contexts.

Chromosome-scale genome assembly of areca palm (Areca catechu)

Areca palm (Areca catechu L.; family Arecaceae) is an important tropical medicinal crop and is also used for masticatory and religious purposes in Asia. Improvements to areca properties made by traditional breeding tools have been very slow, and further advances in its cultivation and practical use require genomic information, which is still unavailable. Here, we present a chromosome-scale reference genome assembly for areca by combining Illumina and PacBio data with Hi-C mapping technologies, covering the predicted A. catechu genome length (2.59 Gb, variety "Reyan#1") to an estimated 240× read depth. The assembly was 2.51 Gb in length with a scaffold N50 of 1.7Mb. The scaffolds were then further assembled into 16 pseudochromosomes, with an N50 of 172 Mb. Transposable elements comprised 80.37% of the areca genome, and 68.68% of them were long-terminal repeat retrotransposon elements. The areca palm genome was predicted to harbour 31,571 protein-coding genes and overall, 92.92% of genes were functionally annotated, including enriched and expanded families of genes responsible for biosynthesis of flavonoid, anthocyanin, monoterpenoid and their derivatives. Comparative analyses indicated that A. catechu probably diverged from its close relatives Elaeis guineensis and Cocos nucifera approximately 50.3 million years ago (Ma). Two whole genome duplication events in areca palm were found to be shared by palms and monocots, respectively. This genome assembly and associated resources represents an important addition to the palm genomics community and will be a valuable resource that will facilitate areca palm breeding and improve our understanding of areca palm biology and evolution.

Gaining Acceptance of Novel Plant Breeding Technologies

Ensuring the sustainability of agriculture under climate change has led to a surge in alternative strategies for crop improvement. Advances in integrated crop breeding, social acceptance, and farm-level adoption are crucial to address future challenges to food security. Societal acceptance can be slow when consumers do not see the need for innovation or immediate benefits. We consider how best to address the issue of social licence and harmonised governance for novel gene technologies in plant breeding. In addition, we highlight optimised breeding strategies that will enable long-term genetic gains to be achieved. Promoted by harmonised global policy change, innovative plant breeding can realise high and sustainable productivity together with enhanced nutritional traits.

Genetic Variation of the Serine Acetyltransferase Gene Family for Sulfur Assimilation in Maize

Improving sulfur assimilation in maize kernels is essential due to humans and animals' inability to synthesize methionine. Serine acetyltransferase (SAT) is a critical enzyme that controls cystine biosynthesis in plants. In this study, all SAT gene members were genome-wide characterized by using a sequence homology search. The RNA-seq quantification indicates that they are highly expressed in leaves, other than root and seeds, consistent with their biological functions in sulfur assimilation. With the recently released 25 genomes of nested association mapping (NAM) founders representing the diverse maize stock, we had the opportunity to investigate the SAT genetic variation comprehensively. The abundant transposon insertions into SAT genes indicate their driving power in terms of gene structure and genome evolution. We found that the transposon insertion into exons could change SAT gene transcription, whereas there was no significant correlation between transposable element (TE) insertion into introns and their gene expression, indicating that other regulatory elements such as promoters could also be involved. Understanding the SAT gene structure, gene expression and genetic variation involved in natural selection and species adaption could precisely guide genetic engineering to manipulate sulfur assimilation in maize and to improve nutritional quality.

Clustered and dispersed chromosomal distribution of the two classes of Revolver transposon family in rye (Secale cereale)

The chromosomal locations of a new class of Revolver transposon-like elements were analyzed by using FISH method on the metaphase chromosome in somatic cell division of the rye cultivar Petkus. First, the Revolver standard element probe λ2 was weakly hybridized throughout the rye chromosome, and comparatively large interstitial signals spotted with a dot shape were detected together with several telomeric regions. The dot shape interstitial signal was stably detected at one site on Chromosome (Chr) 1R (middle part of the interstitial region of the short arm), three sites on Chr 2R (distal part of the interstitial region and adjacent to the centromere on the short arm, middle part of the interstitial region of the long arm), and two sites on Chr 5R (middle part of the interstitial region and adjacent to the centromere on the long arm). The Revolver λ2 probe was effective for identification of 1R, 2R, and 5R chromosomes. On the other hand, Revolver nonautonomous element-specific L626-BARE-100 probe was strongly distributed throughout the rye chromosomes, and considerable numbers and diverse lengths of transcripts were detected by RT-PCR. Although the standard elements were found in localized clusters, the nonautonomous elements tended to be dispersed throughout the genome. Clustered nature of Revolver is a significantly rare case in genomics.

Transcriptional activity of transposable elements along an elevational gradient in Arabidopsis arenosa

Background

Plant genomes can respond rapidly to environmental changes and transposable elements (TEs) arise as important drivers contributing to genome dynamics. Although some elements were reported to be induced by various abiotic or biotic factors, there is a lack of general understanding on how environment influences the activity and diversity of TEs. Here, we combined common garden experiment with short-read sequencing to investigate genomic abundance and expression of 2245 consensus TE sequences (containing retrotransposons and DNA transposons) in an alpine environment in Arabidopsis arenosa. To disentangle general trends from local differentiation, we leveraged four foothill-alpine population pairs from different mountain regions. Seeds of each of the eight populations were raised under four treatments that differed in temperature and irradiance, two factors varying with elevation. RNA-seq analysis was performed on leaves of young plants to test for the effect of elevation and subsequently of temperature and irradiance on expression of TE sequences.

Results

Genomic abundance of the 2245 consensus TE sequences varied greatly between the mountain regions in line with neutral divergence among the regions, representing distinct genetic lineages of A. arenosa. Accounting for intraspecific variation in abundance, we found consistent transcriptomic response for some TE sequences across the different pairs of foothill-alpine populations suggesting parallelism in TE expression. In particular expression of retrotransposon LTR Copia (e.g. Ivana and Ale clades) and LTR Gypsy (e.g. Athila and CRM clades) but also non-LTR LINE or DNA transposon TIR MuDR consistently varied with elevation of origin. TE sequences responding specifically to temperature and irradiance belonged to the same classes as well as additional TE clades containing potentially stress-responsive elements (e.g. LTR Copia Sire and Tar, LTR Gypsy Reina).

Conclusions

Our study demonstrated that the A. arenosa genome harbours a considerable diversity of TE sequences whose abundance and expression response varies across its native range. Some TE clades may contain transcriptionally active elements responding to a natural environmental gradient. This may further contribute to genetic variation between populations and may ultimately provide new regulatory mechanisms to face environmental challenges.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00236-0.

The ecology of the genome and the dynamics of the biological dark matter

Transposable elements (TEs) are essential components of the eukaryotic genomes. While mostly deleterious, evidence is mounting that TEs provide the host with beneficial adaptations. How 'selfish' or 'parasitic' DNA persists until it helps species evolution is emerging as a major evolutionary puzzle, especially in asexual taxa where the lack of sex strongly impede the spread of TEs. Since occasional but unchecked TE proliferations would ultimately drive host lineages toward extinction, asexual genomes are typically predicted to be free of TEs, which contrasts with their persistence in asexual taxa. We designed innovative 'Eco-genomic' models that account for both host demography and within-host molecular mechanisms of transposition and silencing to analyze their impact on TE dynamics in asexual genome populations. We unraveled that the spread of TEs can be limited to a stable level by density-dependent purifying selection when TE copies are over-dispersed among lineages and the host demographic turn-over is fast. We also showed that TE silencing can protect host populations in two ways; by preventing TEs with weak effects to accumulate or by favoring the elimination of TEs with large effects. Our predictions may explain TE persistence in known asexual taxa that typically show fast demography and where TE copy number variation between lineages is expected. Such TE persistence in asexual taxa potentially has important implications for their evolvability and the preservation of sexual reproduction.

Transposable Element Landscape in Drosophila Populations Selected for Longevity

Transposable elements (TEs) inflict numerous negative effects on health and fitness as they replicate by integrating into new regions of the host genome. Even though organisms employ powerful mechanisms to demobilize TEs, transposons gradually lose repression during aging. The rising TE activity causes genomic instability and was implicated in age-dependent neurodegenerative diseases, inflammation, and the determination of lifespan. It is therefore conceivable that long-lived individuals have improved TE silencing mechanisms resulting in reduced TE expression relative to their shorter-lived counterparts and fewer genomic insertions. Here, we test this hypothesis by performing the first genome-wide analysis of TE insertions and expression in populations of Drosophila melanogaster selected for longevity through late-life reproduction for 50–170 generations from four independent studies. Contrary to our expectation, TE families were generally more abundant in long-lived populations compared with nonselected controls. Although simulations showed that this was not expected under neutrality, we found little evidence for selection driving TE abundance differences. Additional RNA-seq analysis revealed a tendency for reducing TE expression in selected populations, which might be more important for lifespan than regulating genomic insertions. We further find limited evidence of parallel selection on genes related to TE regulation and transposition. However, telomeric TEs were genomically and transcriptionally more abundant in long-lived flies, suggesting improved telomere maintenance as a promising TE-mediated mechanism for prolonging lifespan. Our results provide a novel viewpoint indicating that reproduction at old age increases the opportunity of TEs to be passed on to the next generation with little impact on longevity.

Coconut genome assembly enables evolutionary analysis of palms and highlights signaling pathways involved in salt tolerance

Coconut (Cocos nucifera) is the emblematic palm of tropical coastal areas all around the globe. It provides vital resources to millions of farmers. In an effort to better understand its evolutionary history and to develop genomic tools for its improvement, a sequence draft was recently released. Here, we present a dense linkage map (8402 SNPs) aiming to assemble the large genome of coconut (2.42 Gbp, 2n = 32) into 16 pseudomolecules. As a result, 47% of the sequences (representing 77% of the genes) were assigned to 16 linkage groups and ordered. We observed segregation distortion in chromosome Cn15, which is a signature of strong selection among pollen grains, favouring the maternal allele. Comparing our results with the genome of the oil palm Elaeis guineensis allowed us to identify major events in the evolutionary history of palms. We find that coconut underwent a massive transposable element invasion in the last million years, which could be related to the fluctuations of sea level during the glaciations at Pleistocene that would have triggered a population bottleneck. Finally, to better understand the facultative halophyte trait of coconut, we conducted an RNA-seq experiment on leaves to identify key players of signaling pathways involved in salt stress response. Altogether, our findings represent a valuable resource for the coconut breeding community.

DNA profiling and assessment of genetic diversity of relict species Allium altaicum Pall. on the territory of Altai

Analysis of the genetic diversity of natural populations of threatened and endangered species of plants is a main aspect of conservation strategy. The endangered species Allium altaicum is a relict plant of the Ice Age and natural populations are located in extreme climatic conditions of Kazakstan's Altai Mountains. Mobile genetic elements and other interspersed repeats are basic components of a eukaryote genome, which can activate under stress conditions and indirectly promote the survival of an organism against environmental stresses. Detections of chromosomal changes related to recombination processes of mobile genetic elements are performed by various PCR methods. These methods are based on interspersed repeat sequences and are an effective tool for research of biological diversity of plants and their variability. In our research, we used conservative sequences of tRNA primer binding sites (PBS) when initializing the retrotransposon replication as PCR primers to research the genetic diversity of 12 natural populations of A. altaicum found in various ecogeographic conditions of the Kazakhstani Altai. High efficiency of the PBS amplification method used was observed already at the intrapopulation level. Unique amplicons representative of a certain population were found at the intrapopulation level. Analysis of molecular dispersion revealed that the biodiversity of populations of mountainous and lowland A. altaicum is due to intrapopulation differences for climatic zones of habitation. This is likely conditional upon predominance of vegetative reproduction over seed reproduction in some populations. In the case of vegetative reproduction, somatic recombination related to the activity of mobile genetic elements are preserved in subsequent generations. This leads to an increase of intrapopulation genetic diversity. Thus, high genetic diversity was observed in populations such as A. altaicum located in the territory of the Kalbinskii Altai, whereas the minimum diversity was observed in the populations of the Leninororsk ecogeographic group. Distinctions between these populations were also identified depending on the areas of their distribution. Low-land and mid-mountain living environments are characterized by a great variety of shapes and plasticity. This work allowed us to obtain new genetic data on the structure of A. altaicum populations on the territory of the Kazakhstan Altai for the subsequent development of preservation and reproduction strategies for this relict species.

Evolutionary dynamics of transposable elements and satellite DNAs in polyploid Spartina species

Repeated sequences and polyploidy play a central role in plant genome dynamics. Here, we analyze the evolutionary dynamics of repeats in tetraploid and hexaploid Spartina species that diverged during the last 10 million years within the Chloridoideae, one of the poorest investigated grass lineages. From high-throughput genome sequencing, we annotated Spartina repeats and determined what sequence types account for the genome size variation among species. We examined whether differential genome size evolution correlated with ploidy levels and phylogenetic relationships. We also examined the tempo of repeat sequence dynamics associated with allopatric speciation over the last 3-6 million years between hexaploid species that diverged on the American and European Atlantic coasts and tetraploid species from North and South America. The tetraploid S. spartinae, whose phylogenetic placement has been debated, exhibits a similar repeat content as hexaploid species, suggesting common ancestry. Genome expansion or contraction resulting from repeat dynamics seems to be explained mostly by the contrasting divergence times between species, rather than by genome changes triggered by ploidy level change per se. One 370 bp satellite may be exhibiting 'meiotic drive' and driving chromosome evolution in S. alterniflora. Our results provide crucial insights for investigating the genetic and epigenetic consequences of such differential repeat dynamics on the ecology and distribution of the meso- and neopolyploid Spartina species.

Satellite DNA-Mediated Gene Expression Regulation: Physiological and Evolutionary Implication

Satellite DNAs are tandemly repeated sequences organized in large clusters within (peri)centromeric and/or subtelomeric heterochromatin. However, in many species, satellite DNAs are not restricted to heterochromatin but are also dispersed as short arrays within euchromatin. Such genomic organization together with transcriptional activity seems to be a prerequisite for the gene-modulatory effect of satellite DNAs which was first demonstrated in the beetle Tribolium castaneum upon heat stress. Namely, enrichment of a silent histone mark at euchromatic repeats of a major beetle satellite DNA results in epigenetic silencing of neighboring genes. In addition, human satellite III transcripts induced by heat shock contribute to genome-wide gene silencing, providing protection against stress-induced cell death. Gene silencing mediated by satellite RNA was also shown to be fundamental for the early embryonic development of the mosquito Aedes aegypti. Apart from a physiological role during embryogenesis and heat stress response, activation of satellite DNAs in terms of transcription and proliferation can have an evolutionary impact. Spreading of satellite repeats throughout euchromatin promotes the variation of epigenetic landscapes and gene expression diversity, contributing to the evolution of gene regulatory networks and to genome adaptation in fluctuating environmental conditions.

Retrotransposon-based genetic diversity of Deschampsia antarctica Desv. from King George Island (Maritime Antarctic)

Deschampsia antarctica Desv. can be found in diverse Antarctic habitats which may vary considerably in terms of environmental conditions and soil properties. As a result, the species is characterized by wide ecotypic variation in terms of both morphological and anatomical traits. The species is a unique example of an organism that can successfully colonize inhospitable regions due to its phenomenal ability to adapt to both the local mosaic of microhabitats and to general climatic fluctuations. For this reason, D. antarctica has been widely investigated in studies analyzing morphophysiological and biochemical responses to various abiotic stresses (frost, drought, salinity, increased UV radiation). However, there is little evidence to indicate whether the observed polymorphism is accompanied by the corresponding genetic variation. In the present study, retrotransposon-based iPBS markers were used to trace the genetic variation of D. antarctica collected in nine sites of the Arctowski oasis on King George Island (Western Antarctic). The genotyping of 165 individuals from nine populations with seven iPBS primers revealed 125 amplification products, 15 of which (12%) were polymorphic, with an average of 5.6% polymorphic fragments per population. Only one of the polymorphic fragments, observed in population 6, was represented as a private band. The analyzed specimens were characterized by low genetic diversity (uHe = 0.021, I = 0.030) and high population differentiation (F ST = 0.4874). An analysis of Fu's F S statistics and mismatch distribution in most populations (excluding population 2, 6 and 9) revealed demographic/spatial expansion, whereas significant traces of reduction in effective population size were found in three populations (1, 3 and 5). The iPBS markers revealed genetic polymorphism of D. antarctica, which could be attributed to the mobilization of random transposable elements, unique features of reproductive biology, and/or geographic location of the examined populations.

Evolutionary Dynamics of Transposable Elements Following a Shared Polyploidization Event in the Tribe Andropogoneae

Both polyploidization and transposable element (TE) activity are known to be major drivers of plant genome evolution. Here, we utilize the Zea-Tripsacum clade to investigate TE activity and accumulation after a shared polyploidization event. Comparisons of TE evolutionary dynamics in various Zea and Tripsacum species, in addition to two closely related diploid species, Urelytrum digitatum and Sorghum bicolor, revealed variation in repeat content among all taxa included in the study. The repeat composition of Urelytrum is more similar to that of Zea and Tripsacum compared to Sorghum, despite the similarity in genome size with the latter. Although LTR-retrotransposons were abundant in all species, we observed an expansion of the copia superfamily, specifically in Z. mays and T. dactyloides, species that have adapted to more temperate environments. Additional analyses of the genomic distribution of these retroelements provided evidence of biased insertions near genes involved in various biological processes including plant development, defense, and macromolecule biosynthesis. Specifically, copia insertions in Zea and T. dactyloides were significantly enriched near genes involved in abiotic stress response, suggesting independent evolution post Zea-Tripsacum divergence. The lack of copia insertions near the orthologous genes in S. bicolor suggests that duplicate gene copies generated during polyploidization may offer novel neutral sites for TEs to insert, thereby providing an avenue for subfunctionalization via TE insertional mutagenesis.

Comparative Study of Pine Reference Genomes Reveals Transposable Element Interconnected Gene Networks

Sequencing the giga-genomes of several pine species has enabled comparative genomic analyses of these outcrossing tree species. Previous studies have revealed the wide distribution and extraordinary diversity of transposable elements (TEs) that occupy the large intergenic spaces in conifer genomes. In this study, we analyzed the distribution of TEs in gene regions of the assembled genomes of Pinus taeda and Pinus lambertiana using high-performance computing resources. The quality of draft genomes and the genome annotation have significant consequences for the investigation of TEs and these aspects are discussed. Several TE families frequently inserted into genes or their flanks were identified in both species’ genomes. Potentially important sequence motifs were identified in TEs that could bind additional regulatory factors, promoting gene network formation with faster or enhanced transcription initiation. Node genes that contain many TEs were observed in multiple potential transposable element-associated networks. This study demonstrated the increased accumulation of TEs in the introns of stress-responsive genes of pines and suggests the possibility of rewiring them into responsive networks and sub-networks interconnected with node genes containing multiple TEs. Many such regulatory influences could lead to the adaptive environmental response clines that are characteristic of naturally spread pine populations.

Adaptation to Extreme Antarctic Environments Revealed by the Genome of a Sea Ice Green Alga

The unicellular green alga Chlamydomonas sp. ICE-L thrives in polar sea ice, where it tolerates extreme low temperatures, high salinity, and broad seasonal fluctuations in light conditions. Despite the high interest in biotechnological uses of this species, little is known about the adaptations that allow it to thrive in this harsh and complex environment. Here, we assembled a high-quality genome sequence of ∼542 Mb and found that retrotransposon proliferation contributed to the relatively large genome size of ICE-L when compared to other chlorophytes. Genomic features that may support the extremophilic lifestyle of this sea ice alga include massively expanded gene families involved in unsaturated fatty acid biosynthesis, DNA repair, photoprotection, ionic homeostasis, osmotic homeostasis, and reactive oxygen species detoxification. The acquisition of multiple ice binding proteins through putative horizontal gene transfer likely contributed to the origin of the psychrophilic lifestyle in ICE-L. Additional innovations include the significant upregulation under abiotic stress of several expanded ICE-L gene families, likely reflecting adaptive changes among diverse metabolic processes. Our analyses of the genome, transcriptome, and functional assays advance general understanding of the Antarctic green algae and offer potential explanations for how green plants adapt to extreme environments.

The reference genome of the selfing fish Kryptolebias hermaphroditus: Identification of phases I and II detoxification genes

The selfing fish Kryptolebias hermaphroditus has unique reproductive system for self-fertilization, making genetically homozygous offsprings. Here, we report on high density genetic map-based genome assembly for the K. hermaphroditus Panama line (PanRS). The numbers of scaffolds were 5212 and the genome was 683,992,224 bp (N50 = 27.45 Mb). The length of anchored scaffold onto 24 linkage groups was 652,231,070 bp (95.3% of genome) with 0.01% of the gap and 39.33% of GC content and complete Benchmarking Universal Single-Copy Orthologs value was 96.6%. The numbers of annotated genes were 36,756 (average gene length 1368 bp) with the GC content of 54.1%. To examine the difference between the two sister species in the genus Kryptolebias, we compared the genomes of K. hermaphroditus PanRS and Kryptolebias marmoratus PAN line on the composition of transposable elements. To demonstrate applications of genome library, phase I and II detoxification related gene families have been analyzed, and compared the syntenies containing loci of CYP and GST genes on linkage groups. This K. hermaphroditus genome information will be helpful for a better understanding on genome-wide mechanistic view of detoxification and antioxidant-related genes over evolution in the view of fish environmental ecotoxicology.

Inter-retrotransposon amplified polymorphism markers revealed long terminal repeat retrotransposon insertion polymorphism in flax cultivated on the experimental fields around Chernobyl

Ionizing radiation in environment comes from various natural and anthropogenic sources. The effect of radioactivity released after the CNPP (Chernobyl Nuclear Power Plant) on plant systems remains of great interest. Even now, more than three decades after the nuclear accident, the long-lived radionuclides represent a strong stress factor. Herein, the emphasis has been placed on analysis of genetic variability represented by activation of LTR (Long Terminal Repeat)-retrotransposons. Polymorphism in LTR-retrotransposon insertions has been investigated throughout the genome of two flax varieties, Kyivskyi and Bethune. For this purpose, two retrotransposon-based marker techniques, IRAP (Inter-Retrotransposon Amplified Polymorphism) and iPBS (inter-Primer Binding Site), have been employed. The hypothesis that chronic radioactive stress may induce mechanism of retransposition has been supported by the activation of FL9, FL11 and FL12 LTR-retrotransposons in flax seeds harvested from radioactive environment. Out of two retrotransposon-based approaches, IRAP appears to be more suitable for identification of LTR-retrotransposon polymorphism. Even though the LTR-retrotransposon polymorphism was identified, the results suggest the high level of plant adaptation in the radioactive Chernobyl area. However, it is not really surprising that plants developed an effective strategy to survive in radio-contaminated environment over the past 30 years.

Genetic diversity of Alternaria species associated with black point in wheat grains

The genus Alternaria is a widely distributed major plant pathogen that can act as a saprophyte in plant debris. Fungi of this genus frequently infect cereal crops and cause such diseases as black point and wheat leaf blight, which decrease the yield and quality of cereal products. A total of 25 Alternaria sp. isolates were collected from germ grains of various wheat cultivars from different geographic regions in Kazakhstan. We investigated the genetic relationships of the main Alternaria species related to black point disease of wheat in Kazakhstan, using the inter-primer binding site (iPBS) DNA profiling technique. We used 25 retrotransposon-based iPBS primers to identify the differences among and within Alternaria species populations, and analyzed the variation using clustering (UPGMA) and statistical approaches (AMOVA). Isolates of Alternaria species clustered into two main genetic groups, with species of A.alternata and A.tennuissima forming one cluster, and isolates of A. infectoria forming another. The genetic diversity found using retrotransposon profiles was strongly correlated with geographic data. Overall, the iPBS fingerprinting technique is highly informative and useful for the evaluation of genetic diversity and relationships of Alternaria species.

Screening diversity and distribution of Copia retrotransposons reveals a specific amplification of BARE1 elements in genomes of the polyploid Hordeum murinum complex

We explored diversity, distribution and evolutionary dynamics of Ty1-Copia retrotransposons in the genomes of the Hordeum murinum polyploid complex and related taxa. Phylogenetic and fluorescent in situ hybridization (FISH) analyses of reverse transcriptase sequences identified four Copia families in these genomes: the predominant BARE1 (including three groups or subfamilies, A, B and C), and the less represented RIRE1, IKYA and TAR-1. Within the BARE1 family, BARE1-A elements and a subgroup of BARE1-B elements (named B1) have proliferated in the allopolyploid members of the H. murinum complex (H. murinum and H. leporinum), and in their extant diploid progenitor, subsp. glaucum. Moreover, we found a specific amplification of BARE1-B elements within each Hordeum species surveyed. The low occurrence of RIRE1, IKYA and TAR-1 elements in the allopolyploid cytotypes suggests that they are either weakly represented or highly degenerated in their diploid progenitors. The results demonstrate that BARE1-A and BARE1-B1 Copia elements are particularly well represented in the genomes of the H. murinum complex and constitute its genomic hallmark. No BARE1-A and -B1 homologs were detected in the reference barley genome. The similar distribution of RT-Copia probes across chromosomes of diploid, tetraploid and hexaploid taxa of the murinum complex shows no evidence of proliferation following polyploidization.

Long Tandem Arrays of Cassandra Retroelements and Their Role in Genome Dynamics in Plants

Retrotransposable elements are widely distributed and diverse in eukaryotes. Their copy number increases through reverse-transcription-mediated propagation, while they can be lost through recombinational processes, generating genomic rearrangements. We previously identified extensive structurally uniform retrotransposon groups in which no member contains the gag, pol, or env internal domains. Because of the lack of protein-coding capacity, these groups are non-autonomous in replication, even if transcriptionally active. The Cassandra element belongs to the non-autonomous group called terminal-repeat retrotransposons in miniature (TRIM). It carries 5S RNA sequences with conserved RNA polymerase (pol) III promoters and terminators in its long terminal repeats (LTRs). Here, we identified multiple extended tandem arrays of Cassandra retrotransposons within different plant species, including ferns. At least 12 copies of repeated LTRs (as the tandem unit) and internal domain (as a spacer), giving a pattern that resembles the cellular 5S rRNA genes, were identified. A cytogenetic analysis revealed the specific chromosomal pattern of the Cassandra retrotransposon with prominent clustering at and around 5S rDNA loci. The secondary structure of the Cassandra retroelement RNA is predicted to form super-loops, in which the two LTRs are complementary to each other and can initiate local recombination, leading to the tandem arrays of Cassandra elements. The array structures are conserved for Cassandra retroelements of different species. We speculate that recombination events similar to those of 5S rRNA genes may explain the wide variation in Cassandra copy number. Likewise, the organization of 5S rRNA gene sequences is very variable in flowering plants; part of what is taken for 5S gene copy variation may be variation in Cassandra number. The role of the Cassandra 5S sequences remains to be established.

Genome Size Variation and Comparative Genomics Reveal Intraspecific Diversity in Brassica rapa

Traditionally, reference genomes in crop species rely on the assembly of one accession, thus occulting most of intraspecific diversity. However, rearrangements, gene duplications, and transposable element content may have a large impact on the genomic structure, which could generate new phenotypic traits. Comparing two Brassica rapa genomes recently sequenced and assembled using long-read technology and optical mapping, we investigated structural variants and repetitive content between the two accessions and genome size variation among a core collection. We explored the structural consequences of the presence of large repeated sequences in B. rapa 'Z1' genome vs. the B. rapa 'Chiifu' genome, using comparative genomics and cytogenetic approaches. First, we showed that large genomic variants on chromosomes A05, A06, A09, and A10 are due to large insertions and inversions when comparing B. rapa 'Z1' and B. rapa 'Chiifu' at the origin of important length differences in some chromosomes. For instance, lengths of 'Z1' and 'Chiifu' A06 chromosomes were estimated in silico to be 55 and 29 Mb, respectively. To validate these observations, we compared using fluorescent in situ hybridization (FISH) the two A06 chromosomes present in an F1 hybrid produced by crossing these two varieties. We confirmed a length difference of 17.6% between the A06 chromosomes of 'Z1' compared to 'Chiifu.' Alternatively, using a copy number variation approach, we were able to quantify the presence of a higher number of rDNA and gypsy elements in 'Z1' genome compared to 'Chiifu' on different chromosomes including A06. Using flow cytometry, the total genome size of 12 Brassica accessions corresponding to a B. rapa available core collection was estimated and revealed a genome size variation of up to 16% between these accessions as well as some shared inversions. This study revealed the contribution of long-read sequencing of new accessions belonging to different cultigroups of B. rapa and highlighted the potential impact of differential insertion of repeat elements and inversions of large genomic regions in genome size intraspecific variability.

Early Diverging and Core Bromelioideae (Bromeliaceae) Reveal Contrasting Patterns of Genome Size Evolution and Polyploidy

The subfamily Bromelioideae is one of the most diverse groups among the neotropical Bromeliaceae. Previously, key innovations have been identified which account for the extraordinary radiation and species richness of this subfamily, especially in the so-called core Bromelioideae. However, in order to extend our understanding of the evolutionary mechanisms, the genomic mechanisms (e.g. polyploidy, dysploidy) that potentially underlie this accelerated speciation also need to be tested. Here, using PI and DAPI staining and flow cytometry we estimated genome size and GC content of 231 plants covering 30 genera and 165 species and combined it with published data. The evolutionary and ecological significance of all three genomic characters was tested within a previously generated dated phylogenetic framework using ancestral state reconstructions, comparative phylogenetic methods, and multiple regressions with climatic variables. The absolute genome size (2C) of Bromelioideae varied between 0.59 and 4.11 pg, and the GC content ranged between 36.73 and 41.43%. The monoploid genome sizes (Cx) differed significantly between core and early diverging lineages. The occurrence of dysploidy and polyploidy was, with few exceptions, limited to the phylogenetically isolated early diverging tank-less lineages. For Cx and GC content Ornstein-Uhlenbeck models outperformed the Brownian motion models suggesting adaptive potential linked to the temperature conditions. 2C-values revealed different rates of evolution in core and early diverging lineages also related to climatic conditions. Our results suggest that polyploidy is not associated with higher net diversification and fast radiation in core bromelioids. On the other hand, although coupled with higher extinction rates, dysploidy, polyploidy, and resulting genomic reorganizations might have played a role in the survival of the early diverging bromelioids in hot and arid environments.

Developing Transposable Element Marker System for Molecular Breeding

Transposable element (TE) marker system was developed considering the useful properties of the transposable elements such as their large number in the animal and plant genomes, high rate of insertion polymorphism, and ease of detection. Various methods have been employed for developing a large number of TE markers in several crop plants for genomics studies. Here we describe some of these methods including the recent whole genome search. We also review the application of TE markers in molecular breeding.

Copy Number Variation of Transposable Elements in Thinopyrum intermedium and Its Diploid Relative Species

Diploid and polyploid wild species of Triticeae have complex relationships, and the understanding of their evolution and speciation could help to increase the usability of them in wheat breeding as a source of genetic diversity. The diploid species Pseudoroegneria spicata (St), Thinopyrum bessarabicum (J^b), Dasypyrum villosum (V) derived from a hypothetical common ancestor are considered to be possible subgenome donors in hexaploid species Th. intermedium (J^rJ^vsSt, where indices r, v, and s stand for the partial relation to the genomes of Secale, Dasypyrum, and Pseudoroegneria, respectively). We quantified 10 families of transposable elements (TEs) in P. spicata, Th. bessarabicum, D. villosum (per one genome), and Th. intermedium (per one average subgenome) using the quantitative real time PCR assay and compared their abundance within the studied genomes as well as between them. Sabrina was the most abundant among all studied elements in P. spicata, D. villosum, and Th. intermedium, and among Ty3/Gypsy elements in all studied species. Among Ty1/Copia elements, Angela-A and WIS-A showed the highest and close abundance with the exception of D. villosum, and comprised the majority of all studied elements in Th. bessarabicum. Sabrina, BAGY2, and Angela-A showed similar abundance among diploids and in Th. intermedium hexaploid; Latidu and Barbara demonstrated sharp differences between diploid genomes. The relationships between genomes of Triticeae species based on the studied TE abundance and the role of TEs in speciation and polyploidization in the light of the current phylogenetic models is discussed.

TERAD: Extraction of transposable element composition from RADseq data

Transposable elements (TEs) - selfish DNA sequences that can move within the genome - comprise a large proportion of the genomes of many organisms. Although low-coverage whole-genome sequencing can be used to survey TE composition, it is noneconomical for species with large quantities of DNA. Here, we utilize restriction-site associated DNA sequencing (RADSeq) as an alternative method to survey TE composition. First, we demonstrate in silico that double digest restriction-site associated DNA sequencing (ddRADseq) markers contain the same TE compositions as whole genome assemblies across arthropods. Next, we show empirically using eight Synalpheus snapping shrimp species with large genomes that TE compositions from ddRADseq and low-coverage whole-genome sequencing are comparable within and across species. Finally, we develop a new bioinformatic pipeline, TERAD, to extract TE compositions from RADseq data. Our study expands the utility of RADseq to study the repeatome, making comparative studies of genome structure for species with large genomes more tractable and affordable.

Genome-wide variation in nucleotides and retrotransposons in alpine populations of Arabis alpina (Brassicaceae)

Advances in high-throughput sequencing have promoted the collection of reference genomes and genome-wide diversity. However, the assessment of genomic variation among populations has hitherto mainly been surveyed through single-nucleotide polymorphisms (SNPs) and largely ignored the often major fraction of genomes represented by transposable elements (TEs). Despite accumulating evidence supporting the evolutionary significance of TEs, comprehensive surveys remain scarce. Here, we sequenced the full genomes of 304 individuals of Arabis alpina sampled from four nearby natural populations to genotype SNPs as well as polymorphic long terminal repeat retrotransposons (polymorphic TEs; i.e., presence/absence of TE insertions at specific loci). We identified 291,396 SNPs and 20,548 polymorphic TEs, comparing their contributions to genomic diversity and divergence across populations. Few SNPs were shared among populations and overall showed high population-specific variation, whereas most polymorphic TEs segregated among populations. The genomic context of these two classes of variants further highlighted candidate adaptive loci having a putative impact on functional genes. In particular, 4.96% of the SNPs were identified as nonsynonymous or affecting start/stop codons. In contrast, 43% of the polymorphic TEs were present next to Arabis genes enriched in functional categories related to the regulation of reproduction and responses to biotic as well as abiotic stresses. This unprecedented data set, mapping variation gained from SNPs and complementary polymorphic TEs within and among populations, will serve as a rich resource for addressing microevolutionary processes shaping genome variation.

On the Population Dynamics of Junk: A Review on the Population Genomics of Transposable Elements

Transposable elements (TEs) play an important role in shaping genomic organization and structure, and may cause dramatic changes in phenotypes. Despite the genetic load they may impose on their host and their importance in microevolutionary processes such as adaptation and speciation, the number of population genetics studies focused on TEs has been rather limited so far compared to single nucleotide polymorphisms (SNPs). Here, we review the current knowledge about the dynamics of transposable elements at recent evolutionary time scales, and discuss the mechanisms that condition their abundance and frequency. We first discuss non-adaptive mechanisms such as purifying selection and the variable rates of transposition and elimination, and then focus on positive and balancing selection, to finally conclude on the potential role of TEs in causing genomic incompatibilities and eventually speciation. We also suggest possible ways to better model TEs dynamics in a population genomics context by incorporating recent advances in TEs into the rich information provided by SNPs about the demography, selection, and intrinsic properties of genomes.

In silico mining and FISH mapping of a chromosome-specific satellite DNA in Capsicum annuum L

BACKGROUND

A large proportion of eukaryote nuclear genomes is composed of repetitive DNA. Tracing the dynamics of repetitive elements in the genomes of related taxa can reveal important information about their phylogenic relationships as well as traits that have become distinct to a lineage.

OBJECTIVE

Study the genomic abundance and chromosomal location of repetitive DNA in Capsicum annuum L. to understand the repeat dynamics.

METHOD

We quantified repeated DNA content in the C. annuum genome using the RepeatExplorer pipeline.

RESULTS

About 42% of the C. annuum genome dataset comprised repetitive elements. Of these, 0.011, 0.98, 3.09, and 0.024% represented high and low confidence satellite repeats, putative long-terminal repeats (LTRs), and rDNA sequences, respectively. One novel high confidence 167-bp satellite repeat with a genomic proportion of 0.011%, Ca167TR, was identified. Furthermore, FISH with Ca167TR on metaphase chromosomes of C. annuum revealed signals in the subtelomeric regions of the short and long arms of chromosome 3 and 4, respectively.

CONCLUSION

Further understanding of the origin and associated functions of Ca167TR and other repeats in Capsicum will give us insights into the genomic relationships and functions of the genome.

Retrotransposon-based genetic variation of Poa annua populations from contrasting climate conditions

Background

Poa annua L. is an example of a plant characterized by abundant, worldwide distribution from polar to equatorial regions. Due to its high plasticity and extraordinary expansiveness, P. annua is considered an invasive species capable of occupying and surviving in a wide range of habitats including pioneer zones, areas intensively transformed by human activities, remote subarctic meadows and even the Antarctic Peninsula region.

Methods

In the present study, we evaluated the utility of inter-primer binding site (iPBS) markers for assessing the genetic variation of P. annua populations representing contrasting environments from the worldwide range of this species. The electrophoretic patterns of polymerase chain reaction products obtained for each individual were used to estimate the genetic diversity and differentiation between populations.

Results

iPBS genotyping revealed a pattern of genetic variation differentiating the six studied P. annua populations characterized by their different climatic conditions. According to the analysis of molecular variance, the greatest genetic variation was recorded among populations, whereas 41.75% was observed between individuals within populations. The results of principal coordinates analysis (PCoA) and model-based clustering analysis showed a clear subdivision of analyzed populations. According to PCoA, populations from Siberia and the Kola Peninsula were the most different from each other and showed the lowest genetic variability. The application of STRUCTURE software confirmed the unique character of the population from the Kola Peninsula.

Discussion

The lowest variability of the Siberia population suggested that it was subjected to genetic drift. However, although demographic expansion was indicated by negative values of Fu’s F_S statistic and analysis of mismatch distribution, it was not followed by significant traces of a bottleneck or a founder effect. For the Antarctic population, the observed level of genetic variation was surprisingly high, despite the observed significant traces of bottleneck/founder effect following demographic expansion, and was similar to that observed in populations from Poland and the Balkans. For the Antarctic population, the multiple introduction events from different sources are considered to be responsible for such an observation. Moreover, the results of STRUCTURE and PCoA showed that the P. annua from Antarctica has the highest genetic similarity to populations from Europe.

Conclusions

The observed polymorphism should be considered as a consequence of the joint influence of external abiotic stress and the selection process. Environmental changes, due to their ability to induce transposon activation, lead to the acceleration of evolutionary processes through the production of genetic variability.

Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude

Crucihimalaya himalaica, a close relative of Arabidopsis and Capsella, grows on the Qinghai-Tibet Plateau (QTP) about 4,000 m above sea level and represents an attractive model system for studying speciation and ecological adaptation in extreme environments. We assembled a draft genome sequence of 234.72 Mb encoding 27,019 genes and investigated its origin and adaptive evolutionary mechanisms. Phylogenomic analyses based on 4,586 single-copy genes revealed that C. himalaica is most closely related to Capsella (estimated divergence 8.8 to 12.2 Mya), whereas both species form a sister clade to Arabidopsis thaliana and Arabidopsis lyrata, from which they diverged between 12.7 and 17.2 Mya. LTR retrotransposons in C. himalaica proliferated shortly after the dramatic uplift and climatic change of the Himalayas from the Late Pliocene to Pleistocene. Compared with closely related species, C. himalaica showed significant contraction and pseudogenization in gene families associated with disease resistance and also significant expansion in gene families associated with ubiquitin-mediated proteolysis and DNA repair. We identified hundreds of genes involved in DNA repair, ubiquitin-mediated proteolysis, and reproductive processes with signs of positive selection. Gene families showing dramatic changes in size and genes showing signs of positive selection are likely candidates for C. himalaica's adaptation to intense radiation, low temperature, and pathogen-depauperate environments in the QTP. Loss of function at the S-locus, the reason for the transition to self-fertilization of C. himalaica, might have enabled its QTP occupation. Overall, the genome sequence of C. himalaica provides insights into the mechanisms of plant adaptation to extreme environments.

Population-specific dynamics and selection patterns of transposable element insertions in European natural populations

Transposable elements (TEs) are ubiquitous sequences in genomes of virtually all species. While TEs have been investigated for several decades, only recently we have the opportunity to study their genome-wide population dynamics. Most of the studies so far have been restricted either to the analysis of the insertions annotated in the reference genome or to the analysis of a limited number of populations. Taking advantage of the European Drosophila population genomics consortium (DrosEU) sequencing data set, we have identified and measured the dynamics of TEs in a large sample of European Drosophila melanogaster natural populations. We showed that the mobilome landscape is population-specific and highly diverse depending on the TE family. In contrast with previous studies based on SNP variants, no geographical structure was observed for TE abundance or TE divergence in European populations. We further identified de novo individual insertions using two available programs and, as expected, most of the insertions were present at low frequencies. Nevertheless, we identified a subset of TEs present at high frequencies and located in genomic regions with a high recombination rate. These TEs are candidates for being the target of positive selection, although neutral processes should be discarded before reaching any conclusion on the type of selection acting on them. Finally, parallel patterns of association between the frequency of TE insertions and several geographical and temporal variables were found between European and North American populations, suggesting that TEs can be potentially implicated in the adaptation of populations across continents.