The slowdown of Y chromosome expansion in dioecious Silene latifolia due to DNA loss and male-specific silencing of retrotransposons

Sexy ways: approaches to studying plant sex chromosomes

Sex chromosomes have evolved in many plant species with separate sexes. Current plant research is shifting from examining the structure of sex chromosomes to exploring their functional aspects. New studies are progressively unveiling the specific genetic and epigenetic mechanisms responsible for shaping distinct sexes in plants. While the fundamental methods of molecular biology and genomics are generally employed for the analysis of sex chromosomes, it is often necessary to modify classical procedures not only to simplify and expedite analyses but sometimes to make them possible at all. In this review, we demonstrate how, at the level of structural and functional genetics, cytogenetics, and bioinformatics, it is essential to adapt established procedures for sex chromosome analysis.

Non-canonical bases differentially represented in the sex chromosomes of the dioecious plant Silene latifolia

The oxidation of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC), known as oxi-mCs, garners significant interest in plants as potential epigenetic marks. While research in mammals has established a role in cell reprogramming, carcinogenesis, and gene regulation, their functions in plants remain unclear. In rice, 5hmC has been associated with transposable elements (TEs) and heterochromatin. This study utilizes Silene latifolia, a dioecious plant with heteromorphic sex chromosomes and a genome with a large proportion of TEs, which provides a favourable environment for the study of oxi-mCs in individual sexes. Notably, we detected surprisingly high levels of oxi-mCs in S. latifolia comparable with mammals. Nuclei showed enrichment in heterochromatic regions, except for 5hmC whose signal was homogeneously distributed. Intriguingly, the same X chromosome in females displayed overall enrichment of 5hmC and 5fC compared with its counterpart. This fact is shared with 5mC, resembling dosage compensation. Co-localization showed higher correlation between 5mC and 5fC than with 5hmC, indicating no potential relationship between 5hmC and 5fC. Additionally, the promoter of several sex-linked genes and sex-biased TEs clustered in a clear sex-dependent way. Together, these findings unveil a hypothetical role for oxi-mCs in S. latifolia sex chromosome development, warranting further exploration.

Why should we study plant sex chromosomes?

Understanding plant sex chromosomes involves studying interactions between developmental and physiological genetics, genome evolution, and evolutionary ecology. We focus on areas of overlap between these. Ideas about how species with separate sexes (dioecious species, in plant terminology) can evolve are even more relevant to plants than to most animal taxa because dioecy has evolved many times from ancestral functionally hermaphroditic populations, often recently. One aim of studying plant sex chromosomes is to discover how separate males and females evolved from ancestors with no such genetic sex-determining polymorphism, and the diversity in the genetic control of maleness vs femaleness. Different systems share some interesting features, and their differences help to understand why completely sex-linked regions may evolve. In some dioecious plants, the sex-determining genome regions are physically small. In others, regions without crossing over have evolved sometimes extensive regions with properties very similar to those of the familiar animal sex chromosomes. The differences also affect the evolutionary changes possible when the environment (or pollination environment, for angiosperms) changes, as dioecy is an ecologically risky strategy for sessile organisms. Dioecious plants have repeatedly reverted to cosexuality, and hermaphroditic strains of fruit crops such as papaya and grapes are desired by plant breeders. Sex-linked regions are predicted to become enriched in genes with sex differences in expression, especially when higher expression benefits one sex function but harms the other. Such trade-offs may be important for understanding other plant developmental and physiological processes and have direct applications in plant breeding.

Comparative analysis of repetitive DNA in dysploid and non-dysploid Phaseolus beans

Structural karyotype changes result from ectopic recombination events frequently associated with repetitive DNA. Although most Phaseolus species present relatively stable karyotypes with 2n = 22 chromosomes, the karyotypes of species of the Leptostachyus group show high rates of structural rearrangements, including a nested chromosome fusion that led to the dysploid chromosome number of the group (2n = 20). We examined the roles of repetitive landscapes in the rearrangements of species of the Leptostachyus group using genome-skimming data to characterize the repeatome in a range of Phaseolus species and compared them to species of that group (P. leptostachyus and P. macvaughii). LTR retrotransposons, especially the Ty3/gypsy lineage Chromovirus, were the most abundant elements in the genomes. Differences in the abundance of Tekay, Retand, and SIRE elements between P. macvaughii and P. leptostachyus were reflected in their total amounts of Ty3/gypsy and Ty1/copia. The satellite DNA fraction was the most divergent among the species, varying both in abundance and distribution, even between P. leptostachyus and P. macvaughii. The rapid turnover of repeats in the Leptostachyus group may be associated with the several rearrangements observed.

Polygenic risk modeling of tumor stage and survival in bladder cancer

Introduction

Bladder cancer assessment with non-invasive gene expression signatures facilitates the detection of patients at risk and surveillance of their status, bypassing the discomforts given by cystoscopy. To achieve accurate cancer estimation, analysis pipelines for gene expression data (GED) may integrate a sequence of several machine learning and bio-statistical techniques to model complex characteristics of pathological patterns.

Methods

Numerical experiments tested the combination of GED preprocessing by discretization with tree ensemble embeddings and nonlinear dimensionality reductions to categorize oncological patients comprehensively. Modeling aimed to identify tumor stage and distinguish survival outcomes in two situations: complete and partial data embedding. This latter experimental condition simulates the addition of new patients to an existing model for rapid monitoring of disease progression. Machine learning procedures were employed to identify the most relevant genes involved in patient prognosis and test the performance of preprocessed GED compared to untransformed data in predicting patient conditions.

Results

Data embedding paired with dimensionality reduction produced prognostic maps with well-defined clusters of patients, suitable for medical decision support. A second experiment simulated the addition of new patients to an existing model (partial data embedding): Uniform Manifold Approximation and Projection (UMAP) methodology with uniform data discretization led to better outcomes than other analyzed pipelines. Further exploration of parameter space for UMAP and t-distributed stochastic neighbor embedding (t-SNE) underlined the importance of tuning a higher number of parameters for UMAP rather than t-SNE. Moreover, two different machine learning experiments identified a group of genes valuable for partitioning patients (gene relevance analysis) and showed the higher precision obtained by preprocessed data in predicting tumor outcomes for cancer stage and survival rate (six classes prediction).

Conclusions

The present investigation proposed new analysis pipelines for disease outcome modeling from bladder cancer-related biomarkers. Complete and partial data embedding experiments suggested that pipelines employing UMAP had a more accurate predictive ability, supporting the recent literature trends on this methodology. However, it was also found that several UMAP parameters influence experimental results, therefore deriving a recommendation for researchers to pay attention to this aspect of the UMAP technique. Machine learning procedures further demonstrated the effectiveness of the proposed preprocessing in predicting patients’ conditions and determined a sub-group of biomarkers significant for forecasting bladder cancer prognosis.

Modeling binding specificities of transcription factor pairs with random forests

Background

Transcription factors (TFs) bind regulatory DNA regions with sequence specificity, form complexes and regulate gene expression. In cooperative TF-TF binding, two transcription factors bind onto a shared DNA binding site as a pair. Previous work has demonstrated pairwise TF-TF-DNA interactions with position weight matrices (PWMs), which may however not sufficiently take into account the complexity and flexibility of pairwise binding.

Results

We propose two random forest (RF) methods for joint TF-TF binding site prediction: and . We train models with previously published large-scale CAP-SELEX DNA libraries, which comprise DNA sequences enriched for binding of a selected TF pair. builds a random forest with sub-sequences selected from CAP-SELEX DNA reads with previously proposed pairwise PWM. outperforms (area under receiver operating characteristics curve, AUROC, 0.75) the current state-of-the-art method i.e. orientation and spacing specific pairwise PWMs (AUROC 0.59). Thus, may be utilized to improve prediction accuracy for pre-determined binding preferences. However, pairwise TF binding is currently considered flexible; a pair may bind DNA with different orientations and amounts of dinucleotide gaps or overlap between the two motifs. Thus, we developed , which utilizes random forests by considering simultaneously multiple orientations and spacings of the two factors. Our approach outperforms (AUROC 0.78) PWMs, as well as (p<0.00195). provides an approach for predicting TF-TF binding sites without prior knowledge on pairwise binding preferences. However, more research is needed to assess eligibility for practical applications.

Conclusions

Random forest is well suited for modeling pairwise TF-TF-DNA binding specificities, and provides an improvement to pairwise binding site prediction accuracy.

Chemical genetics in Silene latifolia elucidate regulatory pathways involved in gynoecium development

Dioecious plants possess diverse sex determination systems and unique mechanisms of reproductive organ development; however, little is known about how sex-linked genes shape the expression of regulatory cascades that lead to developmental differences between sexes. In Silene latifolia, a dioecious plant with stable dimorphism in floral traits, early experiments suggested that female-regulator genes act on the factors that determine the boundaries of the flower whorls. To identify these regulators, we sequenced the transcriptome of male flowers with fully developed gynoecia, induced by rapid demethylation in the parental generation. Eight candidates were found to have a positive role in gynoecium promotion, floral organ size, and whorl boundary, and affect the expression of class B MADS-box flower genes. To complement our transcriptome analysis, we closely examined the floral organs in their native state using field emission environmental scanning electron microscopy, and examined the differences between females and androhermaphrodites in their placenta and ovule organization. Our results reveal the regulatory pathways potentially involved in sex-specific flower development in the classical model of dioecy, S. latifolia. These pathways include previously hypothesized and unknown female-regulator genes that act on the factors that determine the flower boundaries, and a negative regulator of anther development, SUPERMAN-like (SlSUP).

Do Ty3/Gypsy Transposable Elements Play Preferential Roles in Sex Chromosome Differentiation?

Transposable elements (TEs) comprise a substantial portion of eukaryotic genomes. They have the unique ability to integrate into new locations and serve as the main source of genomic novelties by mediating chromosomal rearrangements and regulating portions of functional genes. Recent studies have revealed that TEs are abundant in sex chromosomes. In this review, we propose evolutionary relationships between specific TEs, such as Ty3/Gypsy, and sex chromosomes in different lineages based on the hypothesis that these elements contributed to sex chromosome differentiation processes. We highlight how TEs can drive the dynamics of sex-determining regions via suppression recombination under a selective force to affect the organization and structural evolution of sex chromosomes. The abundance of TEs in the sex-determining regions originates from TE-poor genomic regions, suggesting a link between TE accumulation and the emergence of the sex-determining regions. TEs are generally considered to be a hallmark of chromosome degeneration. Finally, we outline recent approaches to identify TEs and study their sex-related roles and effects in the differentiation and evolution of sex chromosomes.

The Sister Chromatid Division of the Heteromorphic Sex Chromosomes in Silene Species and Their Transmissibility towards the Mitosis

Young sex chromosomes possess unique and ongoing dynamics that allow us to understand processes that have an impact on their evolution and divergence. The genus Silene includes species with evolutionarily young sex chromosomes, and two species of section Melandrium, namely Silene latifolia (24, XY) and Silene dioica (24, XY), are well-established models of sex chromosome evolution, Y chromosome degeneration, and sex determination. In both species, the X and Y chromosomes are strongly heteromorphic and differ in the genomic composition compared to the autosomes. It is generally accepted that for proper cell division, the longest chromosomal arm must not exceed half of the average length of the spindle axis at telophase. Yet, it is not clear what are the dynamics between males and females during mitosis and how the cell compensates for the presence of the large Y chromosome in one sex. Using hydroxyurea cell synchronization and 2D/3D microscopy, we determined the position of the sex chromosomes during the mitotic cell cycle and determined the upper limit for the expansion of sex chromosome non-recombining region. Using 3D specimen preparations, we found that the velocity of the large chromosomes is compensated by the distant positioning from the central interpolar axis, confirming previous mathematical modulations.

Fundamentally different repetitive element composition of sex chromosomes in Rumex acetosa

BACKGROUND AND AIMS

Dioecious species with well-established sex chromosomes are rare in the plant kingdom. Most sex chromosomes increase in size but no comprehensive analysis of the kind of sequences that drive this expansion has been presented. Here we analyse sex chromosome structure in common sorrel (Rumex acetosa), a dioecious plant with XY1Y2 sex determination, and we provide the first chromosome-specific repeatome analysis for a plant species possessing sex chromosomes.

METHODS

We flow-sorted and separately sequenced sex chromosomes and autosomes in R. acetosa using the two-dimensional fluorescence in situ hybridization in suspension (FISHIS) method and Illumina sequencing. We identified and quantified individual repeats using RepeatExplorer, Tandem Repeat Finder and the Tandem Repeats Analysis Program. We employed fluorescence in situ hybridization (FISH) to analyse the chromosomal localization of satellites and transposons.

KEY RESULTS

We identified a number of novel satellites, which have, in a fashion similar to previously known satellites, significantly expanded on the Y chromosome but not as much on the X or on autosomes. Additionally, the size increase of Y chromosomes is caused by non-long terminal repeat (LTR) and LTR retrotransposons, while only the latter contribute to the enlargement of the X chromosome. However, the X chromosome is populated by different LTR retrotransposon lineages than those on Y chromosomes.

CONCLUSIONS

The X and Y chromosomes have significantly diverged in terms of repeat composition. The lack of recombination probably contributed to the expansion of diverse satellites and microsatellites and faster fixation of newly inserted transposable elements (TEs) on the Y chromosomes. In addition, the X and Y chromosomes, despite similar total counts of TEs, differ significantly in the representation of individual TE lineages, which indicates that transposons proliferate preferentially in either the paternal or the maternal lineage.

Histone post-translational modifications in Silene latifolia X and Y chromosomes suggest a mammal-like dosage compensation system

Silene latifolia is a model organism to study evolutionary young heteromorphic sex chromosome evolution in plants. Previous research indicates a Y-allele gene degeneration and a dosage compensation system already operating. Here, we propose an epigenetic approach based on analysis of several histone post-translational modifications (PTMs) to find the first epigenetic hints of the X:Y sex chromosome system regulation in S. latifolia. Through chromatin immunoprecipitation we interrogated six genes from X and Y alleles. Several histone PTMS linked to DNA methylation and transcriptional repression (H3K27me3, H3K23me, H3K9me2 and H3K9me3) and to transcriptional activation (H3K4me3 and H4K5, 8, 12, 16ac) were used. DNA enrichment (Immunoprecipitated DNA/input DNA) was analyzed and showed three main results: (i) promoters of the Y allele are associated with heterochromatin marks, (ii) promoters of the X allele in males are associated with activation of transcription marks and finally, (iii) promoters of X alleles in females are associated with active and repressive marks. Our finding indicates a transcription activation of X allele and transcription repression of Y allele in males. In females we found a possible differential regulation (up X₁, down X₂) of each female X allele. These results agree with the mammal-like epigenetic dosage compensation regulation.

Epigenetic aspects of floral homeotic genes in relation to sexual dimorphism in the dioecious plant Mercurialis annua

In plants, dioecy characterizes species that carry male and female flowers on separate plants and it occurs in about 6% of angiosperms; however, the molecular mechanisms that underlie dioecy are essentially unknown. The ability for sex-reversal by hormone application raises the hypothesis that the genes required for the expression of both sexes are potentially functional but are regulated by epigenetic means. In this study, proteomic analysis of nuclear proteins isolated from flower buds of females, males, and feminized males of the dioecious plant Mercurialis annua revealed differential expression of nuclear proteins that are implicated in chromatin structure and function, including floral homeotic proteins. Focusing on floral genes, we found that class B genes were mainly expressed in male flowers, while class D genes, as well as SUPERMAN-like genes, were mainly expressed in female flowers. Cytokinin-induced feminization of male plants was associated with down-regulation of male-specific genes concomitantly with up-regulation of female-specific genes. No correlation was found between the expression of class B and D genes and the changes in DNA methylation or chromatin conformation of these genes. Thus, we could not confirm DNA methylation or chromatin conformation of floral genes to be the major determinant regulating sexual dimorphisms. Instead, determination of sex in M. annua might be controlled upstream of floral genes by one or more sex-specific factors that affect hormonal homeostasis. A comprehensive model is proposed for sex-determination in M. annua.

Retrotransposons in Plant Genomes: Structure, Identification, and Classification through Bioinformatics and Machine Learning

Transposable elements (TEs) are genomic units able to move within the genome of virtually all organisms. Due to their natural repetitive numbers and their high structural diversity, the identification and classification of TEs remain a challenge in sequenced genomes. Although TEs were initially regarded as "junk DNA", it has been demonstrated that they play key roles in chromosome structures, gene expression, and regulation, as well as adaptation and evolution. A highly reliable annotation of these elements is, therefore, crucial to better understand genome functions and their evolution. To date, much bioinformatics software has been developed to address TE detection and classification processes, but many problematic aspects remain, such as the reliability, precision, and speed of the analyses. Machine learning and deep learning are algorithms that can make automatic predictions and decisions in a wide variety of scientific applications. They have been tested in bioinformatics and, more specifically for TEs, classification with encouraging results. In this review, we will discuss important aspects of TEs, such as their structure, importance in the evolution and architecture of the host, and their current classifications and nomenclatures. We will also address current methods and their limitations in identifying and classifying TEs.

The distribution of epigenetic histone marks differs between the X and Y chromosomes in Silene latifolia

Contrasting patterns of histone modifications between the X and Y chromosome in Silene latifolia show euchromatic histone mark depletion on the Y chromosome and indicate hyperactivation of one X chromosome in females. Silene latifolia (white campion) is a dioecious plant with heteromorphic sex chromosomes (24, XX in females and 24, XY in males), and a genetically degenerated Y chromosome that is 1.4 times larger than the X chromosome. Although the two sex chromosomes differ in their DNA content, information about epigenetic histone marks and evidence of their function are scarce. We performed immunolabeling experiments using antibodies specific for active and suppressive histone modifications as well as pericentromere-specific histone modifications. We show that the Y chromosome is partially depleted of histone modifications important for transcriptionally active chromatin, and carries these marks only in the pseudo-autosomal region, but that it is not enriched for suppressive and pericentromere histone marks. We also show that two of the active marks are specifically enriched in one of the X chromosomes in females and in the X chromosome in males. Our data support recent findings that genetic imprinting mediates dosage compensation of sex chromosomes in S. latifolia.

Early Sex-Chromosome Evolution in the Diploid Dioecious Plant Mercurialis annua

Suppressed recombination allows divergence between homologous sex chromosomes and the functionality of their genes. Here, we reveal patterns of the earliest stages of sex-chromosome evolution in the diploid dioecious herb Mercurialis annua on the basis of cytological analysis, de novo genome assembly and annotation, genetic mapping, exome resequencing of natural populations, and transcriptome analysis. The genome assembly contained 34,105 expressed genes, of which 10,076 were assigned to linkage groups. Genetic mapping and exome resequencing of individuals across the species range both identified the largest linkage group, LG1, as the sex chromosome. Although the sex chromosomes of M. annua are karyotypically homomorphic, we estimate that about one-third of the Y chromosome, containing 568 transcripts and spanning 22.3 cM in the corresponding female map, has ceased recombining. Nevertheless, we found limited evidence for Y-chromosome degeneration in terms of gene loss and pseudogenization, and most X- and Y-linked genes appear to have diverged in the period subsequent to speciation between M. annua and its sister species M. huetii, which shares the same sex-determining region. Taken together, our results suggest that the M. annua Y chromosome has at least two evolutionary strata: a small old stratum shared with M. huetii, and a more recent larger stratum that is probably unique to M. annua and that stopped recombining ∼1 MYA. Patterns of gene expression within the nonrecombining region are consistent with the idea that sexually antagonistic selection may have played a role in favoring suppressed recombination.

Degenerative Expansion of a Young Supergene

Long-term suppression of recombination ultimately leads to gene loss, as demonstrated by the depauperate Y and W chromosomes of long-established pairs of XY and ZW chromosomes. The young social supergene of the Solenopsis invicta red fire ant provides a powerful system to examine the effects of suppressed recombination over a shorter timescale. The two variants of this supergene are carried by a pair of heteromorphic chromosomes, referred to as the social B and social b (SB and Sb) chromosomes. The Sb variant of this supergene changes colony social organization and has an inheritance pattern similar to a Y or W chromosome because it is unable to recombine. We used high-resolution optical mapping, k-mer distribution analysis, and quantification of repetitive elements on haploid ants carrying alternate variants of this young supergene region. We find that instead of shrinking, the Sb variant of the supergene has increased in length by more than 30%. Surprisingly, only a portion of this length increase is due to consistent increases in the frequency of particular classes of repetitive elements. Instead, haplotypes of this supergene variant differ dramatically in the amounts of other repetitive elements, indicating that the accumulation of repetitive elements is a heterogeneous and dynamic process. This is the first comprehensive demonstration of degenerative expansion in an animal and shows that it occurs through nonlinear processes during the early evolution of a region of suppressed recombination.

Sex and the flower - developmental aspects of sex chromosome evolution

Background

The evolution of dioecious plants is occasionally accompanied by the establishment of sex chromosomes: both XY and ZW systems have been found in plants. Structural studies of sex chromosomes are now being followed up by functional studies that are gradually shedding light on the specific genetic and epigenetic processes that shape the development of separate sexes in plants.

Scope

This review describes sex determination diversity in plants and the genetic background of dioecy, summarizes recent progress in the investigation of both classical and emerging model dioecious plants and discusses novel findings. The advantages of interspecies hybrids in studies focused on sex determination and the role of epigenetic processes in sexual development are also overviewed.

Conclusions

We integrate the genic, genomic and epigenetic levels of sex determination and stress the impact of sex chromosome evolution on structural and functional aspects of plant sexual development. We also discuss the impact of dioecy and sex chromosomes on genome structure and expression.