A deep learning-based method enables the automatic and accurate assembly of chromosome-level genomes

The application of high-throughput chromosome conformation capture (Hi-C) technology enables the construction of chromosome-level assemblies. However, the correction of errors and the anchoring of sequences to chromosomes in the assembly remain significant challenges. In this study, we developed a deep learning-based method, AutoHiC, to address the challenges in chromosome-level genome assembly by enhancing contiguity and accuracy. Conventional Hi-C-aided scaffolding often requires manual refinement, but AutoHiC instead utilizes Hi-C data for automated workflows and iterative error correction. When trained on data from 300+ species, AutoHiC demonstrated a robust average error detection accuracy exceeding 90%. The benchmarking results confirmed its significant impact on genome contiguity and error correction. The innovative approach and comprehensive results of AutoHiC constitute a breakthrough in automated error detection, promising more accurate genome assemblies for advancing genomics research.

Fundamental Patterns of Structural Evolution Revealed by Chromosome-Length Genomes of Cactophilic Drosophila

A thorough understanding of adaptation and speciation requires model organisms with both a history of ecological and phenotypic study as well as a complete set of genomic resources. In particular, high-quality genome assemblies of ecological model organisms are needed to assess the evolution of genome structure and its role in adaptation and speciation. Here, we generate new genomes of cactophilic Drosophila, a crucial model clade for understanding speciation and ecological adaptation in xeric environments. We generated chromosome-level genome assemblies and complete annotations for seven populations across Drosophila mojavensis, Drosophila arizonae, and Drosophila navojoa. We use these data first to establish the most robust phylogeny for this clade to date, and to assess patterns of molecular evolution across the phylogeny, showing concordance with a priori hypotheses regarding adaptive genes in this system. We then show that structural evolution occurs at constant rate across the phylogeny, varies by chromosome, and is correlated with molecular evolution. These results advance the understanding of the D. mojavensis clade by demonstrating core evolutionary genetic patterns and integrating those patterns to generate new gene-level hypotheses regarding adaptation. Our data are presented in a new public database (cactusflybase.arizona.edu), providing one of the most in-depth resources for the analysis of inter- and intraspecific evolutionary genomic data. Furthermore, we anticipate that the patterns of structural evolution identified here will serve as a baseline for future comparative studies to identify the factors that influence the evolution of genome structure across taxa.

An integrative taxonomy approach reveals Saccharomyces chiloensis sp. nov. as a newly discovered species from Coastal Patagonia

Species delineation in microorganisms is challenging due to the limited markers available for accurate species assignment. Here, we applied an integrative taxonomy approach, combining extensive sampling, whole-genome sequence-based classification, phenotypic profiling, and assessment of interspecific reproductive isolation. Our work reveals the presence of a distinct Saccharomyces lineage in Nothofagus forests of coastal Patagonia. This lineage, designated Saccharomyces chiloensis sp. nov., exhibits 7% genetic divergence from its sister species S. uvarum, as revealed by whole-genome sequencing and population analyses. The South America-C (SA-C) coastal Patagonia population forms a unique clade closely related to a previously described divergent S. uvarum population from Oceania (AUS, found in Australia and New Zealand). Our species reclassification is supported by a low Ortho Average Nucleotide Identity (OANI) of 93% in SA-C and AUS relative to S. uvarum, which falls below the suggested species delineation threshold of 95%, indicating an independent evolutionary lineage. Hybrid spore viability assessment provided compelling evidence that SA-C and AUS are reproductively isolated from S. uvarum. In addition, we found unique structural variants between S. chiloensis sp. nov. lineages, including large-scale chromosomal translocations and inversions, together with a distinct phenotypic profile, emphasizing their intraspecies genetic distinctiveness. We suggest that S. chiloensis sp. nov diverged from S. uvarum in allopatry due to glaciation, followed by post-glacial dispersal, resulting in distinct lineages on opposite sides of the Pacific Ocean. The discovery of S. chiloensis sp. nov. illustrates the uniqueness of Patagonia's coastal biodiversity and underscores the importance of adopting an integrative taxonomic approach in species delineation to unveil cryptic microbial species. The holotype of S. chiloensis sp. nov. is CBS 18620T.

Long-read sequencing reveals extensive gut phageome structural variations driven by genetic exchange with bacterial hosts

Genetic variations are instrumental for unraveling phage evolution and deciphering their functional implications. Here, we explore the underlying fine-scale genetic variations in the gut phageome, especially structural variations (SVs). By using virome-enriched long-read metagenomic sequencing across 91 individuals, we identified a total of 14,438 nonredundant phage SVs and revealed their prevalence within the human gut phageome. These SVs are mainly enriched in genes involved in recombination, DNA methylation, and antibiotic resistance. Notably, a substantial fraction of phage SV sequences share close homology with bacterial fragments, with most SVs enriched for horizontal gene transfer (HGT) mechanism. Further investigations showed that these SV sequences were genetic exchanged between specific phage-bacteria pairs, particularly between phages and their respective bacterial hosts. Temperate phages exhibit a higher frequency of genetic exchange with bacterial chromosomes and then virulent phages. Collectively, our findings provide insights into the genetic landscape of the human gut phageome.

Reference-free structural variant detection in microbiomes via long-read co-assembly graphs

MOTIVATION

The study of bacterial genome dynamics is vital for understanding the mechanisms underlying microbial adaptation, growth, and their impact on host phenotype. Structural variants (SVs), genomic alterations of 50 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to the absence of clear reference genomes and the presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing all metagenomic samples in a series (time or other metric) into a single co-assembly graph. The log fold change in graph coverage between successive samples is then calculated to call SVs that are thriving or declining.

RESULTS

We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, particularly as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between successive time and temperature samples, suggesting host advantage. Our approach leverages previous work in assembly graph structural and coverage patterns to provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial gene flux.

AVAILABILITY AND IMPLEMENTATION

rhea is open source and available at: https://github.com/treangenlab/rhea.

Haplotype-resolved gapless genome and chromosome segment substitution lines facilitate gene identification in wild rice

The abundant genetic variation harbored by wild rice (Oryza rufipogon) has provided a reservoir of useful genes for rice breeding. However, the genome of wild rice has not yet been comprehensively assessed. Here, we report the haplotype-resolved gapless genome assembly and annotation of wild rice Y476. In addition, we develop two sets of chromosome segment substitution lines (CSSLs) using Y476 as the donor parent and cultivated rice as the recurrent parents. By analyzing the gapless reference genome and CSSL population, we identify 254 QTLs associated with agronomic traits, biotic and abiotic stresses. We clone a receptor-like kinase gene associated with rice blast resistance and confirm its wild rice allele improves rice blast resistance. Collectively, our study provides a haplotype-resolved gapless reference genome and demonstrates a highly efficient platform for gene identification from wild rice.

Diallel panel reveals a significant impact of low-frequency genetic variants on gene expression variation in yeast

Unraveling the genetic sources of gene expression variation is essential to better understand the origins of phenotypic diversity in natural populations. Genome-wide association studies identified thousands of variants involved in gene expression variation, however, variants detected only explain part of the heritability. In fact, variants such as low-frequency and structural variants (SVs) are poorly captured in association studies. To assess the impact of these variants on gene expression variation, we explored a half-diallel panel composed of 323 hybrids originated from pairwise crosses of 26 natural Saccharomyces cerevisiae isolates. Using short- and long-read sequencing strategies, we established an exhaustive catalog of single nucleotide polymorphisms (SNPs) and SVs for this panel. Combining this dataset with the transcriptomes of all hybrids, we comprehensively mapped SNPs and SVs associated with gene expression variation. While SVs impact gene expression variation, SNPs exhibit a higher effect size with an overrepresentation of low-frequency variants compared to common ones. These results reinforce the importance of dissecting the heritability of complex traits with a comprehensive catalog of genetic variants at the population level.

Major impacts of widespread structural variation on sorghum

Genetic diversity is critical to crop breeding and improvement, and dissection of the genomic variation underlying agronomic traits can both assist breeding and give insight into basic biological mechanisms. Although recent genome analyses in plants reveal many structural variants (SVs), most current studies of crop genetic variation are dominated by single-nucleotide polymorphisms (SNPs). The extent of the impact of SVs on global trait variation, as well as their utility in genome-wide selection, is not yet understood. In this study, we built an SV data set based on whole-genome resequencing of diverse sorghum lines (n = 363), validated the correlation of photoperiod sensitivity and variety type, and identified SV hotspots underlying the divergent evolution of cellulosic and sweet sorghum. In addition, we showed the complementary contribution of SVs for heritability of traits related to sorghum adaptation. Importantly, inclusion of SV polymorphisms in association studies revealed genotype-phenotype associations not observed with SNPs alone. Three-way genome-wide association studies (GWAS) based on whole-genome SNP, SV, and integrated SNP + SV data sets showed substantial associations between SVs and sorghum traits. The addition of SVs to GWAS substantially increased heritability estimates for some traits, indicating their important contribution to functional allelic variation at the genome level. Our discovery of the widespread impacts of SVs on heritable gene expression variation could render a plausible mechanism for their disproportionate impact on phenotypic variation. This study expands our knowledge of SVs and emphasizes the extensive impacts of SVs on sorghum.

Reference-free Structural Variant Detection in Microbiomes via Long-read Coassembly Graphs

Bacterial genome dynamics are vital for understanding the mechanisms underlying microbial adaptation, growth, and their broader impact on host phenotype. Structural variants (SVs), genomic alterations of 10 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to absence of clear reference genomes and presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing a single metagenome coassembly graph constructed from all samples in a series. The log fold change in graph coverage between subsequent samples is then calculated to call SVs that are thriving or declining throughout the series. We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, which is particularly noticeable as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between subsequent time and temperature samples, suggesting host advantage. Our innovative approach leverages raw read patterns rather than references or MAGs to include all sequencing reads in analysis, and thus provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial genome dynamics.

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.

Two chromosome-level genome assemblies of Rhodiola shed new light on genome evolution in rapid radiation and evolution of the biosynthetic pathway of salidroside

Rhodiola L. is a genus that has undergone rapid radiation in the mid-Miocene and may represent a typic case of adaptive radiation. Many species of Rhodiola have also been widely used as an important adaptogen in traditional medicines for centuries. However, a lack of high-quality chromosome-level genomes hinders in-depth study of its evolution and biosynthetic pathway of secondary metabolites. Here, we assembled two chromosome-level genomes for two Rhodiola species with different chromosome number and sexual system. The assembled genome size of R. chrysanthemifolia (2n = 14; hermaphrodite) and R. kirilowii (2n = 22; dioecious) were of 402.67 and 653.62 Mb, respectively, with approximately 57.60% and 69.22% of transposable elements (TEs). The size difference between the two genomes was mostly due to proliferation of long terminal repeat-retrotransposons (LTR-RTs) in the R. kirilowii genome. Comparative genomic analysis revealed possible gene families responsible for high-altitude adaptation of Rhodiola, including a homolog of plant cysteine oxidase 2 gene of Arabidopsis thaliana (AtPCO2), which is part of the core molecular reaction to hypoxia and contributes to the stability of Group VII ethylene response factors (ERF-VII). We found extensive chromosome fusion/fission events and structural variations between the two genomes, which might have facilitated the initial rapid radiation of Rhodiola. We also identified candidate genes in the biosynthetic pathway of salidroside. Overall, our results provide important insights into genome evolution in plant rapid radiations, and possible roles of chromosome fusion/fission and structure variation played in rapid speciation.

Chromosome-level genome assembly of the kiang (Equus kiang) illuminates genomic basis for its high-altitude adaptation

The kiang (Equus kiang) can only be observed in the Qinghai-Tibet Plateau (QTP). The kiang displayed excellent athletic performance in the high-altitude environment, which attracted wide interest in the investigation of the potential adaptive mechanisms to the extreme environment. Here, we assembled a chromosome-level genome of the kiang based on Hi-C sequencing technology. A total of 324.14 Gb clean data were generated, and the chromosome-level genome with 26 chromosomes (25 + X) and scaffold N50 of 101.77 Mb was obtained for the kiang. The genomic synteny analysis revealed large-scale chromosomal rearrangement during the evolution process of Equus species. Phylogenetic and divergence analyses revealed that the kiang was the sister branch to the ass and diverged from a common ancestor at approximately 13.5 Mya. The expanded gene families were mainly related to the hypoxia response, metabolism, and immunity. The kiang suffered a significant loss of olfaction-related genes, which might indicate decreased olfactory sensibility. Positively selected genes (PSGs) detected in the kiang were mainly associated with hypoxia response. Especially, there were two species-specific missense amino acid mutations in the PSG STAT3 annotated in the hypoxia-inducible factor 1 signal pathway, which may play an important role in the high-altitude adaptation of the kiang. Moreover, structure variations in the kiang genome were also identified, which possibly contributed to the high-altitude adaptation of the kiang. Comparative analysis revealed a lot of species-specific insertions and deletions in the kiang genome, such as PIK3CB and AKT with 3258 and 189 bp insertions in the intron region, respectively, possibly affecting the expression and regulation of hypoxia-related downstream pathways. This study provided valuable genomic resources, and our findings help a better understanding of the underlying adaptive strategies to the high-altitude environment in the kiang.

Influence of genomic variations on glanders serodiagnostic antigens using integrative genomic and transcriptomic approaches

Glanders is a highly contagious and life-threatening zoonotic disease caused by Burkholderia mallei (B. mallei). Without an effective vaccine or treatment, early diagnosis has been regarded as the most effective method to prevent glanders transmission. Currently, the diagnosis of glanders is heavily reliant on serological tests. However, given that markedly different host immune responses can be elicited by genetically different strains of the same bacterial species, infection by B. mallei, whose genome is unstable and plastic, may result in various immune responses. This variability can make the serodiagnosis of glanders challenging. Therefore, there is a need for a comprehensive understanding and assessment of how B. mallei genomic variations impact the appropriateness of specific target antigens for glanders serodiagnosis. In this study, we investigated how genomic variations in the B. mallei genome affect gene content (gene presence/absence) and expression, with a special focus on antigens used or potentially used in serodiagnosis. In all the genome sequences of B. mallei isolates available in NCBI's RefSeq database (accessed in July 2023) and in-house sequenced samples, extensive small and large variations were observed when compared to the type strain ATCC 23344. Further pan-genome analysis of those assemblies revealed variations of gene content among all available genomes of B. mallei. Specifically, differences in gene content ranging from 31 to 715 genes with an average of 334 gene presence-absence variations were found in strains with complete or chromosome-level genome assemblies, using the ATCC 23344 strain as a reference. The affected genes included some encoded proteins used as serodiagnostic antigens, which were lost due mainly to structural variations. Additionally, a transcriptomic analysis was performed using the type strain ATCC 23344 and strain Zagreb which has been widely utilized to produce glanders antigens. In total, 388 significant differentially expressed genes were identified between these two strains, including genes related to bacterial pathogenesis and virulence, some of which were associated with genomic variations, particularly structural variations. To our knowledge, this is the first comprehensive study to uncover the impacts of genetic variations of B. mallei on its gene content and expression. These differences would have significant impacts on host innate and adaptive immunity, including antibody production, during infection. This study provides novel insights into B. mallei genetic variants, knowledge which will help to improve glanders serodiagnosis.

Gene mining and genomics-assisted breeding empowered by the pangenome of tea plant Camellia sinensis

Tea is one of the world's oldest crops and is cultivated to produce beverages with various flavours. Despite advances in sequencing technologies, the genetic mechanisms underlying key agronomic traits of tea remain unclear. In this study, we present a high-quality pangenome of 22 elite cultivars, representing broad genetic diversity in the species. Our analysis reveals that a recent long terminal repeat burst contributed nearly 20% of gene copies, introducing functional genetic variants that affect phenotypes such as leaf colour. Our graphical pangenome improves the efficiency of genome-wide association studies and allows the identification of key genes controlling bud flush timing. We also identified strong correlations between allelic variants and flavour-related chemistries. These findings deepen our understanding of the genetic basis of tea quality and provide valuable genomic resources to facilitate its genomics-assisted breeding.

Multiomics analysis reveals the genetic and metabolic characteristics associated with the low prevalence of dental caries

Background

Despite poor oral hygiene, the Baiku Yao (BKY) ethnic group in China presents a low prevalence of dental caries, which may be related to genetic susceptibility. Due to strict intra-ethnic marriage rule, this ethnic has an advantage in studying the interaction between genetic factors and other regulatory factors related to dental caries.

Methods

Peripheral blood from a caries-free adult male was used for whole genome sequencing, and the BKY assembled genome was compared to the Han Chinese genome. Oral saliva samples were collected from 51 subjects for metabolomic and metagenomic analysis. Multiomics data were integrated for combined analysis using bioinformatics approaches.

Results

Comparative genomic analysis revealed the presence of structural variations in several genes associated with dental caries. Metabolomic and metagenomic sequencing demonstrated the caries-free group had significantly higher concentration of antimicrobials and higher abundance of core oral health-related microbiota. The functional analysis indicated that cationic antimicrobial peptide resistance and the lipopolysaccharide biosynthesis pathway were enriched in the caries-free group.

Conclusions

Our study provided new insights into the specific regulatory mechanisms that contribute to the low prevalence of dental caries in the specific population and may provide new evidence for the genetic diagnosis and control of dental caries.

Genome instability footprint under rapamycin and hydroxyurea treatments

The mutational processes dictating the accumulation of mutations in genomes are shaped by genetic background, environment and their interactions. Accurate quantification of mutation rates and spectra under drugs has important implications in disease treatment. Here, we used whole-genome sequencing and time-resolved growth phenotyping of yeast mutation accumulation lines to give a detailed view of the mutagenic effects of rapamycin and hydroxyurea on the genome and cell growth. Mutation rates depended on the genetic backgrounds but were only marginally affected by rapamycin. As a remarkable exception, rapamycin treatment was associated with frequent chromosome XII amplifications, which compensated for rapamycin induced rDNA repeat contraction on this chromosome and served to maintain rDNA content homeostasis and fitness. In hydroxyurea, a wide range of mutation rates were elevated regardless of the genetic backgrounds, with a particularly high occurrence of aneuploidy that associated with dramatic fitness loss. Hydroxyurea also induced a high T-to-G and low C-to-A transversion rate that reversed the common G/C-to-A/T bias in yeast and gave rise to a broad range of structural variants, including mtDNA deletions. The hydroxyurea mutation footprint was consistent with the activation of error-prone DNA polymerase activities and non-homologues end joining repair pathways. Taken together, our study provides an in-depth view of mutation rates and signatures in rapamycin and hydroxyurea and their impact on cell fitness, which brings insights for assessing their chronic effects on genome integrity.

MetaSVs: A pipeline combining long and short reads for analysis and visualization of structural variants in metagenomes

Structural variants (SVs, including large-scale insertions, deletions, inversions, and translocations) significantly impact the functions of genes in the microbial genome, and SVs in the microbiome are associated with diverse biological processes and human diseases. With the advancements in sequencing and bioinformatics technologies, increasingly, sequencing data and analysis tools are already being extensively utilized for microbiome SV analyses, leading to a higher demand for more dedicated SV analysis workflows. Moreover, due to the unique detection biases of various sequencing technologies, including short-read sequencing (such as Illumina platforms) and long-read sequencing (e.g., Oxford Nanopore and PacBio), SV discovery based on multiple platforms is necessary to comprehensively identify the wide variety of SVs. Here, we establish an integrated pipeline MetaSVs combining Nanopore long reads and Illumina short reads to analyze SVs in the microbial genomes from gut microbiome and further identify differential SVs that can be reflective of metabolic differences. Our pipeline provides researchers easy access to SVs and relevant metabolites in the microbial genomes without the requirement of specific technical expertise, which is particularly useful to researchers interested in metagenomic SVs but lacking sophisticated bioinformatic knowledge.

Comprehensive regulatory networks for tomato organ development based on the genome and RNAome of MicroTom tomato

MicroTom has a short growth cycle and high transformation efficiency, and is a prospective model plant for studying organ development, metabolism, and plant-microbe interactions. Here, with a newly assembled reference genome for this tomato cultivar and abundant RNA-seq data derived from tissues of different organs/developmental stages/treatments, we constructed multiple gene co-expression networks, which will provide valuable clues for the identification of important genes involved in diverse regulatory pathways during plant growth, e.g. arbuscular mycorrhizal symbiosis and fruit development. Additionally, non-coding RNAs, including miRNAs, lncRNAs, and circRNAs were also identified, together with their potential targets. Interacting networks between different types of non-coding RNAs (miRNA-lncRNA), and non-coding RNAs and genes (miRNA-mRNA and lncRNA-mRNA) were constructed as well. Our results and data will provide valuable information for the study of organ differentiation and development of this important fruit. Lastly, we established a database (http://eplant.njau.edu.cn/microTomBase/) with genomic and transcriptomic data, as well as details of gene co-expression and interacting networks on MicroTom, and this database should be of great value to those who want to adopt MicroTom as a model plant for research.

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.

Diallel panel reveals a significant impact of low-frequency genetic variants on gene expression variation in yeast

Unraveling the genetic sources of gene expression variation is essential to better understand the origins of phenotypic diversity in natural populations. Genome-wide association studies identified thousands of variants involved in gene expression variation, however, variants detected only explain part of the heritability. In fact, variants such as low-frequency and structural variants (SVs) are poorly captured in association studies. To assess the impact of these variants on gene expression variation, we explored a half-diallel panel composed of 323 hybrids originated from pairwise crosses of 26 natural Saccharomyces cerevisiae isolates. Using short- and long-read sequencing strategies, we established an exhaustive catalog of single nucleotide polymorphisms (SNPs) and SVs for this panel. Combining this dataset with the transcriptomes of all hybrids, we comprehensively mapped SNPs and SVs associated with gene expression variation. While SVs impact gene expression variation, SNPs exhibit a higher effect size with an overrepresentation of low-frequency variants compared to common ones. These results reinforce the importance of dissecting the heritability of complex traits with a comprehensive catalog of genetic variants at the population level.

Natural genetic variation in the pheromone production of C. elegans

From bacterial quorum sensing to human language, communication is essential for social interactions. Nematodes produce and sense pheromones to communicate among individuals and respond to environmental changes. These signals are encoded by different types and mixtures of ascarosides, whose modular structures further enhance the diversity of this nematode pheromone language. Interspecific and intraspecific differences in this ascaroside pheromone language have been described previously, but the genetic basis and molecular mechanisms underlying the variation remain largely unknown. Here, we analyzed natural variation in the production of 44 ascarosides across 95 wild Caenorhabditis elegans strains using high-performance liquid chromatography coupled to high-resolution mass spectrometry. We discovered wild strains defective in the production of specific subsets of ascarosides (e.g., the aggregation pheromone icas#9) or short- and medium-chain ascarosides, as well as inversely correlated patterns between the production of two major classes of ascarosides. We investigated genetic variants that are significantly associated with the natural differences in the composition of the pheromone bouquet, including rare genetic variants in key enzymes participating in ascaroside biosynthesis, such as the peroxisomal 3-ketoacyl-CoA thiolase, daf-22, and the carboxylesterase cest-3. Genome-wide association mappings revealed genomic loci harboring common variants that affect ascaroside profiles. Our study yields a valuable dataset for investigating the genetic mechanisms underlying the evolution of chemical communication.

Highly complete long-read genomes reveal pangenomic variation underlying yeast phenotypic diversity

Understanding the genetic causes of trait variation is a primary goal of genetic research. One way that individuals can vary genetically is through variable pangenomic genes: genes that are only present in some individuals in a population. The presence or absence of entire genes could have large effects on trait variation. However, variable pangenomic genes can be missed in standard genotyping workflows, owing to reliance on aligning short-read sequencing to reference genomes. A popular method for studying the genetic basis of trait variation is linkage mapping, which identifies quantitative trait loci (QTLs), regions of the genome that harbor causative genetic variants. Large-scale linkage mapping in the budding yeast Saccharomyces cerevisiae has found thousands of QTLs affecting myriad yeast phenotypes. To enable the resolution of QTLs caused by variable pangenomic genes, we used long-read sequencing to generate highly complete de novo genome assemblies of 16 diverse yeast isolates. With these assemblies, we resolved QTLs for growth on maltose, sucrose, raffinose, and oxidative stress to specific genes that are absent from the reference genome but present in the broader yeast population at appreciable frequency. Copies of genes also duplicate onto chromosomes where they are absent in the reference genome, and we found that these copies generate additional QTLs whose resolution requires pangenome characterization. Our findings show the need for highly complete genome assemblies to identify the genetic basis of trait variation.

High-quality genome assembly and comparative genomic profiling of yellowhorn (Xanthoceras sorbifolia) revealed environmental adaptation footprints and seed oil contents variations

Yellowhorn (Xanthoceras sorbifolia) is a species of deciduous tree that is native to Northern and Central China, including Loess Plateau. The yellowhorn tree is a hardy plant, tolerating a wide range of growing conditions, and is often grown for ornamental purposes in parks, gardens, and other landscaped areas. The seeds of yellowhorn are edible and contain rich oil and fatty acid contents, making it an ideal plant for oil production. However, the mechanism of its ability to adapt to extreme environments and the genetic basis of oil synthesis remains to be elucidated. In this study, we reported a high-quality and near gap-less yellowhorn genome assembly, containing the highest genome continuity with a contig N50 of 32.5 Mb. Comparative genomics analysis showed that 1,237 and 231 gene families under expansion and the yellowhorn-specific gene family NB-ARC were enriched in photosynthesis and root cap development, which may contribute to the environmental adaption and abiotic stress resistance of yellowhorn. A 3-ketoacyl-CoA thiolase (KAT) gene (Xso_LG02_00600) was identified under positive selection, which may be associated with variations of seed oil content among different yellowhorn cultivars. This study provided insights into environmental adaptation and seed oil content variations of yellowhorn to accelerate its genetic improvement.

The Chlamydomonas Genome Project, version 6: Reference assemblies for mating-type plus and minus strains reveal extensive structural mutation in the laboratory

Five versions of the Chlamydomonas reinhardtii reference genome have been produced over the last two decades. Here we present version 6, bringing significant advances in assembly quality and structural annotations. PacBio-based chromosome-level assemblies for two laboratory strains, CC-503 and CC-4532, provide resources for the plus and minus mating-type alleles. We corrected major misassemblies in previous versions and validated our assemblies via linkage analyses. Contiguity increased over ten-fold and >80% of filled gaps are within genes. We used Iso-Seq and deep RNA-seq datasets to improve structural annotations, and updated gene symbols and textual annotation of functionally characterized genes via extensive manual curation. We discovered that the cell wall-less classical reference strain CC-503 exhibits genomic instability potentially caused by deletion of the helicase RECQ3, with major structural mutations identified that affect >100 genes. We therefore present the CC-4532 assembly as the primary reference, although this strain also carries unique structural mutations and is experiencing rapid proliferation of a Gypsy retrotransposon. We expect all laboratory strains to harbor gene-disrupting mutations, which should be considered when interpreting and comparing experimental results. Collectively, the resources presented here herald a new era of Chlamydomonas genomics and will provide the foundation for continued research in this important reference organism.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

Rates and spectra of de novo structural mutations in Chlamydomonas reinhardtii

Genetic variation originates from several types of spontaneous mutation, including single-nucleotide substitutions, short insertions and deletions (indels), and larger structural changes. Structural mutations (SMs) drive genome evolution and are thought to play major roles in evolutionary adaptation, speciation, and genetic disease, including cancers. Sequencing of mutation accumulation (MA) lines has provided estimates of rates and spectra of single-nucleotide and indel mutations in many species, yet the rate of new SMs is largely unknown. Here, we use long-read sequencing to determine the full mutation spectrum in MA lines derived from two strains (CC-1952 and CC-2931) of the green alga Chlamydomonas reinhardtii The SM rate is highly variable between strains and between MA lines, and SMs represent a substantial proportion of all mutations in both strains (CC-1952 6%; CC-2931 12%). The SM spectra differ considerably between the two strains, with almost all inversions and translocations occurring in CC-2931 MA lines. This variation is associated with heterogeneity in the number and type of active transposable elements (TEs), which comprise major proportions of SMs in both strains (CC-1952 22%; CC-2931 38%). In CC-2931, a Crypton and a previously undescribed type of DNA element have caused 71% of chromosomal rearrangements, whereas in CC-1952, a Dualen LINE is associated with 87% of duplications. Other SMs, notably large duplications in CC-2931, are likely products of various double-strand break repair pathways. Our results show that diverse types of SMs occur at substantial rates, and support prominent roles for SMs and TEs in evolution.

Different structural variant prediction tools yield considerably different results in Caenorhabditis elegans

The accurate characterization of structural variation is crucial for our understanding of how large chromosomal alterations affect phenotypic differences and contribute to genome evolution. Whole-genome sequencing is a popular approach for identifying structural variants, but the accuracy of popular tools remains unclear due to the limitations of existing benchmarks. Moreover, the performance of these tools for predicting variants in non-human genomes is less certain, as most tools were developed and benchmarked using data from the human genome. To evaluate the use of long-read data for the validation of short-read structural variant calls, the agreement between predictions from a short-read ensemble learning method and long-read tools were compared using real and simulated data from Caenorhabditis elegans. The results obtained from simulated data indicate that the best performing tool is contingent on the type and size of the variant, as well as the sequencing depth of coverage. These results also highlight the need for reference datasets generated from real data that can be used as 'ground truth' in benchmarks.

Short- and long-read metagenomics expand individualized structural variations in gut microbiomes

In-depth profiling of genetic variations in the gut microbiome is highly desired for understanding its functionality and impacts on host health and disease. Here, by harnessing the long read advantage provided by Oxford Nanopore Technology (ONT), we characterize fine-scale genetic variations of structural variations (SVs) in hundreds of gut microbiomes from healthy humans. ONT long reads dramatically improve the quality of metagenomic assemblies, enable reliable detection of a large, expanded set of structural variation types (notably including large insertions and inversions). We find SVs are highly distinct between individuals and stable within an individual, representing gut microbiome fingerprints that shape strain-level differentiations in function within species, complicating the associations to metabolites and host phenotypes such as blood glucose. In summary, our study strongly emphasizes that incorporating ONT reads into metagenomic analyses expands the detection scope of genetic variations, enables profiling strain-level variations in gut microbiome, and their intricate correlations with metabolome.

The chromosome-level genome assembly of the Japanese yellowtail jack Seriola aureovittata provides insights into genome evolution and efficient oxygen transport

Fishes of the genus Seriola are widely farmed and highly valued in global aquaculture production. To further understand their economically important traits and help improve aquaculture product quality and sustainability, we performed a chromosome-level genome construction for Seriola aureovittata. Combining two technologies, PacBio and BGISEQ-500, we assembled 649.86 Mb S. aureovittata genome sequences with a contig N50 of 22.21 Mb, and 98% of BUSCO genes were detected in total. The initial assembly was then further scaffolded into 24 pseudochromosomes using Hi-C data, indicating the high quality of the genome. Genome evolution analysis showed that many genes related to fatty acid metabolism and oxygen binding, or transport were expanded, which provided insights into the metabolic characteristics of fatty acids and efficient oxygen transport. Based on the genome data, we confirmed the evolutionary relationship of S. aureovittata, S. dorsalis and S. lalandi and identified chr12 as the putative sex chromosome of S. aureovittata. Our chromosome-level genome assembly provides a genetic foundation for the phylogenetic and taxonomic investigation of different Seriola species. Moreover, the genome will provide an important genomic resource for further biological and aquaculture studies of S. aureovittata.

plotsr: visualizing structural similarities and rearrangements between multiple genomes

SUMMARY

Third-generation genome sequencing technologies have led to a sharp increase in the number of high-quality genome assemblies. This allows the comparison of multiple assembled genomes of individual species and demands new tools for visualizing their structural properties. Here, we present plotsr, an efficient tool to visualize structural similarities and rearrangements between genomes. It can be used to compare genomes on chromosome level or to zoom in on any selected region. In addition, plotsr can augment the visualization with regional identifiers (e.g. genes or genomic markers) or histogram tracks for continuous features (e.g. GC content or polymorphism density).

AVAILABILITY AND IMPLEMENTATION

plotsr is implemented as a python package and uses the standard matplotlib library for plotting. It is freely available under the MIT license at GitHub (https://github.com/schneebergerlab/plotsr) and bioconda (https://anaconda.org/bioconda/plotsr).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Saccharomyces eubayanus haploid resource for research studies

Since its identification, Saccharomyces eubayanus has been recognized as the missing parent of the lager hybrid, S. pastorianus. This wild yeast has never been isolated from fermentation environments, thus representing an interesting candidate for evolutionary, ecological and genetic studies. However, it is imperative to develop additional molecular genetics tools to ease manipulation and thus facilitate future studies. With this in mind, we generated a collection of stable haploid strains representative of three main lineages described in S. eubayanus (PB-1, PB-2 and PB-3), by deleting the HO gene using CRISPR-Cas9 and tetrad micromanipulation. Phenotypic characterization under different conditions demonstrated that the haploid derivates were extremely similar to their parental strains. Genomic analysis in three strains highlighted a likely low frequency of off-targets, and sequencing of a single tetrad evidenced no structural variants in any of the haploid spores. Finally, we demonstrate the utilization of the haploid set by challenging the strains under mass-mating conditions. In this way, we found that S. eubayanus under liquid conditions has a preference to remain in a haploid state, unlike S. cerevisiae that mates rapidly. This haploid resource is a novel set of strains for future yeast molecular genetics studies.

Genome assembly of the Chinese maize elite inbred line RP125 and its EMS mutant collection provide new resources for maize genetics research and crop improvement

Maize is an important crop worldwide, as well as a valuable model with vast genetic diversity. Accurate genome and annotation information for a wide range of inbred lines would provide valuable resources for crop improvement and pan-genome characterization. In this study, we generated a high-quality de novo genome assembly (contig N50 of 15.43 Mb) of the Chinese elite inbred line RP125 using Nanopore long-read sequencing and Hi-C scaffolding, which yield highly contiguous, chromosome-length scaffolds. Global comparison of the RP125 genome with those of B73, W22, and Mo17 revealed a large number of structural variations. To create new germplasm for maize research and crop improvement, we carried out an EMS mutagenesis screen on RP125. In total, we obtained 5818 independent M2 families, with 946 mutants showing heritable phenotypes. Taking advantage of the high-quality RP125 genome, we successfully cloned 10 mutants from the EMS library, including the novel kernel mutant qk1 (quekou: "missing a small part" in Chinese), which exhibited partial loss of endosperm and a starch accumulation defect. QK1 encodes a predicted metal tolerance protein, which is specifically required for Fe transport. Increased accumulation of Fe and reactive oxygen species as well as ferroptosis-like cell death were detected in qk1 endosperm. Our study provides the community with a high-quality genome sequence and a large collection of mutant germplasm.

Genomic structural variants constrain and facilitate adaptation in natural populations of Theobroma cacao, the chocolate tree

Genomic structural variants (SVs) can play important roles in adaptation and speciation. Yet the overall fitness effects of SVs are poorly understood, partly because accurate population-level identification of SVs requires multiple high-quality genome assemblies. Here, we use 31 chromosome-scale, haplotype-resolved genome assemblies of Theobroma cacao-an outcrossing, long-lived tree species that is the source of chocolate-to investigate the fitness consequences of SVs in natural populations. Among the 31 accessions, we find over 160,000 SVs, which together cover eight times more of the genome than single-nucleotide polymorphisms and short indels (125 versus 15 Mb). Our results indicate that a vast majority of these SVs are deleterious: they segregate at low frequencies and are depleted from functional regions of the genome. We show that SVs influence gene expression, which likely impairs gene function and contributes to the detrimental effects of SVs. We also provide empirical support for a theoretical prediction that SVs, particularly inversions, increase genetic load through the accumulation of deleterious nucleotide variants as a result of suppressed recombination. Despite the overall detrimental effects, we identify individual SVs bearing signatures of local adaptation, several of which are associated with genes differentially expressed between populations. Genes involved in pathogen resistance are strongly enriched among these candidates, highlighting the contribution of SVs to this important local adaptation trait. Beyond revealing empirical evidence for the evolutionary importance of SVs, these 31 de novo assemblies provide a valuable resource for genetic and breeding studies in T cacao.

Rapid selection response to ethanol in Saccharomyces eubayanus emulates the domestication process under brewing conditions

Although the typical genomic and phenotypic changes that characterize the evolution of organisms under the human domestication syndrome represent textbook examples of rapid evolution, the molecular processes that underpin such changes are still poorly understood. Domesticated yeasts for brewing, where short generation times and large phenotypic and genomic plasticity were attained in a few generations under selection, are prime examples. To experimentally emulate the lager yeast domestication process, we created a genetically complex (panmictic) artificial population of multiple Saccharomyces eubayanus genotypes, one of the parents of lager yeast. Then, we imposed a constant selection regime under a high ethanol concentration in 10 replicated populations during 260 generations (6 months) and compared them with propagated controls exposed solely to glucose. Propagated populations exhibited a selection differential of 60% in growth rate in ethanol, mostly explained by the proliferation of a single lineage (CL248.1) that competitively displaced all other clones. Interestingly, the outcome does not require the entire time‐course of adaptation, as four lineages monopolized the culture at generation 120. Sequencing demonstrated that de novo genetic variants were produced in all propagated lines, including SNPs, aneuploidies, INDELs and translocations. In addition, the different propagated populations showed correlated responses resembling the domestication syndrome: genomic rearrangements, faster fermentation rates, lower production of phenolic off‐flavours and lower volatile compound complexity. Expression profiling in beer wort revealed altered expression levels of genes related to methionine metabolism, flocculation, stress tolerance and diauxic shift, likely contributing to higher ethanol and fermentation stress tolerance in the evolved populations. Our study shows that experimental evolution can rebuild the brewing domestication process in ‘fast motion’ in wild yeast, and also provides a powerful tool for studying the genetics of the adaptation process in complex populations.

Genomic Evidence of an Ancient East Asian Divergence Event in Wild Saccharomyces cerevisiae

Comparative genome analyses have suggested East Asia to be the cradle of the domesticated microbe Brewer's yeast (Saccharomyces cerevisiae), used in the food and biotechnology industry worldwide. Here, we provide seven new, high-quality long-read genomes of nondomesticated yeast strains isolated from primeval forests and other natural environments in China and Taiwan. In a comprehensive analysis of our new genome assemblies, along with other long-read Saccharomycetes genomes available, we show that the newly sequenced East Asian strains are among the closest living relatives of the ancestors of the global diversity of Brewer's yeast, confirming predictions made from short-read genomic data. Three of these strains (termed the East Asian Clade IX Complex here) share a recent ancestry and evolutionary history suggesting an early divergence from other S. cerevisiae strains before the larger radiation of the species, and prior to its domestication. Our genomic analyses reveal that the wild East Asian strains contain elevated levels of structural variations. The new genomic resources provided here contribute to our understanding of the natural diversity of S. cerevisiae, expand the intraspecific genetic variation found in this heavily domesticated microbe, and provide a foundation for understanding its origin and global colonization history.