Insect mitochondrial genomics : implications for evolution and phylogeny

The mitochondrial (mt) genome is, to date, the most extensively studied genomic system in insects, outnumbering nuclear genomes tenfold and representing all orders versus very few. Phylogenomic analysis methods have been tested extensively, identifying compositional bias and rate variation, both within and between lineages, as the principal issues confronting accurate analyses. Major studies at both inter- and intraordinal levels have contributed to our understanding of phylogenetic relationships within many groups. Genome rearrangements are an additional data type for defining relationships, with rearrangement synapomorphies identified across multiple orders and at many different taxonomic levels. Hymenoptera and Psocodea have greatly elevated rates of rearrangement offering both opportunities and pitfalls for identifying rearrangement synapomorphies in each group. Finally, insects are model systems for studying aberrant mt genomes, including truncated tRNAs and multichromosomal genomes. Greater integration of nuclear and mt genomic studies is necessary to further our understanding of insect genomic evolution.

Genomic organisation, activity and distribution analysis of the microbial putrescine oxidase degradation pathway

The catalytic action of putrescine specific amine oxidases acting in tandem with 4-aminobutyraldehyde dehydrogenase is explored as a degradative pathway in Rhodococcus opacus. By limiting the nitrogen source, increased catalytic activity was induced leading to a coordinated response in the oxidative deamination of putrescine to 4-aminobutyraldehyde and subsequent dehydrogenation to 4-aminobutyrate. Isolating the dehydrogenase by ion exchange chromatography and gel filtration revealed that the enzyme acts principally on linear aliphatic aldehydes possessing an amino moiety. Michaelis-Menten kinetic analysis delivered a Michaelis constant (KM=0.014mM) and maximum rate (Vmax=11.2μmol/min/mg) for the conversion of 4-aminobutyraldehyde to 4-aminobutyrate. The dehydrogenase identified by MALDI-TOF mass spectrometric analysis (E value=0.031, 23% coverage) belongs to a functionally related genomic cluster that includes the amine oxidase, suggesting their association in a directed cell response. Key regulatory, stress and transport encoding genes have been identified, along with candidate dehydrogenases and transaminases for the further conversion of 4-aminobutyrate to succinate. Genomic analysis has revealed highly similar metabolic gene clustering among members of Actinobacteria, providing insight into putrescine degradation notably among Micrococcaceae, Rhodococci and Corynebacterium by a pathway that was previously uncharacterised in bacteria.

Genomic resources notes accepted 1 December 2013 - 31 January 2014

This article documents the public availability of (i) transcriptome sequence data, assembled and annotated contigs and unigenes, and BLAST hits from the Queensland fruit fly, Bactrocera tryoni; (ii) 75 single-nucleotide variants (SNVs) from 454 sequencing of reduced representation libraries for Phalangiidae harvestmen, Megabunus armatus, Megabunus vignai, Megabunus lesserti, and Rilaena triangularis; and (iii) expressed sequence tags from 454 sequencing of the lepidopterans Lymantria dispar and Lymantria monacha.

Promotion of Homologous Recombination and Genomic Stability by RAD51AP1 via RAD51 Recombinase Enhancement

Homologous recombination (HR) repairs chromosome damage and is indispensable for tumor suppression in humans. RAD51 mediates the DNA strand-pairing step in HR. RAD51 associated protein 1 (RAD51AP1) is a RAD51-interacting protein whose function has remained elusive. Knockdown of RAD51AP1 in human cells by RNA interference engenders sensitivity to different types of genotoxic stress, and RAD51AP1 is epistatic to the HR protein XRCC3. Moreover, RAD51AP1-depleted cells are impaired for the recombinational repair of a DNA double-strand break and exhibit chromatid breaks both spontaneously and upon DNA-damaging treatment. Purified RAD51AP1 binds both dsDNA and a D loop structure and, only when able to interact with RAD51, greatly stimulates the RAD51-mediated D loop reaction. Biochemical and cytological results show that RAD51AP1 functions at a step subsequent to the assembly of the RAD51-ssDNA nucleoprotein filament. Our findings provide evidence that RAD51AP1 helps maintain genomic integrity via RAD51 recombinase enhancement.

Integrated analysis of epigenomic and genomic changes by DNA methylation dependent mechanisms provides potential novel biomarkers for prostate cancer

Epigenetic silencing mediated by CpG methylation is a common feature of many cancers. Characterizing aberrant DNA methylation changes associated with tumor progression may identify potential prognostic markers for prostate cancer (PCa). We treated two PCa cell lines, 22Rv1 and DU-145 with the demethylating agent 5-Aza 2’–deoxycitidine (DAC) and global methylation status was analyzed by performing methylation-sensitive restriction enzyme based differential methylation hybridization strategy followed by genome-wide CpG methylation array profiling. In addition, we examined gene expression changes using a custom microarray. Gene Set Enrichment Analysis (GSEA) identified the most significantly dysregulated pathways. In addition, we assessed methylation status of candidate genes that showed reduced CpG methylation and increased gene expression after DAC treatment, in Gleason score (GS) 8 vs. GS6 patients using three independent cohorts of patients; the publically available The Cancer Genome Atlas (TCGA) dataset, and two separate patient cohorts. Our analysis, by integrating methylation and gene expression in PCa cell lines, combined with patient tumor data, identified novel potential biomarkers for PCa patients. These markers may help elucidate the pathogenesis of PCa and represent potential prognostic markers for PCa patients.

Characterization of 67 mitochondrial tRNA gene rearrangements in the Hymenoptera suggests that mitochondrial tRNA gene position is selectively neutral

We present entire sequences of two hymenopteran mitochondrial genomes and the major portion of three others. We combined these data with nine previously sequenced hymenopteran mitochondrial genomes. This allowed us to infer and analyze the evolution of the 67 mitochondrial gene rearrangements so far found in this order. All of these involve tRNA genes, whereas four also involve larger (protein-coding or ribosomal RNA) genes. We find that the vast majority of mitochondrial gene rearrangements are independently derived. A maximum of four of these rearrangements represent shared, derived organizations, whereas three are convergently derived. The remaining mitochondrial gene rearrangements represent new mitochondrial genome organizations. These data are consistent with the proposal that there are an enormous number of alternative mitochondrial genome organizations possible and that mitochondrial genome organization is, for the most part, selectively neutral. Nevertheless, some mitochondrial genes appear less mobile than others. Genes close to the noncoding region are generally more mobile but only marginally so. Some mitochondrial genes rearrange in a pattern consistent with the duplication/random loss model, but more mitochondrial genes move in a pattern inconsistent with this model. An increased rate of mitochondrial gene rearrangement is not tightly associated with the evolution of parasitism. Although parasitic lineages tend to have more mitochondrial gene rearrangements than nonparasitic lineages, there are exceptions (e.g., Orussus and Schlettererius). It is likely that only a small proportion of the total number of mitochondrial gene rearrangements that have occurred during the evolution of the Hymenoptera have been sampled in the present study.

Locality-sensitive hashing for protein classification

Determination of sequence similarity is a central issue in computational biology, a problem addressed primarily through BLAST, an alignment based heuristic which has underpinned much of the analysis and annotation of the genomic era. Despite their success, alignment-based approaches scale poorly with increasing data set size, and are not robust under structural sequence rearrangements. Successive waves of innovation in sequencing technologies – so-called Next Generation Sequencing (NGS) approaches – have led to an explosion in data availability, challenging existing methods and motivating novel approaches to sequence representation and similarity scoring, including adaptation of existing methods from other domains such as information retrieval. In this work, we investigate locality-sensitive hashing of sequences through binary document signatures, applying the method to a bacterial protein classification task. Here, the goal is to predict the gene family to which a given query protein belongs. Experiments carried out on a pair of small but biologically realistic datasets (the full protein repertoires of families of Chlamydia and Staphylococcus aureus genomes respectively) show that a measure of similarity obtained by locality sensitive hashing gives highly accurate results while offering a number of avenues which will lead to substantial performance improvements over BLAST..

Methods for extracting genomic DNA from whole blood samples : current perspectives

Deoxyribonucleic acid (DNA) extraction has considerably evolved since it was initially performed back in 1869. It is the first step required for many of the available downstream applications used in the field of molecular biology. Whole blood samples are one of the main sources used to obtain DNA, and there are many different protocols available to perform nucleic acid extraction on such samples. These methods vary from very basic manual protocols to more sophisticated methods included in automated DNA extraction protocols. Based on the wide range of available options, it would be ideal to determine the ones that perform best in terms of cost-effectiveness and time efficiency. We have reviewed DNA extraction history and the most commonly used methods for DNA extraction from whole blood samples, highlighting their individual advantages and disadvantages. We also searched current scientific literature to find studies comparing different nucleic acid extraction methods, to determine the best available choice. Based on our research, we have determined that there is not enough scientific evidence to support one particular DNA extraction method from whole blood samples. Choosing a suitable method is still a process that requires consideration of many different factors, and more research is needed to validate choices made at facilities around the world.

Sub-genomic RNA of defective interfering (D.I.) dengue viral particles is replicated in the same manner as full length genomes

The predicted secondary structure of sub-genomic RNA in dengue virus defective interfering (D.I.) particles from patients, or generated in vitro, resembled that of the 3′ and 5′ regions of wild type dengue virus (DENV) genomes. While these structures in the sub-genomic RNA were found to be essential for its replication, their nucleotide sequences were not, so long as any new sequences maintained wild type RNA secondary structure. These observations suggested that these sub-genomic fragments of RNA from dengue viruses were replicated in the same manner as the full length genomes of their wild type, “helper”, viruses and that they probably represent the smallest fragments of DENV RNA that can be replicated during a natural infection. While D.I. particles containing sub-genomic RNA are completely parasitic, the relationship between wild type and D.I. DENV may be symbiotic, with the D.I. particles enhancing the transmission of infectious DENV.

BRCA1 deficiency exacerbates estrogen-induced DNA damage and genomic instability

Germline mutations in BRCA1 predispose carriers to a high incidence of breast and ovarian cancers. BRCA1 functions to maintain genomic stability through critical roles in DNA repair, cell cycle arrest and transcriptional control. A major question has been why BRCA1 loss or mutation leads to tumors mainly in estrogen-regulated tissues, given that BRCA1 has essential functions in all cell types. Here we report that estrogen and estrogen metabolites can cause DNA double strand breaks (DSB) in estrogen receptor-α negative breast cells and that BRCA1 is required to repair these DSBs to prevent metabolite-induced genomic instability. We found that BRCA1 also regulates estrogen metabolism and metabolite-mediated DNA damage by repressing the transcription of estrogen-metabolising enzymes, such as CYP1A1, in breast cells. Lastly, we used a knock-in human cell model with a heterozygous BRCA1 pathogenic mutation to show how BRCA1 haploinsufficiency affects these processes. Our findings provide pivotal new insights into why BRCA1 mutation drives the formation of tumours in estrogen-regulated tissues, despite the general role of BRCA1 in DNA repair in all cell types.

Identification of a BRCA1-mRNA splicing complex required for efficient DNA repair and maintenance of genomic stability

Mutations within BRCA1 predispose carriers to a high risk of breast and ovarian cancers. BRCA1 functions to maintain genomic stability through the assembly of multiple protein complexes involved in DNA repair, cell-cycle arrest, and transcriptional regulation. Here, we report the identification of a DNA damage-induced BRCA1 protein complex containing BCLAF1 and other key components of the mRNA-splicing machinery. In response to DNA damage, this complex regulates pre-mRNA splicing of a number of genes involved in DNA damage signaling and repair, thereby promoting the stability of these transcripts/proteins. Further, we show that abrogation of this complex results in sensitivity to DNA damage, defective DNA repair, and genomic instability. Interestingly, mutations in a number of proteins found within this complex have been identified in numerous cancer types. These data suggest that regulation of splicing by the BRCA1-mRNA splicing complex plays an important role in the cellular response to DNA damage.

Characterization of novel microsatellite markers derived from Korean rose bitterling (Rhodeus uyekii) genomic library

Korean rose bitterling (Rhodeus uyekii) is a freshwater fish endemic to Korea. Natural populations of this species have experienced severe declines as a result of habitat fragmentation and water pollution. To conserve and restore R. uyekii, the genetic diversity of this species needs to be assessed at the population level. Eighteen novel polymorphic microsatellite loci for R. uyekii were developed using an enriched partial genomic library. Polymorphisms at these loci were studied in 150 individuals collected from three populations. The number of alleles at each locus ranged from 3 to 47 (mean = 17.1). Within the populations, the observed heterozygosity ranged from 0.032 to 1.000, expected heterozygosity from 0.082 to 0.967, and polymorphism information content from 0.078 to 0.950. Six loci showed significant deviation from Hardy-Weinberg equilibrium after Bonferroni’s correction, and no significant linkage disequilibrium was detected between most locus pairs, except in three cases. These highly informative microsatellite markers should be useful for genetic population structure analyses of R. uyekii.

Optimizing de novo transcriptome assembly and extending genomic resources for striped catfish (Pangasianodon hypophthalmus)

Striped catfish (Pangasianodon hypophthalmus) is a commercially important freshwater fish used in inland aquaculture in the Mekong Delta, Vietnam. The culture industry is facing a significant challenge however from saltwater intrusion into many low topographical coastal provinces across the Mekong Delta as a result of predicted climate change impacts. Developing genomic resources for this species can facilitate the production of improved culture lines that can withstand raised salinity conditions, and so we have applied high-throughput Ion Torrent sequencing of transcriptome libraries from six target osmoregulatory organs from striped catfish as a genomic resource for use in future selection strategies. We obtained 12,177,770 reads after trimming and processing with an average length of 97 bp. De novo assemblies were generated using CLC Genomic Workbench, Trinity and Velvet/Oases with the best overall contig performance resulting from the CLC assembly. De novo assembly using CLC yielded 66,451 contigs with an average length of 478 bp and N50 length of 506 bp. A total of 37,969 contigs (57%) possessed significant similarity with proteins in the non-redundant database. Comparative analyses revealed that a significant number of contigs matched sequences reported in other teleost fishes, ranging in similarity from 45.2% with Atlantic cod to 52% with zebrafish. In addition, 28,879 simple sequence repeats (SSRs) and 55,721 single nucleotide polymorphisms (SNPs) were detected in the striped catfish transcriptome. The sequence collection generated in the current study represents the most comprehensive genomic resource for P. hypophthalmus available to date. Our results illustrate the utility of next-generation sequencing as an efficient tool for constructing a large genomic database for marker development in non-model species.

Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes

Determining the genetic bases of adaptations and their roles in speciation is a prominent issue in evolutionary biology. Cichlid fish species flocks are a prime example of recent rapid radiations, often associated with adaptive phenotypic divergence from a common ancestor within a short period of time. In several radiations of freshwater fishes, divergence in ecomorphological traits - including body shape, colour, lips and jaws - is thought to underlie their ecological differentiation, specialization and, ultimately, speciation. The Midas cichlid species complex (Amphilophus spp.) of Nicaragua provides one of the few known examples of sympatric speciation where species have rapidly evolved different but parallel morphologies in young crater lakes. This study identified significant QTL for body shape using SNPs generated via ddRAD sequencing and geometric morphometric analyses of a cross between two ecologically and morphologically divergent, sympatric cichlid species endemic to crater Lake Apoyo: an elongated limnetic species (Amphilophus zaliosus) and a high-bodied benthic species (Amphilophus astorquii). A total of 453 genome-wide informative SNPs were identified in 240 F-2 hybrids. These markers were used to construct a genetic map in which 25 linkage groups were resolved. Seventy-two segregating SNPs were linked to 11 QTL. By annotating the two most highly supported QTL-linked genomic regions, genes that might contribute to divergence in body shape along the benthic-limnetic axis in Midas cichlid sympatric adaptive radiations were identified. These results suggest that few genomic regions of large effect contribute to early stage divergence in Midas cichlids.