Gene Conversion and Functional Divergence in the β-Globin Gene Family

Different models of gene family evolution have been proposed to explain the mechanism whereby gene copies created by gene duplications are maintained and diverge in function. Ohta proposed a model which predicts a burst of nonsynonymous substitutions following gene duplication and the preservation of duplicates through positive selection. An alternative model, the duplication–degeneration–complementation (DDC) model, does not explicitly require the action of positive Darwinian selection for the maintenance of duplicated gene copies, although purifying selection is assumed to continue to act on both copies. A potential outcome of the DDC model is heterogeneity in purifying selection among the gene copies, due to partitioning of subfunctions which complement each other. By using the d_N/d_S () rate ratio to measure selection pressure, we can distinguish between these two very different evolutionary scenarios. In this study we investigated these scenarios in the -globin family of genes, a textbook example of evolution by gene duplication. We assembled a comprehensive dataset of 72 vertebrate -globin sequences. The estimated phylogeny suggested multiple gene duplication and gene conversion events. By using different programs to detect recombination, we confirmed several cases of gene conversion and detected two new cases. We tested evolutionary scenarios derived from Ohtas model and the DDC model by examining selective pressures along lineages in a phylogeny of -globin genes in eutherian mammals. We did not find significant evidence for an increase in the ratio following major duplication events in this family. However, one exception to this pattern was the duplication of -globin in simian primates, after which a few sites were identified to be under positive selection. Overall, our results suggest that following gene duplications, paralogous copies of -globin genes evolved under a nonepisodic process of functional divergence.[Reviewing Editor: Martin Kreitman]     

Why Is the Correlation between Gene Importance and Gene Evolutionary Rate So Weak?

One of the few commonly believed principles of molecular evolution is that functionally more important genes (or DNA sequences) evolve more slowly than less important ones. This principle is widely used by molecular biologists in daily practice. However, recent genomic analysis of a diverse array of organisms found only weak, negative correlations between the evolutionary rate of a gene and its functional importance, typically measured under a single benign lab condition. A frequently suggested cause of the above finding is that gene importance determined in the lab differs from that in an organism's natural environment. Here, we test this hypothesis in yeast using gene importance values experimentally determined in 418 lab conditions or computationally predicted for 10,000 nutritional conditions. In no single condition or combination of conditions did we find a much stronger negative correlation, which is explainable by our subsequent finding that always-essential (enzyme) genes do not evolve significantly more slowly than sometimes-essential or always-nonessential ones. Furthermore, we verified that functional density, approximated by the fraction of amino acid sites within protein domains, is uncorrelated with gene importance. Thus, neither the lab-nature mismatch nor a potentially biased among-gene distribution of functional density explains the observed weakness of the correlation between gene importance and evolutionary rate. We conclude that the weakness is factual, rather than artifactual. In addition to being weakened by population genetic reasons, the correlation is likely to have been further weakened by the presence of multiple nontrivial rate determinants that are independent from gene importance. These findings notwithstanding, we show that the principle of slower evolution of more important genes does have some predictive power when genes with vastly different evolutionary rates are compared, explaining why the principle can be practically useful despite the weakness of the correlation.     

Interaction between the Sbcc Gene of Escherichia Coli and the Gam Gene of Phage λ

gam mutants of phage lambda carrying long palindromes fail to form plaques on wild-type Escherichia coli but do grow on strains that are mutant in the sbcC gene. gam + lambda carrying the same palindrome grow on both hosts and on a host deleted for the recB, C and D genes. These results suggest that the Gam protein of lambda, known to interact also with E. coli's recBCD protein, can interact with the product of the sbcC gene.     

An Immune-Related Gene Evolved into the Master Sex-Determining Gene in Rainbow Trout, Oncorhynchus mykiss

Since the discovery of Sry in mammals [1, 2], few other master sex-determining genes have been identified in vertebrates [3-7]. To date, all of these genes have been characterized as well-known factors in the sex differentiation pathway, suggesting that the same subset of genes have been repeatedly and independently selected throughout evolution as master sex determinants [8, 9]. Here, we characterized in rainbow trout an unknown gene expressed only in the testis, with a predominant expression during testicular differentiation. This gene is a male-specific genomic sequence that is colocalized along with the sex-determining locus. This gene, named sdY for sexually dimorphic on the Y?chromosome, encodes a protein that displays similarity to the C-terminal domain of interferon regulatory factor 9. The targeted inactivation of sdY in males using zinc-finger nuclease induces ovarian differentiation, and the overexpression of sdY in females using additive transgenesis induces testicular differentiation. Together, these results demonstrate that sdY is a novel vertebrate master sex-determining gene not related to any known sex-differentiating gene. These findings highlight an unexpected evolutionary plasticity in vertebrate sex determination through the demonstration that master sex determinants can arise from the de novo evolution of genes that have not been previously implicated in sex differentiation.     

A Bayesian Approach to Pathway Analysis by Integrating Gene–Gene Functional Directions and Microarray Data

Many statistical methods have been developed to screen for differentially expressed genes associated with specific phenotypes in the microarray data. However, it remains a major challenge to synthesize the observed expression patterns with abundant biological knowledge for more complete understanding of the biological functions among genes. Various methods including clustering analysis on genes, neural network, Bayesian network and pathway analysis have been developed toward this goal. In most of these procedures, the activation and inhibition relationships among genes have hardly been utilized in the modeling steps. We propose two novel Bayesian models to integrate the microarray data with the putative pathway structures obtained from the KEGG database and the directional gene–gene interactions in the medical literature. We define the symmetric Kullback–Leibler divergence of a pathway, and use it to identify the pathway(s) most supported by the microarray data. Monte Carlo Markov Chain sampling algorithm is given for posterior computation in the hierarchical model. The proposed method is shown to select the most supported pathway in an illustrative example. Finally, we apply the methodology to a real microarray data set to understand the gene expression profile of osteoblast lineage at defined stages of differentiation. We observe that our method correctly identifies the pathways that are reported to play essential roles in modulating bone mass.     

Transformation of Pea Lectin Gene and Parasponia Haemoglobin Gene into Rice and Their Expressions

Lectin and leghemoglobin in legumes play the important roles, respectively, in recognition of host plants to their rhizobial bacteria, and lowering the oxygen partial pressure around bacteroids and protecting nitrogenase from oxygen in symbiotic nitrogen-fixing nodules. In order to extend the host range of the rhizobial bacteria and to make them fix nitrogen in non-legumes, pea lectin gene ( pl ) and Parasponia hemoglobin gene ( phb ) have been constructed into a plant expression vector (pCBHUL) and the vector pCBHUL was introduced into rice calli from immature young embryos by particle bombardment. After the calli were regenerated into plantlets on the resistant-selecting media containing hygromycin, they were identified by PCR and Southern blot hybridization. It was indicated that the pl and phb genes were integrated into nucleic genome of the transformed rice plants. GUS activity and the product of the pl gene were determined by GUS staining, Western blot and in situ hybridization at translational level. Eighteen out of 40 plants resistant to hygromycin were positively identified by PCR analysis with the rate of 45%. The pl gene was expressed in 3 out of 18 plants with 17% and 7.5%in 40 plants. The results may provide a clue for exploring whether Rhizobium leguminosarum bv. viceae could extend its host range and make the transgenic rice plants have the possibility of being symbiotic, or associative to nitrogen fixation.     

Gene–Environment Interactions at Nucleotide Resolution

Interactions among genes and the environment are a common source of phenotypic variation. To characterize the interplay between genetics and the environment at single nucleotide resolution, we quantified the genetic and environmental interactions of four quantitative trait nucleotides (QTN) that govern yeast sporulation efficiency. We first constructed a panel of strains that together carry all 32 possible combinations of the 4 QTN genotypes in 2 distinct genetic backgrounds. We then measured the sporulation efficiencies of these 32 strains across 8 controlled environments. This dataset shows that variation in sporulation efficiency is shaped largely by genetic and environmental interactions. We find clear examples of QTN:environment, QTN: background, and environment:background interactions. However, we find no QTN:QTN interactions that occur consistently across the entire dataset. Instead, interactions between QTN only occur under specific combinations of environment and genetic background. Thus, what might appear to be a QTN:QTN interaction in one background and environment becomes a more complex QTN:QTN:environment:background interaction when we consider the entire dataset as a whole. As a result, the phenotypic impact of a set of QTN alleles cannot be predicted from genotype alone. Our results instead demonstrate that the effects of QTN and their interactions are inextricably linked both to genetic background and to environmental variation.     

Gene therapy in Parkinson’s disease

Gene therapy in Parkinsons disease appears to be at the brink of the clinical study phase. Future gene therapy protocols will be based on a substantial amount of preclinical data regarding the use of ex vivo and in vivo genetic modifications with the help of viral or non-viral vectors. To date, the supplementation of neurotrophic factors and substitution for the dopaminergic deficit have formed the focus of trials to achieve relief in animal models of Parkinsons disease. Newer approaches include attempts to influence detrimental cell signalling pathways and to inhibit overactive basal ganglia structures. Nevertheless, current models of Parkinsons disease do not mirror all aspects of the human disease, and important issues with respect to long-term protein expression, choice of target structures and transgenes and safety remain to be solved. Here, we thoroughly review available animal data of gene transfer in models of Parkinsons disease.     

Small RNAs and Gene Network in a Durable Disease Resistance Gene—Mediated Defense Responses in Rice

Accumulating data have suggested that small RNAs (sRNAs) have important functions in plant responses to pathogen invasion. However, it is largely unknown whether and how sRNAs are involved in the regulation of rice responses to the invasion of Xanthomonas oryzae pv. oryzae (Xoo), which causes bacterial blight, the most devastating bacterial disease of rice worldwide. We performed simultaneous genome-wide analyses of the expression of sRNAs and genes during early defense responses of rice to Xoo mediated by a major disease resistance gene, Xa3/Xa26, which confers durable and race-specific qualitative resistance. A large number of sRNAs and genes showed differential expression in Xa3/Xa26-mediated resistance. These differentially expressed sRNAs include known microRNAs (miRNAs), unreported miRNAs, and small interfering RNAs. The candidate genes, with expression that was negatively correlated with the expression of sRNAs, were identified, indicating that these genes may be regulated by sRNAs in disease resistance in rice. These results provide a new perspective regarding the putative roles of sRNA candidates and their putative target genes in durable disease resistance in rice.     

Ultrasound Enhances Retrovirus‐Mediated Gene Transfer

Abstract

Viral vector systems are efficient for transfection of foreign genes into many tissues. Especially, retrovirus based vectors integrate the transgene into the genome of the target cells, which can sustain long term expression. However, it has been demonstrated that the transduction efficiency using retrovirus is relatively lower than those of other viruses. Ultrasound was recently reported to increase gene expression using plasmid DNA, with or without, a delivery vehicle. However, there are no reports, which show an ultrasound effect to retrovirus‐mediated gene transfer efficiency.

Retrovirus‐mediated gene transfer systems were used for transfection of 293T cells, bovine aortic endothelial cells (BAECs), rat aortic smooth muscle cells (RASMCs), and rat skeletal muscle myoblasts (L6 cells) with β‐galactosidase (β‐Gal) genes. Transduction efficiency and cell viability assay were performed on 293T cells that were exposed to varying durations (5 to 30 seconds) and power levels (1.0 watts/cm² to 4.0 watts/cm²) of ultrasound after being transduced by a retrovirus. Effects of ultrasound to the retrovirus itself was evaluated by transduction efficiency of 293T cells. After exposure to varying power levels of ultrasound to a retrovirus for 5 seconds, 293T cells were transduced by a retrovirus, and transduction efficiency was evaluated.

Below 1.0 watts/cm² and 5 seconds exposure, ultrasound showed increased transduction efficiency and no cytotoxicity to 293T cells transduced by a retrovirus. Also, ultrasound showed no toxicity to the virus itself at the same condition. Exposure of 5 seconds at the power of 1.0 watts/cm² of an ultrasound resulted in significant increases in retrovirus‐mediated gene expression in all four cell types tested in this experiment. Transduction efficiencies by ultrasound were enhanced 6.6‐fold, 4.8‐fold, 2.3‐fold, and 3.2‐fold in 293T cells, BAECs, RASMCs, and L6 cells, respectively. Furthermore, β‐Gal activities were also increased by the retrovirus with ultrasound exposure in these cells.

Adjunctive ultrasound exposure was associated with enhanced retrovirus‐mediated transgene expression in vitro. Ultrasound associated local gene therapy has potential for not only plasmid‐DNA‐, but also retrovirus‐mediated gene transfer.     

Deletion Mapping of the λ REX Gene

Deletion mapping has been used to order 12 λ rex^- mutants. Correlation of recombination data with physically-determined positions of deletion end-points (Szybalski 1971; Blattner et al. 1972) suggests that the left-most rex^- mutation, rex209, is located about 260-300 nucleotide pairs from the p_L mutation sex1 and about 475 nucleotide pairs from the left end-point of the region of nonhomology with λimm434.     

Gene expression profiling – Clusters of possibilities

Advances in qPCR technology allow studies of increasingly large systems comprising many genes and samples. The increasing data sizes allow expression profiling both in the gene and the samples dimension while also putting higher demands on sound statistical analysis and expertise to handle and interpret its results. We distinguish between exploratory and confirmatory statistical studies. In this paper we demonstrate several techniques available for exploratory studies on a system of Xenopus laevis development from egg to tadpole. Techniques include hierarchical clustering, heatmap, principal component analysis and self-organizing maps. We stress that even though exploratory studies are excellent for generating hypotheses, results have not been proven statistically significant until an independent confirmatory study has been performed. An exploratory study may certainly be valuable in its own right, and there are often not enough resources to report both an exploratory and a confirmatory study at the same time. However, exploratory and confirmatory studies are intimately connected and we would like to raise that awareness among qPCR practitioners. We suggest that scientific reports should always have a hypothesis focus. Reports are either hypothesis generating, from an exploratory study, or hypothesis validating, from a confirmatory study, or both. In either case, we suggest the generated or validated hypotheses be specifically stated.     

Comparison of β-Adrenergic and Glucocorticoid Signaling on Clock Gene and Osteoblast-Related Gene Expressions in Human Osteoblast

Most living organisms exhibit circadian rhythms that are generated by endogenous circadian clocks, the master one being present in the suprachiasmatic nuclei (SCN). Output signals from the SCN are believed to transmit standard circadian time to peripheral tissue through sympathetic nervous system and humoral routes. Therefore, the authors examined the expression of clock genes following treatment with the β-adrenergic receptor agonist, isoprenaline, or the synthetic glucocorticoid, dexamethasone, in cultured human osteoblast SaM-1 cells. Cells were treated with 10^?6 M isoprenaline or 10^?7 M dexamethasone for 2?h and gene expressions were determined using real-time polymerase chain reaction (PCR) analysis. Treatment with isoprenaline or dexamethasone induced the circadian expression of clock genes human period 1 (hPer1), hPer2, hPer3, and human brain and muscle Arnt-like protein 1 (hBMAL1). Isoprenaline or dexamethasone treatment immediately increased hPer1 and hPer2 and caused circadian oscillation of hPer1 and hPer2 with three peaks within 48?h. hPer3 expression had one peak after isoprenaline or dexamethasone treatment. hBMAL expression had two peaks after isoprenaline or dexamethasone treatment, the temporal pattern being in antiphase to that of the other clock genes. Dexamethasone treatment delayed the oscillation of all clock genes for 2–6?h compared with isoprenaline treatment. The authors also examined the expression of osteoblast-related genes hα-1 type I collagen (hCol1a1), halkaline phosphatase (hALP), and hosteocalcin (hOC). Isoprenaline induced oscillation of hCol1a1, but not hALP and hOC. On the other hand, dexamethasone induced oscillation of hCol1a1 and hALP, but not hOC. Isoprenaline up-regulated hCol1a1 expression, but dexamethasone down-regulated hCol1a1 and hALP expression in the first phase. (Author correspondence: togariaf@dpc.agu.ac.jp)     

Vitamin D Receptor Gene and Aggrecan Gene Polymorphisms and the Risk of Intervertebral Disc Degeneration — A Meta-Analysis

Background

A series of studies have been conducted to evaluate the associations between vitamin D receptor (VDR) and aggrecan variable numbers of tandem repeat (VNTR) polymorphisms and the risk of intervertebral disc degeneration (IDD), but produced conflicting results.

Objective

we performed a meta-analysis to address a more accurate estimation of the associations between the above gene polymorphisms and the risk of IDD.

Methods

A comprehensive literature search was conducted to identify all the relevant studies. The fixed or random effect model was selected based on the heterogeneity test among studies evaluated using the I ². Publication bias was estimated using Begg''s funnel plots and Egger''s regression test.

Results

A total of 9, 5, 3, and 7 studies were finally included in the analyses for the associations between the VDR TaqI (rs731236), FokI (rs2228570), ApaI (rs7975232), or aggrecan VNTR polymorphisms and the risk of IDD, respectively. The combined results showed that none of the VDR (TaqI, FokI, ApaI) polymorphisms were significantly associated with the risk of IDD. In contrast, the alleles with shorter VNTR length was found to significantly increase the risk of IDD (≦25 vs. >25: OR = 1.850, 95%CI 1.477–2.318; ≦23 vs. >23: OR = 1.955, 95%CI 1.41–2.703). Subgroup analysis confirmed the above results. After excluding studies deviated from Hardy-Weinberg equilibrium (HWE) in controls, no other studies were found to significantly influence the pooled effects in each genetic model. No potential publication bias was detected.

Conclusion

This meta-analysis suggested that the alleles with shorter VNTR length significantly increased the risk of IDD, while the VDR (TaqI, FokI, ApaI) gene polymorphisms were not significantly associated with the risk of IDD. Since potential confounders could not be ruled out completely, further studies are needed to confirm these results.