The mating-type chromosome in the filamentous ascomycete Neurospora tetrasperma represents a model for early evolution of sex chromosomes

We combined gene divergence data, classical genetics, and phylogenetics to study the evolution of the mating-type chromosome in the filamentous ascomycete Neurospora tetrasperma. In this species, a large non-recombining region of the mating-type chromosome is associated with a unique fungal life cycle where self-fertility is enforced by maintenance of a constant state of heterokaryosis. Sequence divergence between alleles of 35 genes from the two single mating-type component strains (i.e. the homokaryotic mat A or mat a-strains), derived from one N. tetrasperma heterokaryon (mat A+mat a), was analyzed. By this approach we were able to identify the boundaries and size of the non-recombining region, and reveal insight into the history of recombination cessation. The non-recombining region covers almost 7 Mbp, over 75% of the chromosome, and we hypothesize that the evolution of the mating-type chromosome in this lineage involved two successive events. The first event was contemporaneous with the split of N. tetrasperma from a common ancestor with its outcrossing relative N. crassa and suppressed recombination over at least 6.6 Mbp, and the second was confined to a smaller region in which recombination ceased more recently. In spite of the early origin of the first "evolutionary stratum", genealogies of five genes from strains belonging to an additional N. tetrasperma lineage indicate independent initiations of suppressed recombination in different phylogenetic lineages. This study highlights the shared features between the sex chromosomes found in the animal and plant kingdoms and the fungal mating-type chromosome, despite fungi having no separate sexes. As is often found in sex chromosomes of plants and animals, recombination suppression of the mating-type chromosome of N. tetrasperma involved more than one evolutionary event, covers the majority of the mating-type chromosome and is flanked by distal regions with obligate crossovers.     

Degeneration in codon usage within the region of suppressed recombination in the mating-type chromosomes of Neurospora tetrasperma

The origin and early evolution of sex chromosomes are currently poorly understood. The Neurospora tetrasperma mating-type (mat) chromosomes have recently emerged as a model system for the study of early sex chromosome evolution, since they contain a young (<6 million years ago [Mya]), large (>6.6-Mb) region of suppressed recombination. Here we examined preferred-codon usage in 290 genes (121,831 codon positions) in order to test for early signs of genomic degeneration in N. tetrasperma mat chromosomes. We report several key findings about codon usage in the region of recombination suppression, including the following: (i) this region has been subjected to marked and largely independent degeneration among gene alleles; (ii) the level of degeneration is magnified over longer periods of recombination suppression; and (iii) both mat a and mat A chromosomes have been subjected to deterioration. The frequency of shifts from preferred codons to nonpreferred codons is greater for shorter genes than for longer genes, suggesting that short genes play an especially significant role in early sex chromosome evolution. Furthermore, we show that these degenerative changes in codon usage are best explained by altered selection efficiency in the recombinationally suppressed region. These findings demonstrate that the fungus N. tetrasperma provides an effective system for the study of degenerative genomic changes in young regions of recombination suppression in sex-regulating chromosomes.     

Evidence of the accumulation of allele-specific non-synonymous substitutions in the young region of recombination suppression within the mating-type chromosomes of Neurospora tetrasperma

Currently, little is known about the origin and early evolution of sex chromosomes. This is largely due to the fact that ancient non-recombining sex chromosomes are highly degenerated, and thus provide little information about the early genomic events in their evolution. The Neurospora tetrasperma mating-type (mat) chromosomes contain a young (<6 Mya) and large region (>6.6 Mb) of suppressed recombination, thereby providing a model system to study early stages of sex chromosome evolution. Here, we examined alleles of 207 genes located on the N. tetrasperma mat a and mat A chromosomes to test for signs of genomic alterations at the protein level in the young region of recombination suppression. We report that the N. tetrasperma mat a and mat A chromosomes have each independently accumulated allele-specific non-synonymous codon substitutions in a time-dependent, and gene-specific manner in the recombinationally suppressed region. In addition, examination of the ratio (ω) of non-synonymous substitutions (dN) to synonymous substitutions (dS) using maximum likelihood analyses, indicates that such changes are associated with relaxed purifying selection, a finding consistent with genomic degeneration. We also reveal that sex specific biases in mutation rates or selection pressures are not necessary for genomic alterations in sex chromosomes, and that recombination suppression in itself is sufficient to explain these results. The present findings extend our current understanding of genomic events associated within the young region of recombination suppression in these fungal sex-regulating chromosomes.     

An HMM-Based Comparative Genomic Framework for Detecting Introgression in Eukaryotes

One outcome of interspecific hybridization and subsequent effects of evolutionary forces is introgression, which is the integration of genetic material from one species into the genome of an individual in another species. The evolution of several groups of eukaryotic species has involved hybridization, and cases of adaptation through introgression have been already established. In this work, we report on PhyloNet-HMM—a new comparative genomic framework for detecting introgression in genomes. PhyloNet-HMM combines phylogenetic networks with hidden Markov models (HMMs) to simultaneously capture the (potentially reticulate) evolutionary history of the genomes and dependencies within genomes. A novel aspect of our work is that it also accounts for incomplete lineage sorting and dependence across loci. Application of our model to variation data from chromosome 7 in the mouse (Mus musculus domesticus) genome detected a recently reported adaptive introgression event involving the rodent poison resistance gene Vkorc1, in addition to other newly detected introgressed genomic regions. Based on our analysis, it is estimated that about 9% of all sites within chromosome 7 are of introgressive origin (these cover about 13 Mbp of chromosome 7, and over 300 genes). Further, our model detected no introgression in a negative control data set. We also found that our model accurately detected introgression and other evolutionary processes from synthetic data sets simulated under the coalescent model with recombination, isolation, and migration. Our work provides a powerful framework for systematic analysis of introgression while simultaneously accounting for dependence across sites, point mutations, recombination, and ancestral polymorphism.     

Gene genealogies indicates abundant gene conversions and independent evolutionary histories of the mating-type chromosomes in the evolutionary history of <Emphasis Type="Italic">Neurospora tetrasperma</Emphasis>

Background  

The self-fertile filamentous ascomycete Neurospora tetrasperma contains a large (~7 Mbp) and young (< 6 MYA) region of suppressed recombination within its mating-type (mat) chromosomes. The objective of the present study is to reveal the evolutionary history, including key genomic events, associated with the various regions of the mat chromosomes among ten strains representing all the nine known species (lineages) contained within the N. tetrasperma species complex.     

Recent and Massive Expansion of the Mating-Type-Specific Region in the Smut Fungus Microbotryum

The presence of large genomic regions with suppressed recombination (SR) is a key shared property of some sex- and mating-type determining (mat) chromosomes identified to date in animals, plants, and fungi. Why such regions form and how they evolve remain central questions in evolutionary genetics. The smut fungus Microbotryum lychnis-dioicae is a basidiomycete fungus in which dimorphic mat chromosomes have been reported, but the size, age, and evolutionary dynamics of the SR region remains unresolved. To identify the SR region in M. lychnis-dioicae and to study its evolution, we sequenced 12 genomes (6 per mating type) of this species and identified the genomic contigs that show fixed sequence differences between the mating types. We report that the SR region spans more than half of the mat chromosome (>2.3 Mbp) and that it is of very recent origin (∼2 × 10⁶ years) as the average sequence divergence between mating types was only 2% in the SR region. This contrasts with a much higher divergence in and around the mating-type determining pheromone receptor locus in the SR, suggesting a recent and massive expansion of the SR region. Our results comprise the first reported case of recent massive SR expansion documented in a basidiomycete fungus.     

Molecular karyotype of the white rot fungus Pleurotus ostreatus.

The white rot fungus Pleurotus ostreatus is an edible basidiomycete with increasing agricultural and biotechnological importance. Genetic manipulation and breeding of this organism are restricted because of the lack of knowledge about its genomic structure. In this study, we analyzed the genomic constitution of P. ostreatus by using pulsed-field gel electrophoresis optimized for the separation of its chromosomes. We have determined that it contains 11 pairs of chromosomes with sizes ranging from 1.4 to 4.7 Mbp. In addition to chromosome separation, the use of single-copy DNA probes allowed us to resolve the ambiguities caused by chromosome comigration. When the two nuclei present in the dikaryon were separated by protoplasting, analysis of their karyotypes revealed length polymorphisms affecting various chromosomes. This is, to our knowledge, the clearest chromosome separation available for this species.     

Molecular Karyotype of the White Rot Fungus Pleurotus ostreatus

The white rot fungus Pleurotus ostreatus is an edible basidiomycete with increasing agricultural and biotechnological importance. Genetic manipulation and breeding of this organism are restricted because of the lack of knowledge about its genomic structure. In this study, we analyzed the genomic constitution of P. ostreatus by using pulsed-field gel electrophoresis optimized for the separation of its chromosomes. We have determined that it contains 11 pairs of chromosomes with sizes ranging from 1.4 to 4.7 Mbp. In addition to chromosome separation, the use of single-copy DNA probes allowed us to resolve the ambiguities caused by chromosome comigration. When the two nuclei present in the dikaryon were separated by protoplasting, analysis of their karyotypes revealed length polymorphisms affecting various chromosomes. This is, to our knowledge, the clearest chromosome separation available for this species.     

Multilocus nuclear sequences reveal intra- and interspecific relationships among chromosomally polymorphic species of cactophilic Drosophila

Drosophila mojavensis and Drosophila arizonae, a pair of sibling species endemic to North America, constitute an important model system to study ecological genetics and the evolution of reproductive isolation. This species pair can produce fertile hybrids in some crosses and are sympatric in a large part of their ranges. Despite the potential for hybridization in nature, however, evidence of introgression has not been rigorously sought. Further, the evolutionary relationships within and among the geographically distant populations of the two species have not been characterized in detail using high-resolution molecular studies. Both species have six chromosomes: five large acrocentrics and one 'dot' chromosome. Fixed inversion differences between the species exist in three chromosomes (X, 2 and 3) while three are colinear (4, 5 and 6), suggesting that were introgression to occur, it would be most likely in the colinear chromosomes. We utilized nucleotide sequence variation at multiple loci on five chromosomes to test for evidence of introgression, and to test various scenarios for the evolutionary relationships of these two species and their populations. While we do not find evidence of recent introgression, loci in the colinear chromosomes appear to have participated in exchange in the past. We also found considerable population structure within both species. The level of differentiation discovered among D. arizonae populations was unexpectedly high and suggests that its populations, as well as those of D. mojavensis, may be themselves undergoing incipient speciation and merit further attention.     

Chromosome polymorphisms track trans‐Atlantic divergence and secondary contact in Atlantic salmon

Pleistocene glaciations drove repeated range contractions and expansions shaping contemporary intraspecific diversity. Atlantic salmon (Salmo salar) in the western and eastern Atlantic diverged >600,000 years before present, with the two lineages isolated in different southern refugia during glacial maxima, driving trans‐Atlantic genomic and karyotypic divergence. Here, we investigate the genomic consequences of glacial isolation and trans‐Atlantic secondary contact using 108,870 single nucleotide polymorphisms genotyped in 80 North American and European populations. Throughout North America, we identified extensive interindividual variation and discrete linkage blocks within and between chromosomes with known trans‐Atlantic differences in rearrangements: Ssa01/Ssa23 translocation and Ssa08/Ssa29 fusion. Spatial genetic analyses suggest independence of rearrangements, with Ssa01/Ssa23 showing high European introgression (>50%) in northern populations indicative of post‐glacial trans‐Atlantic secondary contact, contrasting with low European ancestry genome‐wide (3%). Ssa08/Ssa29 showed greater intrapopulation diversity, suggesting a derived chromosome fusion polymorphism that evolved within North America. Evidence of potential selection on both genomic regions suggests that the adaptive role of rearrangements warrants further investigation in Atlantic salmon. Our study highlights how Pleistocene glaciations can influence large‐scale intraspecific variation in genomic architecture of northern species.     

Genome sizes and chromosomes in the basal metazoan Hydra

Hydras belong to one of the earliest eumetazoan animal groups, but to date very little is known about their genome sizes, gene numbers, and chromosomes. Here we provide genome size estimates and corresponding karyotypes for five Hydra species. Nuclear DNA contents were assessed by slide-based Feulgen microphotometry. Hydra oligactis possesses the largest genome of 1450 Mbp, followed by similar 1 C capacities in H. carnea (1350 Mbp), H. vulgaris (1250 Mpb) and H. circumcincta (1150 Mbp). The smallest genome of 380 Mbp was determined in H. viridissima. While the number of chromosomes is identical in all five Hydra species (2n = 30), the size of the chromosomes is strictly correlated to the size of the genome, with H. viridissima having conspicuously small chromosomes. The taxonomic and evolutionary significance of the C-value and chromosomal size variation in this ancient group of metazoans as well as its impact on genomic organization and forthcoming genome projects are discussed.     

Genomic consequences of a recent three‐way admixture in supplemented wild brown trout populations revealed by local ancestry tracts

Understanding the evolutionary consequences of human‐mediated introductions of domesticated strains into the wild and their subsequent admixture with natural populations is of major concern in conservation biology. However, the genomic impacts of stocking from distinct sources (locally derived vs. divergent) on the genetic integrity of wild populations remain poorly understood. We designed an approach based on estimating local ancestry along individual chromosomes to provide a detailed picture of genomic admixture in supplemented populations. We used this approach to document admixture consequences in the brown trout Salmo trutta, for which decades of stocking practices have profoundly impacted the genetic make‐up of wild populations. In southern France, small local Mediterranean populations have been subject to successive introductions of domestic strains derived from the Atlantic and Mediterranean lineages. To address the impact of stocking, we evaluate the extent of admixture from both domestic strains within populations, using 75,684 mapped SNPs obtained from double‐digested restriction site‐associated DNA sequencing. Then, the chromosomal ancestry profiles of admixed individuals reveal a wider diversity of hybrid and introgressed genotypes than estimated using classical methods for inferring ancestry and hybrid pedigrees. In addition, the length distribution of introgressed tracts retained different timings of introgression between the two domestic strains. We finally reveal opposite consequences of admixture on the level of polymorphism of the recipient populations between domestic strains. Our study illustrates the potential of using the information contained in the genomic mosaic of ancestry tracts in combination with classical methods based on allele frequencies for analysing multiple‐way admixture with population genomic data.     

植物基因组研究与利用的新型工具——异源单体附加系

在高等植物中, 以种间杂交和回交把有益基因从一个物种转移到另一个物种为目的育种项目中, 单个外源染色体常常被附加到含有受体细胞完整一套染色体中, 形成异源单体附加系。这种异源单体附加系是阐明基因组结构和转移基因的有效工具。它可以通过回交形成覆盖整个基因组的渗入系重叠群, 用于建立以受体物种基因组为载体的外源物种基因组文库。另外, 一套完整的异源单体附加系也可看作是一个拥有分散供体基因组成为单个染色体单位的文库, 便于精确高通量地将标记分配到单个供体染色体上, 从而可以比较供体染色体和各自的直向同源受体染色体之间的标记位置和同线性关系。同时, 也便于研究同源染色体的渗入机制和配对状态。文中介绍了异源单体附加系的培育和特性, 并着重阐明了它在遗传育种和基础研究中的应用。     

Blocked recombination along the mating-type chromosomes of Neurospora tetrasperma involves both structural heterozygosity and autosomal genes

The Neurospora tetrasperma mating-type chromosomes have been shown to be structurally heterozygous by reciprocal introgression of these chromosomes between N. tetrasperma and N. crassa. This structural heterozygosity correlates with both a previously described recombination block and cytologically visible unpaired chromosomes at pachytene. Genes on the autosomes are also implicated in blocking recombination.     

Population genetics models of local ancestry

Migrations have played an important role in shaping the genetic diversity of human populations. Understanding genomic data thus requires careful modeling of historical gene flow. Here we consider the effect of relatively recent population structure and gene flow and interpret genomes of individuals that have ancestry from multiple source populations as mosaics of segments originating from each population. This article describes general and tractable models for local ancestry patterns with a focus on the length distribution of continuous ancestry tracts and the variance in total ancestry proportions among individuals. The models offer improved agreement with Wright-Fisher simulation data when compared to the state-of-the art and can be used to infer time-dependent migration rates from multiple populations. Considering HapMap African-American (ASW) data, we find that a model with two distinct phases of "European" gene flow significantly improves the modeling of both tract lengths and ancestry variances.     

Selective sorting of ancestral introgression in maize and teosinte along an elevational cline

While often deleterious, hybridization can also be a key source of genetic variation and pre-adapted haplotypes, enabling rapid evolution and niche expansion. Here we evaluate these opposing selection forces on introgressed ancestry between maize (Zea mays ssp. mays) and its wild teosinte relative, mexicana (Zea mays ssp. mexicana). Introgression from ecologically diverse teosinte may have facilitated maize’s global range expansion, in particular to challenging high elevation regions (> 1500 m). We generated low-coverage genome sequencing data for 348 maize and mexicana individuals to evaluate patterns of introgression in 14 sympatric population pairs, spanning the elevational range of mexicana, a teosinte endemic to the mountains of Mexico. While recent hybrids are commonly observed in sympatric populations and mexicana demonstrates fine-scale local adaptation, we find that the majority of mexicana ancestry tracts introgressed into maize over 1000 generations ago. This mexicana ancestry seems to have maintained much of its diversity and likely came from a common ancestral source, rather than contemporary sympatric populations, resulting in relatively low F_ST between mexicana ancestry tracts sampled from geographically distant maize populations.Introgressed mexicana ancestry in maize is reduced in lower-recombination rate quintiles of the genome and around domestication genes, consistent with pervasive selection against introgression. However, we also find mexicana ancestry increases across the sampled elevational gradient and that high introgression peaks are most commonly shared among high-elevation maize populations, consistent with introgression from mexicana facilitating adaptation to the highland environment. In the other direction, we find patterns consistent with adaptive and clinal introgression of maize ancestry into sympatric mexicana at many loci across the genome, suggesting that maize also contributes to adaptation in mexicana, especially at the lower end of its elevational range. In sympatric maize, in addition to high introgression regions we find many genomic regions where selection for local adaptation maintains steep gradients in introgressed mexicana ancestry across elevation, including at least two inversions: the well-characterized 14 Mb Inv4m on chromosome 4 and a novel 3 Mb inversion Inv9f surrounding the macrohairless1 locus on chromosome 9. Most outlier loci with high mexicana introgression show no signals of sweeps or local sourcing from sympatric populations and so likely represent ancestral introgression sorted by selection, resulting in correlated but distinct outcomes of introgression in different contemporary maize populations.     

Massive changes in genome architecture accompany the transition to self-fertility in the filamentous fungus Neurospora tetrasperma

A large region of suppressed recombination surrounds the sex-determining locus of the self-fertile fungus Neurospora tetrasperma. This region encompasses nearly one-fifth of the N. tetrasperma genome and suppression of recombination is necessary for self-fertility. The similarity of the N. tetrasperma mating chromosome to plant and animal sex chromosomes and its recent origin (<5 MYA), combined with a long history of genetic and cytological research, make this fungus an ideal model for studying the evolutionary consequences of suppressed recombination. Here we compare genome sequences from two N. tetrasperma strains of opposite mating type to determine whether structural rearrangements are associated with the nonrecombining region and to examine the effect of suppressed recombination for the evolution of the genes within it. We find a series of three inversions encompassing the majority of the region of suppressed recombination and provide evidence for two different types of rearrangement mechanisms: the recently proposed mechanism of inversion via staggered single-strand breaks as well as ectopic recombination between transposable elements. In addition, we show that the N. tetrasperma mat a mating-type region appears to be accumulating deleterious substitutions at a faster rate than the other mating type (mat A) and thus may be in the early stages of degeneration.     

Evaluation of the genomic extent of effects of fixed inversion differences on intraspecific variation and interspecific gene flow in Drosophila pseudoobscura and D. persimilis

There is increasing evidence that chromosomal inversions may facilitate the formation or persistence of new species by allowing genetic factors conferring species-specific adaptations or reproductive isolation to be inherited together and by reducing or eliminating introgression. However, the genomic domain of influence of the inverted regions on introgression has not been carefully studied. Here, we present a detailed study on the consequences that distance from inversion breakpoints has had on the inferred level of gene flow and divergence between Drosophila pseudoobscura and D. persimilis. We identified the locations of the inversion breakpoints distinguishing D. pseudoobscura and D. persimilis in chromosomes 2, XR, and XL. Population genetic data were collected at specific distances from the inversion breakpoints of the second chromosome and at two loci inside the XR and XL inverted regions. For loci outside the inverted regions, we found that distance from the nearest inversion breakpoint had a significant effect on several measures of divergence and gene flow between D. pseudoobscura and D. persimilis. The data fitted a logarithmic relationship, showing that the suppression of crossovers in inversion heterozygotes also extends to loci located outside the inversion but close to it (within 1-2 Mb). Further, we detected a significant reduction in nucleotide variation inside the inverted second chromosome region of D. persimilis and near one breakpoint, consistent with a scenario in which this inversion arose and was fixed in this species by natural selection.