Unnexin is a protein subunit of a large-pore channel expressed by unicellular organisms.

Cells of vertebrate and invertebrate organisms express proteins specialized in membrane channel-based cell-cell communication that are absent in unicellular organisms. We recently described the prediction of some members of the large-pore channel family in kinetoplastids, consisting of proteins called unnexins, which share several structural features with innexin and pannexin proteins. Here, we demonstrated that the unnexin1 protein (Unx1) is delivered to the cell membrane, displaying a topology consisting of four transmembrane domains with C and N termini on the cytoplasmic side and form large-pore channels that are permeable to small molecules. Low extracellular Ca2+/Mg2+ levels or extracellular alkalinization, but not mechanical stretching, increases channel activity. The Unx1 channel mediates the influx of Ca2+ and does not form intercellular dye coupling between HeLa Unx1 transfected cells. Unx1 channel function was further evidenced by its ability to mediate ionic currents when expressed in Xenopus oocytes. Downregulation of Unx1 mRNA with morpholine contains Trypanosoma cruzi invasion. Phylogenetic analysis revealed the presence of Unx1 homologs in other protozoan parasites, suggesting a conserved function for these channel parasites in other protists. Our data demonstrate that Unx1 forms large-pore membrane channels, which may serve as a diffusional pathway for ions and small molecules that are likely to be metabolic substrates or waste products, and signaling autocrine and paracrine molecules that could be involved in cell invasion. As morpholinos-induced downregulation of Unx1 reduces the infectivity of trypomastigotes, the Unx1 channels might be an attractive target for developing trypanocide drugs.

Characterizing subgenome recombination and chromosomal imbalances in banana varietal lineages.

BACKGROUND: Bananas and plantains (Musa spp.) are among the most important crops worldwide. The cultivated varieties are vegetatively propagated, so their genetic diversity is essentially fixed over time. Musa acuminata, M. balbisiana and M. schizocarpa have provided the named A, B and S subgenomes that predominantly constitute these varieties. Here we aimed to characterize intergenetic recombination and chromosomal imbalances between these A/B/S subgenomes, which often result in copy-number variants (CNVs) leading to changes in gene dosage and phenotype, in a diverse panel of bananas and plantains. This will allow us to characterize varietal lineages better and identify sources of genetic variation. METHODS: We delimited population structure and clonal lineages in a diverse panel of 188 banana and plantain accessions from the most common cultivars using admixture, principal component and phylogenetic analyses. We used new scalable alignment-based methods, Relative Averaged Alignment (RAA) and Relative Coverage, to infer subgenome composition (AA, AAB, etc.) and interspecific recombination. RESULTS: In our panel, we identified ten varietal lineages composed of somatic clones, plus three groups of tetraploid accessions. We identified chromosomal exchanges resulting in gains/losses in chromosomal segments (CNVs), particularly in AAB and ABB varieties. CONCLUSIONS: We demonstrated alignment-based RAA and Relative Coverage can identify subgenome composition and introgressions with similar results to more complex approaches based on single nucleotide polymorphism (SNP) databases. These ab initio species-agnostic methods can be used without sequencing a panel of wild ancestors to find private SNPs, or in recently diverged pools where private SNPs are uncommon. The extensive A/B/S exchanges and the variation in the length of some introgressions between lineages further support multiple foundational events of hybridization and residual backcrossing. Imbalances between A/B/S may have resulted in CNVs and gene dosage variation. Since most edible banana genomes are fixed on time, these CNVs are stable genetic variations probably associated with phenotypic variation for future genetic studies.

Alternative CNDOL Fockians for fast and accurate description of molecular exciton properties.

CNDOL is an a priori, approximate Fockian for molecular wave functions. In this study, we employ several modes of singly excited configuration interaction (CIS) to model molecular excitation properties by using four combinations of the one electron operator terms. Those options are compared to the experimental and theoretical data for a carefully selected set of molecules. The resulting excitons are represented by CIS wave functions that encompass all valence electrons in the system for each excited state energy. The Coulomb-exchange term associated to the calculated excitation energies is rationalized to evaluate theoretical exciton binding energies. This property is shown to be useful for discriminating the charge donation ability of molecular and supermolecular systems. Multielectronic 3D maps of exciton formal charges are showcased, demonstrating the applicability of these approximate wave functions for modeling properties of large molecules and clusters at nanoscales. This modeling proves useful in designing molecular photovoltaic devices. Our methodology holds potential applications in systematic evaluations of such systems and the development of fundamental artificial intelligence databases for predicting related properties.

A physiologic rise in cytoplasmic calcium ion signal increases pannexin1 channel activity via a C-terminus phosphorylation by CaMKII.

Pannexin1 (Panx1) channels are ubiquitously expressed in vertebrate cells and are widely accepted as adenosine triphosphate (ATP)-releasing membrane channels. Activation of Panx1 has been associated with phosphorylation in a specific tyrosine residue or cleavage of its C-terminal domains. In the present work, we identified a residue (S394) as a putative phosphorylation site by Ca2+/calmodulin-dependent kinase II (CaMKII). In HeLa cells transfected with rat Panx1 (rPanx1), membrane stretch (MS)-induced activation-measured by changes in DAPI uptake rate-was drastically reduced by either knockdown of Piezo1 or pharmacological inhibition of calmodulin or CaMKII. By site-directed mutagenesis we generated rPanx1S394A-EGFP (enhanced green fluorescent protein), which lost its sensitivity to MS, and rPanx1S394D-EGFP, mimicking phosphorylation, which shows high DAPI uptake rate without MS stimulation or cleavage of the C terminus. Using whole-cell patch-clamp and outside-out excised patch configurations, we found that rPanx1-EGFP and rPanx1S394D-EGFP channels showed current at all voltages between ±100 mV, similar single channel currents with outward rectification, and unitary conductance (∼30 to 70 pS). However, using cell-attached configuration we found that rPanx1S394D-EGFP channels show increased spontaneous unitary events independent of MS stimulation. In silico studies revealed that phosphorylation of S394 caused conformational changes in the selectivity filter and increased the average volume of lateral tunnels, allowing ATP to be released via these conduits and DAPI uptake directly from the channel mouth to the cytoplasmic space. These results could explain one possible mechanism for activation of rPanx1 upon increase in cytoplasmic Ca2+ signal elicited by diverse physiological conditions in which the C-terminal domain is not cleaved.

Lung Inflammatory Phenotype in Mice Deficient in Fibulin-2 and ADAMTS-12.

Interaction between extracellular matrix (ECM) components plays an important role in the regulation of cellular behavior and hence in tissue function. Consequently, characterization of new interactions within ECM opens the possibility of studying not only the functional but also the pathological consequences derived from those interactions. We have previously described the interaction between fibulin2 and ADAMTS-12 in vitro and the effects of that interaction using cellular models of cancer. Now, we generate a mouse deficient in both ECM components and evaluate functional consequences of their absence using different cancer and inflammation murine models. The main findings indicate that mice deficient in both fibulin2 and ADAMTS12 markedly increase the development of lung tumors following intraperitoneal urethane injections. Moreover, inflammatory phenotype is exacerbated in the lung after LPS treatment as can be inferred from the accumulation of active immune cells in lung parenchyma. Overall, our results suggest that protective effects in cancer or inflammation shown by fibulin2 and ADAMTS12 as interactive partners in vitro are also shown in a more realistic in vivo context.

Axonal and Glial PIEZO1 and PIEZO2 Immunoreactivity in Human Clitoral Krause's Corpuscles.

Krause's corpuscles are typical of cutaneous mucous epithelia, like the lip vermillion or the glans clitoridis, and are associated with rapidly adapting low-threshold mechanoreceptors involved in gentle touch or vibration. PIEZO1 and PIEZO2 are transmembrane mechano-gated proteins that form a part of the cationic ion channels required for mechanosensitivity in mammalian cells. They are involved in somatosensitivity, especially in the different qualities of touch, but also in pain and proprioception. In the present study, immunohistochemistry and immunofluorescence were used to analyze the occurrence and cellular location of PIEZO1 and PIEZO2 in human clitoral Krause's corpuscles. Both PIEZO1 and PIEZO2 were detected in Krause's corpuscles in both the axon and the terminal glial cells. The presence of PIEZOs in the terminal glial cells of Kraus's corpuscles is reported here for the first time. Based on the distribution of PIEZO1 and PIEZO2, it may be assumed they could be involved in mechanical stimuli, sexual behavior, and sexual pleasure.

Developing customized fuel models for shrub and bracken communities in Galicia (NW Spain).

Geospatial fire behaviour and fire hazard simulators, fire effects models and smoke emission software commonly use standard fuel models in order to simplify data collection and the inclusion of complex fuel scenarios. These fuel models are often mapped using remotely sensed data. However, given the great complexity of fuelbeds, with properties that vary widely in both time and space, the use of these standard fuel models can greatly limit accurate fuel mapping. This affects fuel hazard assessment, fuel reduction treatment plans, fire management decision-making and evaluation of the environmental impact of wildfire. In this study, we developed unique customized fire behaviour fuel models for shrub and bracken communities, by using k-medoids clustering analysis based on both fuel structural characteristics and potential fire behaviour. We used an original database of 722 destructive sample plots in nine different shrub and bracken communities covering the entire distribution area in Galicia (NW Spain), one of the regions in Europe most affected by forest fires. Measurements of cover, height and fuel fractions loads differentiated by size and vegetative state (live or dead) were used to estimate the potential rate of fire spread with five different models including fireline intensity, heat per unit area and the flame length for each sampling site and considering extreme environmental conditions. The optimal number of clusters was established by combining practical knowledge about the shrubland communities under study and their associated fire behaviour, with maximization of the mean value of the silhouette variable and minimization of the within-cluster sum of squares. The structural characteristics of the medoids derived from the analysis were associated with each of the proposed customized fuel models. Finally, a simple dichotomous classification based only on shrub height was developed to enable construction of spatially explicit fuel model maps based on remotely sensed data. Thus, the methodology applied allows generation of a more realistic representation of fuel distribution in the landscape, based on fuel structure measurements of natural regional ecosystems rather than on the use of standard models. We believe that the proposed methodology is generally applicable to communities composed of other shrub and fern species in different biogeographical regions.

Analysis of Myocarditis Among 252 Million mRNA-1273 Recipients Worldwide.

BACKGROUND: Growing evidence indicates a causal relationship between SARS-CoV-2 infection and myocarditis. Post-authorization safety data have also identified myocarditis as a rare safety event following mRNA COVID-19 vaccination, particularly among adolescent and young-adult males after dose 2. We further evaluated the potential risk by querying the Moderna global safety database for myocarditis/myopericarditis reports among mRNA-1273 recipients worldwide. METHODS: Myocarditis/myopericarditis reports from 18 December 2020 to 15 February 2022 were reviewed and classified. The reported rate after any known mRNA-1273 dose was calculated according to age and sex, then compared with a population-based incidence rate to calculate observed-to-expected rate ratios (RRs). RESULTS: During the study period, 3017 myocarditis/myopericarditis cases among 252 million mRNA-1273 recipients who received at least 1 dose were reported to the Moderna global safety database. The overall reporting rate was 9.23 per 100 000 person-years, which was similar to the expected reference rate (9.0 cases per 100 000 person-years; RR [95% confidence interval (CI)], 1.03 [.97-1.08]). When stratified by sex and age, observed rates were highest for males aged <40 years, particularly those 18-24 years (53.76 per 100 000 person-years), which was higher than expected (RR [95% CI], 3.10 [2.68-3.58]). When considering only cases occurring within 7 days of a known dose, the observed rate was highest for males aged 18-24 years after dose 2 (4.23 per 100 000 doses administered). CONCLUSIONS: Myocarditis/myopericarditis rates were not higher than expected for the overall population of mRNA-1273 recipients but were higher than expected in males aged 18-24 years, with most cases occurring 7 days after dose 2.

Computational approaches to amino acid side-chain conformation using combined NMR theoretical and experimental results: leucine-67 in Desulfovibrio vulgaris flavodoxin.

In this paper, we show that the combination of NMR theoretical and experimental results can help to solve the molecular structure of peptides, here it is used as an example the residue Leucine-67 in Desulfovibrio vulgaris flavodoxin. We apply a computational protocol based on the leucine amino acid dipeptide, which, using calculated and experimental spin-spin coupling constants, allows us to obtain the conformation of the amino acid side chain. Calculated results show that the best agreement is obtained when three conformers around the lateral chain angle $\chi _1$ are considered or when the dynamic effect in the torsional angles is included. The population of each structure is estimated and analyzed according to the correlation between those two approaches. Independently of the approach, the estimated $\chi _1$ angle in solution is close to the staggered value of -60$^\circ $ and deviates significantly from the average x-ray angle of -90$^\circ $.

The Challenge of ab Initio Calculations in Small Neon Clusters.

Weakly bound neon dimer, trimer and tetramers are studied at HF and CCSD(T) levels using Dunning, ANO and SIGMA-s basis sets. Their ground-state binding energies are studied along with some structural properties. SIGMA-s basis sets have been developed explicitly for this issue but in a manner that can be readily applied to other atoms for the study of larger weakly bound systems. The difficulties for attaining accurate results on these systems are assessed by the computation of total, atomization and correlation energies, as well as equilibrium distances, with several basis sets of increasing size, ranging from non-augmented to double-augmented versions. Extrapolations are proposed to predict stabilization energies and the results are compared with previously published data.

Complex polyploid and hybrid species in an apomictic and sexual tropical forage grass group: genomic composition and evolution in Urochloa (Brachiaria) species.

BACKGROUND AND AIMS: Diploid and polyploid Urochloa (including Brachiaria, Panicum and Megathyrsus species) C4 tropical forage grasses originating from Africa are important for food security and the environment, often being planted in marginal lands worldwide. We aimed to characterize the nature of their genomes, the repetitive DNA and the genome composition of polyploids, leading to a model of the evolutionary pathways within the group including many apomictic species. METHODS: Some 362 forage grass accessions from international germplasm collections were studied, and ploidy was determined using an optimized flow cytometry method. Whole-genome survey sequencing and molecular cytogenetic analysis were used to identify chromosomes and genomes in Urochloa accessions belonging to the 'brizantha' and 'humidicola' agamic complexes and U. maxima. KEY RESULTS: Genome structures are complex and variable, with multiple ploidies and genome compositions within the species, and no clear geographical patterns. Sequence analysis of nine diploid and polyploid accessions enabled identification of abundant genome-specific repetitive DNA motifs. In situ hybridization with a combination of repetitive DNA and genomic DNA probes identified evolutionary divergence and allowed us to discriminate the different genomes present in polyploids. CONCLUSIONS: We suggest a new coherent nomenclature for the genomes present. We develop a model of evolution at the whole-genome level in diploid and polyploid accessions showing processes of grass evolution. We support the retention of narrow species concepts for Urochloa brizantha, U. decumbens and U. ruziziensis, and do not consider diploids and polyploids of single species as cytotypes. The results and model will be valuable in making rational choices of parents for new hybrids, assist in use of the germplasm for breeding and selection of Urochloa with improved sustainability and agronomic potential, and assist in measuring and conserving biodiversity in grasslands.

Serially improved GTOs for molecular applications (SIGMA): Basis sets from H to Ne.

A new approach for generating Gaussian basis sets is reported and tested for atoms from H to Ne. The basis sets thus calculated, named SIGMA basis sets, range from DZ to QZ sizes and have the same composition per shell as Dunning basis sets but with different treatment of the contractions. The standard SIGMA basis sets and their augmented versions have proven to be very suitable for providing good results in atomic and molecular calculations. The performance of the new basis sets is analyzed in terms of total, correlation, and atomization energies, equilibrium distances, and vibrational frequencies in several molecules, and the results are compared at several computational levels with those obtained with the corresponding Dunning and other basis sets.

Acid-Sensing Ion Channels' Immunoreactivity in Nerve Profiles and Glomus Cells of the Human Carotid Body.

The carotid body is a major peripheral chemoreceptor that senses changes in arterial blood oxygen, carbon dioxide, and pH, which is important for the regulation of breathing and cardiovascular function. The mechanisms by which the carotid body senses O2 and CO2 are well known; conversely, the mechanisms by which it senses pH variations are almost unknown. Here, we used immunohistochemistry to investigate how the human carotid body contributes to the detection of acidosis, analyzing whether it expresses acid-sensing ion channels (ASICs) and determining whether these channels are in the chemosensory glomic cells or in the afferent nerves. In ASIC1, ASIC2, and ASIC3, and to a much lesser extent ASIC4, immunoreactivity was detected in subpopulations of type I glomus cells, as well as in the nerves of the carotid body. In addition, immunoreactivity was found for all ASIC subunits in the neurons of the petrosal and superior cervical sympathetic ganglia, where afferent and efferent neurons are located, respectively, innervating the carotid body. This study reports for the first time the occurrence of ASIC proteins in the human carotid body, demonstrating that they are present in glomus chemosensory cells (ASIC1 < ASIC2 > ASIC3 > ASIC4) and nerves, presumably in both the afferent and efferent neurons supplying the organ. These results suggest that the detection of acidosis by the carotid body can be mediated via the ASIC ion channels present in the type I glomus cells or directly via sensory nerve fibers.

Potential Neuroprotective Role of Calretinin-N18 and Calbindin-D28k in the Retina of Adult Zebrafish Exposed to Different Wavelength Lights.

The incidence rates of light-induced retinopathies have increased significantly in the last decades because of continuous exposure to light from different electronic devices. Recent studies showed that exposure to blue light had been related to the pathogenesis of light-induced retinopathies. However, the pathophysiological mechanisms underlying changes induced by light exposure are not fully known yet. In the present study, the effects of exposure to light at different wavelengths with emission peaks in the blue light range (400-500 nm) on the localization of Calretinin-N18 (CaR-N18) and Calbindin-D28K (CaB-D28K) in adult zebrafish retina are studied using double immunofluorescence with confocal laser microscopy. CaB-D28K and CaR-N18 are two homologous cytosolic calcium-binding proteins (CaBPs) implicated in essential process regulation in central and peripheral nervous systems. CaB-D28K and CaR-N18 distributions are investigated to elucidate their potential role in maintaining retinal homeostasis under distinct light conditions and darkness. The results showed that light influences CaB-D28K and CaR-N18 distribution in the retina of adult zebrafish, suggesting that these CaBPs could be involved in the pathophysiology of retinal damage induced by the short-wavelength visible light spectrum.

Development and Pharmacokinetic Evaluation of Two Parenteral Formulations of Albendazole Using Prodrug and Cosolvent Approaches.

Albendazole is a broad-spectrum anthelmintic drug used for parasitic infections. In addition, due to its mechanism of action, it has been studied as an anticancer agent. However, poor and highly variable bioavailability are limiting factors for its use in systemic illnesses. The present study aimed to develop two parenteral formulations of albendazole and to compare its pharmacokinetic profile with the conventional oral administration. Parenteral formulations were developed using two different approaches: a phosphonooxymethylated prodrug and cosolvents. For the albendazole prodrug, once synthetized, its solubility and hydrolysis with alkaline phosphatase were evaluated. A factorial design of experiments was used for the cosolvent formulation. Stability and hemolytic activity were assessed. A pharmacokinetic study was performed on New Zealand rabbits. Both formulations were administered intravenously, and the prodrug was also administered intramuscularly. Results were compared with those obtained after the oral administration of albendazole. A 20,000-fold and 6000-fold increase in albendazole solubility was found with the prodrug and cosolvent formulations, respectively. Both parenteral formulations displayed higher albendazole plasma concentrations for the first 2 h compared with oral administration, even when the oral dose was doubled. The absolute bioavailability of oral albendazole was 15.5% while for the intramuscular administration of the prodrug was 102.6%. Both parenteral formulations showed a significant decrease in the formation of albendazole sulfoxide (ANOVA p<0.05) and allowed greater exposure to albendazole. Albendazole cosolvent parenteral formulation could be a promising option in systemic illnesses considering its ease of preparation and superb pharmacokinetic performance.

Unnexins: Homologs of innexin proteins in Trypanosomatidae parasites.

Large-pore channels, including those formed by connexin, pannexin, innexin proteins, are part of a broad family of plasma membrane channels found in vertebrates and invertebrates, which share topology features. Despite their relevance in parasitic diseases such as Chagas and malaria, it was unknown whether these large-pore channels are present in unicellular organisms. We identified 14 putative proteins in Trypanosomatidae parasites as presumptive homologs of innexin proteins. All proteins possess the canonical motif of the innexin family, a pentapeptide YYQWV, and 10 of them share a classical membrane topology of large-pore channels. A sequence similarity network analysis confirmed their closeness to innexin proteins. A bioinformatic model showed that a homolog of Trypanosoma cruzi (T. cruzi) could presumptively form a stable octamer channel with a highly positive electrostatic potential in the internal cavities and extracellular entrance due to the notable predominance of residues such as Arg or Lys. In vitro dye uptake assays showed that divalent cations-free solution increases YO-PRO-1 uptake and hyperosmotic stress increases DAPI uptake in epimastigotes of T. cruzi. Those effects were sensitive to probenecid. Furthermore, probenecid reduced the proliferation and transformation of T. cruzi. Moreover, probenecid or carbenoxolone increased the parasite sensitivity to antiparasitic drugs commonly used in therapy against Chagas. Our study suggests the existence of innexin homologs in unicellular organisms, which could be protein subunits of new large-pore channels in unicellular organisms.

Lianas decelerate tropical forest thinning during succession.

The well-established pattern of forest thinning during succession predicts an increase in mean tree biomass with decreasing tree density. The forest thinning pattern is commonly assumed to be driven solely by tree-tree competition. The presence of non-tree competitors could alter thinning trajectories, thus altering the rate of forest succession and carbon uptake. We used a large-scale liana removal experiment over 7 years in a 60- to 70-year-old Panamanian forest to test the hypothesis that lianas reduce the rate of forest thinning during succession. We found that lianas slowed forest thinning by reducing tree growth, not by altering tree recruitment or mortality. Without lianas, trees grew and presumably competed more, ultimately reducing tree density while increasing mean tree biomass. Our findings challenge the assumption that forest thinning is driven solely by tree-tree interactions; instead, they demonstrate that competition from other growth forms, such as lianas, slow forest thinning and ultimately delay forest succession.

Vegetative phenologies of lianas and trees in two Neotropical forests with contrasting rainfall regimes.

Among tropical forests, lianas are predicted to have a growth advantage over trees during seasonal drought, with substantial implications for tree and forest dynamics. We tested the hypotheses that lianas maintain higher water status than trees during seasonal drought and that lianas maximize leaf cover to match high, dry-season light conditions, while trees are more limited by moisture availability during the dry season. We monitored the seasonal dynamics of predawn and midday leaf water potentials and leaf phenology for branches of 16 liana and 16 tree species in the canopies of two lowland tropical forests with contrasting rainfall regimes in Panama. In a wet, weakly seasonal forest, lianas maintained higher water balance than trees and maximized their leaf cover during dry-season conditions, when light availability was high, while trees experienced drought stress. In a drier, strongly seasonal forest, lianas and trees displayed similar dry season reductions in leaf cover following strong decreases in soil water availability. Greater soil moisture availability and a higher capacity to maintain water status allow lianas to maintain the turgor potentials that are critical for plant growth in a wet and weakly seasonal forest but not in a dry and strongly seasonal forest.

Sensory innervation of the human palmar aponeurosis in healthy individuals and patients with palmar fibromatosis.

The human palmar aponeurosis is involved in hand proprioception, and it contains different sensory corpuscle morphotypes that serve this role. In palmar fibromatosis (classically referred to as Dupuytren's disease), the palmar aponeurosis undergoes fibrous structural changes that, presumably, also affect the nervous system, causing altered perception. We analysed the various sensory nerve formation morphotypes in the palmar aponeuroses of healthy subjects and patients with palmar fibromatosis. To do this, we used immunohistochemistry for corpuscular constituents and the putative mechanoproteins PIEZO2 and acid-sensing ion channel 2. Free nerve endings and Golgi-Mazzoni, Ruffini, paciniform and Pacinian corpuscles were identified in both the healthy and the pathological conditions. The densities of the free nerve endings and Golgi-Mazzoni corpuscles were slightly increased in the pathological tissues. Furthermore, the Pacinian corpuscles were enlarged and displayed an altered shape. Finally, there was also morphological and immunohistochemical evidence of occasional denervation of the Pacinian corpuscles, although no increase in their number was observed. Both PIEZO2 and acid-sensing ion channel 2 were absent from the altered corpuscles. These results indicate that the human palmar aponeurosis is richly innervated, and the free nerve endings and sensory corpuscles within the palmar aponeurosis undergo quantitative and qualitative changes in patients with palmar fibromatosis, which may explain the sensory alterations occasionally reported for this pathology.

Growth hormone treatment does not to lead to insulin resistance nor excessive rise in IGF-1 levels, while improving height in patients small for gestational age A long-term observational study.

OBJECTIVE: In children born small for gestational age (SGA), the relationship between growth hormone (GH) treatment and insulin resistance (IR) has only been investigated for a short period, necessitating a longer observation period. This study aimed to evaluate the long-term (10 years) effect of GH to SGA-children on IR and safety during treatment. DESIGN: This was a multicenter observational study. PATIENTS: SGA-children who received GH treatment in Spain (stratified by Tanner-stage and age at GH onset [two groups: ≤6 years old or >6 years old]). MEASUREMENTS: The analysed variables (yearly measures) included auxologic, metabolic (insulin-like growth factor-1 (IGF-1), height velocity [HV], weight and homeostatic model assessment-IR [HOMA-IR]) and safety data. Data were collected prospectively (since the study approval: 2007) and retrospectively (since the initiation of GH treatment: 2005-2007). RESULTS: A total of 389 SGA children (369 Tanner-I) were recruited from 27 centres. The mean age (standard deviation) of the children at GH treatment onset was 7.2 (2.8) years old. IGF-1 (standard deviation score [SDS]) and HOMA-IR values tended to increase until the sixth year of GH-treatment, with significant differences being observed only during the first year, while these remained stable in the later years (within normal ranges). Height (SDS) increased significantly (basal: -3.0; tenth year: -1.13), and the maximum HV (SDS) occurred during the first year (2.75 ± 2.39). CONCLUSIONS: HOMA-IR values increased significantly in SGA-children during the first year of GH-treatment, remained stable and were within normal ranges in all cases. Our 10-year data suggests that long-term GH treatment does not promote IR and is well-tolerated, safe and effective.