Comparative analysis of bacterial communities associated with healthy and diseased corals in the Indonesian sea

Coral reef ecosystems are impacted by climate change and human activities, such as increasing coastal development, overfishing, sewage and other pollutant discharge, and consequent eutrophication, which triggers increasing incidents of diseases and deterioration of corals worldwide. In this study, bacterial communities associated with four species of corals: Acropora aspera, Acropora formosa, Cyphastrea sp., and Isopora sp. in the healthy and disease stages with different diseases were compared using tagged 16S rRNA sequencing. In total, 59 bacterial phyla, 190 orders, and 307 genera were assigned in coral metagenomes where Proteobacteria and Firmicutes were predominated followed by Bacteroidetes together with Actinobacteria, Fusobacteria, and Lentisphaerae as minor taxa. Principal Coordinates Analysis (PCoA) showed separated clustering of bacterial diversity in healthy and infected groups for individual coral species. Fusibacter was found as the major bacterial genus across all corals. The lower number of Fusibacter was found in A. aspera infected with white band disease and Isopora sp. with white plaque disease, but marked increases of Vibrio and Acrobacter, respectively, were observed. This was in contrast to A. formosa infected by a black band and Cyphastrea sp. infected by yellow blotch diseases which showed an increasing abundance of Fusibacter but a decrease in WH1-8 bacteria. Overall, infection was shown to result in disturbance in the complexity and structure of the associated bacterial microbiomes which can be relevant to the pathogenicity of the microbes associated with infected corals

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

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

Transcriptome analysis of 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

Uncovering full-length transcript isoforms of sugarcane cultivar Khon Kaen 3 using single-molecule long-read sequencing

Background Sugarcane is an important global food crop and energy resource. To facilitate the sugarcane improvement program, genome and gene information are important for studying traits at the molecular level. Most currently available transcriptome data for sugarcane were generated using second-generation sequencing platforms, which provide short reads. The de novo assembled transcripts from these data are limited in length, and hence may be incomplete and inaccurate, especially for long RNAs. Methods We generated a transcriptome dataset of leaf tissue from a commercial Thai sugarcane cultivar Khon Kaen 3 (KK3) using PacBio RS II single-molecule long-read sequencing by the Iso-Seq method. Short-read RNA-Seq data were generated from the same RNA sample using the Ion Proton platform for reducing base calling errors. Results A total of 119,339 error-corrected transcripts were generated with the N50 length of 3,611 bp, which is on average longer than any previously reported sugarcane transcriptome dataset. 110,253 sequences (92.4%) contain an open reading frame (ORF) of at least 300 bp long with ORF N50 of 1,416 bp. The mean lengths of 5′ and 3′ untranslated regions in 73,795 sequences with complete ORFs are 1,249 and 1,187 bp, respectively. 4,774 transcripts are putatively novel full-length transcripts which do not match with a previous Iso-Seq study of sugarcane. We annotated the functions of 68,962 putative full-length transcripts with at least 90% coverage when compared with homologous protein coding sequences in other plants. Discussion The new catalog of transcripts will be useful for genome annotation, identification of splicing variants, SNP identification, and other research pertaining to the sugarcane improvement program. The putatively novel transcripts suggest unique features of KK3, although more data from different tissues and stages of development are needed to establish a reference transcriptome of this cultivar

A genome scan for quantitative trait loci affecting cyanogenic potential of cassava root in an outbred population

<p>Abstract</p> <p>Background</p> <p>Cassava (<it>Manihot esculenta </it>Crantz) can produce cyanide, a toxic compound, without self-injury. That ability was called the cyanogenic potential (CN). This project aimed to identify quantitative trait loci (QTL) associated with the CN in an outbred population derived from 'Hanatee' × 'Huay Bong 60', two contrasting cultivars. CN was evaluated in 2008 and in 2009 at Rayong province, and in 2009 at Lop Buri province, Thailand. CN was measured using a picrate paper kit. QTL analysis affecting CN was performed with 303 SSR markers.</p> <p>Results</p> <p>The phenotypic values showed continuous variation with transgressive segregation events with more (115 ppm) and less CN (15 ppm) than either parent ('Hanatee' had 33 ppm and 'Huay Bong 60' had 95 ppm). The linkage map consisted of 303 SSR markers, on 27 linkage groups with a map that encompassed 1,328 cM. The average marker interval was 5.8 cM. Five QTL underlying CN were detected. <it>CN08R1</it>from 2008 at Rayong, <it>CN09R1</it>and <it>CN09R2 </it>from 2009 at Rayong, and <it>CN09L1 </it>and <it>CN09L2 </it>from 2009 at Lop Buri were mapped on linkage group 2, 5, 10 and 11, respectively. Among all the identified QTL, <it>CN09R1 </it>was the most significantly associated with the CN trait with LOD score 5.75 and explained the greatest percentage of phenotypic variation (%Expl.) of 26%.</p> <p>Conclusions</p> <p>Five new QTL affecting CN were successfully identified from 4 linkage groups. Discovery of these QTL can provide useful markers to assist in cassava breeding and studying genes affecting the trait.</p

Pins Gene Table v2.0: An Online Genome Database of 37 <i>Pythium insidiosum</i> 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

Microbial communities in the reef water at Kham Island, lower Gulf of Thailand

Coral reefs are among the most biodiverse habitats on Earth, but knowledge of their associated marinemicrobiome remains limited. To increase the understanding of the coral reef ecosystem in the lower Gulf of Thailand, this study utilized 16S and 18S rRNA gene-based pyrosequencing to identify the prokaryotic and eukaryotic microbiota present in the reef water at Kham Island, Trat province, Thailand (N6.97 E100.86). The obtained result was then compared with the published microbiota from different coral reef water and marine sites. The coral reefs at Kham Island are of the fringe type. The reefs remain preserved and abundant. The community similarity indices (i.e., Lennon similarity index, Yue & Clayton similarity index) indicated that the prokaryotic composition of Kham was closely related to that of Kra, another fringing reef site in the lower Gulf of Thailand, followed by coral reef water microbiota at GS048b (Cooks Bay, Fr. Polynesia), Palmyra (Northern Line Islands, United States) and GS108b (Coccos Keeling, Australia), respectively. Additionally, the microbial eukaryotic populations at Kham was analyzed and compared with the available database at Kra. Both eukaryotic microbiota, in summer and winter seasons, were correlated. An abundance of Dinophysis acuminata was noted in the summer season, in accordance with its reported cause of diarrhoeatic shellfish outbreak in the summer season elsewhere. The slightly lower biodiversity in Kham than at Kra might reflect the partly habitat difference due to coastal anthropogenic activities and minor water circulation, as Kham locates close to the mainland and is surrounded by islands (e.g., Chang and Kut islands). The global marine microbiota comparison suggested relatively similar microbial structures among coral sites irrespective of geographical location, supporting the importance of coral-associated marine microbiomes, and Spearman’s correlation analysis between community membership and factors of shore distance and seawater temperature indicated potential correlation of these factors (p-values < 0.05) with Kham, Kra, and some other coral and coastal sites. Together, this study provided the second marine microbial database for the coral reef of the lower Gulf of Thailand, and a comparison of the coral-associated marine microbial diversity among global ocean sites