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.

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.

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.

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

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.

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.

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.

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.

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.

Heat-induced shift in coral microbiome reveals several members of the Rhodobacteraceae family as indicator species for thermal stress in Porites lutea

The coral holobiont is a complex ecosystem consisting of coral animals and a highly diverse consortium of associated microorganisms including algae, fungi, and bacteria. Several studies have highlighted the importance of coral‐associated bacteria and their potential roles in promoting the host fitness and survival. Recently, dynamics of coral‐associated microbiomes have been demonstrated to be linked to patterns of coral heat tolerance. Here, we examined the effect of elevated seawater temperature on the structure and diversity of bacterial populations associated with Porites lutea, using full‐length 16S rRNA sequences obtained from Pacific Biosciences circular consensus sequencing. We observed a significant increase in alpha diversity indices and a distinct shift in microbiome composition during thermal stress. There was a marked decline in the apparent relative abundance of Gammaproteobacteria family Endozoicomonadaceae after P. lutea had been exposed to elevated seawater temperature. Concomitantly, the bacterial community structure shifted toward the predominance of Alphaproteobacteria family Rhodobacteraceae. Interestingly, we did not observe an increase in relative abundance of Vibrio‐related sequences in our heat‐stressed samples even though the appearance of Vibrio spp. has often been detected in parallel with the increase in the relative abundance of Rhodobacteraceae during thermal bleaching in other coral species. The ability of full‐length 16S rRNA sequences in resolving taxonomic uncertainty of associated bacteria at a species level enabled us to identify 24 robust indicator bacterial species for thermally stressed corals. It is worth noting that the majority of those indicator species were members of the family Rhodobacteraceae. The comparison of bacterial community structure and diversity between corals in ambient water temperature and thermally stressed corals may provide a better understanding on how bacteria symbionts contribute to the resilience of their coral hosts to ocean warming.

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.

Transcriptome analysis of oil palm inflorescences revealed candidate genes for an auxin signaling pathway involved in parthenocarpy

Oil palm parthenocarpic fruits, which are produced without fertilization, can be targeted to increase oil content because the majority of the fruit is occupied by mesocarp, the part in which palm oil is stored. Consequently, gaining an understanding of the parthenocarpic mechanism would be instrumental for producing parthenocarpic oil palm. This study aims to determine effects of auxin treatment and analyze differentially expressed genes in oil palm pistils at the pollination/anthesis stage, using an RNA sequencing (RNA seq) approach. The auxin treatment caused 100% parthenocarpy when auxin was sprayed before stigmas opened. The parthenocarpy decreased to 55%, 8% and 5% when the auxin was sprayed 1, 2 and 3 days after the opening of stigmas, respectively. Oil palm plants used for RNA seq were plants untreated with auxin as controls and auxin-treated plants on the day before pollination and 1 day after pollination. The number of raw reads ranged from 8,425,859 to 11,811,166 reads, with an average size ranging from 99 to 137 base pairs (bp). When compared with the oil palm transcriptome, the mapped reads ranged from 8,179,948 to 11,320,799 reads, representing 95.85–98.01% of the oil palm matching. Based on five comparisons between RNA seq of treatments and controls, and confirmation using reverse transcription polymerase chain reaction and quantitative real-time RT-PCR expression, five candidate genes, including probable indole-3-acetic acid (IAA)-amido synthetase GH3.8 (EgGH3.8), IAA-amido synthetase GH3.1 (EgGH3.1), IAA induced ARG7 like (EgARG7), tryptophan amino transferase-related protein 3-like (EgTAA3) and flavin-containing monooxygenase 1 (EgFMO1), were differentially expressed between auxin-treated and untreated samples. This evidence suggests a pathway of parthenocarpic fruit development at the beginning of fruit development. However, more research is needed to identify which genes are definitely involved in parthenocarpy.

Dynamics of coral-associated microbiomes during a thermal bleaching event

Coral‐associated microorganisms play an important role in their host fitness and survival. A number of studies have demonstrated connections between thermal tolerance in corals and the type/relative abundance of Symbiodinium they harbor. More recently, the shifts in coral‐associated bacterial profiles were also shown to be linked to the patterns of coral heat tolerance. Here, we investigated the dynamics of Porites lutea‐associated bacterial and algal communities throughout a natural bleaching event, using full‐length 16S rRNA and internal transcribed spacer sequences (ITS) obtained from PacBio circular consensus sequencing. We provided evidence of significant changes in the structure and diversity of coral‐associated microbiomes during thermal stress. The balance of the symbiosis shifted from a predominant association between corals and Gammaproteobacteria to a predominance of Alphaproteobacteria and to a lesser extent Betaproteobacteria following the bleaching event. On the contrary, the composition and diversity of Symbiodinium communities remained unaltered throughout the bleaching event. It appears that the switching and/or shuffling of Symbiodinium types may not be the primary mechanism used by P. lutea to cope with increasing seawater temperature. The shifts in the structure and diversity of associated bacterial communities may contribute more to the survival of the coral holobiont under heat stress.

The two chromosomes of the mitochondrial genome of a sugarcane cultivar: assembly and recombination analysis using long PacBio reads

Sugarcane accounts for a large portion of the worlds sugar production. Modern commercial cultivars are complex hybrids of S. officinarum and several other Saccharum species. Historical records identify New Guinea as the origin of S. officinarum and that a small number of plants originating from there were used to generate all modern commercial cultivars. The mitochondrial genome can be a useful way to identify the maternal origin of commercial cultivars. We have used the PacBio RSII to sequence and assemble the mitochondrial genome of a South East Asian commercial cultivar, known as Khon Kaen 3. The long read length of this sequencing technology allowed for the mitochondrial genome to be assembled into two distinct circular chromosomes with all repeat sequences spanned by individual reads. Comparison of five commercial hybrids, two S. officinarum and one S. spontaneum to our assembly reveals no structural rearrangements between our assembly, the commercial hybrids and an S. officinarum from New Guinea. The S. spontaneum, from India, and one sample of S. officinarum (unknown origin) are substantially rearranged and have a large number of homozygous variants. This supports the record that S. officinarum plants from New Guinea are the maternal source of all modern commercial hybrids.

Effects of methylation-sensitive enzymes on the enrichment of genic SNPs and the degree of genome complexity reduction in a two-enzyme genotyping-by-sequencing (GBS) approach: a case study in oil palm (Elaeis guineensis)

Advances in next generation sequencing have facilitated a large-scale single nucleotide polymorphism (SNP) discovery in many crop species. Genotyping-by-sequencing (GBS) approach couples next generation sequencing with genome complexity reduction techniques to simultaneously identify and genotype SNPs. Choice of enzymes used in GBS library preparation depends on several factors including the number of markers required, the desired level of multiplexing, and whether the enrichment of genic SNP is preferred. We evaluated various combinations of methylation-sensitive (AatII, PstI, MspI) and methylation-insensitive (SphI, MseI) enzymes for their effectiveness in genome complexity reduction and enrichment of genic SNPs. We discovered that the use of two methylation-sensitive enzymes effectively reduced genome complexity and did not require a size selection step. On the contrary, the genome coverage of libraries constructed with methylation-insensitive enzymes was quite high, and the additional size selection step may be required to increase the overall read depth. We also demonstrated the effectiveness of methylation-sensitive enzymes in enriching for SNPs located in genic regions. When two methylation-insensitive enzymes were used, only 16% of SNPs identified were located in genes and 18% in the vicinity (± 5 kb) of the genic regions, while most SNPs resided in the intergenic regions. In contrast, a remarkable degree of enrichment was observed when two methylation-sensitive enzymes were employed. Almost two thirds of the SNPs were located either inside (32-36%) or in the vicinity (28-31%) of the genic regions. These results provide useful information to help researchers choose appropriate GBS enzymes in oil palm and other crop species.

SNP identification from RNA sequencing and linkage map construction of rubber tree for anchoring the draft genome

Hevea brasiliensis, or rubber tree, is an important crop species that accounts for the majority of natural latex production. The rubber tree nuclear genome consists of 18 chromosomes and is roughly 2.15 Gb. The current rubber tree reference genome assembly consists of 1,150,326 scaffolds ranging from 200 to 531,465 bp and totalling 1.1 Gb. Only 143 scaffolds, totalling 7.6 Mb, have been placed into linkage groups. We have performed RNA-seq on 6 varieties of rubber tree to identify SNPs and InDels and used this information to perform target sequence enrichment and high throughput sequencing to genotype a set of SNPs in 149 rubber tree offspring from a cross between RRIM 600 and RRII 105 rubber tree varieties. We used this information to generate a linkage map allowing for the anchoring of 24,424 contigs from 3,009 scaffolds, totalling 115 Mb or 10.4% of the published sequence, into 18 linkage groups. Each linkage group contains between 319 and 1367 SNPs, or 60 to 194 non-redundant marker positions, and ranges from 156 to 336 cM in length. This linkage map includes 20,143 of the 69,300 predicted genes from rubber tree and will be useful for mapping studies and improving the reference genome assembly.

Genome-wide SNP discovery and identification of QTL associated with agronomic traits in oil palm using genotyping-by-sequencing (GBS)

Oil palm has become one of the most important oil crops in the world. Marker-assisted selections have played a pivotal role in oil palm breeding programs. Here, we report the use of genotyping-by-sequencing (GBS) approach for a large-scale SNP discovery and genotyping of a mapping population. Reduced representation libraries of 108 F2 progeny were sequenced and a total of 524 million reads were obtained. We detected 21,471 single nucleotide substitutions, most of which (62.6%) represented transition events. Of 3417 fully informative SNP markers, we were able to place 1085 on a linkage map, which spanned 1429.6 cM and had an average of one marker every 1.26 cM. Three QTL affecting trunk height were detected on LG 10, 14 and 15, whereas a single QTL associated with fruit bunch weight was identified on LG 3. The use of GBS approach proved to be rapid, cost-effective and highly reproducible in this species.

Construction of a high-density integrated genetic linkage map of rubber tree (Hevea brasiliensis) using genotyping-by-sequencing (GBS)

Construction of linkage maps is crucial for genetic studies and marker-assisted breeding programs. Recent advances in next generation sequencing technologies allow for the generation of high-density linkage maps, especially in non-model species lacking extensive genomic resources. Here, we constructed a high-density integrated genetic linkage map of rubber tree (Hevea brasiliensis), the sole commercial producer of high-quality natural rubber. We applied a genotyping-by-sequencing (GBS) technique to simultaneously discover and genotype single nucleotide polymorphism (SNP) markers in two rubber tree populations. A total of 21,353 single nucleotide substitutions were identified, 55% of which represented transition events. GBS-based genetic maps of populations P and C comprised 1704 and 1719 markers and encompassed 2041 cM and 1874 cM, respectively. The average marker densities of these two maps were one SNP in 1.23-1.25 cM. A total of 1114 shared SNP markers were used to merge the two component maps. An integrated linkage map consisted of 2321 markers and spanned the cumulative length of 2052 cM. The composite map showed a substantial improvement in marker density, with one SNP marker in every 0.89 cM. To our knowledge, this is the most saturated genetic map in rubber tree to date. This integrated map allowed us to anchor 28,965 contigs, covering 135 Mb or 12% of the published rubber tree genome. We demonstrated that GBS is a robust and cost-effective approach for generating a common set of genome-wide SNP data suitable for constructing integrated linkage maps from multiple populations in a highly heterozygous agricultural species.

Large-scale SNP discovery through RNA sequencing and SNP genotyping by targeted enrichment sequencing in cassava (Manihot esculenta Crantz)

Cassava (Manihot esculenta Crantz) is one of the most important crop species being the main source of dietary energy in several countries. Marker-assisted selection has become an essential tool in plant breeding. Single nucleotide polymorphism (SNP) discovery via transcriptome sequencing is an attractive strategy for genome complexity reduction in organisms with large genomes. We sequenced the transcriptome of 16 cassava accessions using the Illumina HiSeq platform and identified 675,559 EST-derived SNP markers. A subset of those markers was subsequently genotyped by capture-based targeted enrichment sequencing in 100 F1 progeny segregating for starch viscosity phenotypes. A total of 2,110 non-redundant SNP markers were used to construct a genetic map. This map encompasses 1,785 cM and consists of 19 linkage groups. A major quantitative trait locus (QTL) controlling starch pasting properties was identified and shown to coincide with the QTL previously reported for this trait. With a high-density SNP-based linkage map presented here, we also uncovered a novel QTL associated with starch pasting time on LG 10.

Characterization of the chloroplast genome sequence of oil palm (Elaeis guineensis Jacq.)

Oil palm (Elaeis guineensis Jacq.) is an economically important crop, which is grown for oil production. To better understand the molecular basis of oil palm chloroplasts, we characterized the complete chloroplast (cp) genome sequence obtained from 454 pyrosequencing. The oil palm cp genome is 156,973 bp in length consisting of a large single-copy region of 85,192 bp flanked on each side by inverted repeats of 27,071 bp with a small single-copy region of 17,639 bp joining the repeats. The genome contains 112 unique genes: 79 protein-coding genes, 4 ribosomal RNA genes and 29 tRNA genes. By aligning the cp genome sequence with oil palm cDNA sequences, we observed 18 non-silent and 10 silent RNA editing events among 19 cp protein-coding genes. Creation of an initiation codon by RNA editing in rpl2 has been reported in several monocots and was also found in the oil palm cp genome. Fifty common chloroplast protein-coding genes from 33 plant taxa were used to construct ML and MP phylogenetic trees. Their topologies are similar and strongly support for the position of E. guineensis as the sister of closely related species Phoenix dactylifera in Arecaceae (palm families) of monocot subtrees.

Transcriptome sequencing of Hevea brasiliensis for development of microsatellite markers and construction of a genetic linkage map

To obtain more information on the Hevea brasiliensis genome, we sequenced the transcriptome from the vegetative shoot apex yielding 2 311 497 reads. Clustering and assembly of the reads produced a total of 113 313 unique sequences, comprising 28 387 isotigs and 84 926 singletons. Also, 17 819 expressed sequence tag (EST)-simple sequence repeats (SSRs) were identified from the data set. To demonstrate the use of this EST resource for marker development, primers were designed for 430 of the EST-SSRs. Three hundred and twenty-three primer pairs were amplifiable in H. brasiliensis clones. Polymorphic information content values of selected 47 SSRs among 20 H. brasiliensis clones ranged from 0.13 to 0.71, with an average of 0.51. A dendrogram of genetic similarities between the 20 H. brasiliensis clones using these 47 EST-SSRs suggested two distinct groups that correlated well with clone pedigree. These novel EST-SSRs together with the published SSRs were used for the construction of an integrated parental linkage map of H. brasiliensis based on 81 lines of an F1 mapping population. The map consisted of 97 loci, consisting of 37 novel EST-SSRs and 60 published SSRs, distributed on 23 linkage groups and covered 842.9 cM with a mean interval of 11.9 cM and ∼4 loci per linkage group. Although the numbers of linkage groups exceed the haploid number (18), but with several common markers between homologous linkage groups with the previous map indicated that the F1 map in this study is appropriate for further study in marker-assisted selection.

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

PREMISE OF THE STUDY

We demonstrated the application of high-throughput 454 sequencing technology in the identification of single nucleotide polymorphism (SNP) markers in Hevea brasiliensis.

METHODS AND RESULTS

A total of 5883 putative SNP positions were discovered in silico, and 10 biallelic SNP markers were validated from 454-derived EST sequences. The polymorphism information content (PIC) and the observed heterozygosity (H(o)) ranged from 0.0963-0.5135 and 0.1071-0.4643, respectively.

CONCLUSIONS

These markers can be useful for the construction of genetic maps, the identification of quantitative trait loci linked to commercially desirable traits, and the study of genetic structure in H. brasiliensis.

Characterization of the complete chloroplast genome of Hevea brasiliensis reveals genome rearrangement, RNA editing sites and phylogenetic relationships

Rubber tree (Hevea brasiliensis) is an economical plant and widely grown for natural rubber production. However, genomic research of rubber tree has lagged behind other species in the Euphorbiaceae family. We report the complete chloroplast genome sequence of rubber tree as being 161,191 bp in length including a pair of inverted repeats of 26,810 bp separated by a small single copy region of 18,362 bp and a large single copy region of 89,209 bp. The chloroplast genome contains 112 unique genes, 16 of which are duplicated in the inverted repeat. Of the 112 unique genes, 78 are predicted protein-coding genes, 4 are ribosomal RNA genes and 30 are tRNA genes. Relative to other plant chloroplast genomes, we observed a unique rearrangement in the rubber tree chloroplast genome: a 30-kb inversion between the trnE(UUC)-trnS(GCU) and the trnT(GGU)-trnR(UCU). A comparison between the rubber tree chloroplast genes and cDNA sequences revealed 51 RNA editing sites in which most (48 sites) were located in 26 protein coding genes and the other 3 sites were in introns. Phylogenetic analysis based on chloroplast genes demonstrated a close relationship between Hevea and Manihot in Euphorbiaceae and provided a strong support for a monophyletic group of the eurosid I.

The chloroplast genome sequence of mungbean (Vigna radiata) determined by high-throughput pyrosequencing: structural organization and phylogenetic relationships

Mungbean is an economically important crop which is grown principally for its protein-rich dry seeds. However, genomic research of mungbean has lagged behind other species in the Fabaceae family. Here, we reported the complete chloroplast (cp) genome sequence of mungbean obtained by the 454 pyrosequencing technology. The mungbean cp genome is 151 271 bp in length which includes a pair of inverted repeats (IRs) of 26 474 bp separated by a small single-copy region of 17 427 bp and a large single-copy region of 80 896 bp. The genome contains 108 unique genes and 19 of these genes are duplicated in the IR. Of these, 75 are predicted protein-coding genes, 4 ribosomal RNA genes and 29 tRNA genes. Relative to other plant cp genomes, we observed two distinct rearrangements: a 50-kb inversion between accD/rps16 and rbcL/trnK-UUU, and a 78-kb rearrangement between trnH/rpl14 and rps19/rps8. We detected sequence length polymorphism in the cp homopolymeric regions at the intra- and inter-specific levels in the Vigna species. Phylogenetic analysis demonstrated a close relationship between Vigna and Phaseolus in the phaseolinae subtribe and provided a strong support for a monophyletic group of the eurosid I.