Comprehensive analysis of alternative splicing in rice and comparative analyses with Arabidopsis

Genomic insights and the conservation potential of captive breeding: The case of Chinese alligator

Despite 40 years of conservation of the critically endangered Chinese alligator (Alligator sinensis), the genomic underpinnings of its status remained uncharted. Genome sequencing data of 244 individuals uncovered relatively low overall genomic diversity/heterozygosity and long runs of homozygosity, with captive populations exhibiting higher heterozygosity and smaller inbreeding coefficients compared to wild individuals. The decreased level of inbreeding in the captive population demonstrates the contribution of the large captive breeding population. The estimated recent effective population size was around a few dozen. To combat challenges of inbreeding depression and reduced adaptability, we used genome-wide SNP-based kinship analysis on captive populations to enable a genome-informed breeding program that minimizes inbreeding. Long-term field monitoring revealed that the Chinese government greatly advanced the conservation of A. sinensis through conservation measures and reintroduction programs. Our research enriches the understanding of the Chinese alligator's genetic landscape, offering invaluable genomic resources for breeding and conservation strategies.

Chromosome-scale genome assembly of Zoysia japonica uncovers cold tolerance candidate genes

Zoysiagrass stands out as a crucial native turfgrass due to its exceptional abiotic stress tolerance, extensive adaptability, and high ornamental value. In this study, we generated a high-quality chromosome-level genome assembly of Compadre (COM) zoysiagrass, leveraging PacBio SMRT sequencing and Hi-C scaffolding technologies. The resulting genome assembly (312.42 Mb) is anchored on 20 chromosomes, with a Scaffold N50 of 18.72 Mb. In total, 49,074 genes and 306,768 repeat sequences were annotated in the assembled genome. The first chromosome-scale genome of Zoysia japonica 'Compadre' provides a critical genetic resource for cold-tolerant turfgrass breeding through identifying stress-responsive candidate genes. Additionally, we have successfully established a cell nucleus extraction and library construction protocol tailored for zoysiagrass ATAC-seq technology, and a total of 80 low temperature tolerance candidate genes were preliminarily identified via ATAC-seq and RNA-seq profiling, thereby initiating the exploration of turfgrass epigenomics.

Jan and mini-Jan, a model system for potato functional genomics

Potato (Solanum tuberosum) is the third-most important food crop in the world. Although the potato genome has been fully sequenced, functional genomics research of potato lags behind that of other major food crops, largely due to the lack of a model experimental potato line. Here, we present a diploid potato line, 'Jan,' which possesses all essential characteristics for facile functional genomics studies. Jan exhibits a high level of homozygosity after seven generations of self-pollination. Jan is vigorous, highly fertile and produces tubers with outstanding traits. Additionally, it demonstrates high regeneration rates and excellent transformation efficiencies. We generated a chromosome-scale genome assembly for Jan, annotated its genes and identified syntelogs relative to the potato reference genome assembly DMv6.1 to facilitate functional genomics. To miniaturize plant architecture, we developed two 'mini-Jan' lines with compact and dwarf plant stature through CRISPR/Cas9-mediated mutagenesis targeting the Dwarf and Erecta genes involved in growth. One mini-Jan mutant, mini-JanE, is fully fertile and will permit higher-throughput studies in limited growth chamber and greenhouse space. Thus, Jan and mini-Jan offer a robust model system that can be leveraged for gene editing and functional genomics research in potato.

Decoding Pecan's Fungal Foe: A Genomic Insight into Colletotrichum plurivorum Isolate W-6

Pecan (Carya illinoinensis) is a world-renowned nut crop that is highly favored by consumers for its high content of healthy nutrients. For a long time, anthracnose has severely threatened the yield and quality of pecan, causing significant economic losses to the global pecan industry. Here, we report the 54.57-Mb gapless chromosome-level assembly of the pathogenic ascomycetes Colletotrichum plurivorum isolate W-6 from pecan plantations in Southeast China. Six of 12 chromosomes contain, at least, telomeric repeats (CCCTAA)n or (TTAGGG)n at one end. A total of 14,343 protein-coding genes were predicted. Pathogenicity- and virulence-related annotations revealed 137 to 4558 genes associated with the TCDB, PHI, Cyt_P450, DFVF, effector, and secretome databases, respectively. A comparative analysis of isolate W-6, together with 51 other Colletotrichum strains, reveled 13 genes unique to the Orchidearum complex to which isolate W-6 belongs, highlighting the major facilitator superfamily transporters. The detailed analyses of MFS transporters associated with secondary metabolite gene clusters in isolate W-6 led to the identification and protein structure analyses of two key virulence factor candidates in DHA1 subclass, prlG and azaK, which were reported as efflux transporters of antibiotics in other pathogenic fungi. The assembly and further functional investigation of two pathogenic genes identified here potentially provide important resources for better understanding the biology and lifestyle of Colletotrichum and pave the way for designing more efficient strategies to control anthracnose in pecan plantations.

Chromosome level genome assembly and annotation of the Swanson rainbow trout homozygous line

The genome of the Swanson doubled haploid (DH) YY male line of rainbow trout was de novo assembled using the Canu pipeline, high-coverage PacBio long-read sequence data, Bionano optical maps, and Hi-C proximity ligation sequence data, resulting in 29 major scaffolds aligning with the karyotype of the Swanson line (2 N = 58). This assembly, totaling 2.3 Gb with an N50 of 52.4 Mb, represents approximately 95% of the genome in 29 chromosome sequences with only 109 gaps between scaffolds. Notably, corrections to previous errors in the Swanson line genome assembly were made, including the identification of a double large inversion on the Omy05 chromosome (~57 Mb), the absence of the Omy20 inversion between the Arlee and Swanson assemblies, and the discovery of a ~6.7 Mb inversion on Omy26. This comprehensive assembly contributes to refining the rainbow trout reference genome and serves as a valuable resource for future genetic studies within this species.

Evidence for evolution of a new sex chromosome within the haploid-dominant Marchantiales plant lineage

Sex chromosomes have evolved independently in numerous lineages across the Tree of Life, in both diploid-dominant species, including many animals and plants, and the less studied haploid-dominant plants and algae. Strict genetic sex determination ensures that individuals reproduce by outcrossing. However, species with separate sexes (termed dioecy in diploid plants, and dioicy in haploid plants) may sometimes evolve different sex systems, and become monoicous, with the ability to self-fertilize. Here, we studied dioicy-monoicy transitions in the ancient liverwort haploid-dominant plant lineage, using three telomere-to-telomere gapless chromosome-scale reference genome assemblies from the Ricciaceae group of Marchantiales. Ancestral liverworts are believed to have been dioicous, with U and V chromosomes (chromosome 9) determining femaleness and maleness, respectively. We confirm the finding that monoicy in Ricciocarpos natans evolved from a dioicous ancestor, and most ancestrally U chromosomal genes have been retained on autosomes in this species. We also describe evidence suggesting the possible re-evolution of dioicy in the genus Riccia, with probable de novo establishment of a sex chromosome from an autosome (chromosome 5), and further translocations of genes from the new sex chromosome to autosomes. Our results also indicated that micro-chromosomes are consistent genomic features, and may have evolved independently from sex chromosomes in Ricciocarpos and Riccia lineages.

Near-complete genome and infection transcriptomes of the maize leaf and sheath spot pathogen Epicoccum sorghinum

Maize leaf and sheath spot disease caused by Epicoccum sorghinum is an emerging disease of maize in China. To disentangle the molecular pathogenesis, we sequenced the genome and infection transcriptomes of the E. sorghinum strain NJC07. The genome was sequenced on Oxford Nanopore GridION and Illumina NovaSeq 6000, producing a near-complete gapless nuclear genome assembly of 32.69 Mb at 285.20-fold depth, comprising 23 contigs (including 12 full-length chromosomes) with an N50 contig number/length of 6/1.66 Mb, and a complete mitochondrial genome assembly of 61.24 kb. The nuclear genome contains 11,779 protein-coding genes, including those predicted to encode potential virulence/pathogenicity factors, such as effectors and carbohydrate-active enzymes. Temporal RNA-Seq analysis revealed that 4,058 of the 11,779 genes were induced during maize infection, with a subset potentially implicated in fungal invasion and colonization of maize plants. Together, the genomic and transcriptomic data generated in the study provide a valuable foundation for the functional analysis of virulence and pathogenicity factors, offering critical insights into the molecular mechanisms driving E. sorghinum pathogenesis on maize.

Jan and mini-Jan, a model system for potato functional genomics

Potato (Solanum tuberosum) is the third most important food crop in the world. Although the potato genome has been fully sequenced, functional genomics research of potato lags relative to other major food crops due primarily to the lack of a model experimental potato line. Here, we present a diploid potato line, 'Jan', which possesses all essential characteristics for facile functional genomics studies. Jan has a high level of homozygosity after seven generations of self-pollination. Jan is vigorous and highly fertile with outstanding tuber traits, high regeneration rates, and excellent transformation efficiencies. We generated a chromosome-scale genome assembly for Jan, annotated genes, and identified syntelogs relative to the potato reference genome assembly DMv6.1 to facilitate functional genomics. To miniaturize plant architecture, we developed two "mini-Jan" lines with compact and dwarf plant stature using CRISPR/Cas9-mediated mutagenesis targeting the Dwarf and Erecta genes related to growth. Mini-Jan mutants are fully fertile and will permit higher-throughput studies in limited growth chamber and greenhouse space. Thus, Jan and mini-Jan provide an outstanding model system that can be leveraged for gene editing and functional genomics research in potato.

A Chromosome-Level Genome Assembly of Chiton Acanthochiton rubrolineatus (Chitonida, Polyplacophora, Mollusca)

(1) Background: Chitons (Mollusca, Polyplacophora) are relatively primitive species in Mollusca that allow the study of biomineralization. Although mitochondrial genomes have been isolated from Polyplacophora, there is no genomic information at the chromosomal level; (2) Methods: Here we report a chromosome-level genome assembly for Acanthochiton rubrolineatus using PacBio (Pacific Biosciences, United States) reads and high-throughput chromosome conformation capture (Hi-C) data; (3) Results: The assembly spans 1.08 Gb with a contig N50 of 3.63 Mb and 99.97% of the genome assigned to eight chromosomes. Among the 32,291 predicted genes, 76.32% had functional predictions. The divergence time of Brachiopoda and Mollusca was ~550.8 Mya (million years ago), and that of A. rubrolineatus and other mollusks was ~548.5 Mya; (4) Conclusions: This study not only offers high-quality reference sequences for the Acanthochiton rubrolineatus genome, but also establishes groundwork for investigating the mechanisms of Polyplacophora biomineralization and its evolutionary history. This research will aid in uncovering the genetic foundations of molluscan adaptations across diverse environments.

The Chromosome-level Genome Provides Insights into the Evolution and Adaptation of Extreme Aggression

Extremely aggressive behavior, as the special pattern, is rare in most species and characteristic as contestants severely injured or killed ending the combat. Current studies of extreme aggression are mainly from the perspectives of behavioral ecology and evolution, while lacked the aspects of molecular evolutionary biology. Here, a high-quality chromosome-level genome of the parasitoid Anastatus disparis was provided, in which the males exhibit extreme mate-competition aggression. The integrated multiomics analysis highlighted that neurotransmitter dopamine overexpression, energy metabolism (especially from lipid), and antibacterial activity are likely major aspects of evolutionary formation and adaptation for extreme aggression in A. disparis. Conclusively, our study provided new perspectives for molecular evolutionary studies of extreme aggression as well as a valuable genomic resource in Hymenoptera.

A Genomic Sequence Resource of Diaporthe mahothocarpus GZU-Y2 Causing Leaf Spot Blight in Camellia oleifera

Diaporthe mahothocarpus GZU-Y2, a new pathogen responsible for leaf spot blight disease, leads to significant damage and economic losses in some Camellia oleifera plantations. The current study annotated the genome of the D. mahothocarpus isolate GZU-Y2 to advance our knowledge of the pathogen and facilitate improved disease management of leaf spot blight. The initial Pacbio-Illumina hybrid draft genome for GZU-Y2 resulted in a high-quality assembly with 62 contigs, characterized by an N50 length of 7.07 Mb. The complete genome of isolate GZU-Y2 was 58.97 Mbp, with a GC content of 50.65%. Importantly, the assembly exhibits remarkable integrity, with 97.93% of complete BUSCO validating genome completeness. The prediction results showed that a total of 15,918 protein-coding genes were annotated using multiple bioinformatics databases. The genome assembly and annotation resource reported here will be useful for the further study of fungal infection mechanisms and pathogen-host interaction.

New insights into decoding the lifestyle of endophytic Fusarium lateritium Fl617 via comparing genomes

Fungal-plant interactions have persisted for 460 million years, and almost all terrestrial plants on Earth have endophytic fungi. However, the mechanism of symbiosis between endophytic fungi and host plants has been inconclusive. In this dissertation, we used a strain of endophytic Fusarium lateritium (Fl617), which was found in the previous stage to promote disease resistance in tomato, and selected the pathogenic Fusarium oxysporum Fo4287 and endophytic Fusarium oxysporum Fo47, which are in the same host and the closest relatives of Fl617, to carry out a comparative genomics analysis of the three systems and to provide a new perspective for the elucidation of the special lifestyle of the fungal endophytes. We found that endophytic F. lateritium has a smaller genome, fewer clusters and genes associated with pathogenicity, and fewer plant cell wall degrading enzymes (PCWDEs). There were also relatively fewer secondary metabolisms and typical Fusarium spp. toxins, and a lack of the key Fusarium spp. pathogenicity factor, secreted in xylem (SIX), but the endophytic fungi may be more sophisticated in their regulation of the colonization process. It is hypothesized that the endophytic fungi may have maintained their symbiosis with plants due to the relatively homogeneous microenvironment in plants for a long period of time, considering only plant interactions and discarding the relevant pathogenicity factors, and that their endophytic evolutionary tendency may tend to be genome streamlining and to enhance the fineness of the regulation of plant interactions, thus maintaining their symbiotic status with plants.

Whole-genome sequence of Sclerotium delphinii, a pathogenic fungus of Dendrobium officinale southern blight

Dendrobium officinale is a rare and precious medicinal plant. Southern blight is a destructive disease in the artificial cultivation of D. officinale, and one of its pathogens is Sclerotium delphinii. S. delphinii is a phytopathogenic fungus with a wide host range with extremely strong pathogenicity. In this study, S. delphinii was isolated from D. officinale with southern blight. Subsequently, this specific strain underwent thorough whole-genome sequencing using the PacBio Sequel II platform, which employed single-molecule real-time (SMRT) technology. Comprehensive annotations were obtained through functional annotation of protein sequences using various publicly available databases. The genome of S. delphinii measures 73.66 Mb, with an N90 contig size of 2,707,110 bp, and it contains 18,506 putative predictive genes. This study represents the first report on the genome size assembly and annotation of S. delphinii, making it the initial species within the Sclerotium genus to undergo whole-genome sequencing, which can provide solid data and a theoretical basis for further research on the pathogenesis, omics of S. delphinii.

Transcriptional Modulation during Photomorphogenesis in Rice Seedlings

Light is one of the most important factors regulating plant gene expression patterns, metabolism, physiology, growth, and development. To explore how light may induce or alter transcript splicing, we conducted RNA-Seq-based transcriptome analyses by comparing the samples harvested as etiolated seedlings grown under continuous dark conditions vs. the light-treated green seedlings. The study aims to reveal differentially regulated protein-coding genes and novel long noncoding RNAs (lncRNAs), their light-induced alternative splicing, and their association with biological pathways. We identified 14,766 differentially expressed genes, of which 4369 genes showed alternative splicing. We observed that genes mapped to the plastid-localized methyl-erythritol-phosphate (MEP) pathway were light-upregulated compared to the cytosolic mevalonate (MVA) pathway genes. Many of these genes also undergo splicing. These pathways provide crucial metabolite precursors for the biosynthesis of secondary metabolic compounds needed for chloroplast biogenesis, the establishment of a successful photosynthetic apparatus, and photomorphogenesis. In the chromosome-wide survey of the light-induced transcriptome, we observed intron retention as the most predominant splicing event. In addition, we identified 1709 novel lncRNA transcripts in our transcriptome data. This study provides insights on light-regulated gene expression and alternative splicing in rice.

Genomic and single-cell analyses reveal genetic signatures of swimming pattern and diapause strategy in jellyfish

Jellyfish exhibit innovative swimming patterns that contribute to exploring the origins of animal locomotion. However, the genetic and cellular basis of these patterns remains unclear. Herein, we generated chromosome-level genome assemblies of two jellyfish species, Turritopsis rubra and Aurelia coerulea, which exhibit straight and free-swimming patterns, respectively. We observe positive selection of numerous genes involved in statolith formation, hair cell ciliogenesis, ciliary motility, and motor neuron function. The lineage-specific absence of otolith morphogenesis- and ciliary movement-related genes in T. rubra may be associated with homeostatic structural statocyst loss and straight swimming pattern. Notably, single-cell transcriptomic analyses covering key developmental stages reveal the enrichment of diapause-related genes in the cyst during reverse development, suggesting that the sustained diapause state favours the development of new polyps under favourable conditions. This study highlights the complex relationship between genetics, locomotion patterns and survival strategies in jellyfish, thereby providing valuable insights into the evolutionary lineages of movement and adaptation in the animal kingdom.

High-quality chromosome-level genome assembly and multi-omics analysis of rosemary (Salvia rosmarinus) reveals new insights into the environmental and genome adaptation

High-quality genome of rosemary (Salvia rosmarinus) represents a valuable resource and tool for understanding genome evolution and environmental adaptation as well as its genetic improvement. However, the existing rosemary genome did not provide insights into the relationship between antioxidant components and environmental adaptability. In this study, by employing Nanopore sequencing and Hi-C technologies, a total of 1.17 Gb (97.96%) genome sequences were mapped to 12 chromosomes with 46 121 protein-coding genes and 1265 non-coding RNA genes. Comparative genome analysis reveals that rosemary had a closely genetic relationship with Salvia splendens and Salvia miltiorrhiza, and it diverged from them approximately 33.7 million years ago (MYA), and one whole-genome duplication occurred around 28.3 MYA in rosemary genome. Among all identified rosemary genes, 1918 gene families were expanded, 35 of which are involved in the biosynthesis of antioxidant components. These expanded gene families enhance the ability of rosemary adaptation to adverse environments. Multi-omics (integrated transcriptome and metabolome) analysis showed the tissue-specific distribution of antioxidant components related to environmental adaptation. During the drought, heat and salt stress treatments, 36 genes in the biosynthesis pathways of carnosic acid, rosmarinic acid and flavonoids were up-regulated, illustrating the important role of these antioxidant components in responding to abiotic stresses by adjusting ROS homeostasis. Moreover, cooperating with the photosynthesis, substance and energy metabolism, protein and ion balance, the collaborative system maintained cell stability and improved the ability of rosemary against harsh environment. This study provides a genomic data platform for gene discovery and precision breeding in rosemary. Our results also provide new insights into the adaptive evolution of rosemary and the contribution of antioxidant components in resistance to harsh environments.

Chromosome-scale genome assembly of medicinal plant Tinospora sagittata (Oliv.) Gagnep. from the Menispermaceae family

Tinospora sagittata (Oliv.) Gagnep. is an important medicinal tetraploid plant in the Menispermaceae family. Its tuber, Radix Tinosporae, used in traditional Chinese medicine, is rich in diterpenoids and benzylisoquinoline alkaloids (BIAs). To enhance our understanding of medicinal compounds' biosynthesis and Menispermaceae's evolution, we herein report assembling a high-quality chromosome-scale genome with both PacBio HiFi and Illumina sequencing technologies. PacBio Sequel II generated 2.5 million circular consensus sequencing (CCS) reads, and a hybrid assembly strategy with Illumina sequencing resulted in 4483 contigs. The assembled genome size was 2.33 Gb, consisting of 4070 scaffolds (N50 = 42.06 Mb), of which 92.05% were assigned to 26 pseudochromosomes. T. sagittata's chromosomal-scale genome assembly, the first species in Menispermaceae, aids Menispermaceae evolution and T. sagittata's secondary metabolites biosynthesis understanding.

Concordant Gene Expression and Alternative Splicing Regulation under Abiotic Stresses in Arabidopsis

The current investigation endeavors to identify differentially expressed alternatively spliced (DAS) genes that exhibit concordant expression with splicing factors (SFs) under diverse multifactorial abiotic stress combinations in Arabidopsis seedlings. SFs serve as the post-transcriptional mechanism governing the spatiotemporal dynamics of gene expression. The different stresses encompass variations in salt concentration, heat, intensive light, and their combinations. Clusters demonstrating consistent expression profiles were surveyed to pinpoint DAS/SF gene pairs exhibiting concordant expression. Through rigorous selection criteria, which incorporate alignment with documented gene functionalities and expression patterns observed in this study, four members of the serine/arginine-rich (SR) gene family were delineated as SFs concordantly expressed with six DAS genes. These regulated SF genes encompass cactin, SR1-like, SR30, and SC35-like. The identified concordantly expressed DAS genes encode diverse proteins such as the 26.5 kDa heat shock protein, chaperone protein DnaJ, potassium channel GORK, calcium-binding EF hand family protein, DEAD-box RNA helicase, and 1-aminocyclopropane-1-carboxylate synthase 6. Among the concordantly expressed DAS/SF gene pairs, SR30/DEAD-box RNA helicase, and SC35-like/1-aminocyclopropane-1-carboxylate synthase 6 emerge as promising candidates, necessitating further examinations to ascertain whether these SFs orchestrate splicing of the respective DAS genes. This study contributes to a deeper comprehension of the varied responses of the splicing machinery to abiotic stresses. Leveraging these DAS/SF associations shows promise for elucidating avenues for augmenting breeding programs aimed at fortifying cultivated plants against heat and intensive light stresses.

Chromosomal-level genome assembly and annotation of the tropical sea cucumber Holothuria scabra

Holothuria scabra, a commercially valuable yet ecologically vulnerable tropical holothuroid, has experienced a severe decline in its wild populations, especially in China. Genomic resources are crucial for the development of effective genomic breeding projects and stock conservation strategies to restore these natural populations. Until now, a high-quality, chromosome-level reference genome for H. scabra has not been available. Here, we employed Oxford Nanopore and Hi-C sequencing technologies to assemble and annotate a high-quality, chromosome-level reference genome of H. scabra. The final genome comprised 31 scaffolds with a total length of 1.19 Gb and a scaffold N50 length of 53.52 Mb. Remarkably, 1,191.67 Mb (99.95%) of the sequences were anchored to 23 pseudo-chromosomes, with the longest one spanning 79.75 Mb. A total of 34,418 protein-coding genes were annotated in the final genome, with BUSCO analysis revealing 98.01% coverage of metazoa_odb10 genes, marking a significant improvement compared to the previous report. These chromosome-level sequences and annotations will provide an essential genomic basis for further investigation into molecular breeding and conservation management of H. scabra.

Genome assemblies of 11 bamboo species highlight diversification induced by dynamic subgenome dominance

Polyploidy (genome duplication) is a pivotal force in evolution. However, the interactions between parental genomes in a polyploid nucleus, frequently involving subgenome dominance, are poorly understood. Here we showcase analyses of a bamboo system (Poaceae: Bambusoideae) comprising a series of lineages from diploid (herbaceous) to tetraploid and hexaploid (woody), with 11 chromosome-level de novo genome assemblies and 476 transcriptome samples. We find that woody bamboo subgenomes exhibit stunning karyotype stability, with parallel subgenome dominance in the two tetraploid clades and a gradual shift of dominance in the hexaploid clade. Allopolyploidization and subgenome dominance have shaped the evolution of tree-like lignified culms, rapid growth and synchronous flowering characteristic of woody bamboos as large grasses. Our work provides insights into genome dominance in a remarkable polyploid system, including its dependence on genomic context and its ability to switch which subgenomes are dominant over evolutionary time.

A haplotype-resolved genome provides insight into allele-specific expression in wild walnut (Juglans regia L.)

Wild germplasm resources are crucial for gene mining and molecular breeding because of their special trait performance. Haplotype-resolved genome is an ideal solution for fully understanding the biology of subgenomes in highly heterozygous species. Here, we surveyed the genome of a wild walnut tree from Gongliu County, Xinjiang, China, and generated a haplotype-resolved reference genome of 562.99 Mb (contig N50 = 34.10 Mb) for one haplotype (hap1) and 561.07 Mb (contig N50 = 33.91 Mb) for another haplotype (hap2) using PacBio high-fidelity (HiFi) reads and Hi-C technology. Approximately 527.20 Mb (93.64%) of hap1 and 526.40 Mb (93.82%) of hap2 were assigned to 16 pseudochromosomes. A total of 41039 and 39744 protein-coding gene models were predicted for hap1 and hap2, respectively. Moreover, 123 structural variations (SVs) were identified between the two haplotype genomes. Allele-specific expression genes (ASEGs) that respond to cold stress were ultimately identified. These datasets can be used to study subgenome evolution, for functional elite gene mining and to discover the transcriptional basis of specific traits related to environmental adaptation in wild walnut.

Genomics insights into flowering and floral pattern formation: regional duplication and seasonal pattern of gene expression in Camellia

BACKGROUND

The formation and domestication of ornamental traits are influenced by various aspects, such as the recognition of esthetic values and cultural traditions. Camellia japonica is widely appreciated and domesticated around the world mainly due to its rich variations in ornamental traits. Ornamental camellias have a diverse range of resources, including different bud variations from Camellia spp. as well as inter- and intra- specific hybridization. Despite research on the formation of ornamental traits, a basic understanding of their genetics and genomics is still lacking.

RESULTS

Here, we report the chromosomal-level reference genome of C. japonica through combining multiple DNA-sequencing technologies and obtain a high-density genetic linkage map of 4255 markers by sequencing 98 interspecific F1 hybrids between C. japonica and C. chekiangoleosa. We identify two whole-genome duplication events in C. japonica: one is a shared ancient γ event, and the other is revealed to be specific to genus Camellia. Based on the micro-collinearity analysis, we find large-scale segmental duplication of chromosome 8, resulting to two copies of the AGAMOUS loci, which may play a key role in the domestication of floral shapes. To explore the regulatory mechanisms of seasonal flowering, we have analyzed year-round gene expression patterns of C. japonica and C. azalea-a sister plant of continuous flowering that has been widely used for cross breeding. Through comparative analyses of gene co-expression networks and annual gene expression patterns, we show that annual expression rhythms of some important regulators of seasonal growth and development, including GIGANTEA and CONSTANS of the photoperiod pathway, have been disrupted in C. azalea. Furthermore, we reveal that the distinctive expression patterns of FLOWERING LOCUS T can be correlated with the seasonal activities of flowering and flushing. We demonstrate that the regulatory module involved in GIGANTEA, CONSTANS, and FLOWERING LOCUS T is central to achieve seasonality.

CONCLUSIONS

Through the genomic and comparative genomics characterizations of ornamental Camellia spp., we propose that duplication of chromosomal segments as well as the establishment of gene expression patterns has played a key role in the formation of ornamental traits (e.g., flower shape, flowering time). This work provides a valuable genomic platform for understanding the molecular basis of ornamental traits.

Genome-wide transcriptome profiling reveals molecular response pathways of Trichoderma harzianum in response to salt stress

Trichoderma harzianum exhibits a strong biological control effect on many important plant pathogens, such as Fusarium oxysporum, Botrytis cinerea, and Meloidogyne. However, its biocontrol effectiveness is weakened or reduced under salt stress. The aim of this study was to investigate the molecular response of T. harzianum to salt stress at the whole-genome level. Here, we present a 44.47 Mb near-complete genome assembly of the T. harzianum qt40003 strain for the first time, which was assembled de novo with 7.59 Gb Nanopore sequencing long reads (~170-fold) and 5.2 Gb Illumina short reads (~116-fold). The assembled qt40003 genome contains 12 contigs, with a contig N50 of 4.81 Mb, in which four of the 12 contigs were entirely reconstructed in a single chromosome from telomere to telomere. The qt40003 genome contains 4.27 Mb of repeat sequences and 12,238 protein-coding genes with a BUSCO completeness of 97.5%, indicating the high accuracy and completeness of our gene annotations. Genome-wide transcriptomic analysis was used to investigate gene expression changes related to salt stress in qt40003 at 0, 2% (T2), and 4% (T4) sodium chloride concentrations. A total of 2,937 and 3,527 differentially expressed genes (DEGs) were obtained under T2 and T4 conditions, respectively. GO enrichment analysis showed that the T2-treatment DEGs were highly enriched in detoxification (p < 0.001), while the T4 DEGs were mainly enriched in cell components, mostly in cellular detoxification, cell surface, and cell wall. KEGG metabolic pathway analysis showed that 91 and 173 DEGs were significantly enriched in the T2 and T4 treatments, respectively (p < 0.01), mainly in the glutathione metabolism pathway. We further experimentally analyzed the differentially expressed glutathione transferase genes in the glutathione metabolic pathway, most of which were downregulated (13/15). In addition, we screened 13 genes related to active oxygen clearance, including six upregulated and seven downregulated genes, alongside five fungal hydrophobic proteins, of which two genes were highly expressed. Our study provides high-quality genome information for the use of T. harzianum for biological control and offers significant insights into the molecular responses of T. harzianum under salt-stress conditions.

A chromosome-level genome assembly of Agave hybrid NO.11648 provides insights into the CAM photosynthesis

The subfamily Agavoideae comprises crassulacean acid metabolism (CAM), C3, and C4 plants with a young age of speciation and slower mutation accumulation, making it a model crop for studying CAM evolution. However, the genetic mechanism underlying CAM evolution remains unclear because of lacking genomic information. This study assembled the genome of Agave hybrid NO.11648, a constitutive CAM plant belonging to subfamily Agavoideae, at the chromosome level using data generated from high-throughput chromosome conformation capture, Nanopore, and Illumina techniques, resulting in 30 pseudo-chromosomes with a size of 4.87 Gb and scaffold N50 of 186.42 Mb. The genome annotation revealed 58 841 protein-coding genes and 76.91% repetitive sequences, with the dominant repetitive sequences being the I-type repeats (Copia and Gypsy accounting for 18.34% and 13.5% of the genome, respectively). Our findings also provide support for a whole genome duplication event in the lineage leading to A. hybrid, which occurred after its divergence from subfamily Asparagoideae. Moreover, we identified a gene duplication event in the phosphoenolpyruvate carboxylase kinase (PEPCK) gene family and revealed that three PEPCK genes (PEPCK3, PEPCK5, and PEPCK12) were involved in the CAM pathway. More importantly, we identified transcription factors enriched in the circadian rhythm, MAPK signaling, and plant hormone signal pathway that regulate the PEPCK3 expression by analysing the transcriptome and using yeast one-hybrid assays. Our results shed light on CAM evolution and offer an essential resource for the molecular breeding program of Agave spp.

Comparative genomics reveals the dynamic evolutionary history of cement protein genes of barnacles from intertidal to deep-sea hydrothermal vents

Thoracican barnacles are a diverse group of marine organisms for which the availability of genome assemblies is currently limited. In this study, we sequenced the genomes of two neolepadoid species (Ashinkailepas kermadecensis, Imbricaverruca yamaguchii) from hydrothermal vents, in addition to two intertidal species. Genome sizes ranged from 481 to 1054 Mb, with repetitive sequence contents of 21.2% to 50.7%. Concordance rates of orthologs and heterozygosity rates were between 82.4% and 91.7% and between 1.0% and 2.1%, respectively, indicating high genetic diversity and heterozygosity. Based on phylogenomic analyses, we revised the nomenclature of cement genes encoding cement proteins that are not homologous to any known proteins. The major cement gene, CP100A, was found in all thoracican species, including vent-associated neolepadoids, and was hypothesised to be essential for thoracican settlement. Duplicated genes, CP100B and CP100C, were found only in balanids, suggesting potential functional redundancy or acquisition of new functions associated with the calcareous base. An ancestor of CP52 genes was duplicated dynamically among lepadids, pollicipedids with multiple copies on a single scaffold, and balanids with multiple sequential repeats of the conserved regions, but no CP52 genes were found in neolepadoids, providing insights into cement gene evolution among thoracican lineages. This study enhances our understanding of the adhesion mechanisms of thoracicans in underwater environments. The newly sequenced genomes provide opportunities for studying their evolution and ecology, shedding light on their adaptation to diverse marine environments, and contributing to our knowledge of barnacle biology with valuable genomic resources for further studies in this field.

Genome Sequencing and Analysis of Nigrospora oryzae, a Rice Leaf Disease Fungus

Nigrospora oryzae is one of several fungal pathogens known to cause brown streaks, leaf spots, and latent infections in rice. In this study, the entire 42.09-Mb genome of N. oryzae was sequenced at a depth of 169× using the Oxford Nanopore Technologies platform. The draft genome sequence was comprised of 26 scaffolds, possessed an average GC content of 58.83%, and contained a total of 10,688 protein-coding genes. Analysis of the complete genome sequence revealed that CAZyme-encoding genes account for 6.11% of all identified genes and that numerous transcription factors (TFs) associated with diverse biological processes belong predominantly to Zn-clus (22.20%) and C2H2 (10.59%) fungal TF classes. In addition, genes encoding 126 transport proteins and 3307 pathogen-host interaction proteins were identified. Comparative analysis of the previously reported N. oryzae reference strain GZL1 genome and the genome of a representative strain ZQ1 obtained here revealed 9722 colinear genes. Collectively, these findings provide valuable insights into N. oryzae genetic mechanisms and phenotypic characteristics.

Integrated Genome Sequencing and Transcriptome Analysis Identifies Candidate Pathogenicity Genes from Ustilago crameri

Ustilago crameri is a pathogenic basidiomycete fungus that causes foxtail millet kernel smut (FMKS), a devastating grain disease in most foxtail-millet-growing regions of the world. Here, we report an assembled high-quality genome sequence of U. crameri strain SCZ-6 isolated from the diseased grains of foxtail millet in Changzhi, Shanxi Province, China. The genome size is 19.55 Mb, consisting of 73 contigs (N50 = 840,209 bp) with a G + C content of 54.09%, and encoding 6576 predicted genes and 6486 genes supported by RNA-seq. Evolutionarily, U. crameri lies close to the barley smut U. hordei, and an obvious co-linearity was observed between these two smut fungi. We annotated the genome of U. crameri strain SCZ-6 using databases, identifying 1827 pathogen-host interaction (PHI)-associated genes, 1324 genes encoding fungal virulence factors, 259 CAZy-related genes, 80 genes encoding transporters, and 206 putative cytochrome P450 genes; their expression profiles at different inoculation time points were also detected. Additionally, 70 candidate pathogen effectors were identified according to their expression patterns and predicted functions. In summary, our results provide important insights into the pathogenic mechanisms of the pathogenesis-related genes of U. crameri and a robust foundation for further investigation.

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-scale genome assembly of the 'Munstead' cultivar of Lavandula angustifolia

OBJECTIVES

Lavandula angustifolia (English lavender) is commercially important not only as an ornamental species but also as a major source of fragrances. To better understand the genomic basis of chemical diversity in lavender, we sequenced, assembled, and annotated the 'Munstead' cultivar of L. angustifolia.

DATA DESCRIPTION

A total of 80 Gb of Oxford Nanopore Technologies reads was used to assemble the 'Munstead' genome using the Canu genome assembler software. Following multiple rounds of error correction and scaffolding using Hi-C data, the final chromosome-scale assembly represents 795,075,733 bp across 25 chromosomes with an N50 scaffold length of 31,371,815 bp. Benchmarking Universal Single Copy Orthologs analysis revealed 98.0% complete orthologs, indicative of a high-quality assembly representative of genic space. Annotation of protein-coding sequences revealed 58,702 high-confidence genes encoding 88,528 gene models. Access to the 'Munstead' genome will permit comparative analyses within and among lavender accessions and provides a pivotal species for comparative analyses within Lamiaceae.

The draft genome of the Temminck's tragopan (Tragopan temminckii) with evolutionary implications

BACKGROUND

High-quality genome data of birds play a significant role in the systematic study of their origin and adaptive evolution. The Temminck's tragopan (Tragopan temminckii) (Galliformes, Phasianidae), a larger pheasant, is one of the most abundant and widely distributed species of the genus Tragopan, and was defined as class II of the list of national key protected wild animals in China. The absence of a sequenced genome has restricted previous evolutionary trait studies of this taxa.

RESULTS

The whole genome of the Temminck's tragopan was sequenced using Illumina and PacBio platform, and then de novo assembled and annotated. The genome size was 1.06 Gb, with a contig N50 of 4.17 Mb. A total of 117.22 Mb (11.00%) repeat sequences were identified. 16,414 genes were predicted using three methods, with 16,099 (98.08%) annotated as functional genes based on five databases. In addition, comparative genome analyses were conducted across 12 Galliformes species. The results indicated that T. temminckii was the first species to branch off from the clade containing Lophura nycthemera, Phasianus colchicus, Chrysolophus pictus, Syrmaticus mikado, Perdix hodgsoniae, and Meleagris gallopavo, with a corresponding divergence time of 31.43 million years ago (MYA). Expanded gene families associated with immune response and energy metabolism were identified. Genes and pathways associated with plumage color and feather development, immune response, and energy metabolism were found in the list of positively selected genes (PSGs).

CONCLUSIONS

A genome draft of the Temminck's tragopan was reported, genome feature and comparative genome analysis were described, and genes and pathways related to plumage color and feather development, immune response, and energy metabolism were identified. The genomic data of the Temminck's tragopan considerably contribute to the genome evolution and phylogeny of the genus Tragopan and the whole Galliformes species underlying ecological adaptation strategies.

Heat shock-inducible clonal analysis reveals the stepwise establishment of cell fate in the rice stem

The stem, consisting of nodes and internodes, is the shoot axis, which supports aboveground organs and connects them to roots. In contrast to other organs, developmental processes of the stem remain elusive, especially those initiating nodes and internodes. By introducing an intron into the Cre recombinase gene, we established a heat shock-inducible clonal analysis system in a single binary vector and applied it to the stem in the flag leaf phytomer of rice (Oryza sativa). With detailed characterizations of stem structure and development, we show that cell fate acquisition for each domain of the stem occurs stepwise. Cell fate for a single phytomer was established in the shoot apical meristem (SAM) by one plastochron before leaf initiation. Cells destined for the foot (nonelongating domain at the stem base) also started emerging before leaf initiation. Cell fate acquisition for the node began just before leaf initiation at the flank of the SAM, separating cell lineages for leaves and stems. Subsequently, cell fates for the axillary bud were established in early leaf primordia. Finally, cells committed to the internode emerged from, at most, a few cell tiers of the 12- to 25-cell stage stem epidermis. Thus, internode cell fate is established last during stem development. This study provides the groundwork to unveil underlying molecular mechanisms in stem development and a valuable tool for clonal analysis, which can be applied to various species.

Cold stress induces differential gene expression of retained homeologs in Camelina sativa cv Suneson

Camelina sativa (L.) Crantz, a member of the Brassicaceae, has potential as a biofuel feedstock which is attributable to the production of fatty acids in its seeds, its fast growth cycle, and low input requirements. While a genome assembly is available for camelina, it was generated from short sequence reads and is thus highly fragmented in nature. Using long read sequences, we generated a chromosome-scale, highly contiguous genome assembly (644,491,969 bp) for the spring biotype cultivar 'Suneson' with an N50 contig length of 12,031,512 bp and a scaffold N50 length of 32,184,682 bp. Annotation of protein-coding genes revealed 91,877 genes that encode 133,355 gene models. We identified a total of 4,467 genes that were significantly up-regulated under cold stress which were enriched in gene ontology terms associated with "response to cold" and "response to abiotic stress". Coexpression analyses revealed multiple coexpression modules that were enriched in genes differentially expressed following cold stress that had putative functions involved in stress adaptation, specifically within the plastid. With access to a highly contiguous genome assembly, comparative analyses with Arabidopsis thaliana revealed 23,625 A. thaliana genes syntenic with 45,453 Suneson genes. Of these, 24,960 Suneson genes were syntenic to 8,320 A. thaliana genes reflecting a 3 camelina homeolog to 1 Arabidopsis gene relationship and retention of all three homeologs. Some of the retained triplicated homeologs showed conserved gene expression patterns under control and cold-stressed conditions whereas other triplicated homeologs displayed diverged expression patterns revealing sub- and neo-functionalization of the homeologs at the transcription level. Access to the chromosome-scale assembly of Suneson will enable both basic and applied research efforts in the improvement of camelina as a sustainable biofuel feedstock.

A chromosome-level genome assembly of the Asian giant softshell turtle Pelochelys cantorii

The Asian giant softshell turtle Pelochelys cantorii is one of the largest aquatic turtles in China and has been designated a First Grade Protected Animal in China. To advance conservation research, a combination of Illumina short-read, PacBio long-read, and Hi-C scaffolding technologies was used to develop a high-quality chromosome-level genome assembly for P. cantorii. A total of 262.77 Gb of clean data were produced (121.6 × depth) and then the genome was assembled into 2.16 Gb with a contig N50 of 41.44 Mb and scaffold N50 length of 120.17 Mb, respectively. Moreover, about 99.98% assembly genome sequences were clustered and ordered onto 33 pseudochromosomes. Genome annotation revealed that 21,833 protein-coding genes were predicted, and 96.40% of them were annotated. This new chromosome-level assembly will be an enabling resource for genetic and genomic studies to support fundamental insight into P. cantorii biology.

Comprehensive Profiling of Alternative Splicing and Alternative Polyadenylation during Fruit Ripening in Watermelon (Citrullus lanatus)

Fruit ripening is a highly complicated process that is accompanied by the formation of fruit quality. In recent years, a series of studies have demonstrated post-transcriptional control play important roles in fruit ripening and fruit quality formation. Till now, the post-transcriptional mechanisms for watermelon fruit ripening have not been comprehensively studied. In this study, we conducted PacBio single-molecule long-read sequencing to identify genome-wide alternative splicing (AS), alternative polyadenylation (APA) and long non-coding RNAs (lncRNAs) in watermelon fruit. In total, 6,921,295 error-corrected and mapped full-length non-chimeric (FLNC) reads were obtained. Notably, more than 42,285 distinct splicing isoforms were derived from 5,891,183 intron-containing full-length FLNC reads, including a large number of AS events associated with fruit ripening. In addition, we characterized 21,506 polyadenylation sites from 11,611 genes, 8703 of which have APA sites. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that fructose and mannose metabolism, starch and sucrose metabolism and carotenoid biosynthesis were both enriched in genes undergoing AS and APA. These results suggest that post-transcriptional regulation might potentially have a key role in regulation of fruit ripening in watermelon. Taken together, our comprehensive PacBio long-read sequencing results offer a valuable resource for watermelon research, and provide new insights into the molecular mechanisms underlying the complex regulatory networks of watermelon fruit ripening.

Multi-omics analysis reveals the evolutionary origin of diterpenoid alkaloid biosynthesis pathways in Aconitum

Diterpenoid alkaloids (DAs) have been often utilized in clinical practice due to their analgesic and anti-inflammatory properties. Natural DAs are prevalent in the family Ranunculaceae, notably in the Aconitum genus. Nevertheless, the evolutionary origin of the biosynthesis pathway responsible for DA production remains unknown. In this study, we successfully assembled a high-quality, pseudochromosome-level genome of the DA-rich species Aconitum vilmorinianum (A. vilmorinianum) (5.76 Gb). An A. vilmorinianum-specific whole-genome duplication event was discovered using comparative genomic analysis, which may aid in the evolution of the DA biosynthesis pathway. We identified several genes involved in DA biosynthesis via integrated genomic, transcriptomic, and metabolomic analyses. These genes included enzymes encoding target ent-kaurene oxidases and aminotransferases, which facilitated the activation of diterpenes and insertion of nitrogen atoms into diterpene skeletons, thereby mediating the transformation of diterpenes into DAs. The divergence periods of these genes in A. vilmorinianum were further assessed, and it was shown that two major types of genes were involved in the establishment of the DA biosynthesis pathway. Our integrated analysis offers fresh insights into the evolutionary origin of DAs in A. vilmorinianum as well as suggestions for engineering the biosynthetic pathways to obtain desired DAs.

Unravelling the genomic features, phylogeny and genetic basis of tooth ontogenesis in Characiformes through analysis of four genomes

Characiformes is a diverse and evolutionarily significant order of freshwater fish encompassing over 2,300 species. Despite its diversity, our understanding of Characiformes' evolutionary relationships and adaptive mechanisms is limited due to insufficient genome sequences. In this study, we sequenced and assembled the genomes of four Characiformes species, three of which were chromosome-level assemblies. Our analyses revealed dynamic changes in gene family evolution, repeat sequences and variations in chromosomal collinearity within these genomes. With the assembled genomes, we were not only able to elucidate the evolutionary relationship of the four main orders in Otophysi but also indicated Characiformes as the paraphyletic group. Comparative genomic analysis with other available fish genomes shed light on the evolution of genes related to tooth development in Characiformes. Notably, variations in the copy number of secretory calcium-binding phosphoproteins (SCPP) genes were observed among different orders of Otophysi, indicating their potential contribution to the diversity of tooth types. Our study offers invaluable genome sequences and novel insights into Characiformes' evolution, paving the way for further genomic and evolutionary research in fish.

MkcDBGAS: a reference-free approach to identify comprehensive alternative splicing events in a transcriptome

Alternative splicing (AS) is an essential post-transcriptional mechanism that regulates many biological processes. However, identifying comprehensive types of AS events without guidance from a reference genome is still a challenge. Here, we proposed a novel method, MkcDBGAS, to identify all seven types of AS events using transcriptome alone, without a reference genome. MkcDBGAS, modeled by full-length transcripts of human and Arabidopsis thaliana, consists of three modules. In the first module, MkcDBGAS, for the first time, uses a colored de Bruijn graph with dynamic- and mixed- kmers to identify bubbles generated by AS with precision higher than 98.17% and detect AS types overlooked by other tools. In the second module, to further classify types of AS, MkcDBGAS added the motifs of exons to construct the feature matrix followed by the XGBoost-based classifier with the accuracy of classification greater than 93.40%, which outperformed other widely used machine learning models and the state-of-the-art methods. Highly scalable, MkcDBGAS performed well when applied to Iso-Seq data of Amborella and transcriptome of mouse. In the third module, MkcDBGAS provides the analysis of differential splicing across multiple biological conditions when RNA-sequencing data is available. MkcDBGAS is the first accurate and scalable method for detecting all seven types of AS events using the transcriptome alone, which will greatly empower the studies of AS in a wider field.

Genome assemblies of the simultaneously hermaphroditic flatworms Macrostomum cliftonense and Macrostomum hystrix

The free-living, simultaneously hermaphroditic flatworms of the genus Macrostomum are increasingly used as model systems in various contexts. In particular, Macrostomum lignano, the only species of this group with a published genome assembly, has emerged as a model for the study of regeneration, reproduction, and stem-cell function. However, challenges have emerged due to M. lignano being a hidden polyploid, having recently undergone whole-genome duplication and chromosome fusion events. This complex genome architecture presents a significant roadblock to the application of many modern genetic tools. Hence, additional genomic resources for this genus are needed. Here, we present such resources for Macrostomum cliftonense and Macrostomum hystrix, which represent the contrasting mating behaviors of reciprocal copulation and hypodermic insemination found in the genus. We use a combination of PacBio long-read sequencing and Illumina shot-gun sequencing, along with several RNA-Seq data sets, to assemble and annotate highly contiguous genomes for both species. The assemblies span ∼227 and ∼220 Mb and are represented by 399 and 42 contigs for M. cliftonense and M. hystrix, respectively. Furthermore, high BUSCO completeness (∼84-85%), low BUSCO duplication rates (8.3-6.2%), and low k-mer multiplicity indicate that these assemblies do not suffer from the same assembly ambiguities of the M. lignano genome assembly, which can be attributed to the complex karyology of this species. We also show that these resources, in combination with the prior resources from M. lignano, offer an excellent foundation for comparative genomic research in this group of organisms.

Chromosome-level assembly of Dictyophora rubrovolvata genome using third-generation DNA sequencing and Hi-C analysis

Dictyophora rubrovolvata, a rare edible mushroom with both nutritional and medicinal values, was regarded as the "queen of the mushroom" for its attractive appearance. Dictyophora rubrovolvata has been widely cultivated in China in recent years, and many researchers were focusing on its nutrition, culture condition, and artificial cultivation. Due to a lack of genomic information, research on bioactive substances, cross breeding, lignocellulose degradation, and molecular biology is limited. In this study, we report a chromosome-level reference genome of D. rubrovolvata using the PacBio single-molecule real-time-sequencing technique and high-throughput chromosome conformation capture (Hi-C) technologies. A total of 1.83 Gb circular consensus sequencing reads representing ∼983.34 coverage of the D. rubrovolvata genome were generated. The final genome was assembled into 136 contigs with a total length of 32.89 Mb. The scaffold and contig N50 length were 2.71 and 2.48 Mb, respectively. After chromosome-level scaffolding, 11 chromosomes with a total length of 28.24 Mb were constructed. Genome annotation further revealed that 9.86% of the genome was composed of repetitive sequences, and a total of 508 noncoding RNA (rRNA: 329, tRNA: 150, ncRNA: 29) were annotated. In addition, 9,725 protein-coding genes were predicted, among which 8,830 (90.79%) genes were predicted using homology or RNA-seq. Benchmarking Universal Single-Copy Orthologs results further revealed that there were 80.34% complete single-copy fungal orthologs. In this study, a total of 360 genes were annotated as belonging to the carbohydrate-active enzymes family. Further analysis also predicted 425 cytochromes P450 genes, which can be classified into 41 families. This highly accurate, chromosome-level reference genome of D. rubrovolvata will provide essential genomic information for understanding the molecular mechanism in its fruiting body formation during morphological development and facilitate the exploitation of medicinal compounds produced by this mushroom.

Genomic Resequencing Unravels the Genetic Basis of Domestication, Expansion, and Trait Improvement in Morus Atropurpurea

Mulberry is an economically important plant in the sericulture industry and traditional medicine. However, the genetic and evolutionary history of mulberry remains largely unknown. Here, this work presents the chromosome-level genome assembly of Morus atropurpurea (M. atropurpurea), originating from south China. Population genomic analysis using 425 mulberry accessions reveal that cultivated mulberry is classified into two species, M. atropurpurea and M. alba, which may have originated from two different mulberry progenitors and have independent and parallel domestication in north and south China, respectively. Extensive gene flow is revealed between different mulberry populations, contributing to genetic diversity in modern hybrid cultivars. This work also identifies the genetic architecture of the flowering time and leaf size. In addition, the genomic structure and evolution of sex-determining regions are identified. This study significantly advances the understanding of the genetic basis and domestication history of mulberry in the north and south, and provides valuable molecular markers of desirable traits for mulberry breeding.

Assembly of Genome and Resequencing Provide Insights into Genetic Differentiation between Parents of Hulong Hybrid Grouper (Epinephelus fuscoguttatus ♀ × E. lanceolatus ♂)

The Hulong hybrid grouper was bred from the brown-marbled grouper (Epinephelus fuscoguttatus) ♀ and the giant grouper (E. lanceolatus) ♂, combining the advantageous traits of both parents. Possessing an excellent performance, this hybrid's cultivation promotes the development of the grouper industry. Its male parent, the giant grouper, possesses the fastest growth and the largest body size among all coral-reef-dwelling fish. This species is not only an economically important species in marine aquaculture, but it is also an ideal male parent in the interspecific crossing of grouper species. In the present study, a high-quality chromosome-level genome of the giant grouper was constructed with a total length of 1.06 Gb, consisting of 24 chromosomes and 69 scaffolds. To analyze the genetic differences between the parents of the Hulong hybrid grouper, the structural variations (SVs) between both parental genomes were detected, and a total of 46,643 SVs were obtained. High-quality SNPs were identified from resequencing data. There were significant differences between the two genomes, and the average FST reached 0.685. A total of 234 highly differentiated regions were detected with an FST > 0.9. The protein-coding genes involved in SVs and highly differentiated regions were significantly enriched in metabolic pathways, including fatty metabolism, carbohydrate metabolism, amino acid metabolism and the TCA cycle. These genes may be related to the differences in feeding preferences and the ability to digest carbohydrates between the two grouper species under natural conditions. In addition, protein-coding genes related to the cell cycle and p53-signaling pathway were also detected. These genes may play important roles in the regulation of body size and growth performance. This research provides genomic resources for further breeding works and evolutionary analyses.

Genome Sequence Resource of a Colletotrichum graminicola Field Strain from China

The maize anthracnose stalk rot and leaf blight diseases caused by the fungal pathogen Colletotrichum graminicola is emerging as an important threat to corn production worldwide. In this work, we provide an improved genome assembly of a C. graminicola strain (TZ-3) by using the PacBio Sequel II and Illumina high-throughput sequencing technologies. The genome of TZ-3 consists of 36 contigs with a length of 59.3 Mb. After correction and evaluation with the Illumina sequencing data and BUSCO, this genome showed a high assembly quality and integrity. Gene annotation of this genome predicted 11,911 protein-coding genes, among which 983 secreted protein-coding genes and 332 effector genes were predicted. Compared with previous genomes of C. graminicola strains, TZ-3 genome is superior in nearly all parameters. The genome assembly and annotation will enhance our knowledge of the genetic makeup of the pathogen and molecular mechanisms underlying its pathogenicity and will provide valuable insights into genome variation across different regions. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

High-quality genome assembly of Verticillium dahliae VD991 allows for screening and validation of pathogenic genes

Verticillium dahliae (V. dahliae) is a notorious soil-borne pathogen causing Verticillium wilt in more than 400 dicotyledonous plants, including a wide range of economically important crops, such as cotton, tomato, lettuce, potato, and romaine lettuce, which can result in extensive economic losses. In the last decade, several studies have been conducted on the physiological and molecular mechanisms of plant resistance to V. dahliae. However, the lack of a complete genome sequence with a high-quality assembly and complete genomic annotations for V. dahliae has limited these studies. In this study, we produced a full genomic assembly for V. dahliae VD991 using Nanopore sequencing technology, consisting of 35.77 Mb across eight pseudochromosomes and with a GC content of 53.41%. Analysis of the genome completeness assessment (BUSCO alignment: 98.62%; Illumina reads alignment: 99.17%) indicated that our efforts resulted in a nearly complete and high-quality genomic assembly. We selected 25 species closely related to V. dahliae for evolutionary analysis, confirming the evolutionary relationship between V. dahliae and related species, and the identification of a possible whole genome duplication event in V. dahliae. The interaction between cotton and V. dahliae was investigated by transcriptome sequencing resulting in the identification of many genes and pathways associated with cotton disease resistance and V. dahliae pathogenesis. These results will provide new insights into the pathogenic mechanisms of V. dahliae and contribute to the cultivation of cotton varieties resistant to Verticillium wilt.

Comparative genome and metabolome analyses uncover the evolution and flavonoid biosynthesis between Apocynum venetum and Apocynum hendersonii

Apocynum species have great application prospects in textile and phytoremediation of saline soil, are rich in flavonoids, and possess medicinal significance. Here, we report the draft genomes of Apocynum venetum and Apocynum hendersonii, and elucidate their evolutionary relationship. The high synteny and collinearity between the two suggested that they have experienced the same WGD event. Comparative analysis revealed that flavone 3-hydroxylase (ApF3H) and differentially evolved flavonoid 3-O-glucosyl transferase (ApUFGT) genes are critical for determining natural variation in flavonoid biosynthesis between the species. Overexpression of ApF3H-1 enhanced the total flavonoid content and promoted the antioxidant capacity of transformed plants compared to the wild-type. ApUFGT5 and 6 explained the diversification of flavonoids or their derivatives. These data provide biochemical insight and knowledge on the genetic regulation of flavonoid biosynthesis, supporting the adoption of these genes in breeding programs aimed at the multipurpose utilization of the plants.

The genetic basis of the leafy seadragon's unique camouflage morphology and avenues for its efficient conservation derived from habitat modeling

The leafy seadragon certainly is among evolution's most "beautiful and wonderful" species aptly named for its extraordinary camouflage mimicking its coastal seaweed habitat. However, limited information is known about the genetic basis of its phenotypes and conspicuous camouflage. Here, we revealed genomic signatures of rapid evolution and positive selection in core genes related to its camouflage, which allowed us to predict population dynamics for this species. Comparative genomic analysis revealed that seadragons have the smallest olfactory repertoires among all ray-finned fishes, suggesting adaptations to the highly specialized habitat. Other positively selected and rapidly evolving genes that serve in bone development and coloration are highly expressed in the leaf-like appendages, supporting a recent adaptive shift in camouflage appendage formation. Knock-out of bmp6 results in dysplastic intermuscular bones with a significantly reduced number in zebrafish, implying its important function in bone formation. Global climate change-induced loss of seagrass beds now severely threatens the continued existence of this enigmatic species. The leafy seadragon has a historically small population size likely due to its specific habitat requirements that further exacerbate its vulnerability to climate change. Therefore, taking climate change-induced range shifts into account while developing future protection strategies.

Alternative splicing shapes the transcriptome complexity in blackgram [Vigna mungo (L.) Hepper]

Vigna mungo, a highly consumed crop in the pan-Asian countries, is vulnerable to several biotic and abiotic stresses. Understanding the post-transcriptional gene regulatory cascades, especially alternative splicing (AS), may underpin large-scale genetic improvements to develop stress-resilient varieties. Herein, a transcriptome based approach was undertaken to decipher the genome-wide AS landscape and splicing dynamics in order to establish the intricacies of their functional interactions in various tissues and stresses. RNA sequencing followed by high-throughput computational analyses identified 54,526 AS events involving 15,506 AS genes that generated 57,405 transcripts isoforms. Enrichment analysis revealed their involvement in diverse regulatory functions and demonstrated that transcription factors are splicing-intensive, splice variants of which are expressed differentially across tissues and environmental cues. Increased expression of a splicing regulator NHP2L1/SNU13 was found to co-occur with lower intron retention events. The host transcriptome is significantly impacted by differential isoform expression of 1172 and 765 AS genes that resulted in 1227 (46.8% up and 53.2% downregulated) and 831 (47.5% up and 52.5% downregulated) transcript isoforms under viral pathogenesis and Fe²⁺ stressed condition, respectively. However, genes experiencing AS operate differently from the differentially expressed genes, suggesting AS is a unique and independent mode of regulatory mechanism. Therefore, it can be inferred that AS mediates a crucial regulatory role across tissues and stressful situations and the results would provide an invaluable resource for future endeavours in V. mungo genomics.

OsCRY2 and OsFBO10 co-regulate photomorphogenesis and photoperiodic flowering in indica rice

Cryptochromes (CRYs) are a class of photoreceptors that perceive blue/ultraviolet-A light of the visible spectrum to mediate a vast number of physiological responses in bacteria, fungi, animals and plants. In the present study, we have characterized OsCRY2 in a photoperiod sensitive indica variety, Basmati 370, by generating and analyzing overexpression (OE) and knock-down (KD) transgenic lines. The OsCRY2^OE lines displayed dwarfism as shown in their reduced plant height and leaf length, attributed largely by an overall reduction in their cell size. The OsCRY2^OE lines flowered significantly earlier and showed shorter and broader seeds with an overall reduced seed weight. The OsCRY2^KD lines showed contrasting phenotypes, such as increased plant height and delayed flowering, however, decreased seed size and weight were also observed in the KD lines, along with reduced spikelet fertility and high seed shattering rate in mature panicles. Novel interactions were confirmed between OsCRY2 and members of ZEITLUPE family of blue/ultraviolet-A light photoreceptors, encoded by OsFBO8, OsFBO9 and OsFBO10 which are orthologous to ZEITLUPE (ZTL), LOV KELCH PROTEIN2 (LKP2) and FLAVIN BINDING, KELCH REPEAT F-BOX1 (FKF1), respectively, of Arabidopsis thaliana. Since FKF1 is known to play a role in regulating photoperiodic flowering, OsFBO10 was chosen for further studies. OsCRY2 and OsFBO10 interacted in the nucleus and cytoplasm of the cell and cross-regulated the expression of each other. They were also found to regulate the expression of several genes involved in photoperiodic flowering in rice. Both OsCRY2 and OsFBO10 played a positive role in photomorphogenic responses in different light conditions. The physical interaction of OsCRY2 with OsFBO10, their involvement in common physiological and developmental pathways and their cross-regulation of each other suggest that the two photoreceptors may regulate common developmental pathways in plants, either jointly or redundantly.

Genome-Wide Identification and Expression Analysis of the Cytochrome P450 Gene Family in Bemisia tabaci MED and Their Roles in the Insecticide Resistance

The whitefly, Bemisia tabaci MED (Hemiptera: Aleyrodidae), is an omnivorous agricultural pest, which causes huge economic losses to agriculture and is highly resistant to many pesticides. The overexpression of cytochrome P450 may play an important role in host adaptation and insecticide resistance in B. tabaci MED. Therefore, the present study systematically analyzed the cytochrome P450 gene family at the genome-wide level to understand its function in B. tabaci MED. Our analysis identified 58 cytochrome P450 genes in B. tabaci MED, among which 24 were novel. Phylogenetic analysis revealed broad functional and species-specific diversification in B. tabaci MED P450, suggesting the role of multiple P450 genes in detoxifying. Reverse transcription-real time quantitative PCR (RT-qPCR) showed that CYP4CS2, CYP4CS5, CYP4CS6, CYP4CS8, CYP6DW4, CYP6DW5, CYP6DW6, CYP6DZ8, and CYP6EN1 genes increased significantly after two days of exposure to imidacloprid. Interestingly, all nine genes belonged to the CYP4 and CYP6 families. A decrease in the expression of five genes (CYP6DW4, CYP6DW5, CYP6DW6, CYP6DZ8, and CYP4CS6) via RNA interference (RNAi) resulted in a significant increase in the mortalities of whiteflies when exposed to imidacloprid. These results indicate that the overexpression of the P450 genes may play an essential role in imidacloprid tolerance of B. tabaci MED. Thus, the present study provides basic information on P450 genes in B. tabaci MED, which will further help elucidate the insecticide resistance mechanism in the agricultural pest whitefly.

Intraspecific Comparative Analysis Reveals Genomic Variation of Didymella arachidicola and Pathogenicity Factors Potentially Related to Lesion Phenotype

Didymella arachidicola is one of the most important fungal pathogens, causing foliar disease and leading to severe yield losses of peanuts (Arachis hypogaea L.) in China. Two main lesion phenotypes of peanut web blotch have been identified as reticulation type (R type) and blotch type (B type). As no satisfactory reference genome is available, the genomic variations and pathogenicity factors of D. arachidicola remain to be revealed. In the present study, we collected 41 D. arachidicola isolates from 26 geographic locations across China (33 for R type and 8 for B type). The chromosome-scale genome of the most virulent isolate (YY187) was assembled as a reference using PacBio and Hi-C technologies. In addition, we re-sequenced 40 isolates from different sampling sites. Genome-wide alignments showed high similarity among the genomic sequences from the 40 isolates, with an average mapping rate of 97.38%. An average of 3242 SNPs and 315 InDels were identified in the genomic variation analysis, which revealed an intraspecific polymorphism in D. arachidicola. The comparative analysis of the most and least virulent isolates generated an integrated gene set containing 512 differential genes. Moreover, 225 genes individually or simultaneously harbored hits in CAZy-base, PHI-base, DFVF, etc. Compared with the R type reference, the differential gene sets from all B type isolates identified 13 shared genes potentially related to lesion phenotype. Our results reveal the intraspecific genomic variation of D. arachidicola isolates and pathogenicity factors potentially related to different lesion phenotypes. This work sets a genomic foundation for understanding the mechanisms behind genomic diversity driving different pathogenic phenotypes of D. arachidicola.