High-quality Reference Genome Research Articles

Multitrait engineering of Hassawi red rice for sustainable cultivation

Sustainable agriculture requires locally adapted varieties that produce nutritious food with limited agricultural inputs. Genome engineering represents a viable approach to develop cultivars that fulfill these criteria. For example, the red Hassawi rice, a native landrace of Saudi Arabia, tolerates local drought and high-salinity conditions and produces grain with diverse health-promoting phytochemicals. However, Hassawi has a long growth cycle, high cultivation costs, low productivity, and susceptibility to lodging. Here, to improve these undesirable traits via genome editing, we established efficient regeneration and Agrobacterium-mediated transformation protocols for Hassawi. In addition, we generated the first high-quality reference genome and targeted the key flowering repressor gene, Hd4, thus shortening the plant's lifecycle and height. Using CRISPR/Cas9 multiplexing, we simultaneously disrupted negative regulators of flowering time (Hd2, Hd4, and Hd5), grain size (GS3), grain number (GN1a), and plant height (Sd1). The resulting homozygous mutant lines flowered extremely early (∼56 days) and had shorter stems (approximately 107 cm), longer grains (by 5.1%), and more grains per plant (by 50.2%), thereby enhancing overall productivity. Furthermore, the awns of grains were 86.4% shorter compared to unedited plants. Moreover, the modified rice grain displayed improved nutritional attributes. As a result, the modified Hassawi rice combines several desirable traits that can incentivize large-scale cultivation and reduce malnutrition.

CpG Island Definition and Methylation Mapping of the T2T-YAO Genome.

Precisely defining and mapping all cytosine (C) positions and their clusters, known as CpG islands (CGIs), as well as their methylation status, are pivotal for genome-wide epigenetic studies, especially when population-centric reference genomes are ready for timely application. Here, we first align the two high-quality reference genomes, T2T-YAO and T2T-CHM13, from different ethnic backgrounds in a base-by-base fashion and compute their genome-wide density-defined and position-defined CGIs. Second, by mapping some representative genome-wide methylation data from selected organs onto the two genomes, we find that there are about 4.7%-5.8% sequence divergency of variable categories depending on quality cutoffs. Genes among the divergent sequences are mostly associated with neurological functions. Moreover, CGIs associated with the divergent sequences are significantly different with respect to CpG density and observed CpG/expected CpG (O/E) ratio between the two genomes. Finally, we find that the T2T-YAO genome not only has a greater CpG coverage than that of the T2T-CHM13 genome when whole-genome bisulfite sequencing (WGBS) data from the European and American populations are mapped to each reference, but also shows more hyper-methylated CpG sites as compared to the T2T-CHM13 genome. Our study suggests that future genome-wide epigenetic studies of the Chinese populations rely on both acquisition of high-quality methylation data and subsequent precision CGI mapping based on the Chinese T2T reference.

The chromosome-level assembly of the wild diploid alfalfa genome provides insights into the full landscape of genomic variations between cultivated and wild alfalfa.

Alfalfa (Medicago sativa L.) is one of the most important forage legumes in the world, including autotetraploid (M. sativa ssp. sativa) and diploid alfalfa (M. sativa ssp. caerulea, progenitor of autotetraploid alfalfa). Here, we reported a high-quality genome of ZW0012 (diploid alfalfa, 769 Mb, contig N50 = 5.5 Mb), which was grouped into the Northern group in population structure analysis, suggesting that our genome assembly filled a major gap among the members of M. sativa complex. During polyploidization, large phenotypic differences occurred between diploids and tetraploids, and the genetic information underlying its massive phenotypic variations remains largely unexplored. Extensive structural variations (SVs) were identified between ZW0012 and XinJiangDaYe (an autotetraploid alfalfa with released genome). We identified 71 ZW0012-specific PAV genes and 1296 XinJiangDaYe-specific PAV genes, mainly involved in defence response, cell growth, and photosynthesis. We have verified the positive roles of MsNCR1 (a XinJiangDaYe-specific PAV gene) in nodulation using an Agrobacterium rhizobia-mediated transgenic method. We also demonstrated that MsSKIP23_1 and MsFBL23_1 (two XinJiangDaYe-specific PAV genes) regulated leaf size by transient overexpression and virus-induced gene silencing analysis. Our study provides a high-quality reference genome of an important diploid alfalfa germplasm and a valuable resource of variation landscape between diploid and autotetraploid, which will facilitate the functional gene discovery and molecular-based breeding for the cultivars in the future.

A high quality reference genome of the brown hare

A high quality reference genome of the brown hare

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

Chromosome-scale Genome assembly of the critically endangered White-eared Night-Heron (Gorsachius magnificus)

The White-eared Night-Heron (Gorsachius magnificus, G. magnificus) is a critically endangered heron that is very poorly known and only found in southern China and northern Vietnam, with an estimated population of 250 to 999 mature individuals. However, the lack of a reference genome has hindered the implementation of conservation management efforts. In this study, we present the first high-quality chromosome-scale reference genome, which was assembled by integrating PacBio long-reads sequencing, Illumina paired-end sequencing, and Hi-C technology. The genome has a total length of 1.176 Gb, with a scaffold N50 of 84.77 Mb and a contig N50 of 18.46 Mb. Utilizing Hi-C data, we anchored 99.89% of the scaffold sequences onto 29 pairs of chromosomes. Additionally, we identified 18,062 protein-coding genes in the genome, with 95.00% of which were functionally annotated. Notably, BUSCO assessment confirmed the presence of 97.2% of highly conserved Aves genes within the genome. This chromosome-level genome assembly and annotation will be valuable for future investigating the G. magnificus’s evolutionary adaptation and conservation.

Metabolic and transcriptomic characterization of summer and winter dormancy in the solitary bee, Osmia lignaria

The solitary bee Osmia lignaria is a native pollinator in North America with growing economic importance. The life cycle of O. lignaria provides a unique opportunity to compare the physiological and molecular mechanisms underlying two ecologically contrasting dormancies within the same species. O. lignaria prepupae become dormant during the summer to avoid high temperatures. Shortly after adult eclosion, they enter a second dormancy and overwinter as diapausing adults. To compare these two dormancies, we measured metabolic rates and gene expression across development as bees initiate, maintain, and terminate both prepupal (summer) and adult (overwintering) dormancies. We observed a moderate temperature-independent decrease in gas exchange during both the prepupal dormancy after cocoon spinning (45 %) and during adult diapause after eclosion (60 %). We sequenced and assembled a high-quality reference genome from a single haploid male bee with a contiguous n50 of 5.5 Mbp to facilitate our transcriptomic analysis. The transcriptomes of dormant prepupae and diapausing adults clustered into distinct groups more closely associated with life stage than dormancy status. Membrane transport, membrane-bound cellular components, oxidoreductase activity, glutathione metabolism, and transcription factor activity increased during adult diapause, relative to prepupal dormancy. Further, the transcriptomes of adults in diapause clustered into two groups, supporting multiple phases of diapause during winter. Late adult diapause was associated with gene expression profiles supporting increased insulin/IGF, juvenile hormone, and ecdysone signaling.

Chromosome-level assembly and gene annotation of Decapterus maruadsi genome using Nanopore and Hi-C technologies

Decapterus maruadsi is one of the representative offshore fish in the Western Pacific. Since the last century, it has become a commercially valuable marine fishery species in the Western Pacific region. Despite its high economic value, there is still a lack of high-quality reference genome of D. maruadsi in germplasm resource evaluation research. Here we report a chromosome-level reference genome of D. maruadsi based on Nanopore sequencing and Hi-C technologies. The whole genome was assembled through 169 contigs with a total length of 723.69 Mb and a contig N50 length of 24.67 Mb. By chromosome scaffolding, 23 chromosomes with a total length of 713.58 Mb were constructed. In addition, a total of 199.49 Mb repetitive elements, 33,515 protein-coding genes, and 6,431 ncRNAs were annotated in the reference genome. This reference genome of D. maruadsi will provide a solid theoretical basis not only for the subsequent development of genomic resources of D. maruadsi but also for the formulation of policies related to the protection of D. maruadsi.

The high-quality genome of Cryptotaenia japonica and comparative genomics analysis reveals anthocyanin biosynthesis in Apiaceae.

Cryptotaenia japonica, a traditional medicinal and edible vegetable crops, is well-known for its attractive flavors and health care functions. As a member of the Apiaceae family, the evolutionary trajectory and biological properties of C. japonica are not clearly understood. Here, we first reported a high-quality genome of C.japonica with a total length of 427 Mb and N50 length 50.76 Mb, was anchored into 10 chromosomes, which confirmed by chromosome (cytogenetic) analysis. Comparative genomic analysis revealed C. japonica exhibited low genetic redundancy, contained a higher percentage of single-cope gene families. The homoeologous blocks, Ks, and collinearity were analyzed among Apiaceae species contributed to the evidence that C. japonica lacked recent species-specific WGD. Through comparative genomic and transcriptomic analyses of Apiaceae species, we revealed the genetic basis of the production of anthocyanins. Several structural genes encoding enzymes and transcription factor genes of the anthocyanin biosynthesis pathway in different species were also identified. The CjANSa, CjDFRb, and CjF3H gene might be the target of Cjaponica_2.2062 (bHLH) and Cjaponica_1.3743 (MYB). Our findings provided a high-quality reference genome of C. japonica and offered new insights into Apiaceae evolution and biology.

Chromosome-level genome assembly of humpback grouper using PacBio HiFi reads and Hi-C technologies

The humpback grouper (Cromileptes altivelis), a medium-sized coral reef teleost, is a naturally rare species distributed in the tropical waters of the Indian and Pacific Oceans. It has high market value, but artificial reproduction and breeding remain limited and need to be improved. Here, we assembled the genome with 1.08 Gb, with a contig N50 of 43.78 Mb. A total of 96.59% of the assembly anchored to 24 pseudochromosomes using Hi-C technology. It contained 24,442 protein-coding sequences, of which 99.3% were functionally annotated. The completeness of the assembly was estimated to be 97.3% using BUSCO. The phylogenomic analysis suggested that humpback grouper should be classified into the genus Epinephelus rather than Cromileptes. The comparative genomic analysis revealed that the gene families related to circadian entrainment were significantly expanded. The high-quality reference genome provides useful genomic tools for exploiting the genomic resource of humpback grouper and supports the functional genomic study of this species in the future.

Multiple Origins of Bioluminescence in Beetles and Evolution of Luciferase Function.

Bioluminescence in beetles has long fascinated biologists, with diverse applications in biotechnology. To date, however, our understanding of its evolutionary origin and functional variation mechanisms remains poor. To address these questions, we obtained high-quality reference genomes of luminous and nonluminous beetles in 6 Elateroidea families. We then reconstructed a robust phylogenetic relationship for all luminous families and related nonluminous families. Comparative genomic analyses and biochemical functional experiments suggested that gene evolution within Elateroidea played a crucial role in the origin of bioluminescence, with multiple parallel origins observed in the luminous beetle families. While most luciferase-like proteins exhibited a conserved nonluminous amino acid pattern (TLA346 to 348) in the luciferin-binding sites, luciferases in the different luminous beetle families showed divergent luminous patterns at these sites (TSA/CCA/CSA/LVA). Comparisons of the structural and enzymatic properties of ancestral, extant, and site-directed mutant luciferases further reinforced the important role of these sites in the trade-off between acyl-CoA synthetase and luciferase activities. Furthermore, the evolution of bioluminescent color demonstrated a tendency toward hypsochromic shifts and variations among the luminous families. Taken together, our results revealed multiple parallel origins of bioluminescence and functional divergence within the beetle bioluminescent system.

Toward genome assemblies for all marine vertebrates: current landscape and challenges.

Marine vertebrate biodiversity is fundamental to ocean ecosystem health but is threatened by climate change, overharvesting, and habitat degradation. High-quality reference genomes are valuable foundational scientific resources that can inform conservation efforts. Consequently, global consortia are striving to produce reference genomes for representatives of all life. Here, we summarize the current landscape of available marine vertebrate reference genomes, including their phylogenetic diversity and geographic hotspots of production. We discuss key logistical and technical challenges that remain to be overcome if we are to realize the vision of a comprehensive reference genome library of all marine vertebrates.

The genome of the black-footed cat: Revealing a rich natural history and urgent conservation priorities for small felids.

Habitat degradation and loss of genetic diversity are common threats faced by almost all of today's wild cats. Big cats, such as tigers and lions, are of great concern and have received considerable conservation attention through policies and international actions. However, knowledge of and conservation actions for small wild cats are lagging considerably behind. The black-footed cat, Felis nigripes, one of the smallest felid species, is experiencing increasing threats with a rapid reduction in population size. However, there is a lack of genetic information to assist in developing effective conservation actions. A de novo assembly of a high-quality chromosome-level reference genome of the black-footed cat was made, and comparative genomics and population genomics analyses were carried out. These analyses revealed that the most significant genetic changes in the evolution of the black-footed cat are the rapid evolution of sensory and metabolic-related genes, reflecting genetic adaptations to its characteristic nocturnal hunting and a high metabolic rate. Genomes of the black-footed cat exhibit a high level of inbreeding, especially for signals of recent inbreeding events, which suggest that they may have experienced severe genetic isolation caused by habitat fragmentation. More importantly, inbreeding associated with two deleterious mutated genes may exacerbate the risk of amyloidosis, the dominant disease that causes mortality of about 70% of captive individuals. Our research provides comprehensive documentation of the evolutionary history of the black-footed cat and suggests that there is an urgent need to investigate genomic variations of small felids worldwide to support effective conservation actions.

The telomere-to-telomere (T2T) genome of Peucedanum praeruptorum Dunn provides insights into the genome evolution and coumarin biosynthesis.

Traditional Chinese medicine has used Peucedanum praeruptorum Dunn (Apiaceae) for a long time. Various coumarins, including the significant constituents praeruptorin (A-E), are the active constituents in the dried roots of P. praeruptorum. Previous transcriptomic and metabolomic studies have attempted to elucidate the distribution and biosynthetic network of these medicinal-valuable compounds. However, the lack of a high-quality reference genome impedes an in-depth understanding of genetic traits and thus the development of better breeding strategies. A telomere-to-telomere (T2T) genome was assembled for P. praeruptorum by combining PacBio HiFi, ONT ultra-long, and Hi-C data. The final genome assembly was approximately 1.798 Gb, assigned to 11 chromosomes with genome completeness >98%. Comparative genomic analysis suggested that P. praeruptorum experienced 2 whole-genome duplication events. By the transcriptomic and metabolomic analysis of the coumarin metabolic pathway, we presented coumarins' spatial and temporal distribution and the expression patterns of critical genes for its biosynthesis. Notably, the COSY and cytochrome P450 genes showed tandem duplications on several chromosomes, which may be responsible for the high accumulation of coumarins. A T2T genome for P. praeruptorum was obtained, providing molecular insights into the chromosomal distribution of the coumarin biosynthetic genes. This high-quality genome is an essential resource for designing engineering strategies for improving the production of these valuable compounds.

A high-quality chromosomal genome assembly of the sea cucumber Chiridota heheva and its hydrothermal adaptation.

Chiridota heheva is a cosmopolitan holothurian well adapted to diverse deep-sea ecosystems, especially chemosynthetic environments. Besides high hydrostatic pressure and limited light, high concentrations of metal ions also represent harsh conditions in hydrothermal environments. Few holothurian species can live in such extreme conditions. Therefore, it is valuable to elucidate the adaptive genetic mechanisms of C. heheva in hydrothermal environments. Herein, we report a high-quality reference genome assembly of C. heheva from the Kairei vent, which is the first chromosome-level genome of Apodida. The chromosome-level genome size was 1.43 Gb, with a scaffold N50 of 53.24 Mb and BUSCO completeness score of 94.5%. Contig sequences were clustered, ordered, and assembled into 19 natural chromosomes. Comparative genome analysis found that the expanded gene families and positively selected genes of C. heheva were involved in the DNA damage repair process. The expanded gene families and the unique genes contributed to maintaining iron homeostasis in an iron-enriched environment. The positively selected gene RFC2 with 10 positively selected sites played an essential role in DNA repair under extreme environments. This first chromosome-level genome assembly of C. heheva reveals the hydrothermal adaptation of holothurians. As the first chromosome-level genome of order Apodida, this genome will provide the resource for investigating the evolution of class Holothuroidea.

High-quality assembly of the T2T genome for Isodon rubescens f. lushanensis reveals genomic structure variations between 2 typical forms of Isodon rubescens.

Rabdosiae rubescentis herba (Isodon rubescens) is widely used as a folk medicine to treat esophageal cancer and sore throat in China. Its germplasm resources are abundant in China, with I. rubescens (Hemsl.) Hara and I. rubescens f. lushanensis as 2 typical forms. I. rubescens (Hemsl.) Hara is featured by biosynthesis of the diterpenoid oridonin with strong anticancer activity, while I. rubescens f. lushanensis produces another diterpenoid with anticancer activity, lushanrubescensin. However, the biosynthetic pathways of both still need to be fully understood. In particular, little is known about the genetic background of I. rubescens f. lushanensis. We used Pacific Biosciences (PacBio) single-molecule real-time and Nanopore Ultra-long sequencing platforms, respectively, and obtained 139.07 Gb of high-quality data, with a sequencing depth of about 328×. We also obtained a high-quality reference genome for I. rubescens f. lushanensis, with a genome size of 349 Mb and a contig N50 of 28.8 Mb. The heterozygosity of the genome is 1.7% and the repeatability is 83.43%. In total, 34,865 protein-coding genes were predicted. Moreover, we found that most of the variant or unique genes in the diterpenoid synthesis pathways of I. rubescens f. lushanensis and I. rubescens (Hemsl.) Hara were enriched in diterpene synthases. We provide the first genome sequence and gene annotation for the I. rubescens f. lushanensis, which provides molecular evidence for understanding the chemotypic differences of I. rubescens.

Chromosome-level genome of Thymus mandschuricus reveals molecular mechanism of aroma compounds biosynthesis.

Thymus mandschuricus is an aromatic and medicinal plant with notable antibacterial and antioxidant properties. However, traditional breeding methods rely on phenotypic selection due to a lack of molecular resources. A high-quality reference genome is crucial for marker-assisted breeding, genome editing, and molecular genetics. We utilized PacBio and Hi-C technologies to generate a high-quality chromosome-level reference genome for T. mandschuricus, with a size of 587.05 Mb and an N50 contig size of 8.41 Mb. The assembled genome contained 29,343 predicted protein-coding genes, and evidence of two distinct whole-genome duplications in T. mandschuricus was discovered. Comparative genomic analysis revealed rapid evolution of genes involved in phenylpropanoid biosynthesis and the CYP450 gene family in T. mandschuricus. Additionally, we reconstructed the gene families of terpenoid biosynthesis structural genes, such as TPS, BAHD, and CYP, and identified regulatory networks controlling the expression of aroma-synthesis genes by integrating transcriptome data from various organs and developmental stages. We discovered that hormones and transcription factors may collaborate in controlling aroma-synthesis gene expression. This study provides the first high-quality genome sequence and gene annotation for T. mandschuricus, an indigenous thyme species unique to China. The genome assembly and the comprehension of the genetic basis of fragrance synthesis acquired from this research could potentially serve as targets for future breeding programs and functional studies.

Chromosome-scale assembly and gene editing of Solanum americanum genome reveals the basis for thermotolerance and fruit anthocyanin composition.

Solanum americanum serves as a promising source of resistance genes against potato late blight and is considered as a leafy vegetable for complementary food and nutrition. The limited availability of high-quality genome assemblies and gene annotations has hindered the exploration and exploitation of stress-resistance genes in S. americanum. Here, we present a chromosome-level genome assembly of a thermotolerant S. americanum ecotype and identify a crucial heat-inducible transcription factor gene, SaHSF17, essential for heat tolerance. The CRISPR/Cas9 system-mediated knockout of SaHSF17 results in remarkably reduced thermotolerance in S. americanum, exhibiting a significant suppression of multiple HSP gene expressions under heat treatment. Furthermore, our transcriptome analysis and anthocyanin component investigation of fruits indicated that delphinidins are the major anthocyanins accumulated in the mature dark-purple fruits. The accumulation of delphinidins and other pigment components during fruit ripening in S. americanum coincides with the transcriptional regulation of key genes, particularly the F3'5'H and F3'H genes, in the anthocyanin biosynthesis pathway. By integrating existing knowledge, the development of this high-quality reference genome for S. americanum will facilitate the identification and utilization of novel abiotic and biotic stress-resistance genes for improvement of Solanaceae and other crops.

Advances in genomics and genome editing for improving strawberry (Fragaria ×ananassa).

The cultivated strawberry, Fragaria ×ananassa, is a recently domesticated fruit species of economic interest worldwide. As such, there is significant interest in continuous varietal improvement. Genomics-assisted improvement, including the use of DNA markers and genomic selection have facilitated significant improvements of numerous key traits during strawberry breeding. CRISPR/Cas-mediated genome editing allows targeted mutations and precision nucleotide substitutions in the target genome, revolutionizing functional genomics and crop improvement. Genome editing is beginning to gain traction in the more challenging polyploid crops, including allo-octoploid strawberry. The release of high-quality reference genomes and comprehensive subgenome-specific genotyping and gene expression profiling data in octoploid strawberry will lead to a surge in trait discovery and modification by using CRISPR/Cas. Genome editing has already been successfully applied for modification of several strawberry genes, including anthocyanin content, fruit firmness and tolerance to post-harvest disease. However, reports on many other important breeding characteristics associated with fruit quality and production are still lacking, indicating a need for streamlined genome editing approaches and tools in Fragaria ×ananassa. In this review, we present an overview of the latest advancements in knowledge and breeding efforts involving CRISPR/Cas genome editing for the enhancement of strawberry varieties. Furthermore, we explore potential applications of this technology for improving other Rosaceous plant species.

High-speed whole-genome sequencing of a Whippet: Rapid chromosome-level assembly and annotation of an extremely fast dog's genome.

The time required for genome sequencing and de novo assembly depends on the interaction between laboratory work, sequencing capacity, and the bioinformatics workflow, often constrained by external sequencing services. Bringing together academic biodiversity institutes and a medical diagnostics company with extensive sequencing capabilities, we aimed at generating a high-quality mammalian de novo genome in minimal time. We present the first chromosome-level genome assembly of the Whippet, using PacBio long-read high-fidelity sequencing and reference-guided scaffolding. The final assembly has a contig N50 of 55 Mbp and a scaffold N50 of 65.7 Mbp. The total assembly length is 2.47 Gbp, of which 2.43 Gpb were scaffolded into 39 chromosome-length scaffolds. Annotation using mammalian genomes and transcriptome data yielded 28,383 transcripts, 90.9% complete BUSCO genes, and identified 36.5% repeat content. Sequencing, assembling, and scaffolding the chromosome-level genome of the Whippet took less than a week, adding another high-quality reference genome to the available sequences of domestic dog breeds.