Trimmomatic: a flexible trimmer for Illumina sequence data

Growth dynamics of gut microbiota in health and disease inferred from single metagenomic samples

Metagenomic sequencing increased our understanding of the role of the microbiome in health and disease, yet it only provides a snapshot of a highly dynamic ecosystem. Here, we show that the pattern of metagenomic sequencing read coverage for different microbial genomes contains a single trough and a single peak, the latter coinciding with the bacterial origin of replication. Furthermore, the ratio of sequencing coverage between the peak and trough provides a quantitative measure of a species' growth rate. We demonstrate this in vitro and in vivo, under different growth conditions, and in complex bacterial communities. For several bacterial species, peak-to-trough coverage ratios, but not relative abundances, correlated with the manifestation of inflammatory bowel disease and type II diabetes.

Comparative Transcriptomes Analysis of Red- and White-Fleshed Apples in an F1 Population of Malus sieversii f. niedzwetzkyana Crossed with M. domestica 'Fuji'

Transcriptome profiles of the red- and white-fleshed apples in an F₁ segregating population of Malus sieversii f.Niedzwetzkyana and M.domestica ‘Fuji’ were generated using the next-generation high-throughput RNA sequencing (RNA-Seq) technology and compared. A total of 114 differentially expressed genes (DEGs) were obtained, of which 88 were up-regulated and 26 were down-regulated in red-fleshed apples. The 88 up-regulated genes were enriched with those related to flavonoid biosynthetic process and stress responses. Further analysis identified 22 genes associated with flavonoid biosynthetic process and 68 genes that may be related to stress responses. Furthermore, the expression of 20 up-regulated candidate genes (10 related to flavonoid biosynthesis, two encoding MYB transcription factors and eight related to stress responses) and 10 down-regulated genes were validated by quantitative real-time PCR. After exploring the possible regulatory network, we speculated that flavonoid metabolism might be involved in stress responses in red-fleshed apple. Our findings provide a theoretical basis for further enriching gene resources associated with flavonoid synthesis and stress responses of fruit trees and for breeding elite apples with high flavonoid content and/or increased stress tolerances.

Virulence characteristics of hcp (+) Campylobacter jejuni and Campylobacter coli isolates from retail chicken

Background

Recently the Type VI secretion system (T6SS), which can play a significant role in bacterial survival and pathogenesis, was reported in Campylobacter spp., having the hcp gene as a key component.

Methods

Campylobacteriosis is associated with the consumption of infected chicken meat. Our study aimed to explore the presence of T6SS in C. jejuni (n = 59) and C. coli (n = 57) isolates, from retail raw chicken and to investigate their pathogenic potential. The hcp gene was used as an indicator for the T6SS presence.

Results

Using multiplex PCR we have identified a significantly higher prevalence of hcp in C. coli isolates (56.1%) than in C. jejuni (28.8%) and AFLP analysis of the isolates showed a high degree of genetic similarity between the isolates carrying the hcp gene. Genome sequencing data showed that 84.3% of the C. coli and 93.7% of the C. jejuni isolates had all 13 T6SS open reading frames. Moreover, the virulence characteristics of hcp + isolates, including motility and the ability to invade human intestinal epithelial cells in vitro, were significantly greater than in the control strain C. jejuni 12502; a human isolate which is hcp positive.

Conclusion

Overall, it was discovered that hcp⁺C. coli and C. jejuni isolated from retail chicken isolates posses genetic and phenotypic properties associated with enhanced virulence. However, since human infections with C. coli are significantly less frequent than those of C. jejuni, the relationship between virulence factors and pathogenesis requires further study.

Vitamin D Modulates Expression of the Airway Smooth Muscle Transcriptome in Fatal Asthma

Globally, asthma is a chronic inflammatory respiratory disease affecting over 300 million people. Some asthma patients remain poorly controlled by conventional therapies and experience more life-threatening exacerbations. Vitamin D, as an adjunct therapy, may improve disease control in severe asthma patients since vitamin D enhances glucocorticoid responsiveness and mitigates airway smooth muscle (ASM) hyperplasia. We sought to characterize differences in transcriptome responsiveness to vitamin D between fatal asthma- and non-asthma-derived ASM by using RNA-Seq to measure ASM transcript expression in five donors with fatal asthma and ten non-asthma-derived donors at baseline and with vitamin D treatment. Based on a Benjamini-Hochberg corrected p-value <0.05, 838 genes were differentially expressed in fatal asthma vs. non-asthma-derived ASM at baseline, and vitamin D treatment compared to baseline conditions induced differential expression of 711 and 867 genes in fatal asthma- and non-asthma-derived ASM, respectively. Functional gene categories that were highly represented in all groups included extracellular matrix, and responses to steroid hormone stimuli and wounding. Genes differentially expressed by vitamin D also included cytokine and chemokine activity categories. Follow-up qPCR and individual analyte ELISA experiments were conducted for four cytokines (i.e. CCL2, CCL13, CXCL12, IL8) to measure TNFα-induced changes by asthma status and vitamin D treatment. Vitamin D inhibited TNFα-induced IL8 protein secretion levels to a comparable degree in fatal asthma- and non-asthma-derived ASM even though IL8 had significantly higher baseline levels in fatal asthma-derived ASM. Our findings identify vitamin D-specific gene targets and provide transcriptomic data to explore differences in the ASM of fatal asthma- and non-asthma-derived donors.

The progesterone antagonist mifepristone/RU486 blocks the negative effect on life span caused by mating in female Drosophila

Mating causes decreased life span in female Drosophila. Here we report that mifepristone blocked this effect, yielding life span increases up to +68%. Drug was fed to females after mating, in the absence of males, demonstrating function in females. Mifepristone did not increase life span of virgin females or males. Mifepristone reduced progeny production but did not reduce food intake. High-throughput RNA sequencing was used to identify genes up-regulated or down-regulated upon mating, and where the change was reduced by mifepristone. Five candidate positive regulators of life span were identified, including dosage compensation regulator Unr and three X-linked genes: multi sex combs (PcG gene), Dopamine 2-like receptor and CG14215. The 37 candidate negative genes included neuropeptide CNMamide and several involved in protein mobilization and immune response. The results inform the interpretation of experiments involving mifepristone, and implicate steroid hormone signaling in regulating the trade-off between reproduction and life span.

Identifying genomic changes associated with insecticide resistance in the dengue mosquito Aedes aegypti by deep targeted sequencing

The capacity of mosquitoes to resist insecticides threatens the control of diseases such as dengue and malaria. Until alternative control tools are implemented, characterizing resistance mechanisms is crucial for managing resistance in natural populations. Insecticide biodegradation by detoxification enzymes is a common resistance mechanism; however, the genomic changes underlying this mechanism have rarely been identified, precluding individual resistance genotyping. In particular, the role of copy number variations (CNVs) and polymorphisms of detoxification enzymes have never been investigated at the genome level, although they can represent robust markers of metabolic resistance. In this context, we combined target enrichment with high-throughput sequencing for conducting the first comprehensive screening of gene amplifications and polymorphisms associated with insecticide resistance in mosquitoes. More than 760 candidate genes were captured and deep sequenced in several populations of the dengue mosquito Ae. aegypti displaying distinct genetic backgrounds and contrasted resistance levels to the insecticide deltamethrin. CNV analysis identified 41 gene amplifications associated with resistance, most affecting cytochrome P450s overtranscribed in resistant populations. Polymorphism analysis detected more than 30,000 variants and strong selection footprints in specific genomic regions. Combining Bayesian and allele frequency filtering approaches identified 55 nonsynonymous variants strongly associated with resistance. Both CNVs and polymorphisms were conserved within regions but differed across continents, confirming that genomic changes underlying metabolic resistance to insecticides are not universal. By identifying novel DNA markers of insecticide resistance, this study opens the way for tracking down metabolic changes developed by mosquitoes to resist insecticides within and among populations.

Development of a Multilocus Sequence Typing Scheme for Molecular Typing of Mycoplasma pneumoniae

Mycoplasma pneumoniae is a major human respiratory pathogen causing both upper and lower respiratory disease in humans of all ages, and it can also result in other serious extrapulmonary sequelae. A multilocus sequence typing (MLST) scheme for M. pneumoniae was developed based on the sequences of eight housekeeping genes (ppa, pgm, gyrB, gmk, glyA, atpA, arcC, and adk) and applied to 55 M. pneumoniae clinical isolates and the two type strains M129 and FH. A total of 12 sequence types (STs) resulted for 57 M. pneumoniae isolates tested, with a discriminatory index of 0.21 STs per isolate. The MLST loci used in this scheme were shown to be stable in 10 strains following 10 sequential subculture passages. Phylogenetic analysis of concatenated sequences of the eight loci indicated two distinct genetic clusters that were directly linked to multilocus variable-number tandem repeat analysis (MLVA) type. Genetic MLST clustering was confirmed by genomic sequence analysis, indicating that the MLST scheme developed in this study is representative of the genome. Furthermore, this MLST scheme was shown to be more discriminatory than both MLVA and P1 typing for the M. pneumoniae isolates examined, providing a method for further and more detailed analysis of observed epidemic peaks of M. pneumoniae infection. This scheme is supported by a public Web-based database (http://pubmlst.org/mpneumoniae).

TAL effectors and activation of predicted host targets distinguish Asian from African strains of the rice pathogen Xanthomonas oryzae pv. oryzicola while strict conservation suggests universal importance of five TAL effectors

Xanthomonas oryzae pv. oryzicola (Xoc) causes the increasingly important disease bacterial leaf streak of rice (BLS) in part by type III delivery of repeat-rich transcription activator-like (TAL) effectors to upregulate host susceptibility genes. By pathogen whole genome, single molecule, real-time sequencing and host RNA sequencing, we compared TAL effector content and rice transcriptional responses across 10 geographically diverse Xoc strains. TAL effector content is surprisingly conserved overall, yet distinguishes Asian from African isolates. Five TAL effectors are conserved across all strains. In a prior laboratory assay in rice cv. Nipponbare, only two contributed to virulence in strain BLS256 but the strict conservation indicates all five may be important, in different rice genotypes or in the field. Concatenated and aligned, TAL effector content across strains largely reflects relationships based on housekeeping genes, suggesting predominantly vertical transmission. Rice transcriptional responses did not reflect these relationships, and on average, only 28% of genes upregulated and 22% of genes downregulated by a strain are up- and down- regulated (respectively) by all strains. However, when only known TAL effector targets were considered, the relationships resembled those of the TAL effectors. Toward identifying new targets, we used the TAL effector-DNA recognition code to predict effector binding elements in promoters of genes upregulated by each strain, but found that for every strain, all upregulated genes had at least one. Filtering with a classifier we developed previously decreases the number of predicted binding elements across the genome, suggesting that it may reduce false positives among upregulated genes. Applying this filter and eliminating genes for which upregulation did not strictly correlate with presence of the corresponding TAL effector, we generated testable numbers of candidate targets for four of the five strictly conserved TAL effectors.

The reduced genomes of Parcubacteria (OD1) contain signatures of a symbiotic lifestyle

Candidate phylum OD1 bacteria (also referred to as Parcubacteria) have been identified in a broad range of anoxic environments through community survey analysis. Although none of these species have been isolated in the laboratory, several genome sequences have been reconstructed from metagenomic sequence data and single-cell sequencing. The organisms have small (generally <1 Mb) genomes with severely reduced metabolic capabilities. We have reconstructed 8 partial to near-complete OD1 genomes from oxic groundwater samples, and compared them against existing genomic data. The conserved core gene set comprises 202 genes, or ~28% of the genomic complement. “Housekeeping” genes and genes for biosynthesis of peptidoglycan and Type IV pilus production are conserved. Gene sets for biosynthesis of cofactors, amino acids, nucleotides, and fatty acids are absent entirely or greatly reduced. The only aspects of energy metabolism conserved are the non-oxidative branch of the pentose-phosphate shunt and central glycolysis. These organisms also lack some activities conserved in almost all other known bacterial genomes, including signal recognition particle, pseudouridine synthase A, and FAD synthase. Pan-genome analysis indicates a broad genotypic diversity and perhaps a highly fluid gene complement, indicating historical adaptation to a wide range of growth environments and a high degree of specialization. The genomes were examined for signatures suggesting either a free-living, streamlined lifestyle, or a symbiotic lifestyle. The lack of biosynthetic capabilities and DNA repair, along with the presence of potential attachment and adhesion proteins suggest that the Parcubacteria are ectosymbionts or parasites of other organisms. The wide diversity of genes that potentially mediate cell-cell contact suggests a broad range of partner/prey organisms across the phylum.

Genomic Analysis of the Emergence and Rapid Global Dissemination of the Clonal Group 258 Klebsiella pneumoniae Pandemic

Multidrug-resistant Klebsiella pneumoniae producing the KPC carbapenemase have rapidly spread throughout the world, causing severe healthcare-associated infections with limited antimicrobial treatment options. Dissemination of KPC-producing K. pneumoniae is largely attributed to expansion of a single dominant strain, ST258. In this study, we explore phylogenetic relationships and evolution within ST258 and its clonal group, CG258, using whole genome sequence analysis of 167 isolates from 20 countries collected over 17 years. Our results show a common ST258 ancestor emerged from its diverse parental clonal group around 1995 and likely acquired blaKPC prior to dissemination. Over the past two decades, ST258 has remained highly clonal despite diversity in accessory elements and divergence in the capsule polysaccharide synthesis locus. Apart from the large recombination event that gave rise to ST258, few mutations set it apart from its clonal group. However, one mutation occurs in a global transcription regulator. Characterization of outer membrane protein sequences revealed a profile in ST258 that includes a truncated OmpK35 and modified OmpK37. Our work illuminates potential genomic contributors to the pathogenic success of ST258, helps us better understand the global dissemination of this strain, and identifies genetic markers unique to ST258.

Transposon Mutagenesis Paired with Deep Sequencing of Caulobacter crescentus under Uranium Stress Reveals Genes Essential for Detoxification and Stress Tolerance

The ubiquitous aquatic bacterium Caulobacter crescentus is highly resistant to uranium (U) and facilitates U biomineralization and thus holds promise as an agent of U bioremediation. To gain an understanding of how C. crescentus tolerates U, we employed transposon (Tn) mutagenesis paired with deep sequencing (Tn-seq) in a global screen for genomic elements required for U resistance. Of the 3,879 annotated genes in the C. crescentus genome, 37 were found to be specifically associated with fitness under U stress, 15 of which were subsequently tested through mutational analysis. Systematic deletion analysis revealed that mutants lacking outer membrane transporters (rsaFa and rsaFb), a stress-responsive transcription factor (cztR), or a ppGpp synthetase/hydrolase (spoT) exhibited a significantly lower survival rate under U stress. RsaFa and RsaFb, which are homologues of TolC in Escherichia coli, have previously been shown to mediate S-layer export. Transcriptional analysis revealed upregulation of rsaFa and rsaFb by 4- and 10-fold, respectively, in the presence of U. We additionally show that rsaFa mutants accumulated higher levels of U than the wild type, with no significant increase in oxidative stress levels. Our results suggest a function for RsaFa and RsaFb in U efflux and/or maintenance of membrane integrity during U stress. In addition, we present data implicating CztR and SpoT in resistance to U stress. Together, our findings reveal novel gene targets that are key to understanding the molecular mechanisms of U resistance in C. crescentus.

IMPORTANCE

Caulobacter crescentus is an aerobic bacterium that is highly resistant to uranium (U) and has great potential to be used in U bioremediation, but its mechanisms of U resistance are poorly understood. We conducted a Tn-seq screen to identify genes specifically required for U resistance in C. crescentus. The genes that we identified have previously remained elusive using other omics approaches and thus provide significant insight into the mechanisms of U resistance by C. crescentus. In particular, we show that outer membrane transporters RsaFa and RsaFb, previously known as part of the S-layer export machinery, may confer U resistance by U efflux and/or by maintaining membrane integrity during U stress.

Root precursors of microRNAs in wild emmer and modern wheats show major differences in response to drought stress

MicroRNAs, small regulatory molecules with significant impacts on the transcriptional network of all living organisms, have been the focus of several studies conducted mostly on modern wheat cultivars. In this study, we investigated miRNA repertoires of modern durum wheat and its wild relatives, with differing degrees of drought tolerance, to identify miRNA candidates and their targets involved in drought stress response. Root transcriptomes of Triticum turgidum ssp. durum variety Kızıltan and two Triticum turgidum ssp. dicoccoides genotypes TR39477 and TTD-22 under control and drought conditions were assembled from individual RNA-Seq reads and used for in silico identification of miRNAs. A total of 66 miRNAs were identified from all species, across all conditions, of which 46 and 38 of the miRNAs identified from modern durum wheat and wild genotypes, respectively, had not been previously reported. Genotype- and/or stress-specific miRNAs provide insights into our understanding of the complex drought response. Particularly, miR1435, miR5024, and miR7714, identified only from drought-stress roots of drought-tolerant genotype TR39477, can be candidates for future studies to explore and exploit the drought response to develop tolerant varieties.

Genomic Epidemiology of a Protracted Hospital Outbreak Caused by a Toxin A-Negative Clostridium difficile Sublineage PCR Ribotype 017 Strain in London, England

Clostridium difficile remains the leading cause of nosocomial diarrhea worldwide, which is largely considered to be due to the production of two potent toxins: TcdA and TcdB. However, PCR ribotype (RT) 017, one of five clonal lineages of human virulent C. difficile, lacks TcdA expression but causes widespread disease. Whole-genome sequencing was applied to 35 isolates from hospitalized patients with C. difficile infection (CDI) and two environmental ward isolates in London, England. The phylogenetic analysis of single nucleotide polymorphisms (SNPs) revealed a clonal cluster of temporally variable isolates from a single hospital ward at University Hospital Lewisham (UHL) that were distinct from other London hospital isolates. De novo assembled genomes revealed a 49-kbp putative conjugative transposon exclusive to this hospital clonal cluster which would not be revealed by current typing methodologies. This study identified three sublineages of C. difficile RT017 that are circulating in London. Similar to the notorious RT027 lineage, which has caused global outbreaks of CDI since 2001, the lineage of toxin-defective RT017 strains appears to be continually evolving. By utilization of WGS technologies to identify SNPs and the evolution of clonal strains, the transmission of outbreaks caused by near-identical isolates can be retraced and identified.

Evolution of Mating Systems in Basidiomycetes and the Genetic Architecture Underlying Mating-Type Determination in the Yeast Leucosporidium scottii

In most fungi, sexual reproduction is bipolar; that is, two alternate sets of genes at a single mating-type (MAT) locus determine two mating types. However, in the Basidiomycota, a unique (tetrapolar) reproductive system emerged in which sexual identity is governed by two unlinked MAT loci, each of which controls independent mechanisms of self/nonself recognition. Tetrapolar-to-bipolar transitions have occurred on multiple occasions in the Basidiomycota, resulting, for example, from linkage of the two MAT loci into a single inheritable unit. Nevertheless, owing to the scarcity of molecular data regarding tetrapolar systems in the earliest-branching lineage of the Basidiomycota (subphylum Pucciniomycotina), it is presently unclear if the last common ancestor was tetrapolar or bipolar. Here, we address this question, by investigating the mating system of the Pucciniomycotina yeast Leucosporidium scottii. Using whole-genome sequencing and chromoblot analysis, we discovered that sexual reproduction is governed by two physically unlinked gene clusters: a multiallelic homeodomain (HD) locus and a pheromone/receptor (P/R) locus that is biallelic, thereby dismissing the existence of a third P/R allele as proposed earlier. Allele distribution of both MAT genes in natural populations showed that the two loci were in strong linkage disequilibrium, but independent assortment of MAT alleles was observed in the meiotic progeny of a test cross. The sexual cycle produces fertile progeny with similar proportions of the four mating types, but approximately 2/3 of the progeny was found to be nonhaploid. Our study adds to others in reinforcing tetrapolarity as the ancestral state of all basidiomycetes.

The genome of the sparganosis tapeworm Spirometra erinaceieuropaei isolated from the biopsy of a migrating brain lesion

BACKGROUND

Sparganosis is an infection with a larval Diphyllobothriidea tapeworm. From a rare cerebral case presented at a clinic in the UK, DNA was recovered from a biopsy sample and used to determine the causative species as Spirometra erinaceieuropaei through sequencing of the cox1 gene. From the same DNA, we have produced a draft genome, the first of its kind for this species, and used it to perform a comparative genomics analysis and to investigate known and potential tapeworm drug targets in this tapeworm.

RESULTS

The 1.26 Gb draft genome of S. erinaceieuropaei is currently the largest reported for any flatworm. Through investigation of β-tubulin genes, we predict that S. erinaceieuropaei larvae are insensitive to the tapeworm drug albendazole. We find that many putative tapeworm drug targets are also present in S. erinaceieuropaei, allowing possible cross application of new drugs. In comparison to other sequenced tapeworm species we observe expansion of protease classes, and of Kuntiz-type protease inhibitors. Expanded gene families in this tapeworm also include those that are involved in processes that add post-translational diversity to the protein landscape, intracellular transport, transcriptional regulation and detoxification.

CONCLUSIONS

The S. erinaceieuropaei genome begins to give us insight into an order of tapeworms previously uncharacterized at the genome-wide level. From a single clinical case we have begun to sketch a picture of the characteristics of these organisms. Finally, our work represents a significant technological achievement as we present a draft genome sequence of a rare tapeworm, and from a small amount of starting material.

Genome Sequence of Vibrio VPAP30, Isolated from an Episode of Massive Mortality of Reared Larvae of the Scallop Argopecten purpuratus

We report here the 5.167-Mbp draft genome sequence of Vibrio VPAP30, isolated from an Argopecten purpuratus larval culture. Vibrio VPAP30 is the etiological agent of a vibriosis outbreak causing a complete collapse of a larval culture of the scallop A. purpuratus, which occurred in a commercial hatchery in Chile.

Comparative Phylodynamics of Rabbit Hemorrhagic Disease Virus in Australia and New Zealand

The introduction of rabbit hemorrhagic disease virus (RHDV) into Australia and New Zealand during the 1990s as a means of controlling feral rabbits is an important case study in viral emergence. Both epidemics are exceptional in that the founder viruses share an origin and the timing of their release is known, providing a unique opportunity to compare the evolution of a single virus in distinct naive populations. We examined the evolution and spread of RHDV in Australia and New Zealand through a genome-wide evolutionary analysis, including data from 28 newly sequenced RHDV field isolates. Following the release of the Australian inoculum strain into New Zealand, no subsequent mixing of the populations occurred, with viruses from both countries forming distinct groups. Strikingly, the rate of evolution in the capsid gene was higher in the Australian viruses than in those from New Zealand, most likely due to the presence of transient deleterious mutations in the former. However, estimates of both substitution rates and population dynamics were strongly sample dependent, such that small changes in sample composition had an important impact on evolutionary parameters. Phylogeographic analysis revealed a clear spatial structure in the Australian RHDV strains, with a major division between those viruses from western and eastern states. Importantly, RHDV sequences from the state where the virus was first released, South Australia, had the greatest diversity and were diffuse throughout both geographic lineages, such that this region was likely a source population for the subsequent spread of the virus across the country.

IMPORTANCE

Most studies of viral emergence lack detailed knowledge about which strains were founders for the outbreak or when these events occurred. Hence, the human-mediated introduction of rabbit hemorrhagic disease virus (RHDV) into Australia and New Zealand from known starting stocks provides a unique opportunity to understand viral evolution and emergence. Within Australia, we revealed a major phylogenetic division between viruses sampled from the east and west of the country, with both regions likely seeded by viruses from South Australia. Despite their common origins, marked differences in evolutionary rates were observed between the Australian and New Zealand RHDV, which led to conflicting conclusions about population growth rates. An analysis of mutational patterns suggested that evolutionary rates have been elevated in the Australian viruses, at least in part due to the presence of low-fitness (deleterious) variants that have yet to be selectively purged.

The Genotoxin Colibactin Is a Determinant of Virulence in Escherichia coli K1 Experimental Neonatal Systemic Infection

Escherichia coli strains expressing the K1 capsule are a major cause of sepsis and meningitis in human neonates. The development of these diseases is dependent on the expression of a range of virulence factors, many of which remain uncharacterized. Here, we show that all but 1 of 34 E. coli K1 neonatal isolates carried clbA and clbP, genes contained within the pks pathogenicity island and required for the synthesis of colibactin, a polyketide-peptide genotoxin that causes genomic instability in eukaryotic cells by induction of double-strand breaks in DNA. Inactivation of clbA and clbP in E. coli A192PP, a virulent strain of serotype O18:K1 that colonizes the gastrointestinal tract and translocates to the blood compartment with very high frequency in experimental infection of the neonatal rat, significantly reduced the capacity of A192PP to colonize the gut, engender double-strand breaks in DNA, and cause invasive, lethal disease. Mutation of clbA, which encodes a pleiotropic enzyme also involved in siderophore synthesis, impacted virulence to a greater extent than mutation of clbP, encoding an enzyme specific to colibactin synthesis. Restoration of colibactin gene function by complementation reestablished the fully virulent phenotype. We conclude that colibactin contributes to the capacity of E. coli K1 to colonize the neonatal gastrointestinal tract and to cause invasive disease in the susceptible neonate.

Transcriptomic profiling of gene expression and RNA processing during Leishmania major differentiation

Protozoan parasites of the genus Leishmania are the etiological agents of leishmaniasis, a group of diseases with a worldwide incidence of 0.9–1.6 million cases per year. We used RNA-seq to conduct a high-resolution transcriptomic analysis of the global changes in gene expression and RNA processing events that occur as L. major transforms from non-infective procyclic promastigotes to infective metacyclic promastigotes. Careful statistical analysis across multiple biological replicates and the removal of batch effects provided a high quality framework for comprehensively analyzing differential gene expression and transcriptome remodeling in this pathogen as it acquires its infectivity. We also identified precise 5′ and 3′ UTR boundaries for a majority of Leishmania genes and detected widespread alternative trans-splicing and polyadenylation. An investigation of possible correlations between stage-specific preferential trans-splicing or polyadenylation sites and differentially expressed genes revealed a lack of systematic association, establishing that differences in expression levels cannot be attributed to stage-regulated alternative RNA processing. Our findings build on and improve existing expression datasets and provide a substantially more detailed view of L. major biology that will inform the field and potentially provide a stronger basis for drug discovery and vaccine development efforts.

IFNγ and IL-12 Restrict Th2 Responses during Helminth/Plasmodium Co-Infection and Promote IFNγ from Th2 Cells

Parasitic helminths establish chronic infections in mammalian hosts. Helminth/Plasmodium co-infections occur frequently in endemic areas. However, it is unclear whether Plasmodium infections compromise anti-helminth immunity, contributing to the chronicity of infection. Immunity to Plasmodium or helminths requires divergent CD4⁺ T cell-driven responses, dominated by IFNγ or IL-4, respectively. Recent literature has indicated that Th cells, including Th2 cells, have phenotypic plasticity with the ability to produce non-lineage associated cytokines. Whether such plasticity occurs during co-infection is unclear. In this study, we observed reduced anti-helminth Th2 cell responses and compromised anti-helminth immunity during Heligmosomoides polygyrus and Plasmodium chabaudi co-infection. Using newly established triple cytokine reporter mice (Il4^gfpIfng^yfpIl17a^FP635), we demonstrated that Il4^gfp+ Th2 cells purified from in vitro cultures or isolated ex vivo from helminth-infected mice up-regulated IFNγ following adoptive transfer into Rag1^–/– mice infected with P. chabaudi. Functionally, Th2 cells that up-regulated IFNγ were transcriptionally re-wired and protected recipient mice from high parasitemia. Mechanistically, TCR stimulation and responsiveness to IL-12 and IFNγ, but not type I IFN, was required for optimal IFNγ production by Th2 cells. Finally, blockade of IL-12 and IFNγ during co-infection partially preserved anti-helminth Th2 responses. In summary, this study demonstrates that Th2 cells retain substantial plasticity with the ability to produce IFNγ during Plasmodium infection. Consequently, co-infection with Plasmodium spp. may contribute to the chronicity of helminth infection by reducing anti-helminth Th2 cells and converting them into IFNγ-secreting cells.

Approximately a third of the world’s population is burdened with chronic intestinal parasitic helminth infections, causing significant morbidities. Identifying the factors that contribute to the chronicity of infection is therefore essential. Co-infection with other pathogens, which is extremely common in helminth endemic areas, may contribute to the chronicity of helminth infections. In this study, we used a mouse model to test whether the immune responses to an intestinal helminth were impaired following malaria co-infection. These two pathogens induce very different immune responses, which, until recently, were thought to be opposing and non-interchangeable. This study identified that the immune cells required for anti-helminth responses are capable of changing their phenotype and providing protection against malaria. By identifying and blocking the factors that drive this change in phenotype, we can preserve anti-helminth immune responses during co-infection. Our studies provide fresh insight into how immune responses are altered during helminth and malaria co-infection.

Hybrid male sterility and genome-wide misexpression of male reproductive proteases

Hybrid male sterility is a common barrier to gene flow between species. Previous studies have posited a link between misregulation of spermatogenesis genes in interspecies hybrids and sterility. However, in the absence of fully fertile control hybrids, it is impossible to differentiate between misregulation associated with sterility vs. fast male gene regulatory evolution. Here, we differentiate between these two possibilities using a D. pseudoobscura species pair that experiences unidirectional hybrid sterility. We identify genes uniquely misexpressed in sterile hybrid male reproductive tracts via RNA-seq. The sterile male hybrids had more misregulated and more over or under expressed genes relative to parental species than the fertile male hybrids. Proteases were the only gene ontology class overrepresented among uniquely misexpressed genes, with four located within a previously identified hybrid male sterility locus. This result highlights the potential role of a previously unexplored class of genes in interspecific hybrid male sterility and speciation.

Chromosomal Copy Number Variation in Saccharomyces pastorianus Is Evidence for Extensive Genome Dynamics in Industrial Lager Brewing Strains

Lager brewing strains of Saccharomyces pastorianus are natural interspecific hybrids originating from the spontaneous hybridization of Saccharomyces cerevisiae and Saccharomyces eubayanus. Over the past 500 years, S. pastorianus has been domesticated to become one of the most important industrial microorganisms. Production of lager-type beers requires a set of essential phenotypes, including the ability to ferment maltose and maltotriose at low temperature, the production of flavors and aromas, and the ability to flocculate. Understanding of the molecular basis of complex brewing-related phenotypic traits is a prerequisite for rational strain improvement. While genome sequences have been reported, the variability and dynamics of S. pastorianus genomes have not been investigated in detail. Here, using deep sequencing and chromosome copy number analysis, we showed that S. pastorianus strain CBS1483 exhibited extensive aneuploidy. This was confirmed by quantitative PCR and by flow cytometry. As a direct consequence of this aneuploidy, a massive number of sequence variants was identified, leading to at least 1,800 additional protein variants in S. pastorianus CBS1483. Analysis of eight additional S. pastorianus strains revealed that the previously defined group I strains showed comparable karyotypes, while group II strains showed large interstrain karyotypic variability. Comparison of three strains with nearly identical genome sequences revealed substantial chromosome copy number variation, which may contribute to strain-specific phenotypic traits. The observed variability of lager yeast genomes demonstrates that systematic linking of genotype to phenotype requires a three-dimensional genome analysis encompassing physical chromosomal structures, the copy number of individual chromosomes or chromosomal regions, and the allelic variation of copies of individual genes.

Genome-wide analysis of human global and transcription-coupled excision repair of UV damage at single-nucleotide resolution

Hu et al. developed a method for genome-wide mapping of DNA excision repair named XR-seq (excision repair sequencing) and used it to produce stranded, nucleotide-resolution maps of repair of two UV-induced DNA damages in human cells. XR-seq and the resulting repair maps will facilitate studies of the effects of genomic location, chromatin context, transcription, and replication on DNA repair in human cells.

We developed a method for genome-wide mapping of DNA excision repair named XR-seq (excision repair sequencing). Human nucleotide excision repair generates two incisions surrounding the site of damage, creating an ∼30-mer. In XR-seq, this fragment is isolated and subjected to high-throughput sequencing. We used XR-seq to produce stranded, nucleotide-resolution maps of repair of two UV-induced DNA damages in human cells: cyclobutane pyrimidine dimers (CPDs) and (6-4) pyrimidine–pyrimidone photoproducts [(6-4)PPs]. In wild-type cells, CPD repair was highly associated with transcription, specifically with the template strand. Experiments in cells defective in either transcription-coupled excision repair or general excision repair isolated the contribution of each pathway to the overall repair pattern and showed that transcription-coupled repair of both photoproducts occurs exclusively on the template strand. XR-seq maps capture transcription-coupled repair at sites of divergent gene promoters and bidirectional enhancer RNA (eRNA) production at enhancers. XR-seq data also uncovered the repair characteristics and novel sequence preferences of CPDs and (6-4)PPs. XR-seq and the resulting repair maps will facilitate studies of the effects of genomic location, chromatin context, transcription, and replication on DNA repair in human cells.

Impact of traumatic brain injury on sleep structure, electrocorticographic activity and transcriptome in mice

Traumatic brain injury (TBI), including mild TBI (mTBI), is importantly associated with vigilance and sleep complaints. Because sleep is required for learning, plasticity and recovery, we here evaluated the bidirectional relationship between mTBI and sleep with two specific objectives: (1) Test that mTBI rapidly impairs sleep-wake architecture and the dynamics of the electrophysiological marker of sleep homeostasis (i.e., non-rapid eye movement sleep delta (1-4Hz) activity); (2) evaluate the impact of sleep loss following mTBI on the expression of plasticity markers that have been linked to sleep homeostasis and on genome-wide gene expression. A closed-head injury model was used to perform a 48h electrocorticographic (ECoG) recording in mice submitted to mTBI or Sham surgery. mTBI was found to immediately decrease the capacity to sustain long bouts of wakefulness as well as the amplitude of the time course of ECoG delta activity during wakefulness. Significant changes in ECoG spectral activity during wakefulness, non-rapid eye movement and rapid eye movement sleep were observed mainly on the second recorded day. A second experiment was performed to measure gene expression in the cerebral cortex and hippocampus after a mTBI followed either by two consecutive days of 6h sleep deprivation (SD) or of undisturbed behavior (quantitative PCR and next-generation sequencing). mTBI modified the expression of genes involved in immunity, inflammation and glial function (e.g., chemokines, glial markers) and SD changed that of genes linked to circadian rhythms, synaptic activity/neuronal plasticity, neuroprotection and cell death and survival. SD appeared to affect gene expression in the cerebral cortex more importantly after mTBI than Sham surgery including that of the astrocytic marker Gfap, which was proposed as a marker of clinical outcome after TBI. Interestingly, SD impacted the hippocampal expression of the plasticity elements Arc and EfnA3 only after mTBI. Overall, our findings reveal alterations in spectral signature across all vigilance states in the first days after mTBI, and show that sleep loss post-mTBI reprograms the transcriptome in a brain area-specific manner and in a way that could be deleterious to brain recovery.

A CASQ1 founder mutation in three Italian families with protein aggregate myopathy and hyperCKaemia

BACKGROUND

Protein aggregate myopathies are increasingly recognised conditions characterised by a surplus of endogenous proteins. The molecular and mutational background for many protein aggregate myopathies has been clarified with the discovery of several underlying mutations. Familial idiopathic hyperCKaemia is a benign genetically heterogeneous condition with autosomal dominant features in a high proportion of cases.

METHODS

In 10 patients from three Italian families with autosomal dominant benign vacuolar myopathy and hyperCKaemia, we performed linkage analysis and exome sequencing as well as morphological and biochemical investigations.

RESULTS AND CONCLUSIONS

We show, by Sanger and exome sequencing, that the protein aggregate myopathy with benign evolution and muscle inclusions composed of excess CASQ1, affecting three Italian families, is due to the D244G heterozygous missense mutation in the CASQ1 gene. Investigation of microsatellite markers revealed a common haplotype in the three families indicating consanguinity and a founder effect. Results from immunocytochemistry, electron microscopy, biochemistry and transfected cell line investigations contribute to our understanding of pathogenetic mechanisms underlining this defect. The mutation is common to other Italian patients and is likely to share a founder effect with them. HyperCKaemia in the CASQ1-related myopathy is common and sometimes the sole overt manifestation. It is likely that CASQ1 mutations may remain undiagnosed if a muscle biopsy is not performed, and the condition could be more common than supposed.

Evolutionary Advantage Conferred by an Eukaryote-to-Eukaryote Gene Transfer Event in Wine Yeasts

Although an increasing number of horizontal gene transfers have been reported in eukaryotes, experimental evidence for their adaptive value is lacking. Here, we report the recent transfer of a 158-kb genomic region between Torulaspora microellipsoides and Saccharomyces cerevisiae wine yeasts or closely related strains. This genomic region has undergone several rearrangements in S. cerevisiae strains, including gene loss and gene conversion between two tandemly duplicated FOT genes encoding oligopeptide transporters. We show that FOT genes confer a strong competitive advantage during grape must fermentation by increasing the number and diversity of oligopeptides that yeast can utilize as a source of nitrogen, thereby improving biomass formation, fermentation efficiency, and cell viability. Thus, the acquisition of FOT genes has favored yeast adaptation to the nitrogen-limited wine fermentation environment. This finding indicates that anthropic environments offer substantial ecological opportunity for evolutionary diversification through gene exchange between distant yeast species.

Basil Downy Mildew (Peronospora belbahrii): Discoveries and Challenges Relative to Its Control

Basil (Ocimum spp.) is one of the most economically important and widely grown herbs in the world. Basil downy mildew, caused by Peronospora belbahrii, has become an important disease in sweet basil (O. basilicum) production worldwide in the past decade. Global sweet basil production is at significant risk to basil downy mildew because of the lack of genetic resistance and the ability of the pathogen to be distributed on infested seed. Controlling the disease is challenging and consequently many crops have been lost. In the past few years, plant breeding efforts have been made to identify germplasm that can be used to introduce downy mildew resistance genes into commercial sweet basils while ensuring that resistant plants have the correct phenotype, aroma, and tastes needed for market acceptability. Fungicide efficacy studies have been conducted to evaluate current and newly developed conventional and organic fungicides for its management with limited success. This review explores the current efforts and progress being made in understanding basil downy mildew and its control.

Identification of Candidate Adherent-Invasive E. coli Signature Transcripts by Genomic/Transcriptomic Analysis

Adherent-invasive Escherichia coli (AIEC) strains are detected more frequently within mucosal lesions of patients with Crohn’s disease (CD). The AIEC phenotype consists of adherence and invasion of intestinal epithelial cells and survival within macrophages of these bacteria in vitro. Our aim was to identify candidate transcripts that distinguish AIEC from non-invasive E. coli (NIEC) strains and might be useful for rapid and accurate identification of AIEC by culture-independent technology. We performed comparative RNA-Sequence (RNASeq) analysis using AIEC strain LF82 and NIEC strain HS during exponential and stationary growth. Differential expression analysis of coding sequences (CDS) homologous to both strains demonstrated 224 and 241 genes with increased and decreased expression, respectively, in LF82 relative to HS. Transition metal transport and siderophore metabolism related pathway genes were up-regulated, while glycogen metabolic and oxidation-reduction related pathway genes were down-regulated, in LF82. Chemotaxis related transcripts were up-regulated in LF82 during the exponential phase, but flagellum-dependent motility pathway genes were down-regulated in LF82 during the stationary phase. CDS that mapped only to the LF82 genome accounted for 747 genes. We applied an in silico subtractive genomics approach to identify CDS specific to AIEC by incorporating the genomes of 10 other previously phenotyped NIEC. From this analysis, 166 CDS mapped to the LF82 genome and lacked homology to any of the 11 human NIEC strains. We compared these CDS across 13 AIEC, but none were homologous in each. Four LF82 gene loci belonging to clustered regularly interspaced short palindromic repeats region (CRISPR)—CRISPR-associated (Cas) genes were identified in 4 to 6 AIEC and absent from all non-pathogenic bacteria. As previously reported, AIEC strains were enriched for pdu operon genes. One CDS, encoding an excisionase, was shared by 9 AIEC strains. Reverse transcription quantitative polymerase chain reaction assays for 6 genes were conducted on fecal and ileal RNA samples from 22 inflammatory bowel disease (IBD), and 32 patients without IBD (non-IBD). The expression of Cas loci was detected in a higher proportion of CD than non-IBD fecal and ileal RNA samples (p <0.05). These results support a comparative genomic/transcriptomic approach towards identifying candidate AIEC signature transcripts.

Identification of Prognostic Groups in High-Grade Serous Ovarian Cancer Treated with Platinum-Taxane Chemotherapy

Disseminated high-grade serous ovarian cancer (HGS-OvCa) is an aggressive disease treated with platinum and taxane combination therapy. While initial response can be favorable, the disease typically relapses and becomes resistant to treatment. As genomic alterations in HGS-OvCa are heterogeneous, identification of clinically meaningful molecular markers for outcome prediction is challenging. We developed a novel computational approach (PSFinder) that fuses transcriptomics and clinical data to identify HGS-OvCa prognostic subgroups for targeted treatment. Application of PSFinder to transcriptomics data from 180 HGS-OvCa patients treated with platinum-taxane therapy revealed 61 transcript isoforms that characterize two poor and one good survival-associated groups (P = 0.007). These groups were validated in eight independent data sets, including a prospectively collected ovarian cancer cohort. Two poor prognostic groups have distinct expression profiles and are characteristic by increased hypermethylation and stroma-related genes. Integration of the PSFinder signature and BRCA1/2 mutation status allowed even better stratification of HGS-OvCa patients' prognosis. The herein introduced novel and generally applicable computational approach can identify outcome-related subgroups and facilitate the development of precision medicine to overcome drug resistance. A limited set of biomarkers divides HGS-OvCa into three prognostic groups and predicts patients in need of targeted therapies.

Caenorhabditis elegans is a useful model for anthelmintic discovery

Parasitic nematodes infect one quarter of the world's population and impact all humans through widespread infection of crops and livestock. Resistance to current anthelmintics has prompted the search for new drugs. Traditional screens that rely on parasitic worms are costly and labour intensive and target-based approaches have failed to yield novel anthelmintics. Here, we present our screen of 67,012 compounds to identify those that kill the non-parasitic nematode Caenorhabditis elegans. We then rescreen our hits in two parasitic nematode species and two vertebrate models (HEK293 cells and zebrafish), and identify 30 structurally distinct anthelmintic lead molecules. Genetic screens of 19 million C. elegans mutants reveal those nematicides for which the generation of resistance is and is not likely. We identify the target of one lead with nematode specificity and nanomolar potency as complex II of the electron transport chain. This work establishes C. elegans as an effective and cost-efficient model system for anthelmintic discovery.

Next-generation sequencing reveals the biological significance of the N(2),3-ethenoguanine lesion in vivo

Etheno DNA adducts are a prevalent type of DNA damage caused by vinyl chloride (VC) exposure and oxidative stress. Etheno adducts are mutagenic and may contribute to the initiation of several pathologies; thus, elucidating the pathways by which they induce cellular transformation is critical. Although N(2),3-ethenoguanine (N(2),3-εG) is the most abundant etheno adduct, its biological consequences have not been well characterized in cells due to its labile glycosidic bond. Here, a stabilized 2'-fluoro-2'-deoxyribose analog of N(2),3-εG was used to quantify directly its genotoxicity and mutagenicity. A multiplex method involving next-generation sequencing enabled a large-scale in vivo analysis, in which both N(2),3-εG and its isomer 1,N(2)-ethenoguanine (1,N(2)-εG) were evaluated in various repair and replication backgrounds. We found that N(2),3-εG potently induces G to A transitions, the same mutation previously observed in VC-associated tumors. By contrast, 1,N(2)-εG induces various substitutions and frameshifts. We also found that N(2),3-εG is the only etheno lesion that cannot be repaired by AlkB, which partially explains its persistence. Both εG lesions are strong replication blocks and DinB, a translesion polymerase, facilitates the mutagenic bypass of both lesions. Collectively, our results indicate that N(2),3-εG is a biologically important lesion and may have a functional role in VC-induced or inflammation-driven carcinogenesis.

Duplication and Adaptive Evolution of a Key Centromeric Protein in Mimulus, a Genus with Female Meiotic Drive

The fundamental asymmetry of female meiosis creates an arena for genetic elements to compete for inclusion in the egg, promoting the selfish evolution of centromere variants that maximize their transmission to the future egg. Such "female meiotic drive" has been hypothesized to explain the paradoxically complex and rapidly evolving nature of centromeric DNA and proteins. Although theoretically widespread, few cases of active drive have been observed, thereby limiting the opportunities to directly assess the impact of centromeric drive on molecular variation at centromeres and binding proteins. Here, we characterize the molecular evolutionary patterns of CENH3, the centromere-defining histone variant, in Mimulus monkeyflowers, a genus with one of the few known cases of active centromere-associated female meiotic drive. First, we identify a novel duplication of CENH3 in diploid Mimulus, including in lineages with actively driving centromeres. Second, we demonstrate long-term adaptive evolution at several sites in the N-terminus of CENH3, a region with some meiosis-specific functions that putatively interacts with centromeric DNA. Finally, we infer that the paralogs evolve under different selective regimes; some sites in the N-terminus evolve under positive selection in the pro-orthologs or only one paralog (CENH3_B) and the paralogs exhibit significantly different patterns of polymorphism within populations. Our finding of long-term, adaptive evolution at CENH3 in the context of centromere-associated meiotic drive supports an antagonistic, coevolutionary battle for evolutionary dominance between centromeric DNA and binding proteins.

Next Generation Sequencing and Transcriptome Analysis Predicts Biosynthetic Pathway of Sennosides from Senna (Cassia angustifolia Vahl.), a Non-Model Plant with Potent Laxative Properties

Senna (Cassia angustifolia Vahl.) is a world's natural laxative medicinal plant. Laxative properties are due to sennosides (anthraquinone glycosides) natural products. However, little genetic information is available for this species, especially concerning the biosynthetic pathways of sennosides. We present here the transcriptome sequencing of young and mature leaf tissue of Cassia angustifolia using Illumina MiSeq platform that resulted in a total of 6.34 Gb of raw nucleotide sequence. The sequence assembly resulted in 42230 and 37174 transcripts with an average length of 1119 bp and 1467 bp for young and mature leaf, respectively. The transcripts were annotated using NCBI BLAST with 'green plant database (txid 33090)', Swiss Prot, Kyoto Encylcopedia of Genes & Genomes (KEGG), Cluster of Orthologous Gene (COG) and Gene Ontology (GO). Out of the total transcripts, 40138 (95.0%) and 36349 (97.7%) from young and mature leaf, respectively, were annotated by BLASTX against green plant database of NCBI. We used InterProscan to see protein similarity at domain level, a total of 34031 (young leaf) and 32077 (mature leaf) transcripts were annotated against the Pfam domains. All transcripts from young and mature leaf were assigned to 191 KEGG pathways. There were 166 and 159 CDS, respectively, from young and mature leaf involved in metabolism of terpenoids and polyketides. Many CDS encoding enzymes leading to biosynthesis of sennosides were identified. A total of 10,763 CDS differentially expressing in both young and mature leaf libraries of which 2,343 (21.7%) CDS were up-regulated in young compared to mature leaf. Several differentially expressed genes found functionally associated with sennoside biosynthesis. CDS encoding for many CYPs and TF families were identified having probable roles in metabolism of primary as well as secondary metabolites. We developed SSR markers for molecular breeding of senna. We have identified a set of putative genes involved in various secondary metabolite pathways, especially those related to the synthesis of sennosides which will serve as an important platform for public information about gene expression, genomics, and functional genomics in senna.

Complete Genome Sequencing of a Multidrug-Resistant and Human-Invasive Salmonella enterica Serovar Typhimurium Strain of the Emerging Sequence Type 213 Genotype

Salmonella enterica subsp. enterica serovar Typhimurium strain YU39 was isolated in 2005 in the state of Yucatán, Mexico, from a human systemic infection. The YU39 strain is representative of the multidrug-resistant emergent sequence type 213 (ST213) genotype. The YU39 complete genome is composed of a chromosome and seven plasmids.

Comparative transcriptome analysis reveals that the extracellular matrix receptor interaction contributes to the venous metastases of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common type of liver cancer in the world. Portal vein tumor thrombus (PVTT) is one of the most serious complications of HCC and is strongly correlated with a poor prognosis for HCC patients. However, the detailed mechanism of PVTT development remains to be explored. In this study, we present a large-scale transcriptome analysis, by RNA sequencing, of 11 patients diagnosed with HCC with PVTT. The dysregulated genes between HCC and PVTT suggested that the extracellular matrix receptor interaction was correlated with the venous metastases of HCC. Among all of the recurrent alternative splicing events, we identified exon 6 skipping of RPS24, which is likely to be a cancer driver. We also identified five common fusion genes between HCC and its corresponding PVTT samples, including ARID1A-GPATCH3, MDM1-NUP107, PTGES3-RARG, PRLR-TERT, and C9orf3-TMC1. All of these findings broaden our knowledge of PVTT development and may also contribute to the diagnosis and treatment of HCC patients with PVTT.

Genetics and genotype-phenotype correlations in Finnish patients with dilated cardiomyopathy

Genetic analysis among patients with dilated cardiomyopathy (DCM) is becoming an important part of clinical assessment, as it is in hypertrophic cardiomyopathy (HCM). The genetics of DCM is complex and therefore next-generation sequencing strategies are essential when providing genetic diagnostics. To achieve maximum yield, the diagnostic approach should include comprehensive clinical phenotyping combined with high-quality, high-coverage deep sequencing of DCM-associated genes and clinical variant classification as a basis for defining true yield in genetic testing. Our study has combined a novel sequencing strategy and clinical interpretation to analyse the yield and genotype–phenotype correlations among well-phenotyped Finnish DCM patients.

Aims

Despite our increased understanding of the genetic basis of dilated cardiomyopathy (DCM), the clinical utility and yield of clinically meaningful findings of comprehensive next-generation sequencing (NGS)-based genetic diagnostics in DCM has been poorly described. We utilized a high-quality oligonucleotide-selective sequencing (OS-Seq)-based targeted sequencing panel to investigate the genetic landscape of DCM in Finnish population and to evaluate the utility of OS-Seq technology as a novel comprehensive diagnostic tool.

Methods and results

Using OS-Seq, we targeted and sequenced the coding regions and splice junctions of 101 genes associated with cardiomyopathies in 145 unrelated Finnish patients with DCM. We developed effective bioinformatic variant filtering strategy and implemented strict variant classification scheme to reveal diagnostic yield and genotype–phenotype correlations. Implemented OS-Seq technology provided high coverage of the target region (median coverage 410× and 99.42% of the nucleotides were sequenced at least 15× read depth). Diagnostic yield was 35.2% (familial 47.6% and sporadic 25.6%, P = 0.004) when both pathogenic and likely pathogenic variants are considered as disease causing. Of these, 20 (53%) were titin (TTN) truncations (non-sense and frameshift) affecting all TTN transcripts. TTN truncations accounted for 20.6% and 14.6% of the familial and sporadic DCM cases, respectively.

Conclusion

Panel-based, high-quality NGS enables high diagnostic yield especially in the familial form of DCM, and bioinformatic variant filtering is a reliable step in the process of interpretation of genomic data in a clinical setting.

Comparative Transcriptome Analysis of Cultivated and Wild Watermelon during Fruit Development

Watermelon [Citrullus lanatus (Thunb.) Matsum. & Nakai] is an important vegetable crop world-wide. Watermelon fruit quality is a complex trait determined by various factors such as sugar content, flesh color and flesh texture. Fruit quality and developmental process of cultivated and wild watermelon are highly different. To systematically understand the molecular basis of these differences, we compared transcriptome profiles of fruit tissues of cultivated watermelon 97103 and wild watermelon PI296341-FR. We identified 2,452, 826 and 322 differentially expressed genes in cultivated flesh, cultivated mesocarp and wild flesh, respectively, during fruit development. Gene ontology enrichment analysis of these genes indicated that biological processes and metabolic pathways related to fruit quality such as sweetness and flavor were significantly changed only in the flesh of 97103 during fruit development, while those related to abiotic stress response were changed mainly in the flesh of PI296341-FR. Our comparative transcriptome profiling analysis identified critical genes potentially involved in controlling fruit quality traits including α-galactosidase, invertase, UDP-galactose/glucose pyrophosphorylase and sugar transporter genes involved in the determination of fruit sugar content, phytoene synthase, β-carotene hydroxylase, 9-cis-epoxycarotenoid dioxygenase and carotenoid cleavage dioxygenase genes involved in carotenoid metabolism, and 4-coumarate:coenzyme A ligase, cellulose synthase, pectinesterase, pectinesterase inhibitor, polygalacturonase inhibitor and α-mannosidase genes involved in the regulation of flesh texture. In addition, we found that genes in the ethylene biosynthesis and signaling pathway including ACC oxidase, ethylene receptor and ethylene responsive factor showed highly ripening-associated expression patterns, indicating a possible role of ethylene in fruit development and ripening of watermelon, a non-climacteric fruit. Our analysis provides novel insights into watermelon fruit quality and ripening biology. Furthermore, the comparative expression profile data we developed provides a valuable resource to accelerate functional studies in watermelon and facilitate watermelon crop improvement.

Impact of ceftiofur injection on gut microbiota and Escherichia coli resistance in pigs

Resistance to extended-spectrum cephalosporins (ESCs) is an important health concern. Here, we studied the impact of the administration of a long-acting form of ceftiofur on the pig gut microbiota and ESC resistance in Escherichia coli. Pigs were orally inoculated with an ESC-resistant E. coli M63 strain harboring a conjugative plasmid carrying a gene conferring resistance, bla CTX-M-1. On the same day, they were given or not a unique injection of ceftiofur. Fecal microbiota were studied using quantitative PCR analysis of the main bacterial groups and quantification of short-chain fatty acids. E. coli and ESC-resistant E. coli were determined by culture methods, and the ESC-resistant E. coli isolates were characterized. The copies of the bla CTX-M-1 gene were quantified. After ceftiofur injection, the main change in gut microbiota was the significant but transitory decrease in the E. coli population. Acetate and butyrate levels were significantly lower in the treated group. In all inoculated groups, E. coli M63 persisted in most pigs, and the bla CTX-M-1 gene was transferred to other E. coli. Culture and PCR results showed that the ceftiofur-treated group shed significantly more resistant strains 1 and 3 days after ESC injection. Thereafter, on most dates, there were no differences between the groups, but notably, one pig in the nontreated group regularly excreted very high numbers of ESC-resistant E. coli, probably leading to a higher contamination level in its pen. In conclusion, the use of ESCs, and also the presence of high-shedding animals, are important features in the spread of ESC resistance.

Complete Dosage Compensation in Anopheles stephensi and the Evolution of Sex-Biased Genes in Mosquitoes

Complete dosage compensation refers to hyperexpression of the entire X or Z chromosome in organisms with heterogametic sex chromosomes (XY male or ZW female) in order to compensate for having only one copy of the X or Z chromosome. Recent analyses suggest that complete dosage compensation, as in Drosophila melanogaster, may not be the norm. There has been no systematic study focusing on dosage compensation in mosquitoes. However, analysis of dosage compensation in Anopheles mosquitoes provides opportunities for evolutionary insights, as the X chromosome of Anopheles and that of its Dipteran relative, D. melanogaster formed independently from the same ancestral chromosome. Furthermore, Culicinae mosquitoes, including the Aedes genus, have homomorphic sex-determining chromosomes, negating the need for dosage compensation. Thus, Culicinae genes provide a rare phylogenetic context to investigate dosage compensation in Anopheles mosquitoes. Here, we performed RNA-seq analysis of male and female samples of the Asian malaria mosquito Anopheles stephensi and the yellow fever mosquito Aedes aegypti. Autosomal and X-linked genes in An. stephensi showed very similar levels of expression in both males and females, indicating complete dosage compensation. The uniformity of average expression levels of autosomal and X-linked genes remained when An. stephensi gene expression was normalized by that of their Ae. aegypti orthologs, strengthening the finding of complete dosage compensation in Anopheles. In addition, we comparatively analyzed the differentially expressed genes between adult males and adult females in both species, investigated sex-biased gene chromosomal distribution patterns in An. stephensi and provided three examples where gene duplications may have enabled the acquisition of sex-specific expression during mosquito evolution.

Draft Genome Sequence of Gerbil-Adapted Carcinogenic Helicobacter pylori Strain 7.13

We report here the draft genome sequence of Helicobacter pylori strain 7.13, a gerbil-adapted strain that causes gastric cancer in gerbils. Strain 7.13 is derived from clinical strain B128, isolated from a patient with a duodenal ulcer. This study reveals genes associated with the virulence of the strain.

Complete Genome Sequence of Listeria monocytogenes Strain DPC6895, a Serotype 1/2b Isolate from Bovine Raw Milk

Listeria monocytogenes is a foodborne pathogen and is the causative agent of listeriosis among humans and animals. The draft genome sequence of L. monocytogenes DPC6895, a serotype 1/2b strain isolated from the raw milk of a cow with subclinical bovine mastitis, is reported.

Determinants of nucleosome positioning and their influence on plant gene expression

Nucleosome positioning influences the access of transcription factors (TFs) to their binding sites and gene expression. Studies in plant, animal, and fungal models demonstrate similar nucleosome positioning patterns along genes and correlations between occupancy and expression. However, the relationships among nucleosome positioning, cis-regulatory element accessibility, and gene expression in plants remain undefined. Here we showed that plant nucleosome depletion occurs on specific 6-mer motifs and this sequence-specific nucleosome depletion is predictive of expression levels. Nucleosome-depleted regions in Arabidopsis thaliana tend to have higher G/C content, unlike yeast, and are centered on specific G/C-rich 6-mers, suggesting that intrinsic sequence properties, such as G/C content, cannot fully explain plant nucleosome positioning. These 6-mer motif sites showed higher DNase I hypersensitivity and are flanked by strongly phased nucleosomes, consistent with known TF binding sites. Intriguingly, this 6-mer-specific nucleosome depletion pattern occurs not only in promoter but also in genic regions and is significantly correlated with higher gene expression level, a phenomenon also found in rice but not in yeast. Among the 6-mer motifs enriched in genes responsive to treatment with the defense hormone jasmonate, there are no significant changes in nucleosome occupancy, suggesting that these sites are potentially preconditioned to enable rapid response without changing chromatin state significantly. Our study provides a global assessment of the joint contribution of nucleosome occupancy and motif sequences that are likely cis-elements to the control of gene expression in plants. Our findings pave the way for further understanding the impact of chromatin state on plant transcriptional regulatory circuits.

Phylogenetically typing bacterial strains from partial SNP genotypes observed from direct sequencing of clinical specimen metagenomic data

We describe an approach for genotyping bacterial strains from low coverage genome datasets, including metagenomic data from complex samples. Sequence reads from unknown samples are aligned to a reference genome where the allele states of known SNPs are determined. The Whole Genome Focused Array SNP Typing (WG-FAST) pipeline can identify unknown strains with much less read data than is needed for genome assembly. To test WG-FAST, we resampled SNPs from real samples to understand the relationship between low coverage metagenomic data and accurate phylogenetic placement. WG-FAST can be downloaded from https://github.com/jasonsahl/wgfast.

Genome-wide analysis of the ATP-binding cassette (ABC) transporter gene family in sea lamprey and Japanese lamprey

Background

Lampreys are extant representatives of the jawless vertebrate lineage that diverged from jawed vertebrates around 500 million years ago. Lamprey genomes contain information crucial for understanding the evolution of gene families in vertebrates. The ATP-binding cassette (ABC) gene family is found from prokaryotes to eukaryotes. The recent availability of two lamprey draft genomes from sea lamprey Petromyzon marinus and Japanese lamprey Lethenteron japonicum presents an opportunity to infer early evolutionary events of ABC genes in vertebrates.

Results

We conducted a genome-wide survey of the ABC gene family in two lamprey draft genomes. A total of 37 ABC transporters were identified and classified into seven subfamilies; namely seven ABCA genes, 10 ABCB genes, 10 ABCC genes, three ABCD genes, one ABCE gene, three ABCF genes, and three ABCG genes. The ABCA subfamily has expanded from three genes in sea squirts, seven and nine in lampreys and zebrafish, to 13 and 16 in human and mouse. Conversely, the multiple copies of ABCB1-, ABCG1-, and ABCG2-like genes found in sea squirts have contracted in the other species examined. ABCB2 and ABCB3 seem to be new additions in gnathostomes (not in sea squirts or lampreys), which coincides with the emergence of the gnathostome-specific adaptive immune system. All the genes in the ABCD, ABCE and ABCF subfamilies were conserved and had undergone limited duplication and loss events. In the sea lamprey transcriptomes, the ABCE and ABCF gene subfamilies were ubiquitously and highly expressed in all tissues while the members in other gene subfamilies were differentially expressed.

Conclusions

Thirteen more lamprey ABC transporter genes were identified in this study compared with a previous study. By concatenating the same gene sequences from the two lampreys, more full length sequences were obtained, which significantly improved both the assignment of gene names and the phylogenetic trees compared with a previous analysis using partial sequences. The ABC gene subfamilies in chordates have undergone obvious expansion or contraction. The ABCA subfamily showed the highest gene expansion rate during chordate evolution. The evolution of ABC transporters in lampreys requires further evaluation because the present results are based on a draft genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1677-z) contains supplementary material, which is available to authorized users.

A primer on precision medicine informatics

In this review, we describe key components of a computational infrastructure for a precision medicine program that is based on clinical-grade genomic sequencing. Specific aspects covered in this review include software components and hardware infrastructure, reporting, integration into Electronic Health Records for routine clinical use and regulatory aspects. We emphasize informatics components related to reproducibility and reliability in genomic testing, regulatory compliance, traceability and documentation of processes, integration into clinical workflows, privacy requirements, prioritization and interpretation of results to report based on clinical needs, rapidly evolving knowledge base of genomic alterations and clinical treatments and return of results in a timely and predictable fashion. We also seek to differentiate between the use of precision medicine in germline and cancer.

Heterologous expression and transcript analysis of gibberellin biosynthetic genes of grasses reveals novel functionality in the GA3ox family

BACKGROUND

The gibberellin (GA) pathway plays a central role in the regulation of plant development, with the 2-oxoglutarate-dependent dioxygenases (2-ODDs: GA20ox, GA3ox, GA2ox) that catalyse the later steps in the biosynthetic pathway of particularly importance in regulating bioactive GA levels. Although GA has important impacts on crop yield and quality, our understanding of the regulation of GA biosynthesis during wheat and barley development remains limited. In this study we identified or assembled genes encoding the GA 2-ODDs of wheat, barley and Brachypodium distachyon and characterised the wheat genes by heterologous expression and transcript analysis.

RESULTS

The wheat, barley and Brachypodium genomes each contain orthologous copies of the GA20ox, GA3ox and GA2ox genes identified in rice, with the exception of OsGA3ox1 and OsGA2ox5 which are absent in these species. Some additional paralogs of 2-ODD genes were identified: notably, a novel gene in the wheat B genome related to GA3ox2 was shown to encode a GA 1-oxidase, named as TaGA1ox-B1. This enzyme is likely to be responsible for the abundant 1β-hydroxylated GAs present in developing wheat grains. We also identified a related gene in barley, located in a syntenic position to TaGA1ox-B1, that encodes a GA 3,18-dihydroxylase which similarly accounts for the accumulation of unusual GAs in barley grains. Transcript analysis showed that some paralogs of the different classes of 2-ODD were expressed mainly in a single tissue or at specific developmental stages. In particular, TaGA20ox3, TaGA1ox1, TaGA3ox3 and TaGA2ox7 were predominantly expressed in developing grain. More detailed analysis of grain-specific gene expression showed that while the transcripts of biosynthetic genes were most abundant in the endosperm, genes encoding inactivation and signalling components were more highly expressed in the seed coat and pericarp.

CONCLUSIONS

The comprehensive expression and functional characterisation of the multigene families encoding the 2-ODD enzymes of the GA pathway in wheat and barley will provide the basis for a better understanding of GA-regulated development in these species. This analysis revealed the existence of a novel, endosperm-specific GA 1-oxidase in wheat and a related GA 3,18-dihydroxylase enzyme in barley that may play important roles during grain expansion and development.