Distinct levels in Pom1 gradients limit Cdr2 activity and localization to time and position division.

Where and when cells divide are fundamental questions. In rod-shaped fission yeast cells, the DYRK-family kinase Pom1 is organized in concentration gradients from cell poles and controls cell division timing and positioning. Pom1 gradients restrict to mid-cell the SAD-like kinase Cdr2, which recruits Mid1/Anillin for medial division. Pom1 also delays mitotic commitment through Cdr2, which inhibits Wee1. Here, we describe quantitatively the distributions of cortical Pom1 and Cdr2. These reveal low profile overlap contrasting with previous whole-cell measurements and Cdr2 levels increase with cell elongation, raising the possibility that Pom1 regulates mitotic commitment by controlling Cdr2 medial levels. However, we show that distinct thresholds of Pom1 activity define the timing and positioning of division. Three conditions-a separation-of-function Pom1 allele, partial downregulation of Pom1 activity, and haploinsufficiency in diploid cells-yield cells that divide early, similar to pom1 deletion, but medially, like wild-type cells. In these cells, Cdr2 is localized correctly at mid-cell. Further, Cdr2 overexpression promotes precocious mitosis only in absence of Pom1. Thus, Pom1 inhibits Cdr2 for mitotic commitment independently of regulating its localization or cortical levels. Indeed, we show Pom1 restricts Cdr2 activity through phosphorylation of a C-terminal self-inhibitory tail. In summary, our results demonstrate that distinct levels in Pom1 gradients delineate a medial Cdr2 domain, for cell division placement, and control its activity, for mitotic commitment.

Using transcription modules to identify expression clusters perturbed in Williams-Beuren syndrome.

The genetic dissection of the phenotypes associated with Williams-Beuren Syndrome (WBS) is advancing thanks to the study of individuals carrying typical or atypical structural rearrangements, as well as in vitro and animal studies. However, little is known about the global dysregulations caused by the WBS deletion. We profiled the transcriptomes of skin fibroblasts from WBS patients and compared them to matched controls. We identified 868 differentially expressed genes that were significantly enriched in extracellular matrix genes, major histocompatibility complex (MHC) genes, as well as genes in which the products localize to the postsynaptic membrane. We then used public expression datasets from human fibroblasts to establish transcription modules, sets of genes coexpressed in this cell type. We identified those sets in which the average gene expression was altered in WBS samples. Dysregulated modules are often interconnected and share multiple common genes, suggesting that intricate regulatory networks connected by a few central genes are disturbed in WBS. This modular approach increases the power to identify pathways dysregulated in WBS patients, thus providing a testable set of additional candidates for genes and their interactions that modulate the WBS phenotypes.

Genome-wide association study identifies eight loci associated with blood pressure.

Elevated blood pressure is a common, heritable cause of cardiovascular disease worldwide. To date, identification of common genetic variants influencing blood pressure has proven challenging. We tested 2.5 million genotyped and imputed SNPs for association with systolic and diastolic blood pressure in 34,433 subjects of European ancestry from the Global BPgen consortium and followed up findings with direct genotyping (N ≤ 71,225 European ancestry, N ≤ 12,889 Indian Asian ancestry) and in silico comparison (CHARGE consortium, N = 29,136). We identified association between systolic or diastolic blood pressure and common variants in eight regions near the CYP17A1 (P = 7 × 10(-24)), CYP1A2 (P = 1 × 10(-23)), FGF5 (P = 1 × 10(-21)), SH2B3 (P = 3 × 10(-18)), MTHFR (P = 2 × 10(-13)), c10orf107 (P = 1 × 10(-9)), ZNF652 (P = 5 × 10(-9)) and PLCD3 (P = 1 × 10(-8)) genes. All variants associated with continuous blood pressure were associated with dichotomous hypertension. These associations between common variants and blood pressure and hypertension offer mechanistic insights into the regulation of blood pressure and may point to novel targets for interventions to prevent cardiovascular disease.

PHYTOCHROME KINASE SUBSTRATE1 regulates root phototropism and gravitropism.

Light promotes the expression of PHYTOCHROME KINASE SUBSTRATE1 (PKS1) in the root of Arabidopsis thaliana, but the function of PKS1 in this organ is unknown. Unilateral blue light induced a negative root phototropic response mediated by phototropin 1 in wild-type seedlings. This response was absent in pks1 mutants. In the wild type, unilateral blue light enhanced PKS1 expression in the subapical region of the root several hours before bending was detectable. The negative phototropism and the enhanced PKS1 expression in response to blue light required phytochrome A (phyA). In addition, the pks1 mutation enhanced the root gravitropic response when vertically oriented seedlings were placed horizontally. The negative regulation of gravitropism by PKS1 occurred even in dark-grown seedlings and did not require phyA. Blue light also failed to induce negative phototropism in pks1 under reduced gravitational stimulation, indicating that the effect of pks1 on phototropism is not simply the consequence of the counteracting effect of enhanced gravitropism. We propose a model where the background level of PKS1 reduces gravitropism. After a phyA-dependent increase in its expression, PKS1 positively affects root phototropism and both effects contribute to negative curvature in response to unilateral blue light.

Common variants in UMOD associate with urinary uromodulin levels: a meta-analysis.

Uromodulin is expressed exclusively in the thick ascending limb and is the most abundant protein excreted in normal urine. Variants in UMOD, which encodes uromodulin, are associated with renal function, and urinary uromodulin levels may be a biomarker for kidney disease. However, the genetic factors regulating uromodulin excretion are unknown. We conducted a meta-analysis of urinary uromodulin levels to identify associated common genetic variants in the general population. We included 10,884 individuals of European descent from three genetic isolates and three urban cohorts. Each study measured uromodulin indexed to creatinine and conducted linear regression analysis of approximately 2.5 million single nucleotide polymorphisms using an additive model. We also tested whether variants in genes expressed in the thick ascending limb associate with uromodulin levels. rs12917707, located near UMOD and previously associated with renal function and CKD, had the strongest association with urinary uromodulin levels (P<0.001). In all cohorts, carriers of a G allele of this variant had higher uromodulin levels than noncarriers did (geometric means 10.24, 14.05, and 17.67 μg/g creatinine for zero, one, or two copies of the G allele). rs12446492 in the adjacent gene PDILT (protein disulfide isomerase-like, testis expressed) also reached genome-wide significance (P<0.001). Regarding genes expressed in the thick ascending limb, variants in KCNJ1, SORL1, and CAB39 associated with urinary uromodulin levels. These data indicate that common variants in the UMOD promoter region may influence urinary uromodulin levels. They also provide insights into uromodulin biology and the association of UMOD variants with renal function.

Genetic polymorphisms of the main transcription factors for adiponectin gene promoter in regulation of adiponectin levels: association analysis in three European cohorts.

Adiponectin serum concentrations are an important biomarker in cardiovascular epidemiology with heritability etimates of 30-70%. However, known genetic variants in the adiponectin gene locus (ADIPOQ) account for only 2%-8% of its variance. As transcription factors are thought to play an under-acknowledged role in carrying functional variants, we hypothesized that genetic polymorphisms in genes coding for the main transcription factors for the ADIPOQ promoter influence adiponectin levels. Single nucleotide polymorphisms (SNPs) at these genes were selected based on the haplotype block structure and previously published evidence to be associated with adiponectin levels. We performed association analyses of the 24 selected SNPs at forkhead box O1 (FOXO1), sterol-regulatory-element-binding transcription factor 1 (SREBF1), sirtuin 1 (SIRT1), peroxisome-proliferator-activated receptor gamma (PPARG) and transcription factor activating enhancer binding protein 2 beta (TFAP2B) gene loci with adiponectin levels in three different European cohorts: SAPHIR (n = 1742), KORA F3 (n = 1636) and CoLaus (n = 5355). In each study population, the association of SNPs with adiponectin levels on log-scale was tested using linear regression adjusted for age, sex and body mass index, applying both an additive and a recessive genetic model. A pooled effect size was obtained by meta-analysis assuming a fixed effects model. We applied a significance threshold of 0.0033 accounting for the multiple testing situation. A significant association was only found for variants within SREBF1 applying an additive genetic model (smallest p-value for rs1889018 on log(adiponectin) = 0.002, β on original scale = -0.217 µg/ml), explaining ∼0.4% of variation of adiponectin levels. Recessive genetic models or haplotype analyses of the FOXO1, SREBF1, SIRT1, TFAPB2B genes or sex-stratified analyses did not reveal additional information on the regulation of adiponectin levels. The role of genetic variations at the SREBF1 gene in regulating adiponectin needs further investigation by functional studies.

Light intensity modulates the regulatory network of the shade avoidance response in Arabidopsis.

Plants such as Arabidopsis thaliana respond to foliar shade and neighbors who may become competitors for light resources by elongation growth to secure access to unfiltered sunlight. Challenges faced during this shade avoidance response (SAR) are different under a light-absorbing canopy and during neighbor detection where light remains abundant. In both situations, elongation growth depends on auxin and transcription factors of the phytochrome interacting factor (PIF) class. Using a computational modeling approach to study the SAR regulatory network, we identify and experimentally validate a previously unidentified role for long hypocotyl in far red 1, a negative regulator of the PIFs. Moreover, we find that during neighbor detection, growth is promoted primarily by the production of auxin. In contrast, in true shade, the system operates with less auxin but with an increased sensitivity to the hormonal signal. Our data suggest that this latter signal is less robust, which may reflect a cost-to-robustness tradeoff, a system trait long recognized by engineers and forming the basis of information theory.

Molecular biology of baculovirus and its use in biological control in Brazil

Baculoviruses are insect viruses found mainly in Lepidoptera. The family Baculoviridae is taxonomically divided in two genera, Nucleopolyhedrovirus and Granulovirus, which differ by occlusion body morphology. NPVs (Nucleopolyhedroviruses) have polyhedrical inclusion bodies (PIBs) containing multiple viral particles, while GVs (Granuloviruses) appear to be generally single particles occluded in oval shaped occlusion bodies. During the life cycle, two different viral progenies are produced: BV (Budded Virus) and PDV (Polyhedra Derived Virus), which are essential for the infectious process and virus propagation in host cells. Baculoviruses are being used for pest control and they are especially safe due to their specificity and invertebrate-restricted host range. Baculoviruses have been used as vectors for high level protein expression ofheterologous genes from prokaryotic and eukaryotic organisms. Also, recombinant DNA techniques have allowed the production of genetically modified viral insecticides. This study is a review on the taxonomy, structure, replication and molecular biology of baculoviruses, as well as their use as bioinsecticides in Brazil.

Common variants near MC4R are associated with fat mass, weight and risk of obesity.

To identify common variants influencing body mass index (BMI), we analyzed genome-wide association data from 16,876 individuals of European descent. After previously reported variants in FTO, the strongest association signal (rs17782313, P = 2.9 x 10(-6)) mapped 188 kb downstream of MC4R (melanocortin-4 receptor), mutations of which are the leading cause of monogenic severe childhood-onset obesity. We confirmed the BMI association in 60,352 adults (per-allele effect = 0.05 Z-score units; P = 2.8 x 10(-15)) and 5,988 children aged 7-11 (0.13 Z-score units; P = 1.5 x 10(-8)). In case-control analyses (n = 10,583), the odds for severe childhood obesity reached 1.30 (P = 8.0 x 10(-11)). Furthermore, we observed overtransmission of the risk allele to obese offspring in 660 families (P (pedigree disequilibrium test average; PDT-avg) = 2.4 x 10(-4)). The SNP location and patterns of phenotypic associations are consistent with effects mediated through altered MC4R function. Our findings establish that common variants near MC4R influence fat mass, weight and obesity risk at the population level and reinforce the need for large-scale data integration to identify variants influencing continuous biomedical traits.

Precision and scaling in morphogen gradient read-out.

Morphogen gradients infer cell fate as a function of cellular position. Experiments in Drosophila embryos have shown that the Bicoid (Bcd) gradient is precise and exhibits some degree of scaling. We present experimental results on the precision of Bcd target genes for embryos with a single, double or quadruple dose of bicoid demonstrating that precision is highest at mid-embryo and position dependent, rather than gene dependent. This confirms that the major contribution to precision is achieved already at the Bcd gradient formation. Modeling this dynamic process, we investigate precision for inter-embryo fluctuations in different parameters affecting gradient formation. Within our modeling framework, the observed precision can only be achieved by a transient Bcd profile. Studying different extensions of our modeling framework reveals that scaling is generally position dependent and decreases toward the posterior pole. Our measurements confirm this trend, indicating almost perfect scaling except for anterior most expression domains, which overcompensate fluctuations in embryo length.

Interactions of dietary whole-grain intake with fasting glucose- and insulin-related genetic loci in individuals of European descent: a meta-analysis of 14 cohort studies.

Whole-grain foods are touted for multiple health benefits, including enhancing insulin sensitivity and reducing type 2 diabetes risk. Recent genome-wide association studies (GWAS) have identified several single nucleotide polymorphisms (SNPs) associated with fasting glucose and insulin concentrations in individuals free of diabetes. We tested the hypothesis that whole-grain food intake and genetic variation interact to influence concentrations of fasting glucose and insulin. Via meta-analysis of data from 14 cohorts comprising ∼ 48,000 participants of European descent, we studied interactions of whole-grain intake with loci previously associated in GWAS with fasting glucose (16 loci) and/or insulin (2 loci) concentrations. For tests of interaction, we considered a P value <0.0028 (0.05 of 18 tests) as statistically significant. Greater whole-grain food intake was associated with lower fasting glucose and insulin concentrations independent of demographics, other dietary and lifestyle factors, and BMI (β [95% CI] per 1-serving-greater whole-grain intake: -0.009 mmol/l glucose [-0.013 to -0.005], P < 0.0001 and -0.011 pmol/l [ln] insulin [-0.015 to -0.007], P = 0.0003). No interactions met our multiple testing-adjusted statistical significance threshold. The strongest SNP interaction with whole-grain intake was rs780094 (GCKR) for fasting insulin (P = 0.006), where greater whole-grain intake was associated with a smaller reduction in fasting insulin concentrations in those with the insulin-raising allele. Our results support the favorable association of whole-grain intake with fasting glucose and insulin and suggest a potential interaction between variation in GCKR and whole-grain intake in influencing fasting insulin concentrations.

Comparative modular analysis of gene expression in vertebrate development

The focus of my PhD research was the concept of modularity. In the last 15 years, modularity has become a classic term in different fields of biology. On the conceptual level, a module is a set of interacting elements that remain mostly independent from the elements outside of the module. I used modular analysis techniques to study gene expression evolution in vertebrates. In particular, I identified ``natural'' modules of gene expression in mouse and human, and I showed that expression of organ-specific and system-specific genes tends to be conserved between such distance vertebrates as mammals and fishes. Also with a modular approach, I studied patterns of developmental constraints on transcriptome evolution. I showed that none of the two commonly accepted models of the evolution of embryonic development (``evo-devo'') are exclusively valid. In particular, I found that the conservation of the sequences of regulatory regions is highest during mid-development of zebrafish, and thus it supports the ``hourglass model''. In contrast, events of gene duplication and new gene introduction are most rare in early development, which supports the ``early conservation model''. In addition to the biological insights on transcriptome evolution, I have also discussed in detail the advantages of modular approaches in large-scale data analysis. Moreover, I re-analyzed several studies (published in high-ranking journals), and showed that their conclusions do not hold out under a detailed analysis. This demonstrates that complex analysis of high-throughput data requires a co-operation between biologists, bioinformaticians, and statisticians.

FastEpistasis: a high performance computing solution for quantitative trait epistasis.

Motivation: Genome-wide association studies have become widely used tools to study effects of genetic variants on complex diseases. While it is of great interest to extend existing analysis methods by considering interaction effects between pairs of loci, the large number of possible tests presents a significant computational challenge. The number of computations is further multiplied in the study of gene expression quantitative trait mapping, in which tests are performed for thousands of gene phenotypes simultaneously. Results: We present FastEpistasis, an efficient parallel solution extending the PLINK epistasis module, designed to test for epistasis effects when analyzing continuous phenotypes. Our results show that the algorithm scales with the number of processors and offers a reduction in computation time when several phenotypes are analyzed simultaneously. FastEpistasis is capable of testing the association of a continuous trait with all single nucleotide polymorphism ( SNP) pairs from 500 000 SNPs, totaling 125 billion tests, in a population of 5000 individuals in 29, 4 or 0.5 days using 8, 64 or 512 processors.

Genome-wide association study identifies new HLA class II haplotypes strongly protective against narcolepsy.

Narcolepsy is a rare sleep disorder with the strongest human leukocyte antigen (HLA) association ever reported. Since the associated HLA-DRB1*1501-DQB1*0602 haplotype is common in the general population (15-25%), it has been suggested that it is almost necessary but not sufficient for developing narcolepsy. To further define the genetic basis of narcolepsy risk, we performed a genome-wide association study (GWAS) in 562 European individuals with narcolepsy (cases) and 702 ethnically matched controls, with independent replication in 370 cases and 495 controls, all heterozygous for DRB1*1501-DQB1*0602. We found association with a protective variant near HLA-DQA2 (rs2858884; P < 3 x 10(-8)). Further analysis revealed that rs2858884 is strongly linked to DRB1*03-DQB1*02 (P < 4 x 10(-43)) and DRB1*1301-DQB1*0603 (P < 3 x 10(-7)). Cases almost never carried a trans DRB1*1301-DQB1*0603 haplotype (odds ratio = 0.02; P < 6 x 10(-14)). This unexpected protective HLA haplotype suggests a virtually causal involvement of the HLA region in narcolepsy susceptibility.