Characterization of mouse profilaggrin: Evidence for nuclear engulfment and translocation of the profilaggrin B-domain during epidermal differentiation

Filaggrin is a keratin filament associated protein that is expressed in granular layer keratinocytes and derived by sequential proteolysis from a polyprotein precursor termed profilaggrin. Depending on the species, each profilaggrin molecule contains between 10 and 20 filaggrin subunits organized as tandem repeats with a calcium-binding domain at the N-terminal end. We now report the characterization of the complete mouse gene. The structural organization of the mouse gene is identical to the human profilaggrin gene and consists of three exons with a 4 kb intron within the 5' noncoding region and a 1.7 kb intron separating the sequences encoding the calcium-binding EF-hand motifs. A processed pseudogene was found embedded within the second intron. The third and largest exon encodes the second EF-hand, a basic domain (designated the B-domain) followed by 12 filaggrin repeats and a unique C-terminal tail domain. A polyclonal anti-body raised against the conceptually translated sequence of the B-domain specifically stained keratohyalin granules and colocalized with a filaggrin antibody in granular layer cells. In upper granular layer cells, B-domain containing keratohyalin granules were in close apposition to the nucleus and, in some cells, appeared to be completely engulfed by the nucleus. In transition layer cells, B-domain staining was evident in the nucleus whereas filaggrin staining remained cytoplasmic. Nuclear staining of the B-domain was also observed in primary mouse keratinocytes induced to differentiate. This study has also revealed significant sequence homology between the mouse and human promoter sequences and in the calcium-binding domain but the remainder of the protein-coding region shows substantial divergence.

Evidence for a QTL on chromosome 19 influencing LDL cholesterol levels in the general population

The genetic basis of cardiovascular disease (CVD) with its complex etiology is still largely elusive. Plasma levels of lipids and apolipoproteins are among the major quantitative risk factors for CVD and are well-established intermediate traits that may be more accessible to genetic dissection than clinical CVD end points. Chromosome 19 harbors multiple genes that have been suggested to play a role in lipid metabolism and previous studies indicated the presence of a quantitative trait locus (QTL) for cholesterol levels in genetic isolates. To establish the relevance of genetic variation at chromosome 19 for plasma levels of lipids and apolipoproteins in the general, out-bred Caucasian population, we performed a linkage study in four independent samples, including adolescent Dutch twins and adult Dutch, Swedish and Australian twins totaling 493 dizygotic twin pairs. The average spacing of short-tandem-repeat markers was 6 - 8 cM. In the three adult twin samples, we found consistent evidence for linkage of chromosome 19 with LDL cholesterol levels ( maximum LOD scores of 4.5, 1.7 and 2.1 in the Dutch, Swedish and Australian sample, respectively); no indication for linkage was observed in the adolescent Dutch twin sample. The QTL effects in the three adult samples were not significantly different and a simultaneous analysis of the samples increased the maximum LOD score to 5.7 at 60 cM pter. Bivariate analyses indicated that the putative LDL-C QTL also contributed to the variance in ApoB levels, consistent with the high genetic correlation between these phenotypes. Our study provides strong evidence for the presence of a QTL on chromosome 19 with a major effect on LDL-C plasma levels in outbred Caucasian populations.

Evaluating plant breeding strategies by simulating gene action and dryland environment effects

Functional genomics is the systematic study of genome-wide effects of gene expression on organism growth and development with the ultimate aim of understanding how networks of genes influence traits. Here, we use a dynamic biophysical cropping systems model (APSIM-Sorg) to generate a state space of genotype performance based on 15 genes controlling four adaptive traits and then search this spice using a quantitative genetics model of a plant breeding program (QU-GENE) to simulate recurrent selection. Complex epistatic and gene X environment effects were generated for yield even though gene action at the trait level had been defined as simple additive effects. Given alternative breeding strategies that restricted either the cultivar maturity type or the drought environment type, the positive (+) alleles for 15 genes associated with the four adaptive traits were accumulated at different rates over cycles of selection. While early maturing genotypes were favored in the Severe-Terminal drought environment type, late genotypes were favored in the Mild-Terminal and Midseason drought environment types. In the Severe-Terminal environment, there was an interaction of the stay-green (SG) trait with other traits: Selection for + alleles of the SG genes was delayed until + alleles for genes associated with the transpiration efficiency and osmotic adjustment traits had been fixed. Given limitations in our current understanding of trait interaction and genetic control, the results are not conclusive. However, they demonstrate how the per se complexity of gene X gene X environment interactions will challenge the application of genomics and marker-assisted selection in crop improvement for dryland adaptation.

Epigenetic and transcriptional signatures of stable versus plastic differentiation of proinflammatory gd T cell subsets

Two distinct subsets of γδ T cells that produce interleukin 17 (IL-17) (CD27(-) γδ T cells) or interferon-γ (IFN-γ) (CD27(+) γδ T cells) develop in the mouse thymus, but the molecular determinants of their functional potential in the periphery remain unknown. Here we conducted a genome-wide characterization of the methylation patterns of histone H3, along with analysis of mRNA encoding transcription factors, to identify the regulatory networks of peripheral IFN-γ-producing or IL-17-producing γδ T cell subsets in vivo. We found that CD27(+) γδ T cells were committed to the expression of Ifng but not Il17, whereas CD27(-) γδ T cells displayed permissive chromatin configurations at loci encoding both cytokines and their regulatory transcription factors and differentiated into cells that produced both IL-17 and IFN-γ in a tumor microenvironment.

Positive Epistasis Drives the Acquisition of Multidrug Resistance

The evolution of multiple antibiotic resistance is an increasing global problem. Resistance mutations are known to impair fitness, and the evolution of resistance to multiple drugs depends both on their costs individually and on how they interact-epistasis. Information on the level of epistasis between antibiotic resistance mutations is of key importance to understanding epistasis amongst deleterious alleles, a key theoretical question, and to improving public health measures. Here we show that in an antibiotic-free environment the cost of multiple resistance is smaller than expected, a signature of pervasive positive epistasis among alleles that confer resistance to antibiotics. Competition assays reveal that the cost of resistance to a given antibiotic is dependent on the presence of resistance alleles for other antibiotics. Surprisingly we find that a significant fraction of resistant mutations can be beneficial in certain resistant genetic backgrounds, that some double resistances entail no measurable cost, and that some allelic combinations are hotspots for rapid compensation. These results provide additional insight as to why multi-resistant bacteria are so prevalent and reveal an extra layer of complexity on epistatic patterns previously unrecognized, since it is hidden in genome-wide studies of genetic interactions using gene knockouts.

The Photoreceptor Cell-Specific Nuclear Receptor Gene (PNR ) Accounts for Retinitis Pigmentosa in the Crypto-Jews from Portugal (Marranos), Survivors from the Spanish Inquisition

The last Crypto-Jews (Marranos) are the survivors of Spanish Jews who were persecuted in the late fifteenth century, escaped to Portugal and were forced to convert to save their lives. Isolated groups still exist in mountainous areas such as Belmonte in the Beira-Baixa province of Portugal. We report here the genetic study of a highly consanguineous endogamic population of Crypto-Jews of Belmonte affected with autosomal recessive retinitis pigmentosa (RP). A genome-wide search for homozygosity allowed us to localize the disease gene to chromosome 15q22-q24 (Zmax=2.95 at θ=0 at the D15S131 locus). Interestingly, the photoreceptor cell-specific nuclear receptor (PNR) gene, the expression of which is restricted to the outer nuclear layer of retinal photoreceptor cells, was found to map to the YAC contig encompassing the disease locus. A search for mutations allowed us to ascribe the RP of Crypto-Jews of Belmonte to a homozygous missense mutation in the PNR gene. Preliminary haplotype studies support the view that this mutation is relatively ancient but probably occurred after the population settled in Belmonte.

Understanding the relationship between central metabolism and virulence in the human pathogen Streptococcus pneumoniae

Dissertation presented to obtain the Ph.D. degree in Biochemistry

Gene Expression Profiling and Association Studies Implicate the Neuregulin Signaling Pathway in Behçet's Disease Susceptibility

Behçet's disease (BD) is a complex disease with genetic and environmental risk factors implicated in its etiology; however, its pathophysiology is poorly understood. To decipher BD's genetic underpinnings, we combined gene expression profiling with pathway analysis and association studies. We compared the gene expression profiles in peripheral blood mononuclear cells (PBMCs) of 15 patients and 14 matched controls using Affymetrix microarrays and found that the neuregulin signaling pathway was over-represented among the differentially expressed genes. The Epiregulin (EREG), Amphiregulin (AREG), and Neuregulin-1 (NRG1) genes of this pathway stand out as they are also among the top differentially expressed genes. Twelve haplotype tagging SNPs at the EREG-AREG locus and 15 SNPs in NRG1 found associated in at least one published BD genome-wide association study were tested for association with BD in a dataset of 976 Iranian patients and 839 controls. We found a novel association with BD for the rs6845297 SNP located downstream of EREG, and replicated three associations at NRG1 (rs4489285, rs383632, and rs1462891). Multifactor dimensionality reduction analysis indicated the existence of epistatic interactions between EREG and NRG1 variants. EREG-AREG and NRG1, which are members of the epidermal growth factor (EGF) family, seem to modulate BD susceptibility through main effects and gene–gene interactions. These association findings support a role for the EGF/ErbB signaling pathway inBD pathogenesis that warrants further investigation and highlight the importance of combining genetic and genomic approaches to dissect the genetic architecture of complex diseases.

Functional genomic analysis of heat stress in Vitis vinifera

Grapevine (Vitis vinifera) is one of most agro-economically important fruit crops worldwide, with a special relevance in Portugal where over 300 varieties are used for wine production. Due to global warming, temperature stress is currently a serious issue affecting crop production especially in temperate climates. Mobile genetic elements such as retrotransposons have been shown to be involved in environmental stress induced genetic and epigenetic modifications. In this study, sequences related to Grapevine Retrotransposon 1 (Gret1) were utilized to determine heat induced genomic and transcriptomic modifications in Touriga Nacional, a traditional Portuguese grapevine variety. For this purpose, growing canes were treated to 42 oC for four hours and leaf genomic DNA and RNA was utilized for various techniques to observe possible genomic alterations and variation in transcription levels of coding and non-coding sequences between non-treated plants and treated plants immediately after heat stress (HS-0 h) or after a 24 hour recovery period (HS-24 h). Heat stress was found to induce a significant decrease in Gret1 related sequences in HS-24 h leaves, indicating an effect of heat stress on genomic structure. In order to identify putative heat induced DNA modifications, genome wide approaches such as Amplified Fragment Length Polymorphism were utilized. This resulted in the identification of a polymorphic DNA fragment in HS-0 h and HS-24 h leaves whose sequence mapped to a genomic region flanking a house keeping gene (NADH) that is represented in multiple copies in the Vitis vinifera genome. Heat stress was also found to affect the transcript levels of various non-coding and gene coding sequences. Accordingly, quantitative real time PCR results established that Gret1 related sequences are up regulated immediately after heat stress whereas the level of transcript of genes involved in identification and repair of double strand breaks are significantly down regulated in HS-0 h plants. Taken together, the results of this work demonstrated heat stress affects both genomic integrity and transcription levels.

Urothelial bladder cancer progression: lessons learned from the bench

Urothelial bladder carcinoma (UBC) is an intricate malignancy with a variable natural history and clinical behavior. Despite developments in diagnosis/prognosis refinement and treatment modalities, the recurrence rate is high, and progression from non-muscle to muscle invasive UBC commonly leads to metastasis. Moreover, patients with muscle-invasive or extra-vesical disease often fail the standard chemotherapy treatment, and overall survival rates are poor. Thus, UBC remains a challenge in the oncology field, representing an ideal candidate for research on biomarkers that could identify patients at increased risk of recurrence, progression, and chemo-refractoriness. However, progress toward personalized medicine has been hampered by the unique genetic complexity of UBC. Recent genome-wide expression and sequencing studies have brought new insights into its molecular features, pathogenesis and clinical diversity, revealing a landscape where classical pathology is intersected by the novel and heterogeneous molecular groups. Hence, it seems plausible to postulate that only an integrated signature of prognostic/predictive biomarkers inherent in different cancer hallmarks will reach clinical validation. In this review, we have summarized ours and others' research into novel putative biomarkers of progression and chemoresistance that encompass several hallmarks of cancer: tumor neovascularization, invasion and metastasis, and energy metabolism reprogramming of the tumor microenvironment.

Extensive admixture and selective pressure across the Sahel Belt

Genome-wide studies of African populations have the potential to reveal powerful insights into the evolution of our species, as these diverse populations have been exposed to intense selective pressures imposed by infectious diseases, diet, and environmental factors. Within Africa, the Sahel Belt extensively overlaps the geographical center of several endemic infections such as malaria, trypanosomiasis, meningitis, and hemorrhagic fevers. We screened 2.5 million single nucleotide polymorphisms in 161 individuals from 13 Sahelian populations, which together with published data cover Western, Central, and Eastern Sahel, and include both nomadic and sedentary groups. We confirmed the role of this Belt as a main corridor for human migrations across the continent. Strong admixture was observed in both Central and Eastern Sahelian populations, with North Africans and Near Eastern/Arabians, respectively, but it was inexistent in Western Sahelian populations. Genome-wide local ancestry inference in admixed Sahelian populations revealed several candidate regions that were significantly enriched for non-autochthonous haplotypes, and many showed to be under positive selection. The DARC gene region in Arabs and Nubians was enriched for African ancestry, whereas the RAB3GAP1/LCT/MCM6 region in Oromo, the TAS2R gene family in Fulani, and the ALMS1/NAT8 in Turkana and Samburu were enriched for non-African ancestry. Signals of positive selection varied in terms of geographic amplitude. Some genomic regions were selected across the Belt, the most striking example being the malaria-related DARC gene. Others were Western-specific (oxytocin, calcium, and heart pathways), Eastern-specific (lipid pathways), or even population-restricted (TAS2R genes in Fulani, which may reflect sexual selection).

Type 2 diabetes-related variants influence the risk of developing multiple myeloma: results from the IMMEnSE consortium

Type 2 diabetes (T2D) has been suggested to be a risk factor for multiple myeloma (MM), but the relationship between the two traits is still not well understood. The aims of this study were to evaluate whether 58 genome-wide-association-studies (GWAS)-identified common variants for T2D influence the risk of developing MM and to determine whether predictive models built with these variants might help to predict the disease risk. We conducted a case–control study including 1420 MM patients and 1858 controls ascertained through the International Multiple Myeloma (IMMEnSE) consortium. Subjects carrying the KCNQ1rs2237892T allele or the CDKN2A-2Brs2383208G/G, IGF1rs35767T/T and MADDrs7944584T/T genotypes had a significantly increased risk of MM (odds ratio (OR)=1.32–2.13) whereas those carrying the KCNJ11rs5215C, KCNJ11rs5219T and THADArs7578597C alleles or the FTOrs8050136A/A and LTArs1041981C/C genotypes showed a significantly decreased risk of developing the disease (OR=0.76–0.85). Interestingly, a prediction model including those T2D-related variants associated with the risk of MM showed a significantly improved discriminatory ability to predict the disease when compared to a model without genetic information (area under the curve (AUC)=0.645 vs AUC=0.629; P=4.05×10-06). A gender-stratified analysis also revealed a significant gender effect modification for ADAM30rs2641348 and NOTCH2rs10923931 variants (Pinteraction=0.001 and 0.0004, respectively). Men carrying the ADAM30rs2641348C and NOTCH2rs10923931T alleles had a significantly decreased risk of MM whereas an opposite but not significant effect was observed in women (ORM=0.71 and ORM=0.66 vs ORW=1.22 and ORW=1.15, respectively). These results suggest that TD2-related variants may influence the risk of developing MM and their genotyping might help to improve MM risk prediction models.

Resolving the ancestry of Austronesian‑speaking populations

There are two very different interpretations of the prehistory of Island Southeast Asia (ISEA), with genetic evidence invoked in support of both. The "out-of-Taiwan" model proposes a major Late Holocene expansion of Neolithic Austronesian speakers from Taiwan. An alternative, proposing that Late Glacial/postglacial sea-level rises triggered largely autochthonous dispersals, accounts for some otherwise enigmatic genetic patterns, but fails to explain the Austronesian language dispersal. Combining mitochondrial DNA (mtDNA), Y-chromosome and genome-wide data, we performed the most comprehensive analysis of the region to date, obtaining highly consistent results across all three systems and allowing us to reconcile the models. We infer a primarily common ancestry for Taiwan/ISEA populations established before the Neolithic, but also detected clear signals of two minor Late Holocene migrations, probably representing Neolithic input from both Mainland Southeast Asia and South China, via Taiwan. This latter may therefore have mediated the Austronesian language dispersal, implying small-scale migration and language shift rather than large-scale expansion.

A genetic perspective on African prehistory

The various genetic systems (mitochondrial DNA, the Y-chromosome and the genome-wide autosomes) indicate that Africa is the most genetically diverse continent in the world and the most likely place of origin for anatomically modern humans. However, where in Africa modern humans arose and how the current genetic makeup within the continent was shaped is still open to debate. Here, we summarize the debate and focus especially on the maternally inherited mitochondrial DNA (mtDNA) and a recently revised chronology for the African mtDNA tree. We discuss the possible origin of modern humans in southern, eastern or Central Africa; the possibility of a migration from southern to eastern Africa more than 100 ka, carrying lineages within mtDNA haplogroup L0; the evidence for a climate-change-mediated population expansion in eastern Africa involving mtDNA haplogroup L3, leading to the “out-of-Africa” migration around 70–60 ka; the re-population of North Africa from the Near East around 40–30 ka suggested by mtDNA haplogroups U6 and M1; the evidence for population expansions and dispersals across the continent at the onset of the Holocene ; and the impact of the Bantu dispersals in Central, eastern and southern Africa within the last few millennia.

CLOCK-controlled polyphonic regulation of circadian rhythms through canonical and noncanonical E-boxes.

In mammalian circadian clockwork, the CLOCK-BMAL1 complex binds to DNA enhancers of target genes and drives circadian oscillation of transcription. Here we identified 7,978 CLOCK-binding sites in mouse liver by chromatin immunoprecipitation-sequencing (ChIP-Seq), and a newly developed bioinformatics method, motif centrality analysis of ChIP-Seq (MOCCS), revealed a genome-wide distribution of previously unappreciated noncanonical E-boxes targeted by CLOCK. In vitro promoter assays showed that CACGNG, CACGTT, and CATG(T/C)G are functional CLOCK-binding motifs. Furthermore, we extensively revealed rhythmically expressed genes by poly(A)-tailed RNA-Seq and identified 1,629 CLOCK target genes within 11,926 genes expressed in the liver. Our analysis also revealed rhythmically expressed genes that have no apparent CLOCK-binding site, indicating the importance of indirect transcriptional and posttranscriptional regulations. Indirect transcriptional regulation is represented by rhythmic expression of CLOCK-regulated transcription factors, such as Krüppel-like factors (KLFs). Indirect posttranscriptional regulation involves rhythmic microRNAs that were identified by small-RNA-Seq. Collectively, CLOCK-dependent direct transactivation through multiple E-boxes and indirect regulations polyphonically orchestrate dynamic circadian outputs.