Search

Abstract Growth hormone (GH) actions are mediated through binding to its cell-surface receptor, the GH receptor (GHR), with consequent activation of downstream signalling. However, nuclear GHR localisation has also been observed and is associated with increased cancer cell proliferation. Here we investigated the functional implications of nuclear translocation of the GHR in the human endometrial cancer cell-line, RL95-2, and human mammary epithelial cell-line, MCF-10A. We found that following GH treatment, the GHR rapidly translocates to the nucleus, with maximal localisation at 5–10 min. Combined immunoprecipitation-mass spectrometry analysis of RL95-2 whole cell lysates identified 40 novel GHR binding partners, including the transcriptional regulator, HMGN1. Moreover, microarray analysis demonstrated that the gene targets of HMGN1 were differentially expressed following GH treatment, and co-immunoprecipitation showed that HMGN1 associates with the GHR in the nucleus. Therefore, our results suggest that GHR nuclear translocation might mediate GH actions via interaction with chromatin factors that then drive changes in specific downstream transcriptional programs. Graphical abstract

SUMMARYRibosomal RNAs (rRNAs) are the most abundant cellular RNAs, and their synthesis from rDNA repeats by RNA Polymerase I accounts for the bulk of all transcription. Despite substantial variation in rRNA transcription rates across cell types, little is known about cell-type-specific factors that bind rDNA and regulate rRNA transcription to meet tissue-specific needs. Using hematopoiesis as a model system, we mapped about 2200 ChIP-Seq datasets for 250 transcription factors (TFs) and chromatin proteins to human and mouse rDNA, and identified robust binding of multiple TF families to canonical TF motifs on rDNA. Using a 47S-FISH-Flow assay developed for nascent rRNA quantification, we demonstrated that targeted degradation of CEBPA (C/EBP alpha), a critical hematopoietic TF with conserved rDNA binding, caused rapid reduction in rRNA transcription due to reduced Pol I occupancy. Our work identifies numerous potential rRNA regulators, and provides a template for dissection of TF roles in rRNA transcription.HIGHLIGHTSMultiple cell-type-specific transcription factors (TFs) bind canonical motifs on rDNA.The hematopoietic TF CEBPA binds to active rDNA alleles at a conserved site.CEBPA promotes Polymerase I occupancy and rRNA transcription in myeloid progenitors.We present ‘47S-FISH-Flow,’ a sensitive assay to quantify nascent rRNA.

Saccharomyces cerevisiae is an exceptional genetic system, with genetic crosses facilitated by its ability to be maintained in haploid and diploid forms. Such crosses are straightforward if the mating type/ploidy of the strains are known. Several techniques can be used to determine mating type (or ploidy), but all have limitations. Here we validate a simple, cheap and robust method to rapidly identify S. cerevisiae mating types. When cells of opposite mating type are mixed in liquid media, they “creep” up the culture vessel sides, a phenotype that can be easily detected visually. In contrast, mixtures of the same mating type or with a diploid simply settle out. The phenotype is robust to different media, cell densities, temperatures and strains, and is observable for several days. Microscopy suggests that cell aggregation during mating is responsible for the phenotype. Yeast knockout collection analysis identified 107 genes required for the creeping phenotype, with these being enriched for mating-specific genes. Surprisingly, the RIM101 signalling pathway was strongly represented. We propose that RIM101 signalling regulates aggregation as part of a wider, previously-unrecognized role in mating. The simplicity and robustness of this method makes it ideal for routine verification of S. cerevisiae mating type, with future studies required to verify its molecular basis.

Genome rearrangements in filamentous fungi are prevalent but little is known about the modalities of their evolution, in part because few complete genomes are available within a single genus. To address this, we have generated and compared 15 complete telomere-to-telomere genomes across the phylogeny of a single genus of filamentous fungi, Epichloë. We find that the strik-ing distinction between gene-rich and repeat-rich regions previously reported for isolated species is ubiquitous across the Epichloë genus. We built a species phylogeny from single-copy gene orthologs to provide a comparative framing to study chromosome composition and structural change through evolutionary time. All Epichloë genomes have exactly seven nuclear chromosomes, but despite this conserved ploidy, analyses reveal low synteny and substantial rearrangement of gene content across the genus. These rearrangements are highly lineage-dependent, with most oc-curring over short evolutionary distances, with long periods of structural stasis. Quantification of chromosomal rearrangements shows they are uncorrelated with numbers of substitutions and evolutionary distances, suggesting that different modes of evolution are acting to create nucleotide and chromosome-scale changes., Journal of Fungi, 8 (7), ISSN:2309-608X

Hybridization is a route to speciation that occurs widely across the eukaryote tree of life. The success of allopolyploids (hybrid species with increased ploidy) and homoploid hybrids (with unchanged ploidy) is well documented. However, their formation and establishment is not straightforward, with a suite of near-instantaneous and longer term biological repercussions faced by the new species. Central to these challenges is the rewiring of gene regulatory networks following the merger of distinct genomes inherited from both parental species. Research on the evolution of hybrid gene expression has largely involved studies on a single hybrid species or a few gene families. Here, we present the first standardized transcriptome-wide study exploring the fates of genes following hybridization across three kingdoms: animals, plants and fungi. Within each kingdom, we pair an allopolyploid system with a closely related homoploid hybrid to decouple the influence of increased ploidy from genome merger. Genome merger, not changes in ploidy, has the greatest effect on posthybridization expression patterns across all study systems. Strikingly, we find that differentially expressed genes in parent species preferentially switch to more similar expression in hybrids across all kingdoms, likely as a consequence of regulatory trans-acting cross-talk within the hybrid nucleus. We also highlight the prevalence of gene loss or silencing among extremely differentially expressed genes in hybrid species across all kingdoms. These shared patterns suggest that the evolutionary process of hybridization leads to common high-level expression outcomes, regardless of the particular species or kingdom.

This study investigated the impact of varying sound conditions (frequency and intensity) on yeast growth, fermentation performance and production of volatile organic compounds (VOCs) in beer. Fermentations were carried out in plastic bags suspended in large water-filled containers fitted with underwater speakers. Ferments were subjected to either 200–800 or 800–2000 Hz at 124 and 140 dB @ 20 µPa. Headspace solid-phase microextraction (HS-SPME) coupled with gas chromatography-mass spectrometry (GC-MS) was used to identify and measure the relative abundance of the VOCs produced. Sound treatment had significant effects on the number of viable yeast cells in suspension at 10 and 24 h (p < 0.05), with control (silence) samples having the highest cell numbers. For wort gravity, there were significant differences between treatments at 24 and 48 h, with the silence control showing the lowest density before all ferments converged to the same final gravity at 140 h. A total of 33 VOCs were identified in the beer samples, including twelve esters, nine alcohols, three acids, three aldehydes, and six hop-derived compounds. Only the abundance of some alcohols showed any consistent response to the sound treatments. These results show that the application of audible sound via underwater transmission to a beer fermentation elicited limited changes to wort gravity and VOCs during fermentation.

The biological effect of sound on microorganisms has been a field of interest for many years, with studies mostly focusing on ultrasonic and infrasonic vibrations. In the audible range (20 Hz to 20 kHz), sound has been shown to both increase colony formation and disrupt microbial growth, depending upon the organism and frequency of sound used. In the brewer’s yeast Saccharomyces cerevisiae, sound has been shown to significantly alter growth, increase alcohol production, and affect the metabolite profile. In this study, S. cerevisiae was exposed to a continuous 90 dB @ 20 μPa tone at different frequencies (0.1 kHz, 10 kHz, and silence). Fermentation characteristics were monitored over a 50-h fermentation in liquid malt extract, with a focus on growth rate and biomass yield. The profile of volatile metabolites at the subsequent stationary phase of the ferment was characterised by headspace gas chromatography–mass spectrometry. Sound treatments resulted in a 23% increase in growth rate compared to that of silence. Subsequent analysis showed significant differences in the volatilomes between all experimental conditions. Specifically, aroma compounds associated with citrus notes were upregulated with the application of sound. Furthermore, there was a pronounced difference in the metabolites produced in high- versus low-frequency sounds. This suggests industrial processes, such as beer brewing, could be modulated by the application of audible sound at specific frequencies during growth.

Ribosomal DNA genes (rDNA) encode the major ribosomal RNAs (rRNA) and in eukaryotic genomes are typically present as one or more arrays of tandem repeats. Species have characteristic rDNA copy numbers, ranging from tens to thousands of copies, with the number thought to be redundant for rRNA production. However, the tandem rDNA repeats are prone to recombination-mediated changes in copy number, resulting in substantial intra-species copy number variation. There is growing evidence that these copy number differences can have phenotypic consequences. However, we lack a comprehensive understanding of what determines rDNA copy number, how it evolves, and what the consequences are, in part because of difficulties in quantifying copy number. Here, we developed a genomic sequence read approach that estimates rDNA copy number from the modal coverage of the rDNA and whole genome to help overcome limitations in quantifying copy number with existing mean coverage-based approaches. We validated our method using strains of the yeast Saccharomyces cerevisiae with previously-determined rDNA copy numbers, and then applied our pipeline to investigate rDNA copy number in a global sample of 788 yeast isolates. We found that wild yeast have a mean copy number of 92, consistent with what is reported for other fungi but much lower than in laboratory strains. We also show that different populations have different rDNA copy numbers. These differences can partially be explained by phylogeny, but other factors such as environment are also likely to contribute to population differences in copy number. Our results demonstrate the utility of the modal coverage method, and highlight the high level of rDNA copy number variation within and between populations.Author summaryThe ribosomal RNA gene repeats (rDNA) form large tandem repeat arrays in most eukaryote genomes. Their tandem arrangement makes the rDNA prone to copy number variation, and there is increasing evidence that this copy number variation has phenotypic consequences. However, difficulties in measuring rDNA copy number hamper investigation into rDNA copy number dynamics and their significance. Here we developed a novel bioinformatics method for measuring rDNA copy number from whole genome sequence data that is based on the modal sequence read coverage. We established parameters for optimal performance of the method and validated it using yeast strains of known rDNA copy numbers. We then applied the method to a dataset of almost 800 global yeast isolates and demonstrate that yeast populations have different rDNA copy numbers that partially correlate with phylogeny. Our work provides a simple and accurate method for determining rDNA copy number that leverages the growing number of whole genome datasets, and highlights the dynamic nature of rDNA copy number.

The three-dimensional organization of the genome contributes to its maintenance and regulation. While chromosomal regions associate with nucleolar ribosomal RNA genes (rDNA), the biological significance of rDNA-genome interactions and whether they are dynamically regulated during disease remain unclear. rDNA chromatin exists in multiple inactive and active states and their transition is regulated by the RNA polymerase I transcription factor UBTF. Here, using a MYC-driven lymphoma model, we demonstrate that during malignant progression the rDNA chromatin converts to the open state, which is required for tumor cell survival. Moreover, this rDNA transition co-occurs with a reorganization of rDNA-genome contacts which correlate with gene expression changes at associated loci, impacting gene ontologies including B-cell differentiation, cell growth and metabolism. We propose that UBTF-mediated conversion to open rDNA chromatin during malignant transformation contributes to the regulation of specific gene pathways that regulate growth and differentiation through reformed long-range physical interactions with the rDNA., Jeannine Diesch et al. report the changes in rDNA chromatin state associated with cell transition into malignancy. They show that a specific transcription factor regulates this transition by altering rDNA chromatin, resulting in the reorganization of contacts between rDNA and the genome.

Endophytes of the genus Epichloë (Clavicipitaceae, Ascomycota) frequently occur within cool-season grasses and form interactions with their hosts that range from mutualistic to antagonistic. Many Epichloë species have arisen via interspecific hybridization, resulting in species with two or three subgenomes that retain all or nearly all of their original parental genomes, a process termed allopolyploidization. Here, we characterize Epichloë hybrida, sp. nov., a mutualistic species that has increasingly become a model system for investigating allopolyploidy in fungi. The Epichloë species so far identified as the closest known relatives of the two progenitors of E. hybrida are E. festucae var. lolii and E. typhina. We confirm that the nuclear genome of E. hybrida contains two homeologs of most protein-coding genes from E. festucae and E. typhina, with genome-wide gene expression analysis indicating a slight bias in overall gene expression from the E. typhina subgenome. Mitochondrial DNA is detectable only from E. festucae, whereas ribosomal DNA is detectable only from E. typhina. Inheriting ribosomal DNA from just one parent might be expected to preferentially favor interactions with ribosomal proteins from the same parent, but we find that ribosomal protein genes from both parental subgenomes are nearly all expressed equally in E. hybrida. Finally, we provide a comprehensive set of resources for this model system that are intended to facilitate further study of fungal hybridization by other researchers.

Koonin argues that CRISPR-Cas systems present the best-known case in point for Lamarckian evolution because they satisfy his proposed criteria for the specific inheritance of acquired adaptive characteristics (IAC). We see two interrelated issues with Koonin’s characterization of CRISPR-Cas systems as Lamarckian. First, at times he appears to confuse an account of the CRISPR-Cas system with an account of the mechanism it employs. We argue there is no evidence for the CRISPR-Cas system being “Lamarckian” in any sense. Second, it is unclear whether the mechanism is more “Lamarckian” than many other forms of genetic change already well-characterized in Darwinian terms. We present three conceptually distinct senses in which the mechanism of IAC may be considered Lamarckian and argue that only the strongest sense of goal-directed IAC would be difficult to accommodate in a Darwinian account. As the CRISPR-Cas mechanism does not qualify as “Lamarckian” in this strong sense, we argue there is no conceptual value in calling it “Lamarckian”. Finally, we suggest that CRISPR-Cas systems do hold the potential for genuinely non-Darwinian, directed evolution in a way that Koonin did not discuss, involving their potential (mis)use as a human gene-editing tool.

The advent of 'next generation sequencing' (NGS) technologies has led to the discovery of many novel mycoviruses, the majority of which are sufficiently different from previously sequenced viruses that there is no appropriate reference sequence on which to base the sequence assembly. Although many new genome sequences are generated by NGS, confirmation of the sequence by Sanger sequencing is still essential for formal classification by the International Committee for the Taxonomy of Viruses (ICTV), although this is currently under review. To empirically test the validity of de novo assembled mycovirus genomes from dsRNA extracts, we compared the results from Illumina sequencing with those from random cloning plus targeted PCR coupled with Sanger sequencing for viruses from five Sclerotinia sclerotiorum isolates. Through Sanger sequencing we detected nine viral genomes while through Illumina sequencing we detected the same nine viruses plus one additional virus from the same samples. Critically, the Illumina derived sequences share >99.3 % identity to those obtained by cloning and Sanger sequencing. Although, there is scope for errors in de novo assembled viral genomes, our results demonstrate that by maximising the proportion of viral sequence in the data and using sufficiently rigorous quality controls, it is possible to generate de novo genome sequences of comparable accuracy from Illumina sequencing to those obtained by Sanger sequencing.

Structural features of genomes, including the three-dimensional arrangement of DNA in the nucleus, are increasingly seen as key contributors to the regulation of gene expression. However, studies on how genome structure and nuclear organisation influence transcription have so far been limited to a handful of model species. This narrow focus limits our ability to draw general conclusions about the ways in which three-dimensional structures are encoded, and to integrate information from three-dimensional data to address a broader gamut of biological questions. Here, we generate a complete and gapless genome sequence for the filamentous fungus, Epichloë festucae. We use Hi-C data to examine the three-dimensional organisation of the genome, and RNA-seq data to investigate how Epichloë genome structure contributes to the suite of transcriptional changes needed to maintain symbiotic relationships with the grass host. Our results reveal a genome in which very repeat-rich blocks of DNA with discrete boundaries are interspersed by gene-rich sequences that are almost repeat-free. In contrast to other species reported to date, the three-dimensional structure of the genome is anchored by these repeat blocks, which act to isolate transcription in neighbouring gene-rich regions. Genes that are differentially expressed in planta are enriched near the boundaries of these repeat-rich blocks, suggesting that their three-dimensional orientation partly encodes and regulates the symbiotic relationship formed by this organism., Author summary How a genome folds within a nucleus can contribute to the regulation of gene expression. There is now considerable interest in understanding how the three-dimensional structure of a genome is maintained and how this structure influences biological processes. However, most studies on these phenomena focus on a few model organisms, which limits our knowledge of how the nuclear organisation is maintained and the range of biological processes to which it contributes. In this paper, we examine the three-dimensional organisation of the genome in Epichloë festucae, a symbiotic fungus. Our results demonstrate that repetitive elements play a key role in shaping the nuclear organisation of E. festucae. These repetitive elements help to divide the genome into distinct regions that have similar gene expression profiles. Our results demonstrate a novel role for repeats in maintaining a three-dimensional organisation that contributes to regulating the symbiotic relationship formed by this species.

The spatial organization of eukaryotic genomes is linked to their functions. However, how individual features of the global spatial structure contribute to nuclear function remains largely unknown. We previously identified a high-frequency interchromosomal interaction within the Saccharomyces cerevisiae genome that occurs between the intergenic spacer of the ribosomal DNA (rDNA) repeats and the intergenic sequence between the locus encoding the second largest RNA polymerase I subunit and a lysine tRNA gene [i.e., RPA135-tK(CUU)P]. Here, we used quantitative chromosome conformation capture in combination with replacement mapping to identify a 75-bp sequence within the RPA135-tK(CUU)P intergenic region that is involved in the interaction. We demonstrate that the RPA135-IGS1 interaction is dependent on the rDNA copy number and the Msn2 protein. Surprisingly, we found that the interaction does not govern RPA135 transcription. Instead, replacement of a 605-bp region within the RPA135-tK(CUU)P intergenic region results in a reduction in the RPA135-IGS1 interaction level and fluctuations in rDNA copy number. We conclude that the chromosomal interaction that occurs between the RPA135-tK(CUU)P and rDNA IGS1 loci stabilizes rDNA repeat number and contributes to the maintenance of nucleolar stability. Our results provide evidence that the DNA loci involved in chromosomal interactions are composite elements, sections of which function in stabilizing the interaction or mediating a functional outcome.

The yeast Saccharomyces cerevisiae is a powerful model system to understand the molecular basis of aging. It has been known for over 50 years that yeast cells have a finite replicative capacity and develop an aging phenotype, and much recent research has focused on the molecular changes that underlie this replicative senescence. A common thread in many yeast replicative aging studies is the involvement of the ribosomal DNA gene repeats (rDNA), beginning with the discovery that the rDNA silencing gene, SIR2, regulates life span. In 2008, a novel aging hypothesis, termed the rDNA theory of aging, was presented where the high level of genomic instability at the rDNA repeats was proposed to dominate global genome stability and determine the life span. Here, we review the rDNA theory of aging and discuss a number of recent studies that provide important new data on the roles of the rDNA in yeast replicative aging. Based on these recent results, we propose an integrative model of the rDNA theory of aging that encompasses genomic instability, chromatin relocalization following DNA repair, and replication stress in a self-reinforcing cyclical pathway that is primarily manifested at the rDNA repeats and results in the aging phenotype.

SummaryGenomic instability is a conserved factor in lifespan reduction, although the molecular mechanism is not known [1]. Studies in the yeast Saccharomyces cerevisiae over the past 20 years have found a connection between the ribosomal RNA gene cluster (rDNA) and lifespan [2]. The highly repetitive rDNA exhibits genomic instability, and the antiaging histone deacetylase gene SIR2 regulates this instability. We previously proposed that SIR2 governs lifespan by repressing rDNA noncoding transcription and rDNA instability, but the extent to which lifespan is affected by SIR2 acting at the rDNA versus other genomic regions, and the relationship between rDNA noncoding transcription/rDNA stability and lifespan have remained controversial. To control rDNA noncoding transcription and rDNA instability, we use a strain in which the rDNA noncoding promoter is replaced with an inducible promoter. Here, we show that repression of noncoding transcription extends lifespan and makes SIR2 dispensable for lifespan extension. These results indicate that Sir2 maintains lifespan through repression of E-pro noncoding transcription in the rDNA cluster, rather than pleiotropically at other loci. The observation of rDNA instability in other organisms, including humans, suggests that this may be a conserved aging pathway.

The ribosomal RNA genes (rDNA) encode the major rRNA species of the ribosome, and thus are essential across life. These genes are highly repetitive in most eukaryotes, forming blocks of tandem repeats that form the core of nucleoli. The primary role of the rDNA in encoding rRNA has been long understood, but more recently the rDNA has been implicated in a number of other important biological phenomena, including genome stability, cell cycle, and epigenetic silencing. Noncoding elements, primarily located in the intergenic spacer region, appear to mediate many of these phenomena. Although sequence information is available for the genomes of many organisms, in almost all cases rDNA repeat sequences are lacking, primarily due to problems in assembling these intriguing regions during whole genome assemblies. Here, we present a method to obtain complete rDNA repeat unit sequences from whole genome assemblies. Limitations of next generation sequencing (NGS) data make them unsuitable for assembling complete rDNA unit sequences; therefore, the method we present relies on the use of Sanger whole genome sequence data. Our method makes use of the Arachne assembler, which can assemble highly repetitive regions such as the rDNA in a memory-efficient way. We provide a detailed step-by-step protocol for generating rDNA sequences from whole genome Sanger sequence data using Arachne, for refining complete rDNA unit sequences, and for validating the sequences obtained. In principle, our method will work for any species where the rDNA is organized into tandem repeats. This will help researchers working on species without a complete rDNA sequence, those working on evolutionary aspects of the rDNA, and those interested in conducting phylogenetic footprinting studies with the rDNA.

The ribosomal RNA genes (rDNA) encode the major rRNA species of the ribosome, and thus are essential across life. These genes are highly repetitive in most eukaryotes, forming blocks of tandem repeats that form the core of nucleoli. The primary role of the rDNA in encoding rRNA has been long understood, but more recently the rDNA has been implicated in a number of other important biological phenomena, including genome stability, cell cycle, and epigenetic silencing. Noncoding elements, primarily located in the intergenic spacer region, appear to mediate many of these phenomena. Although sequence information is available for the genomes of many organisms, in almost all cases rDNA repeat sequences are lacking, primarily due to problems in assembling these intriguing regions during whole genome assemblies. Here, we present a method to obtain complete rDNA repeat unit sequences from whole genome assemblies. Limitations of next generation sequencing (NGS) data make them unsuitable for assembling complete rDNA unit sequences; therefore, the method we present relies on the use of Sanger whole genome sequence data. Our method makes use of the Arachne assembler, which can assemble highly repetitive regions such as the rDNA in a memory-efficient way. We provide a detailed step-by-step protocol for generating rDNA sequences from whole genome Sanger sequence data using Arachne, for refining complete rDNA unit sequences, and for validating the sequences obtained. In principle, our method will work for any species where the rDNA is organized into tandem repeats. This will help researchers working on species without a complete rDNA sequence, those working on evolutionary aspects of the rDNA, and those interested in conducting phylogenetic footprinting studies with the rDNA.

Polyploidy-the increase in the number of whole chromosome sets-is an important evolutionary force in eukaryotes. Polyploidy is well recognized throughout the evolutionary history of plants and animals, where several ancient events have been hypothesized to be drivers of major evolutionary radiations. However, fungi provide a striking contrast: while numerous recent polyploids have been documented, ancient fungal polyploidy is virtually unknown. We present a survey of known fungal polyploids that confirms the absence of ancient fungal polyploidy events. Three hypotheses may explain this finding. First, ancient fungal polyploids are indeed rare, with unique aspects of fungal biology providing similar benefits without genome duplication. Second, fungal polyploids are not successful in the long term, leading to few extant species derived from ancient polyploidy events. Third, ancient fungal polyploids are difficult to detect, causing the real contribution of polyploidy to fungal evolution to be underappreciated. We consider each of these hypotheses in turn and propose that failure to detect ancient events is the most likely reason for the lack of observed ancient fungal polyploids. We examine whether existing data can provide evidence for previously unrecognized ancient fungal polyploidy events but discover that current resources are too limited. We contend that establishing whether unrecognized ancient fungal polyploidy events exist is important to ascertain whether polyploidy has played a key role in the evolution of the extensive complexity and diversity observed in fungi today and, thus, whether polyploidy is a driver of evolutionary diversifications across eukaryotes. Therefore, we conclude by suggesting ways to test the hypothesis that there are unrecognized polyploidy events in the deep evolutionary history of the fungi.

TAR1 (transcript antisense to ribosomal RNA) is a young gene, located antisense to the 25S rRNA gene in Saccharomyces cerevisiae 1. The ribosomal DNA (rDNA) exists as ∼150 tandem repeats 2, making TAR1 the most abundant protein-coding gene in yeast. Oddly, TAR1 is normally silenced by Sir2p, a repressor of RNA polymerase II (pol-II)-transcribed genes 3. Recent reports suggest Tar1p protein is localised to the inner mitochondrial membrane 4, interacts with Coq5p (a protein involved in coenzyme Q synthesis [5]), and can maintain oxidative phosphorylation capacity 5. Direct elucidation of TAR1 function is lacking however, as available observations derive from monitoring a single, modified TAR1 copy 5. This may not be representative of the majority of genomic copies, and some of these results are in conflict with previous reports on rDNA Pol-II transcript expression 3, 6, 7. Given difficulties in probing the function of a multi-copy antisense gene, we examine available data in an effort to better understand the role of TAR1. We propose that TAR1 ameliorates the behaviour of selfish yeast mitochondrial mutants first identified over fifty years ago. The location of TAR1 in the ribosomal DNA (rDNA) repeat array is crucial to our model, as this means it is also present on extra-chromosomal ribosomal circles (ERCs). ERCs are rDNA repeats that have ‘popped out’ of the chromosome by intra-chromatid recombination, and exist in the cell as plasmid-like circular DNA (Fig. 1). ERC generation is thought to curtail replicative (though not chronological) lifespan 8, 9. However, as well as accumulating in old cells, they also accumulate in yeast with defective mitochondria 10. Figure 1 Relationship between the retrograde response and rDNA. A: Simplified view of the retrograde response. The retrograde response enables yeast cells to detect and respond to mitochondrial damage. In brief, some signal, possibly a drop in mitochondrial membrane ... Significantly, ERC accumulation follows activation of the retrograde response pathway in yeast, upon mitochondrial dysfunction (Fig. 1) 11, 12. This enables survival despite diminished respiration capacity. The extent of the response corresponds to the level of mitochondrial dysfunction 13. The downstream effect of the retrograde response is upregulation of mitochondrial damage, nuclear-encoded metabolic, and stress response genes, enabling yeast to grow on fermentable carbon sources 11, 12. It also extends lifespan. Our proposal resolves the seemingly paradoxical outcomes of the retrograde response: on one hand it extends lifespan, yet it also generates lifespan-shortening ERCs that have no known role in this response 11, 12. These lifespan effects are paradoxical only if ERCs serve solely as senescence factors 8. We propose that ERCs have a positive function: suppressing the ability of selfish mitochondrial mutants to overrun populations of sexually reproducing yeast by upregulating TAR1 expression. Available experimental evidence supports this interpretation, which, if correct, indicates that the effect of the retrograde response on lifespan is but a side effect of ERC production, the primary aim of which is preventing the spread of respiration-deficient mitochondria.

Concerted evolution describes the unusual evolutionary pattern exhibited by certain repetitive sequences, whereby all the repeats are maintained in the genome with very similar sequences but differ between related species. The pattern of concerted evolution is thought to result from continual turnover of repeats by recombination, a process known as homogenization. Approaches to studying concerted evolution have largely been observational because of difficulties investigating repeat evolution in an experimental setting with large arrays of identical repeats. Here, we establish an experimental evolution approach to look at the rate and dynamics of concerted evolution in the ribosomal DNA (rDNA) repeats. A small targeted mutation was made in the spacer of a single rDNA unit in Saccharomyces cerevisiae so we could monitor the fate of this unit without the need for a selectable marker. The rate of loss of this single unit was determined, and the frequency of duplication was also estimated. The results show that duplication and deletion events occur at similar rates and are very common: An rDNA unit may be gained or lost as frequently as once every cell division. Investigation of the spatial dynamics of rDNA turnover showed that when the tagged repeat unit was duplicated, the copy predominantly, but not exclusively, ended up near to the tagged repeat. This suggests that variants in the rDNA spread in a semiclustered fashion. Surprisingly, large deletions that remove a significant fraction of total rDNA repeats were frequently found. We propose these large deletions are a driving force of concerted evolution, acting to increase homogenization efficiency over-and-above that afforded by turnover of individual rDNA units. Thus, the results presented here enhance our understanding of concerted evolution by offering insights into both the spatial and temporal dynamics of the homogenization process and suggest an important new aspect in our understanding of concerted evolution.

In eukaryotes, the ribosomal DNA (rDNA) consists of long tandem repeat arrays. These repeated genes are unstable because homologous recombination between them results in copy number loss. To maintain high copy numbers, yeast has an amplification system that works through a pathway involving the replication fork barrier site and unequal sister chromatid recombination. In this study, we show that an active replication origin is essential for amplification, and the amplification rate correlates with origin activity. Moreover, origin activity affects the levels of extrachromosomal rDNA circles (ERC) that are thought to promote aging. Surprisingly, we found that reduction in ERC level results in shorter life span. We instead show that life span correlates with rDNA stability, which is preferentially reduced in mother cells, and that episomes can induce rDNA instability. These data support a model in which rDNA instability itself is a cause of aging in yeast.

Repeat families within genomes are often maintained with similar sequences. Traditionally, this has been explained by concerted evolution, where repeats in an array evolve “in concert” with the same sequence via continual turnover of repeats by recombination. Another form of evolution, birth-and-death evolution, can also explain this pattern, although in this case selection is the critical force maintaining the repeats. The level of intragenomic variation is the key difference between these two forms of evolution. The prohibitive size and repetitive nature of large repeat arrays have made determination of the absolute level of intragenomic repeat variability difficult, thus there is little evidence to support concerted evolution over birth-and-death evolution for many large repeat arrays. Here we use whole-genome shotgun sequence data from the genome projects of five fungal species to reveal absolute levels of sequence variation within the ribosomal RNA gene repeats (rDNA). The level of sequence variation is remarkably low. Furthermore, the polymorphisms that are detected are not functionally constrained and seem to exist beneath the level of selection. These results suggest the rDNA is evolving via concerted evolution. Comparisons with a repeat array undergoing birth-and-death evolution provide a clear contrast in the level of repeat array variation between these two forms of evolution, confirming that the rDNA indeed does evolve via concerted evolution. These low levels of intra-genomic variation are consistent with a model of concerted evolution in which homogenization is very rapid and efficiently maintains highly similar repeat arrays.

This paper proposes a novel approach to comparing cell colony images taken at different times on a Petri dish. The objective is to provide an assistive tool for microbiologists to quantify the loss of cell colonies on two Petri dishes, by mapping cell colonies between a pair of images. This problem is highly non-trivial, as the shape, size and position of the corresponding colonies vary randomly. In addition, the cell colony images for comparison are taken at different times and from slightly different perspectives (i.e. effects of shearing), and this amplifies the complexity of the problem. Experiments show that approaches purely based on SIFT or SURF, or algorithms used in astronomy, do not perform well on the problem domain. We therefore introduce a new approach to addressing these problems. A novel iterative technique that combines triangulation algorithms with the RANSAC alignment algorithm is proposed. Through hundreds of experiments, we demonstrate the efficacy of the new algorithm in comparison to existing ones found in the literature.

Sequences involved in the regulation of genetic and genomic processes are primarily located in noncoding regions. Identifying such cis -acting sequences from sequence data is difficult because their patterns are not readily apparent, and, to date, identification has concentrated on regions associated with gene-coding functions. Here, we used phylogenetic footprinting to look for gene-independent noncoding elements in the ribosomal RNA gene repeats from several Saccharomyces species. Similarity plots of ribosomal intergenic spacer alignments from six closely related Saccharomyces species allowed the identification of previously characterized functional elements, such as the origin-of-replication and replication-fork barrier sites, demonstrating that this method is a powerful predictor of noncoding functional elements. Seventeen previously uncharacterized elements, showing high levels of conservation, were also discovered. The conservation of these elements suggests that they are functional, and we demonstrate the functionality of two classes of these elements, a putative bidirectional noncoding promoter and a series of conserved peaks with matches to the origin-of-replication core consensus. Our results paint a comprehensive picture of the functionality of the Saccharomyces ribosomal intergenic region and demonstrate that functional elements not involved in gene-coding function can be identified by using comparative genomics based on sequence conservation.

We examined ribosomal RNA concerted evolution in an Epichloë endophyte interspecific hybrid (Lp1) and its progenitors (Lp5 and E8). We show that the 5S rrn genes are organized as dispersed copies. Cloned 5S gene sequences revealed two subfamilies exhibiting 12% sequence divergence, with substitutions forming coevolving pairs that maintain secondary structure and presumably function. Observed sequence patterns are not fully consistent with either concerted or classical evolution. The 5S rrn genes are syntenic with the tandemly arranged rDNA genes, despite residing outside the rDNA arrays. We also examined rDNA concerted evolution. Lp1 has rDNA sequence from only one progenitor and contains multiple rDNA arrays. Using 5S rrn genes as chromosomal markers, we propose that interlocus homogenization has replaced all Lp5 rDNA sequence with E8 sequence in the hybrid. This interlocus homogenization appears to have been rapid and efficient and is the first demonstration of hybrid interlocus homogenization in the Fungi.

Background Forming a new species through the merger of two or more divergent parent species is increasingly seen as a key phenomenon in the evolution of many biological systems. However, little is known about how expression of parental gene copies (homeologs) responds following genome merger. High throughput RNA sequencing now makes this analysis technically feasible, but tools to determine homeolog expression are still in their infancy. Results Here we present HyLiTE – a single-step analysis to obtain tables of homeolog expression in a hybrid or allopolyploid and its parent species directly from raw mRNA sequence files. By implementing on-the-fly detection of diagnostic parental polymorphisms, HyLiTE can perform SNP calling and read classification simultaneously, thus allowing HyLiTE to be run as parallelized code. HyLiTE accommodates any number of parent species, multiple data sources (including genomic DNA reads to improve SNP detection), and implements a statistical framework optimized for genes with low to moderate expression. Conclusions HyLiTE is a flexible and easy-to-use program designed for bench biologists to explore patterns of gene expression following genome merger. HyLiTE offers practical advantages over manual methods and existing programs, has been designed to accommodate a wide range of genome merger systems, can identify SNPs that arose following genome merger, and offers accurate performance on non-model organisms. Electronic supplementary material The online version of this article (doi:10.1186/s12859-014-0433-8) contains supplementary material, which is available to authorized users.

Why study cell aging and cell death, in yeast? Most of us have heard these questions before, sometimes from fellow colleagues. Reports from National Institute of Aging (National Institutes of Health), the World Health Organization, and other relevant organizations, outline the substantial challenges related to aging that we will face in the next decade. This concern is also reflected in Horizon2020, the next EU framework for research and innovation, where the challenges of an aging population and the consequences it brings for society, are among the top priorities. One of the key areas of concern is the economic and societal impact of age-related …

The short arms of the five acrocentric human chromosomes harbor sequences that direct the assembly and function of the nucleolus, one of the key functional domains of the nucleus, yet they are absent from the current human genome assembly. Here we describe the genomic architecture of these human nucleolar organizers. Sequences distal and proximal to ribosomal gene arrays are conserved among the acrocentric chromosomes, suggesting they are sites of frequent recombination. Although previously believed to be heterochromatic, characterization of these two flanking regions reveals that they share a complex genomic architecture similar to other euchromatic regions of the genome, but they have distinct genomic characteristics. Proximal sequences are almost entirely segmentally duplicated, similar to the regions bordering centromeres. In contrast, the distal sequence is predominantly unique to the acrocentric short arms and is dominated by a very large inverted repeat. We show that the distal element is localized to the periphery of the nucleolus, where it appears to anchor the ribosomal gene repeats. This, combined with its complex chromatin structure and transcriptional activity, suggests that this region is involved in nucleolar organization. Our results provide a platform for investigating the role of NORs in nucleolar formation and function, and open the door for determining the role of these regions in the well-known empirical association of nucleoli with pathology.

The nucleolus is a prominent nuclear structure that is the site of ribosomal RNA (rRNA) transcription, and hence ribosome biogenesis. Cellular demand for ribo-somes, and hence rRNA, is tightly linked to cell growth and the rRNA makes up the majority of all the RNA within a cell. To fulfill the cellular demand for rRNA, the ribo somal RNA (rDNA) genes are amplified to high copy number and transcribed at very high rates. As such, understanding the rDNA has profound consequences for our compre-hension of genome and transcriptional organization in cells. In this review, we address the question of whether the nucleolus is a raft adrift the sea of nuclear DNA, or actively contributes to genome organization. We present evidence supporting the idea that the nucleolus, and the rDNA contained therein, play more roles in the biology of the cell than simply ribosome biogenesis. We propose that the nucleolus and the rDNA are central factors in the spa-tial organization of the genome, and that rapid alterations in nucleolar structure in response to changing conditions manifest themselves in altered genomic structures that have functional consequences. Finally, we discuss some predictions that result from the nucleolus having a central role in nuclear organization. Keywords: genome architecture; nucleolus; rDNA.

Summary Plant pathogens use a complex arsenal of weapons, such as toxic secondary metabolites, to invade and destroy their hosts. Knowledge of how secondary metabolite pathways evolved is central to understanding the evolution of host specificity. The secondary metabolite dothistromin is structurally similar to aflatoxins and is produced by the fungal pine pathogen Dothistroma septosporum. Our study focused on dothistromin genes, which are widely dispersed across one chromosome, to determine whether this unusual distributed arrangement evolved from an ancestral cluster. We combined comparative genomics and population genetics approaches to elucidate the origins of the dispersed arrangement of dothistromin genes over a broad evolutionary timescale at the phylum, class and species levels. Orthologs of dothistromin genes were found in two major classes of fungi. Their organization is consistent with clustering of core pathway genes in a common ancestor, but with intermediate cluster fragmentation states in the Dothideomycetes fungi. Recombination hotspots in a D. septosporum population matched sites of gene acquisition and cluster fragmentation at higher evolutionary levels. The results suggest that fragmentation of a larger ancestral cluster gave rise to the arrangement seen in D. septosporum. We propose that cluster fragmentation may facilitate metabolic retooling and subsequent host adaptation of plant pathogens.

We sequenced and compared the genomes of the Dothideomycete fungal plant pathogens Cladosporium fulvum (Cfu) (syn. Passalora fulva) and Dothistroma septosporum (Dse) that are closely related phylogenetically, but have different lifestyles and hosts. Although both fungi grow extracellularly in close contact with host mesophyll cells, Cfu is a biotroph infecting tomato, while Dse is a hemibiotroph infecting pine. The genomes of these fungi have a similar set of genes (70% of gene content in both genomes are homologs), but differ significantly in size (Cfu >61.1-Mb; Dse 31.2-Mb), which is mainly due to the difference in repeat content (47.2% in Cfu versus 3.2% in Dse). Recent adaptation to different lifestyles and hosts is suggested by diverged sets of genes. Cfu contains an α-tomatinase gene that we predict might be required for detoxification of tomatine, while this gene is absent in Dse. Many genes encoding secreted proteins are unique to each species and the repeat-rich areas in Cfu are enriched for these species-specific genes. In contrast, conserved genes suggest common host ancestry. Homologs of Cfu effector genes, including Ecp2 and Avr4, are present in Dse and induce a Cf-Ecp2- and Cf-4-mediated hypersensitive response, respectively. Strikingly, genes involved in production of the toxin dothistromin, a likely virulence factor for Dse, are conserved in Cfu, but their expression differs markedly with essentially no expression by Cfu in planta. Likewise, Cfu has a carbohydrate-degrading enzyme catalog that is more similar to that of necrotrophs or hemibiotrophs and a larger pectinolytic gene arsenal than Dse, but many of these genes are not expressed in planta or are pseudogenized. Overall, comparison of their genomes suggests that these closely related plant pathogens had a common ancestral host but since adapted to different hosts and lifestyles by a combination of differentiated gene content, pseudogenization, and gene regulation., Author Summary We compared the genomes of two closely related pathogens with very different lifestyles and hosts: C. fulvum (Cfu), a biotroph of tomato, and D. septosporum (Dse), a hemibiotroph of pine. Some differences in gene content were identified that can be directly related to their different hosts, such as the presence of a gene involved in degradation of a tomato saponin only in Cfu. However, in general the two species share a surprisingly large proportion of genes. Dse has functional homologs of Cfu effector genes, while Cfu has genes for biosynthesis of dothistromin, a toxin probably associated with virulence in Dse. Cfu also has an unexpectedly large content of genes for biosynthesis of other secondary metabolites and degradation of plant cell walls compared to Dse, contrasting with its host preference and lifestyle. However, many of these genes were not expressed in planta or were pseudogenized. These results suggest that evolving species may retain genetic signatures of the host and lifestyle preferences of their ancestor and that evolution of new genes, gene regulation, and pseudogenization are important factors in adaptation.

Bilinski et al . (2012) recently presented a provocative commentary regarding yeast longevity and hypertrophy in this journal. We are compelled, as yeast aging researchers, to present our thinking on some of the critiques presented by Bilinski et al . First, Bilinski et al . argued against the presence of an aging factor because there is little evidence for a ‘soluble, diffusible, and degradable’ factor promoting aging. While perhaps originally defined as such, we consider an ‘aging factor’ to be something that incrementally attenuates lifespan, accumulates in mothers, and is (largely) not passed onto daughters. Thus, yeast hypertrophy, as proposed by Bilinski et al ., fits this definition. Bilinski et al . also note that yeast hypertrophy is usually ignored in articles proposing putative aging factors. While an explicit link between hypertrophy and these aging …

Nadine Hein1,*, Elaine Sanij1,2,*, Jaclyn Quin1,3, Katherine M. Hannan1, Austen Ganley4,# and Ross D. Hannan1,3,5,6,# 1Division of Cancer Research, Peter MacCallum Cancer Centre, Melbourne, 2Department of Pathology, University of Melbourne, Melbourne, 3Department of Biochemistry and Molecular Biology, University of Melbourne, Melbourne, 4Institute of Natural Sciences, Massey University, Auckland, 5Department of Biochemistry and Molecular Biology, Monash University, Melbourne, 6School of Biomedical Sciences, The University of Queensland, Melbourne, 1,2,3,5,6Australia 4New Zealand

This study attempts to resolve relationships among and within the four basal arthropod lineages (Pancrustacea, Myriapoda, Euchelicerata, Pycnogonida) and to assess the widespread expectation that remaining phylogenetic problems will yield to increasing amounts of sequence data. Sixty-eight regions of 62 protein-coding nuclear genes (approximately 41 kilobases (kb)/taxon) were sequenced for 12 taxonomically diverse arthropod taxa and a tardigrade outgroup. Parsimony, likelihood, and Bayesian analyses of total nucleotide data generally strongly supported the monophyly of each of the basal lineages represented by more than one species. Other relationships within the Arthropoda were also supported, with support levels depending on method of analysis and inclusion/exclusion of synonymous changes. Removing third codon positions, where the assumption of base compositional homogeneity was rejected, altered the results. Removing the final class of synonymous mutations--first codon positions encoding leucine and arginine, which were also compositionally heterogeneous--yielded a data set that was consistent with a hypothesis of base compositional homogeneity. Furthermore, under such a data-exclusion regime, all 68 gene regions individually were consistent with base compositional homogeneity. Restricting likelihood analyses to nonsynonymous change recovered trees with strong support for the basal lineages but not for other groups that were variably supported with more inclusive data sets. In a further effort to increase phylogenetic signal, three types of data exploration were undertaken. (1) Individual genes were ranked by their average rate of nonsynonymous change, and three rate categories were assigned--fast, intermediate, and slow. Then, bootstrap analysis of each gene was performed separately to see which taxonomic groups received strong support. Five taxonomic groups were strongly supported independently by two or more genes, and these genes mostly belonged to the slow or intermediate categories, whereas groups supported only by a single gene region tended to be from genes of the fast category, arguing that fast genes provide a less consistent signal. (2) A sensitivity analysis was performed in which increasing numbers of genes were excluded, beginning with the fastest. The number of strongly supported nodes increased up to a point and then decreased slightly. Recovery of Hexapoda required removal of fast genes. Support for Mandibulata (Pancrustacea + Myriapoda) also increased, at times to "strong" levels, with removal of the fastest genes. (3) Concordance selection was evaluated by clustering genes according to their ability to recover Pancrustacea, Euchelicerata, or Myriapoda and analyzing the three clusters separately. All clusters of genes recovered the three concordance clades but were at times inconsistent in the relationships recovered among and within these clades, a result that indicates that the a priori concordance criteria may bias phylogenetic signal in unexpected ways. In a further attempt to increase support of taxonomic relationships, sequence data from 49 additional taxa for three slow genes (i.e., EF-1 alpha, EF-2, and Pol II) were combined with the various 13-taxon data sets. The 62-taxon analyses supported the results of the 13-taxon analyses and provided increased support for additional pancrustacean clades found in an earlier analysis including only EF-1 alpha, EF-2, and Pol II.

Phylogenetic footprinting is powerful technique for finding functional elements from sequence data. Functional elements are thought to have greater sequence constraint than nonfunctional elements, and, thus, undergo a slower rate of sequence change through time. Phylogenetic footprinting uses comparisons of homologous sequences from closely related organisms to identify "phylogenetic footprints," regions with slower rates of sequence change than background. This does not require prior characterization of the sequence in question, therefore, it can be used in a wide range of applications. In particular, it is useful for the identification of functional elements in noncoding DNA, which are traditionally difficult to detect. Here, we describe in detail how to perform a simple yet powerful phylogenetic footprinting analysis. As an example, we use ribosomal DNA repeat sequences from various Saccharomyces yeasts to find functional noncoding DNA elements in the intergenic spacer, and explain critical considerations in performing phylogenetic footprinting analyses, including the number of species and species range, and some of the available software. Our methods are broadly applicable and should appeal to molecular biologists with little experience in bioinformatics.

Phylogenetic footprinting is powerful technique for finding functional elements from sequence data. Functional elements are thought to have greater sequence constraint than nonfunctional elements, and, thus, undergo a slower rate of sequence change through time. Phylogenetic footprinting uses comparisons of homologous sequences from closely related organisms to identify "phylogenetic footprints," regions with slower rates of sequence change than background. This does not require prior characterization of the sequence in question, therefore, it can be used in a wide range of applications. In particular, it is useful for the identification of functional elements in noncoding DNA, which are traditionally difficult to detect. Here, we describe in detail how to perform a simple yet powerful phylogenetic footprinting analysis. As an example, we use ribosomal DNA repeat sequences from various Saccharomyces yeasts to find functional noncoding DNA elements in the intergenic spacer, and explain critical considerations in performing phylogenetic footprinting analyses, including the number of species and species range, and some of the available software. Our methods are broadly applicable and should appeal to molecular biologists with little experience in bioinformatics.

We examined the effect of biochar on the water-soluble arsenic (As) concentration and the extent of organochlorine degradation in a co-contaminated historic sheep-dip soil during a 180-d glasshouse incubation experiment. Soil microbial activity, bacterial community and structure diversity were also investigated. Biochar made from willow feedstock (Salix sp) was pyrolysed at 350 or 550°C and added to soil at rates of 10 g kg-1 and 20 g kg-1 (representing 30 t ha-1 and 60 t ha-1). The isomers of hexachlorocyclohexane (HCH) alpha-HCH and gamma-HCH (lindane), underwent 10-fold and 4-fold reductions in concentration as a function of biochar treatment. Biochar also resulted in a significant reduction in soil DDT levels (P < 0.01), and increased the DDE:DDT ratio. Soil microbial activity was significantly increased (P < 0.01) under all biochar treatments after 60 days of treatment compared to the control. 16S amplicon sequencing revealed that biochar-amended soil contained more members of the Chryseobacterium, Flavobacterium, Dyadobacter and Pseudomonadaceae which are known bioremediators of hydrocarbons. We hypothesise that a recorded short-term reduction in the soluble As concentration due to biochar amendment allowed native soil microbial communities to overcome As-related stress. We propose that increased microbiological activity (dehydrogenase activity) due to biochar amendment was responsible for enhanced degradation of organochlorines in the soil. Biochar therefore partially overcame the co-contaminant effect of As, allowing for enhanced natural attenuation of organochlorines in soil.

Organisms maintain ribosomal RNA gene repeats (rDNA) at stable copy numbers by recombination; the loss of repeats results in gene amplification. Here we report a mechanism of amplification regulation. We show that amplification is dependent on transcription from a noncoding bidirectional promoter (E-pro) within the rDNA spacer. E-pro transcription stimulates the dissociation of cohesin, a DNA binding protein complex that suppresses sister-chromatid-based changes in rDNA copy number. This transcription is regulated by the silencing gene, SIR2 , and by copy number. Transcription-induced cohesin dissociation may be a general mechanism of recombination regulation.

Chytridiomycosis is a recently identified fungal disease associated with global population declines of frogs. Although the fungus, Batrachochytrium dendrobatidis , is considered an emerging pathogen, little is known about its population genetics, including the origin of the current epidemic and how this relates to the dispersal ability of the fungus. In this study, we use multilocus sequence typing to examine genetic diversity and relationships among 35 fungal strains from North America, Africa and Australia. Only five variable nucleotide positions were detected among 10 loci (5918 bp). This low level of genetic variation is consistent with the description of B. dendrobatidis as a recently emerged disease agent. Fixed (i.e. 100%) or nearly fixed frequencies of heterozygous genotypes at two loci suggested that B. dendrobatidis is diploid and primarily reproduces clonally. In contrast to the lack of nucleotide polymorphism, electrophoretic karyotyping of multiple strains demonstrated a number of chromosome length polymorphisms.

Polyploidy, a state in which the chromosome complement has undergone an increase, is a major force in evolution. Understanding the consequences of polyploidy has received much attention, and allopolyploids, which result from the union of two different parental genomes, are of particular interest because they must overcome a suite of biological responses to this merger, known as “genome shock.” A key question is what happens to gene expression of the two gene copies following allopolyploidization, but until recently the tools to answer this question on a genome-wide basis were lacking. Here we utilize high throughput transcriptome sequencing to produce the first genome-wide picture of gene expression response to allopolyploidy in fungi. A novel pipeline for assigning sequence reads to the gene copies was used to quantify their expression in a fungal allopolyploid. We find that the transcriptional response to allopolyploidy is predominantly conservative: both copies of most genes are retained; over half the genes inherit parental gene expression patterns; and parental differential expression is often lost in the allopolyploid. Strikingly, the patterns of gene expression change are highly concordant with the genome-wide expression results of a cotton allopolyploid. The very different nature of these two allopolyploids implies a conserved, eukaryote-wide transcriptional response to genome merger. We provide evidence that the transcriptional responses we observe are mostly driven by intrinsic differences between the regulatory systems in the parent species, and from this propose a mechanistic model in which the cross-kingdom conservation in transcriptional response reflects conservation of the mutational processes underlying eukaryotic gene regulatory evolution. This work provides a platform to develop a universal understanding of gene expression response to allopolyploidy and suggests that allopolyploids are an exceptional system to investigate gene regulatory changes that have evolved in the parental species prior to allopolyploidization., Author Summary Organisms are complex biological systems that must continue to function even as their genomes evolve. While evolution is usually gradual, the formation of new species by the hybridization of different parents—allopolyploidization—occurs nearly instantaneously. A key question is what happens to expression of the two parental gene copies following genome merger. To determine this, we focused on a fungal allopolyploid from a group that dominates many of the world's pastoral economies. To investigate the fate of gene expression in this system, we developed a novel pipeline to assign high throughput RNA sequence reads to the two parental gene copies, thus allowing quantification of expression. We found transcriptional responses to be predominantly conservative: most gene copies either inherit parental expression patterns, or if differentially expressed in the parents, that difference is lost in the hybrid. Moreover, we identified an extraordinary level of concordance in the fate of genome-wide allopolyploid gene expression with that seen in cotton. The very different nature of these two allopolyploids suggests that there is a set of universal rules for the transcriptional response to genome merger. We propose a mechanistic model whereby this conserved response reflects similarities in mutational processes that underlie gene regulatory evolution.

An extraordinary level of length heterogeneity was found in the ribosomal DNA (rDNA) of an asexual hybrid Neotyphodium grass endophyte, isolate Lp1. This hybrid Neotyphodium endophyte is an interspecific hybrid between two grass endophytes, Neotyphodium lolii, and a sexual form, Epichlöe typhina, and the length heterogeneity was not found in either of these progenitor species. The length heterogeneity in the hybrid is localized to the intergenic spacer (IGS) and is the result of copy-number variation of a tandemly repeated subrepeat class within the IGS, the 111-/119-bp subrepeats. Copy number variation of this subrepeat class appears to be a consequence of mitotic unequal crossing over that occurs between these subrepeats. This implies that unequal crossing over plays a role in the concerted evolution of the whole rDNA. Changes in the pattern of IGS length variants occurred in just two rounds of single-spore purification. Analysis of the IGS length heterogeneity revealed features that are unexpected in a simple model of unequal crossing over. Potential refinements of the molecular details of unequal crossing over are presented, and we also discuss evidence for a combination of homogenization mechanisms that drive the concerted evolution of the Lp1 rDNA.