Genome- and transcriptome-assisted development of nuclear insertion/deletion markers for Calanus species (Copepoda: Calanoida) identification

Copepods of the genus Calanus are key zooplankton species in temperate to arctic marine ecosystems. Despite their ecological importance, species identification remains challenging. Furthermore, the recent report of hybrids among Calanus species highlights the need for diagnostic nuclear markers to efficiently identify parental species and hybrids. Using next-generation sequencing analysis of both the genome and transcriptome from two sibling species, Calanus finmarchicus and Calanus glacialis, we developed a panel of 12 nuclear insertion/deletion markers. All the markers showed species-specific amplicon length. Furthermore, most of the markers were successfully amplified in other Calanus species, allowing the molecular identification of Calanus helgolandicus, Calanus hyperboreus and Calanus marshallae.

Parsing parallel evolution:Ecological divergence and differential gene expression in the adaptive radiations of thick-lipped Midas cichlid fishes from Nicaragua

The study of parallel evolution facilitates the discovery of common rules of diversification. Here, we examine the repeated evolution of thick lips in Midas cichlid fishes (the Amphilophus citrinellus species complex) - from two Great Lakes and two crater lakes in Nicaragua - to assess whether similar changes in ecology, phenotypic trophic traits and gene expression accompany parallel trait evolution. Using next-generation sequencing technology, we characterize transcriptome-wide differential gene expression in the lips of wild-caught sympatric thick- and thin-lipped cichlids from all four instances of repeated thick-lip evolution. Six genes (apolipoprotein D, myelin-associated glycoprotein precursor, four-and-a-half LIM domain protein 2, calpain-9, GTPase IMAP family member 8-like and one hypothetical protein) are significantly underexpressed in the thick-lipped morph across all four lakes. However, other aspects of lips' gene expression in sympatric morphs differ in a lake-specific pattern, including the magnitude of differentially expressed genes (97-510). Generally, fewer genes are differentially expressed among morphs in the younger crater lakes than in those from the older Great Lakes. Body shape, lower pharyngeal jaw size and shape, and stable isotopes (dC and dN) differ between all sympatric morphs, with the greatest differentiation in the Great Lake Nicaragua. Some ecological traits evolve in parallel (those related to foraging ecology; e.g. lip size, body and head shape) but others, somewhat surprisingly, do not (those related to diet and food processing; e.g. jaw size and shape, stable isotopes). Taken together, this case of parallelism among thick- and thin-lipped cichlids shows a mosaic pattern of parallel and nonparallel evolution. © 2012 Blackwell Publishing Ltd.

Exploring the Transcriptome of Atlantic Salmon (Salmo salar) Skin, a Major Defense Organ

The skin of fish is the first line of defense against pathogens and parasites. The skin transcriptome of the Atlantic salmon is poorly characterized, and currently only 2,089 expressed sequence tags (ESTs) out of a total of half a million sequences are generated from skin-derived cDNA libraries. The primary aim of this study was to enhance the transcriptomic knowledge of salmon skin by using next-generation sequencing (NGS) technology, namely the Roche-454 platform. An equimolar mixture of high-quality RNA from skin and epidermal samples of salmon reared in either freshwater or seawater was used for 454-sequencing. This technique yielded over 600,000 reads, which were assembled into 34,696 isotigs using Newbler. Of these isotigs, 12 % had not been sequenced in Atlantic salmon, hence representing previously unreported salmon mRNAs that can potentially be skin-specific. Many full-length genes have been acquired, representing numerous biological processes. Mucin proteins are the main structural component of mucus and we examined in greater detail the sequences we obtained for these genes. Several isotigs exhibited homology to mammalian mucins (MUC2, MUC5AC and MUC5B). Mucin mRNAs are generally > 10 kbp and contain large repetitive units, which pose a challenge towards full-length sequence discovery. To date, we have not unearthed any full-length salmon mucin genes with this dataset, but have both N- and C-terminal regions of a mucin type 5. This highlights the fact that, while NGS is indeed a formidable tool for sequence data mining of non-model species, it must be complemented with additional experimental and bioinformatic work to characterize some mRNA sequences with complex features.

Alternate wetting and drying irrigation for rice in Bangladesh: Is it sustainable and has plant breeding something to offer?

The crop management practice of alternate wetting and drying (AWD) is being promoted by IRRI and the national research and extension program in Bangladesh and other parts of the world as a water-saving irrigation practice that reduces the environmental impact of dry season rice production through decreased water usage, and potentially increases yield. Evidence is growing that AWD will dramatically reduce the concentration of arsenic in harvested rice grains conferring a third major advantage over permanently flooded dry season rice production. AWD may also increase the concentration of essential dietary micronutrients in the grain. However, three crucial aspects of AWD irrigation require further investigation. First, why is yield generally altered in AWD? Second, is AWD sustainable economically (viability of farmers' livelihoods) and environmentally (aquifer water table heights) over long-term use? Third, are current cultivars optimized for this irrigation system? This paper describes a multidisciplinary research project that could be conceived which would answer these questions by combining advanced soil biogeochemistry with crop physiology, genomics, and systems biology. The description attempts to show how the breakthroughs in next generation sequencing could be exploited to better utilize local collections of germplasm and identify the molecular mechanisms underlying biological adaptation to the environment within the context of soil chemistry and plant physiology.

Arbuscular mycorrhizal fungal communities are phylogenetically clustered at small scales

Next-generation sequencing technologies with markers covering the full Glomeromycota phylum were used to uncover phylogenetic community structure of arbuscular mycorrhizal fungi (AMF) associated with Festuca brevipila. The study system was a semi-arid grassland with high plant diversity and a steep environmental gradient in pH, C, N, P and soil water content. The AMF community in roots and rhizosphere soil were analyzed separately and consisted of 74 distinct operational taxonomic units (OTUs) in total. Community-level variance partitioning showed that the role of environmental factors in determining AM species composition was marginal when controlling for spatial autocorrelation at multiple scales. Instead, phylogenetic distance and spatial distance were major correlates of AMF communities: OTUs that were more closely related (and which therefore may have similar traits) were more likely to co-occur. This pattern was insensitive to phylogenetic sampling breadth. Given the minor effects of the environment, we propose that at small scales closely related AMF positively associate through biotic factors such as plant-AMF filtering and interactions within the soil biota.

Comparison of Microbiomes from Different Niches of Upper and Lower Airways in Children and Adolescents with Cystic Fibrosis

Changes in the airway microbiome may be important in the pathophysiology of chronic lung disease in patients with cystic fibrosis. However, little is known about the microbiome in early cystic fibrosis lung disease and the relationship between the microbiomes from different niches in the upper and lower airways. Therefore, in this cross-sectional study, we examined the relationship between the microbiome in the upper (nose and throat) and lower (sputum) airways from children with cystic fibrosis using next generation sequencing. Our results demonstrate a significant difference in both α and β-diversity between the nose and the two other sampling sites. The nasal microbiome was characterized by a polymicrobial community while the throat and sputum communities were less diverse and dominated by a few operational taxonomic units. Moreover, sputum and throat microbiomes were closely related especially in patients with clinically stable lung disease. There was a high inter-individual variability in sputum samples primarily due to a decrease in evenness linked to increased abundance of potential respiratory pathogens such as Pseudomonas aeruginosa. Patients with chronic Pseudomonas aeruginosa infection exhibited a less diverse sputum microbiome. A high concordance was found between pediatric and adult sputum microbiomes except that Burkholderia was only observed in the adult cohort. These results indicate that an adult-like lower airways microbiome is established early in life and that throat swabs may be a good surrogate in clinically stable children with cystic fibrosis without chronic Pseudomonas aeruginosa infection in whom sputum sampling is often not feasible.

Computational cancer biology: education is a natural key to many locks

Background: Oncology is a field that profits tremendously from the genomic data generated by high-throughput technologies, including next-generation sequencing. However, in order to exploit, integrate, visualize and interpret such high-dimensional data efficiently, non-trivial computational and statistical analysis methods are required that need to be developed in a problem-directed manner.

Discussion: For this reason, computational cancer biology aims to fill this gap. Unfortunately, computational cancer biology is not yet fully recognized as a coequal field in oncology, leading to a delay in its maturation and, as an immediate consequence, an under-exploration of high-throughput data for translational research.

Summary: Here we argue that this imbalance, favoring 'wet lab-based activities', will be naturally rectified over time, if the next generation of scientists receives an academic education that provides a fair and competent introduction to computational biology and its manifold capabilities. Furthermore, we discuss a number of local educational provisions that can be implemented on university level to help in facilitating the process of harmonization.

The fungal microbiota of de-novo paediatric inflammatory bowel disease

Inflammatory bowel disease (IBD) is characterised by an inappropriate chronic immune response against resident gut microbes. This may be on account of distinct changes in the gut microbiota termed as dysbiosis. The role of fungi in this altered luminal environment has been scarcely reported. We studied the fungal microbiome in de-novo paediatric IBD patients utilising next generation sequencing and compared with adult disease and normal controls. We report a distinct difference in fungal species with Ascomycota predominating in control subjects compared to Basidiomycota dominance in children with IBD, which could be as a result of altered tolerance in these patients. 

Molecular diagnostics of myeloproliferative neoplasms

Since the discovery of the JAK2 V617F mutation in the majority of the myeloproliferative neoplasms (MPN) of polycythemia vera, essential thrombocythemia and primary myelofibrosis ten years ago, further MPN-specific mutational events, notably in JAK2 exon 12, MPL exon 10 and CALR exon 9 have been identified. These discoveries have been rapidly incorporated into evolving molecular diagnostic algorithms. While many of these mutations appear to have prognostic implications, establishing MPN diagnosis is of immediate clinical importance with selection, implementation and the continual evaluation of the appropriate laboratory methodology to achieve this diagnosis similarly vital. The advantages and limitations of these approaches in identifying and quantitating the common MPN-associated mutations is considered herein with particular regard to their clinical utility. The evolution of molecular diagnostic applications and platforms has occurred in parallel with the discovery of MPN-associated mutations and it therefore appears likely that emerging technologies such as next-generation sequencing and digital PCR will in the future, play an increasing role in the molecular diagnosis of MPN. 

Novel tools for conservation genomics: comparing two high-throughput approaches for SNP discovery in the transcriptome of the European hake

The growing accessibility to genomic resources using next-generation sequencing (NGS) technologies has revolutionized the application of molecular genetic tools to ecology and evolutionary studies in non-model organisms. Here we present the case study of the European hake (Merluccius merluccius), one of the most important demersal resources of European fisheries. Two sequencing platforms, the Roche 454 FLX (454) and the Illumina Genome Analyzer (GAII), were used for Single Nucleotide Polymorphisms (SNPs) discovery in the hake muscle transcriptome. De novo transcriptome assembly into unique contigs, annotation, and in silico SNP detection were carried out in parallel for 454 and GAII sequence data. High-throughput genotyping using the Illumina GoldenGate assay was performed for validating 1,536 putative SNPs. Validation results were analysed to compare the performances of 454 and GAII methods and to evaluate the role of several variables (e.g. sequencing depth, intron-exon structure, sequence quality and annotation). Despite well-known differences in sequence length and throughput, the two approaches showed similar assay conversion rates (approximately 43%) and percentages of polymorphic loci (67.5% and 63.3% for GAII and 454, respectively). Both NGS platforms therefore demonstrated to be suitable for large scale identification of SNPs in transcribed regions of non-model species, although the lack of a reference genome profoundly affects the genotyping success rate. The overall efficiency, however, can be improved using strict quality and filtering criteria for SNP selection (sequence quality, intron-exon structure, target region score).

Application of SNPs for population genetics of nonmodel organisms:new opportunities and challenges

Recent improvements in the speed, cost and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs). SNPs are increasingly being used as an addition to the molecular ecology toolkit in nonmodel organisms, but their efficient use remains challenging. Here, we discuss common issues when employing SNP markers, including the high numbers of markers typically employed, the effects of ascertainment bias and the inclusion of nonneutral loci in a marker panel. We provide a critique of considerations specifically associated with the application and population genetic analysis of SNPs in nonmodel taxa, focusing specifically on some of the most commonly applied methods.

Isolation and characterization of thirty microsatellite loci for Atlantic mackerel (Scomber scombrus L.)

Atlantic mackerel (Scomber scombrus L.) is a widely distributed commercially important pelagic species. Little is known about the stock structure of this species, but it is thought to be undergoing a range extension due to environmental changes. Knowledge of the stock structure under these changing conditions is fundamental for effective management. In this paper, 30 highly polymorphic microsatellite loci developed with next generation sequencing are described. The number of alleles per locus ranged from 4 to 39 in two geographically distant populations, observed and expected heterozygosities ranged between 0. 370-0. 978 and 0. 426-0. 962, respectively. These loci are an important resource that will allow assessment of the current population genetic structure of this species, and enable monitoring of climate related changes in the species range and distribution. 

Twenty-two novel microsatellite loci for lumpfish (Cyclopterus lumpus)

Lumpfish (Cyclopterus lumpus L. 1758) are widely distributed on both sides of the North Atlantic. They are a commercially important species, but stock size estimates have declined since the mid-1980s in Canada, Norway and Iceland. Little is known about the biology of this species, in particular the breeding migrations and population structure which are fundamental for effective management. This paper describes the development and characterization of twenty-two polymorphic microsatellite loci using next generation sequencing. The number of alleles per locus ranged from 3 to 27 in two geographically distant North Atlantic populations, with observed and expected heterozygosities ranging between 0. 0625-0. 979 and 0. 0618-0. 946, respectively. These loci are an important resource that will allow assessment of the population genetic structure of this species, and contribute to its appropriate management.

Multiplexing with three-primer PCR for rapid and economical microsatellite validation.

The next generation sequencing revolution has enabled rapid discovery of genetic markers, however, development of fully functioning new markers still requires a long and costly process of marker validation. This study reports a rapid and economical approach for the validation and deployment of polymorphic microsatellite markers obtained from a 454 pyrosequencing library of Atlantic cod, Gadus morhua, Linnaeus 1758. Primers were designed from raw reads to amplify specific amplicon size ranges, allowing effective PCR multiplexing. Multiplexing was combined with a three-primer PCR approach using four universal tails to label amplicons with separate fluorochromes. A total of 192 primer pairs were tested, resulting in 73 polymorphic markers. Of these, 55 loci were combined in six multiplex panels each containing between six and eleven markers. Variability of the loci was assessed on G. morhua from the Celtic Sea (n 46) and the Scotian Shelf (n 46), two locations that have shown genetic differentiation in previous studies. Multilocus FST between the two samples was estimated at 0.067 (P 0.001). After three loci potentially under selection were excluded, the global FST was estimated at 0.043 (P 0.001). Our technique combines three- primer and multiplex PCR techniques, allowing simultaneous screening and validation of relatively large numbers of microsatellite loci.

Comprehensive molecular pathology analysis of small bowel adenocarcinoma reveals novel targets with potential for clinical utility

Small bowel accounts for only 0.5% of cancer cases in the US but incidence rates have been rising at 2.4% per year over the past decade. One-third of these are adenocarcinomas but little is known about their molecular pathology and no molecular markers are available for clinical use. Using a retrospective 28 patient matched normal-tumor cohort, next-generation sequencing, gene expression arrays and CpG methylation arrays were used for molecular profiling. Next-generation sequencing identified novel mutations in IDH1, CDH1, KIT, FGFR2, FLT3, NPM1, PTEN, MET, AKT1, RET, NOTCH1 and ERBB4. Array data revealed 17% of CpGs and 5% of RNA transcripts assayed to be differentially methylated and expressed respectively (p < 0.01). Merging gene expression and DNA methylation data revealed CHN2 as consistently hypermethylated and downregulated in this disease (Spearman -0.71, p < 0.001). Mutations in TP53 which were found in more than half of the cohort (15/28) and Kazald1 hypomethylation were both were indicative of poor survival (p = 0.03, HR = 3.2 and p = 0.01, HR = 4.9 respectively). By integrating high-throughput mutational, gene expression and DNA methylation data, this study reveals for the first time the distinct molecular profile of small bowel adenocarcinoma and highlights potential clinically exploitable markers.