Genome-wide detection of CNVs and their association with meat tenderness in Nelore cattle.

Brazil is one of the largest beef producers and exporters in the world with the Nelore breed representing the vast majority of Brazilian cattle (Bos taurus indicus). Despite the great adaptability of the Nelore breed to tropical climate, meat tenderness (MT) remains to be improved. Several factors including genetic composition can influence MT. In this article, we report a genome-wide analysis of copy number variation (CNV) inferred from Illumina1 High Density SNP-chip data for a Nelore population of 723 males. We detected >2,600 CNV regions (CNVRs) representing 6.5% of the genome. Comparing our results with previous studies revealed an overlap in 1400 CNVRs (>50%). A total of 1,155 CNVRs (43.6%) overlapped 2,750 genes. They were enriched for processes involving guanosine triphosphate (GTP), previously reported to influence skeletal muscle physiology and morphology. Nelore CNVRs also overlapped QTLs for MT reported in other breeds (8.9%, 236 CNVRs) and from a previous study with this population (4.1%, 109 CNVRs). Two CNVRs were also proximal to glutathione metabolism genes that were previously associated with MT. Genome-wide association study of CN state with estimated breeding values derived from meat shear force identified 6 regions, including a region on BTA3 that contains genes of the cAMP and cGMP pathway. Ten CNVRs that overlapped regions associated with MT were successfully validated by qPCR. Our results represent the first comprehensive CNV study in Bos taurus indicus cattle and identify regions in which copy number changes are potentially of importance for the MT phenotype.

Addressing tumor heterogeneity in esophageal adenocarcinoma through different molecular approaches

Several studies have shown epidemiologic, clinical, immune-histochemical and molecular differences among esophageal adenocarcinomas (EAC). Since pathogenesis and biology of this tumor are far to be well defined, our study aimed to examine intra- and inter-tumor heterogeneity and to solve crucial controversies through different molecular approaches. Target sequencing was performed for sorted cancer subpopulations from formalin embedded material obtained from 38 EACs, not treated with neoadjuvant therapy. 35 out 38 cases carried at least one somatic mutation, not present in the corresponding sorted stromal cells. 73.7% of cases carried mutations in TP53 and 10.5% in CDKN2A. Mutations in other genes occurred at lower frequency, including HNF1A, not previously associated with EAC. Sorting allowed us to isolate clones with different mutational loads and/or additional copy number amplifications, confirming the high intra-tumor heterogeneity of these cancers. In our cohort TP53 gene abnormalities correlated with a better survival (P = 0.028); conversely, loss of SMAD4 protein expression was associated with a higher recurrence rate (P = 0.015). Shifting the focus on the epigenetic characterization of EAC, miR-221 and miR-483-3p resulted upregulated from the MicroRNA Array card analysis and confirmed with further testing. The up-regulation of both miRNAs correlated with clinical outcomes, in particular with a reduced cancer-specific survival (miR483-3p P=0.0293; miR221 P=0.0059). In vitro analyses demonstrated an increase for miR-483-3p (fold-change=2.7) that appear to be inversely correlated with SMAD4 expression in FLO-1 cell-line. In conclusion, selective sorting allowed to define the real mutation status and to isolate different cancer subclones. MiRNA expression analysis revealed a significant up-regulation of miR-221 and miR-483-3p, which correlated with worst prognosis, implying that they can be considered oncogenic factors in EAC. Therefore, cell sorting technologies, coupled with next generation sequencing, and the analysis of microRNA profiles seem to be promising strategies to guide treatment and help classify cancer prognosis.

Harnessing the power of Drosophila model to obtain new insights on MYC and ABCC1 cellular functions: from neuroblastoma to autism spectrum disorder

MYCN amplification is a genetic hallmark of the childhood tumour neuroblastoma. MYCN-MAX dimers activate the expression of genes promoting cell proliferation. Moreover, MYCN seems to transcriptionally repress cell differentiation even in absence of MAX. We adopted the Drosophila eye as model to investigate the effect of high MYC to MAX expression ratio on cells. We found that dMyc overexpression in eye cell precursors inhibits cell differentiation and induces the ectopic expression of Antennapedia (the wing Hox gene). The further increase of MYC/MAX ratio results in an eye-to-wing homeotic transformation. Notably, dMyc overexpression phenotype is suppressed by low levels of transcriptional co-repressors and MYCN associates to the promoter of Deformed (the eye Hox gene) in proximity to repressive sites. Hence, we envisage that, in presence of high MYC/MAX ratio, the “free MYC” might inhibit Deformed expression, leading in turn to the ectopic expression of Antennapedia. This suggests that MYCN might reinforce its oncogenic role by affecting the physiological homeotic program. Furthermore, poor neuroblastoma outcome associates with a high level of the MRP1 protein, encoded by the ABCC1 gene and known to promote drug efflux in cancer cells. Intriguingly, this correlation persists regardless of chemotherapy and ABCC1 overexpression enhances neuroblastoma cell motility. We found that Drosophila dMRP contributes to the adhesion between the dorsal and ventral epithelia of the wing by inhibiting the function of integrin receptors, well known regulators of cell adhesion and migration. Besides, integrins play a crucial role during synaptogenesis and ABCC1 locus is included in a copy number variable region of the human genome (16p13.11) involved in neuropsychiatric diseases. Interestingly, we found that the altered dMRP/MRP1 level affects nervous system development in Drosophila embryos. These preliminary findings point out novel ABCC1 functions possibly defining ABCC1 contribution to neuroblastoma and to the pathogenicity of 16p13.11 deletion/duplication

Integrated molecular analysis of endometrial carcinoma: translation of biomarker profiles into the clinical practice

The Cancer Genome Atlas (TCGA) collaborative project identified four distinct prognostic groups of endometrial carcinoma (EC) based on molecular alterations: (i) the ultramutated subtype that encompassed POLE mutated (POLE) cases; (ii) the hypermutated subtype, characterized by MisMatch Repair deficiency (MMRd); (iii) the copy-number high subtype, with p53 abnormal/mutated features (p53abn); (iv) the copy-number low subtype, known as No Specific Molecular Profile (NSMP). Although the prognostic value of TCGA molecular classification, NSMP tumors present a wide variability in molecular alterations and biological aggressiveness. This study aims to investigate the impact of ARID1A and CTNNB1/β-catenin alterations by targeted Next-generation sequencing (NGS) and immunohistochemistry (IHC) in a consecutive series of 125 molecularly classified ECs. NGS and IHC were used to assign surrogate TCGA groups and to identify molecular alterations of multiple target genes including POLE, PTEN, ARID1A, CTNNB1, TP53. Associations with clinicopathologic parameters, molecular subtypes, and outcomes identified NSMP category as the most heterogeneous group in terms of clinicopathologic features and outcome. Integration of surrogate TCGA molecular classification with ARID1A and β-catenin analysis showed NSMP cases with ARID1A mutation characterized by the worst outcome with early recurrence, while NSMP tumors with ARID1A wild-type and β-catenin alteration had indolent clinicopathologic features and no recurrence. This study indicates how the identification of ARID1A and β-catenin alterations in EC represents a simple and effective way to characterize NSMP tumor aggressiveness and metastatic potential.

Tumor heterogeneity in Esophageal Adenocarcinoma: a genomic approach

INTRODUCTION: Esophageal adenocarcinoma (EAC) is a severe malignancy in terms of prognosis and mortality rate. Because its great genetic heterogeneity, disputes regarding classification, prevention and treatments are still unsolved. AIM: We investigated intra- and inter-EAC heterogeneity by defining EAC’s somatic mutational profile and the role of candidate microRNAs, to correlate the molecular profile of tumors to clinical outcomes and to identify biomarkers for classification. METHODS: 38 EAC cases were analyzed via high-throughput cell sorting technology combined with targeted sequencing and whole genome low-pass sequencing. Targeted sequencing of further 169 cases was performed to widen the study. miR221 and miR483-3p expression was profiled via qPCR in 112 EACs and correlation with clinical outcomes was investigated. RESULTS: 35/38 EACs carried at least one somatic mutation absent in stromal cells. TP53 was found mutated in 73.7% of cases. Selective sorting revealed tumor subclones with different mutational loads and copy number alterations, confirming the high intra-tumor heterogeneity of EAC. Mutations were in most cases at homozygous state, and we identified alterations that were missed with the whole-tumor analysis. Mutations in HNF1A gene, not previously associated with EAC, were identified in both cohorts. Higher expression of miR483-3p and miR221 was associated with poorer cancer specific survival (P=0.0293 and P=0.0059), and recurrence in the Lauren intestinal subtype (P=0.0459 and P=0.0002). Median expression levels of miRNAs were higher in patients with advanced tumor stages. The loss of SMAD4 immunoreactivity was significantly associated with poorer cancer specific survival and recurrence (P=0.0452; P=0.022 respectively). CONCLUSION: Combining selective sorting technology and next generation sequencing allowed to better define EAC inter- and intra-tumor heterogeneity. We identified HNF1A as a new mutated gene associated to EAC that could be involved in tumor progression and promising biomarkers such as SMAD4, miR221 and miR483-3p to identify patients at higher risk for more aggressive tumors.

Identification of molecular biomarkers in esophageal adenocarcinoma

Esophageal adenocarcinoma (EAC) is a severe cancer that has been on the rise in Western nations over the past few decades. It has a high mortality rate and the 5-year survival rate is only 35%–45%. EAC has been included in a group of tumors with one of the highest rates of copy number alterations (CNAs), somatic structural rearrangements, high mutation frequency, with different mutational signatures, and with epigenetic mechanisms. The vast heterogeneity of EAC mutations makes it challenging to comprehend the biology that underlies tumor onset and development, identify prognostic biomarkers, and define a molecular classification to stratify patients. The only way to resolve the current disagreements is through an exhaustive molecular analysis of EAC. We examined the genetic profile of 164 patients' esophageal adenocarcinoma samples (without chemo-radiotherapy). The included patients did not receive neoadjuvant therapies, which can change the genetic and molecular composition of the tumor. Using next-generation sequencing technologies (NGS) at high coverage, we examined a custom panel of 26 cancer-related genes. Over the entire cohort, 337 variants were found, with the TP53 gene showing the most frequent alteration (67.27%). Poorer cancer-specific survival was associated with missense mutations in the TP53 gene (Log Rank P=0.0197). We discovered HNF1alpha gene disruptive mutations in 7 cases that were also affected by other gene changes. We started to investigate its role in EAC cell lines by silencing HNF1alpha to mimic our EAC cohort and we use Seahorse technique to analyze its role in the metabolism in esophageal cell. No significant changes were found in transfected cell lines. We conclude by finding that a particular class of TP53 mutations (missense changes) adversely impacted cancer-specific survival in EAC. HNF1alpha, a new EAC-mutated gene, was found, but more research is required to fully understand its function as a tumor suppressor gene.

Phenotypic and molecular investigation of circulating tumor cells (CTCs) isolated from breast cancer patients at single cell resolution

Despite the paramount advances in cancer research, breast cancer (BC) still ranks one of the leading causes of cancer-related death worldwide. Thanks to the screening campaign started in developed countries, BC is often diagnosed at early stages (non-metastatic BC, nmBC), but disease relapse occurrence even after decades and at distant sites is not an uncommon phenomenon. Conversely, metastatic BC (mBC) is considered an incurable disease. The major perpetrators of tumor spread to secondary organs are circulating tumor cells (CTCs), a rare population of cells detectable in the peripheral blood of oncologic patients. In this study, CTCs from patients diagnosed with luminal nmBC and mBC (hormone receptor positive, Human Epidermal Growth Factor Receptor 2 (HER2) negative) were characterized at both phenotypic and molecular levels. To better understand the molecular mechanisms underlying their biology and their metastatic potential, next-generation sequencing (NGS) analyses were performed at single-cell resolution to assess copy number aberrations (CNAs), single nucleotide variants (SNVs) and gene expression profiling. The findings of this study arise hints in CTC detection, and pave the way to new application in CTC research.

The additive genetic variance after bottlenecks is affected by the number of loci involved in epistatic interactions.

We investigated the role of the number of loci coding for a neutral trait on the release of additive variance for this trait after population bottlenecks. Different bottleneck sizes and durations were tested for various matrices of genotypic values, with initial conditions covering the allele frequency space. We used three different types of matrices. First, we extended Cheverud and Routman's model by defining matrices of "pure" epistasis for three and four independent loci; second, we used genotypic values drawn randomly from uniform, normal, and exponential distributions; and third we used two models of simple metabolic pathways leading to physiological epistasis. For all these matrices of genotypic values except the dominant metabolic pathway, we find that, as the number of loci increases from two to three and four, an increase in the release of additive variance is occurring. The amount of additive variance released for a given set of genotypic values is a function of the inbreeding coefficient, independently of the size and duration of the bottleneck. The level of inbreeding necessary to achieve maximum release in additive variance increases with the number of loci. We find that additive-by-additive epistasis is the type of epistasis most easily converted into additive variance. For a wide range of models, our results show that epistasis, rather than dominance, plays a significant role in the increase of additive variance following bottlenecks.

Genetic Control of Residual Variance for Teat Number in Pigs

The genetic improvement in litter size in pigs has been substantial during the last 10-15 years. The number of teats on the sow must increase as well to meet the needs of the piglets, because each piglet needs access to its own teat. We applied a genetic heterogeneity model on teat numberin sows, and estimated medium-high heritability for teat number (0.5), but low heritability for residual variance (0.05), indicating that selection for reduced variance might have very limited effect. A numerically positive correlation (0.8) between additive genetic breeding values for mean and for variance was found, but because of the low heritability for residual variance, the variance will increase very slowly with the mean.

Within- and between-person variations as determinant factors to calculate the number of observations to estimate usual dietary intake of adolescents

This cross-sectional study describes the variation in the dietary intake of energy and macronutrients based on estimated coefficients of within- and between-person variation and intends to calculate the number of days required to evaluate the dietary intake of adolescents in Piracicaba, Brazil. The sample was made up by ninety-two adolescents aged between eleven and sixteen years. Interviews were performed to assess demographic, anthropometric and dietary variables, which were evaluated by the answers obtained through two 24-hour recalls. Descriptive analyses were performed and variances and Coefficients of variation within- and between-person were obtained by results of one-way ANOVA. The mean energy intake observed was 2,326.2 kilocalories; mean macronutrient intake came to 89.0 grams of fatty acids, 305.3 grams of carbohydrates and 82.2 grams of proteins. Coefficients of within-person variation ranged from 36% to 45%, while between-person coefficients varied from 26% to 31%. Variance ratios ranged from 1.35 (carbohydrate) to 2.62 (protein). The lowest number of observations required to correctly evaluate the usual diet calculated (r = 0.90) was six days, for carbohydrates while the highest was eleven, for proteins. Coefficients of variation were similar to those observed in international studies. We conclude that, in this study, two observations were sufficient to obtain reasonable correlations. However, six evaluations are necessary for these adolescents in order to obtain excellent correlations.

A comparative study of speed expressed by the number of throws between heavier and lighter categories in judo

Objective. - The objective of this work was to verify if there was a difference in throwing speed performance between heavier and lighter weight categories in judo. Methods and subjects. - Sixteen (16) judoists 18 +/- 3 years old, eight considered in the lightweight category (< 66 kg) and eight considered in the heavyweight (> 73 kg) category, participated in the study after signing a term of informed consent. A force-velocity test was used to determine the anaerobic power, strength, and pedal speed for each subject. In addition, three trials of Nage-komi exercise, each comprised of a set of Osoto-gari (15s), Uchi-mata (15s) and Seoi-nage (15s) throws were performed by each subject to ascertain throwing speed. Throws within the sets were intersected by one period of three minutes passive rest, while the trials were separated by one period of 10 minutes passive rest. Heart rate and the greatest number of throws within each set were measured for three trials. One-way analysis of variance (Anova) was used to compare the number of throws between the two weight categories and a ""Student"" test when the difference was significant. A correlation was used to examine the link between the different parameters. Results. - The force-velocity test did not show a significant difference in pedal speed between the two categories. However, there was a significant difference between the two categories when throwing speed was measured by the number of throws (p < 0.05) executed during the Seoi-nage (p < 0.01) and Uchi-mata (p <0.05) techniques. There was however, no significant difference between the two categories in Osoto-gari technique. Conclusion. - The throwing speed of judoists represented by the number of throws is significantly different between the two categories. The lighter category has more speed than the heavier category using the arm technique (Seoi-nage), while the heavier category has more speed using the leg technique with half turn of the attacker`s body (Uchi-mata). As a result, throwing speed is related to the type of technique used and not weight category. (C) 2007 Elsevier Masson SAS. All rights reserved.

Patchy populations in stochastic environments: Critical number of patches for persistence

We introduce a model for the dynamics of a patchy population in a stochastic environment and derive a criterion for its persistence. This criterion is based on the geometric mean (GM) through time of the spatial-arithmetic mean of growth rates. For the population to persist, the GM has to be greater than or equal to1. The GM increases with the number of patches (because the sampling error is reduced) and decreases with both the variance and the spatial covariance of growth rates. We derive analytical expressions for the minimum number of patches (and the maximum harvesting rate) required for the persistence of the population. As the magnitude of environmental fluctuations increases, the number of patches required for persistence increases, and the fraction of individuals that can be harvested decreases. The novelty of our approach is that we focus on Malthusian local population dynamics with high dispersal and strong environmental variability from year to year. Unlike previous models of patchy populations that assume an infinite number of patches, we focus specifically on the effect that the number of patches has on population persistence. Our work is therefore directly relevant to patchily distributed organisms that are restricted to a small number of habitat patches.

Effect of changing data collection parameters on statistical motor unit number estimates

The effect of number of samples and selection of data for analysis on the calculation of surface motor unit potential (SMUP) size in the statistical method of motor unit number estimates (MUNE) was determined in 10 normal subjects and 10 with amyotrophic lateral sclerosis (ALS). We recorded 500 sequential compound muscle action potentials (CMAPs) at three different stable stimulus intensities (10–50% of maximal CMAP). Estimated mean SMUP sizes were calculated using Poisson statistical assumptions from the variance of 500 sequential CMAP obtained at each stimulus intensity. The results with the 500 data points were compared with smaller subsets from the same data set. The results using a range of 50–80% of the 500 data points were compared with the full 500. The effect of restricting analysis to data between 5–20% of the CMAP and to standard deviation limits was also assessed. No differences in mean SMUP size were found with stimulus intensity or use of different ranges of data. Consistency was improved with a greater sample number. Data within 5% of CMAP size gave both increased consistency and reduced mean SMUP size in many subjects, but excluded valid responses present at that stimulus intensity. These changes were more prominent in ALS patients in whom the presence of isolated SMUP responses was a striking difference from normal subjects. Noise, spurious data, and large SMUP limited the Poisson assumptions. When these factors are considered, consistent statistical MUNE can be calculated from a continuous sequence of data points. A 2 to 2.5 SD or 10% window are reasonable methods of limiting data for analysis. Muscle Nerve 27: 320–331, 2003

Calculation of the variance in surveys of the economic climate

Public opinion surveys have become progressively incorporated into systems of official statistics. Surveys of the economic climate are usually qualitative because they collect opinions of businesspeople and/or experts about the long-term indicators described by a number of variables. In such cases the responses are expressed in ordinal numbers, that is, the respondents verbally report, for example, whether during a given trimester the sales or the new orders have increased, decreased or remained the same as in the previous trimester. These data allow to calculate the percent of respondents in the total population (results are extrapolated), who select every one of the three options. Data are often presented in the form of an index calculated as the difference between the percent of those who claim that a given variable has improved in value and of those who claim that it has deteriorated. As in any survey conducted on a sample the question of the measurement of the sample error of the results has to be addressed, since the error influences both the reliability of the results and the calculation of the sample size adequate for a desired confidence interval. The results presented here are based on data from the Survey of the Business Climate (Encuesta de Clima Empresarial) developed through the collaboration of the Statistical Institute of Catalonia (Institut d’Estadística de Catalunya) with the Chambers of Commerce (Cámaras de Comercio) of Sabadell and Terrassa.