The GP problem: quantifying gene-to-phenotype relationships

In this paper we refer to the gene-to-phenotype modeling challenge as the GP problem. Integrating information across levels of organization within a genotype-environment system is a major challenge in computational biology. However, resolving the GP problem is a fundamental requirement if we are to understand and predict phenotypes given knowledge of the genome and model dynamic properties of biological systems. Organisms are consequences of this integration, and it is a major property of biological systems that underlies the responses we observe. We discuss the E(NK) model as a framework for investigation of the GP problem and the prediction of system properties at different levels of organization. We apply this quantitative framework to an investigation of the processes involved in genetic improvement of plants for agriculture. In our analysis, N genes determine the genetic variation for a set of traits that are responsible for plant adaptation to E environment-types within a target population of environments. The N genes can interact in epistatic NK gene-networks through the way that they influence plant growth and development processes within a dynamic crop growth model. We use a sorghum crop growth model, available within the APSIM agricultural production systems simulation model, to integrate the gene-environment interactions that occur during growth and development and to predict genotype-to-phenotype relationships for a given E(NK) model. Directional selection is then applied to the population of genotypes, based on their predicted phenotypes, to simulate the dynamic aspects of genetic improvement by a plant-breeding program. The outcomes of the simulated breeding are evaluated across cycles of selection in terms of the changes in allele frequencies for the N genes and the genotypic and phenotypic values of the populations of genotypes.

Restriction-modification gene complexes as selfish gene entities: Roles of a regulatory system in their establishment, maintenance, and apoptotic mutual exclusion

We have reported some type II restriction-modification (RM) gene complexes on plasmids resist displacement by an incompatible plasmid through postsegregational host killing. Such selfish behavior may have contributed to the spread and maintenance of RM systems. Here we analyze the role of regulatory genes (C), often found linked to RM gene complexes, in their interaction with the host and the other RM gene complexes. We identified the C gene of EcoRV as a positive regulator of restriction. A C mutation eliminated postsegregational killing by EcoRV. The C system has been proposed to allow establishment of RM systems in new hosts by delaying the appearance of restriction activity. Consistent with this proposal, bacteria preexpressing ecoRVC were transformed at a reduced efficiency by plasmids carrying the EcoRV RM gene complex. Cells carrying the BamHI RM gene complex were transformed at a reduced efficiency by a plasmid carrying a PvuII RM gene complex, which shares the same C specificity. The reduction most likely was caused by chromosome cleavage at unmodified PvuII sites by prematurely expressed PvuII restriction enzyme. Therefore, association of the C genes of the same specificity with RM gene complexes of different sequence specificities can confer on a resident RM gene complex the capacity to abort establishment of a second, incoming RM gene complex. This phenomenon, termed “apoptotic mutual exclusion,” is reminiscent of suicidal defense against virus infection programmed by other selfish elements. pvuIIC and bamHIC genes define one incompatibility group of exclusion whereas ecoRVC gene defines another.

Identification of planarian homeobox sequences indicates the antiquity of most Hox/homeotic gene subclasses.

The homeotic gene complex (HOM-C) is a cluster of genes involved in the anteroposterior axial patterning of animal embryos. It is composed of homeobox genes belonging to the Hox/HOM superclass. Originally discovered in Drosophila, Hox/HOM genes have been identified in organisms as distantly related as arthropods, vertebrates, nematodes, and cnidarians. Data obtained in parallel from the organization of the complex, the domains of gene expression during embryogenesis, and phylogenetic relationships allow the subdivision of the Hox/HOM superclass into five classes (lab, pb/Hox3, Dfd, Antp, and Abd-B) that appeared early during metazoan evolution. We describe a search for homologues of these genes in platyhelminths, triploblast metazoans emerging as an outgroup to the great coelomate ensemble. A degenerate PCR screening for Hox/HOM homeoboxes in three species of triclad planarians has revealed 10 types of Antennapedia-like genes. The homeobox-containing sequences of these PCR fragments allowed the amplification of the homeobox-coding exons for five of these genes in the species Polycelis nigra. A phylogenetic analysis shows that two genes are clear orthologues of Drosophila labial, four others are members of a Dfd/Antp superclass, and a seventh gene, although more difficult to classify with certainty, may be related to the pb/Hox3 class. Together with previously identified Hox/HOM genes in other flatworms, our analyses demonstrate the existence of an elaborate family of Hox/HOM genes in the ancestor of all triploblast animals.

Non-major-histocompatibility-complex genetics of ankylosing spondylitis

There is strong evidence from twin and family studies indicating that a substantial proportion of the heritability of susceptibility to ankylosing spondylitis (AS) and its clinical manifestations is encoded by non-major-histocompatibility-complex genes. Efforts to identify these genes have included genomewide linkage studies and candidate gene association studies. One region, the interleukin (IL)-I gene complex on chromosome 2, has been repeatedly associated with AS in both Caucasians and Asians. It is likely that more than one gene in this complex is involved in AS, with the strongest evidence to date implicating IL-IA. Identifying the genes underlying other linkage regions has been difficult due to the lack of obvious candidates and the low power of most studies to date to identify genes of the small to moderate magnitude that are likely to be involved. The field is moving towards genomewide association analysis, involving much larger datasets of unrelated cases and controls. Early successes using this approach in other diseases indicates that it is likely to identify genes in common diseases like AS, but there remains the risk that the common-variant, common-disease hypothesis will not hold true in AS. Nonetheless, it is appropriate for the field to be cautiously optimistic that the next few years will bring great advances in our understanding of the genetics of this condition.

Staphylococcus hominis subsp novobiosepticus strains causing nosocomial bloodstream infection in Brazil

To report the isolation of six Staphylococcus hominis subsp. novobiosepticus (SHN) strains from hospitalized patients with bloodstream infections in two Brazilian hospitals and to characterize their susceptibility profile to several antimicrobials. Species identification was performed by biochemical methods and sodA gene sequencing. The MICs of antimicrobials were determined by broth and agar dilution methods and by Etest. Isolates were typed by PFGE and PCR amplification was used to detect the ccr gene complex and the mec class. Morphometric evaluation of cell wall was performed by transmission electron microscopy (TEM). Susceptibility profiles indicated that the majority of isolates (five) were multidrug-resistant. Overlapping and multiplex PCR showed that five out of the six strains harboured SCCmec type III with class A mec and type 3 ccr. The initial vancomycin MIC value of 4 mg/L for these strains increased to 16-32 mg/L after growth for 10 days in BHI broth supplemented with this antimicrobial. TEM indicated that vancomycin resistance was associated with cell wall thickening and to another mechanism not fully elucidated. Only one SHN strain was oxacillin- and vancomycin-susceptible. The nosocomial infections in at least five of the patients from both hospitals were caused by a single clone of SHN. It is very important to consider SHN strains as the cause of nosocomial infections. The clinical implications resulting from the pattern of multidrug resistance in these strains may be complicated by the emergence of vancomycin resistance.

Suppression of Myc oncogenic activity by ribosomal protein haploinsufficiency

The Myc oncogene regulates the expression of several components of the protein synthetic machinery, including ribosomal proteins, initiation factors of translation, RNA polymerase III and ribosomal DNA(1,2). Whether and how increasing the cellular protein synthesis capacity affects the multistep process leading to cancer remains to be addressed. Here we use ribosomal protein heterozygote mice as a genetic tool to restore increased protein synthesis in E mu-Myc/+ transgenic mice to normal levels, and show that the oncogenic potential of Myc in this context is suppressed. Our findings demonstrate that the ability of Myc to increase protein synthesis directly augments cell size and is sufficient to accelerate cell cycle progression independently of known cell cycle targets transcriptionally regulated by Myc. In addition, when protein synthesis is restored to normal levels, Myc- overexpressing precancerous cells are more efficiently eliminated by programmed cell death. Our findings reveal a new mechanism that links increases in general protein synthesis rates downstream of an oncogenic signal to a specific molecular impairment in the modality of translation initiation used to regulate the expression of selective messenger RNAs. We show that an aberrant increase in cap- dependent translation downstream of Myc hyperactivation specifically impairs the translational switch to internal ribosomal entry site ( IRES)- dependent translation that is required for accurate mitotic progression. Failure of this translational switch results in reduced mitotic- specific expression of the endogenous IRES- dependent form of Cdk11 ( also known as Cdc21 and PITSLRE)(3-5), which leads to cytokinesis defects and is associated with increased centrosome numbers and genome instability in E mu-Myc/+ mice. When accurate translational control is re- established in E mu-Myc/+ mice, genome instability is suppressed. Our findings demonstrate how perturbations in translational control provide a highly specific outcome for gene expression, genome stability and cancer initiation that have important implications for understanding the molecular mechanism of cancer formation at the post- genomic level.

Homeotic genes influence the axonal pathway of a Drosophila embryonic sensory neuron

Each abdominal hemisegment of the Drosophila embryo has two sensory neurons intimately associated with a tracheal branch. During embryogenesis, the axons of these sensory neurons, termed the v'td2 neurons, enter the CNS and grow toward the brain with a distinctive pathway change in the third thoracic neuromere. We show that the axons use guidance cues that are under control of the bithorax gene complex (BX-C). Pathway defects in mutants suggest that a drop in Ultrabithorax expression permits the pathway change in the T3 neuromere, while combined Ultrabithorax and abdominal-A expression represses it in the abdominal neuromeres. We propose that the axons do not respond to a particular segmental identity in forming the pathway change; rather they respond to pathfinding cues that come about as a result of a drop in BX-C expression along the antero-posterior axis of the CNS.

Unilateral giant renal angiomyolipoma and pulmonary lymphangioleiomyomatosis

Angiomyolipomas (AMLs) are mesenchymal neoplasms, named so because of the complex tissue composition represented by variable proportions of mature adipose tissue, smooth muscle cells, and dysmorphic blood vessels. Although AMLs may rise in different sites of the body, they are mostly observed in the kidney and liver. In the case of renal AMLs, they are described in two types: isolated AMLs and AMLs associated with tuberous sclerosis (TS). While most cases of AMLs are found incidentally during imaging examinations and are asymptomatic, others may reach huge proportions causing symptoms. Pulmonary lymphangioleiomyomatosis (LAM) is a rare benign disease characterized by cystic changes in the pulmonary parenchyma and smooth muscle proliferation, leading to a mixed picture of interstitial and obstructive disease. AML and LAM constitute major features of tuberous sclerosis complex (TSC), a multisystem autosomal dominant tumor-suppressor gene complex diagnosis. The authors report the case of a young female patient who presented a huge abdominal tumor, which at computed tomography (CT) show a fat predominance. The tumor displaced the right kidney and remaining abdominal viscera to the left. Chest CT also disclosed pulmonary lesions compatible with lymphangioleiomyomatosis. Because of sudden abdominal pain accompanied by a fall in the hemoglobin level, the patient underwent an urgent laparotomy. The excised tumor was shown to be a giant renal AML with signs of bleeding in its interior. The authors call attention to the diagnosis of AML and the huge proportions that the tumor can reach, as well as for ruling out the TSC diagnosis, once it may impose genetic counseling implications.

Novel pseudo-staphylococcal cassette chromosome mec element (ψSCCmec57395) in methicillin-resistant Staphylococcus pseudintermedius CC45

Genetic characterization of methicillin-resistant Staphylococcus pseudintermedius (MRSP) from Thailand and Israel revealed the presence of a predominant atypical clonal lineage which was not typeable by SmaI-PFGE and SCCmec typing. All the atypical isolates (n = 34) belonged to CC45 (30 ST45 and 2 ST179 isolates, 1 ST57 isolate, and 1 ST85 isolate). The isolates originated from healthy and diseased dogs and cats, as well as from the environment of one clinic. Cfr9I-pulsed-field gel electrophoresis (Cfr9I-PFGE) and dru typing permitted the further distinction of CC45 isolates from the two different countries. Microarray analysis identified genes that confer resistance to β-lactams (mecA; blaZ), aminoglycosides [aac(6')-Ie-aph(2')-Ia; aph(3')-III; ant(6)-Ia], macrolides and lincosamides [erm(B)], tetracyclines [tet(M)], trimethoprim [dfr(G)], streptothricin (sat4), and chloramphenicol (catpC221). Fluoroquinolone resistance was attributed to specific amino acid substitutions, i.e., Ser84Leu in GyrA and Ser80Ile and Asp84Asn in GrlA. A novel pseudo-staphylococcal cassette chromosome (ΨSCCmec57395) element was identified in MRSP strain 57395 (sequence type ST45) by whole-genome sequencing. The 12,282-bp ΨSCCmec57395 element contained a class C1 mec gene complex but no ccr genes. In addition to the methicillin resistance gene mecA, ΨSCCmec57395 also carried determinants of resistance to heavy metals, such as arsenic, cadmium, and copper. Bsu36I restriction analysis of the ΨSCCmec57395 element amplified by long-range PCR revealed the presence of ΨSCCmec57395 in the 33 additional isolates of MRSP CC45. The ΨSCCmec57395 element represents a new class of SCCmec and has been identified in MRSP of CC45, which is a predominant clonal lineage in Israel and Thailand.

The mating-type and pathogenicity locus of the fungus Ustilago hordei spans a 500-kb region

The fungal pathogen Ustilago hordei causes the covered smut disease of barley and oats. Mating and pathogenicity in this fungus are controlled by the MAT locus, which contains two distinct gene complexes, a and b. In this study, we tagged the a and b regions with the recognition sequence for the restriction enzyme I-SceI and determined that the distance between the complexes is 500 kb in a MAT-1 strain and 430 kb in a MAT-2 strain. Characterization of the organization of the known genes within the a and b gene complexes provided evidence for nonhomology and sequence inversion between MAT-1 and MAT-2. Antibiotic-resistance markers also were used to tag the a gene complex in MAT-1 strains (phleomycin) and the b gene complex in MAT-2 strains (hygromycin). Crosses were performed with these strains and progeny resistant to both antibiotics were recovered at a very low frequency, suggesting that recombination is suppressed within the MAT region. Overall, the chromosome homologues carrying the MAT locus of U. hordei share features with primitive sex chromosomes, with the added twist that the MAT locus also controls pathogenicity.

Immune-type receptor genes in zebrafish share genetic and functional properties with genes encoded by the mammalian leukocyte receptor cluster

An extensive, highly diversified multigene family of novel immune-type receptor (nitr) genes has been defined in Danio rerio (zebrafish). The genes are predicted to encode type I transmembrane glycoproteins consisting of extracellular variable (V) and V-like C2 (V/C2) domains, a transmembrane region and a cytoplasmic tail. All of the genes examined encode immunoreceptor tyrosine-based inhibition motifs in the cytoplasmic tail. Radiation hybrid panel mapping and analysis of a deletion mutant line (b240) indicate that a minimum of ≈40 nitr genes are contiguous in the genome and span ≈0.6 Mb near the top of zebrafish linkage group 7. One flanking region of the nitr gene complex shares conserved synteny with a region of mouse chromosome 7, which shares conserved synteny with human 19q13.3-q13.4 that encodes the leukocyte receptor cluster. Antibody-induced crosslinking of Nitrs that have been introduced into a human natural killer cell line inhibits the phosphorylation of mitogen-activated protein kinase that is triggered by natural killer-sensitive tumor target cells. Nitrs likely represent intermediates in the evolution of the leukocyte receptor cluster.

The estrogen receptor 1 G594A polymorphism is associated with migraine susceptibility in two independent case/control groups

Migraine is a painful and debilitating disorder with a significant genetic component. Steroid hormones, in particular estrogen, have long been considered to play a role in migraine, as variations in hormone levels are associated with migraine onset in many sufferers of the disorder. Steroid hormones mediate their activity via hormone receptors, which have a wide tissue distribution. Estrogen receptors have been localized to the brain in regions considered to be involved in migraine pathogenesis. Hence it is possible that genetic variation in the estrogen receptor gene may play a role in migraine susceptibility. This study thus examined the estrogen receptor 1 (ESRalpha) gene for a potential role in migraine pathogenesis and susceptibility. A population-based cohort of 224 migraine sufferers and 224 matched controls were genotyped for the G594A polymorphism located in exon 8 of the ESR1 gene. Statistical analysis indicated a significant difference between migraineurs and non-migraineurs in both the allele frequencies (P=0.003) and genotype distributions (P=0.008) in this sample. An independent follow-up study was then undertaken using this marker in an additional population-based cohort of 260 migraine sufferers and 260 matched controls. This resulted in a significant association between the two groups with regard to allele frequencies (P=8x10(-6)) and genotype distributions (P=4x10(-5)). Our findings support the hypothesis that genetic variation in hormone receptors, in particular the ESR1 gene, may play a role in migraine.

The effect of pH on PAMAM dendrimer-siRNA complexation:endosomal considerations as determined by molecular dynamics simulation

Intracellular degradation of genes, most notably within the endo-lysosomal compartment is considered a significant barrier to (non-viral) gene delivery in vivo. Previous reports based on in vitro studies claim that carriers possessing a mixture of primary, secondary and tertiary amines are able to buffer the acidic environment within the endosome, allowing for timely release of their contents, leading to higher transfection rates. In this report, we adopt an atomistic molecular dynamics (MD) simulation approach, comparing the complexation of 21-bp siRNA with low-generation polyamidoamine (PAMAM) dendrimers (G0 and G1) at both neutral and acidic pHs, the latter of which mimics the degradative environment within maturing 'late-endosomes'. Our simulations reveal that the time taken for the dendrimer-gene complex (dendriplex) to reach equilibrium is appreciably longer at low pH and this is accompanied by more compact packaging of the dendriplex, as compared to simulations performed at neutral pH. We also note larger absolute values of calculated binding free energies of the dendriplex at low pH, indicating a higher dendrimer-nucleic acid affinity in comparison with neutral pH. These novel simulations provide a more detailed understanding of low molecular-weight polymer-siRNA behavior, mimicking the endosomal environment and provide input of direct relevance to the "proton sponge theory", thereby advancing the rational design of non-viral gene delivery systems.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Differential tumor surveillance by natural killer (NK) and NKT cells

Natural tumor surveillance capabilities of the host were investigated in six different mouse tumor models where endogenous interleukin (IL)-12. does or does not dictate the efficiency of the innate immune response. Gene-targeted and lymphocyte subset-depleted mice were used to establish the relative importance of natural killer (NK) and NK1.1(+) T (NKT) cells in protection from tumor initiation and metastasis. In the models examined, CD3(-) NK cells were responsible for tumor rejection and protection from metastasis in models where control of major histocompatibility complex class I-deficient tumors was independent of IL-12, A protective role for NKT cells was only observed when tumor rejection required endogenous IL-12 activity. In particular, T cell receptor J alpha 281 gene-targeted mice confirmed a critical function for NKT cells in protection from spontaneous tumors initiated by the chemical carcinogen, methylcholanthrene. This is the first description of an antitumor function for NKT cells in the absence of exogenously administered potent stimulators such as IL-12 or alpha-galactosylceramide.