Chon-Kit Kenneth Chan

Analysis of RNA-Seq Data Using TopHat and Cufflinks

The recent advances in high throughput RNA sequencing (RNA-Seq) have generated huge amounts of data in a very short span of time for a single sample. These data have required the parallel advancement of computing tools to organize and interpret them meaningfully in terms of biological implications, at the same time using minimum computing resources to reduce computation costs. Here we describe the method of analyzing RNA-seq data using the set of open source software programs of the Tuxedo suite: TopHat and Cufflinks. TopHat is designed to align RNA-seq reads to a reference genome, while Cufflinks assembles these mapped reads into possible transcripts and then generates a final transcriptome assembly. Cufflinks also includes Cuffdiff, which accepts the reads assembled from two or more biological conditions and analyzes their differential expression of genes and transcripts, thus aiding in the investigation of their transcriptional and post transcriptional regulation under different conditions. We also describe the use of an accessory tool called CummeRbund, which processes the output files of Cuffdiff and gives an output of publication quality plots and figures of the users choice. We demonstrate the effectiveness of the Tuxedo suite by analyzing RNA-Seq datasets of Arabidopsis thaliana root subjected to two different conditions.

Transcriptomic changes in coral holobionts provide insights into physiological challenges of future climate and ocean change

Tropical reef-building coral stress levels will intensify with the predicted rising atmospheric CO2 resulting in ocean temperature and acidification increase. Most studies to date have focused on the destabilization of coral-dinoflagellate symbioses due to warming oceans, or declining calcification due to ocean acidification. In our study, pH and temperature conditions consistent with the end-of-century scenarios of the Intergovernmental Panel on Climate Change (IPCC) caused major changes in photosynthesis and respiration, in addition to decreased calcification rates in the coral Acropora millepora. Population density of symbiotic dinoflagellates (Symbiodinium) under high levels of ocean acidification and temperature (Representative Concentration Pathway, RCP8.5) decreased to half of that found under present day conditions, with photosynthetic and respiratory rates also being reduced by 40%. These physiological changes were accompanied by evidence for gene regulation of calcium and bicarbonate transporters along with components of the organic matrix. Metatranscriptomic RNA-Seq data analyses showed an overall down regulation of metabolic transcripts, and an increased abundance of transcripts involved in circadian clock control, controlling the damage of oxidative stress, calcium signaling/homeostasis, cytoskeletal interactions, transcription regulation, DNA repair, Wnt signaling and apoptosis/immunity/ toxins. We suggest that increased maintenance costs under ocean acidification and warming, and diversion of cellular ATP to pH homeostasis, oxidative stress response, UPR and DNA repair, along with metabolic suppression, may underpin why Acroporid species tend not to thrive under future environmental stress. Our study highlights the potential increased energy demand when the coral holobiont is exposed to high levels of ocean warming and acidification.

A chromosomal genomics approach to assess and validate the desi and kabuli draft chickpea genome assemblies

With the expansion of next-generation sequencing technology and advanced bioinformatics, there has been a rapid growth of genome sequencing projects. However, while this technology enables the rapid and cost-effective assembly of draft genomes, the quality of these assemblies usually falls short of gold standard genome assemblies produced using the more traditional BAC by BAC and Sanger sequencing approaches. Assembly validation is often performed by the physical anchoring of genetically mapped markers, but this is prone to errors and the resolution is usually low, especially towards centromeric regions where recombination is limited. New approaches are required to validate reference genome assemblies. The ability to isolate individual chromosomes combined with next-generation sequencing permits the validation of genome assemblies at the chromosome level. We demonstrate this approach by the assessment of the recently published chickpea kabuli and desi genomes. While previous genetic analysis suggests that these genomes should be very similar, a comparison of their chromosome sizes and published assemblies highlights significant differences. Our chromosomal genomics analysis highlights short defined regions that appear to have been misassembled in the kabuli genome and identifies large-scale misassembly in the draft desi genome. The integration of chromosomal genomics tools within genome sequencing projects has the potential to significantly improve the construction and validation of genome assemblies. The approach could be applied both for new genome assemblies as well as published assemblies, and complements currently applied genome assembly strategies.

Unfolding the secrets of coral–algal symbiosis

Dinoflagellates from the genus Symbiodinium form a mutualistic symbiotic relationship with reef-building corals. Here we applied massively parallel Illumina sequencing to assess genetic similarity and diversity among four phylogenetically diverse dinoflagellate clades (A, B, C and D) that are commonly associated with corals. We obtained more than 30 000 predicted genes for each Symbiodinium clade, with a majority of the aligned transcripts corresponding to sequence data sets of symbiotic dinoflagellates and <2% of sequences having bacterial or other foreign origin. We report 1053 genes, orthologous among four Symbiodinium clades, that share a high level of sequence identity to known proteins from the SwissProt (SP) database. Approximately 80% of the transcripts aligning to the 1053 SP genes were unique to Symbiodinium species and did not align to other dinoflagellates and unrelated eukaryotic transcriptomes/genomes. Six pathways were common to all four Symbiodinium clades including the phosphatidylinositol signaling system and inositol phosphate metabolism pathways. The list of Symbiodinium transcripts common to all four clades included conserved genes such as heat shock proteins (Hsp70 and Hsp90), calmodulin, actin and tubulin, several ribosomal, photosynthetic and cytochrome genes and chloroplast-based heme-containing cytochrome P450, involved in the biosynthesis of xanthophylls. Antioxidant genes, which are important in stress responses, were also preserved, as were a number of calcium-dependent and calcium/calmodulin-dependent protein kinases that may play a role in the establishment of symbiosis. Our findings disclose new knowledge about the genetic uniqueness of symbiotic dinoflagellates and provide a list of homologous genes important for the foundation of coral–algal symbiosis.

Genome-wide survey of the seagrass Zostera muelleri suggests modification of the ethylene signalling network

Highlight An aquatic life genome sequencing suggests a complete loss of genes for ethylene biosynthesis and signalling pathways in the seagrasses, Zostera muelleri and Zostera marina, a new model for hormone studies.

High-resolution skim genotyping by sequencing reveals the distribution of crossovers and gene conversions in Cicer arietinum and Brassica napus

Key messageWe characterise the distribution of crossover and non-crossover recombination inBrassica napusandCicer arietinumusing a low-coverage genotyping by sequencing pipeline SkimGBS.AbstractThe growth of next-generation DNA sequencing technologies has led to a rapid increase in sequence-based genotyping for applications including diversity assessment, genome structure validation and gene–trait association. We have established a skim-based genotyping by sequencing method for crop plants and applied this approach to genotype-segregating populations of Brassica napus and Cicer arietinum. Comparison of progeny genotypes with those of the parental individuals allowed the identification of crossover and non-crossover (gene conversion) events. Our results identify the positions of recombination events with high resolution, permitting the mapping and frequency assessment of recombination in segregating populations.

Identification and characterization of more than 4 million intervarietal SNPs across the group 7 chromosomes of bread wheat

Despite being a major international crop, our understanding of the wheat genome is relatively poor due to its large size and complexity. To gain a greater understanding of wheat genome diversity, we have identified single nucleotide polymorphisms between 16 Australian bread wheat varieties. Whole-genome shotgun Illumina paired read sequence data were mapped to the draft assemblies of chromosomes 7A, 7B and 7D to identify more than 4 million intervarietal SNPs. SNP density varied between the three genomes, with much greater density observed on the A and B genomes than the D genome. This variation may be a result of substantial gene flow from the tetraploid Triticum turgidum, which possesses A and B genomes, during early co-cultivation of tetraploid and hexaploid wheat. In addition, we examined SNP density variation along the chromosome syntenic builds and identified genes in low-density regions which may have been selected during domestication and breeding. This study highlights the impact of evolution and breeding on the bread wheat genome and provides a substantial resource for trait association and crop improvement. All SNP data are publically available on a generic genome browser GBrowse at www.wheatgenome.info.

Assembly and comparison of two closely related Brassica napus genomes

Summary As an increasing number of plant genome sequences become available, it is clear that gene content varies between individuals, and the challenge arises to predict the gene content of a species. However, genome comparison is often confounded by variation in assembly and annotation. Differentiating between true gene absence and variation in assembly or annotation is essential for the accurate identification of conserved and variable genes in a species. Here, we present the de novo assembly of the B. napus cultivar Tapidor and comparison with an improved assembly of the Brassica napus cultivar Darmor‐bzh. Both cultivars were annotated using the same method to allow comparison of gene content. We identified genes unique to each cultivar and differentiate these from artefacts due to variation in the assembly and annotation. We demonstrate that using a common annotation pipeline can result in different gene predictions, even for closely related cultivars, and repeat regions which collapse during assembly impact whole genome comparison. After accounting for differences in assembly and annotation, we demonstrate that the genome of Darmor‐bzh contains a greater number of genes than the genome of Tapidor. Our results are the first step towards comparison of the true differences between B. napus genomes and highlight the potential sources of error in future production of a B. napus pangenome.

The Genome of a Southern Hemisphere Seagrass Species ( Zostera muelleri )

Adaptation of seagrass to the marine environment involves modification and loss of conserved plant genes. Seagrasses are marine angiosperms that evolved from land plants but returned to the sea around 140 million years ago during the early evolution of monocotyledonous plants. They successfully adapted to abiotic stresses associated with growth in the marine environment, and today, seagrasses are distributed in coastal waters worldwide. Seagrass meadows are an important oceanic carbon sink and provide food and breeding grounds for diverse marine species. Here, we report the assembly and characterization of the Zostera muelleri genome, a southern hemisphere temperate species. Multiple genes were lost or modified in Z. muelleri compared with terrestrial or floating aquatic plants that are associated with their adaptation to life in the ocean. These include genes for hormone biosynthesis and signaling and cell wall catabolism. There is evidence of whole-genome duplication in Z. muelleri; however, an ancient pan-commelinid duplication event is absent, highlighting the early divergence of this species from the main monocot lineages.

Homoeologous exchange is a major cause of gene presence/absence variation in the amphidiploid Brassica napus

Summary Homoeologous exchanges (HEs) have been shown to generate novel gene combinations and phenotypes in a range of polyploid species. Gene presence/absence variation (PAV) is also a major contributor to genetic diversity. In this study, we show that there is an association between these two events, particularly in recent Brassica napus synthetic accessions, and that these represent a novel source of genetic diversity, which can be captured for the improvement of this important crop species. By assembling the pangenome of B. napus, we show that 38% of the genes display PAV behaviour, with some of these variable genes predicted to be involved in important agronomic traits including flowering time, disease resistance, acyl lipid metabolism and glucosinolate metabolism. This study is a first and provides a detailed characterization of the association between HEs and PAVs in B. napus at the pangenome level.