The Microbiome in Bronchial Biopsies from Smokers and Ex-Smokers with Stable COPD - A Metatranscriptomic Approach.

Current knowledge about the respiratory microbiome is mainly based on 16S ribosomal RNA gene sequencing. Newer sequencing approaches, such as metatranscriptomics, offer the technical ability to measure the viable microbiome response to environmental conditions such as smoking as well as to explore its functional role by investigating host-microbiome interactions. However, knowledge about its feasibility in respiratory microbiome research, especially in lung biopsies, is still very limited. RNA sequencing was performed in bronchial biopsies from clinically stable smokers (n = 5) and ex-smokers (n = 6) with COPD not using (inhaled) steroids. The Trinity assembler was used to assemble non-human reads in order to allow unbiased taxonomical and microbial transcriptional analyses. Subsequently, host-microbiome interactions were analyzed based on associations with host transcriptomic data. Ultra-low levels of microbial mass (0.009%) were identified in the RNA-seq data. Overall, no differences were identified in microbiome diversity or transcriptional profiles of microbial communities or individual microbes between COPD smokers and ex-smokers in the initial test dataset as well as a larger replication dataset. We identified an upregulated host gene set, related to the simultaneous presence of Bradyrhizobium, Roseomonas, Brevibacterium.spp., which were related to PERK-mediated unfolded protein response (UPR) and expression of the microRNA-155-5p. Our results show that metatranscriptomic profiling in bronchial biopsy samples from stable COPD patients yields ultra-low levels of microbial mass. Further, this study illustrates the potential of using transcriptional profiling of the host and microbiome to gain more insight into their interaction in the airways.

Population specific analysis of Yakut exomes.

We studied the genetic diversity of the Yakut population using exome sequencing. We performed comparative analysis of the Yakut population and the populations that are included in the "1000 Genomes" project and we identified the alleles specific to the Yakut population. We showed, that the Yakuts population is a separate cluster between Europeans and East Asians.

Using a priori knowledge to align sequencing reads to their exact genomic position.

The use of a priori knowledge in the alignment of targeted sequencing data is investigated using computational experiments. Adapting a Needleman-Wunsch algorithm to incorporate the genomic position information from the targeted capture, we demonstrate that alignment can be done to just the target region of interest. When in addition use is made of direct string comparison, an improvement of up to a factor of 8 in alignment speed compared to the fastest conventional aligner (Bowtie) is obtained. This results in a total alignment time in targeted sequencing of around 7 min for aligning approximately 56 million captured reads. For conventional aligners such as Bowtie, BWA or MAQ, alignment to just the target region is not feasible as experiments show that this leads to an additional 88% SNP calls, the vast majority of which are false positives (≈ 92%).

Systematic generation of in vivo G protein-coupled receptor mutants in the rat.

G-protein-coupled receptors (GPCRs) constitute a large family of cell surface receptors that are involved in a wide range of physiological and pathological processes, and are targets for many therapeutic interventions. However, genetic models in the rat, one of the most widely used model organisms in physiological and pharmacological research, are largely lacking. Here, we applied N-ethyl-N-nitrosourea (ENU)-driven target-selected mutagenesis to generate an in vivo GPCR mutant collection in the rat. A pre-selected panel of 250 human GPCR homologs was screened for mutations in 813 rats, resulting in the identification of 131 non-synonymous mutations. From these, seven novel potential rat gene knockouts were established as well as 45 lines carrying missense mutations in various genes associated with or involved in human diseases. We provide extensive in silico modeling results of the missense mutations and show experimental data, suggesting loss-of-function phenotypes for several models, including Mc4r and Lpar1. Taken together, the approach used resulted not only in a set of novel gene knockouts, but also in allelic series of more subtle amino acid variants, similar as commonly observed in human disease. The mutants presented here may greatly benefit studies to understand specific GPCR function and support the development of novel therapeutic strategies.

Age-related gene and miRNA expression changes in airways of healthy individuals.

Knowledge on age-related miRNA changes in healthy individuals and their interaction with mRNAs is lacking. We studied age-related mRNA and miRNA expression changes and their interactions in normal airways. RNA and small RNA sequencing was performed on bronchial biopsies of 86 healthy individuals (age: 18-73) to determine age-related expression changes. Per age-related miRNA we determined the enrichment of age-related predicted targets and their correlation. We identified 285 age-related genes and 27 age-related miRNAs. Pathway enrichment showed that genes higher expressed with age were involved in synapse-related processes. Genes lower expressed with age were involved in cell cycle regulation, the immune system and DNA damage/repair. MiR-146a-5p, miR-146b-5p and miR-142-5p were lower expressed with increasing age and we found a significant enrichment for predicted targets of these miRNAs among genes that were higher expressed with age. The expression levels of the enriched predicted targets RIMS2 and IGSF1 were negatively correlated with both miR-146a-5p and miR-146b-5p. RIMS2 was present in the enriched process, i.e. positive regulation of synaptic transmission. In conclusion, genes decreased with ageing are involved in several of the ageing hallmarks. Genes higher expressed with ageing were involved in synapse-related processes, of which RIMS2 is potentially regulated by two age-related miRNAs.

Insights in dynamic kinome reprogramming as a consequence of MEK inhibition in MLL-rearranged AML.

Single kinase-targeted cancer therapies often failed prolonged responses because cancer cells bypass through alternative routes. In this study, high-throughput kinomic and proteomic approaches enabled to identify aberrant activity profiles in mixed lineage leukemia (MLL)-rearranged acute myeloid leukemia (AML) that defined druggable targets. This approach revealed impaired activity of proteins belonging to the mitogen-activated protein kinase (MAPK) and phosphatidylinositol 3-kinase (PI3K) pathway. Pharmacological druggable MAPK pathway targets tested in primary MLL-rearranged AML included MAPKK1/2 (MEK), cyclic AMP-responsive element-binding protein (CREB) and MAPK8/9 (JNK). MEK inhibition showed to severely decrease MLL-rearranged AML cell survival without showing cytotoxicity in normal controls, whereas inhibition of CREB and JNK failed to exhibit MLL selectivity. Exploring the working mechanism of MEK inhibition, we assessed proteome activity in response to MEK inhibition in THP-1. MAPK1/3 (Erk) phosphorylation was instantly decreased in concurrence with a sustained Akt/mammalian target of rapamycin (mTOR) phosphorylation that enabled a subpopulation of cells to survive MEK inhibition. After exhaustion of MEK inhibition the AML cells recovered via increased activity of vascular endothelial growth factor receptor-2 (VEGFR-2) and Erk proteins to resume their proliferative state. Combined MEK and VEGFR-2 inhibition strengthened the reduction in MLL-rearranged AML cell survival by blocking the Akt/mTOR and MAPK pathways simultaneously. The generation of insights in cancerous altered activity profiles and alternative escape mechanisms upon targeted therapy allows the rational design of novel combination strategies.

Population variability in Chironomus (Camptochironomus) species (Diptera, Nematocera) with a Holarctic distribution: evidence of mitochondrial gene flow.

Phylogenetic analysis of DNA sequences from mitochondrial (mt) genes (Cytochrome b and Cytochrome oxidase I) and one nuclear gene (globin 2b) was used for the investigation of Nearctic and Palearctic populations representing four Chironomus species of the subgenus Camptochironomus, namely C. biwaprimus, C. pallidivittatus, C. tentans sensu stricto and C. dilutus (the last two species constitute Holarctic C. tentans sensu lato). Phenograms constructed on the basis of mt sequences were not congruent with trees based on nuclear genes, or with morphological and cytological data. The mt tree divided the populations by continental region, rather than by the species groupings recognized by the other data sets. The incongruence is explained by mt gene flow resulting from hybridization between the sympatric species on each continent. Calculation of divergence times, based on the sequence data, suggest that C. tentans (s.l.) and C. pallidivittatus have both been in North America for about 2.5 My.

Phylogeny of the genus Chironomus (Diptera) inferred from DNA sequences of mitochondrial cytochrome b and cytochrome oxidase I.

Two mitochondrial genes, Cytochrome b (Cytb) and Cytochrome c oxidase subunit I (COI), have been used as phylogenetic markers in Chironomids. The nucleotide sequences of 685 bp from Cytb and 596 bp from COI have been determined for 36 Chironomus species from the Palearctic, or Holarctic, and Australasia. The concatenated sequence of 1281 bp from both genes was used to investigate the phylogenetic relationships among these species. The nucleotide sequence alignments were used for construction of phylogenetic trees based on maximum-parsimony and neighbor-joining methods. Both techniques produced similar phylogenies. Monophyly of the genus Chironomus is supported by a bootstrap value of 100% at the basal branch. Six clusters of species have been revealed with high bootstrap values supporting both monophyly of each cluster and the validity of the branching order within each cluster. Four species, C. circumdatus, C. nepeanensis, C. dorsalis, and C. crassiforceps, cannot be placed into any cluster. Cytological phylogenies were constructed using the same set of species, except for C. biwaprimus. These trees showed many similarities to that obtained from the mitochondrial (mt) sequence analysis, but also a number of significant differences. When compared with the tree constructed from the sequence of 23 species available for one of the globin genes, globin 2b (gb2b), there was better support for the mt tree than for the cytological trees. An intron, which varies in its occurrence and position in gb2b, was also investigated and the distribution of the introns supports the phylogenetic history of the genus Chironomus obtained with mt data. The differences observed in the cytological trees seem to be attributable more to the retention of the same chromosome banding sequence across several species, rather than convergent evolutionary events. An important question is the determination of the position of the subgenus Camptochironomus in relation to the representatives of the nominal subgenus Chironomus, since it has been suggested that this is a separate genus. The Camptochironomus species are internal to the trees and have arisen more recently than some of the species of the subgenus Chironomus, indicating that they are not sufficiently differentiated to be considered more than a subgenus.

Insertional polymorphism of a non-LTR mobile element (NLRCth1) in European populations of Chironomus riparius (Diptera, Chironomidae) as detected by transposon insertion display.

The midge Chironomus riparius is distributed all over the Palearctic region and is well characterized both at the morphological and cytogenetic levels. Here we describe a population study based on the insertional polymorphism of the retroposon NLRCth1, by means of a S-SAP (sequence-specific amplification polymorphism) derived technique (transposon insertion display; TID). While a previous study of allozyme polymorphism in Russian samples showed little variability, all the amplicons we identified are polymorphic. Genetic distances between 6 natural populations were calculated according to Nei and did not show a positive correlation with geographic distances. The genetic diversity detected among individuals of a given population was one order of magnitude higher than that among populations. However, the value of phi(ST) was significant (p < 0.001) and indicates that natural populations are more genetically differentiated than random samples of individuals.