Interleukin-16 is increased in obesity and alters adipogenesis and inflammation in vitro

Introduction

Obesity is a chronic condition associated with low-grade inflammation mainly due to immune cell infiltration of white adipose tissue (WAT). WAT is distributed into two main depots: subcutaneous WAT (sWAT) and visceral WAT (vWAT), each with different biochemical features and metabolic roles. Proinflammatory cytokines including interleukin (IL)-16 are secreted by both adipocytes and infiltrated immune cells to upregulate inflammation. IL-16 has been widely studied in the peripheral proinflammatory immune response; however, little is known about its role in adipocytes in the context of obesity.

Aim & Methods

We aimed to study the levels of IL-16 in WAT derived from sWAT and vWAT depots of humans with obesity and the role of this cytokine in palmitate-exposed 3T3-L1 adipocytes.

Results

The results demonstrated that IL-16 expression was higher in vWAT compared with sWAT in individuals with obesity. In addition, IL-16 serum levels were higher in patients with obesity compared with normal-weight individuals, increased at 6 months after bariatric surgery, and at 12 months after surgery decreased to levels similar to before the intervention. Our in vitro models showed that IL-16 could modulate markers of adipogenesis (Pref1), lipid metabolism (Plin1, Cd36, and Glut4), fibrosis (Hif1a, Col4a, Col6a, and Vegf), and inflammatory signaling (IL6) during adipogenesis and in mature adipocytes. In addition, lipid accumulation and glycerol release assays suggested lipolysis alteration.

Discussion

Our results suggest a potential role of IL-16 in adipogenesis, lipid and glucose homeostasis, fibrosis, and inflammation in an obesity context.

Claudin-1 as a novel target gene induced in obesity and associated to inflammation, fibrosis, and cell differentiation

OBJECTIVE

T lymphocytes from visceral and subcutaneous white adipose tissues (vWAT and sWAT, respectively) can have opposing roles in the systemic metabolic changes associated with obesity. However, few studies have focused on this subject. Claudin-1 (CLDN1) is a protein involved canonically in tight junctions and tissue paracellular permeability. We evaluated T-lymphocyte gene expression in vWAT and sWAT and in the whole adipose depots in human samples.

METHODS

A Clariom D-based transcriptomic analysis was performed on T lymphocytes magnetically separated from vWAT and sWAT from patients with obesity (Cohort 1; N = 11). Expression of candidate genes resulting from that analysis was determined in whole WAT from individuals with and without obesity (Cohort 2; patients with obesity: N = 13; patients without obesity: N = 14).

RESULTS

We observed transcriptional differences between T lymphocytes from sWAT compared with vWAT. Specifically, CLDN1 expression was found to be dramatically induced in vWAT T cells relative to those isolated from sWAT in patients with obesity. CLDN1 was also induced in obesity in vWAT and its expression correlates with genes involved in inflammation, fibrosis, and adipogenesis.

CONCLUSION

These results suggest that CLDN1 is a novel marker induced in obesity and differentially expressed in T lymphocytes infiltrated in human vWAT as compared with sWAT. This protein may have a crucial role in the crosstalk between T lymphocytes and other adipose tissue cells and may contribute to inflammation, fibrosis, and alter homeostasis and promote metabolic disease in obesity.

Seminal extracellular vesicle sncRNA sequencing reveals altered miRNA/isomiR profiles as sperm retrieval biomarkers for azoospermia

BACKGROUND

Non-invasive molecular biomarkers for classifying azoospermia by origin into either obstructive or non-obstructive/secretory azoospermia, as well as for inferring the spermatogenic reserve of the testis of non-obstructive/secretory azoospermia patients, are of great interest for testicular sperm retrieval outcome prediction for assisted reproduction. Prior analyses of semen small non-coding RNA expression in azoospermia have focused on microRNAs, but there has been a lack of attention on other regulatory small RNA species. In this regard, studying more in-depth expression changes of small non-coding RNA subtypes in small extracellular vesicles from semen of azoospermic individuals could be useful to select additional non-invasive biomarkers with diagnostic/prognostic purposes.

MATERIAL AND METHODS

A high-throughput small RNA profiling analysis to determine the expression pattern of seminal small extracellular vesicle microRNAs (analyzed at the isomiR level), PIWI-interacting RNAs, and transfer RNA-derived small RNAs in normozoospermic (n = 4) and azoospermic (obstructive azoospermia because of pathological occurring obstruction in the genital tract, n = 4; secretory azoospermic individuals with positive testicular sperm extraction value, n = 5; secretory azoospermic individuals with negative testicular sperm extraction value, n = 4) individuals was carried out. Reverse transcriptase-quantitative real-time polymerase chain reaction validation analysis of selected microRNAs was additionally performed in a larger number of individuals.

RESULTS AND DISCUSSION

Clinically relevant quantitative changes in the small non-coding RNA levels contained in semen small extracellular vesicles can be used as biomarkers for the origin of azoospermia and for predicting the presence of residual spermatogenesis. In this regard, canonical isoform microRNAs (n = 185) but also other isomiR variants (n = 238) stand out in terms of numbers and fold-change differences in expression, underlining the need to consider isomiRs when investigating microRNA-based regulation. Conversely, although transfer RNA-derived small RNAs are shown in our study to represent a high proportion of small non-coding RNA sequences in seminal small extracellular vesicle samples, they are not able to discriminate the origin of azoospermia. PIWI-interacting RNA cluster profiles and individual PIWI-interacting RNAs with significant differential expression were also not able to discriminate. Our study demonstrated that expression values of individual and/or combined canonical isoform microRNAs (miR-10a-5p, miR-146a-5p, miR-31-5p, miR-181b-5p; area under the receiver operating characteristic curve >0.8) in small extracellular vesicles provide considerable clinical value in identifying samples with a high likelihood of sperm retrieval while discriminating azoospermia by origin. Although no individual microRNA showed sufficient discriminating power on its own to identify severe spermatogenic disorders with focal spermatogenesis, multivariate microRNA models in semen small extracellular vesicles have the potential to identify those individuals with residual spermatogenesis. Availability and adoption of such non-invasive molecular biomarkers would represent a great improvement in reproductive treatment decision protocols for azoospermia in clinical practice.

3'IsomiR Species Composition Affects Reliable Quantification of miRNA/isomiR Variants by Poly(A) RT-qPCR: Impact on Small RNA-Seq Profiling Validation

Small RNA-sequencing (small RNA-seq) has revealed the presence of small RNA-naturally occurring variants such as microRNA (miRNA) isoforms or isomiRs. Due to their small size and the sequence similarity among miRNA isoforms, their validation by RT-qPCR is challenging. We previously identified two miR-31-5p isomiRs-the canonical and a 3'isomiR variant (3' G addition)-which were differentially expressed between individuals with azoospermia of different origin. Here, we sought to determine the discriminatory capacity between these two closely-related miRNA isoforms of three alternative poly(A) based-RT-qPCR strategies in both synthetic and real biological context. We found that these poly(A) RT-qPCR strategies exhibit a significant cross-reactivity between these miR-31-5p isomiRs which differ by a single nucleotide, compromising the reliable quantification of individual miRNA isoforms. Fortunately, in the biological context, given that the two miRNA variants show changes in the same direction, RT-qPCR results were consistent with the findings of small RNA-seq study. We suggest that miRNA selection for RT-qPCR validation should be performed with care, prioritizing those canonical miRNAs that, in small RNA-seq, show parallel/homogeneous expression behavior with their most prevalent isomiRs, to avoid confounding RT-qPCR-based results. This is suggested as the current best strategy for robust biomarker selection to develop clinically useful tests.

Monitoring SARS-CoV-2 variant transitions using differences in diagnostic cycle threshold values of target genes

Monitoring the emergence of new SARS-CoV-2 variants is important to detect potential risks of increased transmission or disease severity. We investigated the identification of SARS-CoV-2 variants from real-time reverse transcriptase polymerase chain reaction (RT-PCR) routine diagnostics data. Cycle threshold (Ct) values of positive samples were collected from April 2021 to January 2022 in the Northern Metropolitan Area of Barcelona (n = 15,254). Viral lineage identification from whole genome sequencing (WGS) was available for 4618 (30.3%) of these samples. Pairwise differences in the Ct values between gene targets (ΔCt) were analyzed for variants of concern or interest circulating in our area. A specific delay in the Ct of the N-gene compared to the RdRp-gene (ΔCt_NR) was observed for Alpha, Delta, Eta and Omicron. Temporal differences in ΔCt_NR correlated with the dynamics of viral replacement of Alpha by Delta and of Delta by Omicron according to WGS results. Using ΔCt_NR, prediction of new variants of concern at early stages of circulation was achieved with high sensitivity and specificity (91.1% and 97.8% for Delta; 98.5% and 90.8% for Omicron). Thus, tracking population-wide trends in ΔCt values obtained from routine diagnostics testing in combination with WGS could be useful for real-time management and response to local epidemics.

Detailed stratified GWAS analysis for severe COVID-19 in four European populations

Given the highly variable clinical phenotype of Coronavirus disease 2019 (COVID-19), a deeper analysis of the host genetic contribution to severe COVID-19 is important to improve our understanding of underlying disease mechanisms. Here, we describe an extended genome-wide association meta-analysis of a well-characterized cohort of 3255 COVID-19 patients with respiratory failure and 12 488 population controls from Italy, Spain, Norway and Germany/Austria, including stratified analyses based on age, sex and disease severity, as well as targeted analyses of chromosome Y haplotypes, the human leukocyte antigen region and the SARS-CoV-2 peptidome. By inversion imputation, we traced a reported association at 17q21.31 to a ~0.9-Mb inversion polymorphism that creates two highly differentiated haplotypes and characterized the potential effects of the inversion in detail. Our data, together with the 5th release of summary statistics from the COVID-19 Host Genetics Initiative including non-Caucasian individuals, also identified a new locus at 19q13.33, including NAPSA, a gene which is expressed primarily in alveolar cells responsible for gas exchange in the lung.

A Polygenic Risk Score Based on a Cardioembolic Stroke Multitrait Analysis Improves a Clinical Prediction Model for This Stroke Subtype

Background

Occult atrial fibrillation (AF) is one of the major causes of embolic stroke of undetermined source (ESUS). Knowing the underlying etiology of an ESUS will reduce stroke recurrence and/or unnecessary use of anticoagulants. Understanding cardioembolic strokes (CES), whose main cause is AF, will provide tools to select patients who would benefit from anticoagulants among those with ESUS or AF. We aimed to discover novel loci associated with CES and create a polygenetic risk score (PRS) for a more efficient CES risk stratification.

Methods

Multitrait analysis of GWAS (MTAG) was performed with MEGASTROKE-CES cohort (n = 362,661) and AF cohort (n = 1,030,836). We considered significant variants and replicated those variants with MTAG p-value < 5 × 10⁻⁸ influencing both traits (GWAS-pairwise) with a p-value < 0.05 in the original GWAS and in an independent cohort (n = 9,105). The PRS was created with PRSice-2 and evaluated in the independent cohort.

Results

We found and replicated eleven loci associated with CES. Eight were novel loci. Seven of them had been previously associated with AF, namely, CAV1, ESR2, GORAB, IGF1R, NEURL1, WIPF1, and ZEB2. KIAA1755 locus had never been associated with CES/AF, leading its index variant to a missense change (R1045W). The PRS generated has been significantly associated with CES improving discrimination and patient reclassification of a model with age, sex, and hypertension.

Conclusion

The loci found significantly associated with CES in the MTAG, together with the creation of a PRS that improves the predictive clinical models of CES, might help guide future clinical trials of anticoagulant therapy in patients with ESUS or AF.

GCAT|Panel, a comprehensive structural variant haplotype map of the Iberian population from high-coverage whole-genome sequencing

The combined analysis of haplotype panels with phenotype clinical cohorts is a common approach to explore the genetic architecture of human diseases. However, genetic studies are mainly based on single nucleotide variants (SNVs) and small insertions and deletions (indels). Here, we contribute to fill this gap by generating a dense haplotype map focused on the identification, characterization, and phasing of structural variants (SVs). By integrating multiple variant identification methods and Logistic Regression Models (LRMs), we present a catalogue of 35 431 441 variants, including 89 178 SVs (≥50 bp), 30 325 064 SNVs and 5 017 199 indels, across 785 Illumina high coverage (30x) whole-genomes from the Iberian GCAT Cohort, containing a median of 3.52M SNVs, 606 336 indels and 6393 SVs per individual. The haplotype panel is able to impute up to 14 360 728 SNVs/indels and 23 179 SVs, showing a 2.7-fold increase for SVs compared with available genetic variation panels. The value of this panel for SVs analysis is shown through an imputed rare Alu element located in a new locus associated with Mononeuritis of lower limb, a rare neuromuscular disease. This study represents the first deep characterization of genetic variation within the Iberian population and the first operational haplotype panel to systematically include the SVs into genome-wide genetic studies.

Characterization of Acute HCV Infection and Transmission Networks in People Who Currently Inject Drugs in Catalonia: Usefulness of Dried Blood Spots

BACKGROUND AND AIMS

Accurate identification of recent HCV infections is critical for tracing the extent and mechanisms of ongoing transmission. We aimed to validate dried blood spot (DBS) samples for the assessment of Hepatitis C virus (HCV) genetic diversity and to determine epidemiological parameters including incidence, determinants of acute infection, and phylogenetic clustering in people who inject drugs (PWID).

APPROACH AND RESULTS

HCV nonstructural protein 5B next-generation sequencing was performed from plasma and/or DBS in 220 viremic PWID from the HepCdetect II study. No significant differences were found in consensus sequences or Shannon entropy (SE) intrahost diversity estimate between paired plasma/DBS specimens. SE values were used to identify acute infections with 93.3% sensitivity (95% CI, 0.81-1.06) and 95.0% specificity (95% CI, 0.88-1.02) in a set of well-defined controls. An acute HCV infection (either primary infection or reinfection) was detected in 13.5% of viremic participants and was associated with age ≤30 years (OR, 8.09), injecting less than daily (OR, 4.35), ≤5 years of injected drug use (OR, 3.43), sharing cocaine snorting straws (OR, 2.89), and being unaware of their HCV status (OR, 3.62). Annualized HCV incidence was estimated between 31 and 59/100 person-years. On phylogenetic analysis, 46.8% of viremic cases were part of a transmission pair or cluster; age ≤30 years (OR, 6.16), acute infection (OR, 5.73), and infection with subtype 1a (OR, 4.78) were independently associated with this condition.

CONCLUSIONS

The results obtained from plasma and DBS characterize PWID with acute infection and those involved in ongoing HCV transmission and allow estimating incidence from cross-sectional data. This information is critical for the design and assessment of targeted harm reduction programs and test-and-treat interventions and to facilitate monitoring of HCV elimination in this key population.

Paired-end small RNA sequencing reveals a possible overestimation in the isomiR sequence repertoire previously reported from conventional single read data analysis

BACKGROUND

Next generation sequencing has allowed the discovery of miRNA isoforms, termed isomiRs. Some isomiRs are derived from imprecise processing of pre-miRNA precursors, leading to length variants. Additional variability is introduced by non-templated addition of bases at the ends or editing of internal bases, resulting in base differences relative to the template DNA sequence. We hypothesized that some component of the isomiR variation reported so far could be due to systematic technical noise and not real.

RESULTS

We have developed the XICRA pipeline to analyze small RNA sequencing data at the isomiR level. We exploited its ability to use single or merged reads to compare isomiR results derived from paired-end (PE) reads with those from single reads (SR) to address whether detectable sequence differences relative to canonical miRNAs found in isomiRs are true biological variations or the result of errors in sequencing. We have detected non-negligible systematic differences between SR and PE data which primarily affect putative internally edited isomiRs, and at a much smaller frequency terminal length changing isomiRs. This is relevant for the identification of true isomiRs in small RNA sequencing datasets.

CONCLUSIONS

We conclude that potential artifacts derived from sequencing errors and/or data processing could result in an overestimation of abundance and diversity of miRNA isoforms. Efforts in annotating the isomiRnome should take this into account.

Loss of HDAC11 accelerates skeletal muscle regeneration in mice

Histone deacetylase 11 (HDAC11) is the latest identified member of the histone deacetylase family of enzymes. It is highly expressed in brain, heart, testis, kidney, and skeletal muscle, although its role in these tissues is poorly understood. Here, we investigate for the first time the consequences of HDAC11 genetic impairment on skeletal muscle regeneration, a process principally dependent on its resident stem cells (satellite cells) in coordination with infiltrating immune cells and stromal cells. Our results show that HDAC11 is dispensable for adult muscle growth and establishment of the satellite cell population, while HDAC11 deficiency advances the regeneration process in response to muscle injury. This effect is not caused by differences in satellite cell activation or proliferation upon injury, but rather by an enhanced capacity of satellite cells to differentiate at early regeneration stages in the absence of HDAC11. Infiltrating HDAC11-deficient macrophages could also contribute to this accelerated muscle regenerative process by prematurely producing high levels of IL-10, a cytokine known to promote myoblast differentiation. Altogether, our results show that HDAC11 depletion advances skeletal muscle regeneration and this finding may have potential implications for designing new strategies for muscle pathologies coursing with chronic damage. DATABASE: Data were deposited in NCBI's Gene Expression Omnibus accessible through GEO Series accession number GSE147423.

A likelihood ratio approach for identifying three-quarter siblings in genetic databases

The detection of family relationships in genetic databases is of interest in various scientific disciplines such as genetic epidemiology, population and conservation genetics, forensic science, and genealogical research. Nowadays, screening genetic databases for related individuals forms an important aspect of standard quality control procedures. Relatedness research is usually based on an allele sharing analysis of identity by state (IBS) or identity by descent (IBD) alleles. Existing IBS/IBD methods mainly aim to identify first-degree relationships (parent-offspring or full siblings) and second degree (half-siblings, avuncular, or grandparent-grandchild) pairs. Little attention has been paid to the detection of in-between first and second-degree relationships such as three-quarter siblings (3/4S) who share fewer alleles than first-degree relationships but more alleles than second-degree relationships. With the progressively increasing sample sizes used in genetic research, it becomes more likely that such relationships are present in the database under study. In this paper, we extend existing likelihood ratio (LR) methodology to accurately infer the existence of 3/4S, distinguishing them from full siblings and second-degree relatives. We use bootstrap confidence intervals to express uncertainty in the LRs. Our proposal accounts for linkage disequilibrium (LD) by using marker pruning, and we validate our methodology with a pedigree-based simulation study accounting for both LD and recombination. An empirical genome-wide array data set from the GCAT Genomes for Life cohort project is used to illustrate the method.

PSA Kinetics as Prognostic Markers of Overall Survival in Patients with Metastatic Castration-Resistant Prostate Cancer Treated with Abiraterone Acetate

Background

Abiraterone acetate (AA) is widely used in the treatment of patients with metastatic castration-resistant prostate cancer (mCRPC). However, a significant percentage of patients will still progress, highlighting the need to identify patients more likely to benefit from AA. Parameters linked to prostate-specific antigen (PSA) kinetics are promising prognostic markers. We have examined clinical and PSA-related factors potentially associated with overall survival (OS) in patients treated with AA.

Methods

Between 2011 and 2014, 104 patients with mCRPC treated with AA after progression to docetaxel at centers of the Catalan Institute of Oncology were included in this retrospective study. Patients were assessed monthly. Baseline characteristics and variables related to PSA kinetics were included in univariate and multivariate analyses of OS.

Results

Median OS was 16.4 months (range 12.4-20.6) for all patients. The univariate analysis identified the following baseline characteristics as significantly associated with OS: ECOG PS, location of metastases, time between starting androgen deprivation therapy and starting AA, time between stopping docetaxel treatment and starting AA, neutrophil-lymphocyte ratio (NLR), alkaline phosphatase levels, and PSA levels. Factors related to PSA kinetics associated with longer OS were PSA response >50%, early PSA response (>30% decline at four weeks), PSA decline >50% at week 12, PSA nadir <2.4ng/mL, time to PSA nadir >140 days, the combination of PSA nadir and time to PSA nadir, and low end-of-treatment PSA levels. The multivariate analysis identified ECOG PS (HR 37.46; p<0.001), NLR (HR 3.7; p<0.001), early PSA response (HR 1.22; p=0.002), and time to PSA nadir (HR 0.39; p=0.002) as independent prognostic markers.

Conclusion

Our results indicate an association between PSA kinetics, especially early PSA response, and outcome to AA after progression to docetaxel. Taken together with other factors, lack of an early PSA response could identify patients who are unlikely to benefit from AA and who could be closely monitored with a view to offering alternative therapies.

Genomewide Association Study of Severe Covid-19 with Respiratory Failure

BACKGROUND

There is considerable variation in disease behavior among patients infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes coronavirus disease 2019 (Covid-19). Genomewide association analysis may allow for the identification of potential genetic factors involved in the development of Covid-19.

METHODS

We conducted a genomewide association study involving 1980 patients with Covid-19 and severe disease (defined as respiratory failure) at seven hospitals in the Italian and Spanish epicenters of the SARS-CoV-2 pandemic in Europe. After quality control and the exclusion of population outliers, 835 patients and 1255 control participants from Italy and 775 patients and 950 control participants from Spain were included in the final analysis. In total, we analyzed 8,582,968 single-nucleotide polymorphisms and conducted a meta-analysis of the two case-control panels.

RESULTS

We detected cross-replicating associations with rs11385942 at locus 3p21.31 and with rs657152 at locus 9q34.2, which were significant at the genomewide level (P<5×10-8) in the meta-analysis of the two case-control panels (odds ratio, 1.77; 95% confidence interval [CI], 1.48 to 2.11; P = 1.15×10-10; and odds ratio, 1.32; 95% CI, 1.20 to 1.47; P = 4.95×10-8, respectively). At locus 3p21.31, the association signal spanned the genes SLC6A20, LZTFL1, CCR9, FYCO1, CXCR6 and XCR1. The association signal at locus 9q34.2 coincided with the ABO blood group locus; in this cohort, a blood-group-specific analysis showed a higher risk in blood group A than in other blood groups (odds ratio, 1.45; 95% CI, 1.20 to 1.75; P = 1.48×10-4) and a protective effect in blood group O as compared with other blood groups (odds ratio, 0.65; 95% CI, 0.53 to 0.79; P = 1.06×10-5).

CONCLUSIONS

We identified a 3p21.31 gene cluster as a genetic susceptibility locus in patients with Covid-19 with respiratory failure and confirmed a potential involvement of the ABO blood-group system. (Funded by Stein Erik Hagen and others.).

Morphological and Transcriptional Changes in Human Bone Marrow During Natural Plasmodium vivax Malaria Infections

Background

The presence of Plasmodium vivax malaria parasites in the human bone marrow (BM) is still controversial. However, recent data from a clinical case and experimental infections in splenectomized nonhuman primates unequivocally demonstrated the presence of parasites in this tissue.

Methods

In the current study, we analyzed BM aspirates of 7 patients during the acute attack and 42 days after drug treatment. RNA extracted from CD71⁺ cell suspensions was used for sequencing and transcriptomic analysis.

Results

We demonstrated the presence of parasites in all patients during acute infections. To provide further insights, we purified CD71⁺ BM cells and demonstrated dyserythropoiesis and inefficient erythropoiesis in all patients. In addition, RNA sequencing from 3 patients showed that genes related to erythroid maturation were down-regulated during acute infections, whereas immune response genes were up-regulated.

Conclusions

This study thus shows that during P. vivax infections, parasites are always present in the BM and that such infections induced dyserythropoiesis and ineffective erythropoiesis. Moreover, infections induce transcriptional changes associated with such altered erythropoietic response, thus highlighting the importance of this hidden niche during natural infections.

Multitrait genome association analysis identifies new susceptibility genes for human anthropometric variation in the GCAT cohort

Background

Heritability estimates have revealed an important contribution of SNP variants for most common traits; however, SNP analysis by single-trait genome-wide association studies (GWAS) has failed to uncover their impact. In this study, we applied a multitrait GWAS approach to discover additional factor of the missing heritability of human anthropometric variation.

Methods

We analysed 205 traits, including diseases identified at baseline in the GCAT cohort (Genomes For Life- Cohort study of the Genomes of Catalonia) (n=4988), a Mediterranean adult population-based cohort study from the south of Europe. We estimated SNP heritability contribution and single-trait GWAS for all traits from 15 million SNP variants. Then, we applied a multitrait-related approach to study genome-wide association to anthropometric measures in a two-stage meta-analysis with the UK Biobank cohort (n=336 107).

Results

Heritability estimates (eg, skin colour, alcohol consumption, smoking habit, body mass index, educational level or height) revealed an important contribution of SNP variants, ranging from 18% to 77%. Single-trait analysis identified 1785 SNPs with genome-wide significance threshold. From these, several previously reported single-trait hits were confirmed in our sample with LINC01432 (p=1.9×10⁻⁹) variants associated with male baldness, LDLR variants with hyperlipidaemia (ICD-9:272) (p=9.4×10⁻¹⁰) and variants in IRF4 (p=2.8×10⁻⁵⁷), SLC45A2 (p=2.2×10⁻¹³⁰), HERC2 (p=2.8×10⁻¹⁷⁶), OCA2 (p=2.4×10⁻¹²¹) and MC1R (p=7.7×10⁻²²) associated with hair, eye and skin colour, freckling, tanning capacity and sun burning sensitivity and the Fitzpatrick phototype score, all highly correlated cross-phenotypes. Multitrait meta-analysis of anthropometric variation validated 27 loci in a two-stage meta-analysis with a large British ancestry cohort, six of which are newly reported here (p value threshold <5×10⁻⁹) at ZRANB2-AS2, PIK3R1, EPHA7, MAD1L1, CACUL1 and MAP3K9.

Conclusion

Considering multiple-related genetic phenotypes improve associated genome signal detection. These results indicate the potential value of data-driven multivariate phenotyping for genetic studies in large population-based cohorts to contribute to knowledge of complex traits.

GCAT|Genomes for life: a prospective cohort study of the genomes of Catalonia

PURPOSE

The prevalence of chronic non-communicable diseases (NCDs) is increasing worldwide. NCDs are the leading cause of both morbidity and mortality, and it is estimated that by 2030, they will be responsible for 80% of deaths across the world. The Genomes for Life (GCAT) project is a long-term prospective cohort study that was designed to integrate and assess the role of epidemiological, genomic and epigenomic factors in the development of major chronic diseases in Catalonia, a north-east region of Spain.

PARTICIPANTS

At the end of 2017, the GCAT Study will have recruited 20 000 participants aged 40-65 years. Participants who agreed to take part in the study completed a self-administered computer-driven questionnaire, and underwent blood pressure, cardiac frequency and anthropometry measurements. For each participant, blood plasma, blood serum and white blood cells are collected at baseline. The GCAT Study has access to the electronic health records of the Catalan Public Healthcare System. Participants will be followed biannually at least 20 years after recruitment.

FINDINGS TO DATE

Among all GCAT participants, 59.2% are women and 83.3% of the cohort identified themselves as Caucasian/white. More than half of the participants have higher education levels, 72.2% are current workers and 42.1% are classified as overweight (body mass index ≥25 and <30 kg/m²). We have genotyped 5459 participants, of which 5000 have metabolome data. Further, the whole genome of 808 participants will be sequenced by the end of 2017.

FUTURE PLANS

The first follow-up study started in December 2017 and will end by March 2018. Residences of all subjects will be geocoded during the following year. Several genomic analyses are ongoing, and metabolomic and genomic integrations will be performed to identify underlying genetic variants, as well as environmental factors that influence metabolites.

Genomic profiling in advanced stage non-small-cell lung cancer patients with platinum-based chemotherapy identifies germline variants with prognostic value in SMYD2

OBJECTIVE

The aim of the study was to investigate the relationship between germline variations as a prognosis biomarker in patients with advanced Non-Small-Cell-Lung-Cancer (NSCLC) subjected to first-line platinum-based treatment.

MATERIALS AND METHODS

We carried out a two-stage genome-wide-association study in non-small-cell lung cancer patients with platinum-based chemotherapy in an exploratory sample of 181 NSCLC patients from Caucasian origin, followed by a validation on 356 NSCLC patients from the same ancestry (Valencia, Spain).

RESULTS

We identified germline variants in SMYD2 as a prognostic factor for survival in patients with advanced NSCLC receiving chemotherapy. SMYD2 alleles are associated to a decreased overall survival and with a reduced Time to Progression. In addition, enrichment pathway analysis identified 361 variants in 40 genes to be involved in poorer outcome in advanced-stage NSCLC patients.

CONCLUSION

Germline SMYD2 alleles are associated with bad clinical outcome of first-line platinum-based treatment in advanced NSCLC patients. This result supports the role of SMYD2 in the carcinogenic process, and might be used as prognostic signature directing patient stratification and the choice of therapy.

MICROABSTRACT

A two-Stage Genome wide association study in Caucasian population reveals germline genetic variation in SMYD2 associated to progression disease in first-line platinum-based treatment in advanced NSCLC patients. SMYD2 profiling might have prognostic / predictive value directing choice of therapy and enlighten current knowledge on pathways involved in human carcinogenesis as well in resistance to chemotherapy.

Gene expression changes in blood RNA after swimming in a chlorinated pool

Exposure to disinfection by-products (DBP) such as trihalomethanes (THM) in swimming pools has been linked to adverse health effects in humans, but their biological mechanisms are unclear. We evaluated short-term changes in blood gene expression of adult recreational swimmers after swimming in a chlorinated pool. Volunteers swam 40min in an indoor chlorinated pool. Blood samples were drawn and four THM (chloroform, bromodichloromethane, dibromochloromethane and bromoform) were measured in exhaled breath before and after swimming. Intensity of physical activity was measured as metabolic equivalents (METs). Gene expression in whole blood mRNA was evaluated using IlluminaHumanHT-12v3 Expression-BeadChip. Linear mixed models were used to evaluate the relationship between gene expression changes and THM exposure. Thirty-seven before-after pairs were analyzed. The median increase from baseline to after swimming were: 0.7 to 2.3 for MET, and 1.4 to 7.1μg/m³ for exhaled total THM (sum of the four THM). Exhaled THM increased on average 0.94μg/m³ per 1 MET. While 1643 probes were differentially expressed post-exposure. Of them, 189 were also associated with exhaled levels of individual/total THM or MET after False Discovery Rate. The observed associations with the exhaled THM were low to moderate (Log-fold change range: -0.17 to 0.15). In conclusion, we identified short-term gene expression changes associated with swimming in a pool that were minor in magnitude and their biological meaning was unspecific. The high collinearity between exhaled THM levels and intensity of physical activity precluded mutually adjusted models with both covariates. These exploratory results should be validated in future studies.

Adhesion Molecules Associated with Female Genital Tract Infection

Efforts to develop vaccines that can elicit mucosal immune responses in the female genital tract against sexually transmitted infections have been hampered by an inability to measure immune responses in these tissues. The differential expression of adhesion molecules is known to confer site-dependent homing of circulating effector T cells to mucosal tissues. Specific homing molecules have been defined that can be measured in blood as surrogate markers of local immunity (e.g. α4β7 for gut). Here we analyzed the expression pattern of adhesion molecules by circulating effector T cells following mucosal infection of the female genital tract in mice and during a symptomatic episode of vaginosis in women. While CCR2, CCR5, CXCR6 and CD11c were preferentially expressed in a mouse model of Chlamydia infection, only CCR5 and CD11c were clearly expressed by effector T cells during bacterial vaginosis in women. Other homing molecules previously suggested as required for homing to the genital mucosa such as α4β1 and α4β7 were also differentially expressed in these patients. However, CD11c expression, an integrin chain rarely analyzed in the context of T cell immunity, was the most consistently elevated in all activated effector CD8⁺ T cell subsets analyzed. This molecule was also induced after systemic infection in mice, suggesting that CD11c is not exclusive of genital tract infection. Still, its increase in response to genital tract disorders may represent a novel surrogate marker of mucosal immunity in women, and warrants further exploration for diagnostic and therapeutic purposes.

Cyclic AMP signaling restricts activation and promotes maturation and antioxidant defenses in astrocytes

Background

cAMP signaling produces dramatic changes in astrocyte morphology and physiology. However, its involvement in phenotype acquisition and the transcriptionally mediated mechanisms of action are largely unknown.

Results

Here we analyzed the global transcriptome of cultured astroglial cells incubated with activators of cAMP pathways. A bulk of astroglial transcripts, 6221 annotated genes, were differentially regulated by cAMP signaling. cAMP analogs strongly upregulated genes involved in typical functions of mature astrocytes, such as homeostatic control, metabolic and structural support to neurons, antioxidant defense and communication, whereas they downregulated a considerable number of proliferating and immaturity-related transcripts. Moreover, numerous genes typically activated in reactive cells, such as scar components and immunological mediators, were repressed by cAMP. GSEA analysis contrasting gene expression profiles with transcriptome signatures of acutely isolated astrocytes and in situ evaluation of protein levels in these cells showed that cAMP signaling conferred mature and in vivo–like transcriptional features to cultured astrocytes.

Conclusions

These results indicate that cAMP signaling is a key pathway promoting astrocyte maturation and restricting their developmental and activation features. Therefore, a positive modulation of cAMP signaling may promote the normal state of differentiated astrocytes and favor the protection and function of neuronal networks.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2623-4) contains supplementary material, which is available to authorized users.

The yeast ζ-crystallin/NADPH:quinone oxidoreductase (Zta1p) is under nutritional control by the target of rapamycin pathway and is involved in the regulation of argininosuccinate lyase mRNA half-life

The yeast ζ-crystallin (Zta1p) is a quinone oxidoreductase belonging to the ζ-crystallin family, with activity in the reduction of alkenal/alkenone compounds. Various biological functions have been ascribed to the members of this protein family, such as their ability to interact specifically with AU-rich sequences in mRNA, and thus they have been proposed to act as AU-rich element-binding proteins (AREBPs). In this study, we evaluated the specificity of Zta1p for RNA versus DNA by means of a novel nonisotopic method for the in vitro quantitative detection of protein · RNA complexes. Through comparative transcriptomic analysis, we found that the lack of Zta1p negatively affects the expression of a group of genes involved in amino acid biosynthesis, the argininosuccinate lyase (ARG4) gene being one of them. Here, we propose that Zta1p participates in the post-transcriptional regulation of ARG4 expression by increasing the ARG4 mRNA half-life. In addition, expression of the ζ-crystallin gene (ZTA1) is itself regulated by nutrient availability through the general amino acid control and target of rapamycin pathways. Our results shed new light on the ζ-crystallin family members from yeast to humans as stress response proteins with a bifunctional role in the detoxification of alkenal and alkenone compounds, and the regulation of gene expression.

Microarray and deep sequencing cross-platform analysis of the mirRNome and isomiR variation in response to epidermal growth factor

Background

Epidermal Growth Factor (EGF) plays an important function in the regulation of cell growth, proliferation, and differentiation by binding to its receptor (EGFR) and providing cancer cells with increased survival responsiveness. Signal transduction carried out by EGF has been extensively studied at both transcriptional and post-transcriptional levels. Little is known about the involvement of microRNAs (miRNAs) in the EGF signaling pathway. miRNAs have emerged as major players in the complex networks of gene regulation, and cancer miRNA expression studies have evidenced a direct involvement of miRNAs in cancer progression.

Results

In this study, we have used an integrative high content analysis approach to identify the specific miRNAs implicated in EGF signaling in HeLa cells as potential mediators of cancer mediated functions. We have used microarray and deep-sequencing technologies in order to obtain a global view of the EGF miRNA transcriptome with a robust experimental cross-validation. By applying a procedure based on Rankprod tests, we have delimited a solid set of EGF-regulated miRNAs. After validating regulated miRNAs by reverse transcription quantitative PCR, we have derived protein networks and biological functions from the predicted targets of the regulated miRNAs to gain insight into the potential role of miRNAs in EGF-treated cells. In addition, we have analyzed sequence heterogeneity due to editing relative to the reference sequence (isomiRs) among regulated miRNAs.

Conclusions

We propose that the use of global genomic miRNA cross-validation derived from high throughput technologies can be used to generate more reliable datasets inferring more robust networks of co-regulated predicted miRNA target genes.

A myriad of miRNA variants in control and Huntington's disease brain regions detected by massively parallel sequencing

Huntington disease (HD) is a neurodegenerative disorder that predominantly affects neurons of the forebrain. We have applied the Illumina massively parallel sequencing to deeply analyze the small RNA populations of two different forebrain areas, the frontal cortex (FC) and the striatum (ST) of healthy individuals and individuals with HD. More than 80% of the small-RNAs were annotated as microRNAs (miRNAs) in all samples. Deep sequencing revealed length and sequence heterogeneity (IsomiRs) for the vast majority of miRNAs. Around 80-90% of the miRNAs presented modifications in the 3'-terminus mainly in the form of trimming and/or as nucleotide addition variants, while the 5'-terminus of the miRNAs was specially protected from changes. Expression profiling showed strong miRNA and isomiR expression deregulation in HD, most being common to both FC and ST. The analysis of the upstream regulatory regions in co-regulated miRNAs suggests a role for RE1-Silencing Transcription Factor (REST) and P53 in miRNAs downregulation in HD. The putative targets of deregulated miRNAs and seed-region IsomiRs strongly suggest that their altered expression contributes to the aberrant gene expression in HD. Our results show that miRNA variability is a ubiquitous phenomenon in the adult human brain, which may influence gene expression in physiological and pathological conditions.

Developmental expression of the oligodendrocyte myelin glycoprotein in the mouse telencephalon

The oligodendrocyte myelin glycoprotein is a glycosylphosphatidylinositol-anchored protein expressed by neurons and oligodendrocytes in the central nervous system. Attempts have been made to identify the functions of the myelin-associated inhibitory proteins (MAIPs) after axonal lesion or in neurodegeneration. However, the developmental roles of some of these proteins and their receptors remain elusive. Recent studies indicate that NgR1 and the recently discovered receptor PirB restrict cortical synaptic plasticity. However, the putative factors that trigger these effects are unknown. Because Nogo-A is mostly associated with the endoplasmic reticulum and myelin associated glycoprotein appears late during development, the putative participation of OMgp should be considered. Here, we examine the pattern of development of OMgp immunoreactive elements during mouse telencephalic development. OMgp immunoreactivity in the developing cortex follows the establishment of the thalamo-cortical barrel field. At the cellular level, we located OMgp neuronal membranes in dendrites and axons as well as in brain synaptosome fractions and axon varicosities. Lastly, the analysis of the barrel field in OMgp-deficient mice revealed that although thalamo-cortical connections were formed, their targeting in layer IV was altered, and numerous axons ectopically invaded layers II-III. Our data support the idea that early expressed MAIPs play an active role during development and point to OMgp participating in thalamo-cortical connections.

Regulation of GABA(A) and glutamate receptor expression, synaptic facilitation and long-term potentiation in the hippocampus of prion mutant mice

Background

Prionopathies are characterized by spongiform brain degeneration, myoclonia, dementia, and periodic electroencephalographic (EEG) disturbances. The hallmark of prioniopathies is the presence of an abnormal conformational isoform (PrP^sc) of the natural cellular prion protein (PrP^c) encoded by the Prnp gene. Although several roles have been attributed to PrP^c, its putative functions in neuronal excitability are unknown. Although early studies of the behavior of Prnp knockout mice described minor changes, later studies report altered behavior. To date, most functional PrP^c studies on synaptic plasticity have been performed in vitro. To our knowledge, only one electrophysiological study has been performed in vivo in anesthetized mice, by Curtis and coworkers. They reported no significant differences in paired-pulse facilitation or LTP in the CA1 region after Schaffer collateral/commissural pathway stimulation.

Methodology/Principal Findings

Here we explore the role of PrP^c expression in neurotransmission and neural excitability using wild-type, Prnp −/− and PrP^c-overexpressing mice (Tg20 strain). By correlating histopathology with electrophysiology in living behaving mice, we demonstrate that both Prnp −/− mice but, more relevantly Tg20 mice show increased susceptibility to KA, leading to significant cell death in the hippocampus. This finding correlates with enhanced synaptic facilitation in paired-pulse experiments and hippocampal LTP in living behaving mutant mice. Gene expression profiling using Illumina™ microarrays and Ingenuity pathways analysis showed that 129 genes involved in canonical pathways such as Ubiquitination or Neurotransmission were co-regulated in Prnp −/− and Tg20 mice. Lastly, RT-qPCR of neurotransmission-related genes indicated that subunits of GABA_A and AMPA-kainate receptors are co-regulated in both Prnp −/− and Tg20 mice.

Conclusions/Significance

Present results demonstrate that PrP^c is necessary for the proper homeostatic functioning of hippocampal circuits, because of its relationships with GABA_A and AMPA-Kainate neurotransmission. New PrP^c functions have recently been described, which point to PrP^c as a target for putative therapies in Alzheimer's disease. However, our results indicate that a “gain of function” strategy in Alzheimer's disease, or a “loss of function” in prionopathies, may impair PrP^c function, with devastating effects. In conclusion, we believe that present data should be taken into account in the development of future therapies.

Hnf1alpha (MODY3) controls tissue-specific transcriptional programs and exerts opposed effects on cell growth in pancreatic islets and liver

Heterozygous HNF1A mutations cause pancreatic-islet beta-cell dysfunction and monogenic diabetes (MODY3). Hnf1alpha is known to regulate numerous hepatic genes, yet knowledge of its function in pancreatic islets is more limited. We now show that Hnf1a deficiency in mice leads to highly tissue-specific changes in the expression of genes involved in key functions of both islets and liver. To gain insights into the mechanisms of tissue-specific Hnf1alpha regulation, we integrated expression studies of Hnf1a-deficient mice with identification of direct Hnf1alpha targets. We demonstrate that Hnf1alpha can bind in a tissue-selective manner to genes that are expressed only in liver or islets. We also show that Hnf1alpha is essential only for the transcription of a minor fraction of its direct-target genes. Even among genes that were expressed in both liver and islets, the subset of targets showing functional dependence on Hnf1alpha was highly tissue specific. This was partly explained by the compensatory occupancy by the paralog Hnf1beta at selected genes in Hnf1a-deficient liver. In keeping with these findings, the biological consequences of Hnf1a deficiency were markedly different in islets and liver. Notably, Hnf1a deficiency led to impaired large-T-antigen-induced growth and oncogenesis in beta cells yet enhanced proliferation in hepatocytes. Collectively, these findings show that Hnf1alpha governs broad, highly tissue-specific genetic programs in pancreatic islets and liver and reveal key consequences of Hnf1a deficiency relevant to the pathophysiology of monogenic diabetes.

Hypothalamus transcriptome profile suggests an anorexia-cachexia syndrome in the anx/anx mouse model

The anx/anx mouse displays poor appetite and lean appearance and is considered a good model for the study of anorexia nervosa. To identify new genes involved in feeding behavior and body weight regulation we performed an expression profiling in the hypothalamus of the anx/anx mice. Using commercial microarrays we detected 156 differentially expressed genes and validated 92 of those using TaqMan low-density arrays. The expression of a set of 87 candidate genes selected based on literature evidences was also quantified by TaqMan low-density arrays. Our results showed enrichment in deregulated genes involved in cell death, cell morphology, and cancer, as well as an alteration of several signaling circuits involved in energy balance including neuropeptide Y and melanocortin signaling. The expression profile along with the phenotype led us to conclude that anx/anx mice resemble the anorexia-cachexia syndrome typically observed in cancer, infection with human immunodeficiency virus or chronic diseases, rather than starvation, and that anx/anx mice could be considered a good model for the treatment and investigation of this condition.

Gene expression profiling distinguishes JAK2V617F-negative from JAK2V617F-positive patients in essential thrombocythemia

To explore the gene expression signature in essential thrombocythemia (ET) patients in relation to JAK2V617F mutational status, expression profiling in circulating granulocytes was performed. Twenty ET were studied by microarray analysis and the results were confirmed by real-time quantitative RT-PCR in 40 ET patients, not receiving cytoreductive treatment. A heterogeneous molecular signature characterized by two main gene expression patterns was found: one with an upregulation of inflammatory genes related to neutrophil activation and thrombosis, and the other with significantly lower expression of these genes. Supervised clustering analysis showed 30 genes differentially expressed between JAK2V617F-negative and JAK2V617F-positive ET patients. Among the JAK2V617F-negative, a set of 14 genes (CISH, C13orf18, CCL3, PIM1, MAFF, SOCS3, ID2, GADD45B, KLF5, TNF, LAMB3, HRH4, TAGAP and TRIB1) showed an abnormal expression pattern. In this group of patients, CISH, SOCS2, SOCS3 and PIM1 genes, all involved in JAK-STAT signalling pathway, presented a lower expression. A two-gene predictor model was built comprising FOSB and CISH genes, which were the best discriminators of JAK2V617F status. In conclusion, JAK2V617F-negative ET patients present a characteristic gene expression profile, different from JAK2V617F-positive patients. Other pathways, besides JAK-STAT, might be implicated in the pathophysiology of JAK2V617F-negative ET patients.

Transcriptional regulation by poly(ADP-ribose) polymerase-1 during T cell activation

Background

Accumulating evidence suggests an important role for the enzyme poly(ADP-ribose) polymerase-1 (PARP-1) as an integral part of the gene expression regulatory machinery during development and in response to specific cellular signals. PARP-1 might modulate gene expression through its catalytic activity leading to poly(ADP-ribosyl)ation of nuclear proteins or by its physical association with relevant proteins. Recently, we have shown that PARP-1 is activated during T cell activation. However, the proposed role of PARP-1 in reprogramming T cell gene expression upon activation remains largely unexplored.

Results

In the present study we use oligonucleotide microarray analysis to gain more insight into the role played by PARP-1 during the gene expression reprogramming that takes place in T cells upon activation with anti-CD3 stimulation alone, or in combination with anti-CD28 co-stimulation. We have identified several groups of genes with expression modulated by PARP-1. The expression of 129 early-response genes to anti-CD3 seems to be regulated by PARP-1 either in a positive (45 genes) or in a negative manner (84 genes). Likewise, in the presence of co-stimulation (anti-CD3 + anti-CD28 stimulation), the expression of 203 genes is also regulated by PARP-1 either up (173 genes) or down (30 genes). Interestingly, PARP-1 deficiency significantly alters expression of genes associated with the immune response such as chemokines and genes involved in the Th1/Th2 balance.

Conclusion

This study provides new insights into changes in gene expression mediated by PARP-1 upon T cell activation. Pathway analysis of PARP-1 as a nuclear signalling molecule in T cells would be of relevance for the future development of new therapeutic approaches targeting PARP-1 in the acquired immune response.

Genomic imbalances in Schistosoma-associated and non–Schistosoma-associated bladder carcinoma. An array comparative genomic hybridization analysis

Carcinoma of the urinary bladder is the most common malignancy in many tropical and subtropical countries due to endemic infection by Schistosoma hematobium (bilharzia). In the current study, we performed a high-resolution analysis of gene copy number amplifications using array comparative genomic hybridization to compare DNA copy number changes in pools of Schistosoma-associated (SA) and non-Schistosoma-associated (NSA) bladder cancer (BC). Many DNA copy number changes were detected in all studies, with multiple gains and losses of genetic material. The most frequent alterations were gains on 5p15.2 approximately p15.33, 8q13.1, and 11q13, and losses on 8p21.3 approximately p22 and 22q13. Even when SA pools showed no Schistosoma-specific gene copy number profiling as compared to NSA pools, some genes seemed to be gained (ELN on 7q11.23) and some lost (PRKAG3 on 2q35 and PRDM6 on 5q23.2) in SA-SCC. The following genes were gained in all histopathologic categories: SRC (20q11.23), CEBPB (20q13.13), and GPR9 (Xq13.1). Our study did not provide clear evidence of differences in carcinogenesis of SA-BC and NSA-BC.

Partially Degraded RNA from Bladder Washing is a Suitable Sample for Studying Gene Expression Profiles in Bladder Cancer

OBJECTIVES

To determine the impact of different levels of RNA degradation on gene expression measurements and to ascertain if the gene expression profile obtained from bladder washing (BW) correlates to that obtained from the related bladder tumour (BT).

METHODS

BT and BW RNAs from the same patient were heat shocked to obtain three RNA degradation states, which were compared with intact RNAs from healthy bladders by using complementary DNA (cDNA) microarrays. All samples were amplified by means of a T3N9-based transcription method. In addition, four of the differentially expressed genes in microarrays related to bladder cancer (KRT20, IGF2, GSN, and CCL2) were analyzed in 36 tumoural and 14 control BW samples by quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR).

RESULTS

A high percentage of overlapping differentially expressed genes were detected between BT arrays (85-91%) and between BW arrays (78-93%). Furthermore, the similarity between BW and BT arrays was relatively high and independent of the RNA degradation state (52-60%). Finally, expression differences for the four selected genes were confirmed in the vast majority of extended BW samples tested by qRT-PCR.

CONCLUSIONS

Our results showed that partially degraded RNA samples analyzed by cDNA microarrays yielded gene expression profiles comparable to those obtained using intact RNA. Moreover, BW RNA exhibited gene expression patterns similar to those identified in the BT, indicating that BW is an appropriate sample for studying gene expression profiles of BT using cDNA microarrays. In addition, qRT-PCR results further support the suitability of BW for gene expression profiling and its potential use for routine diagnostics.

Genome-wide differences between microsatellite stable and unstable colorectal tumors

Genomic copy number changes are frequently found in cancers and they have been demonstrated to contribute to carcinogenesis; and it is widely accepted that tumors with microsatellite instability (MSI) are genetically stable and mostly diploid. In the present study we compared the copy number alterations and the gene-expression profiles of microsatellite stable (MSS) and MSI colorectal tumors. A total number of 31 fresh-frozen primary tumors (16 MSS and 15 MSI) were used. Twenty-eight samples (15 MSS and 13 MSI) were analyzed with metaphase comparative genomic hybridization (CGH), nine of which plus one additional sample (4 MSS and 6 MSI) were further analyzed by cDNA-based array-CGH. Gene expression analysis was performed with six samples [3 MSS and 3 MSI, four of these used in metaphase CGH (mCGH) analysis] to identify differentially expressed genes possibly located in the lost or amplified regions found by CGH, stressing the biological significance of copy number changes. Metaphase and array-CGH analysis of two colon cancer cell lines (HTC116 and SW480, reported as MSI and MSS archetypes) gave comparable results. Alterations found by mCGH in MSS tumors were +20, +8q, -8p and -18q. Interestingly, 1p22, 4q26 and 15q21 were found deleted preferentially in MSS tumors, while 22q13 was found gained in MSI tumors. The regions of alterations identified by array-CGH were gains at 8q24, 16q24.3 and 20q13, and the loss of 5q21, appearing in the both types of tumors. Gene expression analysis revealed genes with specific associations with the copy number changes of the corresponding genomic regions. As a conclusion, colorectal cancer is a heterogeneous disease, demonstrated by the genomic profiles of individual samples. However, our data shows that copy number changes do not occur exclusively in the MSS phenotypes.

Human vitamin K-dependentGAS6: Gene structure, allelic variation, and association with stroke

The product of the growth arrest-specific gene 6 (GAS6), a ligand for the Axl, Sky, and Mer tyrosine kinase receptors, is a vitamin K-dependent protein, structurally related to anticoagulant protein S. Gas6-deficient mice are protected against thrombosis, demonstrating the importance of this protein in the cardiovascular system. The present study was aimed at determining the human GAS6 intron-exon structure and analyzing the gene for the presence of allelic variants that could be associated with atherothrombotic disease. Online analyses allowed us to localize 15 GAS6 exons and to determine the sequence of their intron-flanking regions, in a chromosome 13 region spanning 43.8 kb of DNA. SSCP analysis of PCR-amplified GAS6 exons with their intron-flanking regions from a minimum of 12 control DNA samples, revealed the presence of eight different variants, which were confirmed to be single nucleotide polymorphisms (SNPs). Three of them (c.1263G>C, c.1332C>T, and c.1869T>C) are localized in exons 11, 12, and 14, and appear to be neutral since they do not modify the encoded amino acid. The other SNPs (c.280+170C>G, c.712+26G>A, c.713-155C>T, c.834+7G>A, and c.1478-94C>G) are in introns 3, 7, 8, and 12. A preliminary analysis of five of these SNPs in a group of 110 healthy controls and 188 patients with atherothrombotic disease has revealed statistically significant differences between controls and stroke patients in the allelic distributions of one of these variants (c.834+7G>A in intron 8). The SNP identification in GAS6 reported here would be very useful in future association studies aimed at determining the physiologic role of GAS6 in stroke and other human diseases.

LRRN6A/LERN1 (leucine-rich repeat neuronal protein 1), a novel gene with enriched expression in limbic system and neocortex

Human chromosome 15q24-q26 is a very complex genomic region containing several blocks of segmental duplications to which susceptibility to anxiety disorders has been mapped (Gratacos et al., 2001, Cell, 106, 367-379; Pujana et al., 2001, Genome Res., 11, 98-111). Through an in silico gene content analysis of the 15q24-q26 region we have identifie1d a novel gene, LRRN6A (leucine-rich repeat neuronal 6A), and confirmed its location to the centromeric end of this complex region. LRRN6A encodes a transmembrane leucine-rich repeat protein, LERN1 (leucine-rich repeat neuronal protein 1), with similarity to proteins involved in axonal guidance and migration, nervous system development and regeneration processes. The identification of homologous genes to LRRN6A on chromosomes 9 and 19 and the orthologous genes in the mouse genome and other organisms suggests that LERN proteins constitute a novel subfamily of LRR (leucine-rich repeat)-containing proteins. The LRRN6A expression pattern is specific to the central nervous system, highly and broadly expressed during early stages of development and gradually restricted to forebrain structures as development proceeds. Expression level in adulthood is lower in general but remains stable and significantly enriched in the limbic system and cerebral cortex. Taken together, the confirmation of LRRN6A's expression profile, its predicted protein structure and its similarity to nervous system-expressed LRR proteins with essential roles in nervous system development and maintenance suggest that LRRN6A is a novel gene of relevance in the molecular and cellular neurobiology of vertebrates.

The MAPK Hog1 recruits Rpd3 histone deacetylase to activate osmoresponsive genes

Regulation of gene expression by mitogen-activated protein kinases (MAPKs) is essential for proper cell adaptation to extracellular stimuli. Exposure of yeast cells to high osmolarity results in rapid activation of the MAPK Hog1, which coordinates the transcriptional programme required for cell survival on osmostress. The mechanisms by which Hog1 and MAPKs in general regulate gene expression are not completely understood, although Hog1 can modify some transcription factors. Here we propose that Hog1 induces gene expression by a mechanism that involves recruiting a specific histone deacetylase complex to the promoters of genes regulated by osmostress. Cells lacking the Rpd3-Sin3 histone deacetylase complex are sensitive to high osmolarity and show compromised expression of osmostress genes. Hog1 interacts physically with Rpd3 in vivo and in vitro and, on stress, targets the deacetylase to specific osmostress-responsive genes. Binding of the Rpd3-Sin3 complex to specific promoters leads to histone deacetylation, entry of RNA polymerase II and induction of gene expression. Together, our data indicate that targeting of the Rpd3 histone deacetylase to osmoresponsive promoters by the MAPK Hog1 is required to induce gene expression on stress.

Local tissue interactions across the dorsal midline of the forebrain establish CNS laterality

The mechanisms that establish behavioral, cognitive, and neuroanatomical asymmetries are poorly understood. In this study, we analyze the events that regulate development of asymmetric nuclei in the dorsal forebrain. The unilateral parapineal organ has a bilateral origin, and some parapineal precursors migrate across the midline to form this left-sided nucleus. The parapineal subsequently innervates the left habenula, which derives from ventral epithalamic cells adjacent to the parapineal precursors. Ablation of cells in the left ventral epithalamus can reverse laterality in wild-type embryos and impose the direction of CNS asymmetry in embryos in which laterality is usually randomized. Unilateral modulation of Nodal activity by Lefty1 can also impose the direction of CNS laterality in embryos with bilateral expression of Nodal pathway genes. From these data, we propose that laterality is determined by a competitive interaction between the left and right epithalamus and that Nodal signaling biases the outcome of this competition.

Cloning, mapping and expression analysis of C15orf4, a novel human gene with homology to the yeast mitochondrial ribosomal protein Ym130 gene

We have identified C15orf4, a novel human gene showing homology to the yeast mitochondrial ribosomal protein YmL30. C15orf4 encodes a transcript of 1,006 nt with an ORF of 279 amino acids and a predicted protein size of 31.7 kDa. Expression of C15orf4 is enriched in testis. C15orf4 was positioned to chromosome 15q24 by radiation hybrid mapping. We have identified the C15orf4 mouse orthologue as well as homologues in other species.

MSX2 expression in the apical ectoderm ridge is regulated by an MSX2 and Dlx5 binding site

The apical ectodermal ridge (AER) is a specialized ectodermal region essential for limb outgrowth. Msx2 expression patterns in limb development strongly suggest an important role for Msx2 in the AER. Our previous studies identified a 348-bp fragment of the chicken Msx2 gene with AER enhancer activity. In this study, the functions of four potential homeodomain binding TAAT sites in this enhancer were studied using transgenic mice and in vitro protein-DNA interactions. Transgenic studies indicate that the four TAAT sites are not redundant and that only the B-TAAT site is critical for AER enhancer activity. The expression patterns of Msx2 and Dlx5 genes in the AER suggest that they might be involved in the regulation of Msx2. In support of this hypothesis, we found that Msx2 and Dlx5 can bind to the B-TAAT site as well as to a fragment containing the D- and E-TAAT sites in the Msx2 AER enhancer sequences. (c)2002 Elsevier Science (USA).

ATB(0)/SLC1A5 gene. Fine localisation and exclusion of association with the intestinal phenotype of cystic fibrosis

The Na+-dependent amino acid transporter named ATB(0) was previously found to be located in 19q13.3 by fluorescence in situ hybridisation. Genetic heterogeneity in the 19q13.2-13.4 region, syntenic to the Cystic Fibrosis Modulator Locus 1 (CFM1) in mouse, seemed to be associated to the intestinal phenotypic variation of cystic fibrosis (CF). We performed fine chromosomal mapping of ATB(0) on radiation hybrid (RH) panels G3 and TNG. Based on the most accurate location results from TNG-RH panel, mapping analysis evidenced that ATB(0) is localised between STS SHGC-13875 (D19S995) and STS SHGC-6138 in 19q13.3, that corresponds with the immediately telomeric/distal segment of the strongest linkage region within the human CFM1 (hCFM1) syntenic region. Regarding to the genomic structure and exon organisation, our results show that the ATB(0) gene is organised into eight exons. The knowledge of the genomic structure allowed us to perform an exhaustive mutational analysis of the gene. Evaluation of the possible implication of ATB(0) in the intestinal phenotype of CF was performed on the basis of the functional characteristics of the encoded protein, its apparent relevance to meconium ileus (MI) and position in relation to the hCFM1 syntenic region. We have analysed this gene in samples from CF patients with and without MI. Several sequence variations in the ATB(0) gene were identified, although none of them seemed to be related to the intestinal phenotype of CF. Even though no particular allele or haplotype in ATB(0) appears to be associated to CF-MI disease, new SNPs identified should be useful in segregation and linkage disequilibrium analyses in families affected by other disorders caused by the impairment of neutral amino acid transport.

PALML, a novel paralemmin-related gene mapping on human chromosome 1p21

We describe PALML, a novel gene encoding a 551 amino acid protein with similarity to paralemmin and the paralemmin-like amino terminal domain of AKAP2, a protein kinase A anchor protein. PALML mRNA is expressed in many tissues and is most abundant in cardiac and skeletal muscle, while absent from brain and blood. Exogenously expressed PALML fusion protein has a widespread cytoplasmic localization, and it is excluded from the nucleus. Human PALML maps on human chromosome 1p21 (between D1S2767 and D1S223). SSCP-HD analysis of exonic sequences in patients with VUR (familial non-syndromic vesicoureteral reflux syndrome) excluded mutations in the PALML gene from causing this disease. PALML, paralemmin and AKAP2 share the presence of a conserved coiled coil region that may mediate protein interactions with shared partners. Based on its resemblance to paralemmin and AKAP2, PALML is hypothesized to be involved in regulating intracellular signaling and membrane-cytoskeletal interactions.

Identification of C15orf5, a heart-enriched transcript on chromosome 15q23-q24

We have isolated C15orf5, a novel human gene lacking homology to any known protein. The C15orf5 gene encodes a transcript of 1,519 nt with an ORF of 94 amino acids and a predicted protein size of 11.5 kDa. Northern blot analysis showed enhanced expression of C15orf5 in heart. C15orf5 was mapped to chromosome 15q23-q24 using the Stanford TNG4 Radiation Hybrid panel.

Identification and expression analysis of C3orf1, a novel human gene homologous to the Drosophila RP140-upstream gene

We have isolated C3orf1, a novel gene on human chromosome 3 showing homology to the Drosophila RP140-upstream gene. When mutated, RP140-upstream causes lethality in flies through an unknown mechanism, perhaps by interfering with transcription of the RP140 RNA polymerase subunit. The human C3orf1 gene encodes a predicted membrane protein of 32.2 kDa with four transmembrane domains without any other known motifs. Northern blot analysis showed generalized expression of C3orf1, enhanced in heart and skeletal muscle. EST database searching revealed the existence of a homologue gene in mouse. Thus, the C3orf1 gene is conserved and may perform an essential function in all tissues in mammals.

Characterization of human FSD1, a novel brain specific gene on chromosome 19 with paralogy to 9q31

We have characterized a novel human gene, FSD1, on chromosome 19. FSD1 has a BBC, FN3 and SPRY domain, it is distantly related to the midline 1 gene and is expressed only in the brain. We have established its exon-intron structure and we have identified the corresponding orthologous genes in other species. In addition, the identification of FSD1 has led us to identify a homologous counterpart sequence on chromosome 9.

Cloning of the novel gene TM6SF1 reveals conservation of clusters of paralogous genes between human chromosomes 15q24-->q26 and 19p13.3-->p12

As the result of the EUROIMAGE Consortium sequencing project, we have isolated and characterized a novel gene on chromosome 15, TM6SF1. It encodes a 370 amino acid product with enhanced expression in spleen, testis and peripheral blood leukocytes. We have identified another gene, paralogous to TM6SF1 on chromosome 19p12, TM6SF2, with an overall similarity of 68% and 52% identity at the protein level. This conservation has led us to uncover a series of eleven genes in 19p13.3-->p12 with close homology to genes in 15q24--> q26. The percentage of sequence similarity between each paralogous pair of genes at the protein level ranges between 43 and 89%. A partial conservation of synteny with mouse chromosomes 7, 8 and 9 is also observed. The corresponding orthologous genes in mouse of human TM6SF1 and TM6SF2 show a high degree of amino acid sequence conservation.