Eucaryotic genome evolution through the spontaneous duplication of large chromosomal segments

High-quality genome assembly and comparative analysis reveal extensive genomic variation in Talaromyces marneffei

Talaromyces marneffei is a dimorphic fungus that transitions from a filamentous form at 25 °C to a pathogenic yeast form at 37 °C, demonstrating pathogenicity mostly in immunocompromised individuals, such as those with human immunodeficiency virus/AIDS. Though it is one of the most severe infectious fungi in Southeast Asia, the lack of comprehensive genomic analysis has hindered advancement in strain differentiation, diagnosis and treatment. In this study, we assembled a high-quality genome of T. marneffei ATCC 18224, resulting in a 28.9 Mb genome distributed across 11 contigs, using third-generation Oxford Nanopore Technologies sequencing reads. Notably, we identified a strain-specific 740-kb segmental duplication in strain ATCC 18224, potentially mediated by inserting a Ty1/Copia long terminal repeat (LTR) retrotransposon. This segmental duplication includes various functional genes, with 75 differentially expressed during its dimorphic transition. Comparative genomic analysis revealed large-scale rearrangements in strains PM1 and 11CN-20-091, which were inconsistent with the phylogenomic trees of six T. marneffei strains and required further investigation. Additionally, we observed substantial genetic structural variations in LTR retrotransposons, particularly within the Ty1/Copia family, including two significant recent expansions in strain ATCC 18224. In summary, the identification and characterization of these extensive genomic structural variations in T. marneffei contribute to a deep understanding of its genetic diversity and will facilitate improvements in genotyping, classification and genomic surveillance.

Sequence-dependent activity and compartmentalization of foreign DNA in a eukaryotic nucleus

In eukaryotes, DNA-associated protein complexes coevolve with genomic sequences to orchestrate chromatin folding. We investigate the relationship between DNA sequence and the spontaneous loading and activity of chromatin components in the absence of coevolution. Using bacterial genomes integrated into Saccharomyces cerevisiae, which diverged from yeast more than 2 billion years ago, we show that nucleosomes, cohesins, and associated transcriptional machinery can lead to the formation of two different chromatin archetypes, one transcribed and the other silent, independently of heterochromatin formation. These two archetypes also form on eukaryotic exogenous sequences, depend on sequence composition, and can be predicted using neural networks trained on the native genome. They do not mix in the nucleus, leading to a bipartite nuclear compartmentalization, reminiscent of the organization of vertebrate nuclei.

Mechanisms of tandem duplication in the cancer genome

Tandem duplications (TD) are among the most frequent type of structural variant (SV) in the cancer genome. They are characterized by a single breakpoint junction that defines the boundaries and the size of the duplicated segment. Cancer-associated TDs often increase oncogene copy number or disrupt tumor suppressor gene function, and thus have important roles in tumor evolution. TDs in cancer genomes fall into three classes, defined by the size of duplications, and are associated with distinct genetic drivers. In this review, we survey key features of cancer-related TDs and consider possible underlying mechanisms in relation to stressed DNA replication and the 3D organization of the S phase genome.

Mlg1, a yeast acyltransferase located in ER membranes associated with mitochondria (MAMs), is involved in de novo synthesis and remodelling of phospholipids

In cells, phospholipids contain acyl chains of variable lengths and saturation, features that affect their functions. Their de novo synthesis in the endoplasmic reticulum takes place via the cytidine diphosphate diacylglycerol (CDP-DAG) and Kennedy pathways, which are conserved in eukaryotes. PA is a key intermediate for all phospholipids (PI, PIPs, PS, PE, PC, PG and CL). The de novo synthesis of PA occurs by acylation of glycerophosphate leading to the synthesis of 1-acyl lysoPA and subsequent acylation of 1-acyl lysoPA at the sn-2 position. Using membranes from Escherichia coli overexpressing MLG1, we showed that the yeast gene MLG1 encodes an acyltransferase, leading specifically to the synthesis of PA from 1-acyl lysoPA. Moreover, after their de novo synthesis, phospholipids can be remodelled by acyl exchange with one and/or two acyl chains exchanged at the sn-1 and/or sn-2 position. Based on shotgun lipidomics of the reference and mlg1Δ strains, as well as biochemical assays for acyltransferase activities, we identified an additional remodelling activity for Mlg1p, namely, incorporation of palmitic acid into the sn-1 position of PS and PE. By using confocal microscopy and subcellular fractionation, we also found that this acyltransferase is located in ER membranes associated with mitochondria, a finding that highlights the importance of these organelles in the global cellular metabolism of lipids.

The road less travelled? Exploring the nuanced evolutionary consequences of duplicated genes

Duplicated genes have long been appreciated as both substrates and catalysts of evolutionary processes. From even the simplest cell to complex multicellular animals and plants, duplicated genes have made immeasurable contributions to the phenotypic evolution of all life on Earth. Not merely drivers of morphological innovation and speciation events, however, gene duplications sculpt the evolution of genetic architecture in ways we are only just coming to understand now we have the experimental tools to do so. As such, the present article revisits our understanding of the ways in which duplicated genes evolve, examining closely the various fates they can adopt in light of recent work that yields insights from studies of paralogues from across the tree of life that challenge the classical framework.

Megasatellite formation and evolution in vertebrate genes

Since formation of the first proto-eukaryotes, gene repertoire and genome complexity have significantly increased. Among genetic elements responsible for this increase are tandem repeats. Here we describe a genome-wide analysis of large tandem repeats, called megasatellites, in 58 vertebrate genomes. Two bursts occurred, one after the radiation between Agnatha and Gnathostomata fishes and the second one in therian mammals. Megasatellites are enriched in subtelomeric regions and frequently encoded in genes involved in transcription regulation, intracellular trafficking, and cell membrane metabolism, reminiscent of what is observed in fungus genomes. The presence of many introns within young megasatellites suggests that an exon-intron DNA segment is first duplicated and amplified before accumulation of mutations in intronic parts partially erases the megasatellite in such a way that it becomes detectable only in exons. Our results suggest that megasatellite formation and evolution is a dynamic and still ongoing process in vertebrate genomes.

Imaging Intron Evolution

Intron evolution may be readily imaged through the combined use of the “dot plot” function of the NCBI BLAST, aligning two sequences at a time, and the Vertebrate “Multiz” alignment and conservation tool of the UCSC Genome Browser. With the NCBI BLAST, an ideal alignment of two highly conserved sequences generates a diagonal straight line in the plot from the lower left corner to the upper right corner. Gaps in this line correspond to non-conserved sections. In addition, the dot plot of the alignment of a sequence with the same sequence after the removal of the Transposable Elements (TEs) can be observed along the diagonal gaps that correspond to the sites of TE insertion. The UCSC Genome Browser can graph, along the entire sequence of a single gene, the level of overall conservation in vertebrates. This level can be compared with the conservation level of the gene in one or more selected vertebrate species. As an example, we show the graphic analysis of the intron conservation in two genes: the mitochondrial solute carrier 21 (SLC25A21) and the growth hormone receptor (GHR), whose coding sequences are conserved through vertebrates, while their introns show dramatic changes in nucleotide composition and even length. In the SLC25A21, a few short but significant nucleotide sequences are conserved in zebrafish, Xenopus and humans, and the rate of conservation steadily increases from chicken/human to mouse/human alignments. In the GHR, a less conserved gene, the earlier indication of intron conservation is a small signal in chicken/human alignment. The UCSC tool may simultaneously display the conservation level of a gene in different vertebrates, with reference to the level of overall conservation in Vertebrates. It is shown that, at least in SLC25A21, the sites of higher conservation are not always coincident in chicken and zebrafish nor are the sites of higher vertebrate conservation.

Comprehensive analysis of cis- and trans-acting factors affecting ectopic Break-Induced Replication

Break-induced replication (BIR) is a highly mutagenic eukaryotic homologous DNA recombination pathway that repairs one-ended DNA double strand breaks such as broken DNA replication forks and eroded telomeres. While searching for cis-acting factors regulating ectopic BIR efficiency, we found that ectopic BIR efficiency is the highest close to chromosome ends. The variations of ectopic BIR efficiency as a function of the length of DNA to replicate can be described as a combination of two decreasing exponential functions, a property in line with repeated cycles of strand invasion, elongation and dissociation that characterize BIR. Interestingly, the apparent processivity of ectopic BIR depends on the length of DNA already synthesized. Ectopic BIR is more susceptible to disruption during the synthesis of the first ~35–40 kb of DNA than later, notably when the template chromatid is being transcribed or heterochromatic. Finally, we show that the Srs2 helicase promotes ectopic BIR from both telomere proximal and telomere distal regions in diploid cells but only from telomere proximal sites in haploid cells. Altogether, we bring new light on the factors impacting a last resort DNA repair pathway.

DNA is a long molecule composed of two anti-parallel strands that can undergo breaks that need to be efficiently repaired to ensure genomic stability, hence preventing genetic diseases such as cancer. Homologous recombination is a major DNA repair pathway that copies DNA from intact homologous templates to seal DNA double strand breaks. Short DNA repair tracts are favored when homologous sequences for the two extremities of the broken molecule are present. However, when homologous sequences are present for only one extremity of the broken molecule, DNA repair synthesis can proceed up to the end of the chromosome, the telomere. This notably occurs at eroded telomeres when telomerase, the enzyme normally responsible for telomere elongation, is inactive, and at broken DNA replication intermediates. However, this Break-Induced Replication or BIR pathway is highly mutagenic. By initiating BIR at various distances from the telomere, we found that the length of DNA to synthesize significantly reduces BIR efficiency. Interestingly, our findings support two DNA synthesis phases, the first one being much less processive than the second one. Ultimately, this tends to restrain the use of this last resort DNA repair pathway to chromosome extremities notably when it takes place between non-allelic homologous sequences.

Can Introns Stabilize Gene Duplication?

Gene duplication is considered one of the most important events that determine the evolution of genomes. However, the neo-duplication condition of a given gene is particularly unstable due to recombination events. Several mechanisms have been proposed to justify this step. In this "opinion article" we propose a role for intron sequences in stabilizing gene duplication by limiting and reducing the identity of the gene sequence between the two duplicated copies. A review of the topic and a detailed hypothesis are presented.

Genome-wide effect of non-optimal temperatures under anaerobic conditions on gene expression in Saccharomyces cerevisiae

Understanding of thermal adaptation mechanisms in yeast is crucial to develop better-adapted strains to industrial processes, providing more economical and sustainable products. We have analyzed the transcriptomic responses of three Saccharomyces cerevisiae strains, a commercial wine strain, ADY5, a laboratory strain, CEN.PK113-7D and a commercial bioethanol strain, Ethanol Red, grown at non-optimal temperatures under anaerobic chemostat conditions. Transcriptomic analysis of the three strains revealed a huge complexity of cellular mechanisms and responses. Overall, cold exerted a stronger transcriptional response in the three strains comparing with heat conditions, with a higher number of down-regulating genes than of up-regulating genes regardless the strain analyzed. The comparison of the transcriptome at both sub- and supra-optimal temperatures showed the presence of common genes up- or down-regulated in both conditions, but also the presence of common genes up- or down-regulated in the three studied strains. More specifically, we have identified and validated three up-regulated genes at sub-optimal temperature in the three strains, OPI3, EFM6 and YOL014W. Finally, the comparison of the transcriptomic data with a previous proteomic study with the same strains revealed a good correlation between gene activity and protein abundance, mainly at low temperature. Our work provides a global insight into the specific mechanisms involved in temperature adaptation regarding both transcriptome and proteome, which can be a step forward in the comprehension and improvement of yeast thermotolerance.

Functional variability in adhesion and flocculation of yeast megasatellite genes

Megasatellites are large tandem repeats found in all fungal genomes but especially abundant in the opportunistic pathogen Candida glabrata. They are encoded in genes involved in cell-cell interactions, either between yeasts or between yeast and human cells. In the present work, we have been using an iterative genetic system to delete several Candida glabrata megasatellite-containing genes and found that 2 of them were positively involved in adhesion to epithelial cells, whereas 3 genes negatively controlled adhesion. Two of the latter, CAGL0B05061g or CAGL0A04851g, were also negative regulators of yeast-to-yeast adhesion, making them central players in controlling Candida glabrata adherence properties. Using a series of synthetic Saccharomyces cerevisiae strains in which the FLO1 megasatellite was replaced by other tandem repeats of similar length but different sequences, we showed that the capacity of a strain to flocculate in liquid culture was unrelated to its capacity to adhere to epithelial cells or to invade agar. Finally, to understand how megasatellites were initially created and subsequently expanded, an experimental evolution system was set up, in which modified yeast strains containing different megasatellite seeds were grown in bioreactors for more than 200 generations and selected for their ability to sediment at the bottom of the culture tube. Several flocculation-positive mutants were isolated. Functionally relevant mutations included general transcription factors as well as a 230-kbp segmental duplication.

Retention of duplicated genes in evolution

Gene duplication is a prevalent phenomenon across the tree of life. The processes that lead to the retention of duplicated genes are not well understood. Functional genomics approaches in model organisms, such as yeast, provide useful tools to test the mechanisms underlying retention with functional redundancy and divergence of duplicated genes, including fates associated with neofunctionalization, subfunctionalization, back-up compensation, and dosage amplification. Duplicated genes may also be retained as a consequence of structural and functional entanglement. Advances in human gene editing have enabled the interrogation of duplicated genes in the human genome, providing new tools to evaluate the relative contributions of each of these factors to duplicate gene retention and the evolution of genome structure.

Cohesin regulates homology search during recombinational DNA repair

Homologous recombination repairs DNA double-strand breaks (DSB) using an intact dsDNA molecule as a template. It entails a homology search step, carried out along a conserved RecA/Rad51-ssDNA filament assembled on each DSB end. Whether, how and to what extent a DSB impacts chromatin folding, and how this (re)organization in turns influences the homology search process, remain ill-defined. Here we characterize two layers of spatial chromatin reorganization following DSB formation in Saccharomyces cerevisiae. Although cohesin folds chromosomes into cohesive arrays of ~20-kb-long chromatin loops as cells arrest in G2/M, the DSB-flanking regions interact locally in a resection- and 9-1-1 clamp-dependent manner, independently of cohesin, Mec1ATR, Rad52 and Rad51. This local structure blocks cohesin progression, constraining the DSB region at the base of a loop. Functionally, cohesin promotes DSB-dsDNA interactions and donor identification in cis, while inhibiting them in trans. This study identifies multiple direct and indirect ways by which cohesin regulates homology search during recombinational DNA repair.

Quantum Patterns of Genome Size Variation in Angiosperms

Aims:

The discontinuous pattern of genome size variation in angiosperms is an unsolved problem related to genome evolution. In this study, we introduced a genome evolution operator and solved the related eigenvalue equation to deduce the discontinuous pattern.

Background:

Genome is a well-defined system for studying the evolution of species. One of the basic problems is the genome size evolution. The DNA amounts for angiosperm species are highly variable, differing over 1000-fold. One big surprise is the discovery of the discontinuous distribution of nuclear DNA amounts in many angiosperm genera.

Objective:

The discontinuous distribution of nuclear DNA amounts has certain regularity, much like a group of quantum states in atomic physics. The quantum pattern has not been explained by all the evolutionary theories so far and we shall interpret it through the quantum simulation of genome evolution.

Methods:

We introduced a genome evolution operator H to deduce the distribution of DNA amount. The nuclear DNA amount in angiosperms is studied from the eigenvalue equation of the genome evolution operator H. The operator H is introduced by physical simulation and it is defined as a function of the genome size N and the derivative with respect to the size.

Results:

The discontinuity of DNA size distribution and its synergetic occurrence in related angiosperms species are successfully deduced from the solution of the equation. The results agree well with the existing experimental data of Aloe, Clarkia, Nicotiana, Lathyrus, Allium and other genera.

Conclusion:

The success of our approach may infer the existence of a set of genomic evolutionary equations satisfying classical-quantum duality. The classical phase of evolution means it obeys the classical deterministic law, while the quantum phase means it obeys the quantum stochastic law. The discontinuity of DNA size distribution provides novel evidences on the quantum evolution of angiosperms. It has been realized that the discontinuous pattern is due to the existence of some unknown evolutionary constraints. However, our study indicates that these constraints on the angiosperm genome essentially originate from quantum.

Thermo-adaptive evolution to generate improved Saccharomyces cerevisiae strains for cocoa pulp fermentations

Cocoa pulp fermentation is a consequence of the succession of indigenous yeasts, lactic acid bacteria and acetic acid bacteria that not only produce a diversity of metabolites, but also cause the production of flavour precursors. However, as such spontaneous fermentations are less reproducible and contribute to produce variability, interest in a microbial starter culture is growing that could be used to inoculate cocoa pulp fermentations. This study aimed to generate robust S. cerevisiae strains by thermo-adaptive evolution that could be used in cocoa fermentation. We evolved a cocoa strain in a sugary defined medium at high temperature to improve both fermentation and growth capacity. Moreover, adaptive evolution at high temperature (40 °C) also enabled us to unveil the molecular basis underlying the improved phenotype by analysing the whole genome sequence of the evolved strain. Adaptation to high-temperature conditions occurred at different genomic levels, and promoted aneuploidies, segmental duplication, and SNVs in the evolved strain. The lipid profile analysis of the evolved strain also evidenced changes in the membrane composition that contribute to maintain an appropriate cell membrane state at high temperature. Our work demonstrates that experimental evolution is an effective approach to generate better-adapted yeast strains at high temperature for industrial processes.

Transposon-mediated telomere destabilization: a driver of genome evolution in the blast fungus

The fungus Magnaporthe oryzae causes devastating diseases of crops, including rice and wheat, and in various grasses. Strains from ryegrasses have highly unstable chromosome ends that undergo frequent rearrangements, and this has been associated with the presence of retrotransposons (Magnaporthe oryzae Telomeric Retrotransposons-MoTeRs) inserted in the telomeres. The objective of the present study was to determine the mechanisms by which MoTeRs promote telomere instability. Targeted cloning, mapping, and sequencing of parental and novel telomeric restriction fragments (TRFs), along with MinION sequencing of genomic DNA allowed us to document the precise molecular alterations underlying 109 newly-formed TRFs. These included truncations of subterminal rDNA sequences; acquisition of MoTeR insertions by 'plain' telomeres; insertion of the MAGGY retrotransposons into MoTeR arrays; MoTeR-independent expansion and contraction of subtelomeric tandem repeats; and a variety of rearrangements initiated through breaks in interstitial telomere tracts that are generated during MoTeR integration. Overall, we estimate that alterations occurred in approximately sixty percent of chromosomes (one in three telomeres) analyzed. Most importantly, we describe an entirely new mechanism by which transposons can promote genomic alterations at exceptionally high frequencies, and in a manner that can promote genome evolution while minimizing collateral damage to overall chromosome architecture and function.

An Overview of Duplicated Gene Detection Methods: Why the Duplication Mechanism Has to Be Accounted for in Their Choice

Gene duplication is an important evolutionary mechanism allowing to provide new genetic material and thus opportunities to acquire new gene functions for an organism, with major implications such as speciation events. Various processes are known to allow a gene to be duplicated and different models explain how duplicated genes can be maintained in genomes. Due to their particular importance, the identification of duplicated genes is essential when studying genome evolution but it can still be a challenge due to the various fates duplicated genes can encounter. In this review, we first describe the evolutionary processes allowing the formation of duplicated genes but also describe the various bioinformatic approaches that can be used to identify them in genome sequences. Indeed, these bioinformatic approaches differ according to the underlying duplication mechanism. Hence, understanding the specificity of the duplicated genes of interest is a great asset for tool selection and should be taken into account when exploring a biological question.

Genomic Multiplication and Drug Efflux Influence Ketoconazole Resistance in Malassezia restricta

Malassezia restricta is an opportunistic fungal pathogen on human skin; it is associated with various skin diseases, including seborrheic dermatitis and dandruff, which are usually treated using ketoconazole. In this study, we clinically isolated ketoconazole-resistant M. restricta strains (KCTC 27529 and KCTC 27550) from patients with dandruff. To understand the mechanisms of ketoconazole resistance in the isolates, their genomes were sequenced and compared with the susceptible reference strain M. restricta KCTC 27527. Using comparative genome analysis, we identified tandem multiplications of the genomic loci containing ATM1 and ERG11 homologs in M. restricta KCTC 27529 and KCTC 27550, respectively. Additionally, we found that the copy number increase of ATM1 and ERG11 is reflected in the increased expression of these genes; moreover, we observed that overexpression of these homologs caused ketoconazole resistance in a genetically tractable fungal pathogen, Cryptococcus neoformans. In addition to tandem multiplications of the genomic region containing the ATM1 homolog, the PDR5 homolog, which encodes the drug efflux pump protein was upregulated in M. restricta KCTC 27529 compared to the reference strain. Biochemical analysis confirmed that drug efflux was highly activated in M. restricta KCTC 27529, implying that upregulation of the PDR5 homolog may also contribute to ketoconazole resistance in the strain. Overall, our results suggest that multiplication of the genomic loci encoding genes involved in ergosterol synthesis, mitochondrial iron metabolism, and oxidative stress response and overexpression of the drug efflux pumps are the mechanisms underlying ketoconazole resistance in M. restricta.

Aneuploidy Enables Cross-Adaptation to Unrelated Drugs

Aneuploidy is common both in tumor cells responding to chemotherapeutic agents and in fungal cells adapting to antifungal drugs. Because aneuploidy simultaneously affects many genes, it has the potential to confer multiple phenotypes to the same cells. Here, we analyzed the mechanisms by which Candida albicans, the most prevalent human fungal pathogen, acquires the ability to survive both chemotherapeutic agents and antifungal drugs. Strikingly, adaptation to both types of drugs was accompanied by the acquisition of specific whole-chromosome aneuploidies, with some aneuploid karyotypes recovered independently and repeatedly from very different drug conditions. Specifically, strains selected for survival in hydroxyurea, an anticancer drug, acquired cross-adaptation to caspofungin, a first-line antifungal drug, and both acquired traits were attributable to trisomy of the same chromosome: loss of trisomy was accompanied by loss of adaptation to both drugs. Mechanistically, aneuploidy simultaneously altered the copy number of most genes on chromosome 2, yet survival in hydroxyurea or caspofungin required different genes and stress response pathways. Similarly, chromosome 5 monosomy conferred increased tolerance to both fluconazole and to caspofungin, antifungals with different mechanisms of action. Thus, the potential for cross-adaptation is not a feature of aneuploidy per se; rather, it is dependent on specific genes harbored on given aneuploid chromosomes. Furthermore, pre-exposure to hydroxyurea increased the frequency of appearance of caspofungin survivors, and hydroxyurea-adapted C. albicans cells were refractory to antifungal drug treatment in a mouse model of systemic candidiasis. This highlights the potential clinical consequences for the management of cancer chemotherapy patients at risk of fungal infections.

Experimental Evolution of Yeast for High-Temperature Tolerance

Thermotolerance is a polygenic trait that contributes to cell survival and growth under unusually high temperatures. Although some genes associated with high-temperature growth (Htg+) have been identified, how cells accumulate mutations to achieve prolonged thermotolerance is still mysterious. Here, we conducted experimental evolution of a Saccharomyces cerevisiae laboratory strain with stepwise temperature increases for it to grow at 42 °C. Whole genome resequencing of 14 evolved strains and the parental strain revealed a total of 153 mutations in the evolved strains, including single nucleotide variants, small INDELs, and segmental duplication/deletion events. Some mutations persisted from an intermediate temperature to 42 °C, so they might be Htg+ mutations. Functional categorization of mutations revealed enrichment of exonic mutations in the SWI/SNF complex and F-type ATPase, pointing to their involvement in high-temperature tolerance. In addition, multiple mutations were found in a general stress-associated signal transduction network consisting of Hog1 mediated pathway, RAS-cAMP pathway, and Rho1-Pkc1 mediated cell wall integrity pathway, implying that cells can achieve Htg+ partly through modifying existing stress regulatory mechanisms. Using pooled segregant analysis of five Htg+ phenotype-orientated pools, we inferred causative mutations for growth at 42 °C and identified those mutations with stronger impacts on the phenotype. Finally, we experimentally validated a number of the candidate Htg+ mutations. This study increased our understanding of the genetic basis of yeast tolerance to high temperature.

An evolving view of copy number variants

Copy number variants (CNVs) are regions of the genome that vary in integer copy number. CNVs, which comprise both amplifications and deletions of DNA sequence, have been identified across all domains of life, from bacteria and archaea to plants and animals. CNVs are an important source of genetic diversity, and can drive rapid adaptive evolution and progression of heritable and somatic human diseases, such as cancer. However, despite their evolutionary importance and clinical relevance, CNVs remain understudied compared to single-nucleotide variants (SNVs). This is a consequence of the inherent difficulties in detecting CNVs at low-to-intermediate frequencies in heterogeneous populations of cells. Here, we discuss molecular methods used to detect CNVs, the limitations associated with using these techniques, and the application of new and emerging technologies that present solutions to these challenges. The goal of this short review and perspective is to highlight aspects of CNV biology that are understudied and define avenues for further research that address specific gaps in our knowledge of these complex alleles. We describe our recently developed method for CNV detection in which a fluorescent gene functions as a single-cell CNV reporter and present key findings from our evolution experiments in Saccharomyces cerevisiae. Using a CNV reporter, we found that CNVs are generated at a high rate and undergo selection with predictable dynamics across independently evolving replicate populations. Many CNVs appear to be generated through DNA replication-based processes that are mediated by the presence of short, interrupted, inverted-repeat sequences. Our results have important implications for the role of CNVs in evolutionary processes and the molecular mechanisms that underlie CNV formation. We discuss the possible extension of our method to other applications, including tracking the dynamics of CNVs in models of human tumors.

Homologous Recombination and the Formation of Complex Genomic Rearrangements

The maintenance of genome integrity involves multiple independent DNA damage avoidance and repair mechanisms. However, the origin and pathways of the focal chromosomal reshuffling phenomena collectively referred to as chromothripsis remain mechanistically obscure. We discuss here the role, mechanisms, and regulation of homologous recombination (HR) in the formation of simple and complex chromosomal rearrangements. We emphasize features of the recently characterized multi-invasion (MI)-induced rearrangement (MIR) pathway which uniquely amplifies the initial DNA damage. HR intermediates and cellular contexts that endanger genomic stability are discussed as well as the emerging roles of various classes of nucleases in the formation of genome rearrangements. Long-read sequencing and improved mapping of repeats should enable better appreciation of the significance of recombination in generating genomic rearrangements.

Adaptive deletion in resistance gene duplications in the malaria vector Anopheles gambiae

While gene copy-number variations play major roles in long-term evolution, their early dynamics remains largely unknown. However, examples of their role in short-term adaptation are accumulating: identical repetitions of a locus (homogeneous duplications) can provide a quantitative advantage, while the association of differing alleles (heterogeneous duplications) allows carrying two functions simultaneously. Such duplications often result from rearrangements of sometimes relatively large chromosome fragments, and even when adaptive, they can be associated with deleterious side effects that should, however, be reduced by subsequent evolution. Here, we took advantage of the unique model provided by the malaria mosquito Anopheles gambiae s.l. to investigate the early evolution of several duplications, heterogeneous and homogeneous, segregating in natural populations from West Africa. These duplications encompass ~200 kb and 11 genes, including the adaptive insecticide resistance ace-1 locus. Through the survey of several populations from three countries over 3-4 years, we showed that an internal deletion of all coamplified genes except ace-1 is currently spreading in West Africa and introgressing from An. gambiae s.s. to An. coluzzii. Both observations provide evidences of its selection, most likely due to reducing the gene-dosage disturbances caused by the excessive copies of the nonadaptive genes. Our study thus provides a unique example of the early adaptive trajectory of duplications and underlines the role of the environmental conditions (insecticide treatment practices and species ecology). It also emphasizes the striking diversity of adaptive responses in these mosquitoes and reveals a worrisome process of resistance/cost trade-off evolution that could impact the control of malaria vectors in Africa.

Border collies of the genome: domestication of an autonomous retrovirus-like transposon

Retrotransposons often spread rapidly through eukaryotic genomes until they are neutralized by host-mediated silencing mechanisms, reduced by recombination and mutation, and lost or transformed into benevolent entities. But the Ty1 retrotransposon appears to have been domesticated to guard the genome of Saccharomyces cerevisiae.

Population size changes and selection drive patterns of parallel evolution in a host-virus system

Predicting the repeatability of evolution remains elusive. Theory and empirical studies suggest that strong selection and large population sizes increase the probability for parallel evolution at the phenotypic and genotypic levels. However, selection and population sizes are not constant, but rather change continuously and directly affect each other even on short time scales. Here, we examine the degree of parallel evolution shaped through eco-evolutionary dynamics in an algal host population coevolving with a virus. We find high degrees of parallelism at the level of population size changes (ecology) and at the phenotypic level between replicated populations. At the genomic level, we find evidence for parallelism, as the same large genomic region was duplicated in all replicated populations, but also substantial novel sequence divergence between replicates. These patterns of genome evolution can be explained by considering population size changes as an important driver of rapid evolution.

Stoichiometry and compositional plasticity of the yeast nuclear pore complex revealed by quantitative fluorescence microscopy

The nuclear pore complex (NPC) is an eightfold symmetrical channel providing selective transport of biomolecules across the nuclear envelope. Each NPC consists of ∼30 different nuclear pore proteins (Nups) all present in multiple copies per NPC. Significant progress has recently been made in the characterization of the vertebrate NPC structure. However, because of the estimated size differences between the vertebrate and yeast NPC, it has been unclear whether the NPC architecture is conserved between species. Here, we have developed a quantitative image analysis pipeline, termed nuclear rim intensity measurement (NuRIM), to precisely determine copy numbers for almost all Nups within native NPCs of budding yeast cells. Our analysis demonstrates that the majority of yeast Nups are present at most in 16 copies per NPC. This reveals a dramatic difference to the stoichiometry determined for the human NPC, suggesting that despite a high degree of individual Nup conservation, the yeast and human NPC architecture is significantly different. Furthermore, using NuRIM, we examined the effects of mutations on NPC stoichiometry. We demonstrate for two paralog pairs of key scaffold Nups, Nup170/Nup157 and Nup192/Nup188, that their altered expression leads to significant changes in the NPC stoichiometry inducing either voids in the NPC structure or substitution of one paralog by the other. Thus, our results not only provide accurate stoichiometry information for the intact yeast NPC but also reveal an intriguing compositional plasticity of the NPC architecture, which may explain how differences in NPC composition could arise in the course of evolution.

Analyzing Genome Rearrangements in Saccharomyces cerevisiae

Genome rearrangements underlie different human diseases including many cancers. Determining the rates at which genome rearrangements arise and isolating unique, independent genome rearrangements is critical to understanding the genes and pathways that prevent or promote genome rearrangements. Here, we describe quantitative S. cerevisiae genetic assays for measuring the rates of accumulating genome rearrangements including deletions, translocations, and broken chromosomes healed by de novo telomere addition that result in the deletion of two counter-selectable genes, CAN1 and URA3, placed in the nonessential regions of the S. cerevisiae genome. The assays also allow for the isolation of individual genome rearrangements for structural studies, and a method for analyzing genome rearrangements by next-generation DNA sequencing is provided.

Pathways and Mechanisms that Prevent Genome Instability in Saccharomyces cerevisiae

Genome rearrangements result in mutations that underlie many human diseases, and ongoing genome instability likely contributes to the development of many cancers. The tools for studying genome instability in mammalian cells are limited, whereas model organisms such as Saccharomyces cerevisiae are more amenable to these studies. Here, we discuss the many genetic assays developed to measure the rate of occurrence of Gross Chromosomal Rearrangements (called GCRs) in S. cerevisiae. These genetic assays have been used to identify many types of GCRs, including translocations, interstitial deletions, and broken chromosomes healed by de novo telomere addition, and have identified genes that act in the suppression and formation of GCRs. Insights from these studies have contributed to the understanding of pathways and mechanisms that suppress genome instability and how these pathways cooperate with each other. Integrated models for the formation and suppression of GCRs are discussed.

3D organization of synthetic and scrambled chromosomes

Although the design of the synthetic yeast genome Sc2.0 is highly conservative with respect to gene content, the deletion of several classes of repeated sequences and the introduction of thousands of designer changes may affect genome organization and potentially alter cellular functions. We report here the Hi-C-determined three-dimensional (3D) conformations of Sc2.0 chromosomes. The absence of repeats leads to a smoother contact pattern and more precisely tractable chromosome conformations, and the large-scale genomic organization is globally unaffected by the presence of synthetic chromosome(s). Two exceptions are synIII, which lacks the silent mating-type cassettes, and synXII, specifically when the ribosomal DNA is moved to another chromosome. We also exploit the contact maps to detect rearrangements induced in SCRaMbLE (synthetic chromosome rearrangement and modification by loxP-mediated evolution) strains.

Sources of Fungal Genetic Variation and Associating It with Phenotypic Diversity

The first eukaryotic genome to be sequenced was fungal, and there continue to be more sequenced genomes in the kingdom Fungi than in any other eukaryotic kingdom. Comparison of these genomes reveals many sources of genetic variation, from single nucleotide polymorphisms to horizontal gene transfer and on to changes in the arrangement and number of chromosomes, not to mention endofungal bacteria and viruses. Population genomics shows that all sources generate variation all the time and implicate natural selection as the force maintaining genome stability. Variation in wild populations is a rich resource for associating genetic variation with phenotypic variation, whether through quantitative trait locus mapping, genome-wide association studies, or reverse ecology. Subjects of studies associating genetic and phenotypic variation include model fungi, e.g., Saccharomyces and Neurospora, but pioneering studies have also been made with fungi pathogenic to plants, e.g., Pyricularia (= Magnaporthe), Zymoseptoria, and Fusarium, and to humans, e.g., Coccidioides, Cryptococcus, and Candida.

Genome Diversity and Evolution in the Budding Yeasts (Saccharomycotina)

Considerable progress in our understanding of yeast genomes and their evolution has been made over the last decade with the sequencing, analysis, and comparisons of numerous species, strains, or isolates of diverse origins. The role played by yeasts in natural environments as well as in artificial manufactures, combined with the importance of some species as model experimental systems sustained this effort. At the same time, their enormous evolutionary diversity (there are yeast species in every subphylum of Dikarya) sparked curiosity but necessitated further efforts to obtain appropriate reference genomes. Today, yeast genomes have been very informative about basic mechanisms of evolution, speciation, hybridization, domestication, as well as about the molecular machineries underlying them. They are also irreplaceable to investigate in detail the complex relationship between genotypes and phenotypes with both theoretical and practical implications. This review examines these questions at two distinct levels offered by the broad evolutionary range of yeasts: inside the best-studied Saccharomyces species complex, and across the entire and diversified subphylum of Saccharomycotina. While obviously revealing evolutionary histories at different scales, data converge to a remarkably coherent picture in which one can estimate the relative importance of intrinsic genome dynamics, including gene birth and loss, vs. horizontal genetic accidents in the making of populations. The facility with which novel yeast genomes can now be studied, combined with the already numerous available reference genomes, offer privileged perspectives to further examine these fundamental biological questions using yeasts both as eukaryotic models and as fungi of practical importance.

The rise of genomics

A brief history of the development of genomics is provided. Complete sequencing of genomes of uni- and multicellular organisms is based on important progress in sequencing and bioinformatics. Evolution of these methods is ongoing and has triggered an explosion in data production and analysis. Initial analyses focused on the inventory of genes encoding proteins. Completeness and quality of gene prediction remains crucial. Genome analyses profoundly modified our views on evolution, biodiversity and contributed to the detection of new functions, yet to be fully elucidated, such as those fulfilled by non-coding RNAs. Genomics has become the basis for the study of biology and provides the molecular support for a bunch of large-scale studies, the omics.

Massive Amplification at an Unselected Locus Accompanies Complex Chromosomal Rearrangements in Yeast

Gene amplification has been observed in different organisms in response to environmental constraints, such as limited nutrients or exposure to a variety of toxic compounds, conferring them with specific phenotypic adaptations via increased expression levels. However, the presence of multiple gene copies in natural genomes has generally not been found in the absence of specific functional selection. Here, we show that the massive amplification of a chromosomal locus (up to 880 copies per cell) occurs in the absence of any direct selection, and is associated with low-order amplifications of flanking segments in complex chromosomal alterations. These results were obtained from mutants with restored phenotypes that spontaneously appeared from genetically engineered strains of the yeast Saccharomyces cerevisiae suffering from severe fitness reduction. Grossly extended chromosomes (macrotene) were formed, with complex structural alterations but sufficient stability to propagate unchanged over successive generations. Their detailed molecular analysis, including complete genome sequencing, identification of sequence breakpoints, and comparisons between mutants, revealed novel mechanisms causing their formation, whose combined action underlies the astonishing dynamics of eukaryotic chromosomes and their consequences.

The tandem duplicator phenotype as a distinct genomic configuration in cancer

Next-generation sequencing studies have revealed genome-wide structural variation patterns in cancer, such as chromothripsis and chromoplexy, that do not engage a single discernable driver mutation, and whose clinical relevance is unclear. We devised a robust genomic metric able to identify cancers with a chromotype called tandem duplicator phenotype (TDP) characterized by frequent and distributed tandem duplications (TDs). Enriched only in triple-negative breast cancer (TNBC) and in ovarian, endometrial, and liver cancers, TDP tumors conjointly exhibit tumor protein p53 (TP53) mutations, disruption of breast cancer 1 (BRCA1), and increased expression of DNA replication genes pointing at rereplication in a defective checkpoint environment as a plausible causal mechanism. The resultant TDs in TDP augment global oncogene expression and disrupt tumor suppressor genes. Importantly, the TDP strongly correlates with cisplatin sensitivity in both TNBC cell lines and primary patient-derived xenografts. We conclude that the TDP is a common cancer chromotype that coordinately alters oncogene/tumor suppressor expression with potential as a marker for chemotherapeutic response.

Diversification and coevolution of the ghrelin/growth hormone secretagogue receptor system in vertebrates

The gut hormone ghrelin is involved in numerous metabolic functions, such as the stimulation of growth hormone secretion, gastric motility, and food intake. Ghrelin is modified by ghrelin O-acyltransferase (GOAT) or membrane-bound O-acyltransferase domain-containing 4 (MBOAT4) enabling action through the growth hormone secretagogue receptors (GHS-R). During the course of evolution, initially strong ligand/receptor specificities can be disrupted by genomic changes, potentially modifying physiological roles of the ligand/receptor system. Here, we investigated the coevolution of ghrelin, GOAT, and GHS-R in vertebrates. We combined similarity search, conserved synteny analyses, phylogenetic reconstructions, and protein structure comparisons to reconstruct the evolutionary history of the ghrelin system. Ghrelin remained a single-gene locus in all vertebrate species, and accordingly, a single GHS-R isoform was identified in all tetrapods. Similar patterns of the nonsynonymous (dN) and synonymous (dS) ratio (dN/dS) in the vertebrate lineage strongly suggest coevolution of the ghrelin and GHS-R genes, supporting specific functional interactions and common physiological pathways. The selection profiles do not allow confirmation as to whether ghrelin binds specifically to GOAT, but the ghrelin dN/dS patterns are more similar to those of GOAT compared to MBOAT1 and MBOAT2 isoforms. Four GHS-R isoforms were identified in teleost genomes. This diversification of GHS-R resulted from successive rounds of duplications, some of which remained specific to the teleost lineage. Coevolution signals are lost in teleosts, presumably due to the diversification of GHS-R but not the ghrelin gene. The identification of the GHS-R diversity in teleosts provides a molecular basis for comparative studies on ghrelin's physiological roles and regulation, while the comparative sequence and structure analyses will assist translational medicine to determine structure-function relationships of the ghrelin/GHS-R system.

Molecular Combing of Single DNA Molecules on the 10 Megabase Scale

DNA combing allows the investigation of DNA replication on genomic single DNA molecules, but the lengths that can be analysed have been restricted to molecules of 200-500 kb. We have improved the DNA combing procedure so that DNA molecules can be analysed up to the length of entire chromosomes in fission yeast and up to 12 Mb fragments in human cells. Combing multi-Mb-scale DNA molecules revealed previously undetected origin clusters in fission yeast and shows that in human cells replication origins fire stochastically forming clusters of fired origins with an average size of 370 kb. We estimate that a single human cell forms around 3200 clusters at mid S-phase and fires approximately 100,000 origins to complete genome duplication. The procedure presented here will be adaptable to other organisms and experimental conditions.

Deciphering the Code of the Cancer Genome: Mechanisms of Chromosome Rearrangement

Chromosome rearrangement plays a causal role in tumorigenesis by contributing to the inactivation of tumor suppressor genes, the dysregulated expression or amplification of oncogenes and the generation of novel gene fusions. Chromosome breaks are important intermediates in this process. How, when and where these breaks arise and the specific mechanisms engaged in their repair strongly influence the resulting patterns of chromosome rearrangement. Here, we review recent progress in understanding how certain distinctive features of the cancer genome, including clustered mutagenesis, tandem segmental duplications, complex breakpoints, chromothripsis, chromoplexy and chromoanasynthesis may arise.

Elucidating the molecular architecture of adaptation via evolve and resequence experiments

Evolve and resequence (E&R) experiments use experimental evolution to adapt populations to a novel environment, then next-generation sequencing to analyse genetic changes. They enable molecular evolution to be monitored in real time on a genome-wide scale. Here, we review the field of E&R experiments across diverse systems, ranging from simple non-living RNA to bacteria, yeast and the complex multicellular organism Drosophila melanogaster. We explore how different evolutionary outcomes in these systems are largely consistent with common population genetics principles. Differences in outcomes across systems are largely explained by different starting population sizes, levels of pre-existing genetic variation, recombination rates and adaptive landscapes. We highlight emerging themes and inconsistencies that future experiments must address.

Dosage compensation can buffer copy-number variation in wild yeast

Aneuploidy is linked to myriad diseases but also facilitates organismal evolution. It remains unclear how cells overcome the deleterious effects of aneuploidy until new phenotypes evolve. Although laboratory strains are extremely sensitive to aneuploidy, we show here that aneuploidy is common in wild yeast isolates, which show lower-than-expected expression at many amplified genes. We generated diploid strain panels in which cells carried two, three, or four copies of the affected chromosomes, to show that gene-dosage compensation functions at >30% of amplified genes. Genes subject to dosage compensation are under higher expression constraint in wild populations—but they show elevated rates of gene amplification, suggesting that copy-number variation is buffered at these genes. We find that aneuploidy provides a clear ecological advantage to oak strain YPS1009, by amplifying a causal gene that escapes dosage compensation. Our work presents a model in which dosage compensation buffers gene amplification through aneuploidy to provide a natural, but likely transient, route to rapid phenotypic evolution.

DOI:http://dx.doi.org/10.7554/eLife.05462.001

Evolution is driven by changes to the genes and other genetic information found in the DNA of an organism. These changes might, for example, alter the physical characteristics of the organism, or change how efficiently crucial tasks are carried out inside cells. Whatever the change, if it makes it easier for the organism to survive and reproduce, it is more likely to be passed on to future generations.

DNA is organized inside cells in structures called chromosomes. Most of the cells in animals, plants, and fungi contain two copies of each chromosome. However, sometimes mistakes happen during cell division and extra copies of a chromosome—and hence the genes contained within it—may end up in a cell. These extra copies of genes might help to speed up the rate at which a species evolves, as the ‘spare’ copies are free to adapt to new roles. However, having extra copies of genes can also often be harmful, and in humans can cause genetic disorders such as Down syndrome.

In the laboratory, chromosomes are commonly studied in a species of yeast called Saccharomyces cerevisiae. This species consists of several groups—or strains—that are genetically distinct from each other. Over the years, breeding the yeast for experiments has created laboratory strains that have lost some of the characteristics seen in wild strains. Earlier studies suggested that these cells fail to grow properly if they contain extra copies of chromosomes. Now, Hose et al. have studied nearly 50 wild strains of Saccharomyces cerevisiae. In these, extra copies of chromosomes are commonplace, and seemingly have no detrimental effect on growth. Instead, Hose et al. found that cells with too many copies of a gene use many of those genes less often than would be expected. This process is known as ‘dosage compensation’. This dosage compensation has not been observed in laboratory strains, in part because the extra gene copies make them sickly and hard to study.

Together, the results provide examples of how dosage compensation could help new traits to evolve in a species by reducing the negative effects of duplicated genes. This knowledge may have broad application, from suggesting methods to alleviate human disorders to implicating new ways to engineer useful traits in yeast and other microbes.

DOI:http://dx.doi.org/10.7554/eLife.05462.002

The diversity of fungal genome

The genome size of an organism varies from species to species. The C-value paradox enigma is a very complex puzzle with regards to vast diversity in genome sizes in eukaryotes. Here we reported the detailed genomic information of 172 fungal species among different fungal genomes and found that fungal genomes are very diverse in nature. In fungi, the diversity of genomes varies from 8.97 Mb to 177.57 Mb. The average genome sizes of Ascomycota and Basidiomycota fungi are 36.91 and 46.48 Mb respectively. But higher genome size is observed in Oomycota (74.85 Mb) species, a lineage of fungus-like eukaryotic microorganisms. The average coding genes of Oomycota species are almost doubled than that of Acomycota and Basidiomycota fungus.

The Ty1 LTR-retrotransposon of budding yeast, Saccharomyces cerevisiae

Long-terminal repeat (LTR)-retrotransposons generate a copy of their DNA (cDNA) by reverse transcription of their RNA genome in cytoplasmic nucleocapsids. They are widespread in the eukaryotic kingdom and are the evolutionary progenitors of retroviruses [1]. The Ty1 element of the budding yeast Saccharomyces cerevisiae was the first LTR-retrotransposon demonstrated to mobilize through an RNA intermediate, and not surprisingly, is the best studied. The depth of our knowledge of Ty1 biology stems not only from the predominance of active Ty1 elements in the S. cerevisiae genome but also the ease and breadth of genomic, biochemical and cell biology approaches available to study cellular processes in yeast. This review describes the basic structure of Ty1 and its gene products, the replication cycle, the rapidly expanding compendium of host co-factors known to influence retrotransposition and the nature of Ty1's elaborate symbiosis with its host. Our goal is to illuminate the value of Ty1 as a paradigm to explore the biology of LTR-retrotransposons in multicellular organisms, where the low frequency of retrotransposition events presents a formidable barrier to investigations of retrotransposon biology.

Macrotene chromosomes provide insights to a new mechanism of high-order gene amplification in eukaryotes

Copy number variation of chromosomal segments is now recognized as a major source of genetic polymorphism within natural populations of eukaryotes, as well as a possible cause of genetic diseases in humans, including cancer, but its molecular bases remain incompletely understood. In the baker's yeast Saccharomyces cerevisiae, a variety of low-order amplifications (segmental duplications) were observed after adaptation to limiting environmental conditions or recovery from gene dosage imbalance, and interpreted in terms of replication-based mechanisms associated or not with homologous recombination. Here we show the emergence of novel high-order amplification structures, with corresponding overexpression of embedded genes, during evolution under favourable growth conditions of severely unfit yeast cells bearing genetically disabled genomes. Such events form massively extended chromosomes, which we propose to call macrotene, whose characteristics suggest the products of intrachromosomal rolling-circle type of replication structures, probably initiated by increased accidental template switches under important cellular stress conditions.

Metagenomic chromosome conformation capture (meta3C) unveils the diversity of chromosome organization in microorganisms

Genomic analyses of microbial populations in their natural environment remain limited by the difficulty to assemble full genomes of individual species. Consequently, the chromosome organization of microorganisms has been investigated in a few model species, but the extent to which the features described can be generalized to other taxa remains unknown. Using controlled mixes of bacterial and yeast species, we developed meta3C, a metagenomic chromosome conformation capture approach that allows characterizing individual genomes and their average organization within a mix of organisms. Not only can meta3C be applied to species already sequenced, but a single meta3C library can be used for assembling, scaffolding and characterizing the tridimensional organization of unknown genomes. By applying meta3C to a semi-complex environmental sample, we confirmed its promising potential. Overall, this first meta3C study highlights the remarkable diversity of microorganisms chromosome organization, while providing an elegant and integrated approach to metagenomic analysis.

DOI:http://dx.doi.org/10.7554/eLife.03318.001

Microbial communities play vital roles in the environment and sustain animal and plant life. Marine microbes are part of the ocean's food chain; soil microbes support the turnover of major nutrients and facilitate plant growth; and the microbial communities residing in the human gut support digestion and the immune system, among other roles. These communities are very complex systems, often containing 1000s of different species engaged in co-dependent relationships, and are therefore very difficult to study.

The entire DNA sequence of an organism constitutes its genome, and much of this genetic information is stored in large structures called chromosomes. Examining the genome of a species can provide important clues about its lifestyle and how it evolved. To do this, DNA is extracted from cells and is then usually cut into smaller fragments, amplified, and sequenced. The small stretches of sequence obtained, called reads, are finally assembled, yielding ideally the complete genome of the organism under study.

Metagenomics attempts to interpret the combined genome of all the different species in a microbial community and has been instrumental in deciphering how the different species interact with each other. Metagenomics involves sequencing stretches of the community's DNA and matching these pieces to individual species to ultimately assemble whole genomes. While this may be a relatively straightforward task for communities that contain only a handful of members, the metagenomes derived from complex microbial communities are huge, fragmented, and incomplete. This often makes it very difficult or even nearly impossible to match the inferred DNA stretches to individual species.

A method called chromosome conformation capture (or ‘3C’ for short) can reveal the physical contacts between different regions of a chromosome and between the different chromosomes of a cell. How often each of these chromosomal contacts occurs provides a kind of physical signature to each genome and each individual chromosome within it.

Marbouty et al. took advantage of these interactions to develop a technique that combines metagenomics and chromosome conformation capture—called meta3C—that can analyze the DNA of many different species mixed together. Testing meta3C on artificial mixtures of a few species of yeast or bacteria showed that meta3C can separate the genomes of the different species without any prior knowledge of the composition of the mix. In a single experiment, meta3C can identify individual chromosomes, match each of them to its species of origin, and reveal the three-dimensional structure of each genome in the mix. Further tests showed that meta3C can also interpret more complex communities where the number and types of the species present are not known.

Meta3C holds great promise for understanding how microbial communities work and how the genomes of the species within a community are organized. However, further developments of the technique will be required to investigate communities as diverse as those present in most natural environments.

DOI:http://dx.doi.org/10.7554/eLife.03318.002

High-quality genome (re)assembly using chromosomal contact data

Closing gaps in draft genome assemblies can be costly and time-consuming, and published genomes are therefore often left ‘unfinished.’ Here we show that genome-wide chromosome conformation capture (3C) data can be used to overcome these limitations, and present a computational approach rooted in polymer physics that determines the most likely genome structure using chromosomal contact data. This algorithm—named GRAAL—generates high-quality assemblies of genomes in which repeated and duplicated regions are accurately represented and offers a direct probabilistic interpretation of the computed structures. We first validated GRAAL on the reference genome of Saccharomyces cerevisiae, as well as other yeast isolates, where GRAAL recovered both known and unknown complex chromosomal structural variations. We then applied GRAAL to the finishing of the assembly of Trichoderma reesei and obtained a number of contigs congruent with the know karyotype of this species. Finally, we showed that GRAAL can accurately reconstruct human chromosomes from either fragments generated in silico or contigs obtained from de novo assembly. In all these applications, GRAAL compared favourably to recently published programmes implementing related approaches.

The correct assembly of genomes from sequencing data remains a challenge due to difficulties in correctly assigning the location of repeated DNA elements. Here the authors describe GRAAL, an algorithm that utilizes genome-wide chromosome contact data within a probabilistic framework to produce accurate genome assemblies.

Evolutionary Perspectives on Human Fungal Pathogens

Owing to their small size and paucity of phenotypic characters, progress in the evolutionary biology of microbes in general, and human pathogenic fungi in particular, has been linked to a series of advances in DNA sequencing over the past quarter century. Phylogenetics was the first area to benefit, with the achievement of a basic understanding of fungal phylogeny. Population genetics was the next advance, finding cryptic species everywhere, and recombination in species previously thought to be asexual. Comparative genomics saw the next advance, in which variation in gene content and changes in gene family size were found to be important sources of variation. Fungal population genomics is showing that gene flow among closely related populations and species provides yet another source of adaptive, genetic variation. Now, two means to associate genetic variation with phenotypic variation, "reverse ecology" for adaptive phenotypes, and genome-wide association of any phenotype, are letting evolutionary biology make a profound contribution to molecular developmental biology of pathogenic fungi.

Structures of naturally evolved CUP1 tandem arrays in yeast indicate that these arrays are generated by unequal nonhomologous recombination

An important issue in genome evolution is the mechanism by which tandem duplications are generated from single-copy genes. In the yeast Saccharomyces cerevisiae, most strains contain tandemly duplicated copies of CUP1, a gene that encodes a copper-binding metallothionein. By screening 101 natural isolates of S. cerevisiae, we identified five different types of CUP1-containing repeats, as well as strains that only had one copy of CUP1. A comparison of the DNA sequences of these strains indicates that the CUP1 tandem arrays were generated by unequal nonhomologous recombination events from strains that had one CUP1 gene.