De novo assemblies of Luffa acutangula and Luffa cylindrica genomes reveal an expansion associated with substantial accumulation of transposable elements

Genome Assembly and Structural Variation Analysis of Luffa acutangula Provide Insights on Flowering Time and Ridge Development

Luffa spp. is an important worldwide cultivated vegetable and medicinal plant from the Cucurbitaceae family. In this study, we report a high-quality chromosome-level genome of the high-generation inbred line SG261 of Luffa acutangula. The genomic sequence was determined by PacBio long reads, Hi-C sequencing reads, and 10× Genomics sequencing, with an assembly size of 739.82 Mb, contig N50 of 18.38 Mb, and scaffold N50 of 56.08 Mb. The genome of L. acutangula SG261 was predicted to contain 27,312 protein-coding genes and 72.56% repetitive sequences, of which long terminal repeats (LTRs) were an important form of repetitive sequences, accounting for 67.84% of the genome. Phylogenetic analysis reveals that L. acutangula evolved later than Luffa cylindrica, and Luffa is closely related to Momodica charantia. Comparing the genome of L. acutangula SG261 and L. cylindrica with PacBio data, 67,128 high-quality structural variations (SVs) and 55,978 presence-absence variations (PAVs) were identified in SG261, resulting in 2424 and 1094 genes with variation in the CDS region, respectively, and there are 287 identical genes affected by two different structural variation analyses. In addition, we found that the transcription factor FY (FLOWERING LOCUS Y) families had a large expansion in L. acutangula SG261 (flowering in the morning) compared to L. cylindrica (flowering in the afternoon), which may result in the early flowering time in L. acutangula SG261. This study provides valuable reference for the breeding of and pan-genome research into Luffa species.

Analyzing genetic diversity in luffa and developing a Fusarium wilt-susceptible linked SNP marker through a single plant genome-wide association (sp-GWAS) study

BACKGROUND

Luffa (Luffa spp.) is an economically important crop of the Cucurbitaceae family, commonly known as sponge gourd or vegetable gourd. It is an annual cross-pollinated crop primarily found in the subtropical and tropical regions of Asia, Australia, Africa, and the Americas. Luffa serves not only as a vegetable but also exhibits medicinal properties, including anti-inflammatory, antidiabetic, and anticancer effects. Moreover, the fiber derived from luffa finds extensive applications in various fields such as biotechnology and construction. However, luffa Fusarium wilt poses a severe threat to its production, and existing control methods have proven ineffective in terms of cost-effectiveness and environmental considerations. Therefore, there is an urgent need to develop luffa varieties resistant to Fusarium wilt. Single-plant GWAS (sp-GWAS) has been demonstrated as a promising tool for the rapid and efficient identification of quantitative trait loci (QTLs) associated with target traits, as well as closely linked molecular markers.

RESULTS

In this study, a collection of 97 individuals from 73 luffa accessions including two major luffa species underwent single-plant GWAS to investigate luffa Fusarium wilt resistance. Utilizing the double digest restriction site associated DNA (ddRAD) method, a total of 8,919 high-quality single nucleotide polymorphisms (SNPs) were identified. The analysis revealed the potential for Fusarium wilt resistance in accessions from both luffa species. There are 6 QTLs identified from 3 traits, including the area under the disease progress curve (AUDPC), a putative disease-resistant QTL, was identified on the second chromosome of luffa. Within the region of linkage disequilibrium, a candidate gene homologous to LOC111009722, which encodes peroxidase 40 and is associated with disease resistance in Cucumis melo, was identified. Furthermore, to validate the applicability of the marker associated with resistance from sp-GWAS, an additional set of 21 individual luffa plants were tested, exhibiting 93.75% accuracy in detecting susceptible of luffa species L. aegyptiaca Mill.

CONCLUSION

In summary, these findings give a hint of genome position that may contribute to luffa wild resistance to Fusarium and can be utilized in the future luffa wilt resistant breeding programs aimed at developing wilt-resistant varieties by using the susceptible-linked SNP marker.

Establishment of virus-induced gene silencing (VIGS) system in Luffa acutangula using Phytoene desaturase (PDS) and tendril synthesis related gene (TEN)

BACKGROUND

Virus-induced gene silencing (VIGS) is a reverse genetics technology that can efficiently and rapidly identify plant gene functions. Although a variety of VIGS vectors have been successfully used in plants, only a few reports on VIGS technology in Luffa exist.

RESULTS

In the present study, a new cucumber green mottle mosaic virus (CGMMV)-based VIGS vector, pV190, was applied to establish the CGMMV-VIGS to investigate the feasibility of the silencing system for Luffa. Phytoene desaturase (PDS) gene was initially selected as a VIGS marker gene to construct a recombinant vector. Plants infected with Agrobacterium harboring pV190-PDS successfully induced effective silencing in Luffa, and an effective gene silencing phenotype with obvious photobleaching was observed. To further validate the efficiency, we selected TEN for gene-silencing, which encodes a CYC/TB1-like transcription factor and is involved in tendril development. Luffa plants inoculated with the pV190-TEN exhibited shorter tendril length and nodal positions where tendrils appear are higher compared to those of non-inoculated plants. RT-qPCR showed that the expression levels of PDS and TEN were significantly reduced in the CGMMV-VIGS plants. Moreover, we evaluated the CGMMV-VIGS efficiency in three cucurbits, including cucumber, ridge gourd, and bottle gourd.

CONCLUSION

We successfully established a CGMMV-based VIGS system on ridge gourd and used marker genes to identify the feasibility of the silencing system in Luffa leaves and stems.

Inheritance and molecular mapping of solitary/cluster fruit-bearing habit in Luffa

Fruiting behaviour and sex form are important goals for Luffa breeders and this study aimed to shed light upon inheritance patterns for both these traits. The hermaphrodite form of Luffa acutangula (known as Satputia) is an underutilized vegetable with a unique clustered fruiting habit. Its desirable traits, such as plant architecture, earliness, as well as contrasting traits like unique clustered fruiting, bisexual flower, and crossability with Luffa acutangula (monoecious ridge gourd with solitary fruits), make it a potential source for trait improvement and mapping of desirable traits in Luffa. In the present study, we have elucidated the inheritance pattern of fruiting behaviour in Luffa using F₂ mapping population generated from a cross between Pusa Nutan (Luffa acutangula, monoecious, solitary fruiting) × DSat-116 (Luffa acutangula, hermaphrodite, cluster fruiting). In F₂ generation, the observed distribution of plant phenotypes fitted in the expected ratio of 3:1 (solitary vs cluster) for fruit-bearing habit. This is the first report of monogenic recessive control for cluster fruit-bearing habit in Luffa. Herein, we designate for the first time the gene symbol cl for cluster fruit bearing in Luffa. Linkage analysis revealed that SRAP marker ME10 EM4-280 was linked to the fruiting trait at the distance of 4.6 cM from the Cl locus. In addition, the inheritance pattern of hermaphrodite sex form in Luffa was also studied in the F₂ population of Pusa Nutan × DSat-116 that segregated into 9:3:3:1 ratio (monoecious:andromonoecious:gynoecious:hermaphrodite), suggesting a digenic recessive control of hermaphrodite sex form in Luffa, which was further confirmed by the test cross. The inheritance and identification of molecular marker for cluster fruiting trait provides a basis for breeding in Luffa species.

Lateral Root Initiation in Cucumber (Cucumis sativus): What Does the Expression Pattern of Rapid Alkalinization Factor 34 (RALF34) Tell Us?

In Arabidopsis, the small signaling peptide (peptide hormone) RALF34 is involved in the gene regulatory network of lateral root initiation. In this study, we aimed to understand the nature of the signals induced by RALF34 in the non-model plant cucumber (Cucumis sativus), where lateral root primordia are induced in the apical meristem of the parental root. The RALF family members of cucumber were identified using phylogenetic analysis. The sequence of events involved in the initiation and development of lateral root primordia in cucumber was examined in detail. To elucidate the role of the small signaling peptide CsRALF34 and its receptor CsTHESEUS1 in the initial stages of lateral root formation in the parental root meristem in cucumber, we studied the expression patterns of both genes, as well as the localization and transport of the CsRALF34 peptide. CsRALF34 is expressed in all plant organs. CsRALF34 seems to differ from AtRALF34 in that its expression is not regulated by auxin. The expression of AtRALF34, as well as CsRALF34, is regulated in part by ethylene. CsTHESEUS1 is expressed constitutively in cucumber root tissues. Our data suggest that CsRALF34 acts in a non-cell-autonomous manner and is not involved in lateral root initiation in cucumber.

Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus)

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management. Here, we generated 55 Gb of short Illumina reads (92× coverage) and 73 Gb of long Nanopore reads (122×) to produce the first genome assembly of the marine teleost tarakihi [Nemadactylus macropterus (Forster, 1801)], a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq reads were obtained to assist in gene annotation. The final genome assembly was 568 Mb long with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii Benchmarking Universal Single-Copy Orthologs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared with the same values reported for other fishes. Iso-Seq analysis recovered 91,313 unique transcripts from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly and the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics and comparative eco-evolutionary studies in teleosts and related organisms.

QTL mapping reveals candidate genes for main agronomic traits in Luffa based on a high-resolution genetic map

Luffa is an important medicinal and edible vegetable crop of Cucurbitaceae. Strong heterosis effects and strikingly complementary characteristics were found between the two domesticated Luffa cultivars, Luffa acutangula and Luffa cylindrica. To explore the genetic basis underlying their important agronomic traits, we constructed the first interspecific high-density genetic linkage map using a BC₁ population of 110 lines derived from a cross between S1174 (Luffa acutangula) and P93075 (Luffa cylindrica). The map spanned a total of 2246.74 cM with an average distance of 0.48 cM between adjacent markers. Thereafter, a large-scale field-based quantitative trait loci (QTLs) mapping was conducted for 25 important agronomic traits and 40 significant genetic loci distributed across 11 chromosomes were detected. Notably, a vital QTL (qID2) located on chromosome 9 with a minimum distance of 23 kb was identified to be responsible for the internode diameter and explained 11% of the phenotypic variation. Lac09g006860 (LacCRWN3), encoding a nuclear lamina protein involved in the control of nuclear morphology, was the only gene harbored in qID2. Sequence alignment showed completely different promoter sequences between the two parental alleles of LacCRWN3 except for some nonsynonymous single nucleotide polymorphisms (SNPs) in exons, and the expression level in thick-stem P93075 was distinctively higher than that in thin-stem S1174. According to the natural variation analysis of a population of 183 inbred lines, two main haplotypes were found for LacCRWN3: the P93075-like and S1174-like, with the former haplotype lines exhibiting significantly thicker internode diameters than those of the latter haplotype lines. It showed that LacCRWN3, as the only CRWN3 gene in Cucurbitaceae, was the most likely candidate gene regulating the internode diameter of Luffa. Our findings will be beneficial for deciphering the molecular mechanism of key phenotypic traits and promoting maker-assisted breeding in Luffa.

CuGenDBv2: an updated database for cucurbit genomics

The Cucurbitaceae (cucurbit) family consists of about 1,000 species in 95 genera, including many economically important and popular fruit and vegetable crops. During the past several years, reference genomes have been generated for >20 cucurbit species, and variome and transcriptome profiling data have been rapidly accumulated for cucurbits. To efficiently mine, analyze and disseminate these large-scale datasets, we have developed an updated version of Cucurbit Genomics Database. The updated database, CuGenDBv2 (http://cucurbitgenomics.org/v2), currently hosts 34 reference genomes from 27 cucurbit species/subspecies belonging to 10 different genera. Protein-coding genes from these genomes have been comprehensively annotated by comparing their protein sequences to various public protein and domain databases. A novel 'Genotype' module has been implemented to facilitate mining and analysis of the functionally annotated variome data including SNPs and small indels from large-scale genome sequencing projects. An updated 'Expression' module has been developed to provide a comprehensive gene expression atlas for cucurbits. Furthermore, synteny blocks between any two and within each of the 34 genomes, representing a total of 595 pair-wise genome comparisons, have been identified and can be explored and visualized in the database.

Assessment of the Genetic Diversity and Population Structure of Rhizophora apiculata Blume (Rhizophoraceae) in Thailand

Rhizophora apiculata is one of the most widespread and economically important mangrove trees in the Indo-West Pacific region. Knowledge of the genetic variation of R. apiculata in Thailand is limited. Here, we generated a whole-genome sequence of R. apiculata using the 10× Genomics technology. R. apiculata genome assembly was 230.47 Mb. Based on its genome, 2640 loci of high-quality biallelic SNPs were identified from 82 R. apiculata accessions collected from 17 natural mangrove forests in Thailand to assess the genetic diversity and population structure among them. A moderate level of genetic diversity of R. apiculata was observed. The average observed heterozygosity (Ho = 0.48) was higher than the average expected heterozygosity (He = 0.36). Two subpopulations were observed and confirmed from three approaches: population structure, PCA, and phylogenetic analyses. They corresponded to the Gulf of Thailand and the Andaman Sea separated by the Malay Peninsula. AMOVA analyses indicated that genetic variation was attributable to 76.22% within populations and 23.78% among populations. A high level of genetic differentiation between the two subpopulations (FST = 0.24, p < 0.001) was observed. This study evaluated the genetic diversity and population structure of R. apiculata, providing useful information for sustainable mangrove management in Thailand.

A pH- and redox-stimulated responsive hollow mesoporous silica for triggered delivery of fungicides to control downy mildew of Luffa cylindrica

BACKGROUND

Downy mildew, a devastating disease of cucurbitaceous crops caused by Pseudoperonospora cubensis. Although a variety of fungicides are used to control downy mildew, choosing an effective product can be challenging. Environmental stimulus-responsive pesticide delivery systems have great potential to improve the effectiveness of disease and pest control and reduce the impact on environmentally beneficial organisms.

RESULTS

In this work, a disulfide bond (SS)-modified and chitosan oligosaccharide (COS)-capped hollow mesoporous silica (HMS) pesticide delivery system was synthesized using a hard template method for the control of downy mildew in cucurbit crops. The synthesized nanoparticles were loaded with dimethomorph (DMM), denoted as DMM@HMS-SS-COS, and the developmental toxicity of these nanoparticles to zebrafish embryos were evaluated. The results showed that the prepared DMM@HMS-SS-COS exhibited excellent dual response properties to pH and glutathione (GSH), with an encapsulation rate of up to 24.36%. DMM@HMS-SS-COS has good ultraviolet (UV) radiation stability and adhesion properties. Compared with dimethomorph suspension concentrate (SC), DMM@HMS-SS-COS was more effective against downy mildew for up to 21 days. Toxicity tests showed that DMM@HMS-SS-COS significantly reduced the effect of DMM on the hatching rate and survival rate of zebrafish embryos.

CONCLUSIONS

This work not only demonstrates that DMM@HMS-SS-COS could be used as a nanodelivery system for intelligent control of downy mildew but also emphasizes the necessity of increasing the acute toxicity of nanoformulations to non-target organisms in environmental risk assessment. © 2022 Society of Chemical Industry.

A de novo reference assembly of the yellow mangrove Ceriops zippeliana genome

Mangroves are of great ecological and economical importance, providing shelters for a wide range of species and nursery habitats for commercially important marine species. Ceriops zippeliana (yellow mangrove) belongs to Rhizophoraceae family and is commonly distributed in the tropical and subtropical coastal communities. In this study, we present a high-quality assembly of the C. zippeliana genome. We constructed an initial draft assembly of 240,139,412 bases with an N50 contig length of 564,761 bases using the 10x Genomics linked-read technology. This assembly was further scaffolded with RagTag using a chromosome-scale assembly of a closely related Ceriops species as a reference. The final assembly contained 243,228,612 bases with an N50 scaffold length of 10,559,178 Mb. The size of the final assembly was close to those estimated using DNA flow cytometry (248 Mb) and the k-mer distribution analysis (246 Mb). We predicted a total of 23,474 gene models and 21,724 protein-coding genes in the C. zippeliana genome, of which 16,002 were assigned gene ontology terms. We recovered 97.1% of the highly conserved orthologs based on the Benchmarking Universal Single-Copy Orthologs analysis. The phylogenetic analysis based on single-copy orthologous genes illustrated that C. zippeliana and Ceriops tagal diverged approximately 10.2 million years ago (MYA), and their last common ancestor and Kandelia obovata diverged approximately 29.9 MYA. The high-quality assembly of C. zippeliana presented in this work provides a useful genomic resource for studying mangroves' unique adaptations to stressful intertidal habitats and for developing sustainable mangrove forest restoration and conservation programs.

De novo reference assembly of the upriver orange mangrove (Bruguiera sexangula) genome

Upriver orange mangrove (Bruguiera sexangula) is a member of the most mangrove-rich taxon (Rhizophoraceae family) and is commonly distributed in the intertidal zones in tropical and subtropical latitudes. In this study, we employed the 10× Genomics linked-read technology to obtain a preliminary de novo assembly of the B. sexangula genome, which was further scaffolded to a pseudomolecule level using the Bruguiera parviflora genome as a reference. The final assembly of the B. sexangula genome contained 260 Mb with an N50 scaffold length of 11,020,310 bases. The assembly comprised 18 pseudomolecules (corresponding to the haploid chromosome number in B. sexangula), covering 204,645,832 bases or 78.6% of the 260-Mb assembly. We predicted a total of 23,978 protein-coding sequences, 17,598 of which were associated with gene ontology terms. Our gene prediction recovered 96.6% of the highly conserved orthologs based on the Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis. The chromosome-level assembly presented in this work provides a valuable genetic resource to help strengthen our understanding of mangroves’ physiological and morphological adaptations to the intertidal zones.

Cucurbitaceae genome evolution, gene function and molecular breeding

The Cucurbitaceae is one of the most genetically diverse plant families in the world. Many of them are important vegetables or medicinal plants and are widely distributed worldwide. The rapid development of sequencing technologies and bioinformatic algorithms has enabled the generation of genome sequences of numerous important Cucurbitaceae species. This has greatly facilitated research on gene identification, genome evolution, genetic variation and molecular breeding of cucurbit crops. So far, genome sequences of 18 different cucurbit species belonging to tribes Benincaseae, Cucurbiteae, Sicyoeae, Momordiceae and Siraitieae have been deciphered. This review summarizes the genome sequence information, evolutionary relationship, and functional genes associated with important agronomic traits (e.g., fruit quality). The progress of molecular breeding in cucurbit crops and prospects for future applications of Cucurbitaceae genome information are also discussed.

Genetic diversity and population structure of ridge gourd (Luffa acutangula) accessions in a Thailand collection using SNP markers

This study explored a germplasm collection consisting of 112 Luffa acutangula (ridge gourd) accessions, mainly from Thailand. A total of 2834 SNPs were used to establish population structure and underlying genetic diversity while exploring the fruit characteristics together with genetic information which would help in the selection of parental lines for a breeding program. The study found that the average polymorphism information content value of 0.288 which indicates a moderate genetic diversity for this L. acutangula germplasm. STRUCTURE analysis (ΔK at K = 6) allowed us to group the accessions into six subpopulations that corresponded well with the unrooted phylogenetic tree and principal coordinate analyses. When plotted, the STRUCTURE bars to the area of collection, we observed an admixed genotype from surrounding accessions and a geneflow confirmed by the value of F_ST = 0.137. AMOVA based on STRUCTURE clustering showed a low 12.83% variation between subpopulations that correspond well with the negative inbreeding coefficient value (F_IS =  − 0.092) and low total fixation index (F_IT = 0.057). There were distinguishing fruit shapes and length characteristics in specific accessions for each subpopulation. The genetic diversity and different fruit shapes in the L. acutangula germplasm could benefit the ridge gourd breeding programs to meet the demands and needs of consumers, farmers, and vegetable exporters such as increasing the yield of fruit by the fruit width but not by the fruit length to solve the problem of fruit breakage during exportation.

Review on the Development and Applications of Medicinal Plant Genomes

With the development of sequencing technology, the research on medicinal plants is no longer limited to the aspects of chemistry, pharmacology, and pharmacodynamics, but reveals them from the genetic level. As the price of next-generation sequencing technology becomes affordable, and the long-read sequencing technology is established, the medicinal plant genomes with large sizes have been sequenced and assembled more easily. Although the review of plant genomes has been reported several times, there is no review giving a systematic and comprehensive introduction about the development and application of medicinal plant genomes that have been reported until now. Here, we provide a historical perspective on the current situation of genomes in medicinal plant biology, highlight the use of the rapidly developing sequencing technologies, and conduct a comprehensive summary on how the genomes apply to solve the practical problems in medicinal plants, like genomics-assisted herb breeding, evolution history revelation, herbal synthetic biology study, and geoherbal research, which are important for effective utilization, rational use and sustainable protection of medicinal plants.