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.

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.

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.

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.

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.

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.

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.

The mutational profile of the yeast genome is shaped by replication.

Despite the scrutiny that has been directed for years at the yeast genome, relatively little is known about the impact of replication on the substitution dynamics in Saccharomyces cerevisiae. Here, we show that the mutation rate increases with the replication timing by more than 30% between the earliest and the latest replicating regions. In addition, we found a mutational asymmetry associated with the polarity of replication resulting in higher rates of substitutions toward C and A than toward G and T in leading strands (reciprocally more substitutions toward G and T in lagging strands). Such mutational asymmetries applied over long evolutionary periods should generate compositional skews between the two DNA strands. Thus, we show that the leading replicating strands present an excess of C over G and of A over T in the genome of S. cerevisiae (reciprocally an excess of G + T over C + A in lagging strands). We also show that the nucleotide frequencies at mutational equilibrium predict a compositional skew at equilibrium very close to the observed skew between leading and lagging strands, suggesting that compositional equilibrium has been nearly attained in the present day genome of S. cerevisiae. Surprisingly, the direction of this skew is inverted compared with the one in the human genome.

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.

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.

Functional study of genes essential for autogamy and nuclear reorganization in Paramecium.

Like all ciliates, Paramecium tetraurelia is a unicellular eukaryote that harbors two kinds of nuclei within its cytoplasm. At each sexual cycle, a new somatic macronucleus (MAC) develops from the germ line micronucleus (MIC) through a sequence of complex events, which includes meiosis, karyogamy, and assembly of the MAC genome from MIC sequences. The latter process involves developmentally programmed genome rearrangements controlled by noncoding RNAs and a specialized RNA interference machinery. We describe our first attempts to identify genes and biological processes that contribute to the progression of the sexual cycle. Given the high percentage of unknown genes annotated in the P. tetraurelia genome, we applied a global strategy to monitor gene expression profiles during autogamy, a self-fertilization process. We focused this pilot study on the genes carried by the largest somatic chromosome and designed dedicated DNA arrays covering 484 genes from this chromosome (1.2% of all genes annotated in the genome). Transcriptome analysis revealed four major patterns of gene expression, including two successive waves of gene induction. Functional analysis of 15 upregulated genes revealed four that are essential for vegetative growth, one of which is involved in the maintenance of MAC integrity and another in cell division or membrane trafficking. Two additional genes, encoding a MIC-specific protein and a putative RNA helicase localizing to the old and then to the new MAC, are specifically required during sexual processes. Our work provides a proof of principle that genes essential for meiosis and nuclear reorganization can be uncovered following genome-wide transcriptome analysis.

Assessing the impact of transgenerational epigenetic variation on complex traits.

Loss or gain of DNA methylation can affect gene expression and is sometimes transmitted across generations. Such epigenetic alterations are thus a possible source of heritable phenotypic variation in the absence of DNA sequence change. However, attempts to assess the prevalence of stable epigenetic variation in natural and experimental populations and to quantify its impact on complex traits have been hampered by the confounding effects of DNA sequence polymorphisms. To overcome this problem as much as possible, two parents with little DNA sequence differences, but contrasting DNA methylation profiles, were used to derive a panel of epigenetic Recombinant Inbred Lines (epiRILs) in the reference plant Arabidopsis thaliana. The epiRILs showed variation and high heritability for flowering time and plant height ( approximately 30%), as well as stable inheritance of multiple parental DNA methylation variants (epialleles) over at least eight generations. These findings provide a first rationale to identify epiallelic variants that contribute to heritable variation in complex traits using linkage or association studies. More generally, the demonstration that numerous epialleles across the genome can be stable over many generations in the absence of selection or extensive DNA sequence variation highlights the need to integrate epigenetic information into population genetics studies.

DYRK1A interacts with the REST/NRSF-SWI/SNF chromatin remodelling complex to deregulate gene clusters involved in the neuronal phenotypic traits of Down syndrome.

The molecular mechanisms that lead to the cognitive defects characteristic of Down syndrome (DS), the most frequent cause of mental retardation, have remained elusive. Here we use a transgenic DS mouse model (152F7 line) to show that DYRK1A gene dosage imbalance deregulates chromosomal clusters of genes located near neuron-restrictive silencer factor (REST/NRSF) binding sites. We found that Dyrk1a binds the SWI/SNF complex known to interact with REST/NRSF. The mutation of a REST/NRSF binding site in the promoter of the REST/NRSF target gene L1cam modifies the transcriptional effect of Dyrk1a-dosage imbalance on L1cam. Dyrk1a dosage imbalance perturbs Rest/Nrsf levels with decreased Rest/Nrsf expression in embryonic neurons and increased expression in adult neurons. Using transcriptome analysis of embryonic brain subregions of transgenic 152F7 mouse line, we identified a coordinated deregulation of multiple genes that are responsible for dendritic growth impairment present in DS. Similarly, Dyrk1a overexpression in primary mouse cortical neurons induced severe reduction of the dendritic growth and dendritic complexity. We propose that DYRK1A overexpression-related neuronal gene deregulation via disturbance of REST/NRSF levels, and the REST/NRSF-SWI/SNF chromatin remodelling complex, significantly contributes to the neural phenotypic changes that characterize DS.

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).

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.

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.

Epigenetic natural variation in Arabidopsis thaliana.

Cytosine methylation of repetitive sequences is widespread in plant genomes, occurring in both symmetric (CpG and CpNpG) as well as asymmetric sequence contexts. We used the methylation-dependent restriction enzyme McrBC to profile methylated DNA using tiling microarrays of Arabidopsis Chromosome 4 in two distinct ecotypes, Columbia and Landsberg erecta. We also used comparative genome hybridization to profile copy number polymorphisms. Repeated sequences and transposable elements (TEs), especially long terminal repeat retrotransposons, are densely methylated, but one third of genes also have low but detectable methylation in their transcribed regions. While TEs are almost always methylated, genic methylation is highly polymorphic, with half of all methylated genes being methylated in only one of the two ecotypes. A survey of loci in 96 Arabidopsis accessions revealed a similar degree of methylation polymorphism. Within-gene methylation is heritable, but is lost at a high frequency in segregating F(2) families. Promoter methylation is rare, and gene expression is not generally affected by differences in DNA methylation. Small interfering RNA are preferentially associated with methylated TEs, but not with methylated genes, indicating that most genic methylation is not guided by small interfering RNA. This may account for the instability of gene methylation, if occasional failure of maintenance methylation cannot be restored by other means.