Genes controlling hydrolysate toxin tolerance identified by QTL analysis of the natural Saccharomyces cerevisiae BCC39850

Lignocellulosic material can be converted to valorized products such as fuels. Pretreatment is an essential step in conversion, which is needed to increase the digestibility of the raw material for microbial fermentation. However, pretreatment generates by-products (hydrolysate toxins) that are detrimental to microbial growth. In this study, natural Saccharomyces strains isolated from habitats in Thailand were screened for their tolerance to synthetic hydrolysate toxins (synHTs). The Saccharomyces cerevisiae natural strain BCC39850 (toxin-tolerant) was crossed with the laboratory strain CEN.PK2-1C (toxin-sensitive), and quantitative trait locus (QTL) analysis was performed on the segregants using phenotypic scores of growth (OD600) and glucose consumption. VMS1, DET1, KCS1, MRH1, YOS9, SYO1, and YDR042C were identified from QTLs as candidate genes associated with the tolerance trait. CEN.PK2-1C knockouts of the VMS1, YOS9, KCS1, and MRH1 genes exhibited significantly greater hydrolysate toxin sensitivity to growth, whereas CEN.PK2-1C knock-ins with replacement of VMS1 and MRH1 genes from the BCC39850 alleles showed significant increased ethanol production titers compared with the CEN.PK2-1C parental strain in the presence of synHTs. The discovery of VMS1, YOS9, MRH1, and KCS1 genes associated with hydrolysate toxin tolerance in S. cerevisiae indicates the roles of the endoplasmic-reticulum-associated protein degradation pathway, plasma membrane protein association, and the phosphatidylinositol signaling system in this trait. KEY POINTS: • QTL analysis was conducted using a hydrolysate toxin-tolerant S. cerevisiae natural strain • Deletion of VMS1, YOS9, MRH1, and KCS1 genes associated with hydrolysate toxin-sensitivity • Replacement of VMS1 and MRH1 with natural strain alleles increased ethanol production titers in the presence of hydrolysate toxins.

De novo assembly and analysis of Sonneratia ovata genome and population analysis

Mangroves are an important part of coastal and estuarine ecosystems where they serve as nurseries for marine species and prevent coastal erosion. Here we report the genome of Sonneratia ovata, which is a true mangrove that grows in estuarine environments and can tolerate moderate salt exposure. We sequenced the S. ovata genome and assembled it into chromosome-level scaffolds through the use of Hi-C. The genome is 212.3 Mb and contains 12 chromosomes that range in size from 12.2 to 23.2 Mb. Annotation identified 29,829 genes with a BUSCO completeness of 95.9%. We identified salt genes and found copy number expansion of salt genes such as ADP-ribosylation factor 1, and elongation factor 1-alpha. Population analysis identified a low level of genetic variation and a lack of population structure within S. ovata.

Comprehensive genome assembly reveals genetic diversity and carcass consumption insights in critically endangered Asian king vultures

The Asian king vulture (AKV), a vital forest scavenger, is facing globally critical endangerment. This study aimed to construct a reference genome to unveil the mechanisms underlying its scavenger abilities and to assess the genetic relatedness of the captive population in Thailand. A reference genome of a female AKV was assembled from sequencing reads obtained from both PacBio long-read and MGI short-read sequencing platforms. Comparative genomics with New World vultures (NWVs) and other birds in the Family Accipitridae revealed unique gene families in AKV associated with retroviral genome integration and feather keratin, contrasting with NWVs' genes related to olfactory reception. Expanded gene families in AKV were linked to inflammatory response, iron regulation and spermatogenesis. Positively selected genes included those associated with anti-apoptosis, immune response and muscle cell development, shedding light on adaptations for carcass consumption and high-altitude soaring. Using restriction site-associated DNA sequencing (RADseq)-based genome-wide single nucleotide polymorphisms (SNPs), genetic relatedness and inbreeding status of five captive AKVs were determined, revealing high genomic inbreeding in two females. In conclusion, the AKV reference genome was established, providing insights into its unique characteristics. Additionally, the potential of RADseq-based genome-wide SNPs for selecting AKV breeders was demonstrated.

Determination of Morpho-Physiological Traits for Assessing Drought Tolerance in Sugarcane

Drought is a significant constraint to sugarcane productivity. Therefore, understanding how different varieties of sugarcane respond to drought stress can facilitate breeding programs and set up criteria for selecting drought-tolerant varieties. In the present study, we examined eight morpho-physiological traits to distinguish 40 sugarcane genotypes categorized into four groups based on significant differences in cane yield under non-stressed conditions and reduction of cane yield under drought-stressed conditions. The study was conducted during the formative stage in a greenhouse, encompassing both control and drought conditions. Drought treatments resulted in significant changes and differences in the mean values of various morpho-physiological traits. The hierarchical clustering analysis, utilizing stay-green traits such as higher chlorophyll fluorescence ratio (Fv/Fm), leaf chlorophyll content (SPAD), leaf relative water content (RWC), and lower leaf rolling score (LR), leaf drying score (LD), and drought recovery score (DR), successfully grouped 40 sugarcane genotypes into four major clusters, similar to the previously categorized groups. Correlation analysis showed significant relationships among cane yield, reduction of cane yield under drought conditions, and the stay-green traits. Our results demonstrated that morpho-physiological traits contributing to the "stay-green" phenotypes could be useful as selection criteria for drought tolerance in sugarcane.

Adaptation of African swine fever virus to MA-104 cells: Implications of unique genetic variations

African swine fever virus (ASFV) is a large, double-stranded DNA virus that causes a fatal, contagious disease specifically in pigs. However, prevention and control of ASFV outbreaks have been hampered by the lack of an effective vaccine or antiviral treatment for ASFV. Although ASFV has been reported to adapt to a variety of continuous cell lines, the phenotypic and genetic changes associated with ASFV adaptation to MA-104 cells remain poorly understood. Here, we adapted ASFV field isolates to efficiently propagate through serial viral passages in MA-104 cells. The adapted ASFV strain developed a pronounced cytopathic effect and robust infection in MA-104 cells. Interestingly, the adapted variant maintained its tropism in primary porcine kidney macrophages. Whole genome analysis of the adapted virus revealed unique gene deletions in the left and right variable regions of the viral genome compared to other previously reported cell culture-adapted ASFV strains. Notably, gene duplications at the 5' and 3' ends of the viral genome were in reverse complementary alignment with their paralogs. Single point mutations in protein-coding genes and intergenic regions were also observed in the viral genome. Collectively, our results shed light on the significance of these unique genetic changes during adaptation, which facilitate the growth of ASFV in MA-104 cells.

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.

Pins Gene Table v2.0: An Online Genome Database of 37 Pythium insidiosum Strains for Gene Content Exploration and Phylogenomic Analysis

Unlike most pathogenic oomycetes, Pythium insidiosum infects humans and animals instead of plants. P. insidiosum has three clinically relevant genotypes/clades that cause a severe disease called pythiosis. To develop strategies for infection control, it is necessary to understand the biology and pathogenesis of this pathogen. Investigating the evolutionary mechanisms behind the host-specific adaptation is vital, and comparative genomic analysis can help with this. To facilitate genomic analysis, an online bioinformatics tool called P. insidiosum (Pins) Gene Table v2.0 was developed. This tool includes genomic data from 37 genetically diverse P. insidiosum strains and four related species. The database contains 732,686 genes, grouped into 80,061 unique clusters and further divided into core and variable categories at genus, species, and genotype levels. A high-resolution phylogenomic relationship among P. insidiosum strains and other oomycetes was projected through hierarchical clustering and core gene analyses. 3156 P. insidiosum-specific genes were shared among all genotypes and may be responsible for causing disease in humans and animals. After comparing these species-specific genes to the MvirDB database, 112 had significant matches with 66 known virulence proteins, some of which might be involved in vascular occlusion, which is a pathological feature of pythiosis. The correlation of genotypes, geographic origins, and affected hosts of P. insidiosum suggests that clade-I strains are more specific to animals, while clade-II/III strains are more specific to humans. The clade-specific genes might link to host preference. In summary, Pins Gene Table v2.0 is a comprehensive genome database accessible to users with minimal bioinformatics experience for the analysis of P. insidiosum genomes.

A SNP variation in the Sucrose synthase (SoSUS) gene associated with sugar-related traits in sugarcane

Background

Sugarcane (Saccharum spp.) is an economically significant crop for both the sugar and biofuel industries. Breeding sugarcane cultivars with high-performance agronomic traits is the most effective approach for meeting the rising demand for sugar and biofuels. Molecular markers associated with relevant agronomic traits could drastically reduce the time and resources required to develop new sugarcane varieties. Previous sugarcane candidate gene association analyses have found single nucleotide polymorphism (SNP) markers associated with sugar-related traits. This study aims to validate these associated SNP markers of six genes, including Lesion simulating disease 1 (LSD), Calreticulin (CALR), Sucrose synthase 1 (SUS1), DEAD-box ATP-dependent RNA helicase (RH), KANADI1 (KAN1), and Sodium/hydrogen exchanger 7 (NHX7), in a diverse population in 2-year and two-location evaluations.

Methods

After genotyping of seven targeted SNP markers was performed by PCR Allelic Competitive Extension (PACE) SNP genotyping, the association with sugar-related traits and important cane yield component traits was determined on a set of 159 sugarcane genotypes. The marker-trait relationships were validated and identified by both t-test analysis and an association analysis based on the general linear model.

Results

The mSoSUS1_SNPCh10.T/C and mSoKAN1_SNPCh7.T/C markers that were designed from the SUS1 and KAN1 genes, respectively, showed significant associations with different amounts of sugar-related traits and yield components. The mSoSUS1_SNPCh10.T/C marker was found to have more significant association with sugar-related traits, including pol, CCS, brix, fiber and sugar yield, with p values of 6.08 × 10-6 to 4.35 × 10-2, as well as some cane yield component traits with p values of 1.61 × 10-4 to 3.35 × 10-2. The significant association is consistent across four environments.

Conclusion

Sucrose synthase (SUS) is considered a crucial enzyme involved in sucrose metabolism. This marker is a high potential functional marker that may be used in sugarcane breeding programs to select superior sugarcane with good fiber and high sugar contents.

Genetic diversity of the dengue virus population in dengue fever and dengue hemorrhagic fever patients

BACKGROUND

The error-prone replication of dengue virus (DENV) in host results in the highly diverse viral population. Together with the host factor, intra-host diversity may influence the disease severity. Therefore, it is worth investigating whether there is a correlation between intra-host genetic diversity and disease severity.

OBJECTIVE

To investigate the genetic diversity in DENV for four serotypes of the dengue population from patients with dengue fever (DF) and dengue hemorrhagic fever (DHF) using next-generation sequencing (NGS) technology.

METHODS

Forty RNA samples categorized into eight groups by severity and serotypes were sequenced and analyzed for genetic variation. Analysis on the hot-cold genomic regions, selection pressure and correlation between genotype and disease severity were performed in this study.

RESULTS

Comparison between the NGS data of the DF and DHF specimens showed conservation between their major populations with the consensus sequences for DF and DHF sharing 99% similarity. However, the minor populations in DF and DHF were more diverse. Many genes in DF had an #NS/#S ratio higher than in DHF. Only NS4B of DENV1 DF has #NS/#S ratio higher than one. Hot regions of the DF were detected in NS3 of DENV1, DENV2 and Envelope of DENV3, whereas the hot regions of the DHF samples were detected in the small region in 3'UTR of DENV2 and DENV3.

CONCLUSIONS

Various explorations of the variations of DF and DHF were performed in this study. However, we have not yet found any specific characteristics of intra-host diversity associated with disease severity.

Comparative genomics and genome-wide SNPs of endangered Eld's deer provide breeder selection for inbreeding avoidance

Eld's deer, a conserved wildlife species of Thailand, is facing inbreeding depression, particularly in the captive Siamese Eld's deer (SED) subspecies. In this study, we constructed genomes of a male SED and a male Burmese Eld's deer (BED), and used genome-wide single nucleotide polymorphisms to evaluate the genetic purity and the inbreeding status of 35 SED and 49 BED with limited pedigree information. The results show that these subspecies diverged approximately 1.26 million years ago. All SED were found to be purebred. A low proportion of admixed SED genetic material was observed in some BED individuals. Six potential breeders from male SED with no genetic relation to any female SED and three purebred male BED with no relation to more than 10 purebred female BED were identified. This study provides valuable insights about Eld's deer populations and appropriate breeder selection in efforts to repopulate this endangered species while avoiding inbreeding.

MGI short-read genome assemblies of Pythium insidiosum (reclassified as Pythium periculosum) strains Pi057C3 and Pi050C3

OBJECTIVES

Pythium insidiosum causes a difficult-to-treat infectious condition called pythiosis, with high morbidity and mortality. So far, genome data of at least 10 strains of P. insidiosum, primarily classified in the phylogenetic clades I and II, have been sequenced using various next-generation sequencing platforms. The MGI short-read platform was employed to obtain genome data of 2 clade-III strains of P. insidiosum (recently reclassified as Pythium periculosum) from patients in Thailand and the United States. This work is a part of our attempt to generate a comprehensive genome database from diverse pathogen strains.

DATA DESCRIPTION

A 150-bp paired-end library was prepared from a gDNA sample of P. insidiosum (P. periculosum) strains Pi057C3 and Pi050C3 (also known as ATCC90586) to generate draft genome sequences using an MGISEQ-2000RS sequencer. As a result, for the strain Pi057C3, we obtained a 42.5-Mb assembled genome (164x coverage) comprising 14,134 contigs, L50 of 241, N50 of 45,748, 57.6% CG content, and 12,147 ORFs. For the strain Pi050C3, we received a 43.3-Mb draft genome (230x coverage) containing 14,511 contigs, L50 of 245, N50 of 45,208, 57.7% CG content, and 12,249 ORFs. The genome sequences have been deposited in the NCBI/DDBJ databases under the accession numbers JAKCXM000000000.1 (strain Pi057C3) and JAKCXL000000000.1 (strain Pi050C3).

PacBio long read-assembled draft genome of Pythium insidiosum strain Pi-S isolated from a Thai patient with pythiosis

OBJECTIVES

Pythium insidiosum is the causative agent of pythiosis, a difficult-to-treat condition, in humans and animals worldwide. Biological information about this filamentous microorganism is sparse. Genomes of several P. insidiosum strains were sequenced using the Illumina short-read NGS platform, producing incomplete genome sequence data. PacBio long-read platform was employed to obtain a better-quality genome of Pythium insidiosum. The obtained genome data could promote basic research on the pathogen's biology and pathogenicity.

DATA DESCRIPTION

gDNA sample was extracted from the P. insidiosum strain Pi-S for whole-genome sequencing by PacBio long-read NGS platform. Raw reads were assembled using CANU (v2.1), polished using ARROW (SMRT link version 5.0.1), aligned with the original raw PacBio reads using pbmm2 (v1.2.1), consensus sequence checked using ARROW, and gene predicted using Funannotate pipeline (v1.7.4). The genome completion was assessed using BUSCO (v4.0.2). As a result, 840 contigs (maximum length: 1.3 Mb; N50: 229.9 Kb; L50: 70) were obtained. Sequence assembly showed a genome size of 66.7 Mb (178x coverage; 57.2% G-C content) that contained 20,375 ORFs. A BUSCO-based assessment revealed 85.5% genome completion. All assembled contig sequences have been deposited in the NCBI database under the accession numbers BBXB02000001 - BBXB02000840.

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.

A Large-Scale Candidate-Gene Association Mapping for Drought Tolerance and Agronomic Traits in Sugarcane

Dissection of the genetic loci controlling drought tolerance traits with a complex genetic inheritance is important for drought-tolerant sugarcane improvement. In this study, we conducted a large-scale candidate gene association study of 649 candidate genes in a sugarcane diversity panel to identify genetic variants underlying agronomic traits and drought tolerance indices evaluated in plant cane and ratoon cane under water-stressed (WS) and non-stressed (NS) environments. We identified 197 significant marker-trait associations (MTAs) in 141 candidate genes associated with 18 evaluated traits with the Bonferroni correction threshold (α = 0.05). Out of the total, 95 MTAs in 78 candidate genes and 62 MTAs in 58 candidate genes were detected under NS and WS conditions, respectively. Most MTAs were found only in specific water regimes and crop seasons. These MTAs explained 7.93-30.52% of phenotypic variation. Association mapping results revealed that 34, 59, and 104 MTAs involved physiological and molecular adaptation, phytohormone metabolism, and drought-inducible genes. They identified 19 pleiotropic genes associated with more than one trait and many genes related to drought tolerance indices. The genetic and genomic resources identified in this study will enable the combining of yield-related traits and sugar-related traits with agronomic value to optimize the yield of sugarcane cultivars grown under drought-stressed and non-stressed environments.

Assessment of the Genetic Diversity and Population Structure of Rhizophora mucronata along Coastal Areas in Thailand

Unique and biodiverse, mangrove ecosystems provide humans with benefits and contribute to coastal protection. Rhizophora mucronata, a member of the Rhizophoraceae family, is prevalent in the mangrove forests of Thailand. R. mucronata's population structure and genetic diversity have received scant attention. Here, we sequenced the entire genome of R. mucronata using 10× Genomics technology and obtained an assembly size of 219 Mb with the N50 length of 542,540 bases. Using 2857 single nucleotide polymorphism (SNP) markers, this study investigated the genetic diversity and population structure of 80 R. mucronata accessions obtained from the mangrove forests in Thailand. The genetic diversity of R. mucronata was moderate (I = 0.573, Ho = 0.619, He = 0.391). Two subpopulations were observed and confirmed from both population structure and principal component analysis (PCA). Analysis of molecular variance (AMOVA) showed that there was more variation within populations than between them. Mean pairwise genetic differentiation (FST = 0.09) showed that there was not much genetic difference between populations. Intriguingly, the predominant clustering pattern in the R. mucronata population did not correspond to the Gulf of Thailand and the Andaman Sea, which are separated by the Malay Peninsula. Several factors could have influenced the R. mucronata genetic pattern, such as hybridization and anthropogenic factors. This research will provide important information for the future conservation and management of R. mucronata in Thailand.

Resequencing of durian genomes reveals large genetic variations among different cultivars

Durian (Durio zibethinus), which yields the fruit known as the "King of Fruits," is an important economic crop in Southeast Asia. Several durian cultivars have been developed in this region. In this study, we resequenced the genomes of three popular durian cultivars in Thailand, including Kradumthong (KD), Monthong (MT), and Puangmanee (PM) to investigate genetic diversities of cultivated durians. KD, MT, and PM genome assemblies were 832.7, 762.6, and 821.6 Mb, and their annotations covered 95.7, 92.4, and 92.7% of the embryophyta core proteins, respectively. We constructed the draft durian pangenome and analyzed comparative genomes with related species in Malvales. Long terminal repeat (LTR) sequences and protein families in durian genomes had slower evolution rates than that in cotton genomes. However, protein families with transcriptional regulation function and protein phosphorylation function involved in abiotic and biotic stress responses appeared to evolve faster in durians. The analyses of phylogenetic relationships, copy number variations (CNVs), and presence/absence variations (PAVs) suggested that the genome evolution of Thai durians was different from that of the Malaysian durian, Musang King (MK). Among the three newly sequenced genomes, the PAV and CNV profiles of disease resistance genes and the expressions of methylesterase inhibitor domain containing genes involved in flowering and fruit maturation in MT were different from those in KD and PM. These genome assemblies and their analyses provide valuable resources to gain a better understanding of the genetic diversity of cultivated durians, which may be useful for the future development of new durian cultivars.

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.

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.

The complete chloroplast genome sequence of Intsia bijuga (Colebr.) Kuntze (Fabaceae: Detaroideae: Afzelieae)

Intsia bijuga (Colebr.) Kuntze. (1891) is a threatened mangrove species, belonging to the Fabaceae family and is native to the western Pacific coast and Southeast Asia. Here, we applied short-read Illumina technology to sequence and assemble its chloroplast genome. The complete chloroplast genome is 158,363 bp in length, composed of one large single-copy (LSC) region of 87,489 bp, one small single-copy (SSC) region of 19,438 bp, and a pair of inverted repeats (IRs) of 25,719 bp. A total of 129 unique genes were annotated, comprising 84 protein-coding genes, eight rRNA genes, and 37 tRNA genes. Our phylogenetic analysis showed the placement of I. bijuga (OL699920.1) with Afzelia species within Fabaceae family.

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.

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.

The complete chloroplast genome of Sonneratia griffithii Kurz (Lythraceae)

Sonneratia griffithii Kurz is a critically endangered mangrove species that can be found along the western coast of Thailand. In this study, we reported the complete chloroplast genome of S. griffithii. The chloroplast genome is 152,730 bp, consisting of one large single-copy (LSC) region, one small single-copy (SSC) region and a pair of inverted repeats (IRs). The LSC, SSC, and IR lengths are 87,226, 17,764, and 23,870 bp, respectively. The genome contains 113 unique genes, including 79 protein-coding, 30 tRNA, and 4 rRNA genes. The GC content of the chloroplast genome is 37.31%. The phylogenetic analysis based on 76 protein-coding genes showed a monophyletic group of S. griffithii and other Sonneratia species.

Mass spectrometry and synchrotron-FTIR microspectroscopy reveal the anti-inflammatory activity of Bua Bok extracts

INTRODUCTION

Bua Bok or Centella asiatica (CA) is an Asian vegetable with anti-inflammatory benefits. Asiaticoside, asiatic acid, madecassoside and madecassic have been characterised as major active ingredients with a wide range of pharmacological advantages. In manufacturing processes, high-performance liquid chromatography (HPLC) and liquid chromatography-mass spectrometry (LCMS) are used to routinely determine the active compounds in raw materials.

OBJECTIVES

This research aims to explore anti-inflammatory properties, characterise metabolites and observe the biochemical changes of the inflammatory induced macrophages after pretreatment with the potential extracted fractions.

METHODS

Bua Bok leaf extracts were prepared. Macrophages were pretreated with non-toxic fractions to determine the anti-inflammatory action. Tentative metabolites of effective fractions were identified by LC-MS. Synchrotron fourier-transform infrared (S-FTIR) microspectroscopy was utilised to observe the biochemical change of the lipopolysaccharide (LPS)-induced cells after pretreatment with potential fractions.

RESULTS

Fractions of ethyl acetate, 30% and 100% ethanol highly increased the nitrile scavenging and suppressed the function of phospholipase A₂ . Fractions of 70% and 100% ethanol strongly decreased nitric oxide production. The comparison of 39 chemical compounds was presented. The change of proteins was improved after pretreatment of macrophages with fraction 70% ethanol. Fraction of 100% ethanol revealed the lipid accumulation was lower than 70% ethanol and diclofenac.

CONCLUSION

While the anti-inflammatory actions of 70% and 100% ethanol were similar. S-FTIR expressed they inhibited inflammatory response with the distinct features of biomolecules. The S-FTIR, LC-MS and biological assay confidently provided the efficient strategies to inform the advantage of herbal extract on cellular organisation instead of a single compound.

Chromosome-level genome assembly of Indian mangrove (Ceriops tagal) revealed a genome-wide duplication event predating the divergence of Rhizophoraceae mangrove species

Mangrove ecosystems are unique, highly diverse, provide benefits to humans, and aid in coastal protection. The Indian mangrove, or spurred mangrove, [Ceriops tagal (Perr.) C. B. Rob.] is a member of the Rhizophoraceae family and is commonly found along the intertidal zones in tropical regions in Southeast Asia, southern Asia, and Africa. Here, we present the first high-quality reference genome assembly of the Ceriops species. A preliminary draft assembly, generated from the 10× Genomics linked-read library, was scaffolded using the proximity ligation chromatin contact mapping technique (Hi-C) to obtain a chromosome-scale assembly of 231,919,005 bases with an N50 length of 11,408,429 bases. The benchmarking universal single-copy orthologs (BUSCO) analysis revealed that C. tagal gene predictions recovered 95.8% of the highly conserved orthologs. Phylogenetic analyses suggested that C. tagal diverged from the last common ancestor of flat-leaf spurred mangrove [C. decandra (Griff.) Ding Hou] and C. zippeliana Blume ∼10.4 million yr ago (MYA), and the last common ancestor of genera Ceriops, Kandelia, and Rhizophora diverged from that of genus Bruguiera ∼49.4 MYA. In addition, our analysis of the transversion rate at fourfold-degenerate sites from orthologous gene pairs provided evidence supporting a recent whole-genome duplication in C. tagal. The STRUCTURE and principal component analyses illustrated that C. tagal individuals investigated in this study were the admixture of two subpopulations, the genetic background of which was influenced primarily by location. The availability of genomic and transcriptomic resources and biodiversity data reported in this work will be useful for future studies that may shed light on adaptive evolutions of mangrove species.

Quantitative analysis of methoxyflavones discriminates between the two types of Kaempferia parviflora

INTRODUCTION

Kaempferia parviflora or black ginger is abundantly cultivated because its rhizomes contain methoxyflavones that have many pharmacological properties. K. parviflora can be divided into two types, based on morphological characteristics, but differences in their chemical compositions have never been explored.

OBJECTIVES

This research aims to find chemical markers that can be used to differentiate between the two types of K. parviflora, the red-leaf and green-leaf types, by quantifying the amounts of methoxyflavones.

MATERIAL AND METHODS

K. parviflora samples were collected from 39 locations in Thailand. Their genetic diversity was assessed by a genotyping-by-sequencing (GBS) technique to construct the population structure. Their chemical compositions were analyzed by high performance liquid chromatography-photodiode array detection to determine the methoxyflavone contents.

RESULTS

The population structure based on >3,000 single nucleotide polymorphism (SNP) markers showed that the samples can be divided into two groups, which were consistent with the classification by leaf margin color (red-leaf and green-leaf types). HPLC analysis revealed 3,5,7,3',4'-pentamethoxyflavone (PMF), 5,7-dimethoxyflavone (DMF), 5,7,4'-trimethoxyflavone (TMF), 3,5,7-trimethoxyflavone and 3,5,7,4'-tetramethoxyflavone as major methoxyflavones that can be used as chemical markers. The red-leaf type showed higher amounts of PMF, TMF and 3,5,7,4'-tetramethoxyflavone than the green-leaf type, while the green-leaf type showed higher amounts of DMF and 3,5,7-trimethoxyflavone than the red-leaf type.

CONCLUSION

These results provide another approach to discriminate the two types of K. parviflora using chemical profiles alongside genetic and morphological analyses. Therefore, a specific type of K. parviflora can be selected over the other based on preferences for a certain methoxyflavone.

The complete mitochondrial genome sequence of the karst-dwelling crab, Terrapotamon thungwa (Crustacea: Brachyura: Potamidae)

Terrapotamon thungwa is a new species of terrestrial long-legged crab discovered in a karst landscape of southern Thailand. Here, we report the first complete mitochondrial genome of this crab species. The mitochondrial genome size is 16,156 base-pairs (bp), including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA), and two ribosomal RNA (rRNA) genes. The AT and GC content of the mitochondrial genome sequence is 73.2% and 26.8%, respectively. Phylogenetic analysis with 26 crustacean species, based on 13 mitochondrial conserve genes, showed that T. thungwa was closely related to other freshwater crab species in the family Potamidae.

Assembly of a hybrid mangrove, Bruguiera hainesii, and its two ancestral contributors, Bruguiera cylindrica and Bruguiera gymnorhiza

Mangroves are plants that live in tropical and subtropical coastal regions of the world, they are adapted to high salt environments and cyclic tidal flooding. Mangroves play important ecological roles, including acting as breeding grounds for many fish species and to prevent coastal erosion. The genomes of three mangrove species, Bruguiera gymnorhiza, Bruguiera cylindrica, and a hybrid of the two, Bruguiera hainesii were sequenced, assembled and annotated. The two progenitor species, B. gymnorhiza and B. cylindrica, were found to be highly similar to each other and sufficiently similar to B. parviflora to allow it to be used for reference based scaffolding to generate chromosome level scaffolds. The two subgenomes of B. hainesii were independently assembled and scaffolded. Analysis of B. hainesii confirms that it is a hybrid and the hybridisation event was estimated at 2.4 to 3.5 million years ago using a Bayesian Relaxed Molecular Clock approach.

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.

A SNP variation in an expansin (EgExp4) gene affects height in oil palm

Oil palm (Elaeis guineensis Jacq.), an Aracaceae family plant, is utilized for both consumable and non-consumable products, including cooking oil, cosmetics and biodiesel production. Oil palm is a perennial tree with 25 years of optimal harvesting time and a height of up to 18 m. However, harvesting of oil palm fruit bunches with heights of more than 2-3 meters is challenging for oil palm farmers. Thus, understanding the genetic control of height would be beneficial for using gene-based markers to speed up oil palm breeding programs to select semi-dwarf oil palm varieties. This study aims to identify Insertion/Deletions (InDels) and single nucleotide polymorphisms (SNPs) of five height-related genes, including EgDELLA1, EgGRF1, EgGA20ox1, EgAPG1 and EgExp4, in short and tall oil palm groups by PacBio SMRT sequencing technology. Then, the SNP variation's association with height was validated in the Golden Tenera (GT) population. All targeted genes were successfully amplified by two rounds of PCR amplification with expected sizes that ranged from 2,516 to 3,015 base pair (bp), covering 5' UTR, gene sequences and 3' UTR from 20 short and 20 tall oil palm trees. As a result, 1,166, 909, 1,494, 387 and 5,384 full-length genomic DNA sequences were revealed by PacBio SMRT sequencing technology, from EgDELLA1, EgGRF1, EgGA20ox1, EgAPG1 and EgExp4 genes, respectively. Twelve variations, including eight InDels and four SNPs, were identified from EgDELLA1, EgGRF1, EgGA20ox1 and EgExp4. No variation was found for EgAPG1. After SNP through-put genotyping of 4 targeted SNP markers was done by PACE™ SNP genotyping, the association with height was determined in the GT population. Only the mEgExp4_SNP118 marker, designed from EgExp4 gene, was found to associate with height in 2 of 4 height-recordings, with p values of 0.0383 for height (HT)-1 and 0.0263 for HT-4. In conclusion, this marker is a potential gene-based marker that may be used in oil palm breeding programs for selecting semi-dwarf oil palm varieties in the near future.

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.

The First Genetic Linkage Map of Winged Bean [Psophocarpus tetragonolobus (L.) DC.] and QTL Mapping for Flower-, Pod-, and Seed-Related Traits

Winged bean [Psophocarpus tetragonolobus (L.) DC.] (2n = 2× = 18) is a tropical legume crop with multipurpose usages. Recently, the winged bean has regained attention from scientists as a food protein source. Currently, there is no breeding program for winged bean cultivars. All winged bean cultivars are landraces or selections from landraces. Molecular markers and genetic linkage maps are pre-requisites for molecular plant breeding. The aim of this study was to develop a high-density linkage map and identify quantitative trait loci (QTLs) for pod and seed-related traits of the winged bean. An F₂ population of 86 plants was developed from a cross between winged bean accessions W054 and TPT9 showing contrasting pod length, and pod, flower and seed colors. A genetic linkage map of 1384 single nucleotide polymorphism (SNP) markers generated from restriction site-associated DNA sequencing was constructed. The map resolved nine haploid chromosomes of the winged bean and spanned the cumulative length of 4552.8 cM with the number of SNPs per linkage ranging from 36 to 218 with an average of 153.78. QTL analysis in the F₂ population revealed 31 QTLs controlling pod length, pod color, pod anthocyanin content, flower color, and seed color. The number of QTLs per trait varied between 1 (seed length) to 7 (banner color). Interestingly, the major QTLs for pod color, anthocyanin content, and calyx color, and for seed color and flower wing color were located at the same position. The high-density linkage map QTLs reported in this study will be useful for molecular breeding of winged beans.

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.

The complete mitochondrial genome of the Hipposideros pendleburyi (Pendlebury's leaf-nosed bat) an endemic species in Thailand

This study presents the first complete mitochondrial genome of the Hipposideros pendleburyi (Pendlebury's leaf-nosed bat), an endemic species in Thailand. The mitochondrial genome was 16,820 bp in length and contains 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. The overall base composition was 31.5% A, 26.2% T, 28.3% C, and 14.0% G. A maximum-likelihood tree revealed that H. pendleburyi was grouped with Hipposideros armiger within the Hipposideridae clade, which has Rhinolophidae as a sister clade.

The complete mitogenome of the Thai soldier crab Mictyris thailandensis Davie, Wisespongpand & Shih, 2013 (Crustacea: Decapoda: Mictyridae)

Mictyris thailandensis has been described recently in the family Mictyridae which is found only in the Andaman Sea, west coast of Thailand. In this study, we performed shotgun genome sequencing of a male M. thailandensis using a paired-end (150 bp) sequencing chemistry on MGISEQ-2000RS and report the complete mitochondrial genome of M. thailandensis (15,557 bp). A total of 37 genes have been annotated: 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a control region. Comparative phylogenetic analysis with 29 crustaceans based on 13 conserved genes demonstrated that M. thailandensis is closely related to other soldier crabs in the family Mictyridae. The mitogenome of M. thailandensis presented here provides useful genetic information to help understand the evolutionary relationships among the Mictyridae family members.

Complete chloroplast genome sequences of five Bruguiera species (Rhizophoraceae): comparative analysis and phylogenetic relationships

Bruguiera is a genus of true mangroves that are mostly distributed in the Indo-West Pacific region. However, the number of published whole chloroplast genome sequences of Bruguiera species are limited. Here, the complete chloroplast sequences of five Bruguiera species were sequenced and assembled using Illumina data. The chloroplast genomes of B. gymnorhiza, B. hainesii, B. cylindrica, B. parviflora and B. sexangula were assembled into 161,195, 164,295, 164,297, 163,228 and 164,170 bp, respectively. All chloroplast genomes contain 37 tRNA and eight rRNA genes, with either 84 or 85 protein-coding genes. A comparative analysis of these genomes revealed high similarity in gene structure, gene order and boundary position of the LSC, SSC and two IR regions. Interestingly, B. gymnorhiza lost a rpl32 gene in the SSC region. In addition, a ndhF gene in B. parviflora straddles both the SSC and IRB boundary regions. These genes reveal differences in chloroplast evolution among Bruguiera species. Repeats and SSRs in the chloroplast genome sequences were found to be highly conserved between B. cylindrica and B. hainesii as well as B. gymnorhiza and B. sexangula indicating close genetic relationships based on maternal inheritance. Notably, B. hainesii, which is considered a hybrid between B. gymnorhiza and B. cylindrica, appears to have inherited the chloroplast from B. cylindrica. Investigating the effects of selection events on shared protein-coding genes showed a positive selection in rps7 and rpl36 genes in all species compared to land-plant species. A phylogenetic analysis, based on 59 conserved chloroplast protein-coding genes, showed strong support that all Bruguiera species are in the clade Rhizophoraceae. This study provides valuable genetic information for the study of evolutionary relationships and population genetics in Bruguiera and other mangrove species.

Establishment of Human-Induced Pluripotent Stem Cell-Derived Neurons-A Promising In Vitro Model for a Molecular Study of Rabies Virus and Host Interaction

Rabies is a deadly viral disease caused by the rabies virus (RABV), transmitted through a bite of an infected host, resulting in irreversible neurological symptoms and a 100% fatality rate in humans. Despite many aspects describing rabies neuropathogenesis, numerous hypotheses remain unanswered and concealed. Observations obtained from infected primary neurons or mouse brain samples are more relevant to human clinical rabies than permissive cell lines; however, limitations regarding the ethical issue and sample accessibility become a hurdle for discovering new insights into virus-host interplays. To better understand RABV pathogenesis in humans, we generated human-induced pluripotent stem cell (hiPSC)-derived neurons to offer the opportunity for an inimitable study of RABV infection at a molecular level in a pathologically relevant cell type. This study describes the characteristics and detailed proteomic changes of hiPSC-derived neurons in response to RABV infection using LC-MS/MS quantitative analysis. Gene ontology (GO) enrichment of differentially expressed proteins (DEPs) reveals temporal changes of proteins related to metabolic process, immune response, neurotransmitter transport/synaptic vesicle cycle, cytoskeleton organization, and cell stress response, demonstrating fundamental underlying mechanisms of neuropathogenesis in a time-course dependence. Lastly, we highlighted plausible functions of heat shock cognate protein 70 (HSC70 or HSPA8) that might play a pivotal role in regulating RABV replication and pathogenesis. Our findings acquired from this hiPSC-derived neuron platform help to define novel cellular mechanisms during RABV infection, which could be applicable to further studies to widen views of RABV-host interaction.

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.

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.

Multi-scale comparative transcriptome analysis reveals key genes and metabolic reprogramming processes associated with oil palm fruit abscission

BACKGROUND

Fruit abscission depends on cell separation that occurs within specialized cell layers that constitute an abscission zone (AZ). To determine the mechanisms of fleshy fruit abscission of the monocot oil palm (Elaeis guineensis Jacq.) compared with other abscission systems, we performed multi-scale comparative transcriptome analyses on fruit targeting the developing primary AZ and adjacent tissues.

RESULTS

Combining between-tissue developmental comparisons with exogenous ethylene treatments, and naturally occurring abscission in the field, RNAseq analysis revealed a robust core set of 168 genes with differentially regulated expression, spatially associated with the ripe fruit AZ, and temporally restricted to the abscission timing. The expression of a set of candidate genes was validated by qRT-PCR in the fruit AZ of a natural oil palm variant with blocked fruit abscission, which provides evidence for their functions during abscission. Our results substantiate the conservation of gene function between dicot dry fruit dehiscence and monocot fleshy fruit abscission. The study also revealed major metabolic transitions occur in the AZ during abscission, including key senescence marker genes and transcriptional regulators, in addition to genes involved in nutrient recycling and reallocation, alternative routes for energy supply and adaptation to oxidative stress.

CONCLUSIONS

The study provides the first reference transcriptome of a monocot fleshy fruit abscission zone and provides insight into the mechanisms underlying abscission by identifying key genes with functional roles and processes, including metabolic transitions, cell wall modifications, signalling, stress adaptations and transcriptional regulation, that occur during ripe fruit abscission of the monocot oil palm. The transcriptome data comprises an original reference and resource useful towards understanding the evolutionary basis of this fundamental plant process.

The complete mitochondrial genome sequence of the mountain crab Indochinamon bhumibol

Indochinamon bhumibol has been found as the biggest freshwater crab in Thailand. In this study, we report the first complete sequence of mitochondrial genome from I. bhumibol encoding 13 protein-coding genes, 22 transfer RNAs, and 2 ribosomal RNAs. The nucleotide composition of I. bhumibol mitogenome showed a strong AT bias (70.4%) with a low GC content (29.6%). Comparative phylogenetic analysis with 28 crustaceans based on nine conserved genes demonstrated that I. bhumibol was closely related to members of the Potamidae family.

Candidate gene association of gene expression data in sugarcane contrasting for sucrose content

Association mapping of gene expression data, generated from transcriptome and proteome studies, provides a means of understanding the functional significance and trait association potential of candidate genes. In this study, we applied candidate gene association mapping to validate sugarcane genes, using data from the starch and sucrose metabolism pathway, transcriptome, and proteome. We performed multiplex PCR targeted amplicon sequencing of 109 candidate genes, using NGS technology. A range of statistical models, both single-locus and multi-locus, were compared for minimization of false positives in association mapping of four sugar-related traits with different heritability. The Fixed and random model Circulating Probability Unification model effectively suppressed false positives for both low- and high-heritability traits. We identified favorable alleles of the candidate genes involved in signalling and transcriptional regulation. The results will support genetic improvement of sugarcane and may help clarify the genetic architecture of sugar-related traits.

Assembly of the durian chloroplast genome using long PacBio reads

We have assembled the complete sequence of the Durio zibethinus chloroplast genome using long PacBio reads. Durian is a valuable commercial tree that produces durian fruit, which is popular in Southeast Asia. The chloroplast genome assembled into a single 143 kb cyclic contig that contained 111 genes. There were 46 short direct repeats (45 to 586 bp) and five short inverted repeats (63 to 169 bp). The long reads that were used for the assembly span the entire chloroplast with > 10 kb overlaps and multiple long reads join the start of the contig to the end of the contig. The durian chloroplast was found to lack the large inverted repeat that is common in chloroplast genomes. An additional 24 durian varieties were sequenced and compared to the assembly and found to also lack the large inverted repeat. There were nine SNPs among the varieties.

The Genome and Transcriptome Analysis of the Vigna mungo Chloroplast

Vigna mungo is cultivated in approximately 5 million hectares worldwide. The chloroplast genome of this species has not been previously reported. In this study, we sequenced the genome and transcriptome of the V. mungo chloroplast. We identified many positively selected genes in the photosynthetic pathway (e.g., rbcL, ndhF, and atpF) and RNA polymerase genes (e.g., rpoC2) from the comparison of the chloroplast genome of V. mungo, temperate legume species, and tropical legume species. Our transcriptome data from PacBio isoform sequencing showed that the 51-kb DNA inversion could affect the transcriptional regulation of accD polycistronic. Using Illumina deep RNA sequencing, we found RNA editing of clpP in the leaf, shoot, flower, fruit, and root tissues of V. mungo. We also found three G-to-A RNA editing events that change guanine to adenine in the transcripts transcribed from the adenine-rich regions of the ycf4 gene. The edited guanine bases were found particularly in the chloroplast genome of the Vigna species. These G-to-A RNA editing events were likely to provide a mechanism for correcting DNA base mutations. The V. mungo chloroplast genome sequence and the analysis results obtained in this study can apply to phylogenetic studies and chloroplast genome engineering.

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.

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.

Differential expression between drought-tolerant and drought-sensitive sugarcane under mild and moderate water stress as revealed by a comparative analysis of leaf transcriptome

Sugarcane contributes 80% of global sugar production and to bioethanol generation for the bioenergy industry. Its productivity is threatened by drought that can cause up to 60% yield loss. This study used RNA-Seq to gain a better understanding of the underlying mechanism by which drought-tolerant sugarcane copes with water stress. We compared gene expression in KPS01-12 (drought-tolerant genotype) and UT12 (drought-sensitive genotype) that have significantly different yield loss rates under drought conditions. We treated KPS01-12 and UT12 with mild and moderate water stress and found differentially expressed genes in various biological processes. KPS01-12 had higher expression of genes that were involved in water retention, antioxidant secondary metabolite biosynthesis, and oxidative and osmotic stress response than UT12. In contrast, the sensitive genotype had more down-regulated genes that were involved in photosynthesis, carbon fixation and Calvin cycle than the tolerant genotype. Our obtained expression profiles suggest that the tolerant sugarcane has a more effective genetic response than the sensitive genotype at the initiation of drought stress. The knowledge gained from this study may be applied in breeding programs to improve sugarcane production in drought conditions.

The complete mitochondrial genome of Luffa acutangula

Based on PacBio de novo assembly, we report the first complete mitochondrial genome of Luffa acutangula (460,333 bp) containing nine large chloroplast-derived sequences (1.9–17.3 kb) across the mitogenome. The base composition of the mitogenome in descending order is A: 28.02%, C: 22.04%, G: 21.83% and T: 28.10%, and the G + C content is 43.87%. There are 63 mitochondrial genes including 40 protein-coding genes, 3 rRNA genes and 20 tRNA genes. Additionally, a total of 288 repeats ranging from 31 to 5,301 bp were identified, accounting for 5.7% of the mitogenome. Two large direct repeats (5,301 and 405 bp) within the mitogenome were found for the formation of four subgenomic molecules. A phylogenetic analysis showed that L. acutangula was closely related to other species in Cucurbiaceae. This mitogenome provides useful genetic information for evolutionary studies.