Differences in transcriptomic responses upon Phytophthora palmivora infection among cultivars reveal potential underlying resistant mechanisms in durian.

BACKGROUND: Phytophthora palmivora is a devastating oomycete pathogen in durian, one of the most economically important crops in Southeast Asia. The use of fungicides in Phytophthora management may not be a long-term solution because of emerging chemical resistance issues. It is crucial to develop Phytophthora-resistant durian cultivars, and information regarding the underlying resistance mechanisms is valuable for smart breeding programs. RESULTS: In this study, we conducted RNA sequencing (RNA-seq) to investigate early gene expression responses (at 8, 24, and 48 h) after the P. palmivora infection in three durian cultivars, which included one resistant cultivar (Puangmanee; PM) and two susceptible cultivars (Monthong; MT and Kradumthong; KD). We performed co-expression and differential gene expression analyses to capture gene expression patterns and identify the differentially expressed genes. The results showed that genes encoding heat shock proteins (HSPs) were upregulated in all infected durians. The expression levels of genes encoding HSPs, such as ERdj3B, were high only in infected PM. A higher level of P. palmivora resistance in PM appeared to be associated with higher expression levels of various genes encoding defense and chitin response proteins, such as lysM domain receptor-like kinases. MT had a lower resistance level than PM, although it possessed more upregulated genes during P. palmivora infection. Many photosynthetic and defense genes were upregulated in the infected MT, although their expression levels were lower than those in the infected PM. KD, the least resistant cultivar, showed downregulation of genes involved in cell wall organization or biogenesis during P. palmivora infection. CONCLUSIONS: Our results showed that the three durian cultivars exhibited significantly different gene expression patterns in response to P. palmivora infection. The upregulation of genes encoding HSPs was common in all studied durians. The high expression of genes encoding chitin response proteins likely contributed to P. palmivora resistance in durians. Durian susceptibility was associated with low basal expression of defense genes and downregulation of several cell wall-related genes. These findings enhance our understanding of durian resistance to Phytophthora infection and could be useful for the development of elite durian cultivars.

De novo chromosome-level assembly of the Centella asiatica genome.

Centella asiatica is a herbaceous, perennial species indigenous to India and Southeast Asia. C. asiatica possesses several medicinal properties: anti-aging, anti-inflammatory, wound healing and memory enhancing. The lack of available genomics resources significantly impedes the improvement of C. asiatica varieties through molecular breeding. Here, we combined the 10× Genomics linked-read technology and the long-range HiC technique to obtain the genome assembly. The final assembly contained nine pseudomolecules, corresponding to the haploid chromosome number in C. asiatica. These nine chromosomes covered 402,536,584 bases or 93.6% of the 430-Mb assembly. Comparative genomics analyses based on single-copy orthologous genes showed that C. asiatica and the common ancestor of Coriandrum sativum (coriander) and Daucus carota (carrot) diverged about 48 million years ago. This assembly provides a valuable reference genome for future molecular studies, varietal development through marker-assisted breeding and comparative genomics studies in C. asiatica.

Taxonomic profiling of Symbiodiniaceae and bacterial communities associated with Indo-Pacific corals in the Gulf of Thailand using PacBio sequencing of full-length ITS and 16S rRNA genes.

Corals live with complex assemblages of microbes including bacteria, the dinoflagellate Symbiodiniaceae, fungi and viruses in a coral holobiont. These coral-associated microorganisms play an important role in their host fitness and survival. Here, we investigated the structure and diversity of algal and bacterial communities associated with five Indo-Pacific coral species, using full-length 16S rRNA and internal transcribed spacer sequences. While the dinoflagellate communities associated with Poriteslutea were dominated with Symbiodiniaceae genus Cladocopium, the other four coral hosts were associated mainly with members of the Durusdinium genus, suggesting that host species was one of the underlying factors influencing the structure and composition of dinoflagellate communities associated with corals in the Gulf of Thailand. Alphaproteobacteria dominated the microbiomes of Pocillopora spp. while Pavonafrondifera and P. lutea were associated primarily with Gammaproteobacteria. Finally, we demonstrated a superior performance of full-length 16S rRNA sequences in achieving species-resolution taxonomic classification of coral-associated microbiota.

ACC oxidase and miRNA 159a, and their involvement in fresh fruit bunch yield (FFB) via sex ratio determination in oil palm.

Oil palm (Elaeis guineesis Jacq.) is the most productive oil-bearing crop, yielding more oil per area than any other oil-bearing crops. However, there are still efforts to improve oil palm yield, in order to serve consumer and manufacturer demand. Oil palm produces female and male inflorescences in an alternating cycle. So, high sex ratio (SR), the ratio of female inflorescences to the total inflorescences, is a favorable trait in term of increasing yields in oil palm. This study aims to understand the genetic control for SR related traits, such as fresh fruit bunch yield (FFB), by characterizing genes at FFB quantitative trait loci (QTLs) on linkage 10 (chromosome 6) and linkage 15 (chromosome 10). Published oil palm sequences at the FFB QTLs were used to develop gene-based and simple sequence repeat (SSR) markers. We used the multiple QTL analysis model (MQM) to characterize the relationship of new markers with the SR traits in the oil palm population. The RNA expression of the most linked QTL genes was also evaluated in various tissues of oil palm. We identified EgACCO1 (encoding aminocyclopropane carboxylate (ACC) oxidase) at chromosome 10 and EgmiR159a (microRNA 159a) at chromosome 6 to be the most linked QTL genes or determinants for FFB yield and/or female inflorescence number with a phenotype variance explained (PVE) from 10.4 to 15 % and suggest that these play the important roles in sex determination and differentiation in oil palm. The strongest expression of EgACCO1 and the predicted precursor of EgmiR159a was found in ovaries and, to a lesser extent, fruit development. In addition, highly normalized expression of EgmiR159a was found in female flowers. In summary, the QTL analysis and the RNA expression reveal that EgACCO1 and EgmiR159a are the potential genetic factors involved in female flower determination and hence would affect yield in oil palm. However, to clarify how these genetic factors regulate female flower determination, more investigation of their down regulation or target may be essential. Additionally, if more sex determination genes controlled by plant hormones are identified, it may possible to reveal a crosstalk of sex determination genes with hormones and environment factors.

Assembly and analysis of a male sterile rubber tree mitochondrial genome reveals DNA rearrangement events and a novel transcript.

BACKGROUND: The rubber tree, Hevea brasiliensis, is an important plant species that is commercially grown to produce latex rubber in many countries. The rubber tree variety BPM 24 exhibits cytoplasmic male sterility, inherited from the variety GT 1. RESULTS: We constructed the rubber tree mitochondrial genome of a cytoplasmic male sterile variety, BPM 24, using 454 sequencing, including 8 kb paired-end libraries, plus Illumina paired-end sequencing. We annotated this mitochondrial genome with the aid of Illumina RNA-seq data and performed comparative analysis. We then compared the sequence of BPM 24 to the contigs of the published rubber tree, variety RRIM 600, and identified a rearrangement that is unique to BPM 24 resulting in a novel transcript containing a portion of atp9. CONCLUSIONS: The novel transcript is consistent with changes that cause cytoplasmic male sterility through a slight reduction to ATP production efficiency. The exhaustive nature of the search rules out alternative causes and supports previous findings of novel transcripts causing cytoplasmic male sterility.

Transcriptome analysis of normal and mantled developing oil palm flower and fruit.

Elaeis guineensis (oil palm) accounts for a large and increasing proportion of the world's cooking oil production. Cloning via somatic embryogenesis results in a somaclonal variant known as mantled which produce fruit with little to no oil yield. The mantled phenotype is believed to be epigenetic in nature. We performed RNA-Seq on developing flower and fruit samples of normal and mantled oil palm to characterize their transcriptomes. We present expression data for all transcripts in normal and mantled flower and fruit samples. Many genes are differentially expressed, including several from pathways that may be the cause of the mantled phenotype if disrupted, such as genes involved in primary hormone responses, DNA replication and repair, chromatin remodeling and a gene involved in RNA mediated DNA methylation. In addition, the gene expression data for developing flower and fruit will serve as a valuable resource for oil palm genetics and genomic studies.

A de novo chromosome-scale assembly of the Lablab purpureus genome.

Introduction: Lablab (Lablab purpureus (L.) Sweet), an underutilized tropical legume crop, plays a crucial role in global food and nutritional security. To enhance our understanding of its genetic makeup towards developing elite cultivars, we sequenced and assembled a draft genome of L. purpureus accession PK2022T020 using a single tube long fragment read (stLFR) technique. Results and discussion: The preliminary assembly encompassed 367 Mb with a scaffold N50 of 4.3 Mb. To improve the contiguity of our draft genome, we employed a chromatin contact mapping (Hi-C) approach to obtain a pseudochromosome-level assembly containing 366 Mb with an N50 length of 31.1 Mb. A total of 327.4 Mb had successfully been anchored into 11 pseudomolecules, corresponding to the haploid chromosome number in lablab. Our gene prediction recovered 98.4% of the highly conserved orthologs based on the Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis. Comparative analyses utilizing sequence information from single-copy orthologous genes demonstrated that L. purpureus diverged from the last common ancestor of the Phaseolus/Vigna species approximately 27.7 million years ago. A gene family expansion analysis revealed a significant expansion of genes involved in responses to biotic and abiotic stresses. Our high-quality chromosome-scale reference assembly provides an invaluable genomic resource for lablab genetic improvement and future comparative genomics studies among legume species.

microRNA profiling of exosomes derived from plasma and their potential as biomarkers for Opisthorchis viverrini-associated cholangiocarcinoma.

Cholangiocarcinoma (CCA) is a life-threatening disease that impacts patients worldwide. In Southeast Asian countries, the liver fluke Opisthorchis viverrini plays a major role in inducing carcinogenesis of the bile ducts. Due to its asymptomatic nature, O. viverrini infections are rarely treated, consequently leading to the development of advanced stages of CCA before diagnosis. Despite the current use of exosomal microRNAs (miRNA) as diagnostic biomarkers for the early detection of many types of cancer, the applications for miRNA remain limited with CCA. Circulating exosomes, membranous vesicles essential for intercellular communication, were found to contain unique miRNA. In this study, we conducted next-generation sequencing (Ion Torrent PGM) and bioinformatics to characterize and compare the contents of exosomal miRNA derived from the plasma of CCA patients, O. viverrini-infected patients, and healthy individuals, as well as to identify and validate key molecules as markers for screening the diagnosis of CCA and O. viverrini infection. The obtained results showed the success of using NGS technology in discovering exosomal miRNAs, specifically miR-194-5p and miR-192-5p, both of which were upregulated in the O. viverrini-infected group. Interestingly, miR-192-5p was upregulated while miR-194-5p was downregulated in CCA, suggesting their potential use as biomarkers for screening CCA and O. viverrini infection, especially in O. viverrini-endemic areas.

The complete chloroplast genome sequence and phylogenetic analysis of Heritiera fomes Buch.-Ham. (Malvales: Sterculiaceae).

Heritiera fomes Buch.-Ham. (1800) is a species of mangrove in the family Malvaceae, widely distributed in the Indo-Pacific and listed as 'endangered' (EN) on the International Union for Conservation of Nature's (IUCN) red list. We reported the complete chloroplast genome sequence of H. fomes. The genome was 168,521 bp in length and included two inverted repeats (IRs) of 34,496 bp, separated by a large single-copy (LSC) region of 88,604 bp and a small single-copy (SSC) region of 10,925 bp, respectively. The genome contained 87 protein-coding genes (PCGs), 8 rRNA genes, and 37 tRNA genes. The maximum-likelihood (ML) phylogenetic tree suggested that H. fomes is closely related to Heritiera angustata and Heritiera parvifolia with relatively high support bootstrap values of 86% and 100% with other species (Heritiera littoralis and Heritiera javanica), suggesting a relatively close genetic relationship between the five Heritiera plants. The chloroplast genome sequence provided a useful resource for conservation genetics studies of H. fomes and for phylogenetic studies of Heritiera.

Genome assembly of the Pendlebury's roundleaf bat, Hipposideros pendleburyi, revealed the expansion of Tc1/Mariner DNA transposons in Rhinolophoidea.

Bats (Chiroptera) constitute the second largest order of mammals and have several distinctive features, such as true self-powered flight and strong immunity. The Pendlebury's roundleaf bat, Hipposideros pendleburyi, is endemic to Thailand and listed as a vulnerable species. We employed the 10× Genomics linked-read technology to obtain a genome assembly of H. pendleburyi. The assembly size was 2.17 Gb with a scaffold N50 length of 15,398,518 bases. Our phylogenetic analysis placed H. pendleburyi within the rhinolophoid clade of the suborder Yinpterochiroptera. A synteny analysis showed that H. pendleburyi shared conserved chromosome segments (up to 105 Mb) with Rhinolophus ferrumequinum and Phyllostomus discolor albeit having different chromosome numbers and belonging different families. We found positive selection signals in genes involved in inflammation, spermatogenesis and Wnt signalling. The analyses of transposable elements suggested the contraction of short interspersed nuclear elements (SINEs) and the accumulation of young mariner DNA transposons in the analysed hipposiderids. Distinct mariners were likely horizontally transferred to hipposiderid genomes over the evolution of this family. The lineage-specific profiles of SINEs and mariners might involve in the evolution of hipposiderids and be associated with the phylogenetic separations of these bats from other bat families.

Genome assemblies of Vigna reflexo-pilosa (créole bean) and its progenitors, Vigna hirtella and Vigna trinervia, revealed homoeolog expression bias and expression-level dominance in the allotetraploid.

Vigna reflexo-pilosa (créole bean) is a wild legume belonging to the subgenus Ceratoropis and is widely distributed in Asia. Créole bean is the only tetraploid species in the genus Vigna, and it has been shown to derive from the hybridization of Vigna hirtella and Vigna trinervia. In this study, we combined the long-read PacBio technology with the chromatin contact mapping (Hi-C) technique to obtain a chromosome-level assembly of V. reflexo-pilosa. The final assembly contained 998,724,903 bases with an N50 length of 42,545,650 bases. Our gene prediction recovered 99.4% of the highly conserved orthologs based on the BUSCO analysis. To investigate homoeolog expression bias and expression level dominance in the tetraploid, we also sequenced and assembled the genomes of its progenitors. Overall, the majority of the homoeolog pairs (72.9%) displayed no expression bias, and among those that exhibited biased expression, 16.3% showed unbalanced homoeolog expression bias toward the V. trinervia subgenome. Moreover, 41.2% and 36.2% of the expressed gene pairs exhibited transgressive expression and expression level dominance, respectively. Interestingly, the genome-wide expression level dominance in the tetraploid was biased toward the V. trinervia subgenome. The analysis of methylation patterns also revealed that the average methylation levels in coding regions were higher in the V. hirtella subgenome than those in the V. trinervia subgenome. The genomic/transcriptomic resources for these three species are useful not only for the development of elite cultivars in Vigna breeding programs but also to researchers studying comparative genomics and investigating genomic/epigenomic changes following polyploid events.

Comparative Analysis and Phylogenetic Relationships of Ceriops Species (Rhizophoraceae) and Avicennia lanata (Acanthaceae): Insight into the Chloroplast Genome Evolution between Middle and Seaward Zones of Mangrove Forests.

Ceriops and Avicennia are true mangroves in the middle and seaward zones of mangrove forests, respectively. The chloroplast genomes of Ceriops decandra, Ceriops zippeliana, and Ceriops tagal were assembled into lengths of 166,650, 166,083 and 164,432 bp, respectively, whereas Avicennia lanata was 148,264 bp in length. The gene content and gene order are highly conserved among these species. The chloroplast genome contains 125 genes in A. lanata and 129 genes in Ceriops species. Three duplicate genes (rpl2, rpl23, and trnM-CAU) were found in the IR regions of the three Ceriops species, resulting in expansion of the IR regions. The rpl32 gene was lost in C. zippeliana, whereas the infA gene was present in A. lanata. Short repeats (<40 bp) and a lower number of SSRs were found in A. lanata but not in Ceriops species. The phylogenetic analysis supports that all Ceriops species are clustered in Rhizophoraceae and A. lanata is in Acanthaceae. In a search for genes under selective pressures of coastal environments, the rps7 gene was under positive selection compared with non-mangrove species. Finally, two specific primer sets were developed for species identification of the three Ceriops species. Thus, this finding provides insightful genetic information for evolutionary relationships and molecular markers in Ceriops and Avicennia species.

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.

A Chromosome-Scale Genome Assembly of Mitragyna speciosa (Kratom) and the Assessment of Its Genetic Diversity in Thailand.

Mitragyna speciosa (Kratom) is a tropical narcotic plant native to Southeast Asia with unique pharmacological properties. Here, we report the first chromosome-scale assembly of the M. speciosa genome. We employed PacBio sequencing to obtain a preliminary assembly, which was subsequently scaffolded using the chromatin contact mapping technique (Hi-C) into 22 pseudomolecules. The final assembly was 692 Mb with a scaffold N50 of 26 Mb. We annotated a total of 39,708 protein-coding genes, and our gene predictions recovered 98.4% of the highly conserved orthologs based on the BUSCO analysis. The phylogenetic analysis revealed that M. speciosa diverged from the last common ancestors of Coffea arabica and Coffea canephora approximately 47.6 million years ago. Our analysis of the sequence divergence at fourfold-degenerate sites from orthologous gene pairs provided evidence supporting a genome-wide duplication in M. speciosa, agreeing with the report that members of the genus Mitragyna are tetraploid. The STRUCTURE and principal component analyses demonstrated that the 85 M. speciosa accessions included in this study were an admixture of two subpopulations. The availability of our high-quality chromosome-level genome assembly and the transcriptomic resources will be useful for future studies on the alkaloid biosynthesis pathway, as well as comparative phylogenetic studies in Mitragyna and related species.

A draft chromosome-scale genome assembly of a commercial sugarcane.

Sugarcane accounts for a large portion of the worlds sugar production. Modern commercial cultivars are complex hybrids of S. officinarum, S. spontaneum, and several other Saccharum species, resulting in an auto-allopolyploid with 8-12 copies of each chromosome. The current genome assembly gold standard is to generate a long read assembly followed by chromatin conformation capture sequencing to scaffold. We used the PacBio RSII and chromatin conformation capture sequencing to sequence and assemble the genome of a South East Asian commercial sugarcane cultivar, known as Khon Kaen 3. The Khon Kaen 3 genome assembled into 104,477 contigs totalling 7 Gb, which scaffolded into 56 pseudochromosomes containing 5.2 Gb of sequence. Genome annotation produced 242,406 genes from 30,927 orthogroups. Aligning the Khon Kaen 3 genome sequence to S. officinarum and S. spontaneum revealed a high level of apparent recombination, indicating a chimeric assembly. This assembly error is explained by high nucleotide identity between S. officinarum and S. spontaneum, where 91.8% of S. spontaneum aligns to S. officinarum at 94% identity. Thus, the subgenomes of commercial sugarcane are so similar that using short reads to correct long PacBio reads produced chimeric long reads. Future attempts to sequence sugarcane must take this information into account.

Single nucleotide polymorphism marker development in the rubber tree, Hevea brasiliensis (Euphorbiaceae).

PREMISE OF THE STUDY: We demonstrated the application of high-throughput 454 sequencing technology in the identification of single nucleotide polymorphism (SNP) markers in Hevea brasiliensis. METHODS AND RESULTS: A total of 5883 putative SNP positions were discovered in silico, and 10 biallelic SNP markers were validated from 454-derived EST sequences. The polymorphism information content (PIC) and the observed heterozygosity (H(o)) ranged from 0.0963-0.5135 and 0.1071-0.4643, respectively. CONCLUSIONS: These markers can be useful for the construction of genetic maps, the identification of quantitative trait loci linked to commercially desirable traits, and the study of genetic structure in H. brasiliensis.

Iron homeostasis in the absence of ferricrocin and its consequences in fungal development and insect virulence in Beauveria bassiana.

The putative ferricrocin synthetase gene ferS in the fungal entomopathogen Beauveria bassiana BCC 2660 was identified and characterized. The 14,445-bp ferS encodes a multimodular nonribosomal siderophore synthetase tightly clustered with Fusarium graminearum ferricrocin synthetase. Functional analysis of this gene was performed by disruption with the bar cassette. ΔferS mutants were verified by Southern and PCR analyses. HPLC and TLC analyses of crude extracts indicated that biosynthesis of ferricrocin was abolished in ΔferS. Insect bioassays surprisingly indicated that ΔferS killed the Spodoptera exigua larvae faster (LT50 59 h) than wild type (66 h). Growth and developmental assays of the mutant and wild type demonstrated that ΔferS had a significant increase in germination under iron depletion and radial growth and a decrease in conidiation. Mitotracker staining showed that the mitochondrial activity was enriched in ΔferS under both iron excess and iron depletion. Comparative transcriptomes between wild type and ΔferS indicated that the mutant was increased in the expression of eight cytochrome P450 genes and those in iron homeostasis, ferroptosis, oxidative stress response, ergosterol biosynthesis, and TCA cycle, compared to wild type. Our data suggested that ΔferS sensed the iron excess and the oxidative stress and, in turn, was up-regulated in the antioxidant-related genes and those in ergosterol biosynthesis and TCA cycle. These increased biological pathways help ΔferS grow and germinate faster than the wild type and caused higher insect mortality than the wild type in the early phase of infection.

Transcriptome sequencing revealed the influence of blue light on the expression levels of light-stress response genes in Centella asiatica.

Centella asiatica is rich in medical and cosmetic properties. While physiological responses of C. asiatica to light have been widely reported, the knowledge of the effects of light on its gene expression is sparse. In this study, we used RNA sequencing (RNA-seq) to investigate the expression of the C. asiatica genes in response to monochromatic red and blue light. Most of the differentially expressed genes (DEGs) under blue light were up-regulated but those under red light were down-regulated. The DEGs encoded for CRY-DASH and UVR3 were among up-regulated genes that play significant roles in responses under blue light. The DEGs involved in the response to photosystem II photodamages and in the biosynthesis of photoprotective xanthophylls were also up-regulated. The expression of flavonoid biosynthetic DEGs under blue light was up-regulated but that under red light was down-regulated. Correspondingly, total flavonoid content under blue light was higher than that under red light. The ABI5, MYB4, and HYH transcription factors appeared as hub nodes in the protein-protein interaction network of the DEGs under blue light while ERF38 was a hub node among the DEGs under red light. In summary, stress-responsive genes were predominantly up-regulated under blue light to respond to stresses that could be induced under high energy light. The information obtained from this study can be useful to better understand the responses of C. asiatica to different light qualities.

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 chromosome-scale assembly of the black gram (Vigna mungo) genome.

Black gram (Vigna mungo) is an important short duration grain legume crop. Black gram seeds provide an inexpensive source of dietary protein. Here, we applied the 10X Genomics linked-read technology to obtain a de novo whole genome assembly of V. mungo cultivated variety Chai Nat 80 (CN80). The preliminary assembly contained 12,228 contigs and had an N50 length of 5.2 Mb. Subsequent scaffolding using the long-range Chicago and HiC techniques yielded the first high-quality, chromosome-level assembly of 499 Mb comprising 11 pseudomolecules. Comparative genomics analyses based on sequence information from single-copy orthologous genes revealed that black gram and mungbean (Vigna radiata) diverged about 2.7 million years ago . The transversion rate (4DTv) analysis in V. mungo revealed no evidence supporting a recent genome-wide duplication event observed in the tetraploid créole bean (Vigna reflexo-pilosa). The proportion of repetitive elements in the black gram genome is slightly lower than the numbers reported for related Vigna species. The majority of long terminal repeat retrotransposons appeared to integrate into the genome within the last five million years. We also examined alternative splicing events in V. mungo using full-length transcript sequences. While intron retention was the most prevalent mode of alternative splicing in several plant species, alternative 3' acceptor site selection represented the majority of events in black gram. Our high-quality genome assembly along with the genomic variation information from the germplasm provides valuable resources for accelerating the development of elite varieties through marker-assisted breeding and for future comparative genomics and phylogenetic studies in legume species.