Combining metagenomics with metaproteomics and stable isotope probing reveals metabolic pathways used by a naturally occurring marine methylotroph

A variety of culture-independent techniques have been developed that can be used in conjunction with culture-dependent physiological and metabolic studies of key microbial organisms, in order to better understand how the activity of natural populations influences and regulates all major biogeochemical cycles. In this study, we combined DNA-stable isotope probing with metagenomics and metaproteomics to characterize an as yet uncultivated marine methylotroph that actively incorporated carbon from 13C-labeled methanol into biomass. By metagenomic sequencing of the heavy DNA, we retrieved virtually the whole genome of this bacterium and determined its metabolic potential. Through protein-stable isotope probing, the RuMP cycle was established as the main carbon assimilation pathway, and the classical methanol dehydrogenase-encoding gene mxaF, as well as three out of four identified xoxF homologues were found to be expressed. This proof-of-concept study is the first in which theculture-independent techniques of DNA- and protein-stable isotope probing have been used to characterize the metabolism of a naturally-ocurring Methylophaga-like bacterium in the marine environment (i.e. M. thiooxydans L4) and thus provides a powerful approach to access the genome and proteome of uncultivated microbes involved in key processes in the environment

Genome-wide characterization reveals complex interplay between TP53 and TP63 in response to genotoxic stress

In response to genotoxic stress the TP53 tumour suppressor activates target gene expression to induce cell cycle arrest or apoptosis depending on the extent of DNA damage. These canonical activities can be repressed by TP63 in normal stratifying epithelia to maintain proliferative capacity or drive proliferation of squamous cell carcinomas, where TP63 is frequently overexpressed/amplified. Here we use ChIP-sequencing, integrated with microarray analysis, to define the genome-wide interplay between TP53 and TP63 in response to genotoxic stress in normal cells. We reveal that TP53 and TP63 bind to overlapping, but distinct cistromes of sites through utilization of distinctive consensus motifs and that TP53 is constitutively bound to a number of sites. We demonstrate that cisplatin and adriamycin elicit distinct effects on TP53 and TP63 binding events, through which TP53 can induce or repress transcription of an extensive network of genes by direct binding and/or modulation of TP63 activity. Collectively, this results in a global TP53-dependent repression of cell cycle progression, mitosis and DNA damage repair concomitant with activation of anti-proliferative and pro-apoptotic canonical target genes. Further analyses reveal that in the absence of genotoxic stress TP63 plays an important role in maintaining expression of DNA repair genes, loss of which results in defective repair.

Characterization of nontypable Haemophilus influenzae isolates recovered from adult patients with underlying chronic lung disease reveals genotypic and phenotypic traits associated with persistent infection

Nontypable Haemophilus influenzae (NTHi) has emerged as an important opportunistic pathogen causing infection in adults suffering obstructive lung diseases. Existing evidence associates chronic infection by NTHi to the progression of the chronic respiratory disease, but specific features of NTHi associated with persistence have not been comprehensively addressed. To provide clues about adaptive strategies adopted by NTHi during persistent infection, we compared sequential persistent isolates with newly acquired isolates in sputa from six patients with chronic obstructive lung disease. Pulse field gel electrophoresis (PFGE) identified three patients with consecutive persistent strains and three with new strains. Phenotypic characterisation included infection of respiratory epithelial cells, bacterial self-aggregation, biofilm formation and resistance to antimicrobial peptides (AMP). Persistent isolates differed from new strains in showing low epithelial adhesion and inability to form biofilms when grown under continuous-flow culture conditions in microfermenters. Self-aggregation clustered the strains by patient, not by persistence. Increasing resistance to AMPs was observed for each series of persistent isolates; this was not associated with lipooligosaccharide decoration with phosphorylcholine or with lipid A acylation. Variation was further analyzed for the series of three persistent isolates recovered from patient 1. These isolates displayed comparable growth rate, natural transformation frequency and murine pulmonary infection. Genome sequencing of these three isolates revealed sequential acquisition of single-nucleotide variants in the AMP permease sapC, the heme acquisition systems hgpB, hgpC, hup and hxuC, the 3-deoxy-D-manno-octulosonic acid kinase kdkA, the long-chain fatty acid transporter ompP1, and the phosphoribosylamine glycine ligase purD. Collectively, we frame a range of pathogenic traits and a repertoire of genetic variants in the context of persistent infection by NTHi.

An unbiased genetic screen reveals the polygenic nature of the influenza virus anti-interferon response

UNLABELLED: Influenza A viruses counteract the cellular innate immune response at several steps, including blocking RIG I-dependent activation of interferon (IFN) transcription, interferon (IFN)-dependent upregulation of IFN-stimulated genes (ISGs), and the activity of various ISG products; the multifunctional NS1 protein is responsible for most of these activities. To determine the importance of other viral genes in the interplay between the virus and the host IFN response, we characterized populations and selected mutants of wild-type viruses selected by passage through non-IFN-responsive cells. We reasoned that, by allowing replication to occur in the absence of the selection pressure exerted by IFN, the virus could mutate at positions that would normally be restricted and could thus find new optimal sequence solutions. Deep sequencing of selected virus populations and individual virus mutants indicated that nonsynonymous mutations occurred at many phylogenetically conserved positions in nearly all virus genes. Most individual mutants selected for further characterization induced IFN and ISGs and were unable to counteract the effects of exogenous IFN, yet only one contained a mutation in NS1. The relevance of these mutations for the virus phenotype was verified by reverse genetics. Of note, several virus mutants expressing intact NS1 proteins exhibited alterations in the M1/M2 proteins and accumulated large amounts of deleted genomic RNAs but nonetheless replicated to high titers. This suggests that the overproduction of IFN inducers by these viruses can override NS1-mediated IFN modulation. Altogether, the results suggest that influenza viruses replicating in IFN-competent cells have tuned their complete genomes to evade the cellular innate immune system and that serial replication in non-IFN-responsive cells allows the virus to relax from these constraints and find a new genome consensus within its sequence space.

IMPORTANCE: In natural virus infections, the production of interferons leads to an antiviral state in cells that effectively limits virus replication. The interferon response places considerable selection pressure on viruses, and they have evolved a variety of ways to evade it. Although the influenza virus NS1 protein is a powerful interferon antagonist, the contributions of other viral genes to interferon evasion have not been well characterized. Here, we examined the effects of alleviating the selection pressure exerted by interferon by serially passaging influenza viruses in cells unable to respond to interferon. Viruses that grew to high titers had mutations at many normally conserved positions in nearly all genes and were not restricted to the NS1 gene. Our results demonstrate that influenza viruses have fine-tuned their entire genomes to evade the interferon response, and by removing interferon-mediated constraints, viruses can mutate at genome positions normally restricted by the interferon response.

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

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

"Genome-Wide Linkage in a Highly Consanguineous Pedigree Reveals Two Novel Loci on Chromosome 7 for Non-Syndromic Familial Premature Ovarian Failure"

Background: The human condition known as Premature Ovarian Failure (POF) is characterized by loss of ovarian function before the age of 40. A majority of POF cases are sporadic, but 10–15% are familial, suggesting a genetic origin of the disease. Although several causal mutations have been identified, the etiology of POF is still unknown for about 90% of the patients. Methodology/Principal Findings: We report a genome-wide linkage and homozygosity analysis in one large consanguineous Middle-Eastern POF-affected family presenting an autosomal recessive pattern of inheritance. We identified two regions with a LODmax of 3.26 on chromosome 7p21.1-15.3 and 7q21.3-22.2, which are supported as candidate regions by homozygosity mapping. Sequencing of the coding exons and known regulatory sequences of three candidate genes (DLX5, DLX6 and DSS1) included within the largest region did not reveal any causal mutations. Conclusions/Significance: We detect two novel POF-associated loci on human chromosome 7, opening the way to the identification of new genes involved in the control of ovarian development and function.

Haplotype analysis of vernalization loci in European barley germplasm reveals novel VRN-H1 alleles and a predominant winter VRN-H1/VRN-H2 multi-locus haplotype

In barley, variation in the requirement for vernalization (an extended period of low temperature before flowering can occur) is determined by the VRN-H1, -H2 and -H3 loci. In European cultivated germplasm, most variation in vernalization requirement is accounted for by alleles at VRN-H1 and VRN-H2 only, but the range of allelic variation is largely unexplored. Here we characterise VRN-H1 and VRN-H2 haplotypes in 429 varieties representing a large portion of the acreage sown to barley in Western Europe over the last 60 years. Analysis of genotype, intron I sequencing data and growth habit tests identified three novel VRN-H1 alleles and determined the most frequent VRN-H1 intron I rearrangements. Combined analysis of VRN-H1 and VRN-H2 alleles resulted in the classification of seventeen VRN-H1/VRN-H2 multi-locus haplotypes, three of which account for 79% of varieties. The molecular markers employed here represent powerful diagnostic tools for prediction of growth habit and assessment of varietal purity. These markers will also allow development of germplasm to test the behaviour of individual alleles with the aim of understanding the relationship between allelic variation and adaptation to specific agri-environments.

Metaproteomics analysis reveals the adaptation process for the chicken gut microbiota

The animal gastrointestinal tract houses a large microbial community, the gut microbiota, that confers many benefits to its host, such as protection from pathogens and provision of essential metabolites. Metagenomic approaches have defined the chicken fecal microbiota in other studies, but here, we wished to assess the correlation between the metagenome and the bacterial proteome in order to better understand the healthy chicken gut microbiota. Here, we performed high-throughput sequencing of 16S rRNA gene amplicons and metaproteomics analysis of fecal samples to determine microbial gut composition and protein expression. 16 rRNA gene sequencing analysis identified Clostridiales, Bacteroidaceae, and Lactobacillaceae species as the most abundant species in the gut. For metaproteomics analysis, peptides were generated by using the Fasp method and subsequently fractionated by strong anion exchanges. Metaproteomics analysis identified 3,673 proteins. Among the most frequently identified proteins, 380 proteins belonged to Lactobacillus spp., 155 belonged to Clostridium spp., and 66 belonged to Streptococcus spp. The most frequently identified proteins were heat shock chaperones, including 349 GroEL proteins, from many bacterial species, whereas the most abundant enzymes were pyruvate kinases, as judged by the number of peptides identified per protein (spectral counting). Gene ontology and KEGG pathway analyses revealed the functions and locations of the identified proteins. The findings of both metaproteomics and 16S rRNA sequencing analyses are discussed.

Genome-wide mapping of 5-hydroxymethylcytosine in three rice cultivars reveals its preferential localization in transcriptionally silent transposable element genes

5-Hydroxymethylcytosine (5hmC), a modified form of cytosine that is considered the sixth nucleobase in DNA, has been detected in mammals and is believed to play an important role in gene regulation. In this study, 5hmC modification was detected in rice by employing a dot-blot assay, and its levels was further quantified in DNA from different rice tissues using liquid chromatography-multistage mass spectrometry (LC-MS/MS/MS). The results showed large intertissue variation in 5hmC levels. The genome-wide profiles of 5hmC modification in three different rice cultivars were also obtained using a sensitive chemical labelling followed by a next-generation sequencing method. Thousands of 5hmC peaks were identified, and a comparison of the distributions of 5hmC among different rice cultivars revealed the specificity and conservation of 5hmC modification. The identified 5hmC peaks were significantly enriched in heterochromatin regions,and mainly located in transposable element (TE) genes, especially around retrotransposons. The correlation analysis of 5hmC and gene expression data revealed a close association between 5hmC and silent TEs. These findings provide a resource for plant DNA 5hmC epigenetic studies and expand our knowledge of 5hmC modification.

Characterisation of monotreme caseins reveals lineage-specific expansion of an ancestral casein locus in mammals

Using a milk-cell cDNA sequencing approach we characterised milk-protein sequences from two monotreme species, platypus (Ornithorhynchus anatinus) and echidna (Tachyglossus aculeatus) and found a full set of caseins and casein variants. The genomic organisation of the platypus casein locus is compared with other mammalian genomes, including the marsupial opossum and several eutherians. Physical linkage of casein genes has been seen in the casein loci of all mammalian genomes examined and we confirm that this is also observed in platypus. However, we show that a recent duplication of β-casein occurred in the monotreme lineage, as opposed to more ancient duplications of α-casein in the eutherian lineage, while marsupials possess only single copies of α- and β-caseins. Despite this variability, the close proximity of the main α- and β-casein genes in an inverted tail-tail orientation and the relative orientation of the more distant kappa-casein genes are similar in all mammalian genome sequences so far available. Overall, the conservation of the genomic organisation of the caseins indicates the early, pre-monotreme development of the fundamental role of caseins during lactation. In contrast, the lineage-specific gene duplications that have occurred within the casein locus of monotremes and eutherians but not marsupials, which may have lost part of the ancestral casein locus, emphasises the independent selection on milk provision strategies to the young, most likely linked to different developmental strategies. The monotremes therefore provide insight into the ancestral drivers for lactation and how these have adapted in different lineages.

H19-DMR allele-specific methylation analysis reveals epigenetic heterogeneity of CTCF binding site 6 but not of site 5 in head-and-neck carcinomas: A pilot case-control analysis

Aberrant methylation of seven potential binding sites of the CTCF factor in the differentially methylated region upstream of the H19 gene (H19-DMR) has been suggested as critical for the regulation of IGF2 and H19 imprinted genes. In this study, we analyzed the allele-specific methylation pattern of CTCF binding sites 5 and 6 using methylationsensitive restriction enzyme PCR followed by RFLP analysis in matched tumoral and lymphocyte DNA from head-and-neck squamous cell carcinoma (HNSCC) patients, as well as in lymphocyte DNA from control individuals who were cancer-free. The monoallelic methylation pattern was maintained in CTCF binding site 5 in 22 heterozygous out of 91 samples analyzed. Nevertheless, a biallelic methylation pattern was detected in CTCF binding site 6 in a subgroup of HNSCC patients as a somatic acquired feature of tumor cells. An atypical biallelic methylation was also observed in both tumor and lymphocyte DNA from two patients, and at a high frequency in the control group (29 out of 64 informative controls). Additionally, we found that the C/T transition detected by HhaI RFLP suppressed one dinucleotide CpG in critical CTCF binding site 6, of a mutation showing polymorphic frequencies. Although a heterogeneous methylation pattern was observed after DNA sequencing modified by sodium bisulfite, the biallelic methylation pattern was confirmed in 9 out of 10 HNSCCs. These findings are likely to be relevant in the epigenetic regulation of the DMR, especially in pathological conditions in which the imprinting of IGF2 and H19 genes is disrupted.

Comparative genomics of two Leptospira interrogans serovars reveals novel insights into physiology and pathogenesis

Leptospira species colonize a significant proportion of rodent populations worldwide and produce life-threatening infections in accidental hosts, including humans. Complete genome sequencing of Leptospira interrogans serovar Copenhageni and comparative analysis with the available Leptospira interrogans serovar Lai genome reveal that despite overall genetic similarity there are significant structural differences, including a large chromosomal inversion and extensive variation in the number and distribution of insertion sequence elements. Genome sequence analysis elucidates many of the novel aspects of leptospiral physiology relating to energy metabolism, oxygen tolerance, two-component signal transduction systems, and mechanisms of pathogenesis. A broad array of transcriptional regulation proteins and two new families of afimbrial adhesins which contribute to host tissue colonization in the early steps of infection were identified. Differences in genes involved in the biosynthesis of lipopolysaccharide 0 side chains between the Copenhageni and Lai serovars were identified, offering an important starting point for the elucidation of the organism's complex polysaccharide surface antigens. Differences in adhesins and in lipopolysaccharide might be associated with the adaptation of serovars Copenhageni and Lai to different animal hosts. Hundreds of genes encoding surface-exposed lipoproteins and transmembrane outer membrane proteins were identified as candidates for development of vaccines for the prevention of leptospirosis.

FISH using a gag-like fragment probe reveals a common Ty3-gypsy-like retrotransposon in genome of Coffea species

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Gene Content Analysis of Sugarcane Public ESTs Reveals Thousands of Missing Coding-Genes and an Unexpected Pool of Grasses Conserved ncRNAs

Sugarcane is the most important crop for sugar industry and raw material for bioethanol. Here we present a quantitative analysis of the gene content from publicly available sugarcane ESTs. The current sugarcane EST collection sampled orthologs for ~58 % of the closely-related sorghum proteome, suggesting that more than 10,000 sugarcane coding-genes remain undiscovered. Moreover the existence of more than 2,000 ncRNAs conserved between sugarcane and sorghum was revealed, among which over 500 are also detected in rice, supporting the existence of hundreds of conserved ncRNAs in grasses. New efforts towards sugarcane transcriptome sequencing were needed to sample the missing coding-genes as well as to expand the catalog of ncRNAs. © 2012 Springer Science+Business Media, LLC.

Genome-Wide Analysis Reveals Selection for Important Traits in Domestic Horse Breeds

Intense selective pressures applied over short evolutionary time have resulted in homogeneity within, but substantial variation among, horse breeds. Utilizing this population structure, 744 individuals from 33 breeds, and a 54,000 SNP genotyping array, breed-specific targets of selection were identified using an FST-based statistic calculated in 500-kb windows across the genome. A 5.5-Mb region of ECA18, in which the myostatin (MSTN) gene was centered, contained the highest signature of selection in both the Paint and Quarter Horse. Gene sequencing and histological analysis of gluteal muscle biopsies showed a promoter variant and intronic SNP of MSTN were each significantly associated with higher Type 2B and lower Type 1 muscle fiber proportions in the Quarter Horse, demonstrating a functional consequence of selection at this locus. Signatures of selection on ECA23 in all gaited breeds in the sample led to the identification of a shared, 186-kb haplotype including two doublesex related mab transcription factor genes (DMRT2 and 3). The recent identification of a DMRT3 mutation within this haplotype, which appears necessary for the ability to perform alternative gaits, provides further evidence for selection at this locus. Finally, putative loci for the determination of size were identified in the draft breeds and the Miniature horse on ECA11, as well as when signatures of selection surrounding candidate genes at other loci were examined. This work provides further evidence of the importance of MSTN in racing breeds, provides strong evidence for selection upon gait and size, and illustrates the potential for population-based techniques to find genomic regions driving important phenotypes in the modern horse. © 2013 Petersen et al.