Regulation of cyclic peptide biosynthesis and pathogenicity in Cochliobolus carbonum by TOXEp, a novel protein with a bZIP basic DNA-binding motif and four ankyrin repeats

HC-toxin is an epoxide-containing cyclic tetrapeptide that is a critical virulence determinant in the pathogenic interaction between the filamentous fungus Cochliobolus carbonum and maize. HC-toxin exerts a potent cytostatic effect on plant and animal cells by inhibiting histone deacetylase. The biosynthesis of HC-toxin by C. carbonum is controlled by a complex genetic locus, TOX2, that contains multiple, duplicated copies of genes encoding export and biosynthetic enzymes. A new gene in the TOX2 complex, TOXE, has now been isolated. Mutation of TOXE by targeted gene disruption has no effect on growth and sporulation but abolishes HC-toxin production and pathogenicity. TOXE is required for the expression of three genes with a known or putative role in HC-toxin production, but is not required for expression of HTS1, which encodes the large, multifunctional peptide synthetase that is the central enzyme in HC-toxin biosynthesis. At its N-terminus, TOXEp has a bZIP basic DNA binding domain, but it does not contain any discernible leucine zipper or helix-loop-helix. At its carboxy terminus, TOXEp contains four ankyrin repeats. In having these two common regulatory motifs in a single polypeptide, TOXEp appears to represent a novel class of regulatory protein. TOXE is present only in HC-toxin-producing (Tox2⁺) isolates of C. carbonum. Most Tox2⁺ isolates have two copies; in strain SB111, one copy of TOXE is on the same 3.5-Mb chromosome that contains all of the other genes known to be involved in HC-toxin biosynthesis, and the second copy of TOXE is on a 0.7-Mb chromosome. Received: 20 April 1998 / Accepted: 21 September 1998     

Chromosomal organization of TOX2, a complex locus controlling host-selective toxin biosynthesis in Cochliobolus carbonum.

Race 1 isolates of the filamentous fungus Cochliobolus carbonum are exceptionally virulent on certain genotypes of maize due to production of a cyclic tetrapeptide, HC-toxin. In crosses between toxin-producing (Tox2+) and toxin-nonproducing (Tox2-) isolates, toxin production segregates in a simple 1:1 pattern, suggesting the involvement of a single genetic locus, which has been named TOX2. Earlier work had shown that in isolate SB111, TOX2 consists in part of two copies of a gene, HTS1, that encodes a 570-kD cyclic peptide synthetase and is lacking in Tox2- isolates. The genomic structure of TOX2 and the relationship between the two copies of HTS1 have now been clarified by using pulsedfield gel electrophoresis and physical mapping. In isolate SB111, both copies of HTS1 are on the largest chromosome (3.5 Mb), which is not present in the related Tox2- strain SB114. Two other genes known or thought to be important for HC-toxin biosynthesis, TOXA and TOXC, are also on the same chromosome in multiple copies. Other independent Tox2+ isolates also have two linked copies of HTS1, but in some isolates the size of the chromosome containing HTS1 is 2.2 Mb. Evidence obtained with Tox2+ -unique and with random probes is consistent with a reciprocal translocation as the cause of the difference in the size of the HTS1-containing chromosome among the Tox2+ isolates studied here. Physical mapping of the 3.5-Mb chromosome of SB111 that contains HTS1 using rare-cutting restriction enzymes and engineered restriction sites was used to map the chromosome location of the two copies of HTS1 and the three copies of TOXC. The results indicate that TOX2 is a complex locus that extends over more than 500 kb. The capacity to produce HC-toxin did not evolve by any single, simple mechanism.     

Reduced virulence caused by meiotic instability of the TOX2 chromosome of the maize pathogen Cochliobolus carbonum

The mechanisms by which pathogenic fungi evolve are poorly understood. Production of the host-selective cyclic peptide HC-toxin is controlled by a complex locus, TOX2, in the plant pathogen Cochliobolus carbonum. Crosses between toxin-producing (Tox2+) and toxin-nonproducing (Tox2-) isolates, as well as crosses between isolates in which the TOX2 genes were on chromosomes of different size, yielded progeny that had lost one or more copies of one or more of the TOX2 genes. Of approximately 200 progeny analyzed, eight (4%) had lost at least one TOX2 gene. All of them still had at least one functional copy of all of the known genes required for HC-toxin production (HTS1, TOXA, TOXC, and TOXE). Most deletion strains could be explained by simple chromosome breaks resulting in the loss of major contiguous portions (0.8 to 1.4 Mb) of the 3.5-Mb TOX2 chromosome, whereas others had more complicated patterns. All deletion strains had normal growth and were fertile, indicating that the 1.4 Mb of DNA contained no essential housekeeping genes. Most strains were also still virulent (Tox2+), but two had a novel phenotype of reduced virulence (RV), characterized by smaller lesions that expanded at a reduced rate and an inability to colonize plants systemically. Although the RV strains made no detectable HC-toxin in culture, the RV phenotype was dependent on the presence of a functional copy of HTS1, which encodes the central enzyme in HC-toxin biosynthesis. We propose that the RV strains still make a low level of HC-toxin, at least in planta, and that this is due to the loss of one or more genes that contribute to, but are not absolutely required for, HC-toxin synthesis.     

An extended physical map of the TOX2 locus of Cochliobolus carbonum required for biosynthesis of HC-toxin

In genetic crosses, HC-toxin production in the filamentous fungus Cochliobolus carbonum appears to be controlled by a single locus, TOX2. At the molecular level, TOX2 is composed of at least seven duplicated and coregulated genes involved in HC-toxin biosynthesis, export, and regulation. All copies of four of the TOX2 genes were previously mapped within a 540-kb stretch of DNA in strain SB111. Subsequently, an additional three TOX2 genes, TOXE, TOXF, and TOXG, have been discovered. In this paper we have mapped all copies of the new genes, a total of seven, and show that except for one of the two copies of TOXE, which was previously shown to be on a chromosome of 0.7 Mb in strain SB111, they are all linked to the previously known TOX2 genes within approximately 600 kb of each other on a chromosome of 3.5 Mb. We show here that this chromosome also contains at least one non-TOX2 gene, EXG2, which encodes an exo-beta1,3-glucanase. EXG2 is still present in strains that have undergone spontaneous deletion of up to approximately 1.4 Mb of the 3.5-Mb chromosome. The results contribute to our understanding of the complex organization of the genes involved in HC-toxin biosynthesis and are consistent with the hypothesis that a reciprocal chromosomal translocation accounts for the pattern of distribution of the TOX2 genes in different C. carbonum isolates.     

A eukaryotic alanine racemase gene involved in cyclic peptide biosynthesis

The cyclic tetrapeptide HC-toxin is an essential virulence determinant for the plant pathogenic fungus Cochliobolus carbonum and an inhibitor of histone deacetylase. The major form of HC-toxin contains the D-isomers of Ala and Pro. The non-ribosomal peptide synthetase that synthesizes HC-toxin has only one epimerizing domain for conversion of L-Pro to D-Pro; the source of D-Ala has remained unknown. Here we present the cloning and characterization of a new gene involved in HC-toxin biosynthesis, TOXG. TOXG is present only in HC-toxin-producing (Tox2(+)) isolates of C. carbonum. TOXG is able to support D-Ala-independent growth of a strain of Escherichia coli defective in D-Ala synthesis. A C. carbonum strain with both of its copies of TOXG mutated grows normally in culture, and although it no longer makes the three forms of HC-toxin that contain D-Ala, it still makes a minor form of HC-toxin that contains Gly in place of D-Ala. The addition of D-Ala to the culture medium restores production of the D-Ala-containing forms of HC-toxin by the toxG mutant. The toxG mutant has only partially reduced virulence. It is concluded that TOXG encodes an alanine racemase whose function is to synthesize D-Ala for incorporation into HC-toxin.     

Structure-Based Computational Study of Two Disease Resistance Gene Homologues (Hm1 and Hm2) in Maize (Zea mays L.) with Implications in Plant-Pathogen Interactions

The NADPH-dependent HC-toxin reductases (HCTR1 and 2) encoded by enzymatic class of disease resistance homologous genes (Hm1 and Hm2) protect maize by detoxifying a cyclic tetrapeptide, HC-toxin, secreted by the fungus Cochliobolus carbonum race 1(CCR1). Unlike the other classes'' resistance (R) genes, HCTR-mediated disease resistance is an inimitable mechanism where the avirulence (Avr) component from CCR1 is not involved in toxin degradation. In this study, we attempted to decipher cofactor (NADPH) recognition and mode of HC-toxin binding to HCTRs through molecular docking, molecular dynamics (MD) simulations and binding free energy calculation methods. The rationality and the stability of docked complexes were validated by 30-ns MD simulation. The binding free energy decomposition of enzyme-cofactor complex was calculated to find the driving force behind cofactor recognition. The overall binding free energies of HCTR1-NADPH and HCTR2-NADPH were found to be −616.989 and −16.9749 kJ mol⁻¹ respectively. The binding free energy decomposition revealed that the binding of NADPH to the HCTR1 is mainly governed by van der Waals and nonpolar interactions, whereas electrostatic terms play dominant role in stabilizing the binding mode between HCTR2 and NADPH. Further, docking analysis of HC-toxin with HCTR-NADPH complexes showed a distinct mode of binding and the complexes were stabilized by a strong network of hydrogen bond and hydrophobic interactions. This study is the first in silico attempt to unravel the biophysical and biochemical basis of cofactor recognition in enzymatic class of R genes in cereal crop maize.     

The cyclic peptide synthetase catalyzing HC-toxin production in the filamentous fungus Cochliobolus carbonum is encoded by a 15.7-kilobase open reading frame.

Race 1 of Cochliobolus carbonum, a fungal plant pathogen, owes its exceptional virulence on certain genotypes of maize to the production of HC-toxin, a cyclic tetrapeptide. Production of HC-toxin is controlled by a single known gene, TOX2. Race 1, but not races that do not make HC-toxin, contains two copies of a 22-kilobase (kb) region of chromosomal DNA that is required for HC-toxin biosynthesis and hence virulence. We have sequenced this 22-kb region and here show that it contains an open reading frame of 15.7 kb that encodes a multifunctional cyclic peptide synthetase of potential M(r)574,620. This gene, called HTS1, apparently contains no introns. The predicted gene product, HC-toxin synthetase (HTS), contains four amino acid-binding (adenylate-forming) domains that are highly similar to those found in other cyclic peptide synthetases and other adenylate-binding enzymes. The DNA sequence encodes tryptic peptides derived from two HC-toxin biosynthetic enzymes, HC-toxin synthetase 1 (HTS-1) and HC-toxin synthetase 2 (HTS-2), indicating that these two enzymes exist in vivo as part of a single polypeptide. Consistent with this, in some enzyme preparations antibodies against the enzyme HTS-2, which was originally purified as a protein with a subunit M(r) of 160,000, recognize a protein with an estimated subunit M(r) greater than 480,000.     

Polymorphic Chromosomes Bearing the Tox2 Locus in Cochliobolus carbonum Behave as Homologs during Meiosis

The HTS1 gene in the Tox2 locus of the fungal pathogen Cochliobolus carbonum race 1 is required for synthesis of a host-selective phytotoxin and for increased virulence on susceptible genotypes of maize. The locus is present in race 1 isolates but absent from isolates of the other races, which do not produce the toxin. By pulsed-field gel electrophoresis and Southern analysis with HTS1 sequences and chromosome-specific markers, the HTS1 gene was detected on a 4-Mb chromosome in one group of isolates and on a 2.3-Mb chromosome in another group, which lacked the 4-Mb chromosome. A chromosome-specific marker from C. heterostrophus hybridized to a 2.3-Mb chromosome in non-toxin-producing isolates and in toxin-producing isolates, including those with a 4-Mb chromosome. A marker from C. carbonum hybridized to the 4-Mb chromosome, but in isolates lacking the 4-Mb chromosome, this marker hybridized to a smaller, 2.0-Mb chromosome. Thus, the Tox2 locus is on different chromosomes in different groups of race 1 isolates. Single ascospore progeny from crosses between isolates having HTS1 on different chromosomes were analyzed for toxin-producing ability, virulence, and the presence and chromosomal location of HTS1. All progeny produced HC toxin in culture, incited race 1-type lesions on susceptible maize genotypes, and contained HTS1 sequences, as determined by PCR amplification with gene-specific primers. Analysis of the chromosomal complements of several progeny indicated that they all had only one Tox2-containing chromosome. Thus, despite their differences in size, these chromosomes behave as homologs during meiosis and may have arisen by a translocation.     

Genome Information of Methylobacterium oryzae,a Plant-Probiotic Methylotroph in the Phyllosphere

Pink-pigmented facultative methylotrophs in the Rhizobiales are widespread in the environment, and many Methylobacterium species associated with plants produce plant growth-promoting substances. To gain insights into the life style at the phyllosphere and the genetic bases of plant growth promotion, we determined and analyzed the complete genome sequence of Methylobacterium oryzae CBMB20^T, a strain isolated from rice stem. The genome consists of a 6.29-Mb chromosome and four plasmids, designated as pMOC1 to pMOC4. Among the 6,274 coding sequences in the chromosome, the bacterium has, besides most of the genes for the central metabolism, all of the essential genes for the assimilation and dissimilation of methanol that are either located in methylotrophy islands or dispersed. M. oryzae is equipped with several kinds of genes for adaptation to plant surfaces such as defense against UV radiation, oxidative stress, desiccation, or nutrient deficiency, as well as high proportion of genes related to motility and signaling. Moreover, it has an array of genes involved in metabolic pathways that may contribute to promotion of plant growth; they include auxin biosynthesis, cytokine biosynthesis, vitamin B₁₂ biosynthesis, urea metabolism, biosorption of heavy metals or decrease of metal toxicity, pyrroloquinoline quinone biosynthesis, 1-aminocyclopropane-1-carboxylate deamination, phosphate solubilization, and thiosulfate oxidation. Through the genome analysis of M. oryzae, we provide information on the full gene complement of M. oryzae that resides in the aerial parts of plants and enhances plant growth. The plant-associated lifestyle of M. oryzae pertaining to methylotrophy and plant growth promotion, and its potential as a candidate for a bioinoculant targeted to the phyllosphere and focused on phytostimulation are illuminated.     

Enzymatic Detoxification of HC-toxin, the Host-Selective Cyclic Peptide from Cochliobolus carbonum

Resistance to the fungal plant pathogen Cochliobolus carbonum race 1 and to its host-selective toxin, HC-toxin, is determined by Hm, a single dominant gene in the host plant maize, (Zea mays L). Radiolabeled HC-toxin of specific activity 70 milliCuries per millimole, prepared by feeding tritiated d,l-alanine to the fungus, was used to study its fate in maize leaf tissues. HC-toxin was converted by resistant leaf segments to a single compound, identified by mass spectrometry and nuclear magnetic resonance as the 8-hydroxy derivative of HC-toxin formed by reduction of the 8-keto group of 2-amino-9, 10-epoxy-8-oxo-decanoic acid, one of the amino acids in HC-toxin. Reduction of HC-toxin occurred in cell-free preparations from etiolated (Hm/hm) maize shoots, and the activity was sensitive to heat and proteolytic digestion, dependent on NADPH, and inhibited by p-hydroxymercuribenzoate and disulfiram. The enzyme (from the Hm/hm genotype) was partially purified by ammonium sulfate precipitation and diethylaminoethyl-ion exchange chromatography. By gel filtration chromatography, the enzyme had a molecular weight of 42,000. NADH was approximately 30% as effective as NADPH as a hydride donor, and flavin-containing cofactors had no effect on activity. When HC-toxin was introduced to maize leaf segments through the transpiration stream, leaf segments from both resistant and susceptible maize inactivated toxin equally well over a time-course of 9 hours. Although these data suggest no relationship between toxin metabolism and host selectivity, we discuss findings in apparent conflict with the current data and describe why the relationship between enzymatic reduction of HC-toxin and Hm remains unresolved.     

Plasmid Copy Number Underlies Adaptive Mutability in Bacteria

The origin of mutations under selection has been intensively studied using the Cairns-Foster system, in which cells of an Escherichia coli lac mutant are plated on lactose and give rise to 100 Lac⁺ revertants over several days. These revertants have been attributed variously to stress-induced mutagenesis of nongrowing cells or to selective improvement of preexisting weakly Lac⁺ cells with no mutagenesis. Most revertant colonies (90%) contain stably Lac⁺ cells, while others (10%) contain cells with an unstable amplification of the leaky mutant lac allele. Evidence is presented that both stable and unstable Lac⁺ revertant colonies are initiated by preexisting cells with multiple copies of the F′lac plasmid, which carries the mutant lac allele. The tetracycline analog anhydrotetracycline (AnTc) inhibits growth of cells with multiple copies of the tetA gene. Populations with tetA on their F′lac plasmid include rare cells with an elevated plasmid copy number and multiple copies of both the tetA and lac genes. Pregrowth of such populations with AnTc reduces the number of cells with multiple F′lac copies and consequently the number of Lac⁺ colonies appearing under selection. Revertant yield is restored rapidly by a few generations of growth without AnTc. We suggest that preexisting cells with multiple F′lac copies divide very little under selection but have enough energy to replicate their F′lac plasmids repeatedly until reversion initiates a stable Lac⁺ colony. Preexisting cells whose high-copy plasmid includes an internal lac duplication grow under selection and produce an unstable Lac⁺ colony. In this model, all revertant colonies are initiated by preexisting cells and cannot be stress induced.     

Regulation of Sex Determination in Mice by a Non-coding Genomic Region

To identify novel genomic regions that regulate sex determination, we utilized the powerful C57BL/6J-Y^POS (B6-Y^POS) model of XY sex reversal where mice with autosomes from the B6 strain and a Y chromosome from a wild-derived strain, Mus domesticus poschiavinus (Y^POS), show complete sex reversal. In B6-Y^POS, the presence of a 55-Mb congenic region on chromosome 11 protects from sex reversal in a dose-dependent manner. Using mouse genetic backcross designs and high-density SNP arrays, we narrowed the congenic region to a 1.62-Mb genomic region on chromosome 11 that confers 80% protection from B6-Y^POS sex reversal when one copy is present and complete protection when two copies are present. It was previously believed that the protective congenic region originated from the 129S1/SviMJ (129) strain. However, genomic analysis revealed that this region is not derived from 129 and most likely is derived from the semi-inbred strain POSA. We show that the small 1.62-Mb congenic region that protects against B6-Y^POS sex reversal is located within the Sox9 promoter and promotes the expression of Sox9, thereby driving testis development within the B6-Y^POS background. Through 30 years of backcrossing, this congenic region was maintained, as it promoted male sex determination and fertility despite the female-promoting B6-Y^POS genetic background. Our findings demonstrate that long-range enhancer regions are critical to developmental processes and can be used to identify the complex interplay between genome variants, epigenetics, and developmental gene regulation.     

Draft Genome Sequence of the Antifungal-Producing Plant-Benefiting Bacterium Burkholderia pyrrocinia CH-67

Burkholderia pyrrocinia CH-67 was isolated from forest soil as a biocontrol agent to be utilized in agriculture. Here, we report the 8.05-Mb draft genome sequence of this bacterium. Its genome contains genes involved in biosynthesis of secondary metabolites and plant growth promotion, which may contribute to probiotic effects on plants.     

Effects of Helminthosporium carbonum Toxin on Absorption of Solutes by Corn Roots

Susceptible corn roots exposed to the host-selective toxin of Helminthosporium carbonum took up and retained more NO₃⁻, Na⁺, Cl⁻, 3-o-methylglucose, and leucine than did control roots. Stimulatory effects on uptake were more pronounced with freshly cut roots than with roots that were washed and aged. Solutes were accumulated against a concentration gradient, and toxin-treated tissues developed a steeper gradient than did control tissues. Toxin affected both the low and high affinity uptake systems for Na⁺ and Cl⁻. Toxin did not affect uptake of Na₂⁻, K⁺, Ca²⁺, phosphate ion (H₂PO₄⁻ and HPO₄⁻), SO₄⁻, and glutamic acid. No toxin-induced leakage of any solute tested was detected within 5 to 6 hr after initial exposure to toxin. The data suggest that toxin from H. carbonum does not cause the general plasma membrane derangement caused by other host-selective toxins. Instead, H. carbonum toxin may cause specific changes in characteristics of the plasmalemma, which result in increased uptake of certain solutes.     

CCL3L1 and CCL4L1 chemokine genes are located in a segmental duplication at chromosome 17q12

Sixteen CC chemokine genes localize to a 2.06-Mb interval at 17q11.2-q12 on genomic contig NT_010799.13. Four of these genes comprise two closely related paralogous pairs: CCL3-CCL3L1 and CCL4-CCL4L1. Members within each pair share 95% sequence identity at both the genomic and the amino acid levels. One BAC clone (AC131056.5) on the contig with substantial internal sequence duplication contains two complete copies of CCL3L1 and CCL4L1 and one truncated copy of CCL3L1, while a partially overlapping clone (AC003976.1) contains one copy each of CCL3 and CCL4. Dot-matrix comparison of the regions of AC131056.5 with those of AC003976.1 containing the four genes reveals 90% sequence similarity over 37 kb. These observations support the idea that the multiple copies of CCL3L1 and CCL4L1 present in a single diploid genome are the result of segmental duplication.     

Finished sequence and assembly of the DUF1220-rich 1q21 region using a haploid human genome

Background

Although the reference human genome sequence was declared finished in 2003, some regions of the genome remain incomplete due to their complex architecture. One such region, 1q21.1-q21.2, is of increasing interest due to its relevance to human disease and evolution. Elucidation of the exact variants behind these associations has been hampered by the repetitive nature of the region and its incomplete assembly. This region also contains 238 of the 270 human DUF1220 protein domains, which are implicated in human brain evolution and neurodevelopment. Additionally, examinations of this protein domain have been challenging due to the incomplete 1q21 build. To address these problems, a single-haplotype hydatidiform mole BAC library (CHORI-17) was used to produce the first complete sequence of the 1q21.1-q21.2 region.

Results

We found and addressed several inaccuracies in the GRCh37sequence of the 1q21 region on large and small scales, including genomic rearrangements and inversions, and incorrect gene copy number estimates and assemblies. The DUF1220-encoding NBPF genes required the most corrections, with 3 genes removed, 2 genes reassigned to the 1p11.2 region, 8 genes requiring assembly corrections for DUF1220 domains (~91 DUF1220 domains were misassigned), and multiple instances of nucleotide changes that reassigned the domain to a different DUF1220 subtype. These corrections resulted in an overall increase in DUF1220 copy number, yielding a haploid total of 289 copies. Approximately 20 of these new DUF1220 copies were the result of a segmental duplication from 1q21.2 to 1p11.2 that included two NBPF genes. Interestingly, this duplication may have been the catalyst for the evolutionarily important human lineage-specific chromosome 1 pericentric inversion.

Conclusions

Through the hydatidiform mole genome sequencing effort, the 1q21.1-q21.2 region is complete and misassemblies involving inter- and intra-region duplications have been resolved. The availability of this single haploid sequence path will aid in the investigation of many genetic diseases linked to 1q21, including several associated with DUF1220 copy number variations. Finally, the corrected sequence identified a recent segmental duplication that added 20 additional DUF1220 copies to the human genome, and may have facilitated the chromosome 1 pericentric inversion that is among the most notable human-specific genomic landmarks.     

Distinct spermiogenic phenotypes underlie sperm elimination in the Segregation Distorter meiotic drive system

Segregation Distorter (SD) is a male meiotic drive system in Drosophila melanogaster. Males heterozygous for a selfish SD chromosome rarely transmit the homologous SD⁺ chromosome. It is well established that distortion results from an interaction between Sd, the primary distorting locus on the SD chromosome and its target, a satellite DNA called Rsp, on the SD⁺ chromosome. However, the molecular and cellular mechanisms leading to post-meiotic SD⁺ sperm elimination remain unclear. Here we show that SD/SD⁺ males of different genotypes but with similarly strong degrees of distortion have distinct spermiogenic phenotypes. In some genotypes, SD⁺ spermatids fail to fully incorporate protamines after the removal of histones, and degenerate during the individualization stage of spermiogenesis. In contrast, in other SD/SD⁺ genotypes, protamine incorporation appears less disturbed, yet spermatid nuclei are abnormally compacted, and mature sperm nuclei are eventually released in the seminal vesicle. Our analyses of different SD⁺ chromosomes suggest that the severity of the spermiogenic defects associates with the copy number of the Rsp satellite. We propose that when Rsp copy number is very high (> 2000), spermatid nuclear compaction defects reach a threshold that triggers a checkpoint controlling sperm chromatin quality to eliminate abnormal spermatids during individualization.     

A new plant sex-linked gene with high sequence diversity and possible introgression of the X copy

We describe patterns of DNA sequence diversity in a newly identified sex-linked gene, SlX9/SlY9, in Silene latifolia (Caryophyllaceae). The copies on both sex chromosomes seem to be functional, and each maps close to the respective X- and Y-linked copy of another sex-linked gene pair, SlCypX/SlCypY. The Y-linked copy has low diversity, similar to what has been found for several other Y-linked genes in S. latifolia, and consistent with the theoretical expectations of hitch-hiking processes occurring on a non-recombining chromosome. However, SlX9 has higher diversity than other genes on the S. latifolia X chromosome. We evaluate the hypothesis of introgression from the closely related species S. dioica as an explanation for the high sequence diversity observed.