Identifying regulational alterations in gene regulatory networks by state space representation of vector autoregressive models and variational annealing

Background: In the analysis of effects by cell treatment such as drug dosing, identifying changes on gene network structures between normal and treated cells is a key task. A possible way for identifying the changes is to compare structures of networks estimated from data on normal and treated cells separately. However, this approach usually fails to estimate accurate gene networks due to the limited length of time series data and measurement noise. Thus, approaches that identify changes on regulations by using time series data on both conditions in an efficient manner are demanded. Methods: We propose a new statistical approach that is based on the state space representation of the vector autoregressive model and estimates gene networks on two different conditions in order to identify changes on regulations between the conditions. In the mathematical model of our approach, hidden binary variables are newly introduced to indicate the presence of regulations on each condition. The use of the hidden binary variables enables an efficient data usage; data on both conditions are used for commonly existing regulations, while for condition specific regulations corresponding data are only applied. Also, the similarity of networks on two conditions is automatically considered from the design of the potential function for the hidden binary variables. For the estimation of the hidden binary variables, we derive a new variational annealing method that searches the configuration of the binary variables maximizing the marginal likelihood. Results: For the performance evaluation, we use time series data from two topologically similar synthetic networks, and confirm that our proposed approach estimates commonly existing regulations as well as changes on regulations with higher coverage and precision than other existing approaches in almost all the experimental settings. For a real data application, our proposed approach is applied to time series data from normal Human lung cells and Human lung cells treated by stimulating EGF-receptors and dosing an anticancer drug termed Gefitinib. In the treated lung cells, a cancer cell condition is simulated by the stimulation of EGF-receptors, but the effect would be counteracted due to the selective inhibition of EGF-receptors by Gefitinib. However, gene expression profiles are actually different between the conditions, and the genes related to the identified changes are considered as possible off-targets of Gefitinib. Conclusions: From the synthetically generated time series data, our proposed approach can identify changes on regulations more accurately than existing methods. By applying the proposed approach to the time series data on normal and treated Human lung cells, candidates of off-target genes of Gefitinib are found. According to the published clinical information, one of the genes can be related to a factor of interstitial pneumonia, which is known as a side effect of Gefitinib.

Global pharmacogenomics: Impact of population diversity on the distribution of polymorphisms in the CYP2C cluster among Brazilians

The impact of biogeographical ancestry, self-reported 'race/color' and geographical origin on the frequency distribution of 10 CYP2C functional polymorphisms (CYP2C8*2, *3, *4, CYP2C9*2, *3, *5, *11, CYP2C19*2, *3 and *17) and their haplotypes was assessed in a representative cohort of the Brazilian population (n = 1034). TaqMan assays were used for allele discrimination at each CYP2C locus investigated. Individual proportions of European, African and Amerindian biogeographical ancestry were estimated using a panel of insertion-deletion polymorphisms. Multinomial log-linear models were applied to infer the statistical association between the CYP2C alleles and haplotypes (response variables), and biogeographical ancestry, self-reported Color and geographical origin (explanatory variables). The results showed that CYP2C19*3, CYP2C9*5 and CYP2C9*11 were rare alleles (<1%), the frequency of other variants ranged from 3.4% (CYP2C8*4) to 17.3% (CYP2C19*17). Two distinct haplotype blocks were identified: block 1 consists of three single nucleotide polymorphisms (SNPs) (CYP2C19*17, CYP2C19*2 and CYP2C9*2) and block 2 of six SNPs (CYP2C9*11, CYP2C9*3, CYP2C9*5, CYP2C8*2, CYP2C8*4 and CYP2C8*3). Diplotype analysis generated 41 haplotypes, of which eight had frequencies greater than 1% and together accounted for 96.4% of the overall genetic diversity. The distribution of CYP2C8 and CYP2C9 (but not CYP2C19) alleles, and of CYP2C haplotypes was significantly associated with self-reported Color and with the individual proportions of European and African genetic ancestry, irrespective of Color self-identification. The individual odds of having alleles CYP2C8*2, CYP2C8*3, CYP2C9*2 and CYP2C9*3, and haplotypes including these alleles, varied continuously as the proportion of European ancestry increased. Collectively, these data strongly suggest that the intrinsic heterogeneity of the Brazilian population must be acknowledged in the design and interpretation of pharmacogenomic studies of the CYP2C cluster in order to avoid spurious conclusions based on improper matching of study cohorts. This conclusion extends to other polymorphic pharmacogenes among Brazilians, and most likely to other admixed populations of the Americas. The Pharmacogenomics Journal (2012) 12, 267-276; doi: 10.1038/tpj.2010.89; published online 21 December 2010

Variable expressivity of osteogenesis imperfecta in a Brazilian family due to p.G1079S mutation in the COL1A1 gene

Osteogenesis imperfecta (OI) is a Mendelian disease with genetic heterogeneity characterized by bone fragility, recurrent fractures, blue sclerae, and short stature, caused mostly by mutations in COL1A1 or COL1A2 genes, which encode the pro-alpha 1(I) and pro-alpha 2(I) chains of type I collagen, respectively. A Brazilian family that showed variable expression of autosomal dominant OI was identified and characterized. Scanning for mutations was carried out using SSCP and DNA sequence analysis. The missense mutation c.3235G>A was identified within exon 45 of the COL1A1 gene in a 16-year-old girl diagnosed as having OI type I; it resulted in substitution of a glycine residue (G) by a serine (S) at codon 1079 (p.G1079S). The proband's mother had the disease signs, but without bone fractures, as did five of nine uncles and aunts of the patient. All of them carried the mutation, which was excluded in four healthy brothers of the patient's mother. This is the first description in a Brazilian family with OI showing variable expression; only one among seven carriers for the c.3235G>A mutation developed bone fractures, the most striking clinical feature of this disease. This finding has a significant implication for prenatal diagnosis in OI disease.

Contribution and perspectives of quantitative genetics to plant breeding in Brazil

The purpose of this article is to show how quantitative genetics has contributed to the huge genetic progress obtained in plant breeding in Brazil in the last forty years. The information obtained through quantitative genetics has given Brazilian breeders the possibility of responding to innumerable questions in their work in a much more informative way, such as the use or not of hybrid cultivars, which segregating population to use, which breeding method to employ, alternatives for improving the efficiency of selection programs, and how to handle the data of progeny and/or cultivars evaluations to identify the most stable ones and thus improve recommendations.

Molecular evolution of a malaria resistance gene (DARC) in primates

Genes involved in host-pathogen interactions are often strongly affected by positive natural selection. The Duffy antigen, coded by the Duffy antigen receptor for chemokines (DARC) gene, serves as a receptor for Plasmodium vivax in humans and for Plasmodium knowlesi in some nonhuman primates. In the majority of sub-Saharan Africans, a nucleic acid variant in GATA-1 of the gene promoter is responsible for the nonexpression of the Duffy antigen on red blood cells and consequently resistance to invasion by P. vivax. The Duffy antigen also acts as a receptor for chemokines and is expressed in red blood cells and many other tissues of the body. Because of this dual role, we sequenced a 3,000-bp region encompassing the entire DARC gene as well as part of its 5' and 3' flanking regions in a phylogenetic sample of primates and used statistical methods to evaluate the nature of selection pressures acting on the gene during its evolution. We analyzed both coding and regulatory regions of the DARC gene. The regulatory analysis showed accelerated rates of substitution at several sites near known motifs. Our tests of positive selection in the coding region using maximum likelihood by branch sites and maximum likelihood by codon sites did not yield statistically significant evidence for the action of positive selection. However, the maximum likelihood test in which the gene was subdivided into different structural regions showed that the known binding region for P. vivax/P. knowlesi is under very different selective pressures than the remainder of the gene. In fact, most of the gene appears to be under strong purifying selection, but this is not evident in the binding region. We suggest that the binding region is under the influence of two opposing selective pressures, positive selection possibly exerted by the parasite and purifying selection exerted by chemokines.

Variable patterns of Y chromosome homology in Akodontini rodents (Sigmodontinae): a phylogenetic signal revealed by chromosome painting

The Akodontini is the second most speciose tribe of sigmodontine rodents, one of the most diverse groups of neotropical mammals. Molecular phylogenetic analyses are discordant regarding the interrelationships of genera, with low support for some clades. However, two clades are concordant, one (clade A) with Akodon sensu strictu (excluding Akodon serrensis), "Akodon" serrensis, Bibimys, Deltamys, Juscelinomys, Necromys, Oxymycterus, Podoxymys, Thalpomys and Thaptomys, and another (clade B) with Blarinomys, Brucepattersonius, Kunsia, Lenoxus and Scapteromys. Here, we present chromosome painting using Akodon paranaensis (APA) Y paint, after suppression of simple repetitive sequences, on ten Akodontini genera. Partial Y chromosome homology, in addition to the homology already reported on the Akodon genus, was detected on the Y chromosomes of "A." serrensis, Thaptomys, Deltamys, Necromys and Thalpomys and on Y and X chromosomes in Oxymycterus. In Blarinomys, Brucepattersonius, Scapteromys and Kunsia, no APA Y signal was observed using different hybridization conditions; APA X paint gave positive signals only on the X chromosome in all genera. The Y chromosome homology was variable in size and positioning among the species studied as follow: (1) whole acrocentric Y chromosome in Akodon and "A." serrensis, (2) Yp and pericentromeric region in submetacentric Y of Necromys and Thaptomys, (3) pericentromeric region in acrocentric Y of Deltamys, (4) distal Yq in the acrocentric Y chromosome of Thalpomys and (5) proximal Yq in the acrocentric Y and Xp in the basal clade A genus Oxymycterus. The results suggest that the homology involves pairing (pseudoautosomal) and additional regions that have undergone rearrangement during divergence. The widespread Y homology represents a phylogenetic signal in Akodontini that provides additional evidence supporting the monophyly of clade A. The findings also raise questions about the evolution of the pseudoautosomal region observed in Oxymycterus. The Y chromosomes of these closely related species seem to have undergone dynamic rearrangements, including restructuring and reduction of homologous segments. Furthermore, the changes observed may indicate progressive attrition of the Y chromosome in more distantly related species.

Is missing maxillary lateral incisor in complete cleft lip and palate a product of genetics or local environment?

Objective: To test the null hypothesis: Subjects with isolated complete unilateral cleft lip and palate (UCLP) show no differences in overall frequency of tooth agenesis (hypodontia), comparing a subsample with cleft-side maxillary lateral incisor (MxI2) agenesis to a subsample without cleftside MxI2 agenesis. Findings could clarify the origins of cleft-side MxI2 agenesis. Materials and Methods: Tooth agenesis was identified from dental radiographs of 141 subjects with UCLP. The UCLP cohort was segregated into four categories according to the status and location of MxI2 in the region of the unilateral cleft: group M: subjects with one tooth, located on the mesial side of the alveolar cleft; group D: subjects with one tooth, located on the distal side of the alveolar cleft; group MD: subjects with two teeth present, one mesial and one distal to the cleft; and group ABS: subjects with lateral incisor absent (agenesis) in the cleft area. Results: The null hypothesis was rejected. Among UCLP subjects, there was a twofold increase (P < .0008) in overall frequency of tooth agenesis outside the cleft region in a subsample with cleftside MxI2 agenesis (ABS), compared to a subsample presenting with no agenesis of the cleft-side MxI2 (M+D+MD). Conclusions: Cleft-side MxI2 agenesis in CLP subjects appears to be largely a genetically controlled anomaly associated with cleft development, rather than a collateral environmental consequence of the adjacent cleft defect, since increased hypodontia involving multiple missing teeth observed remote from a cleft clearly has a significant genetic basis. (Angle Orthod. 2012;82:959-963.)

Duplications of BHLHA9 are associated with ectrodactyly and tibia hemimelia inherited in non-Mendelian fashion

Background Split-hand/foot malformation (SHFM)-also known as ectrodactyly-is a congenital disorder characterised by severe malformations of the distal limbs affecting the central rays of hands and/or feet. A distinct entity termed SHFLD presents with SHFM and long bone deficiency. Mouse models suggest that a defect of the central apical ectodermal ridge leads to the phenotype. Although six different loci/mutations (SHFM1-6) have been associated with SHFM, the underlying cause in a large number of cases is still unresolved. Methods High resolution array comparative genomic hybridisation (CGH) was performed in patients with SHFLD to detect copy number changes. Candidate genes were further evaluated for expression and function during limb development by whole mount in situ hybridisation and morpholino knock-down experiments. Results Array CGH showed microduplications on chromosome 17p13.3, a locus previously associated with SHFLD. Detailed analysis of 17 families revealed that this copy number variation serves as a susceptibility factor for a highly variable phenotype with reduced penetrance, particularly in females. Compared to other known causes for SHFLD 17p duplications appear to be the most frequent cause of SHFLD. A similar to 11.8 kb minimal critical region was identified encompassing a single gene, BHLHA9, a putative basic loop helix transcription factor. Whole mount in situ hybridisation showed expression restricted to the limb bud mesenchyme underlying the apical ectodermal ridge in mouse and zebrafish embryos. Knock down of bhlha9 in zebrafish resulted in shortening of the pectoral fins. Conclusions Genomic duplications encompassing BHLHA9 are associated with SHFLD and non-Mendelian inheritance characterised by a high degree of non-penetrance with sex bias. Knock-down of bhlha9 in zebrafish causes severe reduction defects of the pectoral fin, indicating a role for this gene in limb development.

A familial case with interstitial 2q36 deletion: Variable phenotypic expression in full and mosaic state

Submicroscopic chromosomal anomalies play an important role in the etiology of craniofacial malformations, including midline facial defects with hypertelorism (MFDH). MFDH is a common feature combination in several conditions, of which Frontonasal Dysplasia is the most frequently encountered manifestation; in most cases the etiology remains unknown. We identified a parent to child transmission of a 6.2 Mb interstitial deletion of chromosome region 2q36.1q36.3 by array-CGH and confirmed by FISH and microsatellite analysis. The patient and her mother both presented an MFDH phenotype although the phenotype in the mother was much milder than her daughter. Inspection of haplotype segregation within the family of 2q36.1 region suggests that the deletion arose on a chromosome derived from the maternal grandfather. Evidences based on FISH, microsatellite and array-CGH analysis point to a high frequency mosaicism for presence of a deleted region 2q36 occurring in blood of the mother. The frequency of mosaicism in other tissues could not be determined. We here suggest that the milder phenotype observed in the proband's mother can be explained by the mosaic state of the deletion. This most likely arose by an early embryonic deletion in the maternal embryo resulting in both gonadal and somatic mosaicism of two cell lines, with and without the deleted chromosome. The occurrence of gonadal mosaicism increases the recurrence risk significantly and is often either underestimated or not even taken into account in genetic counseling where new mutation is suspected. (C) 2012 Elsevier Masson SAS. All rights reserved.

Statistical models for experimental model membranes.

Biological membranes are constituted from lipid bilayers and proteins. Investigation of protein-membrane interaction, essential for biological function of cells, must rest upon solid knowledge of lipid bilayer behavior. Thus, extensive studies of an experimental model for membranes, lipid bilayers in water solution, have been undertaken in the last decades. These systems present structural, thermal and electrical properties which depend on temperature, ionic strength or concentration. In this talk, we shall discuss statistical models for lipid bilayers, as well as the relation between their properties and results for properties of lipid dispersions investigated by the laboratories supervised by Teresa Lamy (IF-USP) and Amando Ito (FFCL-USP).