Extraction and selection of high-molecular-weight DNA for long-read sequencing from Chlamydomonas reinhardtii.

Recent advances in long-read sequencing technologies have enabled the complete assembly of eukaryotic genomes from telomere to telomere by allowing repeated regions to be fully sequenced and assembled, thus filling the gaps left by previous short-read sequencing methods. Furthermore, long-read sequencing can also help characterizing structural variants, with applications in the fields of genome evolution or cancer genomics. For many organisms, the main bottleneck to sequence long reads remains the lack of robust methods to obtain high-molecular-weight (HMW) DNA. For this purpose, we developed an optimized protocol to extract DNA suitable for long-read sequencing from the unicellular green alga Chlamydomonas reinhardtii, based on CTAB/phenol extraction followed by a size selection step for long DNA molecules. We provide validation results for the extraction protocol, as well as statistics obtained with Oxford Nanopore Technologies sequencing.

Domestication signatures in the non-conventional yeast Lachancea cidri.

Evaluating domestication signatures beyond model organisms is essential for a thorough understanding of the genotype-phenotype relationship in wild and human-related environments. Structural variations (SVs) can significantly impact phenotypes playing an important role in the physiological adaptation of species to different niches, including during domestication. A detailed characterization of the fitness consequences of these genomic rearrangements, however, is still limited in non-model systems, largely due to the paucity of direct comparisons between domesticated and wild isolates. Here, we used a combination of sequencing strategies to explore major genomic rearrangements in a Lachancea cidri yeast strain isolated from cider (CBS2950) and compared them to those in eight wild isolates from primary forests. Genomic analysis revealed dozens of SVs, including a large reciprocal translocation (~16 kb and 500 kb) present in the cider strain, but absent from all wild strains. Interestingly, the number of SVs was higher relative to single-nucleotide polymorphisms in the cider strain, suggesting a significant role in the strain's phenotypic variation. The set of SVs identified directly impacts dozens of genes and likely underpins the greater fermentation performance in the L. cidri CBS2950. In addition, the large reciprocal translocation affects a proline permease (PUT4) regulatory region, resulting in higher PUT4 transcript levels, which agrees with higher ethanol tolerance, improved cell growth when using proline, and higher amino acid consumption during fermentation. These results suggest that SVs are responsible for the rapid physiological adaptation of yeast to a human-related environment and demonstrate the key contribution of SVs in adaptive fermentative traits in non-model species.IMPORTANCEThe exploration of domestication signatures associated with human-related environments has predominantly focused on studies conducted on model organisms, such as Saccharomyces cerevisiae, overlooking the potential for comparisons across other non-Saccharomyces species. In our research, employing a combination of long- and short-read data, we found domestication signatures in Lachancea cidri, a non-model species recently isolated from fermentative environments in cider in France. The significance of our study lies in the identification of large array of major genomic rearrangements in a cider strain compared to wild isolates, which underly several fermentative traits. These domestication signatures result from structural variants, which are likely responsible for the phenotypic differences between strains, providing a rapid path of adaptation to human-related environments.

Telomerase-independent survival leads to a mosaic of complex subtelomere rearrangements in Chlamydomonas reinhardtii.

Telomeres and subtelomeres, the genomic regions located at chromosome extremities, are essential for genome stability in eukaryotes. In the absence of the canonical maintenance mechanism provided by telomerase, telomere shortening induces genome instability. The landscape of the ensuing genome rearrangements is not accessible by short-read sequencing. Here, we leverage Oxford Nanopore Technologies long-read sequencing to survey the extensive repertoire of genome rearrangements in telomerase mutants of the model green microalga Chlamydomonas reinhardtii In telomerase-mutant strains grown for hundreds of generations, most chromosome extremities were capped by short telomere sequences that were either recruited de novo from other loci or maintained in a telomerase-independent manner. Other extremities did not end with telomeres but only with repeated subtelomeric sequences. The subtelomeric elements, including rDNA, were massively rearranged and involved in breakage-fusion-bridge cycles, translocations, recombinations, and chromosome circularization. These events were established progressively over time and displayed heterogeneity at the subpopulation level. New telomere-capped extremities composed of sequences originating from more internal genomic regions were associated with high DNA methylation, suggesting that de novo heterochromatin formation contributes to the restoration of chromosome end stability in C. reinhardtii The diversity of alternative strategies present in the same organism to maintain chromosome integrity and the variety of rearrangements found in telomerase mutants are remarkable, and illustrate genome plasticity at short timescales.

Telomere-to-telomere assemblies of 142 strains characterize the genome structural landscape in Saccharomyces cerevisiae.

Pangenomes provide access to an accurate representation of the genetic diversity of species, both in terms of sequence polymorphisms and structural variants (SVs). Here we generated the Saccharomyces cerevisiae Reference Assembly Panel (ScRAP) comprising reference-quality genomes for 142 strains representing the species' phylogenetic and ecological diversity. The ScRAP includes phased haplotype assemblies for several heterozygous diploid and polyploid isolates. We identified circa (ca.) 4,800 nonredundant SVs that provide a broad view of the genomic diversity, including the dynamics of telomere length and transposable elements. We uncovered frequent cases of complex aneuploidies where large chromosomes underwent large deletions and translocations. We found that SVs can impact gene expression near the breakpoints and substantially contribute to gene repertoire evolution. We also discovered that horizontally acquired regions insert at chromosome ends and can generate new telomeres. Overall, the ScRAP demonstrates the benefit of a pangenome in understanding genome evolution at population scale.

Late Pleistocene-dated divergence between South Hemisphere populations of the non-conventional yeast L. cidri.

Most organisms belonging to the Saccharomycotina subphylum have high genetic diversity and a vast repertoire of metabolisms and lifestyles. Lachancea cidri is an ideal yeast model for exploring the interplay between genetics, ecological function and evolution. Lachancea cidri diverged from the Saccharomyces lineage before the whole-genome duplication and is distributed across the South Hemisphere, displaying an important ecological success. We applied phylogenomics to investigate the genetic variation of L. cidri isolates obtained from Australia and South America. Our approach revealed the presence of two main lineages according to their geographic distribution (Aus and SoAm). Estimation of the divergence time suggests that SoAm and Aus lineages diverged near the last glacial maximum event during the Pleistocene (64-8 KYA). Interestingly, we found that the French reference strain is closely related to the Australian strains, with a recent divergence (405-51 YA), likely associated to human movements. Additionally, we identified different lineages within the South American population, revealing that Patagonia contains a similar genetic diversity comparable to that of other lineages in S. cerevisiae. These findings support the idea of a Pleistocene-dated divergence between South Hemisphere lineages, where the Nothofagus and Araucaria ecological niches likely favoured the extensive distribution of L. cidri in Patagonia.

An evolutionary model identifies the main evolutionary biases for the evolution of genome-replication profiles.

Recent results comparing the temporal program of genome replication of yeast species belonging to the Lachancea clade support the scenario that the evolution of the replication timing program could be mainly driven by correlated acquisition and loss events of active replication origins. Using these results as a benchmark, we develop an evolutionary model defined as birth-death process for replication origins and use it to identify the evolutionary biases that shape the replication timing profiles. Comparing different evolutionary models with data, we find that replication origin birth and death events are mainly driven by two evolutionary pressures, the first imposes that events leading to higher double-stall probability of replication forks are penalized, while the second makes less efficient origins more prone to evolutionary loss. This analysis provides an empirically grounded predictive framework for quantitative evolutionary studies of the replication timing program.

A Versatile Protocol to Generate Translocations in Yeast Genomes Using CRISPR/Cas9.

Genomic engineering methods represent powerful tools to examine chromosomal modifications and to subsequently study their impacts on cellular phenotypes. However, quantifying the fitness impact of translocations, independently from base substitutions or the insertion of genetic markers, remains a challenge. Here we report a rapid and straightforward protocol for engineering either targeted reciprocal translocations at the base pair level of resolution between two chromosomes or multiple simultaneous rearrangements in the yeast genome, without inserting any marker sequence in the chromosomes. Our CRISPR/Cas9-based method consists of inducing either (1) two double-strand breaks (DSBs) in two different chromosomes with two distinct guide RNAs (gRNAs) while providing specifically designed homologous donor DNA forcing the trans-repair of chromosomal extremities to generate a targeted reciprocal translocation or (2) multiple DSBs with a single gRNA targeting dispersed repeated sequences and leaving endogenous uncut copies of the repeat to be used as donor DNA, thereby generating multiple translocations, often associated with large segmental duplications (Fleiss, et al. PLoS Genet 15:e1008332, 2019).

Evidence for Two Main Domestication Trajectories in Saccharomyces cerevisiae Linked to Distinct Bread-Making Processes.

Production of leavened bread dates to the second millennium BCE. Since then, the art of bread making has developed, yet the evolution of bread-associated microbial species remains largely unknown. Nowadays, leavened bread is made either by using a pure commercial culture of the yeast Saccharomyces cerevisiae or by propagating a sourdough-a mix of flour and water spontaneously fermented by yeasts and bacteria. We studied the domestication of S. cerevisiae originating from industrial sources and artisanal sourdoughs and tested whether different bread-making processes led to population divergence. We found that S. cerevisiae bakery strains are polyphyletic with 67% of strains clustering into two main clades: most industrial strains were tetraploid and clustered with strains having diverse origins, including beer. By contrast, most sourdough strains were diploid and grouped in a second clade of strains having mosaic genomes and diverse origins, including fruits and natural environments. They harbored a higher copy number of genes involved in maltose utilization, and a high level of gene flow from multiple contributors was detected. Bakery strains displayed higher CO2 production than do strains from other domesticated lineages (such as beer and wine), revealing a specific phenotypic signature of domestication. Interestingly, industrial strains had a shorter fermentation onset than sourdough strains, which were better adapted to a sourdough-like environment, suggesting divergent selection by industrial and artisanal processes. Our results reveal that the domestication of bakery yeast has been accompanied by dispersion, hybridization, and divergent selection through industrial and artisanal processes.

Reshuffling yeast chromosomes with CRISPR/Cas9.

Genome engineering is a powerful approach to study how chromosomal architecture impacts phenotypes. However, quantifying the fitness impact of translocations independently from the confounding effect of base substitutions has so far remained challenging. We report a novel application of the CRISPR/Cas9 technology allowing to generate with high efficiency both uniquely targeted and multiple concomitant reciprocal translocations in the yeast genome. Targeted translocations are constructed by inducing two double-strand breaks on different chromosomes and forcing the trans-chromosomal repair through homologous recombination by chimerical donor DNAs. Multiple translocations are generated from the induction of several DSBs in LTR repeated sequences and promoting repair using endogenous uncut LTR copies as template. All engineered translocations are markerless and scarless. Targeted translocations are produced at base pair resolution and can be sequentially generated one after the other. Multiple translocations result in a large diversity of karyotypes and are associated in many instances with the formation of unanticipated segmental duplications. To test the phenotypic impact of translocations, we first recapitulated in a lab strain the SSU1/ECM34 translocation providing increased sulphite resistance to wine isolates. Surprisingly, the same translocation in a laboratory strain resulted in decreased sulphite resistance. However, adding the repeated sequences that are present in the SSU1 promoter of the resistant wine strain induced sulphite resistance in the lab strain, yet to a lower level than that of the wine isolate, implying that additional polymorphisms also contribute to the phenotype. These findings illustrate the advantage brought by our technique to untangle the phenotypic impacts of structural variations from confounding effects of base substitutions. Secondly, we showed that strains with multiple translocations, even those devoid of unanticipated segmental duplications, display large phenotypic diversity in a wide range of environmental conditions, showing that simply reconfiguring chromosome architecture is sufficient to provide fitness advantages in stressful growth conditions.

bHLH-PAS protein RITMO1 regulates diel biological rhythms in the marine diatom Phaeodactylum tricornutum.

Periodic light-dark cycles govern the timing of basic biological processes in organisms inhabiting land as well as the sea, where life evolved. Although prominent marine phytoplanktonic organisms such as diatoms show robust diel rhythms, the mechanisms regulating these processes are still obscure. By characterizing a Phaeodactylum tricornutum bHLH-PAS nuclear protein, hereby named RITMO1, we shed light on the regulation of the daily life of diatoms. Alteration of RITMO1 expression levels and timing by ectopic overexpression results in lines with deregulated diurnal gene expression profiles compared with the wild-type cells. Reduced gene expression oscillations are also observed in these lines in continuous darkness, showing that the regulation of rhythmicity by RITMO1 is not directly dependent on light inputs. We also describe strong diurnal rhythms of cellular fluorescence in wild-type cells, which persist in continuous light conditions, indicating the existence of an endogenous circadian clock in diatoms. The altered rhythmicity observed in RITMO1 overexpression lines in continuous light supports the involvement of this protein in circadian rhythm regulation. Phylogenetic analysis reveals a wide distribution of RITMO1-like proteins in the genomes of diatoms as well as in other marine algae, which may indicate a common function in these phototrophs. This study adds elements to our understanding of diatom biology and offers perspectives to elucidate timekeeping mechanisms in marine organisms belonging to a major, but under-investigated, branch of the tree of life.

Characterizing meiotic chromosomes' structure and pairing using a designer sequence optimized for Hi-C.

In chromosome conformation capture experiments (Hi-C), the accuracy with which contacts are detected varies due to the uneven distribution of restriction sites along genomes. In addition, repeated sequences or homologous regions remain indistinguishable because of the ambiguities they introduce during the alignment of the sequencing reads. We addressed both limitations by designing and engineering 144 kb of a yeast chromosome with regularly spaced restriction sites (Syn-HiC design). In the Syn-HiC region, Hi-C signal-to-noise ratio is enhanced and can be used to measure the shape of an unbiased distribution of contact frequencies, allowing to propose a robust definition of a Hi-C experiment resolution. The redesigned region is also distinguishable from its native homologous counterpart in an otherwise isogenic diploid strain. As a proof of principle, we tracked homologous chromosomes during meiotic prophase in synchronized and pachytene-arrested cells and captured important features of their spatial reorganization, such as chromatin restructuration into arrays of Rec8-delimited loops, centromere declustering, individualization, and pairing. Overall, we illustrate the promises held by redesigning genomic regions to explore complex biological questions.

The evolution of the temporal program of genome replication.

Genome replication is highly regulated in time and space, but the rules governing the remodeling of these programs during evolution remain largely unknown. We generated genome-wide replication timing profiles for ten Lachancea yeasts, covering a continuous evolutionary range from closely related to more divergent species. We show that replication programs primarily evolve through a highly dynamic evolutionary renewal of the cohort of active replication origins. We found that gained origins appear with low activity yet become more efficient and fire earlier as they evolutionarily age. By contrast, origins that are lost comprise the complete range of firing strength. Additionally, they preferentially occur in close vicinity to strong origins. Interestingly, despite high evolutionary turnover, active replication origins remain regularly spaced along chromosomes in all species, suggesting that origin distribution is optimized to limit large inter-origin intervals. We propose a model on the evolutionary birth, death, and conservation of active replication origins.

Reconstruction of ancestral chromosome architecture and gene repertoire reveals principles of genome evolution in a model yeast genus.

Reconstructing genome history is complex but necessary to reveal quantitative principles governing genome evolution. Such reconstruction requires recapitulating into a single evolutionary framework the evolution of genome architecture and gene repertoire. Here, we reconstructed the genome history of the genus Lachancea that appeared to cover a continuous evolutionary range from closely related to more diverged yeast species. Our approach integrated the generation of a high-quality genome data set; the development of AnChro, a new algorithm for reconstructing ancestral genome architecture; and a comprehensive analysis of gene repertoire evolution. We found that the ancestral genome of the genus Lachancea contained eight chromosomes and about 5173 protein-coding genes. Moreover, we characterized 24 horizontal gene transfers and 159 putative gene creation events that punctuated species diversification. We retraced all chromosomal rearrangements, including gene losses, gene duplications, chromosomal inversions and translocations at single gene resolution. Gene duplications outnumbered losses and balanced rearrangements with 1503, 929, and 423 events, respectively. Gene content variations between extant species are mainly driven by differential gene losses, while gene duplications remained globally constant in all lineages. Remarkably, we discovered that balanced chromosomal rearrangements could be responsible for up to 14% of all gene losses by disrupting genes at their breakpoints. Finally, we found that nonsynonymous substitutions reached fixation at a coordinated pace with chromosomal inversions, translocations, and duplications, but not deletions. Overall, we provide a granular view of genome evolution within an entire eukaryotic genus, linking gene content, chromosome rearrangements, and protein divergence into a single evolutionary framework.

A Versatile Procedure to Generate Genome-Wide Spatiotemporal Program of Replication in Yeast Species.

Here, we describe a complete protocol, comprising both the experimental and the analytical procedures, that allows to generate genome-wide spatiotemporal program of replication and to find the location of chromosomally active replication origins in yeast. The first step consists on synchronizing a cell population by physical discrimination of G1 cells according to their sedimentation coefficient. G1 cells are then synchronously released into S-phase and time-point samples are regularly taken until they reach the G2 phase. Progression through the cell cycle is monitored by measuring DNA content variation by flow cytometry. DNA samples, covering the entire S-phase, are then extracted and analyzed using deep sequencing. The gradual change of DNA copy number is measured to determine the mean replication time along the genome. A simple method of peak calling allows to infer from the replication profile the location of replication origins along the chromosomes. Our protocol is versatile enough to be applied to virtually any yeast species of interest and generate its replication profile.

An integrated "omics" approach to the characterization of maize (Zea mays L.) mutants deficient in the expression of two genes encoding cytosolic glutamine synthetase.

BACKGROUND: To identify the key elements controlling grain production in maize, it is essential to have an integrated view of the responses to alterations in the main steps of nitrogen assimilation by modification of gene expression. Two maize mutant lines (gln1.3 and gln1.4), deficient in two genes encoding cytosolic glutamine synthetase, a key enzyme involved in nitrogen assimilation, were previously characterized by a reduction of kernel size in the gln1.4 mutant and by a reduction of kernel number in the gln1.3 mutant. In this work, the differences in leaf gene transcripts, proteins and metabolite accumulation in gln1.3 and gln1.4 mutants were studied at two key stages of plant development, in order to identify putative candidate genes, proteins and metabolic pathways contributing on one hand to the control of plant development and on the other to grain production. RESULTS: The most interesting finding in this study is that a number of key plant processes were altered in the gln1.3 and gln1.4 mutants, including a number of major biological processes such as carbon metabolism and transport, cell wall metabolism, and several metabolic pathways and stress responsive and regulatory elements. We also found that the two mutants share common or specific characteristics across at least two or even three of the "omics" considered at the vegetative stage of plant development, or during the grain filling period. CONCLUSIONS: This is the first comprehensive molecular and physiological characterization of two cytosolic glutamine synthetase maize mutants using a combined transcriptomic, proteomic and metabolomic approach. We find that the integration of the three "omics" procedures is not straight forward, since developmental and mutant-specific levels of regulation seem to occur from gene expression to metabolite accumulation. However, their potential use is discussed with a view to improving our understanding of nitrogen assimilation and partitioning and its impact on grain production.

Total synthesis of a functional designer eukaryotic chromosome.

Rapid advances in DNA synthesis techniques have made it possible to engineer viruses, biochemical pathways and assemble bacterial genomes. Here, we report the synthesis of a functional 272,871-base pair designer eukaryotic chromosome, synIII, which is based on the 316,617-base pair native Saccharomyces cerevisiae chromosome III. Changes to synIII include TAG/TAA stop-codon replacements, deletion of subtelomeric regions, introns, transfer RNAs, transposons, and silent mating loci as well as insertion of loxPsym sites to enable genome scrambling. SynIII is functional in S. cerevisiae. Scrambling of the chromosome in a heterozygous diploid reveals a large increase in a-mater derivatives resulting from loss of the MATα allele on synIII. The complete design and synthesis of synIII establishes S. cerevisiae as the basis for designer eukaryotic genome biology.

The spatiotemporal program of replication in the genome of Lachancea kluyveri.

We generated a genome-wide replication profile in the genome of Lachancea kluyveri and assessed the relationship between replication and base composition. This species diverged from Saccharomyces cerevisiae before the ancestral whole genome duplication. The genome comprises eight chromosomes among which a chromosomal arm of 1 Mb has a G + C-content much higher than the rest of the genome. We identified 252 active replication origins in L. kluyveri and found considerable divergence in origin location with S. cerevisiae and with Lachancea waltii. Although some global features of S. cerevisiae replication are conserved: Centromeres replicate early, whereas telomeres replicate late, we found that replication origins both in L. kluyveri and L. waltii do not behave as evolutionary fragile sites. In L. kluyveri, replication timing along chromosomes alternates between regions of early and late activating origins, except for the 1 Mb GC-rich chromosomal arm. This chromosomal arm contains an origin consensus motif different from other chromosomes and is replicated early during S-phase. We showed that precocious replication results from the specific absence of late firing origins in this chromosomal arm. In addition, we found a correlation between GC-content and distance from replication origins as well as a lack of replication-associated compositional skew between leading and lagging strands specifically in this GC-rich chromosomal arm. These findings suggest that the unusual base composition in the genome of L. kluyveri could be linked to replication.

A head and neck cancer tumor response-specific gene signature for cisplatin, 5-fluorouracil induction chemotherapy fails with added taxanes.

BACKGROUND: It is a major clinical challenge to predict which patients, with advanced stage head and neck squamous cell carcinoma, will not exhibit a reduction in tumor size following induction chemotherapy in order to avoid toxic effects of ineffective chemotherapy and delays for instituting other therapeutic options. Further, it is of interest to know to what extent a gene signature, which identifies patients with tumors that will not respond to a particular induction chemotherapy, is applicable when additional chemotherapeutic agents are added to the regimen. METHODOLOGY/PRINCIPAL FINDINGS: To identify genes that predict tumor resistance to induction with cisplatin/5-fluorouracil (PF) or PF and a taxane, we analyzed patient tumor biopsies with whole genome microarrays and quantitative reverse transcriptase-PCR (TLDA) cards. A leave one out cross-validation procedure allowed evaluation of the prediction tool. A ten-gene microarray signature correctly classified 12/13 responders and 7/10 non-responders to PF (92% specificity, 82.6% accuracy). TLDA analysis (using the same classifier) of the patients correctly classified 12/12 responders and 8/10 non-responders (100% specificity, 90.9% accuracy). Further, TLDA analysis correctly predicted the response of 5 new patients and, overall, 12/12 responders and 13/15 non-responders (100% specificity, 92.6% accuracy). The protein products of the genes constituting the signature physically associate with 27 other proteins, involved in regulating gene expression, constituting an interaction network. In contrast, TLDA-based prediction (with the same gene signature) of responses to induction with PF and either of two taxanes was poor (0% specificity, 25% accuracy and 33.3% specificity, 25% accuracy). CONCLUSIONS/SIGNIFICANCE: Successful transfer of the microarray-based gene signature to an independent, PCR-based technology suggests that TLDA-based signatures could be a useful hospital-based technology for determining therapeutic options. Although highly specific for tumor responses to PF induction, the gene signature is unsuccessful when taxanes are added. The results illustrate the subtlety in developing "personalized medicine".

The use of metabolomics integrated with transcriptomic and proteomic studies for identifying key steps involved in the control of nitrogen metabolism in crops such as maize.

Linking plant phenotype to gene and protein expression and also to metabolite synthesis and accumulation is one of the main challenges for improving agricultural production worldwide. Such a challenge is particularly relevant to crop nitrogen use efficiency (NUE). Here, the differences in leaf gene transcript, protein, and metabolite accumulation in maize subjected to long-term nitrogen (N)-deficient growth conditions at two important stages of plant development have been studied. The impact of N deficiency was examined at the transcriptomic, proteomic, and metabolomic levels. It was found that a number of key plant biological functions were either up- or down-regulated when N was limiting, including major alterations to photosynthesis, carbon (C) metabolism, and, to a lesser extent, downstream metabolic pathways. It was also found that the impact of the N deficiency stress resembled the response of plants to a number of other biotic and abiotic stresses, in terms of transcript, protein, and metabolite accumulation. The genetic and metabolic alterations were different during the N assimilation and the grain-filling period, indicating that plant development is an important component for identifying the key elements involved in the control of plant NUE. It was also found that integration of the three 'omics' studies is not straightforward, since different levels of regulation seem to occur in a stepwise manner from gene expression to metabolite accumulation. The potential use of these 'omics' studies is discussed with a view to improve our understanding of whole plant nitrogen economics, which should have applications in breeding and agronomy.