The Nature of Genetic Engineering and the Uses and Potential Abuses of Biotechnology

Over the past decade the topic of genetic engineering has been has been readily debated in the media, but often these debates consist of political rhetoric and fail to offer objective information on the methods and the potential benefits to human health and their environment. In truth, humans have been manipulating the genomes of organisms for thousands of years, and it has been an evolution of scientific knowledge that has led to the more precise methods of genetic engineering. This paper discusses how scientists utilize natural processes to alter the genetic constituents of both prokaryotic and eukaryotic organisms, benefits to human health and the environment, as well as potential misuses of biotechnology such as bioterrorism.

Whole genome analysis of Fusobacterium nucleatum, Rickettsia typhi and Francisella tularensis

The genomes of Fusobacterium nucleatum subspecies polymorphum strain ATCC 10953, Rickettsia typhi strain Wilmington, and Francisella tularensis subspecies holarctica strain OSU18 were sequenced, annotated, and analyzed. Each genome was then compared to the sequenced genomes of closely related bacteria. The genome of F. nucleatum ATCC 10953 was compared to two additional F. nucleatum subspecies, subspecies nucleatum and subspecies vincentii. This analysis revealed substantial evidence of horizontal gene transfer along with considerable genetic diversity within the species of F. nucleatum. R. typhi was compared to R. prowazekii and R. conorii. This analysis uncovered a hotspot for chromosomal rearrangements in the Spotted Fever Group but not the Typhus Group Rickettsia and revealed the close genetic relationship between the Typhus Group rickettsial species. F. tularensis OSU18 was compared to two additional F. tularensis strains. These comparisons uncovered significant chromosomal rearrangements between F. tularensis subspecies due to recombination between insertion sequence elements. ^

Roles of mismatch repair proteins in the recognition and processing of DNA interstrand crosslinks in human cells

DNA interstrand crosslinks (ICLs) are among the most toxic type of damage to a cell. Many ICL-inducing agents are widely used as therapeutic agents, e.g. cisplatin, psoralen. A bettor understanding of the cellular mechanism that eliminates ICLs is important for the improvement of human health. However, ICL repair is still poorly understood in mammals. Using a triplex-directed site-specific ICL model, we studied the roles of mismatch repair (MMR) proteins in ICL repair in human cells. We are also interested in using psoralen-conjugated triplex-forming oligonucleotides (TFOs) to direct ICLs to a specific site in targeted DNA and in the mammalian genomes. ^ MSH2 protein is the common subunit of two MMR recognition complexes, and MutSα and MutSβ. We showed that MSH2 deficiency renders human cell hypersensitive to psoralen ICLs. MMR recognition complexes bind specifically to triplex-directed psoralen ICLs in vitro. Together with the fact that psoralen ICL-induced repair synthesis is dramatically decreased in MSH2 deficient cell extracts, we demonstrated that MSH2 function is critical for the recognition and processing of psoralen ICLs in human cells. Interestingly, lack of MSH2 does not reduce the level of psoralen ICL-induced mutagenesis in human cells, suggesting that MSH2 does not contribute to error-generating repair of psoralen ICLs, and therefore, may represent a novel error-free mechanism for repairing ICLs. We also studied the role of MLH1, anther key protein in MMR, in the processing of psoralen ICLs. MLH1-deficient human cells are more resistant to psoralen plus UVA treatment. Importantly, MLH1 function is not required for the mutagenic repair of psoralen ICLs, suggesting that it is not involved in the error-generating repair of this type of DNA damage in human cells. ^ These are the first data indicating mismatch repair proteins may participate in a relatively error-free mechanism for processing psoralen ICL in human cells. Enhancement of MMR protein function relative to nucleotide excision repair proteins may reduce the mutagenesis caused by DNA ICLs in humans. ^ In order to specifically target ICLs to mammalian genes, we identified novel TFO target sequences in mouse and human genomes. Using this information, many critical mammalian genes can now be targeted by TFOs.^

The molecular basis for beta-lactamase gene expression in B. anthracis, B. cereus and B. thuringiensis

The susceptibility of most Bacillus anthracis strains to β-lactam antibiotics is intriguing considering that the B. anthracis genome harbors two β-lactamase genes, bla1 and bla2, and closely-related species, Bacillus cereus and Bacillus thuringiensis, typically produce β-lactamases. This work demonstrates that B. anthracis bla expression is affected by two genes, sigP and rsp, predicted to encode an extracytoplasmic function sigma factor and an antisigma factor, respectively. Deletion of the sigP/rsp locus abolished bla expression in a penicillin-resistant clinical isolate and had no effect on bla expression in a prototypical penicillin-susceptible strain. Complementation with sigP/rsp from the penicillin-resistant strain, but not the penicillin-susceptible strain, conferred β-lactamase activity upon both mutants. These results are attributed to a nucleotide deletion near the 5' end of rsp in the penicillin-resistant strain that is predicted to result in a nonfunctional protein. B. cereus and B. thuringiensis sigP and rsp homologues are required for inducible penicillin resistance in those species. Expression of the B. cereus or B. thuringiensis sigP and rsp genes in a B. anthracis sigP/rsp-null mutant confers resistance to β-lactam antibiotics, suggesting that while B. anthracis contains the genes necessary for sensing β-lactam antibiotics, the B. anthracis sigP/rsp gene products are insufficient for bla induction. ^ Because alternative sigma factors recognize unique promoter sequence, direct targets can be elucidated by comparing transcriptional profiling results with an in silico search using the sigma factor binding sequence. Potential σP -10 and -35 promoter elements were identified upstream from bla1 bla2 and sigP. Results obtained from searching the B. anthracis genome with the conserved sequences were evaluated against transcriptional profiling results comparing B. anthracis 32 and an isogenic sigP/rsp -null strain. Results from these analyses indicate that while the absence of the sigP gene significantly affects the transcript levels of 16 genes, only bla1, bla2 and sigP are directly regulated by σP. The genomes of B. cereus and B. thuringiensis strains were also analyzed for the potential σP binding elements. The sequence was located upstream from the sigP and bla genes, and previously unidentified genes predicted to encode a penicillin-binding protein (PBP) and a D-alanyl-D-alanine carboxypeptidase, indicating that the σ P regulon in these species responds to cell-wall stress caused by β-lactam antibiotics. ^ β-lactam antibiotics prevent attachment of new peptidoglycan to the cell wall by blocking the active site of PBPs. A B. cereus and B. thuringiensis pbp-encoding gene located near bla1 contains a potential σP recognition sequence upstream from the annotated translational start. Deletion of this gene abolished β-lactam resistance in both strains. Mutations in the active site of the PBP were detrimental to β-lactam resistance in B. cereus, but not B. thuringiensis, indicating that the transpeptidase activity is only important in B. cereus. I also found that transcript levels of the PBP-encoding gene are not significantly affected by the presence of β-lactam antibiotic. Based on these data I hypothesize that the gene product acts a sensor of β-lactam antibiotic. ^

Exploration of genetic susceptibility to pancreatic cancer: A GWAS approach

Pancreatic cancer is the 4th most common cause for cancer death in the United States, accompanied by less than 5% five-year survival rate based on current treatments, particularly because it is usually detected at a late stage. Identifying a high-risk population to launch an effective preventive strategy and intervention to control this highly lethal disease is desperately needed. The genetic etiology of pancreatic cancer has not been well profiled. We hypothesized that unidentified genetic variants by previous genome-wide association study (GWAS) for pancreatic cancer, due to stringent statistical threshold or missing interaction analysis, may be unveiled using alternative approaches. To achieve this aim, we explored genetic susceptibility to pancreatic cancer in terms of marginal associations of pathway and genes, as well as their interactions with risk factors. We conducted pathway- and gene-based analysis using GWAS data from 3141 pancreatic cancer patients and 3367 controls with European ancestry. Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Using the logistic kernel machine (LKM) test, we analyzed 17906 genes defined by University of California Santa Cruz (UCSC) database. Using the likelihood ratio test (LRT) in a logistic regression model, we analyzed 177 pathways and 17906 genes for interactions with risk factors in 2028 pancreatic cancer patients and 2109 controls with European ancestry. After adjusting for multiple comparisons, six pathways were marginally associated with risk of pancreatic cancer ( P < 0.00025): Fc epsilon RI signaling, maturity onset diabetes of the young, neuroactive ligand-receptor interaction, long-term depression (Ps < 0.0002), and the olfactory transduction and vascular smooth muscle contraction pathways (P = 0.0002; Nine genes were marginally associated with pancreatic cancer risk (P < 2.62 × 10−5), including five reported genes (ABO, HNF1A, CLPTM1L, SHH and MYC), as well as four novel genes (OR13C4, OR 13C3, KCNA6 and HNF4 G); three pathways significantly interacted with risk factors on modifying the risk of pancreatic cancer (P < 2.82 × 10−4): chemokine signaling pathway with obesity ( P < 1.43 × 10−4), calcium signaling pathway (P < 2.27 × 10−4) and MAPK signaling pathway with diabetes (P < 2.77 × 10−4). However, none of the 17906 genes tested for interactions survived the multiple comparisons corrections. In summary, our current GWAS study unveiled unidentified genetic susceptibility to pancreatic cancer using alternative methods. These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer, once confirmed, will shed promising light on the prevention and treatment of this disease. ^

MOLECULAR AND GENOMIC BASED INSIGHTS INTO THE EVOLUTION OF ENTEROCOCCUS FAECIUM FROM COMMENSAL TO HOSPITAL-ADAPTED PATHOGEN

The basis for the recent transition of Enterococcus faecium from a primarily commensal organism to one of the leading causes of hospital-acquired infections in the United States is not yet understood. To address this, the first part of my project assessed isolates from early outbreaks in the USA and South America using sequence analysis, colony hybridizations, and minimal inhibitory concentrations (MICs) which showed clinical isolates possess virulence and antibiotic resistance determinants that are less abundant or lacking in community isolates. I also revealed that the level of ampicillin resistance increased over time in clinical strains. By sequencing the pbp5 gene, I demonstrated an ~5% difference in the pbp5 gene between strains with MICs <4ug/ml and those with MICs >4µg/ml, but no specific sequence changes correlated with increases in MICs within the latter group. A 3-10% nucleotide difference was also seen in three other genes analyzed, which suggested the existence of two distinct subpopulations of E. faecium. This led to the second part of my project analyzing concatenated core gene sequences, SNPs, the 16S rRNA, and phylogenetics of 21 E. faecium genomes confirming two distinct clades; a community-associated (CA) clade and hospital-associated (HA) clade. Molecular clock calculations indicate that these two clades likely diverged ~ 300,000 to > 1 million years ago, long before the modern antibiotic era. Genomic analysis also showed that, in addition to core genomic differences, HA E. faecium harbor specific accessory genetic elements that may confer selection advantages over CA E. faecium. The third part of my project discovered 6 E. faecium genes with the newly identified “WxL” domain. My analyses, using RT-PCR, western blots, patient sera, whole-cell ELISA, and immunogold electron microscopy, indicated that E. faecium WxL genes exist in operons, encode bacterial cell surface localized proteins, that WxL proteins are antigenic in humans, and are more exposed on the surface of clinical isolates versus community isolates (even though they are ubiquitous in both clades). ELISAs and BIAcore analyses also showed that proteins encoded by these operons bind several different host extracellular matrix proteins, as well as to each other, suggesting a novel cell-surface complex. In summary, my studies provide new insights into the evolution of E. faecium by showing that there are two distantly related clades; one being more successful in the hospital setting. My studies also identified operons encoding WxL proteins whose characteristics could also contribute to colonization and virulence within this species.

Molecular analysis of a neurovirulent murine leukemia virus, {\it ts\/}1

ts1 is a neurovirulent spontaneous temperature-sensitive mutant of Moloney murine leukemia virus TB (MoMuLV-TB). MoMuLV-TB causes T-cell lymphoma or lymphoid leukemia in mice after a long latency period whereas ts1 causes a progressive hindlimb paralytic disease after a much shorter latency period. In previous studies, it had been shown that the temperature-sensitive defect resided in the $env$ gene. At the restrictive temperature, the envelope precursor polyprotein, gPr80$\sp{env}$, is inefficiently processed intracellularly into a heterodimer consisting of two cleavage products, gp70 and Prp15E. This inefficient processing is correlated with neurovirulence. In this study, the nucleotide sequences of the env genes for both ts1 and MoMuLV-TB were determined, and the encoded amino acid sequences were deduced from the DNA sequences. There were four unique amino acid substitutions in the gPr80$\sp{env}$ of ts1. In order to determine which unique amino acid was responsible for the phenotypic characteristics of ts1, a set of hybrid genomes was constructed by exchanging restriction fragments between ts1 and MoMuLV-TB. NIH 3T3 cells were transfected with the hybrid genomes to obtain infectious hybrid viruses. Assays of the hybrid viruses showed that a Val-25$\to$Ile substitution in gPr80$\sp{env}$ was responsible for the temperature sensitivity, inefficient processing, and neurovirulence of ts1. In further studies, the Ile-25 in gPr80$\sp{env}$ was substituted with Thr, Ala, Leu, Gly, and Glu by site-directed mutagenesis to generate a new set of mutant viruses, i.e., ts1-T, -A, -L, -G, and -E, respectively. The rank order of the mutants for temperature sensitivity was: ts1-E $>$ ts1-G $>$ ts1-L $>$ ts1-A $>$ ts1 $>$ ts1-T. The degree of temperature sensitivity of each of the mutants also correlated with the degree of inefficient processing of gPr80$\sp{env}$. The mutant viruses were assayed for neurovirulence. ts1-T caused whole body tremor, ts1-A caused hindlimb paralysis, ts1-L caused paraparesis, but ts1-G and -E were not neurovirulent. These results show that inefficient processing of gPr80$\sp{env}$ is correlated with neurovirulence, but if processing of gPr80$\sp{env}$ is too inefficient there is no neurovirulence. Furthermore, the disease profile of each of the neurovirulent viruses depends on the degree of inefficient processing of gPr80$\sp{env}$. ^

Genome -wide and locus -specific instability

The discovery of expanded simple repeated sequences causing or associated with human disease has lead to a new area of research involved in the elucidation of how the expanded repeat causes disease and how the repeat becomes unstable. ^ To study the genetic basis of the (CTG)n repeat instability in the DMPK gene in myotonic dystrophy (DM1) patients, somatic cell hybrids were constructed between the lymphocytes of DM1 patients and a variety of Chinese hamster ovary (CHO) cell DNA repair gene deficient mutants. By using small pool PCR (SP-PCR), the instability of the (CTG)n can be quantitated for both the frequency and sizes of length change mutations. ^ Additional SP-PCR analysis on 2/11 subclones generated from this original hybrid showed a marked increase in large repeat deletions, ∼50%. A bimodal distribution of repeats was seen around the progenitor allele and at a large deleted product (within the normal range) with no intermediate products present. ^ To determine if the repair capacity of the CHO cell led to a mutator phenotype in the hamster and hybrid clones, SP-PCR was also done on 3 hamster microsatellites in a variety of hamster cell backgrounds. No variant alleles were seen in over 2500 genome equivalents screened. ^ Human-hamster hybrids have long been shown to be chromosomally unstable, yet information about the stability of repeated sequences was not known. To test if repeat instability was associated with either intact or non-intact human chromosomes, more than 300 microsatellite repeats on 13 human chromosomes (intact and non-intact) were analyzed in eight hybrid cells. No variants were seen between the hybrid and patient alleles in the hybrids. ^ To identify whether DM1 patients have a previously undetected level of genome wide instability or if the instability is truly locus specific, SP-PCR was done on 6 human microsatellites within the patient used to make the hybrid cells. No variants were seen in over 1000 genomes screened. ^ These studies show that the somatic cell hybrid approach is a genetically stable system that allows for the determination of factors that could lead to changes in microsatellite instability. It also shows that there is something inherent about the DM1 expanded (CTG)n repeat that it is solely targeted by, as of yet, and unknown mechanism that causes the repeat to be unstable. (Abstract shortened by UMI.)^

Cross -packing among retroviruses and characterization of the feline immunodeficiency virus (FIV) packaging signal: Implications for the development of FIV-based gene transfer systems

Feline immunodeficiency virus (FIV)-based gene transfer systems are being seriously considered for human gene therapy as an alternative to vectors based on primate lentiviruses, a genetically complex group of retroviruses capable of infecting non-dividing cells. The greater phylogenetic distance between the feline and primate lentiviruses is thought to reduce chances of the generation of recombinant viruses. However, safety of FIV-based vector systems has not been tested experimentally. Since primate lentiviruses such as human and simian immunodeficiency viruses (HIV/SIV) can cross-package each other's genomes, we tested this trait with respect to FIV. Unexpectedly, both feline and primate lentiviruses were reciprocally able to both cross-package and propagate each other's RNA genomes. This was largely due to the recognition of viral packaging signals by the heterologous proteins. However, a simple retrovirus such as Mason-Pfizer monkey virus (MPMV) was unable to package FIV RNA. Interestingly, FIV could package MPMV RNA, but not propagate it for further steps of replication. These findings suggest that upon co-infection of the same host, cross-packaging may allow distinct retroviruses to generate chimeric variants with unknown pathogenic potential. ^ In order to understand the packaging determinants in FIV, we conducted a detailed mutational analysis of the region thought to contain FIV packaging signal. We show that the first 90–120 nt of the 5′ untranslated region (UTR) and the first 90 nt of gag were simultaneously required for efficient FIV RNA packaging. These results suggest that the primary FIV packaging signal is multipartite and discontinuous, composed of two core elements separated by 150 nt of the 5 ′UTR. ^ The above studies are being used towards the development of safer FIV-based self-inactivating (SIN) vectors. These vectors are being designed to eliminate the ability of FIV transfer vector RNAs to be mobilized by primate lentiviral proteins that may be present in the target cells. Preliminary test of the first generation of these vectors has revealed that they are incapable of being propagated by feline proteins. The inability of FIV transfer vectors to express packageable vector RNA after integration should greatly increase the safety of FIV vectors for human gene therapy. ^

Evolution of cis-regulatory elements in a repetitive sequence adjacent to a sea urchin aboral ectoderm -specific gene

The creation, preservation, and degeneration of cis-regulatory elements controlling developmental gene expression are fundamental genome-level evolutionary processes about which little is known. In this study, critical differences in cis-regulatory elements controlling the expression of the sea urchin aboral ectoderm-specific spec genes were identified and explored. In genomes of species within the Strongylocentrotidae family, multiple copies of a repetitive sequence element termed RSR were present, but RSRs were not detected in genomes of species outside Strongylocentrotidae. RSRs are invariably associated with spec genes, and in Strongylocentrotus purpuratus, the spec2a RSR functioned as a transcriptional enhancer displaying greater activity than RSRs from the spec1 or spec2c paralogs. Single base-pair differences at two cis-regulatory elements within the spec2a RSR greatly increased the binding affinities of four transcription factors: SpCCAAT-binding factor at one element and SpOtx, SpGoosecoid, and SpGATA-E at another. The cis-regulatory elements to which SpCCAAT-binding factor, SpOtx, SpGoosecoid, and SpGATA-E bound were recent evolutionary acquisitions that could act either to activate or repress transcription, depending on the cell type. These elements were found in the spec2a RSR ortholog in Strongylocentrotus pallidus but not in the RSR orthologs of Strongylocentrotus droebachiensis or Hemicentrotus pulcherrimus. These results indicate that spec genes exhibit a dynamic pattern of cis-regulatory element evolution while stabilizing selection preserves their aboral ectoderm expression domain. ^

Populus spp : estabilidad y ganancia genética sobre la altura media dominante en tres ambientes de la pampa ondulada, Argentina

El objetivo de este trabajo fue evaluar la altura media dominante (AMD) y su estabilidad en 15 clones de álamos de 5 años de edad, implantados en tres ambientes diferentes de la pampa ondulada, Argentina. Los sitios fueron: Teodelina (Sitios 1 y 2), Santa Fe (34° 12' LS; 61° 43' W; 90 msnm) y Alberti (Sitio 3), Buenos Aires (34° 50' LS; 60° 30' W; 55 msnm) y se caracterizaron en base a variables de clima y suelo. Se realizaron los análisis de la varianza del conjunto de clones entre sitios y entre los clones por sitio. La comparación de AMD se realizó aplicando el test de Tukey. Se analizó la interacción clon-sitio. Se construyó un criterio de elección basado en parámetros genéticos de estabilidad sobre el carácter (AMD), estimándose la ganancia genética al seleccionar los mejores clones. Se realizaron los cálculos de heredabilidad en sentido amplio. Los posicionamientos en AMD fueron significativos entre clones por sitio y entre sitios. La interacción clon sitio fue significativa. Los valores de AMD y los de ecovalencia permitieron seleccionar genomas con mayor amplitud de adaptación y consecuentemente conciliar la producción con la plasticidad dentro de la disimilitud de los sitios evaluados.

Caracterización genética de germoplasma de maní cultivado (Arachis hypogaea L.) mediante el empleo de marcadores microsatélites

El maní cultivado (Arachis hypogaea L.), es una especie de gran importancia económica, nativo de América del Sur. Se divide en dos subespecies y seis variedades botánicas. Genéticamente es alotetraploide, constituido por dos juegos genómicos duplicados. El empleo de marcadores microsatélites resulta más apropiado para realizar la caracterización genética de esta especie, puesto que permiten detectar un elevado nivel de polimorfismo. El objetivo del presente trabajo fue caracterizar la diversidad genética existente en las entradas de germoplasma de maní cultivado pertenecientes al Banco Activo de Germoplasma de Maní del Instituto Nacional de Tecnología Agropecuaria (INTA). Veinticinco entradas fueron genotipificadas con 23 marcadores microsatélites, de los cuales, 17 resultaron polimórficos. Se observaron 75 fragmentos polimórficos amplificados, con un promedio de 4,41 alelos por locus y un rango de 1 a 9 alelos. El contenido de información polimórfica osciló entre 0,15 y 0,58. El valor de la diversidad genética promedio fue de 0,165. Tanto el análisis de conglomerados como el de coordenadas principales evidenciaron dos grupos, uno formado por los materiales representantes de la subespecie fastigiata y otro por los de la subespecie hypogaea. Los resultados del análisis molecular de la varianza mostraron varianza tanto dentro como entre, las subespecies analizadas.

Experiment on marine fungus Microascus brevicaulis strain

The marine fungus Microascus brevicaulis strain LF580 is a non-model secondary metabolite producer with high yields of the two secondary metabolites scopularide A and B, which exhibit distinct activities against tumour cell lines. A mutant strain was obtained using UV mutagenesis, showing besides higher production levels faster growth and differences in pellet formation. Comparative proteomics were applied to gain deeper understanding of the regulation of production and of the physiology of this fungus and its mutant. For this purpose, an optimised protein extraction protocol was established. Here, we show the first proteome study of a marine fungus. In total, 4759 proteins were identified. The central metabolic pathway of LF580 could be mapped by using KEGG pathway analysis and GO annotation. Using iTRAQ labelling, 318 proteins were shown to be significantly regulated in the mutant strain: 189 were down- and 129 upregulated. Proteomics are a powerful tool for the understanding of regulatory aspects: The differences on proteome level could be attributed to a limited nutrient availability in wild type strain due to a strong pellet formation. This information can be applied to optimisation on strain and process level. The linkage between nutrient limitation and pellet formation in the non-model fungus M. brevicaulis is in consensus with the knowledge on model organisms like Aspergillus niger and Penicillium chrysogenum.

Individual codes, GenBank accession numbers and GenBank accession links for Mediterranean sponges

Metabolomic analysis has shown the chemical richness of the sponge-associated actinomycetes Streptomyces sp. SBT349, Nonomureae sp. SBT364, and Nocardiopsis sp. SBT366. The genomes of these actinomycetes were sequenced and the genomic potential for secondary metabolism was evaluated. Their draft genomes have sizes of 8.0, 10, and 5.8Mb having 687, 367, and 179 contigs with a GC content of 71.6, 70.7, and 72.7%, respectively. Moreover, antiSMASH 3.0 predicted 108, 149, and 75 secondary metabolite gene clusters, respectively which highlight the metabolic capacity of the three actinomycete species to produce diverse classes of natural products.

Individual codes, GenBank accession numbers and GenBank accession links for sponge-associated actinobacterial isolates

Here, we report the draft genome sequences of three actinobacterial isolates, Micromonospora sp. RV43, Rubrobacter sp. RV113, and Nocardiopsis sp. RV163 that had previously been isolated from Mediterranean sponges. The draft genomes were analyzed for the presence of gene clusters indicative of secondary metabolism using antiSMASH 3.0 and NapDos pipelines. Our findings demonstrated the chemical richness of sponge-associated actinomycetes and the efficacy of genome mining in exploring the genomic potential of sponge-derived actinomycetes.