Ecological and genetic interactions between cyanobacteria and viruses in a low-oxygen mat community inferred through metagenomics and metatranscriptomics.

Metagenomic and metatranscriptomic sequencing was conducted on cyanobacterial mats of the Middle Island Sinkhole (MIS), Lake Huron. Metagenomic data from 14 samples collected over 5 years were used to reconstruct genomes of two genotypes of a novel virus, designated PhV1 type A and PhV1 type B. Both viral genotypes encode and express nblA, a gene involved in degrading phycobilisomes, which are complexes of pigmented proteins that harvest light for photosynthesis. Phylogenetic analysis indicated that the viral-encoded nblA is derived from the host cyanobacterium, Phormidium MIS-PhA. The cyanobacterial host also has two complete CRISPR (clustered regularly interspaced short palindromic repeats) systems that serve as defence mechanisms for bacteria and archaea against viruses and plasmids. One 45 bp CRISPR spacer from Phormidium had 100% nucleotide identity to PhV1 type B, but this region was absent from PhV1 type A. Transcripts from PhV1 and the Phormidium CRISPR loci were detected in all six metatranscriptomic data sets (three during the day and three at night), indicating that both are transcriptionally active in the environment. These results reveal ecological and genetic interactions between viruses and cyanobacteria at MIS, highlighting the value of parallel analysis of viruses and hosts in understanding ecological interactions in natural communities.

LINE-1 Mediated Insertion into Poc1a (Protein of Centriole 1 A) Causes Growth Insufficiency and Male Infertility in Mice.

Skeletal dysplasias are a common, genetically heterogeneous cause of short stature that can result from disruptions in many cellular processes. We report the identification of the lesion responsible for skeletal dysplasia and male infertility in the spontaneous, recessive mouse mutant chagun. We determined that Poc1a, encoding protein of the centriole 1a, is disrupted by the insertion of a processed Cenpw cDNA, which is flanked by target site duplications, suggestive of a LINE-1 retrotransposon-mediated event. Mutant fibroblasts have impaired cilia formation and multipolar spindles. Male infertility is caused by defective spermatogenesis early in meiosis and progressive germ cell loss. Spermatogonial stem cell transplantation studies revealed that Poc1a is essential for normal function of both Sertoli cells and germ cells. The proliferative zone of the growth plate is small and disorganized because chondrocytes fail to re-align after cell division and undergo increased apoptosis. Poc1a and several other genes associated with centrosome function can affect the skeleton and lead to skeletal dysplasias and primordial dwarfisms. This mouse mutant reveals how centrosome dysfunction contributes to defects in skeletal growth and male infertility.

The IRE1α/XBP1s Pathway Is Essential for the Glucose Response and Protection of ß Cells.

Although glucose uniquely stimulates proinsulin biosynthesis in ß cells, surprisingly little is known of the underlying mechanism(s). Here, we demonstrate that glucose activates the unfolded protein response transducer inositol-requiring enzyme 1 alpha (IRE1α) to initiate X-box-binding protein 1 (Xbp1) mRNA splicing in adult primary ß cells. Using mRNA sequencing (mRNA-Seq), we show that unconventional Xbp1 mRNA splicing is required to increase and decrease the expression of several hundred mRNAs encoding functions that expand the protein secretory capacity for increased insulin production and protect from oxidative damage, respectively. At 2 wk after tamoxifen-mediated Ire1α deletion, mice develop hyperglycemia and hypoinsulinemia, due to defective ß cell function that was exacerbated upon feeding and glucose stimulation. Although previous reports suggest IRE1α degrades insulin mRNAs, Ire1α deletion did not alter insulin mRNA expression either in the presence or absence of glucose stimulation. Instead, ß cell failure upon Ire1α deletion was primarily due to reduced proinsulin mRNA translation primarily because of defective glucose-stimulated induction of a dozen genes required for the signal recognition particle (SRP), SRP receptors, the translocon, the signal peptidase complex, and over 100 other genes with many other intracellular functions. In contrast, Ire1α deletion in ß cells increased the expression of over 300 mRNAs encoding functions that cause inflammation and oxidative stress, yet only a few of these accumulated during high glucose. Antioxidant treatment significantly reduced glucose intolerance and markers of inflammation and oxidative stress in mice with ß cell-specific Ire1α deletion. The results demonstrate that glucose activates IRE1α-mediated Xbp1 splicing to expand the secretory capacity of the ß cell for increased proinsulin synthesis and to limit oxidative stress that leads to ß cell failure.

Genome-Wide Identification of Klebsiella pneumoniae Fitness Genes during Lung Infection.

UNLABELLED: Klebsiella pneumoniae is an urgent public health threat because of resistance to carbapenems, antibiotics of last resort against Gram-negative bacterial infections. Despite the fact that K. pneumoniae is a leading cause of pneumonia in hospitalized patients, the bacterial factors required to cause disease are poorly understood. Insertion site sequencing combines transposon mutagenesis with high-throughput sequencing to simultaneously screen thousands of insertion mutants for fitness defects during infection. Using the recently sequenced K. pneumoniae strain KPPR1 in a well-established mouse model of pneumonia, insertion site sequencing was performed on a pool of >25,000 transposon mutants. The relative fitness requirement of each gene was ranked based on the ratio of lung to inoculum read counts and concordance between insertions in the same gene. This analysis revealed over 300 mutants with at least a 2-fold fitness defect and 69 with defects ranging from 10- to >2,000-fold. Construction of 6 isogenic mutants for use in competitive infections with the wild type confirmed their requirement for lung fitness. Critical fitness genes included those for the synthesis of branched-chain and aromatic amino acids that are essential in mice and humans, the transcriptional elongation factor RfaH, and the copper efflux pump CopA. The majority of fitness genes were conserved among reference strains representative of diverse pathotypes. These results indicate that regulation of outer membrane components and synthesis of amino acids that are essential to its host are critical for K. pneumoniae fitness in the lung. IMPORTANCE: Klebsiella pneumoniae is a bacterium that commonly causes pneumonia in patients after they are admitted to the hospital. K. pneumoniae is becoming resistant to all available antibiotics, and when these infections spread to the bloodstream, over half of patients die. Since currently available antibiotics are failing, we must discover new ways to treat these infections. In this study, we asked what genes the bacterium needs to cause an infection, since the proteins encoded by these genes could be targets for new antibiotics. We identified over 300 genes that K. pneumoniae requires to grow in a mouse model of pneumonia. Many of the genes that we identified are found in K. pneumoniae isolates from throughout the world, including antibiotic-resistant forms. If new antibiotics could be made against the proteins that these genes encode, they may be broadly effective against K. pneumoniae.

Alternative splicing of human NT5E in cirrhosis and hepatocellular carcinoma produces a negative regulator of ecto-5'-nucleotidase (CD73).

Ecto-5'-nucleotidase (CD73), encoded by NT5E, is the major enzymatic source of extracellular adenosine. CD73 controls numerous pathophysiological responses and is a potential disease target, but its regulation is poorly understood. We examined NT5E regulation by alternative splicing. Genomic database analysis of human transcripts led us to identify NT5E-2, a novel splice variant that was expressed at low abundance in normal human tissues but was significantly up-regulated in cirrhosis and hepatocellular carcinoma (HCC). NT5E-2 encodes a shorter CD73 isoform we named CD73S. The presence of CD73S protein, which lacks 50 amino acids, was detected in HCC using an isoform-specific antibody. A noncanonical mouse mRNA, similar to human CD73S, was observed, but the corresponding protein was undetectable. The two human isoforms exhibited functional differences, such that ectopic expression of canonical CD73 (CD73L) in human HepG2 cells was associated with decreased expression of the proliferation marker Ki67, whereas CD73S expression did not have an effect on Ki67 expression. CD73S was glycosylated, catalytically inactive, unable to dimerize, and complexed intracellularly with the endoplasmic reticulum chaperone calnexin. Furthermore, CD73S complexed with CD73L and promoted proteasome-dependent CD73L degradation. The findings reveal species-specific CD73 regulation, with potential significance to cancer, fibrosis, and other diseases characterized by changes in CD73 expression and function.

Complete Genome Sequence of Klebsiella pneumoniae Strain ATCC 43816 KPPR1, a Rifampin-Resistant Mutant Commonly Used in Animal, Genetic, and Molecular Biology Studies.

Klebsiella pneumoniae is an urgent public health threat due to the spread of carbapenem-resistant strains causing serious, and frequently fatal, infections. To facilitate genetic, molecular, and immunological studies of this pathogen, we report the complete chromosomal sequence of a genetically tractable, prototypical strain used in animal models.

CD1d-unrestricted NKT cells are endowed with a hybrid function far superior than that of iNKT cells.

Invariant natural killer T (iNKT) cells to date represent the best example of cells known to have a hybrid function, representing both innate and adaptive immunity. Shared phenotypic similarities with NK cells together with a rapid response to a cytokine stimulus and a productive TCR engagement are the features that underline the hybrid nature of iNKT cells. Using these criteria, we provide molecular and functional evidence demonstrating that CD1d-independent (CD1d(ind)) NKT cells, a population of CD1d-unrestricted NKT cells, are endowed with a hybrid function far superior to that of iNKT cells: (i) an extensive shared program with NK cells, (ii) a closer Euclidian distance with NK cells, and (iii) the ability to respond to innate stimuli (Poly:IC) with cytotoxic potential in the same manner as NK cells identify a hybrid feature in CD1d(ind)NKT cells that truly fulfills the dual function of an NK and a T cell. Our finding that CD1d(ind)NKT cells are programmed to act like NK cells in response to innate signals while being capable of adaptive responses is unprecedented, and thus might reemphasize CD1d-unrestricted NKT cells as a subset of lymphocytes that could affect biological processes of antimicrobial and tumor immunity in a unique way.

Novel bioinformatics method for identification of genome-wide non-canonical spliced regions using RNA-Seq data.

SETTING: During endoplasmic reticulum (ER) stress, the endoribonuclease (RNase) Ire1α initiates removal of a 26 nt region from the mRNA encoding the transcription factor Xbp1 via an unconventional mechanism (atypically within the cytosol). This causes an open reading frame-shift that leads to altered transcriptional regulation of numerous downstream genes in response to ER stress as part of the unfolded protein response (UPR). Strikingly, other examples of targeted, unconventional splicing of short mRNA regions have yet to be reported. OBJECTIVE: Our goal was to develop an approach to identify non-canonical, possibly very short, splicing regions using RNA-Seq data and apply it to ER stress-induced Ire1α heterozygous and knockout mouse embryonic fibroblast (MEF) cell lines to identify additional Ire1α targets. RESULTS: We developed a bioinformatics approach called the Read-Split-Walk (RSW) pipeline, and evaluated it using two Ire1α heterozygous and two Ire1α-null samples. The 26 nt non-canonical splice site in Xbp1 was detected as the top hit by our RSW pipeline in heterozygous samples but not in the negative control Ire1α knockout samples. We compared the Xbp1 results from our approach with results using the alignment program BWA, Bowtie2, STAR, Exonerate and the Unix "grep" command. We then applied our RSW pipeline to RNA-Seq data from the SKBR3 human breast cancer cell line. RSW reported a large number of non-canonical spliced regions for 108 genes in chromosome 17, which were identified by an independent study. CONCLUSIONS: We conclude that our RSW pipeline is a practical approach for identifying non-canonical splice junction sites on a genome-wide level. We demonstrate that our pipeline can detect novel splice sites in RNA-Seq data generated under similar conditions for multiple species, in our case mouse and human.

Regulation of the human endogenous retrovirus K (HML-2) transcriptome by the HIV-1 Tat protein.

UNLABELLED: Approximately 8% of the human genome is made up of endogenous retroviral sequences. As the HIV-1 Tat protein activates the overall expression of the human endogenous retrovirus type K (HERV-K) (HML-2), we used next-generation sequencing to determine which of the 91 currently annotated HERV-K (HML-2) proviruses are regulated by Tat. Transcriptome sequencing of total RNA isolated from Tat- and vehicle-treated peripheral blood lymphocytes from a healthy donor showed that Tat significantly activates expression of 26 unique HERV-K (HML-2) proviruses, silences 12, and does not significantly alter the expression of the remaining proviruses. Quantitative reverse transcription-PCR validation of the sequencing data was performed on Tat-treated PBLs of seven donors using provirus-specific primers and corroborated the results with a substantial degree of quantitative similarity. IMPORTANCE: The expression of HERV-K (HML-2) is tightly regulated but becomes markedly increased following infection with HIV-1, in part due to the HIV-1 Tat protein. The findings reported here demonstrate the complexity of the genome-wide regulation of HERV-K (HML-2) expression by Tat. This work also demonstrates that although HERV-K (HML-2) proviruses in the human genome are highly similar in terms of DNA sequence, modulation of the expression of specific proviruses in a given biological situation can be ascertained using next-generation sequencing and bioinformatics analysis.

A bioinformatics approach reveals novel interactions of the OVOL transcription factors in the regulation of epithelial - mesenchymal cell reprogramming and cancer progression.

BACKGROUND: Mesenchymal to Epithelial Transition (MET) plasticity is critical to cancer progression, and we recently showed that the OVOL transcription factors (TFs) are critical regulators of MET. Results of that work also posed the hypothesis that the OVOLs impact MET in a range of cancers. We now test this hypothesis by developing a model, OVOL Induced MET (OI-MET), and sub-model (OI-MET-TF), to characterize differential gene expression in MET common to prostate cancer (PC) and breast cancer (BC). RESULTS: In the OI-MET model, we identified 739 genes differentially expressed in both the PC and BC models. For this gene set, we found significant enrichment of annotation for BC, PC, cancer, and MET, as well as regulation of gene expression by AP1, STAT1, STAT3, and NFKB1. Focusing on the target genes for these four TFs plus the OVOLs, we produced the OI-MET-TF sub-model, which shows even greater enrichment for these annotations, plus significant evidence of cooperation among these five TFs. Based on known gene/drug interactions, we prioritized targets in the OI-MET-TF network for follow-on analysis, emphasizing the clinical relevance of this work. Reflecting these results back to the OI-MET model, we found that binding motifs for the TF pair AP1/MYC are more frequent than expected and that the AP1/MYC pair is significantly enriched in binding in cancer models, relative to non-cancer models, in these promoters. This effect is seen in both MET models (solid tumors) and in non-MET models (leukemia). These results are consistent with our hypothesis that the OVOLs impact cancer susceptibility by regulating MET, and extend the hypothesis to include mechanisms not specific to MET. CONCLUSIONS: We find significant evidence of the OVOL, AP1, STAT1, STAT3, and NFKB1 TFs having important roles in MET, and more broadly in cancer. We prioritize known gene/drug targets for follow-up in the clinic, and we show that the AP1/MYC TF pair is a strong candidate for intervention.

H2A.B facilitates transcription elongation at methylated CpG loci.

H2A.B is a unique histone H2A variant that only exists in mammals. Here we found that H2A.B is ubiquitously expressed in major organs. Genome-wide analysis of H2A.B in mouse ES cells shows that H2A.B is associated with methylated DNA in gene body regions. Moreover, H2A.B-enriched gene loci are actively transcribed. One typical example is that H2A.B is enriched in a set of differentially methylated regions at imprinted loci and facilitates transcription elongation. These results suggest that H2A.B positively regulates transcription elongation by overcoming DNA methylation in the transcribed region. It provides a novel mechanism by which transcription is regulated at DNA hypermethylated regions.

SNPAAMapper: An efficient genome-wide SNP variant analysis pipeline for next-generation sequencing data.

UNLABELLED: Many NGS analysis tools focusing on read alignment and variant calling functions for exome sequencing data have been developed in recent years. However, publicly available tools dealing with the downstream analysis of genome-wide variants are fewer and have limited functionality. We developed SNPAAMapper, a novel variant analysis pipeline that can effectively classify variants by region (e.g. CDS, UTRs, intron, upstream, downstream), predict amino acid change type (e.g. synonymous, non-synonymous mutation), and prioritize mutation effects (e.g. CDS versus UTRs). Additional functionality afforded by our pipeline includes: checking variants at exon/intron junctions, customized homozygosity and allele frequency cutoff parameters, and annotation of known variants with dbSNP information, listing original and mutated amino acid sequences containing variants. The final result is reported in a spreadsheet format table containing all variant associated information and prioritized amino acids effects for investigators to examine. AVAILABILITY: Perl scripts and required input files are available on the web at http://www.ccmb.med.umich.edu/ccdu /SNPAAMapper.

Community transcriptomic assembly reveals microbes that contribute to deep-sea carbon and nitrogen cycling.

The deep ocean is an important component of global biogeochemical cycles because it contains one of the largest pools of reactive carbon and nitrogen on earth. However, the microbial communities that drive deep-sea geochemistry are vastly unexplored. Metatranscriptomics offers new windows into these communities, but it has been hampered by reliance on genome databases for interpretation. We reconstructed the transcriptomes of microbial populations from Guaymas Basin, in the deep Gulf of California, through shotgun sequencing and de novo assembly of total community RNA. Many of the resulting messenger RNA (mRNA) contiguous sequences contain multiple genes, reflecting co-transcription of operons, including those from dominant members. Also prevalent were transcripts with only limited representation (2.8 times coverage) in a corresponding metagenome, including a considerable portion (1.2 Mb total assembled mRNA sequence) with similarity (96%) to a marine heterotroph, Alteromonas macleodii. This Alteromonas and euryarchaeal marine group II populations displayed abundant transcripts from amino-acid transporters, suggesting recycling of organic carbon and nitrogen from amino acids. Also among the most abundant mRNAs were catalytic subunits of the nitrite oxidoreductase complex and electron transfer components involved in nitrite oxidation. These and other novel genes are related to novel Nitrospirae and have limited representation in accompanying metagenomic data. High throughput sequencing of 16S ribosomal RNA (rRNA) genes and rRNA read counts confirmed that Nitrospirae are minor yet widespread members of deep-sea communities. These results implicate a novel bacterial group in deep-sea nitrite oxidation, the second step of nitrification. This study highlights metatranscriptomic assembly as a valuable approach to study microbial communities.

Next-generation sequencing identifies the Danforth's short tail mouse mutation as a retrotransposon insertion affecting Ptf1a expression.

The semidominant Danforth's short tail (Sd) mutation arose spontaneously in the 1920s. The homozygous Sd phenotype includes severe malformations of the axial skeleton with an absent tail, kidney agenesis, anal atresia, and persistent cloaca. The Sd mutant phenotype mirrors features seen in human caudal malformation syndromes including urorectal septum malformation, caudal regression, VACTERL association, and persistent cloaca. The Sd mutation was previously mapped to a 0.9 cM region on mouse chromosome 2qA3. We performed Sanger sequencing of exons and intron/exon boundaries mapping to the Sd critical region and did not identify any mutations. We then performed DNA enrichment/capture followed by next-generation sequencing (NGS) of the critical genomic region. Standard bioinformatic analysis of paired-end sequence data did not reveal any causative mutations. Interrogation of reads that had been discarded because only a single end mapped correctly to the Sd locus identified an early transposon (ETn) retroviral insertion at the Sd locus, located 12.5 kb upstream of the Ptf1a gene. We show that Ptf1a expression is significantly upregulated in Sd mutant embryos at E9.5. The identification of the Sd mutation will lead to improved understanding of the developmental pathways that are misregulated in human caudal malformation syndromes.

Identification of cross-species shared transcriptional networks of diabetic nephropathy in human and mouse glomeruli.

Murine models are valuable instruments in defining the pathogenesis of diabetic nephropathy (DN), but they only partially recapitulate disease manifestations of human DN, limiting their utility. To define the molecular similarities and differences between human and murine DN, we performed a cross-species comparison of glomerular transcriptional networks. Glomerular gene expression was profiled in patients with early type 2 DN and in three mouse models (streptozotocin DBA/2, C57BLKS db/db, and eNOS-deficient C57BLKS db/db mice). Species-specific transcriptional networks were generated and compared with a novel network-matching algorithm. Three shared human-mouse cross-species glomerular transcriptional networks containing 143 (Human-DBA STZ), 97 (Human-BKS db/db), and 162 (Human-BKS eNOS(-/-) db/db) gene nodes were generated. Shared nodes across all networks reflected established pathogenic mechanisms of diabetes complications, such as elements of Janus kinase (JAK)/signal transducer and activator of transcription (STAT) and vascular endothelial growth factor receptor (VEGFR) signaling pathways. In addition, novel pathways not previously associated with DN and cross-species gene nodes and pathways unique to each of the human-mouse networks were discovered. The human-mouse shared glomerular transcriptional networks will assist DN researchers in selecting mouse models most relevant to the human disease process of interest. Moreover, they will allow identification of new pathways shared between mice and humans.

Current status and future perspectives for sequencing livestock genomes.

Only in recent years, the draft sequences for several agricultural animals have been assembled. Assembling an individual animal's entire genome sequence or specific region(s) of interest is increasingly important for agricultural researchers to perform genetic comparisons between animals with different performance. We review the current status for several sequenced agricultural species and suggest that next generation sequencing (NGS) technology with decreased sequencing cost and increased speed of sequencing can benefit agricultural researchers. By taking advantage of advanced NGS technologies, genes and chromosomal regions that are more labile to the influence of environmental factors could be pinpointed. A more long term goal would be addressing the question of how animals respond at the molecular and cellular levels to different environmental models (e.g. nutrition). Upon revealing important genes and gene-environment interactions, the rate of genetic improvement can also be accelerated. It is clear that NGS technologies will be able to assist animal scientists to efficiently raise animals and to better prevent infectious diseases so that overall costs of animal production can be decreased.

Disruption of RAB40AL function leads to Martin--Probst syndrome, a rare X-linked multisystem neurodevelopmental human disorder.

BACKGROUND AND AIM: Martin--Probst syndrome (MPS) is a rare X-linked disorder characterised by deafness, cognitive impairment, short stature and distinct craniofacial dysmorphisms, among other features. The authors sought to identify the causative mutation for MPS. METHODS AND RESULTS: Massively parallel sequencing in two affected, related male subjects with MPS identified a RAB40AL (also called RLGP) missense mutation (chrX:102,079,078-102,079,079AC→GA p.D59G; hg18). RAB40AL encodes a small Ras-like GTPase protein with one suppressor of cytokine signalling box. The p.D59G variant is located in a highly conserved region of the GTPase domain between ß-2 and ß-3 strands. Using RT-PCR, the authors show that RAB40AL is expressed in human fetal and adult brain and kidney, and adult lung, heart, liver and skeletal muscle. RAB40AL appears to be a primate innovation, with no orthologues found in mouse, Xenopus or zebrafish. Western analysis and fluorescence microscopy of GFP-tagged RAB40AL constructs from transiently transfected COS7 cells show that the D59G missense change renders RAB40AL unstable and disrupts its cytoplasmic localisation. CONCLUSIONS: This is the first study to show that mutation of RAB40AL is associated with a human disorder. Identification of RAB40AL as the gene mutated in MPS allows for further investigations into the molecular mechanism(s) of RAB40AL and its roles in diverse processes such as cognition, hearing and skeletal development.

Convergence of genetic influences in comorbidity.

BACKGROUND: Predisposition to complex diseases is explained in part by genetic variation, and complex diseases are frequently comorbid, consistent with pleiotropic genetic variation influencing comorbidity. Genome Wide Association (GWA) studies typically assess association between SNPs and a single-disease phenotype. Fisher meta-analysis combines evidence of association from single-disease GWA studies, assuming that each study is an independent test of the same hypothesis. The Rank Product (RP) method overcomes limitations posed by Fisher assumptions, though RP was not designed for GWA data. METHODS: We modified RP to accommodate GWA data, and we call it modRP. Using p-values output from GWA studies, we aggregate evidence for association between SNPs and related phenotypes. To assess significance, RP randomly samples the observed ranks to develop the null distribution of the RP statistic, and then places the observed RPs into the null distribution. ModRP eliminates the effect of linkage disequilibrium and controls for differences in power at tested SNPs, to meet RP assumptions in application to GWA data. RESULTS: After validating modRP based on both positive and negative control studies, we searched for pleiotropic influences on comorbid substance use disorders in a novel study, and found two SNPs to be significantly associated with comorbid cocaine, opium, and nicotine dependence. Placing these SNPs into biological context, we developed a protein network modeling the interaction of cocaine, nicotine, and opium with these variants. CONCLUSIONS: ModRP is a novel approach to identifying pleiotropic genetic influences on comorbid complex diseases. It can be used to assess association for related phenotypes where raw data is unavailable or inappropriate for analysis using other approaches. The method is conceptually simple and produces statistically significant, biologically relevant results.

Decade-long bacterial community dynamics in cystic fibrosis airways.

The structure and dynamics of bacterial communities in the airways of persons with cystic fibrosis (CF) remain largely unknown. We characterized the bacterial communities in 126 sputum samples representing serial collections spanning 8-9 y from six age-matched male CF patients. Sputum DNA was analyzed by bar-coded pyrosequencing of the V3-V5 hypervariable region of the 16S rRNA gene, defining 662 operational taxonomic units (OTUs) from >633,000 sequences. Bacterial community diversity decreased significantly over time in patients with typically progressive lung disease but remained relatively stable in patients with a mild lung disease phenotype. Antibiotic use, rather than patient age or lung function, was the primary driver of decreasing diversity. Interpatient variability in community structure exceeded intrapatient variability in serial samples. Antibiotic treatment was associated with pronounced shifts in community structure, but communities showed both short- and long-term resilience after antibiotic perturbation. There was a positive correlation between OTU occurrence and relative abundance, with a small number of persistent OTUs accounting for the greatest abundance. Significant changes in community structure, diversity, or total bacterial density at the time of pulmonary exacerbation were not observed. Despite decreasing community diversity in patients with progressive disease, total bacterial density remained relatively stable over time. These findings show the critical relationship between airway bacterial community structure, disease stage, and clinical state at the time of sample collection. These features are the key parameters with which to assess the complex ecology of the CF airway.

The NIH National Center for Integrative Biomedical Informatics (NCIBI).

The National Center for Integrative and Biomedical Informatics (NCIBI) is one of the eight NCBCs. NCIBI supports information access and data analysis for biomedical researchers, enabling them to build computational and knowledge models of biological systems to address the Driving Biological Problems (DBPs). The NCIBI DBPs have included prostate cancer progression, organ-specific complications of type 1 and 2 diabetes, bipolar disorder, and metabolic analysis of obesity syndrome. Collaborating with these and other partners, NCIBI has developed a series of software tools for exploratory analysis, concept visualization, and literature searches, as well as core database and web services resources. Many of our training and outreach initiatives have been in collaboration with the Research Centers at Minority Institutions (RCMI), integrating NCIBI and RCMI faculty and students, culminating each year in an annual workshop. Our future directions include focusing on the TranSMART data sharing and analysis initiative.