The gene transformer-2 of Anastrepha fruit flies (Diptera, Tephritidae) and its evolution in insects

Background: In the tephritids Ceratitis, Bactrocera and Anastrepha, the gene transformer provides the memory device for sex determination via its auto-regulation; only in females is functional Tra protein produced. To date, the isolation and characterisation of the gene transformer-2 in the tephritids has only been undertaken in Ceratitis, and it has been shown that its function is required for the female-specific splicing of doublesex and transformer pre-mRNA. It therefore participates in transformer auto-regulatory function. In this work, the characterisation of this gene in eleven tephritid species belonging to the less extensively analysed genus Anastrepha was undertaken in order to throw light on the evolution of transformer-2. Results: The gene transformer-2 produces a protein of 249 amino acids in both sexes, which shows the features of the SR protein family. No significant partially spliced mRNA isoform specific to the male germ line was detected, unlike in Drosophila. It is transcribed in both sexes during development and in adult life, in both the soma and germ line. The injection of Anastrepha transformer-2 dsRNA into Anastrepha embryos caused a change in the splicing pattern of the endogenous transformer and doublesex pre-mRNA of XX females from the female to the male mode. Consequently, these XX females were transformed into pseudomales. The comparison of the eleven Anastrepha Transformer-2 proteins among themselves, and with the Transformer-2 proteins of other insects, suggests the existence of negative selection acting at the protein level to maintain Transformer-2 structural features. Conclusions: These results indicate that transformer-2 is required for sex determination in Anastrepha through its participation in the female-specific splicing of transformer and doublesex pre-mRNAs. It is therefore needed for the auto-regulation of the gene transformer. Thus, the transformer/transfomer-2 > doublesex elements at the bottom of the cascade, and their relationships, probably represent the ancestral state ( which still exists in the Tephritidae, Calliphoridae and Muscidae lineages) of the extant cascade found in the Drosophilidae lineage ( in which tra is just another component of the sex determination gene cascade regulated by Sex-lethal). In the phylogenetic lineage that gave rise to the drosophilids, evolution co-opted for Sex-lethal, modified it, and converted it into the key gene controlling sex determination.

CONTRASTING PHYLOGEOGRAPHIC PATTERNS IN MITOCHONDRIAL DNA AND MICROSATELLITES: EVIDENCE OF FEMALE PHILOPATRY AND MALE-BIASED GENE FLOW AMONG REGIONAL POPULATIONS OF THE BLUE-AND-YELLOW MACAW (PSITTACIFORMES: ARA ARARAUNA) IN BRAZIL

Comparing the patterns of population differentiation among genetic markers with different modes of inheritance call provide insights into patterns of sex-biased dispersal and gene flow. The blue-and-yellow Macaw (Ara ararauna) is a Neotropical parrot with a broad geographic distribution ill South America. However, little is known about the natural history and current status Of remaining wild populations, including levels of genetic variability. The progressive decline and possible fragmentation of populations may endanger this species in the near future. We analyzed mitochondrial DNA (mtDNA) control-region sequences and six microsatellite 106 Of Blue-and-yellow Macaws sampled throughout their geographic range ill Brazil to describe population genetic Structure, to make inferences about historical demography and dispersal behavior, and to provide insight for conservation efforts. Analyses of population genetic structure based on mtDNA showed evidence of two major populations ill western and eastern Brazil that share a few low-frequency haplotypes. This phylogeographic pattern seems to have originated by the historical isolation of Blue-and-yellow Macaw populations similar to 374,000 years ago and has been maintained by restricted gene flow and female philopatry. By contrast, variation ill biparentally inherited microsatellites was not structured geographically, Male-biased dispersal and female philopatry best explain the different patterns observed in these two markers. Because females disperse less than males, the two regional populations with well-differentiated mtDNA haplogroups should be considered two different management units for conservation purposes. Received 4 November 2007 accepted 10 December 2008.

A novel COL1A1 gene-splicing mutation (c. 1875+1G > C) in a Brazilian patient with osteogenesis imperfecta

Osteogenesis imperfecta is a heterogeneous genetic disorder characterized by bone fragility and deformity, recurrent fractures, blue sclera, short stature, and dentinogenesis imperfecta. Most cases are caused by mutations in COL1A1 and COL1A2 genes. We present a novel splicing mutation in the COL1A1 gene (c. 1875+ 1G>C) in a 16-year-old Brazilian boy diagnosed as a type III osteogenesis imperfecta patient. This splicing mutation and its association with clinical phenotypes will be submitted to the reference database of COL1A1 mutations, which has no other description of this mutation.

Deletion of the Basement Membrane Heparan Sulfate Proteoglycan Type XVIII Collagen Causes Hypertriglyceridemia in Mice and Humans

Background: Lipoprotein lipase (Lpl) acts on triglyceride-rich lipoproteins in the peripheral circulation, liberating free fatty acids for energy metabolism or storage. This essential enzyme is synthesized in parenchymal cells of adipose tissue, heart, and skeletal muscle and migrates to the luminal side of the vascular endothelium where it acts upon circulating lipoproteins. Prior studies suggested that Lpl is immobilized by way of heparan sulfate proteoglycans on the endothelium, but genetically altering endothelial cell heparan sulfate had no effect on Lpl localization or lipolysis. The objective of this study was to determine if extracellular matrix proteoglycans affect Lpl distribution and triglyceride metabolism. Methods and Findings: We examined mutant mice defective in collagen XVIII (Col18), a heparan sulfate proteoglycan present in vascular basement membranes. Loss of Col18 reduces plasma levels of Lpl enzyme and activity, which results in mild fasting hypertriglyceridemia and diet-induced hyperchylomicronemia. Humans with Knobloch Syndrome caused by a null mutation in the vascular form of Col18 also present lower than normal plasma Lpl mass and activity and exhibit fasting hypertriglyceridemia. Conclusions: This is the first report demonstrating that Lpl presentation on the lumenal side of the endothelium depends on a basement membrane proteoglycan and demonstrates a previously unrecognized phenotype in patients lacking Col18.

Expression of an activating mutation in the gene encoding the K(ATP) channel subunit Kir6.2 in mouse pancreatic beta cells recapitulates neonatal diabetes

Neonatal diabetes is a rare monogenic form of diabetes that usually presents within the first six months of life. It is commonly caused by gain-of-function mutations in the genes encoding the Kir6.2 and SUR1 subunits of the plasmalemmal ATP-sensitive K(+) (K(ATP)) channel. To better understand this disease, we generated a mouse expressing a Kir6.2 mutation (V59M) that causes neonatal diabetes in humans and we used Cre-lox technology to express the mutation specifically in pancreatic beta cells. These beta-V59M mice developed severe diabetes soon after birth, and by 5 weeks of age, blood glucose levels were markedly increased and insulin was undetectable. Islets isolated from beta-V59M mice secreted substantially less insulin and showed a smaller increase in intracellular calcium in response to glucose. This was due to a reduced sensitivity of K(ATP) channels in pancreatic beta cells to inhibition by ATP or glucose. In contrast, the sulfonylurea tolbutamide, a specific blocker of K(ATP) channels, closed K(ATP) channels, elevated intracellular calcium levels, and stimulated insulin release in beta-V59M beta cells, indicating that events downstream of K(ATP) channel closure remained intact. Expression of the V59M Kir6.2 mutation in pancreatic beta cells alone is thus sufficient to recapitulate the neonatal diabetes observed in humans. beta-V59M islets also displayed a reduced percentage of beta cells, abnormal morphology, lower insulin content, and decreased expression of Kir6.2, SUR1, and insulin mRNA. All these changes are expected to contribute to the diabetes of beta-V59M mice. Their cause requires further investigation.

Brabykinin B1 Receptor Antagonism Is Beneficial in Renal Ischemia-Reperfusion Injury

Previously we have demonstrated that bradykinin B1 receptor deficient mice (B1KO) were protected against renal ischemia and reperfusion injury (IRI). Here, we aimed to analyze the effect of B1 antagonism on renal IRI and to study whether B1R knockout or antagonism could modulate the renal expression of pro and anti-inflammatory molecules. To this end, mice were subjected to 45 minutes ischemia and reperfused at 4, 24, 48 and 120 hours. Wild-type mice were treated intra-peritoneally with antagonists of either B1 (R-954, 200 mg/kg) or B2 receptor (HOE140, 200 mg/kg) 30 minutes prior to ischemia. Blood samples were collected to ascertain serum creatinine level, and kidneys were harvested for gene transcript analyses by real-time PCR. Herein, B1R antagonism ( R-954) was able to decrease serum creatinine levels, whereas B2R antagonism had no effect. The protection seen under B1R deletion or antagonism was associated with an increased expression of GATA-3, IL-4 and IL-10 and a decreased T-bet and IL-1b transcription. Moreover, treatment with R-954 resulted in lower MCP-1, and higher HO-1 expression. Our results demonstrated that bradykinin B1R antagonism is beneficial in renal IRI.

Positive Selection Results in Frequent Reversible Amino Acid Replacements in the G Protein Gene of Human Respiratory Syncytial Virus

Human respiratory syncytial virus (HRSV) is the major cause of lower respiratory tract infections in children under 5 years of age and the elderly, causing annual disease outbreaks during the fall and winter. Multiple lineages of the HRSVA and HRSVB serotypes co-circulate within a single outbreak and display a strongly temporal pattern of genetic variation, with a replacement of dominant genotypes occurring during consecutive years. In the present study we utilized phylogenetic methods to detect and map sites subject to adaptive evolution in the G protein of HRSVA and HRSVB. A total of 29 and 23 amino acid sites were found to be putatively positively selected in HRSVA and HRSVB, respectively. Several of these sites defined genotypes and lineages within genotypes in both groups, and correlated well with epitopes previously described in group A. Remarkably, 18 of these positively selected tended to revert in time to a previous codon state, producing a ""flipflop'' phylogenetic pattern. Such frequent evolutionary reversals in HRSV are indicative of a combination of frequent positive selection, reflecting the changing immune status of the human population, and a limited repertoire of functionally viable amino acids at specific amino acid sites.

Gene Expression Noise in Spatial Patterning: hunchback Promoter Structure Affects Noise Amplitude and Distribution in Drosophila Segmentation

Positional information in developing embryos is specified by spatial gradients of transcriptional regulators. One of the classic systems for studying this is the activation of the hunchback (hb) gene in early fruit fly (Drosophila) segmentation by the maternally-derived gradient of the Bicoid (Bcd) protein. Gene regulation is subject to intrinsic noise which can produce variable expression. This variability must be constrained in the highly reproducible and coordinated events of development. We identify means by which noise is controlled during gene expression by characterizing the dependence of hb mRNA and protein output noise on hb promoter structure and transcriptional dynamics. We use a stochastic model of the hb promoter in which the number and strength of Bcd and Hb (self-regulatory) binding sites can be varied. Model parameters are fit to data from WT embryos, the self-regulation mutant hb(14F), and lacZ reporter constructs using different portions of the hb promoter. We have corroborated model noise predictions experimentally. The results indicate that WT (self-regulatory) Hb output noise is predominantly dependent on the transcription and translation dynamics of its own expression, rather than on Bcd fluctuations. The constructs and mutant, which lack self-regulation, indicate that the multiple Bcd binding sites in the hb promoter (and their strengths) also play a role in buffering noise. The model is robust to the variation in Bcd binding site number across a number of fly species. This study identifies particular ways in which promoter structure and regulatory dynamics reduce hb output noise. Insofar as many of these are common features of genes (e. g. multiple regulatory sites, cooperativity, self-feedback), the current results contribute to the general understanding of the reproducibility and determinacy of spatial patterning in early development.

Molecular determinants of improved cathepsin B inhibition by new cystatins obtained by DNA shuffling

Background: Cystatins are inhibitors of cysteine proteases. The majority are only weak inhibitors of human cathepsin B, which has been associated with cancer, Alzheimer's disease and arthritis. Results: Starting from the sequences of oryzacystatin-1 and canecystatin-1, a shuffling library was designed and a hybrid clone obtained, which presented higher inhibitory activity towards cathepsin B. This clone presented two unanticipated point mutations as well as an N-terminal deletion. Reversing each point mutation independently or both simultaneously abolishes the inhibitory activity towards cathepsin B. Homology modeling together with experimental studies of the reverse mutants revealed the likely molecular determinants of the improved inhibitory activity to be related to decreased protein stability. Conclusion: A combination of experimental approaches including gene shuffling, enzyme assays and reverse mutation allied to molecular modeling has shed light upon the unexpected inhibitory properties of certain cystatin mutants against Cathepsin B. We conclude that mutations disrupting the hydrophobic core of phytocystatins increase the flexibility of the N-terminus, leading to an increase in inhibitory activity. Such mutations need not affect the inhibitory site directly but may be observed distant from it and manifest their effects via an uncoupling of its three components as a result of increased protein flexibility.

Identification of the Schistosoma mansoni TNF-Alpha Receptor Gene and the Effect of Human TNF-Alpha on the Parasite Gene Expression Profile

Background: Schistosoma mansoni is the major causative agent of schistosomiasis. The parasite takes advantage of host signals to complete its development in the human body. Tumor necrosis factor-alpha (TNF-alpha) is a human cytokine involved in skin inflammatory responses, and although its effect on the adult parasite's metabolism and egg-laying process has been previously described, a comprehensive assessment of the TNF-alpha pathway and its downstream molecular effects is lacking. Methodology/Principal Findings: In the present work we describe a possible TNF-alpha receptor (TNFR) homolog gene in S. mansoni (SmTNFR). SmTNFR encodes a complete receptor sequence composed of 599 amino acids, and contains four cysteine-rich domains as described for TNFR members. Real-time RT-PCR experiments revealed that SmTNFR highest expression level is in cercariae, 3.5 (+/- 0.7) times higher than in adult worms. Downstream members of the known human TNF-alpha pathway were identified by an in silico analysis, revealing a possible TNF-alpha signaling pathway in the parasite. In order to simulate parasite's exposure to human cytokine during penetration of the skin, schistosomula were exposed to human TNF-alpha just 3 h after cercariae-to-schistosomula in vitro transformation, and large-scale gene expression measurements were performed with microarrays. A total of 548 genes with significantly altered expression were detected, when compared to control parasites. In addition, treatment of adult worms with TNF-alpha caused a significantly altered expression of 1857 genes. Interestingly, the set of genes altered in adults is different from that of schistosomula, with 58 genes in common, representing 3% of altered genes in adults and 11% in 3 h-old early schistosomula. Conclusions/Significance: We describe the possible molecular elements and targets involved in human TNF-alpha effect on S. mansoni, highlighting the mechanism by which recently transformed schistosomula may sense and respond to this host mediator at the site of cercarial penetration into the skin.

Novel pathogenic mutations and skin biopsy analysis in Knobloch syndrome

Purpose: To facilitate future diagnosis of Knobloch syndrome (KS) and better understand its etiology, we sought to identify not yet described COL18A1 mutations in KS patients. In addition, we tested whether mutations in this gene lead to absence of the COL18A1 gene product and attempted to better characterize the functional effect of a previously reported missense mutation. Methods: Direct sequencing of COL18A1 exons was performed in KS patients from four unrelated pedigrees. We used immunofluorescent histochemistry in skin biopsies to evaluate the presence of type XVIII collagen in four KS patients carrying two already described mutations: c. 3277C>T, a nonsense mutation, and c. 3601G>A, a missense mutation. Furthermore, we determined the binding properties of the mutated endostatin domain p.A1381T (c.3601G>A) to extracellular matrix proteins using ELISA and surface plasmon resonance assays. Results: We identified four novel mutations in COL18A1, including a large deletion involving exon 41. Skin biopsies from KS patients revealed lack of type XVIII collagen in epithelial basement membranes and blood vessels. We also found a reduced affinity of p.A1381T endostatin to some extracellular matrix components. Conclusions: COL18A1 mutations involved in Knobloch syndrome have a distribution bias toward the coding exons of the C-terminal end. Large deletions must also be considered when point mutations are not identified in patients with characteristic KS phenotype. We report, for the first time, lack of type XVIII collagen in KS patients by immunofluorescent histochemistry in skin biopsy samples. As a final point, we suggest the employment of this technique as a preliminary and complementary test for diagnosis of KS in cases when mutation screening either does not detect mutations or reveals mutations of uncertain effect, such as the p.A1381T change.

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.

Identification of COL6A1 as a differentially expressed gene in human astrocytomas

Diffuse infiltrating gliomas are the most common tumors of the central nervous system. Gliomas are classified by the WHO according to their histopathological and clinical characteristics into four classes: grade I (pilocytic astrocytoma), grade II (diffuse astrocytoma), grade III (anaplastic astrocytoma), and grade IV (glioblastoma multiforme). Several genes have already been correlated with astrocytomas, but many others are yet to be uncovered. By analyzing the public SAGE data from 21 patients, comprising low malignant grade astrocytomas and glioblastomas, we found COL6A1 to be differentially expressed, confirming this finding by real time RT-PCR in 66 surgical samples. To the best of our knowledge, COL6A1 has never been described in gliomas. The expression of this gene has significantly different means when normal glia is compared with low-grade astrocytomas (grades I and II) and high-grade astrocytomas (grades III and IV), with a tendency to be greater in higher grade samples, thus rendering it a powerful tumor marker.

Inference of gene regulatory networks from time series by Tsallis entropy

Background: The inference of gene regulatory networks (GRNs) from large-scale expression profiles is one of the most challenging problems of Systems Biology nowadays. Many techniques and models have been proposed for this task. However, it is not generally possible to recover the original topology with great accuracy, mainly due to the short time series data in face of the high complexity of the networks and the intrinsic noise of the expression measurements. In order to improve the accuracy of GRNs inference methods based on entropy (mutual information), a new criterion function is here proposed. Results: In this paper we introduce the use of generalized entropy proposed by Tsallis, for the inference of GRNs from time series expression profiles. The inference process is based on a feature selection approach and the conditional entropy is applied as criterion function. In order to assess the proposed methodology, the algorithm is applied to recover the network topology from temporal expressions generated by an artificial gene network (AGN) model as well as from the DREAM challenge. The adopted AGN is based on theoretical models of complex networks and its gene transference function is obtained from random drawing on the set of possible Boolean functions, thus creating its dynamics. On the other hand, DREAM time series data presents variation of network size and its topologies are based on real networks. The dynamics are generated by continuous differential equations with noise and perturbation. By adopting both data sources, it is possible to estimate the average quality of the inference with respect to different network topologies, transfer functions and network sizes. Conclusions: A remarkable improvement of accuracy was observed in the experimental results by reducing the number of false connections in the inferred topology by the non-Shannon entropy. The obtained best free parameter of the Tsallis entropy was on average in the range 2.5 <= q <= 3.5 (hence, subextensive entropy), which opens new perspectives for GRNs inference methods based on information theory and for investigation of the nonextensivity of such networks. The inference algorithm and criterion function proposed here were implemented and included in the DimReduction software, which is freely available at http://sourceforge.net/projects/dimreduction and http://code.google.com/p/dimreduction/.

Construct and Compare Gene Coexpression Networks with DAPfinder and DAPview

Background: DAPfinder and DAPview are novel BRB-ArrayTools plug-ins to construct gene coexpression networks and identify significant differences in pairwise gene-gene coexpression between two phenotypes. Results: Each significant difference in gene-gene association represents a Differentially Associated Pair (DAP). Our tools include several choices of filtering methods, gene-gene association metrics, statistical testing methods and multiple comparison adjustments. Network results are easily displayed in Cytoscape. Analyses of glioma experiments and microarray simulations demonstrate the utility of these tools. Conclusions: DAPfinder is a new friendly-user tool for reconstruction and comparison of biological networks.