ABSTRACT
The life cycle of genome builds spans interlocking pillars of assembly, annotation, and comparative genomics to drive biological insights. While tools exist to address each pillar separately, there is a growing need for tools to integrate different pillars of a genome project holistically. For example, comparative approaches can provide quality control of assembly or annotation; genome assembly, in turn, can help to identify artifacts that may complicate the interpretation of genome comparisons. The JCVI library is a versatile Python-based library that offers a suite of tools that excel across these pillars. Featuring a modular design, the JCVI library provides high-level utilities for tasks such as format parsing, graphics generation, and manipulation of genome assemblies and annotations. Supporting genomics algorithms like MCscan and ALLMAPS are widely employed in building genome releases, producing publication-ready figures for quality assessment and evolutionary inference. Developed and maintained collaboratively, the JCVI library emphasizes quality and reusability.
ABSTRACT
Scaffolding is crucial for constructing most chromosome-level genomes. The high-throughput chromatin conformation capture (Hi-C) technology has become the primary scaffolding strategy due to its convenience and cost-effectiveness. As sequencing technologies and assembly algorithms advance, constructing haplotype-resolved genomes is increasingly preferred because haplotypes can provide additional genetic information on allelic and non-allelic variations. ALLHiC is a widely used allele-aware scaffolding tool designed for this purpose. However, its dependence on chromosome-level reference genomes and a higher chromosome misassignment rate still impede the unravelling of haplotype-resolved genomes. Here we present HapHiC, a reference-independent allele-aware scaffolding tool with superior performance on chromosome assignment as well as contig ordering and orientation. In addition, we provide new insights into the challenges in allele-aware scaffolding by conducting comprehensive analyses on various adverse factors. Finally, with the help of HapHiC, we constructed the haplotype-resolved allotriploid genome for Miscanthus × giganteus, an important lignocellulosic bioenergy crop.
Subject(s)
Chromosomes, Plant , Genome, Plant , Haplotypes , Chromosomes, Plant/genetics , Chromatin/genetics , Poaceae/genetics , High-Throughput Nucleotide Sequencing/methods , AllelesABSTRACT
Vernicia montana is a dioecious plant widely cultivated for high-quality tung oil production and ornamental purposes in the Euphorbiaceae family. The lack of genomic information has severely hindered molecular breeding for genetic improvement and early sex identification in V. montana. Here, we present a chromosome-level reference genome of a male V. montana with a total size of 1.29 Gb and a contig N50 of 3.69 Mb. Genome analysis revealed that different repeat lineages drove the expansion of genome size. The model of chromosome evolution in the Euphorbiaceae family suggests that polyploidization-induced genomic structural variation reshaped the chromosome structure, giving rise to the diverse modern chromosomes. Based on whole-genome resequencing data and analyses of selective sweep and genetic diversity, several genes associated with stress resistance and flavonoid synthesis such as CYP450 genes and members of the LRR-RLK family, were identified and presumed to have been selected during the evolutionary process. Genome-wide association studies were conducted and a putative sex-linked insertion and deletion (InDel) (Chr 2: 102 799 917-102 799 933 bp) was identified and developed as a polymorphic molecular marker capable of effectively detecting the gender of V. montana. This InDel is located in the second intron of VmBASS4, suggesting a possible role of VmBASS4 in sex determination in V. montana. This study sheds light on the genome evolution and sex identification of V. montana, which will facilitate research on the development of agronomically important traits and genomics-assisted breeding.
ABSTRACT
Camellia oil is valuable as an edible oil and serves as a base material for a range of high-value products. Camellia plants of significant economic importance, such as Camellia sinensis and Camellia oleifera, have been classified into sect. Thea and sect. Oleifera, respectively. Fatty acid desaturases play a crucial role in catalyzing the formation of double bonds at specific positions of fatty acid chains, leading to the production of unsaturated fatty acids and contributing to lipid synthesis. Comparative genomics results have revealed that expanded gene families in oil tea are enriched in functions related to lipid, fatty acid, and seed processes. To explore the function of the FAD gene family, a total of 82 FAD genes were identified in tea and oil tea. Transcriptome data showed the differential expression of the FAD gene family in mature seeds of tea tree and oil tea tree. Furthermore, the structural analysis and clustering of FAD proteins provided insights for the further exploration of the function of the FAD gene family and its role in lipid synthesis. Overall, these findings shed light on the role of the FAD gene family in Camellia plants and their involvement in lipid metabolism, as well as provide a reference for understanding their function in oil synthesis.
ABSTRACT
Sugarcane is a vital crop with significant economic and industrial value. However, the cultivated sugarcane's ultra-complex genome still needs to be resolved due to its high ploidy and extensive recombination between the two subgenomes. Here, we generate a chromosomal-scale, haplotype-resolved genome assembly for a hybrid sugarcane cultivar ZZ1. This assembly contains 10.4 Gb genomic sequences and 68,509 annotated genes with defined alleles in two sub-genomes distributed in 99 original and 15 recombined chromosomes. RNA-seq data analysis shows that sugar accumulation-associated gene families have been primarily expanded from the ZZSO subgenome. However, genes responding to pokkah boeng disease susceptibility have been derived dominantly from the ZZSS subgenome. The region harboring the possible smut resistance genes has expanded significantly. Among them, the expansion of WAK and FLS2 families is proposed to have occurred during the breeding of ZZ1. Our findings provide insights into the complex genome of hybrid sugarcane cultivars and pave the way for future genomics and molecular breeding studies in sugarcane.
Subject(s)
Saccharum , Saccharum/genetics , Plant Breeding , Genomics , Haplotypes/genetics , ChromosomesABSTRACT
Following the milestone success of the Human Genome Project, the 'Encyclopedia of DNA Elements (ENCODE)' initiative was launched in 2003 to unearth information about the numerous functional elements within the genome. This endeavor coincided with the emergence of numerous novel technologies, accompanied by the provision of vast amounts of whole-genome sequences, high-throughput data such as ChIP-Seq and RNA-Seq. Extracting biologically meaningful information from this massive dataset has become a critical aspect of many recent studies, particularly in annotating and predicting the functions of unknown genes. The core idea behind genome annotation is to identify genes and various functional elements within the genome sequence and infer their biological functions. Traditional wet-lab experimental methods still rely on extensive efforts for functional verification. However, early bioinformatics algorithms and software primarily employed shallow learning techniques; thus, the ability to characterize data and features learning was limited. With the widespread adoption of RNA-Seq technology, scientists from the biological community began to harness the potential of machine learning and deep learning approaches for gene structure prediction and functional annotation. In this context, we reviewed both conventional methods and contemporary deep learning frameworks, and highlighted novel perspectives on the challenges arising during annotation underscoring the dynamic nature of this evolving scientific landscape.
Subject(s)
Deep Learning , Humans , Genome , Algorithms , Software , Computational Biology/methods , Molecular Sequence AnnotationABSTRACT
MOTIVATION: Constructing a phylogenetic tree requires calculating the evolutionary distance between samples or species via large-scale resequencing data, a process that is both time-consuming and computationally demanding. Striking the right balance between accuracy and efficiency is a significant challenge. RESULTS: To address this, we introduce a new algorithm, MIKE (MinHash-based k-mer algorithm). This algorithm is designed for the swift calculation of the Jaccard coefficient directly from raw sequencing reads and enables the construction of phylogenetic trees based on the resultant Jaccard coefficient. Simulation results highlight the superior speed of MIKE compared to existing state-of-the-art methods. We used MIKE to reconstruct a phylogenetic tree, incorporating 238 yeast, 303 Zea, 141 Ficus, 67 Oryza, and 43 Saccharum spontaneum samples. MIKE demonstrated accurate performance across varying evolutionary scales, reproductive modes, and ploidy levels, proving itself as a powerful tool for phylogenetic tree construction. AVAILABILITY AND IMPLEMENTATION: MIKE is publicly available on Github at https://github.com/Argonum-Clever2/mike.git.
Subject(s)
Algorithms , Phylogeny , Computer Simulation , Sequence Analysis, DNAABSTRACT
A high-quality reference genome is indispensable for resolving biologically essential traits. Ficus hispida is a dioecious plant. A complete Ficus reference genome will be crucial for understanding their sex evolution and important biological characteristics, such as aerial roots, mutualistic symbiosis with ficus-wasps, and fruiting from old stems. Here, we generated a telomere-to-telomere (T2T) genome for F. hispida using PacBio HiFi and Oxford Nanopore Ultra-long sequencing technologies. The genome contiguity and completeness has shown improvement compared with the previously released genome, with the annotation of six centromeres and 28 telomeres. We have refined our previously reported 2-Mb male-specific region into a 7.2-Mb genomic region containing 51 newly predicted genes and candidate sex-determination genes AG2 and AG3. Many of these genes showed extremely low expression, likely attributed to hypermethylation in the gene body and promoter regions. Gene regulatory networks (GRNs) revealed that AG2 and AG3 are related to the regulation of stamen development in male flowers, while the AG1 gene is responsible for regulating female flowers' defense responses and secondary metabolite processes. Comparative analysis of GRNs showed that the NAC, WRKY, and MYB transcription factor families dominate the female GRN, whereas the MADS and MYB transcription factor families are prevalent in the male GRN.
ABSTRACT
During the oolong tea withering process, abiotic stresses induce significant changes in the content of various flavor substances and jasmonic acid (JA). However, the changes in chromatin accessibility during withering and their potential impact remain poorly understood. By integrating ATAC-seq, RNA-seq, metabolite, and hormone assays, we characterized the withering treatment-induced changes in chromatin accessibility, gene expression levels, important metabolite contents, and JA and JA-ILE contents. Additionally, we analyzed the effects of chromatin accessibility alterations on gene expression changes, content changes of important flavor substances, and JA hyperaccumulation. Our analysis identified a total of 3451 open- and 13 426 close-differentially accessible chromatin regions (DACRs) under withering treatment. Our findings indicate that close-DACRs-mediated down-regulated differentially expressed genes (DEGs) resulted in the reduced accumulation of multiple catechins during withering, whereas open-DACRs-mediated up-regulated DEGs contributed to the increased accumulation of important terpenoids, JA, JA-ILE and short-chain C5/C6 volatiles. We further highlighted important DACRs-mediated DEGs associated with the synthesis of catechins, terpenoids, JA and JA and short-chain C5/C6 volatiles and confirmed the broad effect of close-DACRs on catechin synthesis involving almost all enzymes in the pathway during withering. Importantly, we identified a novel MYB transcription factor (CsMYB83) regulating catechin synthesis and verified the binding of CsMYB83 in the promoter-DACRs regions of key catechin synthesis genes using DAP-seq. Overall, our results not only revealed a landscape of chromatin alters-mediated transcription, flavor substance and hormone changes under oolong tea withering, but also provided target genes for flavor improvement breeding in tea plant.
Subject(s)
Catechin , Cyclopentanes , Isoleucine/analogs & derivatives , Oxylipins , Transcriptome , Catechin/analysis , Catechin/metabolism , Chromatin/genetics , Chromatin/metabolism , Plant Breeding , Tea/chemistry , Tea/metabolism , Hormones/analysis , Hormones/metabolism , Terpenes/metabolism , Plant Leaves/metabolismABSTRACT
Upland rice is a distinctive drought-aerobic ecotype of cultivated rice highly resistant to drought stress. However, the genetic and genomic basis for the drought-aerobic adaptation of upland rice remains largely unclear due to the lack of genomic resources. In this study, we identified 25 typical upland rice accessions and assembled a high-quality genome of one of the typical upland rice varieties, IRAT109, comprising 384 Mb with a contig N50 of 19.6 Mb. Phylogenetic analysis revealed upland and lowland rice have distinct ecotype differentiation within the japonica subgroup. Comparative genomic analyses revealed that adaptive differentiation of lowland and upland rice is likely attributable to the natural variation of many genes in promoter regions, formation of specific genes in upland rice, and expansion of gene families. We revealed differentiated gene expression patterns in the leaves and roots of the two ecotypes and found that lignin synthesis mediated by the phenylpropane pathway plays an important role in the adaptive differentiation of upland and lowland rice. We identified 28 selective sweeps that occurred during domestication and validated that the qRT9 gene in selective regions can positively regulate drought resistance in rice. Eighty key genes closely associated with drought resistance were appraised for their appreciable potential in drought resistance breeding. Our study enhances the understanding of the adaptation of upland rice and provides a genome navigation map of drought resistance breeding, which will facilitate the breeding of drought-resistant rice and the "blue revolution" in agriculture.
Subject(s)
Drought Resistance , Oryza , Oryza/metabolism , Phylogeny , Plant Breeding , Droughts , GenomicsABSTRACT
Actinidia ('Mihoutao' in Chinese) includes species with complex ploidy, among which diploid Actinidia chinensis and hexaploid Actinidia deliciosa are economically and nutritionally important fruit crops. Actinidia deliciosa has been proposed to be an autohexaploid (2n = 174) with diploid A. chinensis (2n = 58) as the putative parent. A CCS-based assembly anchored to a high-resolution linkage map provided a chromosome-resolved genome for hexaploid A. deliciosa yielded a 3.91-Gb assembly of 174 pseudochromosomes comprising 29 homologous groups with 6 members each, which contain 39 854 genes with an average of 4.57 alleles per gene. Here we provide evidence that much of the hexaploid genome matches diploid A. chinensis; 95.5% of homologous gene pairs exhibited >90% similarity. However, intragenome and intergenome comparisons of synteny indicate chromosomal changes. Our data, therefore, indicate that if A. deliciosa is an autoploid, chromosomal rearrangement occurred following autohexaploidy. A highly diversified pattern of gene expression and a history of rapid population expansion after polyploidisation likely facilitated the adaptation and niche differentiation of A. deliciosa in nature. The allele-defined hexaploid genome of A. deliciosa provides new genomic resources to accelerate crop improvement and to understand polyploid genome evolution.
Subject(s)
Actinidia , Actinidia/genetics , Chromosome Mapping , Genome, Plant/genetics , Ploidies , Chromosomes , Fruit/geneticsABSTRACT
Oolong tea has gained great popularity in China due to its pleasant floral and fruity aromas. Although numerous studies have investigated the aroma differences across various tea cultivars, the genetic mechanism is unclear. This study performed multiomics analysis of three varieties suitable for oolong tea and three others with different processing suitability. Our analysis revealed that oolong tea varieties contained higher levels of cadinane sesquiterpenoids. PanTFBS was developed to identify variants of transcription factor binding sites (TFBSs). We found that the CsDCS gene had two TFBS variants in the promoter sequence and a single nucleotide polymorphism (SNP) in the coding sequence. Integrating data on genetic variations, gene expression, and protein-binding sites indicated that CsDCS might be a pivotal gene involved in the biosynthesis of cadinane sesquiterpenoids. These findings advance our understanding of the genetic factors involved in the aroma formation of oolong tea and offer insights into the enhancement of tea aroma.
Subject(s)
Camellia sinensis , Sesquiterpenes , Volatile Organic Compounds , Camellia sinensis/genetics , Camellia sinensis/chemistry , Multiomics , Plant Leaves/chemistry , Volatile Organic Compounds/metabolism , Polycyclic Sesquiterpenes/analysis , Polycyclic Sesquiterpenes/metabolism , Sesquiterpenes/metabolism , Tea/chemistryABSTRACT
Peach (Prunus persica) is an economically important fruit crop globally and an excellent material for genomic studies. While considerable progress has been made in unveiling trait-associated genes within cultivars and wild relatives, certain novel genes controlling valuable traits in peach landraces, such as the red-flowering gene, remained unclear. In this study, we sequenced and assembled the diploid genome of the red-flower landrace 'Yingzui' (abbreviated as 'RedY'). Multi-omics profiling of red petals of 'RedY' revealed the intensified red coloration associated with anthocyanins accumulation and concurrent decline in flavonols. This phenomenon is likely attributed to a natural variant of Flavonol Synthase (FLS) harboring a 9-bp exonic insertion. Intriguingly, the homozygous allelic configurations of this FLS variant were only observed in red-flowered peaches. Furthermore, the 9-bp sequence variation tightly associated with pink/red petal color in genome-wide association studies (GWAS) of collected peach germplasm resources. Functional analyses of the FLS variant, purified from procaryotic expression system, demonstrated its diminished enzymatic activity in flavonols biosynthesis, impeccably aligning with the cardinal trait of red flowers. Therefore, the natural FLS variant was proposed as the best candidate gene for red-flowering trait in peach. The pioneering unveiling of the red-flowered peach genome, coupled with the identification of the candidate gene, expanded the knowledge boundaries of the genetic basis of peach traits and provided valuable insights for future peach breeding efforts.
ABSTRACT
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.
Subject(s)
Camellia sinensis , Camellia sinensis/genetics , Genome-Wide Association Study , Plant Breeding , Genomics , TeaABSTRACT
Ixodes scapularis, the black-legged tick, is the principal vector of the Lyme disease spirochete, Borrelia burgdorferi, and is responsible for most of the â¼470,000 estimated Lyme disease cases annually in the USA. Ixodes scapularis can transmit six additional pathogens of human health significance. Because of its medical importance, I. scapularis was the first tick genome to be sequenced and annotated. However, the first assembly, I. scapularis Wikel (IscaW), was highly fragmented because of the technical challenges posed by the long, repetitive genome sequences characteristic of arthropod genomes and the lack of long-read sequencing techniques. Although I. scapularis has emerged as a model for tick research because of the availability of new tools such as embryo injection and CRISPR-Cas9-mediated gene editing yet the lack of chromosome-scale scaffolds has slowed progress in tick biology and the development of tools for their control. Here we combine diverse technologies to produce the I. scapularis Gulia-Nuss (IscGN) genome assembly and gene set. We used DNA from eggs and male and female adult ticks and took advantage of Hi-C, PacBio HiFi sequencing, and Illumina short-read sequencing technologies to produce a chromosome-level assembly. In this work, we present the predicted pseudochromosomes consisting of 13 autosomes and the sex pseudochromosomes: X and Y, and a markedly improved genome annotation compared with the existing assemblies and annotations.
Subject(s)
Borrelia burgdorferi , Ixodes , Lyme Disease , Animals , Male , Female , Humans , Ixodes/genetics , Lyme Disease/genetics , Borrelia burgdorferi/genetics , Genome/genetics , High-Throughput Nucleotide SequencingABSTRACT
Cyclocarya paliurus, an endemic species in the genus Juglandaceae with the character of heterodichogamy, is one of triterpene-rich medicinal plants in China. To uncover the genetic mechanisms behind the special characteristics, we sequenced the genomes of two diploid (protandry, PA-dip and protogyny, PG-dip) and one auto-tetraploid (PA-tetra) C. paliurus genomes. Based on 134.9 (~225x), 75.5 (~125x) and 271.8 Gb (~226x) subreads of PacBio platform sequencing data, we assembled 586.62 Mb (contig N50 = 1.9 Mb), 583.45 Mb (contig N50 = 1.4 Mb), and 2.38 Gb (contig N50 = 430.9 kb) for PA-dip, PG-dip and PA-tetra genome, respectively. Furthermore, 543.53, 553.87, and 2168.65 Mb in PA-dip, PG-dip, and PA-tetra, were respectively anchored to 16, 16, and 64 pseudo-chromosomes using over 65.4 Gb (~109x), 68 Gb (~113x), and 264 (~220x) Hi-C sequencing data. Annotation of PA-dip, PG-dip, and PA-tetra genome assembly identified 34,699, 35,221, and 34,633 protein-coding genes (90,752 gene models) or allele-defined genes, respectively. In addition, 45 accessions from nine locations were re-sequenced, and more than 10 × coverage reads were generated.
Subject(s)
Genome, Plant , Juglandaceae , Chromosomes , Diploidy , Juglandaceae/genetics , Molecular Sequence Annotation , Phylogeny , TetraploidyABSTRACT
In plants, 5mC DNA methylation is an important and conserved epistatic mark involving genomic stability, gene transcriptional regulation, developmental regulation, abiotic stress response, metabolite synthesis, etc. However, the roles of 5mC DNA methylation modification (5mC methylation) in tea plant growth and development (in pre-harvest processing) and flavor substance synthesis in pre- and post-harvest processing are unknown. We therefore conducted a comprehensive methylation analysis of four key pre-harvest tissues (root, leaf, flower, and fruit) and two processed leaves during oolong tea post-harvest processing. We found that differential 5mC methylation among four key tissues is closely related to tissue functional differentiation and that genes expressed tissue-specifically, responsible for tissue-specific functions, maintain relatively low 5mC methylation levels relative to non-tissue-specifically expressed genes. Importantly, hypomethylation modifications of CsAlaDC and TS/GS genes in roots provided the molecular basis for the dominant synthesis of theanine in roots. In addition, integration of 5mC DNA methylationomics, metabolomics, and transcriptomics of post-harvest leaves revealed that content changes in flavor metabolites during oolong tea processing were closely associated with transcription level changes in corresponding metabolite synthesis genes, and changes in transcript levels of these important synthesis genes were strictly regulated by 5mC methylation. We further report that some key genes during processing are regulated by 5mC methylation, which can effectively explain the content changes of important aroma metabolites, including α-farnesene, nerolidol, lipids, and taste substances such as catechins. Our results not only highlight the key roles of 5mC methylation in important flavor substance synthesis in pre- and post-harvest processing, but also provide epimutation-related gene targets for future improvement of tea quality or breeding of whole-tissue high-theanine varieties.