Non-recombining Region Research Articles

One-step leaping evolution from an autosomal pair to the heteromorphic sex chromosomes.

Sex chromosomes evolve from an autosomal pair after the acquisition of a sex-determining gene. The primary sex chromosomes are homomorphic in both sexes and often undergo heteromorphism in either sex (XY in males or ZW in females) in association with chromosome rearrangements such as inversion, which creates a non-recombining region, called a stratum. Then, multiple strata may form by sequential inversions and extend the non-recombining region, where gene divergence accelerates, and degeneration of the Y or W chromosome progressively occurs. Contrast to the conventional theory, we propose a shortcut in heteromorphic sex chromosome evolution, where an autosomal pair directly evolves into a heteromorphic sex chromosome pair. We illustrate this with two frog cases where Y chromosome or autosome, which is morphologically inverted, was introgressed from another species through interspecific hybridization, instantly forming a new heteromorphic sex chromosome pair. This event resulted in a distinct non-recombining region immediately after hybridization. The introduction of an inverted chromosome from a different species may be associated with benefits in morphology, breeding behavior, hybrid viability, sex determination, and recovery of the sex ratio of the hybrids. We discuss the molecular mechanisms driving preferential mutations in the introduced, inverted chromosome through interspecific hybridization.

Delineating the W sex chromosome in the clam shrimp, Eulimnadia texana

Introduction: Sex chromosomes have evolved independently across various lineages, often showing convergent degradation of the sex-limited chromosome. While extensively studied in model organisms with ancient sex chromosomal systems, the evolution of early-stage sex chromosomes remains poorly understood. Eulimnadia texana, a freshwater crustacean with a unique androdioecious breeding system (ZZ, ZW, and viable WW genotypes), provides a rare opportunity to study early sex chromosome evolution. This study examines E. texana's W chromosome for evidence of a small localized non-recombining region, characterized by a transposable element (TE) "hotspot", low gene density, and low GC content. Methods: Sex-linked markers were mapped onto the W chromosome (Scaffold 1). TEs in the WW genome were identified using RepeatModeler and RepeatMasker. Statistical analyses compared TE distribution between the genome and Scaffold 1, which was then divided into 20 equal-sized "bins" for finer-scale statistical analyses. Gene density and GC content were analyzed across these bins. Results: While no significant TE accumulation was found across the entire W chromosome compared to the remaining genome, a specific region (6.6-8.8 Mb, fourth bin) showed significantly higher TE accumulation. This region also exhibited low gene density and low GC content, indicative of reduced recombination. Conclusion: Our findings suggest that E. texana's W chromosome contains a smaller region of crossover suppression, supporting the hypothesis that it is in a proto-sex chromosome in early evolutionary development. This study provides valuable insights into early sex chromosome evolution and establishes E. texana as an ideal model for further investigation of evolutionary processes driving proto-sex chromosome differentiation.

A million-cow genome-wide association study of productive life in U.S. Holstein cows

BackgroundProductive life (PL) of a cow is the time the cow remains in the milking herd from first calving to exit from the herd due to culling or death and is an important economic trait in U.S. Holstein cattle. The large samples of Holstein genomic evaluation data that have become available recently provided unprecedented statistical power to identify genetic factors affecting PL in Holstein cows using the approach of genome-wide association study (GWAS).MethodsThe GWAS analysis used 1,103,641 Holstein cows with phenotypic observations on PL and genotypes of 75,282 single nucleotide polymorphism (SNP) markers. The statistical tests and estimation of SNP additive and dominance effects used the approximate generalized least squares method implemented by the EPISNPmpi computer program.ResultsThe GWAS detected 5390 significant additive effects of PL distributed over all 29 autosomes and the X–Y nonrecombining region of the X chromosome (Chr31). Two chromosome regions had the most significant and largest cluster of additive effects, the SLC4A4-GC-NPFFR2 (SGN) region of Chr06 with pleiotropic effects for PL, fertility, somatic cell score and milk yield; and the 32–52 Mb region of Chr10 with peak effects for PL in or near RASGRP1 with many important immunity functions. The dominance tests detected 38 significant dominance effects including 12 dominance effects with sharply negative homozygous recessive genotypes on Chr18, Chr05, Chr23 and Chr24.ConclusionsThe GWAS results showed that highly significant genetic effects for PL were in chromosome regions known to have highly significant effects for fertility and health and a chromosome region with multiple genes with reproductive and immunity functions. SNPs with rare but sharply negative homozygous recessive genotypes for PL existed and should be used for eliminating heifers carrying those homozygous recessive genotypes.

The genetic architecture of the load linked to dominant and recessive self-incompatibility alleles in Arabidopsis halleri and Arabidopsis lyrata

The long-term balancing selection acting on mating types or sex-determining genes is expected to lead to the accumulation of deleterious mutations in the tightly linked chromosomal segments that are locally ‘sheltered’ from purifying selection. However, the factors determining the extent of this accumulation are poorly understood. Here, we took advantage of variations in the intensity of balancing selection along a dominance hierarchy formed by alleles at the sporophytic self-incompatibility system of the Brassicaceae to compare the pace at which linked deleterious mutations accumulate among them. We first experimentally measured the phenotypic manifestation of the linked load at three different levels of the dominance hierarchy. We then sequenced and phased polymorphisms in the chromosomal regions linked to 126 distinct copies of S-alleles in two populations of Arabidopsis halleri and three populations of Arabidopsis lyrata. We find that linkage to the S-locus locally distorts phylogenies over about 10–30 kb along the chromosome. The more intense balancing selection on dominant S-alleles results in greater fixation of linked deleterious mutations, while recessive S-alleles accumulate more linked deleterious mutations that are segregating. Hence, the structure rather than the overall magnitude of the linked genetic load differs between dominant and recessive S-alleles. Our results have consequences for the long-term evolution of new S-alleles, the evolution of dominance modifiers between them, and raise the question of why the non-recombining regions of some sex and mating type chromosomes expand over evolutionary times while others, such as the S-locus of the Brassicaceae, remain restricted to small chromosomal regions.

The genetic architecture of the load linked to dominant and recessive self-incompatibility alleles in Arabidopsis halleri and Arabidopsis lyrata.

The long-term balancing selection acting on mating types or sex-determining genes is expected to lead to the accumulation of deleterious mutations in the tightly linked chromosomal segments that are locally 'sheltered' from purifying selection. However, the factors determining the extent of this accumulation are poorly understood. Here, we took advantage of variations in the intensity of balancing selection along a dominance hierarchy formed by alleles at the sporophytic self-incompatibility system of the Brassicaceae to compare the pace at which linked deleterious mutations accumulate among them. We first experimentally measured the phenotypic manifestation of the linked load at three different levels of the dominance hierarchy. We then sequenced and phased polymorphisms in the chromosomal regions linked to 126 distinct copies of S-alleles in two populations of Arabidopsis halleri and three populations of Arabidopsis lyrata. We find that linkage to the S-locus locally distorts phylogenies over about 10-30 kb along the chromosome. The more intense balancing selection on dominant S-alleles results in greater fixation of linked deleterious mutations, while recessive S-alleles accumulate more linked deleterious mutations that are segregating. Hence, the structure rather than the overall magnitude of the linked genetic load differs between dominant and recessive S-alleles. Our results have consequences for the long-term evolution of new S-alleles, the evolution of dominance modifiers between them, and raise the question of why the non-recombining regions of some sex and mating type chromosomes expand over evolutionary times while others, such as the S-locus of the Brassicaceae, remain restricted to small chromosomal regions.

Discovery unbalanced DNA mixtures and evaluation mixing ratio via a droplet digital PCR method.

Small amounts of DNA from a perpetrator collected during crime-scene investigations can be masked by large amounts of DNA from the victim. These samples can provide important information for the perpetrator's conviction. Short tandem repeat (STR) detection system is not sensitive enough to detect trace amounts of minor components in unbalanced mixed DNA. We developed a system using droplet digital polymerase chain reaction (ddPCR) capable of discovering trace components and accurately determining the ratio of mixed DNA in extremely unbalanced mixtures. The non-recombining regions of the X chromosome and Y chromosome were quantified in the DNA of male and female mixtures using duplex ddPCR. Absolute quantification of low-abundance portions of trace samples and unbalanced mixtures was done using different mixing ratios. The ddPCR system could be used to detect low-abundance samples with < 5 copies of DNA components in an extremely unbalanced mixture at a mixing ratio of 10000:1. The high sensitivity and specificity of the system could identify the mixing ratio of mixed DNA accurately. A ddPCR system was developed for evaluation of mixed samples of male DNA and female DNA. Our system could detect DNA quantities as low as 5 copies in extremely unbalanced mixed samples with good specificity and applicability. This method could assist forensic investigators in avoiding the omission of important physical evidence, and evaluating the ratio of mixed male/female trace samples.

Genome-wide association studies using chromosome-scale genomes of male and female lines redefines two sex-linked loci in linkage disequilibrium in Ficus carica L.

Despite its significance in breeding and agricultural practices, the genetic mechanisms that determine sex in the fig (Ficus carica) are not clearly understood, highlighting the need for comprehensive genomic and genetic research to clarify these complex processes. With the aim to accurately identifying the genetic loci involved in sex-related traits of the fig, two chromosome-scale genome sequences for ‘Horaishi’ (a female line) and ‘Caprifig 6085′ (a male line) were determined, and a genome-wide association study (GWAS) was subsequently conducted using single nucleotide polymorphisms (SNPs) obtained from whole genome sequencing across genetically diverse lines and a biparental mapping population derived from a cross between the two lines. As a result, two major GWAS signals were detected. Of these, only the SNPs on FcRAN1 were found to perfectly associate with the sex phenotype, which was suggested as a sex determinant in our previous study. On the other hand, the SNP pattern of FcAG did not perfectly match the sex phenotype, but FcAG was located near and showed another strong GWAS signal, and in contrast, the polymorphisms in the deletion of FcAG did not contradict the sex phenotype or the SNP pattern of FcRAN1. These results suggest that the candidate genes for the sex-linked traits in figs are narrowed down to both FcRAN1 and FcAG, indicating that they are not limited to one or the other. Furthermore, they suggest some potential scenarios, including the possibility that the sex-linked traits are controlled either complementarily or collaboratively by these two genes, which are in linkage disequilibrium. In addition, we confirmed that a non-recombining sex-linked region exists between these two gene loci. The genome information and the finding for sex-related traits would be helpful to understand sex determination mechanisms in Ficus plants as well as the fig breeding programs.

A catalogue of recombination coldspots in interspecific tomato hybrids.

Increasing natural resistance and resilience in plants is key for ensuring food security within a changing climate. Breeders improve these traits by crossing cultivars with their wild relatives and introgressing specific alleles through meiotic recombination. However, some genomic regions are devoid of recombination especially in crosses between divergent genomes, limiting the combinations of desirable alleles. Here, we used pooled-pollen sequencing to build a map of recombinant and non-recombinant regions between tomato and five wild relatives commonly used for introgressive tomato breeding. We detected hybrid-specific recombination coldspots that underscore the role of structural variations in modifying recombination patterns and maintaining genetic linkage in interspecific crosses. Crossover regions and coldspots show strong association with specific TE superfamilies exhibiting differentially accessible chromatin between somatic and meiotic cells. About two-thirds of the genome are conserved coldspots, located mostly in the pericentromeres and enriched with retrotransposons. The coldspots also harbor genes associated with agronomic traits and stress resistance, revealing undesired consequences of linkage drag and possible barriers to breeding. We presented examples of linkage drag that can potentially be resolved by pairing tomato with other wild species. Overall, this catalogue will help breeders better understand crossover localization and make informed decisions on generating new tomato varieties.

Genetic landscape of the people of Bangladesh depicted with 17 Y-Chromosome-Specific microsatellites

Seventeen microsatellite loci from the non-recombining region of the human Y chromosome were typed using AmpFlSTR® Yfiler® PCR amplification systems in 404 male subjects belonging to the three largest ethnic groups in Bangladesh. A total of 150 haplotypes from the Chakma, 144 from the Tripura, and 110 from the Khasia were detected with a corresponding discrimination capacity of 73.885%, 65.563%, and 81.250%, respectively. The highest allele frequency of 0.828 was detected in DYS391 locus in the Tripuras, while the lowest allele frequency of 0.009 was detected at the same locus for the Chakma population. The highest gene diversity (0.964) was observed at DYS385a/b locus in the Khasias, while the lowest gene diversity (0.301) was detected at DYS391 locus in Tripuras population. The overall haplotype diversity for the studied populations was 0.986141. Both the Neighbour-Joining tree and pairwise genetic distances showed that Chakma lies closer to a clade consisting of Tripuras (Khagrachari, Bangladesh) and Tripuri (Tripura, India). In contrast, the Khasias demonstrated a close affinity with the Oraon (Chhattisgarh, India), followed by the Santals. The Y-STR haplotype matching probabilities within and between populations demonstrated that the Chakma, Tripura, and Khasia were 100% genetically distinct. The studied ethnic populations exhibited higher frequency for haplogroups L and Q as opposed to haplogroups R1a, H, and L found in the mainstream Bengali population. The Median-joining networking showed haplogroups L and R1a have the most compact clustering within populations, followed by haplogroups Q and H. The presence of haplogroup R1a suggests that Bengali may have originated through west-to-east migration, whereas haplogroups L and Q distribution in the studied tribes reveal a very significant affinity with the South-East Asian populations and may have shared a common ancestral origin with the Mongoloid stock populations. Bioresearch Commu. 10(2): 1482-1488, 2024 (July)

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae

BackgroundThe Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems.ResultsWe explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19–27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary.ConclusionsTogether, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Evolutionary genomics of socially polymorphic populations of Pogonomyrmex californicus

BackgroundSocial insects vary considerably in their social organization both between and within species. In the California harvester ant, Pogonomyrmex californicus (Buckley 1867), colonies are commonly founded and headed by a single queen (haplometrosis, primary monogyny). However, in some populations in California (USA), unrelated queens cooperate not only during founding (pleometrosis) but also throughout the life of the colony (primary polygyny). The genetic architecture and evolutionary dynamics of this complex social niche polymorphism (haplometrosis vs pleometrosis) have remained unknown. ResultsWe provide a first analysis of its genomic basis and evolutionary history using population genomics comparing individuals from a haplometrotic population to those from a pleometrotic population. We discovered a recently evolved (< 200 k years), 8-Mb non-recombining region segregating with the observed social niche polymorphism. This region shares several characteristics with supergenes underlying social polymorphisms in other socially polymorphic ant species. However, we also find remarkable differences from previously described social supergenes. Particularly, four additional genomic regions not in linkage with the supergene show signatures of a selective sweep in the pleometrotic population. Within these regions, we find for example genes crucial for epigenetic regulation via histone modification (chameau) and DNA methylation (Dnmt1). ConclusionsAltogether, our results suggest that social morph in this species is a polygenic trait involving a potential young supergene. Further studies targeting haplo- and pleometrotic individuals from a single population are however required to conclusively resolve whether these genetic differences underlie the alternative social phenotypes or have emerged through genetic drift.

Identification of a common haplotype in carriers of rob(1;29) in 32 Italian cattle breeds

Robertsonian translocation 1;29 (rob(1;29)), a widespread chromosomal anomaly affecting cattle fertility, appears to have originated from a common ancestor. This study utilizes routine SNP data to investigate the chromosomal region associated with rob(1;29) and confirm the presence of a shared haplotype among carriers in diverse Italian breeds. Three datasets were employed: Dataset 1 included 151 subjects from 5 beef cattle breeds genotyped with the GGP Bovine 33 k SNP chip; Dataset 2 encompassed 800 subjects from 32 Italian breeds genotyped with the Illumina 50 k SNP chip, sourced from the BOVITA dataset; Dataset 3 combined Dataset 2 with 21 karyologically tested subjects from breeds with a high carrier frequency, genotyped using the Affymetrix 65 K SNP chip. FST analysis pinpointed a distinctive genomic region on the first six Mb of BTA29, the centromeric region involved in the translocation. Haplotype comparisons within this non-recombining region revealed a common haplotype shared among all carriers, supporting the theory of a common ancestor. Principal component and haplotype analysis allowed clear differentiation of rob(1;29) homozygous and heterozygous carriers. Expanding to Dataset 2 revealed rob(1;29) carriers in unexpected breeds, all sharing the same ancestral haplotype. Notably, previously untested breeds, including Cinisara, exhibited a high carrier prevalence (nearly 50%), confirmed by karyological analysis. This study validates the presence of a shared haplotype among all identified rob(1;29) carriers, reinforcing the common ancestor theory as the origin of this translocation's spread throughout the cattle population. Furthermore, it underscores the potential of SNP data analysis as a rapid, accurate, and cost-effective tool for broad rob(1;29) screening, given the translocation's consistent nature across all analyzed breeds.

Genomic architecture of resistance to latania scale (H. lataniae) in kiwifruit (A. chinensis var. chinensis)

BackgroundLatania scale (Hemiberlesia lataniae Signoret) is an armoured scale insect known to cause damage to kiwifruit plants and fruit, which ultimately reduces crop values and creates post-harvest export and quarantine issues. Resistance to H. lataniae does exist in some commercial cultivars of kiwifruit. However, some of the commercial cultivars bred in New Zealand have not inherited alleles for resistance to H. lataniae carried by their parents. To elucidate the architecture of resistance in the parents and develop molecular markers to assist breeding, these experiments analysed the inheritance of resistance to H. lataniae from families related to commercial cultivars.ResultsThe first experiment identified a 15.97 Mb genomic region of interest for resistance to H. lataniae in rtGBS data of 3.23 to 19.20 Mb on chromosome 10. A larger population was then QTL mapped, which confirmed the region of interest as the sole locus contributing to H. lataniae resistance. inDel markers mapping the region of low recombination under the QTL peak further narrowed the region associated with H. lataniae resistance to a 5.73 Mb region.ConclusionsThe kiwifruit populations and genomic methods used in this study identify the same non-recombinant region of chromosome 10 which confers resistance of A. chinensis var. chinensis to H. lataniae. The markers developed to target the H. lataniae resistance loci will reduce the amount of costly and time-consuming phenotyping required for breeding H. lataniae scale resistance into new kiwifruit cultivars.

Inka Child Mummy Found in Cerro Aconcagua (Argentina) Traced Back to Populations of the Northern Peruvian Coast through Y-Chromosome Analysis

The mummy of a seven-year-old child that was discovered in 1985 in Cerro Aconcagua (Mendoza, Argentina) was likely part of an Inka sacrificial religious practice known as capacocha. Previous uniparental DNA marker studies conducted by some scholars have suggested that the mummified child may be related to the southern Andean population of Peru. However, autosome genome-wide analysis performed by others has indicated that the child was more closely related to the population along the northern Peruvian coast than to that of the southern Andes. In this study, we aimed to determine possible genealogical connections in the male lineage of the mummified child. To achieve this, we compared the genetic profile of the mummy with an extensive database of contemporary individuals from the northern Peruvian coastal and southern Andean regions. We used single nucleotide polymorphisms and short tandem repeats from the nonrecombining region of the Y-chromosome for our analysis. Our results confirmed that the Inka child mummy was closely related to individuals from the north coast of Peru. This suggests that the child was likely descended from the Muchik–Chimor-speaking people.

Comprehensive analyses of genetic diversities and population structure of the Guizhou Dong group based on 44 Y-markers.

The non-recombining region of the human Y chromosome (NRY) is a strictly paternally inherited genetic marker and the best material to trace the paternal lineages of populations. Y chromosomal short tandem repeat (Y-STR) is characterized by high polymorphism and paternal inheritance pattern, so it has been widely used in forensic medicine and population genetic research. This study aims to understand the genetic distribution of Y-STRs in the Guizhou Dong population, provide reference data for forensic application, and explore the phylogenetic relationships between the Guizhou Dong population and other comparison populations. Based on the allele profile of 44 Y-markers in the Guizhou Dong group, we estimate their allele frequencies and haplotype frequencies. In addition, we also compare the forensic application efficiency of different Y-STR sets in the Guizhou Dong group. Finally, genetic relationships among Guizhou Dong and other reference populations are dissected by the multi-dimensional scaling and the phylogenetic tree. A total of 393 alleles are observed in 312 Guizhou Dong individuals for these Y-markers, with allele frequencies ranging from 0.0032 to 0.9679. The haplotype diversity and discriminatory capacity for these Y-markers in the Guizhou Dong population are 0.99984 and 0.97440, respectively. The population genetic analyses of the Guizhou Dong group and other reference populations show that the Guizhou Dong group has the closest genetic relationship with the Hunan Dong population, and followed by the Guizhou Tujia population. In conclusion, these 44 Y-markers can be used as an effective tool for male differentiation in the Guizhou Dong group. The haplotype data in this study not only enrich the Y-STR data of different ethnic groups in China, but also have important significance for population genetics and forensic research.

Dynamics of transposable element accumulation in the non-recombining regions of mating-type chromosomes in anther-smut fungi

In the absence of recombination, the number of transposable elements (TEs) increases due to less efficient selection, but the dynamics of such TE accumulations are not well characterized. Leveraging a dataset of 21 independent events of recombination cessation of different ages in mating-type chromosomes of Microbotryum fungi, we show that TEs rapidly accumulated in regions lacking recombination, but that TE content reached a plateau at ca. 50% of occupied base pairs by 1.5 million years following recombination suppression. The same TE superfamilies have expanded in independently evolved non-recombining regions, in particular rolling-circle replication elements (Helitrons). Long-terminal repeat (LTR) retrotransposons of the Copia and Ty3 superfamilies also expanded, through transposition bursts (distinguished from gene conversion based on LTR divergence), with both non-recombining regions and autosomes affected, suggesting that non-recombining regions constitute TE reservoirs. This study improves our knowledge of genome evolution by showing that TEs can accumulate through bursts, following non-linear decelerating dynamics.

The fate of sex chromosomes during the evolution of monoicy from dioicy in liverworts

Liverworts comprise one of six primary land plant lineages, with the predicted origin of extant liverwort diversity dating to the Silurian. The ancestral liverwort has been inferred to have been dioicous (unisexual) with chromosomal sex determination in which the U chromosome of females and the V chromosome of males were dimorphic with an extensive non-recombining region. In liverworts, sex is determined by a U chromosomal "feminizer" gene that promotes female development, and in its absence, male development ensues. Monoicy (bisexuality) has independently evolved multiple times within liverworts. Here, we explore the evolution of monoicy, focusing on the monoicous species Ricciocarpos natans, and propose that the evolution of monoicy in R.natans involved the appearance of an aneuploid spore that possessed both U and V chromosomes. Chromosomal rearrangements involving the U chromosome resulted in distribution of essential U chromosome genes, including the feminizer, to several autosomal locations. By contrast, we infer that the ancestral V chromosome was inherited largely intact, probably because it carries numerous dispersed "motility" genes distributed across the chromosome. The genetic networks for sex differentiation in R.natans appear largely unchanged except that the feminizer is developmentally regulated, allowing for temporally separated differentiation of female and male reproductive organs on a single plant. A survey of other monoicous liverworts suggests that similar genomic rearrangements may have occurred repeatedly in lineages transitioning to monoicy from dioicy. These data provide a foundation for understanding how genetic networks controlling sex determination can be subtly rewired to produce profound changes in sexual systems.

Ongoing Rapid Evolution of a Post-Y Region Revealed by Chromosome-Scale Genome Assembly of a Hexaploid Monoecious Persimmon (Diospyros kaki).

Plants have evolved sex chromosomes independently in many lineages, and loss of separate sexes can also occur. In this study, we assembled a monoecious recently hexaploidized persimmon (Diospyros kaki), in which the Y chromosome has lost the maleness-determining function. Comparative genomic analysis of D. kaki and its dioecious relatives uncovered the evolutionary process by which the nonfunctional Y chromosome (or Ymonoecy) was derived, which involved silencing of the sex-determining gene, OGI, approximately 2 million years ago. Analyses of the entire X and Ymonoecy chromosomes suggested that D. kaki's nonfunctional male-specific region of the Y chromosome (MSY), which we call a post-MSY, has conserved some characteristics of the original functional MSY. Specifically, comparing the functional MSY in Diospyros lotus and the nonfunctional "post-MSY" in D. kaki indicated that both have been rapidly rearranged, mainly via ongoing transposable element bursts, resembling structural changes often detected in Y-linked regions, some of which can enlarge the nonrecombining regions. The recent evolution of the post-MSY (and possibly also MSYs in dioecious Diospyros species) therefore probably reflects these regions' ancestral location in a pericentromeric region, rather than the presence of male-determining genes and/or genes controlling sexually dimorphic traits.

The origin and evolution of sex chromosomes, revealed by sequencing of the Silene latifolia female genome

White campion (Silene latifolia, Caryophyllaceae) was the first vascular plant where sex chromosomes were discovered. This species is a classic model for studies on plant sex chromosomes due to presence of large, clearly distinguishable X and Y chromosomes that originated de novo about 11 million years ago (mya), but lack of genomic resources for this relatively large genome (∼2.8 Gb) remains a significant hurdle. Here we report S.latifolia female genome assembly integrated with sex-specific genetic maps of this species, focusing on sex chromosomes and their evolution. The analysis reveals a highly heterogeneous recombination landscape with strong reduction in recombination rate in the central parts of all chromosomes. Recombination on the X chromosome in female meiosis primarily occurs at the very ends, and over 85% of the X chromosome length is located in a massive (∼330 Mb) gene-poor, rarely recombining pericentromeric region (Xpr). The results indicate that the non-recombining region on the Y chromosome (NRY) initially evolved in a relatively small (∼15 Mb), actively recombining region at the end of the q-arm, possibly as a result of inversion on the nascent X chromosome. The NRY expanded about 6 mya via linkage between the Xpr and the sex-determining region, which may have been caused by expanding pericentromeric recombination suppression on the X chromosome. These findings shed light on the origin of sex chromosomes in S.latifolia and yield genomic resources to assist ongoing and future investigations into sex chromosome evolution.

Heterochiasmy and Sex Chromosome Evolution in Silene

The evolution of a non-recombining sex-specific region is a key step in sex chromosome evolution. Suppression of recombination between the (proto-) X- and Y-chromosomes in male meiosis creates a non-recombining Y-linked region (NRY), while the X-chromosome continues to recombine in females. Lack of recombination in the NRY defines its main properties-genetic degeneration and accumulation of repetitive DNA, making X and Y chromosomes very different from each other. How and why recombination suppression on sex chromosomes evolves remains controversial. A strong difference in recombination rates between the sexes (heterochiasmy) can facilitate or even cause recombination suppression. In the extreme case-complete lack of recombination in the heterogametic sex (achiasmy)-the entire sex-specific chromosome is automatically non-recombining. In this study, I analyse sex-specific recombination rates in a dioecious plant Silene latifolia (Caryophyllaceae), which evolved separate sexes and sex chromosomes ~11 million years ago. I reconstruct high-density RNAseq-based genetic maps including over five thousand genic markers for the two sexes separately. The comparison of the male and female maps reveals only modest heterochiasmy across the genome, with the exception of the sex chromosomes, where recombination is suppressed in males. This indicates that heterochiasmy likely played only a minor, if any, role in NRY evolution in S. latifolia, as recombination suppression is specific to NRY rather than to the entire genome in males. Other mechanisms such as structural rearrangements and/or epigenetic modifications were likely involved, and comparative genome analysis and genetic mapping in multiple Silene species will help to shed light on the mechanism(s) of recombination suppression that led to the evolution of sex chromosomes.