Almost all lizard families in the pleurodont clade share the same XY system. This system was meticulously studied in Anolis carolinensis, where it shows a highly degenerated Y chromosome and a male-specific X chromosome dosage compensation mechanism. Corytophanids (casque-headed lizards) have been proposed as the only family in the pleurodont clade to lack the XY system. In this study, we worked with extensive genomic and transcriptomic data from Basiliscus vittatus, a member of the Corytophanidae family that inhabits the tropical rainforests of Mexico. We confirmed that B. vittatus underwent a sex chromosome system turnover, which consisted in the loss of the pleurodont XY system and the gain of a new pair of XY chromosomes that are orthologous to chicken chromosome 17. We estimated the origin of the sex chromosome system to have occurred ∼63 Ma in the ancestor of corytophanids. Moreover, we identified 12 XY gametologues with particular attributes, such as functions related to the membrane and intracellular trafficking, very low expression levels, blood specificity, and incomplete dosage compensation in males.

We used two-dimensional quantitative trait locus analysis to identify interacting genetic loci that contribute to the native airway constrictor hyperresponsiveness to methacholine that characterizes A/J mice, relative to C57BL/6J mice. We quantified airway responsiveness to intravenous methacholine boluses in eighty-eight (C57BL/6J X A/J) F₂ and twenty-seven (A/J X C57BL/6J) F₂ mice as well as ten A/J mice and six C57BL/6J mice; all studies were performed in male mice. Mice were genotyped at 384 SNP markers, and from these data two-QTL analyses disclosed one pair of interacting loci on chromosomes 11 and 18; the homozygous A/J genotype at each locus constituted the genetic interaction linked to the hyperresponsive A/J phenotype. Bioinformatic network analysis of potential interactions among proteins encoded by genes in the linked regions disclosed two high priority subnetworks--Myl7, Rock1, Limk2; and Npc1, Npc1l1. Evidence in the literature supports the possibility that either or both networks could contribute to the regulation of airway constrictor responsiveness. Together, these results should stimulate evaluation of the genetic contribution of these networks in the regulation of airway responsiveness in humans.

Comprehensive identification of the acquired mutations that cause common cancers will require genomic analyses of large sets of tumor samples. Typically, the tissue material available from tumor specimens is limited, which creates a demand for accurate template amplification. We therefore evaluated whether phi29-mediated whole genome amplification introduces false positive structural mutations by massive mate-pair sequencing of a normal human genome before and after such amplification. Multiple displacement amplification led to a decrease in clone coverage and an increase by two orders of magnitude in the prevalence of inversions, but did not increase the prevalence of translocations. While multiple strand displacement amplification may find uses in translocation analyses, it is likely that alternative amplification strategies need to be developed to meet the demands of cancer genomics.

Proteins play a major role in the three-dimensional organization of nuclear genome and its function. While histones arrange DNA into a nucleosome fiber, other proteins contribute to higher-order chromatin structures in interphase nuclei, and mitotic/meiotic chromosomes. Despite the key role of proteins in maintaining genome integrity and transferring hereditary information to daughter cells and progenies, the knowledge about their function remains fragmentary. This is particularly true for the proteins of condensed chromosomes and, in particular, chromosomes of plants. Here, we purified barley mitotic metaphase chromosomes by a flow cytometric sorting and characterized their proteins. Peptides from tryptic protein digests were fractionated either on a cation exchanger or reversed-phase microgradient system before liquid chromatography coupled to tandem mass spectrometry. Chromosomal proteins comprising almost 900 identifications were classified based on a combination of software prediction, available database localization information, sequence homology, and domain representation. A biological context evaluation indicated the presence of several groups of abundant proteins including histones, topoisomerase 2, POLYMERASE 2, condensin subunits, and many proteins with chromatin-related functions. Proteins involved in processes related to DNA replication, transcription, and repair as well as nucleolar proteins were found. We have experimentally validated the presence of FIBRILLARIN 1, one of the nucleolar proteins, on metaphase chromosomes, suggesting that plant chromosomes are coated with proteins during mitosis, similar to those of human and animals. These results improve significantly the knowledge of plant chromosomal proteins and provide a basis for their functional characterization and comparative phylogenetic analyses.

The Sc2.0 project is building a eukaryotic synthetic genome from scratch. A major milestone has been achieved with all individual Sc2.0 chromosomes assembled. Here, we describe the consolidation of multiple synthetic chromosomes using advanced endoreduplication intercrossing with tRNA expression cassettes to generate a strain with 6.5 synthetic chromosomes. The 3D chromosome organization and transcript isoform profiles were evaluated using Hi-C and long-read direct RNA sequencing. We developed CRISPR Directed Biallelic URA3-assisted Genome Scan, or "CRISPR D-BUGS," to map phenotypic variants caused by specific designer modifications, known as "bugs." We first fine-mapped a bug in synthetic chromosome II (synII) and then discovered a combinatorial interaction associated with synIII and synX, revealing an unexpected genetic interaction that links transcriptional regulation, inositol metabolism, and tRNASerCGA abundance. Finally, to expedite consolidation, we employed chromosome substitution to incorporate the largest chromosome (synIV), thereby consolidating >50% of the Sc2.0 genome in one strain.

Near-haploid chromosome numbers have been found in less than 1% of cytogenetically reported tumors, but seem to be more common in certain neoplasms including the malignant cartilage-producing tumor chondrosarcoma. By a literature survey of published karyotypes from chondrosarcomas we could confirm that loss of chromosomes resulting in hyperhaploid-hypodiploid cells is common and that these cells may polyploidize. Sixteen chondrosarcomas were investigated by single nucleotide polymorphism (SNP) array and the majority displayed SNP patterns indicative of a hyperhaploid-hypodiploid origin, with or without subsequent polyploidization. Except for chromosomes 5, 7, 19, 20 and 21, autosomal loss of heterozygosity was commonly found, resulting from chromosome loss and subsequent duplication of monosomic chromosomes giving rise to uniparental disomy. Additional gains, losses and rearrangements of genetic material, and even repeated rounds of polyploidization, may affect chondrosarcoma cells resulting in highly complex karyotypes. Loss of chromosomes and subsequent polyploidization was not restricted to a particular chondrosarcoma subtype and, although commonly found in chondrosarcoma, binucleated cells did not seem to be involved in these events.

Sex chromosomes evolve distinctive types of chromatin from a pair of ancestral autosomes that are usually euchromatic. In Drosophila, the dosage-compensated X becomes enriched for hyperactive chromatin in males (mediated by H4K16ac), while the Y chromosome acquires silencing heterochromatin (enriched for H3K9me2/3). Drosophila autosomes are typically mostly euchromatic but the small dot chromosome has evolved a heterochromatin-like milieu (enriched for H3K9me2/3) that permits the normal expression of dot-linked genes, but which is different from typical pericentric heterochromatin. In Drosophila busckii, the dot chromosomes have fused to the ancestral sex chromosomes, creating a pair of 'neo-sex' chromosomes. Here we collect genomic, transcriptomic and epigenomic data from D. busckii, to investigate the evolutionary trajectory of sex chromosomes from a largely heterochromatic ancestor. We show that the neo-sex chromosomes formed <1 million years ago, but nearly 60% of neo-Y linked genes have already become non-functional. Expression levels are generally lower for the neo-Y alleles relative to their neo-X homologs, and the silencing heterochromatin mark H3K9me2, but not H3K9me3, is significantly enriched on silenced neo-Y genes. Despite rampant neo-Y degeneration, we find that the neo-X is deficient for the canonical histone modification mark of dosage compensation (H4K16ac), relative to autosomes or the compensated ancestral X chromosome, possibly reflecting constraints imposed on evolving hyperactive chromatin in an originally heterochromatic environment. Yet, neo-X genes are transcriptionally more active in males, relative to females, suggesting the evolution of incipient dosage compensation on the neo-X. Our data show that Y degeneration proceeds quickly after sex chromosomes become established through genomic and epigenetic changes, and are consistent with the idea that the evolution of sex-linked chromatin is influenced by its ancestral configuration.

B chromosomes constitute a heterogeneous mixture of genomic parasites that are sometimes derived intraspecifically from the standard genome of the host species, but result from interspecific hybridization in other cases. The mode of origin determines the DNA content, with the B chromosomes showing high similarity with the A genome in the first case, but presenting higher similarity with a different species in the second. The characid fish Moenkhausia sanctaefilomenae harbours highly invasive B chromosomes, which are present in all populations analyzed to date in the Parana and Tietê rivers. To investigate the origin of these B chromosomes, we analyzed two natural populations: one carrying B chromosomes and the other lacking them, using a combination of molecular cytogenetic techniques, nucleotide sequence analysis and high-throughput sequencing (Illumina HiSeq2000). Our results showed that i) B chromosomes have not yet reached the Paranapanema River basin; ii) B chromosomes are mitotically unstable; iii) there are two types of B chromosomes, the most frequent of which is lightly C-banded (similar to euchromatin in A chromosomes) (B1), while the other is darkly C-banded (heterochromatin-like) (B2); iv) the two B types contain the same tandem repeat DNA sequences (18S ribosomal DNA, H3 histone genes, MS3 and MS7 satellite DNA), with a higher content of 18S rDNA in the heterochromatic variant; v) all of these repetitive DNAs are present together only in the paracentromeric region of autosome pair no. 6, suggesting that the B chromosomes are derived from this A chromosome; vi) the two B chromosome variants show MS3 sequences that are highly divergent from each other and from the 0B genome, although the B2-derived sequences exhibit higher similarity with the 0B genome (this suggests an independent origin of the two B variants, with the less frequent, B2 type presumably being younger); and vii) the dN/dS ratio for the H3.2 histone gene is almost 4-6 times higher for B chromosomes than for A chromosome sequences, suggesting that purifying selection is relaxed for the DNA sequences located on the B chromosomes, presumably because they are mostly inactive.

The karyotype of most birds has remained considerably stable during more than 100 million years' evolution, except for some groups, such as parrots. The evolutionary processes and underlying genetic mechanism of chromosomal rearrangements in parrots, however, are poorly understood. Here, using chromosome-level assemblies of four parrot genomes, we uncover frequent chromosome fusions and fissions, with most of them occurring independently among lineages. The increased activities of chromosomal rearrangements in parrots are likely associated with parrot-specific loss of two genes, ALC1 and PARP3, that have known functions in the repair of double-strand breaks and maintenance of genome stability. We further find that the fusion of the ZW sex chromosomes and chromosome 11 has created a pair of neo-sex chromosomes in the ancestor of parrots, and the chromosome 25 has been further added to the sex chromosomes in monk parakeet. Together, the combination of our genomic and cytogenetic analyses characterizes the complex evolutionary history of chromosomal rearrangements and sex chromosomes in parrots.

Pyruvate dehydrogenase complex (PDC) deficiency caused by mutations in the X-linked PDHA1 gene has a broad clinical presentation, and the pattern of X-chromosome inactivation has been proposed as a major factor contributing to its variable expressivity in heterozygous females. Here, we report the first set of monozygotic twin females with PDC deficiency, caused by a novel, de novo heterozygous missense mutation in exon 11 of PDHA1 (NM_000284.3: c.1100A>T). Both twins presented in infancy with a similar clinical phenotype including developmental delay, episodes of hypotonia or encephalopathy, epilepsy, and slowly progressive motor impairment due to pyramidal, extrapyramidal, and cerebellar involvement. However, they exhibited clear differences in disease severity that correlated well with residual PDC activities (approximately 60% and 20% of mean control values, respectively) and levels of immunoreactive E1α subunit in cultured skin fibroblasts. To address whether the observed clinical and biochemical differences could be explained by the pattern of X-chromosome inactivation, we undertook an androgen receptor assay in peripheral blood. In the less severely affected twin, a significant bias in the relative activity of the two X chromosomes with a ratio of approximately 75:25 was detected, while the ratio was close to 50:50 in the other twin. Although it may be difficult to extrapolate these results to other tissues, our observation provides further support to the hypothesis that the pattern of X-chromosome inactivation may influence the phenotypic expression of the same mutation in heterozygous females and broadens the clinical and genetic spectrum of PDC deficiency.

How the avian sex chromosomes first evolved from autosomes remains elusive as 100 million years (My) of divergence and degeneration obscure their evolutionary history. The Sylvioidea group of songbirds is interesting for understanding avian sex chromosome evolution because a chromosome fusion event ∼24 Ma formed "neo-sex chromosomes" consisting of an added (new) and an ancestral (old) part. Here, we report the complete female genome (ZW) of one Sylvioidea species, the great reed warbler (Acrocephalus arundinaceus). Our long-read assembly shows that the added region has been translocated to both Z and W, and whereas the added-Z has retained its gene order the added-W part has been heavily rearranged. Phylogenetic analyses show that recombination between the homologous added-Z and -W regions continued after the fusion event, and that recombination suppression across this region took several million years to be completed. Moreover, recombination suppression was initiated across multiple positions over the added-Z, which is not consistent with a simple linear progression starting from the fusion point. As expected following recombination suppression, the added-W show signs of degeneration including repeat accumulation and gene loss. Finally, we present evidence for nonrandom maintenance of slowly evolving and dosage-sensitive genes on both ancestral- and added-W, a process causing correlated evolution among orthologous genes across broad taxonomic groups, regardless of sex linkage.

The orb web is a remarkable example of animal architecture that is observed in families of spiders that diverged over 200 million years ago. While several genomes exist for araneid orb-weavers, none exist for other orb-weaving families, hampering efforts to investigate the genetic basis of this complex behavior. Here we present a chromosome-level genome assembly for the cribellate orb-weaving spider Uloborus diversus. The assembly reinforces evidence of an ancient arachnid genome duplication and identifies complete open reading frames for every class of spidroin gene, which encode the proteins that are the key structural components of spider silks. We identified the 2 X chromosomes for U. diversus and identify candidate sex-determining loci. This chromosome-level assembly will be a valuable resource for evolutionary research into the origins of orb-weaving, spidroin evolution, chromosomal rearrangement, and chromosomal sex determination in spiders.

Sex chromosomes originated from autosomes but have evolved a highly specialized chromatin structure. Drosophila Y chromosomes are composed entirely of silent heterochromatin, while male X chromosomes have highly accessible chromatin and are hypertranscribed as a result of dosage compensation. Here, we dissect the molecular mechanisms and functional pressures driving heterochromatin formation and dosage compensation of the recently formed neo-sex chromosomes of Drosophila miranda. We show that the onset of heterochromatin formation on the neo-Y is triggered by an accumulation of repetitive DNA. The neo-X has evolved partial dosage compensation and we find that diverse mutational paths have been utilized to establish several dozen novel binding consensus motifs for the dosage compensation complex on the neo-X, including simple point mutations at pre-binding sites, insertion and deletion mutations, microsatellite expansions, or tandem amplification of weak binding sites. Spreading of these silencing or activating chromatin modifications to adjacent regions results in massive mis-expression of neo-sex linked genes, and little correspondence between functionality of genes and their silencing on the neo-Y or dosage compensation on the neo-X. Intriguingly, the genomic regions being targeted by the dosage compensation complex on the neo-X and those becoming heterochromatic on the neo-Y show little overlap, possibly reflecting different propensities along the ancestral chromosome that formed the sex chromosome to adopt active or repressive chromatin configurations. Our findings have broad implications for current models of sex chromosome evolution, and demonstrate how mechanistic constraints can limit evolutionary adaptations. Our study also highlights how evolution can follow predictable genetic trajectories, by repeatedly acquiring the same 21-bp consensus motif for recruitment of the dosage compensation complex, yet utilizing a diverse array of random mutational changes to attain the same phenotypic outcome.

The dioecious genus Spinacia is thought to include two wild relatives (S. turkestanica Ilj. and S. tetrandra Stev.) of cultivated spinach (S. oleracea L.). In this study, nuclear and chloroplast sequences from 21 accessions of Spinacia germplasm and six spinach cultivars or lines were subjected to phylogenetic analysis to define the relationships among the three species. Maximum-likelihood sequence analysis suggested that the Spinacia plant samples could be classified into two monophyletic groups (Group 1 and Group 2): Group 1 consisted of all accessions, cultivars, and lines of S. oleracea L. and S. turkestanica Ilj. and two of five S. tetrandra Stev. accessions, whereas Group 2 was composed of the three remaining S. tetrandra Stev. accessions. By using flow cytometry, we detected a distinct difference in nuclear genome size between the groups. Group 2 also was characterized by a sexual dimorphism in inflorescence structure, which was not observed in Group 1. Interspecific crosses between the groups produced hybrids with drastically reduced pollen fertility and showed that the male is the heterogametic sex (XY) in Group 2, as is the case in S. oleracea L. (Group 1). Cytogenetic and DNA marker analyses suggested that Group 1 and Group 2 have homomorphic and heteromorphic sex chromosome pairs (XY), respectively, and that the sex chromosome pairs of the two groups evolved from a common ancestral pair. Our data suggest that the Spinacia genus may serve as a good model for investigation of evolutionary mechanisms underlying the emergence of heteromorphic sex chromosome pairs from ancestral homomorphic pairs.

Supernumerary (B) chromosomes are dispensable genomic elements occurring frequently among grasshoppers. Most B chromosomes are enriched with repetitive DNAs, including satellite DNAs (satDNAs) that could be implicated in their evolution. Although studied in some species, the specific ancestry of B chromosomes is difficult to ascertain and it was determined in only a few examples. Here we used bioinformatics and cytogenetics to characterize the composition and putative ancestry of B chromosomes in three grasshopper species, Rhammatocerus brasiliensis, Schistocerca rubiginosa, and Xyleus discoideus angulatus. Using the RepeatExplorer pipeline we searched for the most abundant satDNAs in Illumina sequenced reads, and then we generated probes used in fluorescent in situ hybridization (FISH) to determine chromosomal position. We used this information to infer ancestry and the events that likely occurred at the origin of B chromosomes. We found twelve, nine, and eighteen satDNA families in the genomes of R. brasiliensis, S. rubiginosa, and X. d. angulatus, respectively. Some satDNAs revealed clustered organization on A and B chromosomes varying in number of sites and position along chromosomes. We did not find specific satDNA occurring in the B chromosome. The satDNAs shared among A and B chromosomes support the idea of putative intraspecific ancestry from small autosomes in the three species, i.e., pair S11 in R. brasiliensis, pair S9 in S. rubiginosa, and pair S10 in X. d. angulatus. The possibility of involvement of other chromosomal pairs in B chromosome origin is also hypothesized. Finally, we discussed particular aspects in composition, origin, and evolution of the B chromosome for each species.

Genetic predisposition controlled by susceptibility quantitative trait loci (QTLs) contributes to a large proportion of common cancers. Studies of genetics of cancer susceptibility, however, did not address systematically the relationship between susceptibility to cancers in different organs. We present five sets of data on genetic architecture of colon and lung cancer susceptibility in mice, humans and rats. They collectively show that the majority of genes for colon and lung cancer susceptibility are linked pair-wise and are likely identical or related. Four CcS/Dem recombinant congenic strains, each differing from strain BALB/cHeA by a different small random subset of ±12.5% of genes received from strain STS/A, suggestively show either extreme susceptibility or extreme resistance for both colon and lung tumors, which is unlikely if the two tumors were controlled by independent susceptibility genes. Indeed, susceptibility to lung cancer (Sluc) loci underlying the extreme susceptibility or resistance of such CcS/Dem strains, mapped in 226 (CcS-10 x CcS-19)F2 mice, co-localize with susceptibility to colon cancer (Scc) loci. Analysis of additional Sluc loci that were mapped in OcB/Dem strains and Scc loci in CcS/Dem strains, respectively, shows their widespread pair-wise co-localization (P  =  0.0036). Finally, the majority of published human and rat colon cancer susceptibility genes map to chromosomal regions homologous to mouse Sluc loci. 12/12 mouse Scc loci, 9/11 human and 5/7 rat colon cancer susceptibility loci are close to a Sluc locus or its homologous site, forming 21 clusters of lung and colon cancer susceptibility genes from one, two or three species. Our data shows that cancer susceptibility QTLs can have much broader biological effects than presently appreciated. It also demonstrates the power of mouse genetics to predict human susceptibility genes. Comparison of molecular mechanisms of susceptibility genes that are organ-specific and those with trans-organ effects can provide a new dimension in understanding individual cancer susceptibility.

Sex determination (SD) is not conserved among teleost fishes and can even differ between populations of the same species. Across the outstandingly species-rich fish family Cichlidae, more and more SD systems are being discovered. Still, the picture of SD evolution in this group is far from being complete. Lake Tanganyika and its affluent rivers are home to Astatotilapia burtoni, which belongs to the extremely successful East African cichlid lineage Haplochromini. Previously, in different families of an A. burtoni laboratory strain, an XYW system and an XY system have been described. The latter was also found in a second laboratory strain. In a laboratory-reared family descending from a population of the species' southern distribution, a second XY system was discovered. Yet, an analysis of sex chromosomes for the whole species distribution is missing. Here, we examined the genomes of 11 natural populations of A. burtoni, encompassing a wide range of its distribution, for sex-linked regions. We did not detect signs of differentiated sex chromosomes and also not the previously described sex chromosomal systems present in laboratory lines, suggesting different SD systems in the same species under natural and (long-term) artificial conditions. We suggest that SD in A. burtoni is more labile than previously assumed and consists of a combination of non-genetic, polygenic, or poorly differentiated sex chromosomes.

Proteins secreted by pathogens during host colonization largely determine the outcome of pathogen-host interactions and are commonly called 'effectors'. In fungal plant pathogens, coordinated transcriptional up-regulation of effector genes is a key feature of pathogenesis and effectors are often encoded in genomic regions with distinct repeat content, histone code and rate of evolution. In the tomato pathogen Fusarium oxysporum f. sp. lycopersici (Fol), effector genes reside on one of four accessory chromosomes, known as the 'pathogenicity' chromosome, which can be exchanged between strains through horizontal transfer. The three other accessory chromosomes in the Fol reference strain may also be important for virulence towards tomato. Expression of effector genes in Fol is highly up-regulated upon infection and requires Sge1, a transcription factor encoded on the core genome. Interestingly, the pathogenicity chromosome itself contains 13 predicted transcription factor genes and for all except one, there is a homolog on the core genome. We determined DNA binding specificity for nine transcription factors using oligonucleotide arrays. The binding sites for homologous transcription factors were highly similar, suggesting that extensive neofunctionalization of DNA binding specificity has not occurred. Several DNA binding sites are enriched on accessory chromosomes, and expression of FTF1, its core homolog FTF2 and SGE1 from a constitutive promoter can induce expression of effector genes. The DNA binding sites of only these three transcription factors are enriched among genes up-regulated during infection. We further show that Ftf1, Ftf2 and Sge1 can activate transcription from their binding sites in yeast. RNAseq analysis revealed that in strains with constitutive expression of FTF1, FTF2 or SGE1, expression of a similar set of plant-responsive genes on the pathogenicity chromosome is induced, including most effector genes. We conclude that the Fol pathogenicity chromosome may be partially transcriptionally autonomous, but there are also extensive transcriptional connections between core and accessory chromosomes.

Next-generation mate-pair sequencing (MPS) has revealed that many constitutional complex chromosomal rearrangements (CCRs) are associated with local shattering of chromosomal regions (chromothripsis). Although MPS promises to identify the molecular basis of the abnormal phenotypes associated with many CCRs, none of the reported mate-pair sequenced complex rearrangements have been simultaneously studied with state-of-the art molecular cytogenetic techniques. Here, we studied chromothripsis-associated CCR involving chromosomes 2, 5 and 7, associated with global developmental and psychomotor delay and severe speech disorder. We identified three truncated genes: CDH12, DGKB and FOXP2, confirming the role of FOXP2 in severe speech disorder, and suggestive roles of CDH12 and/or DGKB for the global developmental and psychomotor delay. Our study confirmes the power of MPS for detecting breakpoints and truncated genes at near nucleotide resolution in chromothripsis. However, only by combining MPS data with conventional G-banding and extensive fluorescence in situ hybridizations could we delineate the precise structure of the derivative chromosomes.

The significance of introgression as an evolutionary force shaping natural populations is well established, especially in animal and plant systems. However, the abundance and size of introgression tracts, and to what degree interspecific gene flow is the result of adaptive processes, are largely unknown. In this study, we present medium coverage genomic data from species of the filamentous ascomycete Neurospora, and we use comparative genomics to investigate the introgression landscape at the genomic level in this model genus. We revealed one large introgression tract in each of the three investigated phylogenetic lineages of Neurospora tetrasperma (sizes of 5.6 Mbp, 5.2 Mbp, and 4.1 Mbp, respectively). The tract is located on the chromosome containing the locus conferring sexual identity, the mating-type (mat) chromosome. The region of introgression is confined to the region of suppressed recombination and is found on one of the two mat chromosomes (mat a). We used Bayesian concordance analyses to exclude incomplete lineage sorting as the cause for the observed pattern, and multilocus genealogies from additional species of Neurospora show that the introgression likely originates from two closely related, freely recombining, heterothallic species (N. hispaniola and N. crassa/N. perkinsii). Finally, we investigated patterns of molecular evolution of the mat chromosome in Neurospora, and we show that introgression is correlated with reduced level of molecular degeneration, consistent with a shorter time of recombination suppression. The chromosome specific (mat) and allele specific (mat a) introgression reported herein comprise the largest introgression tracts reported to date from natural populations. Furthermore, our data contradicts theoretical predictions that introgression should be less likely on sex-determining chromosomes. Taken together, the data presented herein advance our general understanding of introgression as a force shaping eukaryotic genomes.