Next-generation sequencing and large genome assemblies

Cardiovascular genetics: technological advancements and applicability for dilated cardiomyopathy

Genetics plays an important role in the pathophysiology of cardiovascular diseases, and is increasingly being integrated into clinical practice. Since 2008, both capacity and cost-efficiency of mutation screening of DNA have been increased magnificently due to the technological advancement obtained by next-generation sequencing. Hence, the discovery rate of genetic defects in cardiovascular genetics has grown rapidly and the financial threshold for gene diagnostics has been lowered, making large-scale DNA sequencing broadly accessible. In this review, the genetic variants, mutations and inheritance models are briefly introduced, after which an overview is provided of current clinical and technological applications in gene diagnostics and research for cardiovascular disease and in particular, dilated cardiomyopathy. Finally, a reflection on the future perspectives in cardiogenetics is given.

Challenges and opportunities of airborne metagenomics

Recent metagenomic studies of environments, such as marine and soil, have significantly enhanced our understanding of the diverse microbial communities living in these habitats and their essential roles in sustaining vast ecosystems. The increase in the number of publications related to soil and marine metagenomics is in sharp contrast to those of air, yet airborne microbes are thought to have significant impacts on many aspects of our lives from their potential roles in atmospheric events such as cloud formation, precipitation, and atmospheric chemistry to their major impact on human health. In this review, we will discuss the current progress in airborne metagenomics, with a special focus on exploring the challenges and opportunities of undertaking such studies. The main challenges of conducting metagenomic studies of airborne microbes are as follows: 1) Low density of microorganisms in the air, 2) efficient retrieval of microorganisms from the air, 3) variability in airborne microbial community composition, 4) the lack of standardized protocols and methodologies, and 5) DNA sequencing and bioinformatics-related challenges. Overcoming these challenges could provide the groundwork for comprehensive analysis of airborne microbes and their potential impact on the atmosphere, global climate, and our health. Metagenomic studies offer a unique opportunity to examine viral and bacterial diversity in the air and monitor their spread locally or across the globe, including threats from pathogenic microorganisms. Airborne metagenomic studies could also lead to discoveries of novel genes and metabolic pathways relevant to meteorological and industrial applications, environmental bioremediation, and biogeochemical cycles.

RNA-seq analysis of the response of the halophyte, Mesembryanthemum crystallinum (ice plant) to high salinity

Understanding the molecular mechanisms that convey salt tolerance in plants is a crucial issue for increasing crop yield. The ice plant (Mesembryanthemum crystallinum) is a halophyte that is capable of growing under high salt conditions. For example, the roots of ice plant seedlings continue to grow in 140 mM NaCl, a salt concentration that completely inhibits Arabidopsis thaliana root growth. Identifying the molecular mechanisms responsible for this high level of salt tolerance in a halophyte has the potential of revealing tolerance mechanisms that have been evolutionarily successful. In the present study, deep sequencing (RNAseq) was used to examine gene expression in ice plant roots treated with various concentrations of NaCl. Sequencing resulted in the identification of 53,516 contigs, 10,818 of which were orthologs of Arabidopsis genes. In addition to the expression analysis, a web-based ice plant database was constructed that allows broad public access to the data. The results obtained from an analysis of the RNAseq data were confirmed by RT-qPCR. Novel patterns of gene expression in response to high salinity within 24 hours were identified in the ice plant when the RNAseq data from the ice plant was compared to gene expression data obtained from Arabidopsis plants exposed to high salt. Although ABA responsive genes and a sodium transporter protein (HKT1), are up-regulated and down-regulated respectively in both Arabidopsis and the ice plant; peroxidase genes exhibit opposite responses. The results of this study provide an important first step towards analyzing environmental tolerance mechanisms in a non-model organism and provide a useful dataset for predicting novel gene functions.

Evolutionary histories of transposable elements in the genome of the largest living marsupial carnivore, the Tasmanian devil

The largest living carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii), is the sole survivor of a lineage originating about 12 Ma. We set out to investigate the spectrum of transposable elements found in the Tasmanian devil genome, the first high-coverage genome of an Australian marsupial. Marsupial genomes have been shown to have the highest amount of transposable elements among vertebrates. We analyzed the horizontally transmitted DNA transposons OC1 and hAT-1_MEu in the Tasmanian devil genome. OC1 is present in all carnivorous marsupials, while having a very limited distribution among the remaining Australian marsupial orders. In contrast, hAT-1_MEu is present in all Australian marsupial orders, and has so far only been identified in a few placental mammals. We screened 158 introns for phylogenetically informative retrotransposons in the order Dasyuromorphia, and found that the youngest SINE (Short INterspersed Element), WSINE1, is no longer active in the subfamily Dasyuridae. The lack of detectable WSINE1 activity in this group may be due to a retrotransposon inactivation event approximately 30 Ma. We found that the Tasmanian devil genome contains a relatively low number of continuous full-length LINE-1 (Long INterspersed Element 1, L1) retrotransposons compared with the opossum genome. Furthermore, all L1 elements in the Tasmanian devil appeared to be nonfunctional. Hidden Markov Model approaches suggested that other potential sources of functional reverse transcriptase are absent from the genome. We discuss the issues associated with assembling long, highly similar L1 copies from short read Illumina data and describe how assembly artifacts can potentially lead to erroneous conclusions.

Validation and assessment of variant calling pipelines for next-generation sequencing

Background

The processing and analysis of the large scale data generated by next-generation sequencing (NGS) experiments is challenging and is a burgeoning area of new methods development. Several new bioinformatics tools have been developed for calling sequence variants from NGS data. Here, we validate the variant calling of these tools and compare their relative accuracy to determine which data processing pipeline is optimal.

Results

We developed a unified pipeline for processing NGS data that encompasses four modules: mapping, filtering, realignment and recalibration, and variant calling. We processed 130 subjects from an ongoing whole exome sequencing study through this pipeline. To evaluate the accuracy of each module, we conducted a series of comparisons between the single nucleotide variant (SNV) calls from the NGS data and either gold-standard Sanger sequencing on a total of 700 variants or array genotyping data on a total of 9,935 single-nucleotide polymorphisms. A head to head comparison showed that Genome Analysis Toolkit (GATK) provided more accurate calls than SAMtools (positive predictive value of 92.55% vs. 80.35%, respectively). Realignment of mapped reads and recalibration of base quality scores before SNV calling proved to be crucial to accurate variant calling. GATK HaplotypeCaller algorithm for variant calling outperformed the UnifiedGenotype algorithm. We also showed a relationship between mapping quality, read depth and allele balance, and SNV call accuracy. However, if best practices are used in data processing, then additional filtering based on these metrics provides little gains and accuracies of >99% are achievable.

Conclusions

Our findings will help to determine the best approach for processing NGS data to confidently call variants for downstream analyses. To enable others to implement and replicate our results, all of our codes are freely available at http://metamoodics.org/wes.

(Post-)genomics approaches in fungal research

To date, hundreds of fungal genomes have been sequenced and many more are in progress. This wealth of genomic information has provided new directions to study fungal biodiversity. However, to further dissect and understand the complicated biological mechanisms involved in fungal life styles, functional studies beyond genomes are required. Thanks to the developments of current -omics techniques, it is possible to produce large amounts of fungal functional data in a high-throughput fashion (e.g. transcriptome, proteome, etc.). The increasing ease of creating -omics data has also created a major challenge for downstream data handling and analysis. Numerous databases, tools and software have been created to meet this challenge. Facing such a richness of techniques and information, hereby we provide a brief roadmap on current wet-lab and bioinformatics approaches to study functional genomics in fungi.

Human pharmacogenomic variation of antihypertensive drugs: from population genetics to personalized medicine

Aim: To investigate the human pharmacogenetic variation related to antihypertensive drugs, providing a survey of functional interpopulation differences in hypertension pharmacogenes. Materials & methods: The study was divided into two stages. In the first stage, we analyzed 1249 variants located in 57 hypertension pharmacogenes. This first-stage analysis confirmed that geographic origin strongly affects hypertension pharmacogenomic variation and that 31 pharmacogenes are geographically differentiated. In the second stage, we focused our attention on the ethnic-differentiated pharmacogenes, investigating 55,521 genetic variants. In silico analyses were performed to predict the effect of genetic variation. Results: Our analyses indicated functional interpopulation differences, suggesting insight into the mechanisms of antihypertensive drug response. Moreover, our data suggested that rare variants mainly determine the functionality of genes related to antihypertensive drugs. Conclusion: Our study provided important knowledge about the genetics of the antihypertensive drug response, suggesting that next-generation sequencing technologies may develop reliable pharmacogenetic tests for antihypertensive drugs.

Original submitted 19 September 2013; Revision submitted 14 November 2013

Bioprospecting in the genomic age

The genomic revolution promises great advances in the search for useful biocatalysts. Function-based metagenomic approaches have identified several enzymes with properties that make them useful candidates for a variety of bioprocesses. As DNA sequencing costs continue to decline, the volume of genomic data, along with their corresponding predicted protein sequences, will continue to increase dramatically, necessitating new approaches to leverage this information for gene-based bioprospecting efforts. Additionally, as new functions are discovered and correlated with this sequence information, the knowledge of the often complex relationship between a protein's sequence and function will improve. This in turn will lead to better gene-based bioprospecting approaches and facilitate the tailoring of desired properties through protein engineering projects. In this chapter, we discuss a number of recent advances in bioprospecting within the context of the genomic age.

Deep sequencing of mixed total DNA without barcodes allows efficient assembly of highly plastic ascidian mitochondrial genomes

Ascidians or sea squirts form a diverse group within chordates, which includes a few thousand members of marine sessile filter-feeding animals. Their mitochondrial genomes are characterized by particularly high evolutionary rates and rampant gene rearrangements. This extreme variability complicates standard polymerase chain reaction (PCR) based techniques for molecular characterization studies, and consequently only a few complete Ascidian mitochondrial genome sequences are available. Using the standard PCR and Sanger sequencing approach, we produced the mitochondrial genome of Ascidiella aspersa only after a great effort. In contrast, we produced five additional mitogenomes (Botrylloides aff. leachii, Halocynthia spinosa, Polycarpa mytiligera, Pyura gangelion, and Rhodosoma turcicum) with a novel strategy, consisting in sequencing the pooled total DNA samples of these five species using one Illumina HiSeq 2000 flow cell lane. Each mitogenome was efficiently assembled in a single contig using de novo transcriptome assembly, as de novo genome assembly generally performed poorly for this task. Each of the new six mitogenomes presents a different and novel gene order, showing that no syntenic block has been conserved at the ordinal level (in Stolidobranchia and in Phlebobranchia). Phylogenetic analyses support the paraphyly of both Ascidiacea and Phlebobranchia, with Thaliacea nested inside Phlebobranchia, although the deepest nodes of the Phlebobranchia-Thaliacea clade are not well resolved. The strategy described here thus provides a cost-effective approach to obtain complete mitogenomes characterized by a highly plastic gene order and a fast nucleotide/amino acid substitution rate.

AutoAssemblyD: a graphical user interface system for several genome assemblers

Next-generation sequencing technologies have increased the amount of biological data generated. Thus, bioinformatics has become important because new methods and algorithms are necessary to manipulate and process such data. However, certain challenges have emerged, such as genome assembly using short reads and high-throughput platforms. In this context, several algorithms have been developed, such as Velvet, Abyss, Euler-SR, Mira, Edna, Maq, SHRiMP, Newbler, ALLPATHS, Bowtie and BWA. However, most such assemblers do not have a graphical interface, which makes their use difficult for users without computing experience given the complexity of the assembler syntax. Thus, to make the operation of such assemblers accessible to users without a computing background, we developed AutoAssemblyD, which is a graphical tool for genome assembly submission and remote management by multiple assemblers through XML templates.

AVAILABILITY

AssemblyD is freely available at https://sourceforge.net/projects/autoassemblyd. It requires Sun jdk 6 or higher.

Assemblathon 2: evaluating de novo methods of genome assembly in three vertebrate species

BACKGROUND

The process of generating raw genome sequence data continues to become cheaper, faster, and more accurate. However, assembly of such data into high-quality, finished genome sequences remains challenging. Many genome assembly tools are available, but they differ greatly in terms of their performance (speed, scalability, hardware requirements, acceptance of newer read technologies) and in their final output (composition of assembled sequence). More importantly, it remains largely unclear how to best assess the quality of assembled genome sequences. The Assemblathon competitions are intended to assess current state-of-the-art methods in genome assembly.

RESULTS

In Assemblathon 2, we provided a variety of sequence data to be assembled for three vertebrate species (a bird, a fish, and snake). This resulted in a total of 43 submitted assemblies from 21 participating teams. We evaluated these assemblies using a combination of optical map data, Fosmid sequences, and several statistical methods. From over 100 different metrics, we chose ten key measures by which to assess the overall quality of the assemblies.

CONCLUSIONS

Many current genome assemblers produced useful assemblies, containing a significant representation of their genes and overall genome structure. However, the high degree of variability between the entries suggests that there is still much room for improvement in the field of genome assembly and that approaches which work well in assembling the genome of one species may not necessarily work well for another.

How are comparative genomics and the study of microRNAs changing our views on arthropod endocrinology and adaptations to the environment?

As the last few decades of work has shown, precise regulation of biosynthesis and release of arthropod hormones is essential to cope with environmental stresses and challenges. In crustaceans and insects, the sesquiterpenoids methyl farnesoate (MF), farnesoic acid (FA) and juvenile hormone (JH) regulate many developmental, physiological, and reproductive processes. In this review, we discuss how comparative genomics has and will impact our views on arthropod endocrinology. We will also highlight the current knowledge of regulation of genes involved in arthropod hormone biosynthesis by microRNAs, and describe the potential insights into arthropod endocrinology, evolution, and adaptation that are likely to come from the study of microRNAs.

De novo assembly of the complete organelle genome sequences of azuki bean (Vigna angularis) using next-generation sequencers

Since chloroplasts and mitochondria are maternally inherited and have unique features in evolution, DNA sequences of those organelle genomes have been broadly used in phylogenetic studies. Thanks to recent progress in next-generation sequencer (NGS) technology, whole-genome sequencing can be easily performed. Here, using NGS data generated by Roche GS Titanium and Illumina Hiseq 2000, we performed a hybrid assembly of organelle genome sequences of Vigna angularis (azuki bean). Both the mitochondrial genome (mtDNA) and the chloroplast genome (cpDNA) of V. angularis have very similar size and gene content to those of V. radiata (mungbean). However, in structure, mtDNA sequences have undergone many recombination events after divergence from the common ancestor of V. angularis and V. radiata, whereas cpDNAs are almost identical between the two. The stability of cpDNAs and the variability of mtDNAs was further confirmed by comparative analysis of Vigna organelles with model plants Lotus japonicus and Arabidopsis thaliana.

Genotyping and phenotyping of platelet function disorders

The majority of patients with platelet function disorders (PFDs) have normal platelet counts and mild day-to-day bleeding symptoms, but are at risk of major hemorrhage at times of trauma, surgery, or childbirth. This group is challenging to investigate, because the assays are often time-intensive and labour-intensive, and interpretation is difficult, especially in patients with mild disorders. In addition, interuser variability in performance of the assays, including the currently accepted gold standard, light transmission aggregometry, makes the results difficult to compare between laboratories. Furthermore, a similar pattern of mucocutaneous bleeding is seen in disorders in other components of the hemostatic pathway, including type 1 von Willebrand disease (VWD). We have undertaken an extensive investigation of patients with clinically diagnosed excessive bleeding, using a genotyping and platelet phenotyping approach based on lumi-aggregometry, and other specialist tests of platelet function, in combination with Sanger and next-generation sequencing (NGS). We found a functional defect in ~ 60% of patients, the majority being associated with feedback pathways of platelet activation. Function-disrupting mutations were identified in known and novel genes, and coinheritance with other genetic disorders of hemostasis, including type 1 VWD, was shown. A significant number of mutations are heterozygous and unlikely to cause extensive bleeding in isolation, consistent with incomplete penetrance of inheritance of bleeding disorders and a multifactorial etiology for excessive bleeding in many patients. Mucocutaneous bleeding is a complex trait, and this has important implications for NGS in the assessment of a PFD.

Proteomics of nonmodel plant species

Until recently, large scale proteomic investigations in the plant field have only been possible for a few model species for which the whole genome sequence had been fully determined. In contrast, for many other species with a strong economic interest as sources of human food and animal feed, as well as industrial and pharmacological molecules, little was known about their genome sequence and identifying the proteome in these species was still considered challenging. However, progress has been made as a result of several recent advances in proteomics tools, e.g. in MS technology and data search programs, and the increasing availability of genomic and cDNA sequences from various species. Moreover, next-generation sequencing technologies now make it possible to rapidly determine, at a reasonable cost, the genome or RNA sequence of species not currently considered as models, thus considerably expanding the plant sequence databases. This review will show how these advances make it possible to identify a large set of proteins, even for species for which few sequences are currently available.