"Role of long non-coing RNAs in the stressed heart". Wisper : a novel cardiac fibroblast-enriched IncRNA that modulate fibrosis in the failing heart

LncRNAs are transcripts greater than 200 nucleotides in length with no apparent coding potential. They exert important regulatory functions in the genome. Their role in cardiac fibrosis is however unexplored. To identify IncRNAs that could modulate cardiac fibrosis, we profiled the long non-coding transcriptome in the infarcted mouse heart, and identified 1500 novel IncRNAs. These IncRNAs have unique characteristics such as high tissue and cell type specificity. Their expression is highly correlated with parameters of cardiac dimensions and function. The majority of these novel IncRNAs are conserved in human. Importantly, human IncRNAs appear to be differentially expressed in heart disease. Using a computational pipeline, we identified a super-enhancer-associated IncRNA, which is dynamically expressed after myocardial infarction. We named this particular transcript Wisper for «Wisp2 super-enhancer- derived IncRNA ». Interestingly, Wisper expression is overexpressed in cardiac fibroblasts as compared to cardiomyocytes or to fibroblasts isolated from other organs than the heart. The importance of Wisper in the biology of fibroblasts was demonstrated in knockdown experiments. Differentiation of cardiac fibroblast into myofibroblasts in vitro is significantly impaired upon Wisper knockdown. Wisper downregulation in cardiac fibroblasts results in a dramatic reduction of fibrotic gene expression, a diminished cell proliferation and an increase in apoptotic cell death. In vivo, depletion of Wisper during the acute phase of the response to infarction is detrimental via increasing the risk of cardiac rupture. On the other hand, Wisper knockdown following infarction in a prevention study reduces fibrosis and preserves cardiac function. Since WISPER is detectable in the human heart, where it is associated with severe cardiac fibrosis, these data suggest that Wisper could represent a novel therapeutic target for limiting the extent of the fibrotic response in the heart. -- Les long ARN non-codants (IncRNAs) sont des ARN de plus de 200 nucléotides qui ne codent pas pour des protéines. Ils exercent d'importantes fonctions dans le génome. Par contre, leur importance dans le développement de la fibrose cardiaque n'a pas été étudiée. Pour identifier des IncRNAs jouant un rôle dans ce processus, le transcriptome non-codant a été étudié dans le coeur de'souris après un infarctus du myocarde. Nous avons découverts 1500 nouveaux IncRNAs. Ces transcrits ont d'uniques caractéristiques. En particulier ils sont extrêmement spécifiques de sous-populations de cellules cardiaques. Par ailleurs, leur expression est remarquablement corrélée avec les paramètres définissant les dimensions du coeur et la fonction cardiaque. La majorité de ces IncRNAs sont conservés chez l'humain. Certains sont modulés dans des pathologies cardiaques. En utilisant une approche bioinformatique, nous avons identifié un IncRNA qui est associé à des séquences amplificatrices et qui est particulièrement enrichi dans les fibroblastes cardiaques. Ce transcrit a été nommé Wisper pour «Wisp2 super-enhancer-derived IncRNA ». L'importance de Wisper dans la biologie des fibroblastes cardiaques est démontrée dans des expériences de déplétion. En l'absence de Wisper, l'expression de protéines impliquées dans le développement de la fibrose est dramatiquement réduite dans les fibroblastes cardiaques. Ceux-ci montrent une prolifération réduite. Le niveau d'apoptose est largement augmenté. In vivo, la déplétion de Wisper pendant la phase aiguë de l'infarctus rehausse le risque de rupture cardiaque. Au contraire, la réduction de l'expression de Wisper pendant la phase chronique diminue la fibrose cardiaque et améliore la fonction du coeur. Puisque Wisper est exprimé dans le coeur humain, ce transcrit représente une nouvelle cible thérapeutique pour limiter la réponse fibrotique dans le coeur.

Annotation des ARN non codants du génome de Candida albicans par méthode bioinformatique

La bio-informatique est un champ pluridisciplinaire qui utilise la biologie, l’informatique, la physique et les mathématiques pour résoudre des problèmes posés par la biologie. L’une des thématiques de la bio-informatique est l’analyse des séquences génomiques et la prédiction de gènes d’ARN non codants. Les ARN non codants sont des molécules d’ARN qui sont transcrites mais pas traduites en protéine et qui ont une fonction dans la cellule. Trouver des gènes d’ARN non codants par des techniques de biochimie et de biologie moléculaire est assez difficile et relativement coûteux. Ainsi, la prédiction des gènes d’ARNnc par des méthodes bio-informatiques est un enjeu important. Cette recherche décrit un travail d’analyse informatique pour chercher des nouveaux ARNnc chez le pathogène Candida albicans et d’une validation expérimentale. Nous avons utilisé comme stratégie une analyse informatique combinant plusieurs logiciels d’identification d’ARNnc. Nous avons validé un sous-ensemble des prédictions informatiques avec une expérience de puces à ADN couvrant 1979 régions du génome. Grace à cette expérience nous avons identifié 62 nouveaux transcrits chez Candida albicans. Ce travail aussi permit le développement d’une méthode d’analyse pour des puces à ADN de type tiling array. Ce travail présente également une tentation d’améliorer de la prédiction d’ARNnc avec une méthode se basant sur la recherche de motifs d’ARN dans les séquences.

Rôle du régulateur Fur et des petits ARN non codants RfrA et RfrB dans l’homéostasie du fer et la virulence de Salmonella

La régulation de l’homéostasie du fer est cruciale chez les bactéries. Chez Salmonella, l’expression des gènes d’acquisition et du métabolisme du fer au moment approprié est importante pour sa survie et sa virulence. Cette régulation est effectuée par la protéine Fur et les petits ARN non codants RfrA et RfrB. Le rôle de ces régulateurs est d’assurer que le niveau de fer soit assez élevé pour la survie et le métabolisme de Salmonella, et assez faible pour éviter l’effet toxique du fer en présence d’oxygène. Les connaissances concernant le rôle de ces régulateurs ont été principalement obtenues par des études chez S. Typhimurium, un sérovar généraliste causant une gastro-entérite chez les humains. Très peu d’informations sont connues sur le rôle de ces régulateurs chez S. Typhi, un sérovar humain-spécifique responsable de la fièvre typhoïde. Le but de cette étude était de déterminer les rôles de Fur, RfrA et RfrB dans l’homéostasie du fer et la virulence de Salmonella, et de démontrer qu’ils ont une implication distincte chez les sérovars Typhi et Typhimurium. Premièrement, Fur, RfrA et RfrB régulent l’homéostasie du fer de Salmonella. Les résultats de cette étude ont démontré que Fur est requis pour la résistance au stress oxydatif et pour une croissance optimale dans différentes conditions in vitro. La sensibilité du mutant fur est due à l’expression des petits ARN RfrA et RfrB, et cette sensibilité est beaucoup plus importante chez S. Typhi que chez S. Typhimurium. Également, Fur inhibe la transcription des gènes codant pour les sidérophores en conditions riches en fer, tandis que les petits ARN RfrA et RfrB semblent être importants pour la production d’entérobactine et de salmochélines chez S. Typhi lors de conditions pauvres en fer. Ensuite, ces régulateurs affectent la virulence de Salmonella. Fur est important pour la motilité de Salmonella, particulièrement chez S. Typhi. Fur est nécessaire pour l’invasion des deux sérovars dans les cellules épithéliales, et pour l’entrée et la survie de S. Typhi dans les macrophages. Chez S. Typhimurium, Fur ne semble pas impliqué dans l’interaction avec les macrophages. De plus, les petits ARN RfrA et RfrB sont importants pour la multiplication intracellulaire de Salmonella dans les macrophages pour les deux sérovars. Finalement, la protéine Fur et les petits ARN RfrA et RfrB régulent l’expression de l’opéron fimbriaire tcf, absent du génome de S. Typhimurium. Un site de liaison putatif de la protéine Fur a été identifié dans la région promotrice de tcfA chez S. Typhi, mais une régulation directe n’a pas été confirmée. L’expression de tcf est induite par le fer et par Fur, et est inhibée par les petits ARN RfrA et RfrB. Ainsi, ces régulateurs affectent des gènes de virulence qui sont retrouvés spécifiquement chez S. Typhi. En somme, ce projet a permis de démontrer que les régulateurs de l’homéostasie du fer de Salmonella peuvent affecter la résistance de cette bactérie pathogène à différents stress, notamment le stress oxydatif, la croissance en conditions de carence en fer ainsi que la virulence. Ces régulateurs jouent un rôle distinct chez les sérovars Typhi et Typhimurium.

Origin of avian genome size and structure in non-avian dinosaurs

Avian genomes are small and streamlined compared with those of other amniotes by virtue of having fewer repetitive elements and less non-coding DNA(1,2). This condition has been suggested to represent a key adaptation for flight in birds, by reducing the metabolic costs associated with having large genome and cell sizes(3,4). However, the evolution of genome architecture in birds, or any other lineage, is difficult to study because genomic information is often absent for long-extinct relatives. Here we use a novel bayesian comparative method to show that bone-cell size correlates well with genome size in extant vertebrates, and hence use this relationship to estimate the genome sizes of 31 species of extinct dinosaur, including several species of extinct birds. Our results indicate that the small genomes typically associated with avian flight evolved in the saurischian dinosaur lineage between 230 and 250 million years ago, long before this lineage gave rise to the first birds. By comparison, ornithischian dinosaurs are inferred to have had much larger genomes, which were probably typical for ancestral Dinosauria. Using comparative genomic data, we estimate that genome-wide interspersed mobile elements, a class of repetitive DNA, comprised 5 - 12% of the total genome size in the saurischian dinosaur lineage, but was 7 - 19% of total genome size in ornithischian dinosaurs, suggesting that repetitive elements became less active in the saurischian lineage. These genomic characteristics should be added to the list of attributes previously considered avian but now thought to have arisen in non-avian dinosaurs, such as feathers(5), pulmonary innovations 6, and parental care and nesting

Non-classical HLA-E gene variability in Brazilians: a nearly invariable locus surrounded by the most variable genes in the human genome

The non-classical human leukocyte antigen (HLA) class I genes present a very low rate of variation. So far, only 10 HLA-E alleles encoding three proteins have been described, but only two are frequently found in worldwide populations. Because of its historical background, Brazilians are very suitable for population genetic studies. Therefore, 104 bone marrow donors from Brazil were evaluated for HLA-E exons 14. Seven variation sites were found, including two known single nucleotide polymorphisms (SNPs) at positions +424 and +756 and five new SNPs at positions +170 (intron 1), +1294 (intron 3), +1625, +1645 and +1857 (exon 4). Haplotyping analysis did show eight haplotypes, three of them known as E*01:01:01, E*01:03:01 and E*01:03:02:01 and five HLA-E new alleles that carry the new variation sites. The HLA-E*01:01:01 allele was the predominant haplotype (62.50%), followed by E*01:03:02:01 (24.52%). Selective neutrality tests have disclosed an interesting pattern of selective pressures in which balancing selection is probably shaping allele frequency distributions at an SNP at exon 3 (codon 107), sequence diversity at exon 4 and the non-coding regions is facing significant purifying pressure. Even in an admixed population such as the Brazilian one, the HLA-E locus is very conserved, presenting few polymorphic SNPs in the coding region.

RNA non-codificanti indotti da danno al DNA regolano il fattore di trascrizione HIF-1α

Recenti analisi sull’intero trascrittoma hanno rivelato una estensiva trascrizione di RNA non codificanti (ncRNA), le quali funzioni sono tuttavia in gran parte sconosciute. In questo lavoro è stato dimostrato che alte dosi di camptotecina (CPT), un farmaco antitumorale inibitore della Top1, aumentano la trascrizione di due ncRNA antisenso in 5’ e 3’ (5'aHIF-1α e 3'aHIF-1α rispettivamente) al locus genico di HIF-1α e diminuiscono i livelli dell’mRNA di HIF-1α stesso. Gli effetti del trattamento sono Top1-dipendenti, mentre non dipendono dal danno al DNA alla forca di replicazione o dai checkpoint attivati dal danno al DNA. I ncRNA vengono attivati in risposta a diversi tipi di stress, il 5'aHIF-1α è lungo circa 10 kb e possiede sia il CAP in 5’ sia poliadenilazione in 3’ (in letteratura è noto che il 3'aHIF-1α è un trascritto di 1,7 kb, senza 5’CAP né poliadenilazione). Analisi di localizzazione intracellulare hanno dimostrato che entrambi sono trascritti nucleari. In particolare 5'aHIF-1α co-localizza con proteine del complesso del poro nucleare, suggerendo un suo possibile ruolo come mediatore degli scambi della membrana nucleare. È stata dimostrata inoltre la trascrizione dei due ncRNA in tessuti di tumore umano del rene, evidenziandone possibili ruoli nello sviluppo del cancro. È anche noto in letteratura che basse dosi di CPT in condizioni di ipossia diminuiscono i livelli di proteina di HIF-1α. Dopo aver dimostrato su diverse linee cellulari che i due ncRNA sopracitati non potessero essere implicati in tale effetto, abbiamo studiato le variazioni dell’intero miRnoma alle nuove condizioni sperimentali. In tal modo abbiamo scoperto che il miR-X sembra essere il mediatore molecolare dell’abbattimento di HIF-1α dopo trattamento con basse dosi di CPT in ipossia. Complessivamente, questi risultati suggeriscono che il fattore di trascrizione HIF-1α venga finemente regolato da RNA non-codificanti indotti da danno al DNA.

miR-34a and miR-15a/16 are co-regulated in non-small cell lung cancer and control cell cycle progression in a synergistic and Rb-dependent manner

microRNAs (miRNAs) are small non-coding RNAs that are frequently involved in carcinogenesis. Although many miRNAs form part of integrated networks, little information is available how they interact with each other to control cellular processes. miR-34a and miR-15a/16 are functionally related; they share common targets and control similar processes including G1-S cell cycle progression and apoptosis. The aim of this study was to investigate the combined action of miR-34a and miR-15a/16 in non-small cell lung cancer (NSCLC) cells.

Mapping new non-genetic dependencies in malignant pleural mesothelioma

Malignant Pleural Mesothelioma (MPM) is a very aggressive cancer whose incidence is growing worldwide. MPM escapes the classical models of carcinogenesis and lacks a distinctive genetic fingerprint, keeping obscure the molecular events that lead to tumorigenesis. This severely impacts on the limited therapeutic options and on the lack of specific biomarkers, concurring to make MPM one of the deadliest cancers. Here we combined a functional genome-wide loss of function CRISPR/Cas9 screening with patients’ transcriptomic and clinical data, to identify genes essential for MPM progression. Besides, we explored the role of non-coding RNAs to MPM progression by analysing gene expression profiles and clinical data from the MESO-TCGA dataset. We identified TRIM28 and the lncRNA LINC00941 as new vulnerabilities of MPM, associated with disease aggressiveness and bad outcome of patients. TRIM28 is a multi-domain protein involved in many processes, including transcription regulation. We showed that TRIM28 silencing impairs MPM cells’ growth and clonogenicity by blocking cells in mitosis. RNA-seq profiling showed that TRIM28 loss abolished the expression of major mitotic players. Our data suggest that TRIM28 is part of the B-MYB/FOXM1-MuvB complex that specifically drives the activation of mitotic genes, keeping the time of mitosis. In parallel, we found LINC00941 as strongly associated with reduced survival probability in MPM patients. LINC00941 KD profoundly reduced MPM cells’ growth, migration and invasion. This is accompanied by changes in morphology, cytoskeleton organization and cell-cell adhesion properties. RNA-seq profiling showed that LINC00941 KD impacts crucial functions of MPM, including HIF1α signalling. Collectively these data provided new insights into MPM biology and demonstrated that the integration of functional screening with patients’ clinical data is a powerful tool to highlight new non-genetic cancer dependencies that associate to a bad outcome in vivo, paving the way to new MPM-oriented targeted strategies and prognostic tools to improve patients risk-based stratification.

Crosstalk Between Kinases, Phosphatases And Mirnas In Cancer.

Reversible phosphorylation of proteins, performed by kinases and phosphatases, is the major post translational protein modification in eukaryotic cells. This intracellular event represents a critical regulatory mechanism of several signaling pathways and can be related to a vast array of diseases, including cancer. Cancer research has produced increasing evidence that kinase and phosphatase activity can be compromised by mutations and also by miRNA silencing, performed by small non-coding and endogenously produced RNA molecules that lead to translational repression. miRNAs are believed to target about one-third of human mRNAs while a single miRNA may target about 200 transcripts simultaneously. Regulation of the phosphorylation balance by miRNAs has been a topic of intense research over the last years, spanning topics going as far as cancer aggressiveness and chemotherapy resistance. By addressing recent studies that have shown miRNA expression patterns as phenotypic signatures of cancers and how miRNA influence cellular processes such as apoptosis, cell cycle control, angiogenesis, inflammation and DNA repair, we discuss how kinases, phosphatases and miRNAs cooperatively act in cancer biology.

The Putative Leishmania Telomerase Rna (leishter) Undergoes Trans-splicing And Contains A Conserved Template Sequence.

Telomerase RNAs (TERs) are highly divergent between species, varying in size and sequence composition. Here, we identify a candidate for the telomerase RNA component of Leishmania genus, which includes species that cause leishmaniasis, a neglected tropical disease. Merging a thorough computational screening combined with RNA-seq evidence, we mapped a non-coding RNA gene localized in a syntenic locus on chromosome 25 of five Leishmania species that shares partial synteny with both Trypanosoma brucei TER locus and a putative TER candidate-containing locus of Crithidia fasciculata. Using target-driven molecular biology approaches, we detected a ∼2,100 nt transcript (LeishTER) that contains a 5' spliced leader (SL) cap, a putative 3' polyA tail and a predicted C/D box snoRNA domain. LeishTER is expressed at similar levels in the logarithmic and stationary growth phases of promastigote forms. A 5'SL capped LeishTER co-immunoprecipitated and co-localized with the telomerase protein component (TERT) in a cell cycle-dependent manner. Prediction of its secondary structure strongly suggests the existence of a bona fide single-stranded template sequence and a conserved C[U/C]GUCA motif-containing helix II, representing the template boundary element. This study paves the way for further investigations on the biogenesis of parasite TERT ribonucleoproteins (RNPs) and its role in parasite telomere biology.

Phylogenetic relationships in Solanaceae and related species based on cpDNA sequence from plastid trnE-trnT region

Intergenic spacers of chloroplast DNA (cpDNA) are very useful in phylogenetic and population genetic studies of plant species, to study their potential integration in phylogenetic analysis. The non-coding trnE-trnT intergenic spacer of cpDNA was analyzed to assess the nucleotide sequence polymorphism of 16 Solanaceae species and to estimate its ability to contribute to the resolution of phylogenetic studies of this group. Multiple alignments of DNA sequences of trnE-trnT intergenic spacer made the identification of nucleotide variability in this region possible and the phylogeny was estimated by maximum parsimony and rooted with Convolvulaceae Ipomoea batalas, the most closely related family. Besides, this intergenic spacer was tested for the phylogenetic ability to differentiate taxonomic levels. For this purpose, species from four other families were analyzed and compared with Solanaceae species. Results confirmed polymorphism in the trnE-trnT region at different taxonomic levels.

Differential expression of microRNAs in mouse pain models

Background: MicroRNAs (miRNAs) are short non-coding RNAs that inhibit translation of target genes by binding to their mRNAs. The expression of numerous brain-specific miRNAs with a high degree of temporal and spatial specificity suggests that miRNAs play an important role in gene regulation in health and disease. Here we investigate the time course gene expression profile of miR-1, -16, and -206 in mouse dorsal root ganglion (DRG), and spinal cord dorsal horn under inflammatory and neuropathic pain conditions as well as following acute noxious stimulation. Results: Quantitative real-time polymerase chain reaction analyses showed that the mature form of miR-1, -16 and -206, is expressed in DRG and the dorsal horn of the spinal cord. Moreover, CFA-induced inflammation significantly reduced miRs-1 and -16 expression in DRG whereas miR-206 was downregulated in a time dependent manner. Conversely, in the spinal dorsal horn all three miRNAs monitored were upregulated. After sciatic nerve partial ligation, miR-1 and -206 were downregulated in DRG with no change in the spinal dorsal horn. On the other hand, axotomy increases the relative expression of miR-1, -16, and 206 in a time-dependent fashion while in the dorsal horn there was a significant downregulation of miR-1. Acute noxious stimulation with capsaicin also increased the expression of miR-1 and -16 in DRG cells but, on the other hand, in the spinal dorsal horn only a high dose of capsaicin was able to downregulate miR-206 expression. Conclusions: Our results indicate that miRNAs may participate in the regulatory mechanisms of genes associated with the pathophysiology of chronic pain as well as the nociceptive processing following acute noxious stimulation. We found substantial evidence that miRNAs are differentially regulated in DRG and the dorsal horn of the spinal cord under different pain states. Therefore, miRNA expression in the nociceptive system shows not only temporal and spatial specificity but is also stimulus-dependent.

Genotypic distribution of hepatitis C among hepatitis C and HIV co-infected patients in Brazil

Information on hepatitis C virus (HCV) genotypic distribution among HIV-HCV co-infected patients is lacking in Brazil as well as other Latin American countries. The objective of this study was to evaluate the level of exposure to different risk factors associated with HCV transmission among a group of co-infected patients and to characterize the genotypic distribution of HCV in this cluster. A series of 100 HIV-HCV co-infected patients was analysed. The data to be analysed were collected from specific laboratory tests. Information was collected through a questionnaire. HCV genotyping was carried out by sequencing the 5 ` non-coding region of HCV. Chi-square and Fischer association tests or Kruskal-Wallis test were used to study the association between HCV transmission-related variables and the established genotypes. In conclusion, exposure to multiple risk factors associated with HCV transmission was common among HIV co-infected patients and an association between HCV genotype 3 and intravenous drug user was observed.

Prevalence and distribution of the GBV-C/HGV among HIV-1-infected patients under anti-retroviral therapy

Infection with GB virus C (GBV-C) or hepatitis G virus (HGV) is highly prevalent among HIV/AIDS patients. GBV-C/HGV viremia has not been associated with liver disease and seems to slow HIV disease progression. To study the GBV-C/HGV genotypes prevalence among HIV/AIDS patients and its association with HIV viral load (VL) and CD4+ lymphocyte counts. From February 2003 to February 2004, we analyzed 210 HIV-1-infected subjects who were on anti-retroviral therapy (ART). For 63 of them a PCR-nested to the non-coding 5` (5`NCR) region of the GBV-C/HGV was done, and for 49 a DNA direct sequencing was done. A phylogenetic analysis was performed by PHYLIP program. 63(30%) of the HIV-1-infected patients were co-infected with GBV-C/HGV. The phylogenetic analysis revealed the following genotypes (and respective relative frequencies): 1(10%), 2a (41%), 2b (43%), and 3 (6%). Co-infected patients presented lower HIV-1 VL and higher T CD4+ lymphocyte cells counts as compared with patients negative for GBV-C/HGV sequences (log = 4.52 vs. 4.71, p = 0.036), and T CD4+ lymphocyte counts (cells/mm(3) = 322.6 vs. 273.5, p = 0.081, respectively). T CD4+ cells counts equal to, or higher than, 200/mm(3) were significantly more common among co-infected patients than among HIV-infected-only patients (p = 0.042). The lowest T CD4+ cells counts were associated with genotype 1 and the highest with genotype 2b (p = 0.05). The GBV-C/HGV infection prevalence was 30% among HIV-1-infected subjects, and was associated with lower VL and higher CD4+ lymphocyte counts. GBV-C/HGV genotype 2b may be associated with better immunological response. Published by Elsevier B.V.

Representational Difference Analysis (RDA) reveals differential expression of conserved as well as novel genes during caste-specific development of the honey bee (Apis mellifera L.) ovary

In highly eusocial insects, such as the honey bee, Apis mellifera, the reproductive bias has become embedded in morphological caste differences. These are most expressively denoted in ovary size, with adult queens having large ovaries consisting of 150-200 ovarioles each, while workers typically have only 1-20 ovarioles per ovary. This morphological differentiation is a result of hormonal signals triggered by the diet change in the third larval instar, which eventually generate caste-specific gene expression patterns. To reveal these we produced differential gene expression libraries by Representational Difference Analysis (RDA) for queen and worker ovaries in a developmental stage when cell death is a prominent feature in the ovarioles of workers, whereas all ovarioles are maintained and extend in length in queens. In the queen library, 48% of the gene set represented homologs of known Drosophila genes, whereas in the worker ovary, the largest set (59%) were ESTs evidencing novel genes, not even computationally predicted in the honey bee genome. Differential expression was confirmed by quantitative RT-PCR for a selected gene set, denoting major differences for two queen and two worker library genes. These included two unpredicted genes located in chromosome 11 (Group11.35 and Group11.31, respectively) possibly representing long non-coding RNAs. Being candidates as modulators of ovary development, their expression and functional analysis should be a focal point for future studies. (C) 2011 Elsevier Ltd. All rights reserved.