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.

Rice height QTLs in KDML105 chromosome segment substitution lines.

Rice is an important crop that is consumed by approximately half of the world's population on a regular basis. Plant height is an important characteristic with shorter rice often having higher lodging resistance and better soil nutrient utilization allowing for lower fertilizer use. We used a Chromosome Segment Substitution Line (CSSL) population generated by introgressing segments of CT9993 and IR62266 into KDML 105. We identified height QTLs on chromosomes 1 and 4. We performed whole genome sequencing of the parental lines and found that IR62266 has the deletion in Gibberellin 20-oxidase 2 corresponding to the semi-dwarf 1 locus. However, short height on chromosome 1 came from CT9993 with no mutation in Gibberellin 20-oxidase 2, or any known height genes. The height QTL on chromosome 4 contains mutations in Peroxisome biogenesis protein 6, which has been linked to a reduced growth phenotype in A. thaliana, making this a good candidate height gene.

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.

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.

Clonal relationship of synchronous head and neck cancer and esophageal cancer assessed by single nucleotide polymorphism-based loss of heterozygosity analysis.

BACKGROUND: The prognoses of head and neck squamous cell carcinoma (HNSCC) and esophageal squamous cell carcinoma (ESCC) are poor, especially when both tumors occur at the same time. We examined the clonal relatedness of HNSCCs with synchronous ESCCs to confirm whether the second tumors were metastasis or separate second primary malignancies (SPMs) using loss of heterozygosity (LOH) analysis. METHODS: Twenty-one pairs of formalin-fixed paraffin-embedded tissue from HNSCC patients with synchronous esophageal cancer were analyzed by single nucleotide polymorphism (SNP) array using the Illumina HumanCytoSNP FFPE-12 BeadChip (San Diego, CA), which contains approximately 300,000 probes. LOH was identified using Nexus Copy Number software (El Segundo, CA). RESULTS: Comparing the LOH pattern between HNSCC and paired ESCC, we found that 20 out of 21 paired tissues had a high number of discordant LOHs (LOH identified solely in the primary HNSCC but not in synchronous ESCC at the same genomic location) and a low number of concordant LOHs (LOH at the same genomic location in both HNSCC and ESCC). Only one case fell into the undetermined category. Therefore, these 20 ESCCs were classified as SPMs or second field tumors (SFTs). Moreover, the HNSCC patients with molecularly confirmed esophageal SPM had significantly poorer survival than the other patients. CONCLUSIONS: We propose the use of a genome-wide SNP array as a tool to differentiate metastatic tumors from SPM/SFT. The SNP array offers genome-wide LOH information that earlier microsatellite analysis studies lack. The ability to accurately identify SPM should contribute to a better treatment plan and follow-up care of these patients.

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.

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.

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.

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.

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 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.

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.

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 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.

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.

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.

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

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

A 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 chromosome-scale reference genome assembly of yellow mangrove (Bruguiera parviflora) reveals a whole genome duplication event associated with the Rhizophoraceae lineage.

Mangrove forest ecosystems support a diverse flora and fauna of marine and terrestrial species and have important direct and indirect economic, ecological and social values to mankind. Yellow mangrove (Bruguiera parviflora) belongs to the Rhizophoraceae family and is widely distributed in the intertidal zones along sheltered coastal areas in tropical latitudes. Here, we present a high-quality, chromosome-level assembly of the B. parviflora genome. We employed the 10x Genomics linked-read technology to obtain a preliminary assembly, which was subsequently scaffolded using the long-range chromatin contact mapping technique (HiC) to obtain a final assembly containing 213,026,782 bases in 10,045 scaffolds with an N50 length of 10,906,948 bases. Our gene prediction recovered 96.5% of the highly conserved orthologues in the Embryophyta lineage based on the Benchmarking Universal Single-Copy Orthologues (BUSCO) analysis. We analysed the transversion rate at fourfold-degenerate sites from orthologous gene pairs and discovered evidence supporting a recent whole-genome duplication event in B. parviflora and other Rhizophoreae members. Comparative studies based on single-copy orthologous genes indicated that B. parviflora and Bruguiera gymnorrhiza diverged approximately 24.1 million years ago. The population structure analysis revealed that 63 B. parviflora accessions from different geographical regions in Thailand were an admixture of two subpopulations. The examination of alternative splicing events in B. parviflora showed that the most prevalent splicing mechanism was intron retention. This high-quality genome assembly together with the genetic diversity information obtained from the germplasm provide useful genomic resources for future studies on comparative phylogenetics and evolution of adaptive traits in mangrove species.