BMC Genomics - Latest Articles

Genomics-driven discovery of the pneumocandin biosynthetic gene cluster in the fungus Glarea lozoyensis

Li Chen — 2013-05-20T00:00:00Z

Background: The antifungal therapy caspofungin is a semi-synthetic derivative of pneumocandin B0, a lipohexapeptide produced by the fungus Glarea lozoyensis, and was the first member of the echinocandin class approved for human therapy. The nonribosomal peptide synthetase (NRPS)-polyketide synthases (PKS) gene cluster responsible for pneumocandin biosynthesis from G. lozoyensis has not been elucidated to date. In this study, we report the elucidation of the pneumocandin biosynthetic gene cluster by whole genome sequencing of the G. lozoyensis wild-type strain ATCC 20868. Results: The pneumocandin biosynthetic gene cluster contains a NRPS (GLNRPS4) and a PKS (GLPKS4) arranged in tandem, two cytochrome P450 monooxygenases, seven other modifying enzymes, and genes for L-homotyrosine biosynthesis, a component of the peptide core. Thus, the pneumocandin biosynthetic gene cluster is significantly more autonomous and organized than that of the recently characterized echinocandin B gene cluster. Disruption mutants of GLNRPS4 and GLPKS4 no longer produced the pneumocandins (A0 and B0), and the Deltaglnrps4 and Deltaglpks4 mutants lost antifungal activity against the human pathogenic fungus Candida albicans. In addition to pneumocandins, the G. lozoyensis genome encodes a rich repertoire of natural product-encoding genes including 24 PKSs, six NRPSs, five PKS-NRPS hybrids, two dimethylallyl tryptophan synthases, and 14 terpene synthases. Conclusions: Characterization of the gene cluster provides a blueprint for engineering new pneumocandin derivatives with improved pharmacological properties. Whole genome estimation of the secondary metabolite-encoding genes from G. lozoyensis provides yet another example of the huge potential for drug discovery from natural products from the fungal kingdom.

The role of microRNAs in the pathogenesis of MMPi-induced skin fibrodysplasia

Daniel Tonge — 2013-05-20T00:00:00Z

Background: Matrix metalloproteinases (MMPs) are a family of proteolytic enzymes involved in extracellular matrix (ECM) homeostasis. MMPs have been an attractive pharmacological target for a number of indications. However, development has been hampered by the propensity of compounds targeting these enzymes to cause connective-tissue pathologies. The broad-spectrum MMP-inhibitor (MMPi) AZM551248 has been shown to induce such effects in the dog. Histopathological changes were consistent with fibrodysplasia (FD), characterised by fibroblast proliferation and the deposition of collagen in the subcutaneous tissues. We conducted a time-course study administering 20mg/kg/day AZM551248 between 4 and 17 days. Cervical subcutaneous tissue and plasma were sampled during the time-course. miRNA expression profiles in subcutaneous skin specimens following the administration of AZM551248 were determined by high-throughput-sequencing. Results: An increasing number of miRNAs were differentially expressed compared with vehicle treated control animals as the study progressed. Several of these were members of the miR-200 family and were significantly attenuated in response to MMPi. As the severity of FD increased at the later time-points, other miRNAs associated with TGFbeta synthesis and regulation of the acute inflammatory response were modulated. Evidence indicative of epithelial to mesenchymal transition was present at all study time points. Receiver operator curve (ROC) analysis revealed that miR-21 expression in the cervical subcutaneous tissue was a sensitive and specific biomarker of FD incidence. Conclusions: Our data reveal significant perturbations in canine skin miRNA expression in response to MMPi administration. Furthermore, we have identified dysregulated miRNAs that are associated with processes relevant to the key histopathological events of MMPi-induced FD.

Integrative genomic analysis of CREB defines a critical role for transcription factor networks in mediating the fed/fasted switch in liver

Logan Everett — 2013-05-17T00:00:00Z

Background: Metabolic homeostasis in mammals critically depends on the regulation of fasting-induced genes by CREB in the liver. Previous genome-wide analysis has shown that only a small percentage of CREB target genes are induced in response to fasting-associated signaling pathways. The precise molecular mechanisms by which CREB specifically targets these genes in response to alternating hormonal cues remain to be elucidated. Results: We performed chromatin immunoprecipitation coupled to high-throughput sequencing of CREB in livers from both fasted and re-fed mice. In order to quantitatively compare the extent of CREB-DNA interactions genome-wide between these two physiological conditions we developed a novel, robust analysis method, termed the 'single sample independence' (SSI) test that greatly reduced the number of false-positive peaks. We found that CREB remains constitutively bound to its target genes in the liver regardless of the metabolic state. Integration of the CREB cistrome with expression microarrays of fasted and re-fed mouse livers and ChIP-seq data for additional transcription factors revealed that the gene expression switches between the two metabolic states are associated with co-localization of additional transcription factors at CREB sites. Conclusions: Our results support a model in which CREB is constitutively bound to thousands of target genes, and combinatorial interactions between DNA-binding factors are necessary to achieve the specific transcriptional response of the liver to fasting. Furthermore, our genome-wide analysis identifies thousands of novel CREB target genes in liver, and suggests a previously unknown role for CREB in regulating ER stress genes in response to nutrient influx.

A simple and reproducible breast cancer prognostic test

Luigi Marchionni — 2013-05-17T00:00:00Z

Background: A small number of prognostic and predictive tests based on gene expression are currently offered as reference laboratory tests. In contrast to such success stories, a number of flaws and errors have recently been identified in other genomic-based predictors and the success rate for developing clinically useful genomic signatures is low. These errors have led to widespread concerns about the protocols for conducting and reporting of computational research. As a result, a need has emerged for a template for reproducible development of genomic signatures that incorporates full transparency, data sharing and statistical robustness. Results: Here we present the first fully reproducible analysis of the data used to train and test MammaPrint, an FDA-cleared prognostic test for breast cancer based on a 70-gene expression signature. We provide all the software and documentation necessary for researchers to build and evaluate genomic classifiers based on these data. As an example of the utility of this reproducible research resource, we develop a simple prognostic classifier that uses only 16 genes from the MammaPrint signature and is equally accurate in predicting 5-year disease free survival. Conclusions: Our study provides a prototypic example for reproducible development of computational algorithms for learning prognostic biomarkers in the era of personalized medicine.

Frequent loss of lineages and deficient duplications accounted for low copy number of disease resistance genes in Cucurbitaceae

Xiao Lin — 2013-05-17T00:00:00Z

Background: The sequenced genomes of cucumber, melon and watermelon have relatively few R-genes, with 70, 75 and 55 copies only, respectively. The mechanism for low copy number of R-genes in Cucurbitaceae genomes remains unknown. Results: Manual annotation of R-genes in the sequenced genomes of Cucurbitaceae species showed that approximately half of them are pseudogenes. Comparative analysis of R-genes showed frequent loss of R-gene loci in different Cucurbitaceae species. Phylogenetic analysis, data mining and PCR cloning using degenerate primers indicated that Cucurbitaceae has limited number of R-gene lineages (subfamilies). Comparison between R-genes from Cucurbitaceae and those from poplar and soybean suggested frequent loss of R-gene lineages in Cucurbitaceae. Furthermore, the average number of R-genes per lineage in Cucurbitaceae species is approximately 1/3 that in soybean or poplar. Therefore, both loss of lineages and deficient duplications in extant lineages accounted for the low copy number of R-genes in Cucurbitaceae. No extensive chimeras of R-genes were found in any of the sequenced Cucurbitaceae genomes. Nevertheless, one lineage of R-genes from Trichosanthes kirilowii, a wild Cucurbitaceae species, exhibits chimeric structures caused by gene conversions, and may contain a large number of distinct R-genes in natural populations. Conclusions: Cucurbitaceae species have limited number of R-gene lineages and each genome harbors relatively few R-genes. The scarcity of R-genes in Cucurbitaceae species was due to frequent loss of R-gene lineages and infrequent duplications in extant lineages. The evolutionary mechanisms for large variation of copy number of R-genes in different plant species were discussed.

Genetic parameters and genome-wide association study of hyperpigmentation of the visceral peritoneum in chickens

Chenglong Luo — 2013-05-16T00:00:00Z

Background: Hyperpigmentation of the visceral peritoneum (HVP) has recently garnered much attention in the poultry industry because of the possible risk to the health of affected animals and the damage it causes to the appearance of commercial chicken carcasses. However, the heritable characters of HVP remain unclear. The objective of this study was to investigate the genetic parameters of HVP by genome-wide association study (GWAS) in chickens. Results: HVP was found to be influenced by genetic factors, with a heritability score of 0.33. HVP had positive genetic correlations with growth and carcass traits, such as leg muscle weight (rg = 0.34), but had negative genetic correlations with immune traits, such as the antibody response to Newcastle disease virus (rg = -0.42). The GWAS for HVP using 39,833 single nucleotide polymorphisms indicated the genetic factors associated with HVP displayed an additive effect rather than a dominance effect. In addition, we determined that three genomic regions, involving the 50.5--54.0 Mb region of chicken (Gallus gallus) chromosome 1 (GGA1), the 58.5--60.5 Mb region of GGA1, and the 10.5--12.0 Mb region of GGA20, were strongly associated (P < 6.28 x 10-7) with HVP in chickens. Variants in these regions explained >50% of additive genetic variance for HVP. This study also confirmed that expression of BMP7, which codes for a bone morphogenetic protein and is located in one of the candidate regions, was significantly higher in the visceral peritoneum of Huiyang Beard chickens with HVP than in that of chickens without pigmentation (P < 0.05). Conclusions: HVP is a quantitative trait with moderate heritability. Genomic variants resulting in HVP were identified on GGA1 and GGA20, and expression of the BMP7 gene appears to be upregulated in HVP-affected chickens. Findings from this study should be used as a basis for further functional validation of candidate genes involved in HVP.

Genomic and physiological variability within Group II (non-proteolytic) Clostridium botulinum

Sandra Stringer — 2013-05-16T00:00:00Z

Background: Clostridium botulinum is a group of four physiologically and phylogenetically distinct bacteria that produce botulinum neurotoxin. While studies have characterised variability between strains of Group I (proteolytic) C. botulinum, the genetic and physiological variability and relationships between strains within Group II (non-proteolytic) C. botulinum are not well understood. In this study the genome of Group II strain C. botulinum Eklund 17B (NRP) was sequenced and used to construct a whole genome DNA microarray. This was used in a comparative genomic indexing study to compare the relatedness of 43 strains of Group II C. botulinum (14 type B, 24 type E and 5 type F). These results were compared with characteristics determined from physiological tests. Results: Whole genome indexing showed that strains of Group II C. botulinum isolated from a wide variety of environments over more than 75 years clustered together indicating the genetic background of Group II C. botulinum is stable. Further analysis showed that strains forming type B or type F toxin are closely related with only toxin cluster genes targets being unique to either type. Strains producing type E toxin formed a separate subset. Carbohydrate fermentation tests supported the observation that type B and F strains form a separate subset to type E strains. All the type F strains and most of type B strains produced acid from amylopectin, amylose and glycogen whereas type E strains did not. However, these two subsets did not differ strongly in minimum growth temperature or maximum NaCl concentration for growth. No relationship was found between tellurite resistance and toxin type despite all the tested type B and type F strains carrying tehB, while the sequence was absent or diverged in all type E strains. Conclusions: Although Group II C. botulinum form a tight genetic group, genomic and physiological analysis indicates there are two distinct subsets within this group. All type B strains and type F strains are in one subset and all type E strains in the other.

Structural and functional annotation of the porcine immunome

Harry Dawson — 2013-05-15T00:00:00Z

Background: The domestic pig is known as an excellent model for human immunology and the two species share many pathogens. Susceptibility to infectious disease is one of the major constraints on swine performance, yet the structure and function of genes comprising the pig immunome are not well-characterized. The completion of the pig genome provides the opportunity to annotate the pig immunome, and compare and contrast pig and human immune systems. Results: The Immune Response Annotation Group (IRAG) used computational curation and manual annotation of the swine genome assembly 10.2 (Sscrofa10.2) to refine the currently available automated annotation of 1,369 immunity-related genes through sequence-based comparison to genes in other species. Within these genes, we annotated 3,472 transcripts. Annotation provided evidence for gene expansions in several immune response families, and identified artiodactyl-specific expansions in the cathelicidin and type 1 Interferon families. We found gene duplications for 18 genes, including 13 immune response genes and five non-immune response genes discovered in the annotation process. Manual annotation provided evidence for many new alternative splice variants and 8 gene duplications. Over 1,100 transcripts without porcine sequence evidence were detected using cross-species annotation. We used a functional approach to discover and accurately annotate porcine immune response genes. A co-expression clustering analysis of transcriptomic data from selected experimental infections or immune stimulations of blood, macrophages or lymph nodes identified a large cluster of genes that exhibited a correlated positive response upon infection across multiple pathogens or immune stimuli. Interestingly, this gene cluster (cluster 4) is enriched for known general human immune response genes, yet contains many un-annotated porcine genes. A phylogenetic analysis of the encoded proteins of cluster 4 genes showed that 15% exhibited an accelerated evolution as compared to 4.1% across the entire genome. Conclusions: This extensive annotation dramatically extends the genome-based knowledge of the molecular genetics and structure of a major portion of the porcine immunome. Our complementary functional approach using co-expression during immune response has provided new putative immune response annotation for over 500 porcine genes. Our phylogenetic analysis of this core immunome cluster confirms rapid evolutionary change in this set of genes, and that, as in other species, such genes are important components of the pig’s adaptation to pathogen challenge over evolutionary time. These comprehensive and integrated analyses increase the value of the porcine genome sequence and provide important tools for global analyses and data-mining of the porcine immune response.

Differential transcript isoform usage pre- and post-zygotic genome activation in zebrafish

Håvard Aanes — 2013-05-15T00:00:00Z

Background: Zebrafish embryos are transcriptionally silent until activation of the zygotic genome during the 10th cell cycle. Onset of transcription is followed by cellular and morphological changes involving cell speciation and gastrulation. Previous genome-wide surveys of transcriptional changes only assessed gene expression levels; however, recent studies have shown the necessity to map isoform-specific transcriptional changes. Here, we perform isoform discovery and quantification on transcriptome sequences from before and after zebrafish zygotic genome activation (ZGA). Results: We identify novel isoforms and isoform switches during ZGA for genes related to cell adhesion, pluripotency and DNA methylation. Isoform switching events include alternative splicing and changes in transcriptional start sites and in 3' untranslated regions. New isoforms are identified even for well-characterized genes such as pou5f1, sall4 and dnmt1. Genes involved in cell-cell interactions such as f11r and magi1 display isoform switches with alterations of coding sequences. We also detect over 1000 transcripts that acquire a longer 3' terminal exon when transcribed by the zygote compared to their maternal transcript counterparts. ChIP-sequencing data mapped onto skipped exon events reveal a correlation between histone H3K36 trimethylation peaks and skipped exons, suggesting epigenetic marks being part of alternative splicing regulation. Conclusions: The novel isoforms and isoform switches reported here include regulators of transcriptional, cellular and morphological changes taking place around ZGA. Our data display an array of isoform-related functional changes and represent a valuable resource complementary to existing early embryo transcriptomes.

Characterisation of a transcriptome to find sequence differences between two differentially migrating subspecies of the willow warbler Phylloscopus trochilus

Max Lundberg — 2013-05-14T00:00:00Z

Background: Animal migration requires adaptations in morphological, physiological and behavioural traits. Several of these traits have been shown to possess a strong heritable component in birds, but little is known about their genetic architecture. Here we used 454 sequencing of brain-derived transcriptomes from two differentially migrating subspecies of the willow warbler Phylloscopus trochilus to detect genes potentially underlying traits associated with migration. Results: The transcriptome sequencing resulted in 1.8 million reads following filtering steps. Most of the reads (84%) were successfully mapped to the genome of the zebra finch Taeniopygia gutatta. The mapped reads were situated within at least 12,101 predicted zebra finch genes, with the greatest sequencing depth in exons. Reads that were mapped to intergenic regions were generally located close to predicted genes and possibly located in uncharacterized untranslated regions (UTRs). Out of 85,000 single nucleotide polymorphisms (SNPs) with a minimum sequencing depth of eight reads from each of two subspecies-specific pools, only 55 showed high differentiation, confirming previous studies showing that most of the genetic variation is shared between the subspecies. Validation of a subset of the most highly differentiated SNPs using Sanger sequencing demonstrated that several of them also were differentiated between an independent set of individuals of each subspecies. These SNPs were clustered in two chromosome regions that are likely to be influenced by divergent selection between the subspecies and that could potentially be associated with adaptations to their different migratory strategies. Conclusions: Our study represents the first large-scale sequencing analysis aiming at detecting genes underlying migratory phenotypes in birds and provides new candidates for genes potentially involved in migration.