Genetic Linkage Map Construction and QTL Analysis of Two Interspecific Reproductive Isolation Traits in Sponge Gourd

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.

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.

De novo transcriptome assembly, gene annotation, and EST-SSR marker development of an important medicinal and edible crop, Amomum tsaoko (Zingiberaceae)

BACKGROUND

Amomum tsaoko is a medicinal and food dual-use crop that belongs to the Zingiberaceae family. However, the lack of transcriptomic and genomic information has limited the understanding of the genetic basis of this species. Here, we performed transcriptome sequencing of samples from different A. tsaoko tissues, and identified and characterized the expressed sequence tag-simple sequence repeat (EST-SSR) markers.

RESULTS

A total of 58,278,226 high-quality clean reads were obtained and de novo assembled to generate 146,911 unigenes with an N50 length of 2002 bp. A total of 128,174 unigenes were successfully annotated by searching seven protein databases, and 496 unigenes were identified as annotated as putative terpenoid biosynthesis-related genes. Furthermore, a total of 55,590 EST-SSR loci were detected, and 42,333 primer pairs were successfully designed. We randomly selected 80 primer pairs to validate their polymorphism in A. tsaoko; 18 of these primer pairs produced distinct, clear, and reproducible polymorphisms. A total of 98 bands and 96 polymorphic bands were amplified by 18 pairs of EST-SSR primers for the 72 A. tsaoko accessions. The Shannon's information index (I) ranged from 0.477 (AM208) to 1.701 (AM242) with an average of 1.183, and the polymorphism information content (PIC) ranged from 0.223 (AM208) to 0.779 (AM247) with an average of 0.580, indicating that these markers had a high level of polymorphism. Analysis of molecular variance (AMOVA) indicated relatively low genetic differentiation among the six A. tsaoko populations. Cross-species amplification showed that 14 of the 18 EST-SSR primer pairs have transferability between 11 Zingiberaceae species.

CONCLUSIONS

Our study is the first to provide transcriptome data of this important medicinal and edible crop, and these newly developed EST-SSR markers are a very efficient tool for germplasm evaluation, genetic diversity, and molecular marker-assisted selection in A. tsaoko.

Identification of suitable reference genes for quantitative reverse transcription PCR in Luffa (Luffa cylindrica)

Reverse transcription real-time quantitative PCR is widely used to quantify gene expression. Reference genes are usually used as internal controls to measure the target gene expression level. To date, there is no consensus on the use of systematically validated reference genes in different tissues of Luffa. This study evaluated the expression stability of 11 candidate reference genes in different tissues using five algorithms (BestKeeper, comparative delta-Ct method, GeNorm, NormFinder, and RefFinder). Protein phosphatase 2A was the most stable gene, while alpha Tubulin was the least stable. The relative expression of ethylene-related genes in different tissues was also analyzed to reveal their role in sex determination. This study provides the basis for using suitable reference genes to evaluate targeted gene expression.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01182-8.

Metabolome and Transcriptome Analyses of Cucurbitacin Biosynthesis in Luffa (Luffa acutangula)

Cucurbitacins are extremely bitter compounds mainly present in Cucurbitaceae, where Luffa belongs. However, there is no comprehensive analysis of cucurbitacin biosynthesis in Luffa fruit. Therefore, this study analyzed bitter (WM709) and non-bitter (S1174) genotypes of Luffa to reveal the underlying mechanism of cucurbitacin biosynthesis by integrating metabolome and transcriptome analyses. A total of 422 metabolites were detected, including vitamins, essential amino acids, antioxidants, and antitumor substances. Of these, 131 metabolites showed significant differences between bitter (WM709) and non-bitter (S1174) Luffa fruits. The levels of isocucurbitacin B, cucurbitacin D, 23,24-dihydro cucurbitacin E, cucurbitacin F were significantly higher in bitter than in non-bitter Luffa. Transcriptome analysis showed that Bi, cytochromes P450s (CYP450s), and acyltransferase (ACT) of the cucurbitacin biosynthesis pathway, were significantly up-regulated. Moreover, drought stress and abscisic acid (ABA) activated genes of the cucurbitacin biosynthesis pathway. Furthermore, dual-luciferase reporter and yeast one-hybrid assays demonstrated that ABA-response element binding factor 1 (AREB1) binds to the Bi promoter to activate Bi expression. Comparative analysis of the Luffa and cucumber genomes showed that Bi, CYP450s, and ACT are located in the conserved syntenic loci, and formed a cucurbitacin biosynthesis cluster. This study provides important insights into major genes and metabolites of the cucurbitacin biosynthetic pathway, deepening the understanding of regulatory mechanisms of cucurbitacin biosynthesis in Luffa.

A mutation in LacDWARF1 results in a GA-deficient dwarf phenotype in sponge gourd (Luffa acutangula)

KEY MESSAGE

A dwarfism gene LacDWARF1 was mapped by combined BSA-Seq and comparative genomics analyses to a 65.4 kb physical genomic region on chromosome 05. Dwarf architecture is one of the most important traits utilized in Cucurbitaceae breeding because it saves labor and increases the harvest index. To our knowledge, there has been no prior research about dwarfism in the sponge gourd. This study reports the first dwarf mutant WJ209 with a decrease in cell size and internodes. A genetic analysis revealed that the mutant phenotype was controlled by a single recessive gene, which is designated Lacdwarf1 (Lacd1). Combined with bulked segregate analysis and next-generation sequencing, we quickly mapped a 65.4 kb region on chromosome 5 using F2 segregation population with InDel and SNP polymorphism markers. Gene annotation revealed that Lac05g019500 encodes a gibberellin 3β-hydroxylase (GA3ox) that functions as the most likely candidate gene for Lacd1. DNA sequence analysis showed that there is an approximately 4 kb insertion in the first intron of Lac05g019500 in WJ209. Lac05g019500 is transcribed incorrectly in the dwarf mutant owing to the presence of the insertion. Moreover, the bioactive GAs decreased significantly in WJ209, and the dwarf phenotype could be restored by exogenous GA3 treatment, indicating that WJ209 is a GA-deficient mutant. All these results support the conclusion that Lac05g019500 is the Lacd1 gene. In addition, RNA-Seq revealed that many genes, including those related to plant hormones, cellular process, cell wall, membrane and response to stress, were significantly altered in WJ209 compared with the wild type. This study will aid in the use of molecular marker-assisted breeding in the dwarf sponge gourd.

Progress on Genetic Improvement and Analysis of Bottle Gourd [Lagenaria siceraria (Molina) Standl.] for Agronomic Traits, Nutrient Compositions, and Stress Tolerance: A Review

Bottle gourd [Lagenaria siceraria (Molina) Standl.] is an important multi-purpose cucurbit crop grown for its leaf, fruit, and seed. It is widely cultivated and used for human consumption in sub-Saharan Africa (SSA) providing vital human nutrition and serving as food security crop. There is wide genetic variation among bottle gourd genetic resources in Africa for diverse qualitative and quantitative attributes for effective variety design, product development, and marketing. However, the crop is under- researched and -utilized, and improved varieties are yet to be developed and commercialized in the region. Therefore, the objective of this review is to provide the progress on bottle gourd genetic improvement and genetic analysis targeting agronomic and horticultural attributes, nutritional composition, biotic, and abiotic stress tolerance to guide current and future cultivar development, germplasm access, and conservation in SSA. The first section of the paper presents progress on breeding of bottle gourd for horticultural traits, agronomic performance, nutritional and anti-nutritional composition, and biotic and abiotic stress tolerance. This is followed by important highlights on key genetic resources of cultivated and wild bottle gourd for demand driven breeding. Lastly, the review summaries advances in bottle gourd genomics, genetic engineering and genome editing. Information presented in this paper should aid bottle gourd breeders and agronomists to develop and deploy new generation and promising varieties with farmer- and market -preferred attributes.

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

Luffa spp. (sponge gourd or ridge gourd) is an economically important vegetable crop widely cultivated in China, India and Southeast Asia. Here, we employed PacBio long-read single-molecule real-time (SMRT) sequencing to perform de novo genome assemblies of two commonly cultivated Luffa species, L. acutangula and L. cylindrica. We obtained preliminary draft genomes of 734.6 Mb and 689.8 Mb with scaffold N50 of 786,130 and 578,616 bases for L. acutangula and L. cylindrica, respectively. We also applied long-range Chicago and HiC techniques to obtain the first chromosome-scale whole-genome assembly of L. acutangula. The final assembly contained 13 pseudomolecules, corresponding to the haploid chromosome number in Luffa spp. (1n = 13, 2n = 26). The sizes of the assembled Luffa genomes are approximately twice as large as the genome assemblies of related Cucurbitaceae. A large proportion of L. acutangula (62.17%; 456.69 Mb) and L. cylindrica (56.78%; 391.65 Mb) genome assemblies contained repetitive elements. Phylogenetic analyses revealed that the substantial accumulation of transposable elements likely contributed to the expansion of the Luffa genomes. We also investigated alternative splicing events in Luffa using full-length transcript sequences obtained from PacBio Isoform Sequencing (Iso-seq). While the predominant form of alternative splicing in most plant species examined was intron retention, alternative 3' acceptor site selection appeared to be a major event observed in Luffa. High-quality genome assemblies for L. acutangula and L. cylindrica reported here provide valuable resources for Luffa breeding and future genetics and comparative genomics studies in Cucurbitaceae.

A high-quality sponge gourd (Luffa cylindrica) genome

Sponge gourd (Luffa cylindrica) is an important cultivated vegetable and medicinal plant in the family Cucurbitaceae. In this study, a draft genome sequence of the sponge gourd inbred line P93075 was analyzed. Using Illumina, PacBio, and 10× Genomics sequencing techniques as well as new assembly techniques such as FALCON and chromatin interaction mapping (Hi-C), a chromosome-scale genome of approximately 656.19 Mb, with an N50 scaffold length of 48.76 Mb, was generated. From this assembly, 25,508 protein-coding gene loci were identified, and 63.81% of the whole-genome consisted of transposable elements, which are major contributors to the expansion of the sponge gourd genome. According to a phylogenetic analysis of conserved genes, the sponge gourd lineage diverged from the bitter gourd lineage approximately 41.6 million years ago. Additionally, many genes that respond to biotic and abiotic stresses were found to be lineage specific or expanded in the sponge gourd genome, as demonstrated by the presence of 462 NBS-LRR genes, a much greater number than are found in the genomes of other cucurbit species; these results are consistent with the high stress resistance of sponge gourd. Collectively, our study provides insights into genome evolution and serves as a valuable reference for the genetic improvement of sponge gourd.

Construction of a high-density genetic linkage map and identification of gene controlling resistance to cucumber mosaic virus in Luffa cylindrica (L.) Roem. based on specific length amplified fragment sequencing

Luffa cylindrica L. is a cash crop which has important health, medicinal and industrial value, but no high saturation genetic map has been constructed owing to a lack of efficient markers. Furthermore, no genes were reportedly responsible for CMV resistance in Luffa spp. Specific length amplified fragment sequencing (SLAF-seq) is a valuable tool for large-scale discovery of markers and genetic mapping. The present study reported the construction of a high-density genetic map and the mapping of CMV resistant genes by using an F₂ population of 130 individuals and their two inbred line parents. A total of 271.01 Mb pair-end reads were generated. 100,077 high-quality SLAFs were detected, and 7404 of them were polymorphic. Finally, 3701 of the polymorphic markers were selected for genetic map construction, and 13 linkage groups were generated. The map spanned 1518.56 cM with an average distance of 0.41 cM between adjacent markers. Based on the newly constructed high-density map, one gene located on chromosome 1 (100.072-100.457 cM) was identified to regulate CMV resistance in L. cylindrica. A gag-polypeptide of LTR copia-type retrotransposon was predicted as the candidate gene responsible for CMV resistance in L. cylindrica. The high-density genetic map and the CMV resistant gene mapped and predicted in this study will be useful in future research.

The genetic basis of hybrid male sterility in sympatric Primulina species

BACKGROUND

Sympatric sister species provide an opportunity to investigate the genetic mechanisms and evolutionary forces that maintain species boundaries. The persistence of morphologically and genetically distinct populations in sympatry can only occur if some degree of reproductive isolation exists. A pair of sympatric sister species of Primulina (P. depressa and P. danxiaensis) was used to explore the genetic architecture of hybrid male sterility.

RESULTS

We mapped one major- and seven minor-effect quantitative trait loci (QTLs) that underlie pollen fertility rate (PFR). These loci jointly explained 55.4% of the phenotypic variation in the F2 population. A Bateson-Dobzhansky-Muller (BDM) model involving three loci was observed in this system. We found genotypic correlations between hybrid male sterility and flower morphology, consistent with the weak but significant phenotypic correlations between PFR and floral traits.

CONCLUSIONS

Hybrid male sterility in Primulina is controlled by a polygenic genetic basis with a complex pattern. The genetic incompatibility involves a three-locus BDM model. Hybrid male sterility is genetically correlated with floral morphology and divergence hitchhiking may occur between them.

Genome Survey Sequencing of Luffa Cylindrica L. and Microsatellite High Resolution Melting (SSR-HRM) Analysis for Genetic Relationship of Luffa Genotypes

Luffa cylindrica (L.) Roem. is an economically important vegetable crop in China. However, the genomic information on this species is currently unknown. In this study, for the first time, a genome survey of L. cylindrica was carried out using next-generation sequencing (NGS) technology. In total, 43.40 Gb sequence data of L. cylindrica, about 54.94× coverage of the estimated genome size of 789.97 Mb, were obtained from HiSeq 2500 sequencing, in which the guanine plus cytosine (GC) content was calculated to be 37.90%. The heterozygosity of genome sequences was only 0.24%. In total, 1,913,731 contigs (>200 bp) with 525 bp N₅₀ length and 1,410,117 scaffolds (>200 bp) with 885.01 Mb total length were obtained. From the initial assembled L. cylindrica genome, 431,234 microsatellites (SSRs) (≥5 repeats) were identified. The motif types of SSR repeats included 62.88% di-nucleotide, 31.03% tri-nucleotide, 4.59% tetra-nucleotide, 0.96% penta-nucleotide and 0.54% hexa-nucleotide. Eighty genomic SSR markers were developed, and 51/80 primers could be used in both “Zheda 23” and “Zheda 83”. Nineteen SSRs were used to investigate the genetic diversity among 32 accessions through SSR-HRM analysis. The unweighted pair group method analysis (UPGMA) dendrogram tree was built by calculating the SSR-HRM raw data. SSR-HRM could be effectively used for genotype relationship analysis of Luffa species.