Search

Many treatments against breast cancer decrease the level of estrogen in blood, resulting in bone loss, osteoporosis and fragility fractures in breast cancer patients. This retrospective study aimed to evaluate a novel opportunistic screening for cancer treatment-induced bone loss (CTIBL) in breast cancer patients using CT radiomics. Between 2011 and 2021, a total of 412 female breast cancer patients who received treatment and were followed up in our institution, had post-treatment dual-energy X-ray absorptiometry (DXA) examination of the lumbar vertebrae and had post-treatment chest CT scan that encompassed the L1 vertebra, were included in this study. Results indicated that the T-score of L1 vertebra had a strongly positive correlation with the average T-score of L1-L4 vertebrae derived from DXA (r = 0.91, p < 0.05). On multivariable analysis, four clinical variables (age, body weight, menopause status, aromatase inhibitor exposure duration) and three radiomic features extracted from the region of interest of L1 vertebra (original&#95;firstorder&#95;RootMeanSquared, wavelet.HH&#95;glcm&#95;InverseVariance, and wavelet.LL&#95;glcm&#95;MCC) were selected for building predictive models of L1 T-score and bone health. The predictive model combining clinical and radiomic features showed the greatest adjusted R 2 value (0.557), sensitivity (83.6%), specificity (74.2%) and total accuracy (79.4%) compared to models that relied solely on clinical data, radiomic features, or Hounsfield units. In conclusion, the clinical-radiomic predictive model may be used as an opportunistic screening tool for early identification of breast cancer survivors at high risk of CTIBL based on non-contrast CT images of the L1 vertebra, thereby facilitating early intervention for osteoporosis., (© 2024. The Author(s).)

Wolfberry, also known as goji berry or Lycium barbarum , is a highly valued fruit with significant health benefits and nutritional value. For more efficient and comprehensive usage of published L. barbarum genomic data, we established the Wolfberry database. The utility of the Wolfberry Genome Database (WGDB) is highlighted through the Genome browser, which enables the user to explore the L. barbarum genome, browse specific chromosomes, and access gene sequences. Gene annotation features provide comprehensive information about gene functions, locations, expression profiles, pathway involvement, protein domains, and regulatory transcription factors. The transcriptome feature allows the user to explore gene expression patterns using transcripts per kilobase million (TPM) and fragments per kilobase per million mapped reads (FPKM) metrics. The Metabolism pathway page provides insights into metabolic pathways and the involvement of the selected genes. In addition to the database content, we also introduce six analysis tools developed for the WGDB. These tools offer functionalities for gene function prediction, nucleotide and amino acid BLAST analysis, protein domain analysis, GO annotation, and gene expression pattern analysis. The WGDB is freely accessible at https://cosbi7.ee.ncku.edu.tw/Wolfberry/. Overall, WGDB serves as a valuable resource for researchers interested in the genomics and transcriptomics of L. barbarum . Its user-friendly web interface and comprehensive data facilitate the exploration of gene functions, regulatory mechanisms, and metabolic pathways, ultimately contributing to a deeper understanding of wolfberry and its potential applications in agronomy and nutrition., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Cao, Chen, Li, Li, Lin, Lee, Hsieh, Hsiao, Fan, Luo, Zhao, Yin, An, Shi, Chow, Chang, Huang, Chang, Liu, Wu and Tsai.)

Prior research has shown that the deconvolution of cell-free RNA can uncover the tissue origin. The conventional deconvolution approaches rely on constructing a reference tissue-specific gene panel, which cannot capture the inherent variation present in actual data. To address this, we have developed a novel method that utilizes a neural network framework to leverage the entire training dataset. Our approach involved training a model that incorporated 15 distinct tissue types. Through one semi-independent and two complete independent validations, including deconvolution using a semi in silico dataset, deconvolution with a custom normal tissue mixture RNA-seq data, and deconvolution of longitudinal circulating tumor cell RNA-seq (ctcRNA) data from a cancer patient with metastatic tumors, we demonstrate the efficacy and advantages of the deep-learning approach which were exerted by effectively capturing the inherent variability present in the dataset, thus leading to enhanced accuracy. Sensitivity analyses reveal that neural network models are less susceptible to the presence of missing data, making them more suitable for real-world applications. Moreover, by leveraging the concept of organotropism, we applied our approach to trace the migration of circulating tumor cell-derived RNA (ctcRNA) in a cancer patient with metastatic tumors, thereby highlighting the potential clinical significance of early detection of cancer metastasis., (© 2023. The Author(s).)

Motivation: As an important player in transcriptome regulation, microRNAs may effectively diffuse somatic mutation impacts to broad cellular processes and ultimately manifest disease and dictate prognosis. Previous studies that tried to correlate mutation with gene expression dysregulation neglected to adjust for the disparate multitudes of false positives associated with unequal sample sizes and uneven class balancing scenarios., Results: To properly address this issue, we developed a statistical framework to rigorously assess the extent of mutation impact on microRNAs in relation to a permutation-based null distribution of a matching sample structure. Carrying out the framework in a pan-cancer study, we ascertained 9008 protein-coding genes with statistically significant mutation impacts on miRNAs. Of these, the collective miRNA expression for 83 genes showed significant prognostic power in nine cancer types. For example, in lower-grade glioma, 10 genes' mutations broadly impacted miRNAs, all of which showed prognostic value with the corresponding miRNA expression. Our framework was further validated with functional analysis and augmented with rich features including the ability to analyze miRNA isoforms; aggregative prognostic analysis; advanced annotations such as mutation type, regulator alteration, somatic motif, and disease association; and instructive visualization such as mutation OncoPrint, Ideogram, and interactive mRNA-miRNA network., Availability and Implementation: The data underlying this article are available in MutMix, at http://innovebioinfo.com/Database/TmiEx/MutMix.php., (© The Author(s) 2023. Published by Oxford University Press.)

Background: The metabolome is the best representation of cancer phenotypes. Gene expression can be considered a confounding covariate affecting metabolite levels. Data integration across metabolomics and genomics to establish the biological relevance of cancer metabolism is challenging. This study aimed to eliminate the confounding effect of metabolic gene expression to reflect actual metabolite levels in microsatellite instability (MSI) cancers., Methods: In this study, we propose a new strategy using covariate-adjusted tensor classification in high dimensions (CATCH) models to integrate metabolite and metabolic gene expression data to classify MSI and microsatellite stability (MSS) cancers. We used datasets from the Cancer Cell Line Encyclopedia (CCLE) phase II project and treated metabolomic data as tensor predictors and data on gene expression of metabolic enzymes as confounding covariates., Results: The CATCH model performed well, with high accuracy (0.82), sensitivity (0.66), specificity (0.88), precision (0.65), and F1 score (0.65). Seven metabolite features adjusted for metabolic gene expression, namely, 3-phosphoglycerate, 6-phosphogluconate, cholesterol ester, lysophosphatidylethanolamine (LPE), phosphatidylcholine, reduced glutathione, and sarcosine, were found in MSI cancers. Only one metabolite, Hippurate, was present in MSS cancers. The gene expression of phosphofructokinase 1 (PFKP), which is involved in the glycolytic pathway, was related to 3-phosphoglycerate. ALDH4A1 and GPT2 were associated with sarcosine. LPE was associated with the expression of CHPT1, which is involved in lipid metabolism. The glycolysis, nucleotide, glutamate, and lipid metabolic pathways were enriched in MSI cancers., Conclusions: We propose an effective CATCH model for predicting MSI cancer status. By controlling the confounding effect of metabolic gene expression, we identified cancer metabolic biomarkers and therapeutic targets. In addition, we provided the possible biology and genetics of MSI cancer metabolism., (© 2023. The Author(s).)

Introduction: Cystic periventricular leukomalacia (PVL) is the most common white matter injury and a common cause of cerebral palsy in preterm infants. Postnatal epilepsy may occur after cystic PVL, but their causal relationship remains uncertain. Our aim was to validate the contribution of cystic PVL to postnatal epilepsy in very preterm infants and demonstrate their seizure characteristics., Methods: This prospective cohort study enrolled 1,342 preterm infants (birth weight &lt;1,500 g and gestational age &lt;32 weeks) from 2003 to 2015. Cystic PVL was diagnosed by serial cerebral ultrasound, and other comorbidities were recorded during hospitalization. Neurological developments and consequences, including epilepsy, were serially accessed until the age of 5., Results: A total of 976 preterm infants completed a 5-year neurological follow-up; 47 (4.8%) had cystic PVL. Preterm infants with cystic PVL were commonly associated with other comorbidities, including necrotizing enterocolitis stage III, neonatal seizures, and intraventricular hemorrhage during hospitalization. At age 5, 14 of the 47 (29.8%) preterm infants with cystic PVL had postnatal epilepsy. After adjusting for gender, gestational age, and three common comorbidities, cystic PVL was an independent risk factor for postnatal epilepsy (adjust OR: 16.2; 95% CI: 6.8-38.4; p &lt; 0.001). Postnatal epilepsy after cystic PVL was commonly the generalized type (13 of 14, 92.9%), not intractable and most occurred after 1 year of age., Discussion/conclusion: Cystic PVL would independently lead to postnatal epilepsy. Preterm infants with cystic PVL are at risk of postnatal epilepsy after age 1 in addition to cerebral palsy., (© 2023 S. Karger AG, Basel.)

Containing the largest number of species, the orchid family provides not only materials for studying plant evolution and environmental adaptation, but economically and culturally important ornamental plants for human society. Previously, we collected genome and transcriptome information of Dendrobium catenatum, Phalaenopsis equestris, and Apostasia shenzhenica which belong to two different subfamilies of Orchidaceae, and developed user-friendly tools to explore the orchid genetic sequences in the OrchidBase 4.0. The OrchidBase 4.0 offers the opportunity for plant science community to compare orchid genomes and transcriptomes and retrieve orchid sequences for further study.In the year 2022, two whole-genome sequences of Orchidoideae species, Platanthera zijinensis and Platanthera guangdongensis, were de novo sequenced, assembled and analyzed. In addition, systemic transcriptomes from these two species were also established. Therefore, we included these datasets to develop the new version of OrchidBase 5.0. In addition, three new functions including synteny, gene order, and miRNA information were also developed for orchid genome comparisons and miRNA characterization.OrchidBase 5.0 extended the genetic information to three orchid subfamilies (including five orchid species) and provided new tools for orchid researchers to analyze orchid genomes and transcriptomes. The online resources can be accessed at https://cosbi.ee.ncku.edu.tw/orchidbase5/., (© 2022. The Author(s).)

Background: Improper endotracheal tube (ETT) positioning is frequently observed and potentially hazardous in the intensive care unit. The authors developed a deep learning-based automatic detection algorithm detecting the ETT tip and carina on portable supine chest radiographs to measure the ETT-carina distance. This study investigated the hypothesis that the algorithm might be more accurate than frontline critical care clinicians in ETT tip detection, carina detection, and ETT-carina distance measurement., Methods: A deep learning-based automatic detection algorithm was developed using 1,842 portable supine chest radiographs of 1,842 adult intubated patients, where two board-certified intensivists worked together to annotate the distal ETT end and tracheal bifurcation. The performance of the deep learning-based algorithm was assessed in 4-fold cross-validation (1,842 radiographs), external validation (216 radiographs), and an observer performance test (462 radiographs) involving 11 critical care clinicians. The performance metrics included the errors from the ground truth in ETT tip detection, carina detection, and ETT-carina distance measurement., Results: During 4-fold cross-validation and external validation, the median errors (interquartile range) of the algorithm in ETT-carina distance measurement were 3.9 (1.8 to 7.1) mm and 4.2 (1.7 to 7.8) mm, respectively. During the observer performance test, the median errors (interquartile range) of the algorithm were 2.6 (1.6 to 4.8) mm, 3.6 (2.1 to 5.9) mm, and 4.0 (1.7 to 7.2) mm in ETT tip detection, carina detection, and ETT-carina distance measurement, significantly superior to that of 6, 10, and 7 clinicians (all P < 0.05), respectively. The algorithm outperformed 7, 3, and 0, 9, 6, and 4, and 5, 5, and 3 clinicians (all P < 0.005) regarding the proportions of chest radiographs within 5 mm, 10 mm, and 15 mm error in ETT tip detection, carina detection, and ETT-carina distance measurement, respectively. No clinician was significantly more accurate than the algorithm in any comparison., Conclusions: A deep learning-based algorithm can match or even outperform frontline critical care clinicians in ETT tip detection, carina detection, and ETT-carina distance measurement., (Copyright © 2022, the American Society of Anesthesiologists. All Rights Reserved.)

Mutation density patterns reveal unique biological properties of specific genomic regions and shed light on the mechanisms of carcinogenesis. Although previous studies reported insightful mutation density patterns associated with certain genomic regions such as transcription start sites and DNA replication origins, a tool that can systematically investigate mutational spatial patterns is still lacking. Thus, we developed MutDens, a bioinformatic tool for comprehensive analysis of mutation density patterns around genomic features, namely, genomic positions, in humans and model species. By scanning the bidirectional vicinity regions of given positions, MutDens systematically characterizes the mutation density for single-base substitution mutational classes after adjusting for total mutation burden and local nucleotide proportion. Analysis results using MutDens not only verified the previously reported transcriptional strand bias around transcription start sites and replicative strand bias around DNA replication origins, but also identified novel mutation density patterns around other genomics features, such as enhancers and retrotransposon insertion polymorphism sites. To our knowledge, MutDens is the first tool that systematically calculates, examines, and compares mutation density patterns, thus providing a valuable avenue for investigating the mutational landscapes associated with important genomic features., (© 2022 Yu et al.; Published by Cold Spring Harbor Laboratory Press.)

Patients with thrombocytopenia (platelet count <150 × 10 3 /μL) often develop pulmonary hemorrhage (PH) after Stenotrophomonas maltophilia (SM) respiratory infection, resulting in a high respiratory failure rate and increased mortality. Developing an efficient method for early prediction of PH in these patients may improve survival. This study aimed to evaluate risk factors in PH and to develop an index measuring serial platelet deficit to predict PH in patients with SM respiratory infection. Data of patients with SM respiratory infection and thrombocytopenia treated in a tertiary university hospital during 2018-2020 were retrospectively retrieved from electronic medical records and analyzed. SM respiratory infection was defined as SM isolated from sputum, endotracheal suction, or bronchial alveolar lavage plus acute respiratory symptoms. Between PH and non-PH groups, clinical characteristics and laboratory parameters were collected and compared. The newly developed platelet dissimilarity index (d-index) was calculated by accumulating differences between the actual and the lowest normal level of the platelet count in each patient at different time points. Within 1,039 patients with positive SM culture, 437 cases matched the criteria and were analyzed. A total of 125 (28.6%) patients developed PH and 312 (71.4%) did not. The patients with PH had increased prothrombin time/international normalized ratio (PT/INR), lower platelet count, and higher platelet d-index. Multivariate analysis revealed that extreme thrombocytopenia (platelet count <50 × 10 3 /μL) is a common independent risk factor in PH and mortality. The performance of platelet deficit and d-index varied between patients with different comorbidities. Performance of platelet deficit to predict PH is better in patients with hematology/oncology or liver disease (area under curve, 0.705-0.757), while d-index is better in patients with sepsis/treatment and various other groups (0.711-0.816). Prolonged and extreme thrombocytopenia is a determinant risk factor in PH in patients with SM respiratory infection. Given the complexity of causes of thrombocytopenia and associated comorbidities, different strategies should be applied when assessing the risk for PH., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Wu, Wu, Cheng, Chen, Chen, Li and Shen.)

To evaluate whether adjusted computed tomography (CT) scan image-based radiomics combined with immune genomic expression can achieve accurate stratification of cancer recurrence and identify potential therapeutic targets in stage III colorectal cancer (CRC), this cohort study enrolled 71 patients with postoperative stage III CRC. Based on preoperative CT scans, radiomic features were extracted and selected to build pixel image data using covariate-adjusted tensor classification in the high-dimension (CATCH) model. The differentially expressed RNA genes, as radiomic covariates, were identified by cancer recurrence. Predictive models were built using the pixel image and immune genomic expression factors, and the area under the curve (AUC) and F1 score were used to evaluate their performance. Significantly adjusted radiomic features were selected to predict recurrence. The association between the significantly adjusted radiomic features and immune gene expression was also investigated. Overall, 1037 radiomic features were converted into 33 × 32-pixel image data. Thirty differentially expressed genes were identified. We performed 100 iterations of 3-fold cross-validation to evaluate the performance of the CATCH model, which showed a high sensitivity of 0.66 and an F1 score of 0.69. The area under the curve (AUC) was 0.56. Overall, ten adjusted radiomic features were significantly associated with cancer recurrence in the CATCH model. All of these methods are texture-associated radiomics. Compared with non-adjusted radiomics, 7 out of 10 adjusted radiomic features influenced recurrence-free survival. The adjusted radiomic features were positively associated with PECAM1 , PRDM1 , AIF1 , IL10 , ISG20 , and TLR8 expression. We provide individualized cancer therapeutic strategies based on adjusted radiomic features in recurrent stage III CRC. Adjusted CT scan image-based radiomics with immune genomic expression covariates using the CATCH model can efficiently predict cancer recurrence. The correlation between adjusted radiomic features and immune genomic expression can provide biological relevance and individualized therapeutic targets.

Background: Neonatal hypoxic-ischemic encephalopathy (HIE) is the most common cause of mortality and neurological disability in infancy after perinatal asphyxia. Reliable biomarkers to predict neurological outcomes of neonates after perinatal asphyxia are still not accessible in clinical practice., Methods: A prospective cohort study enrolled neonates with perinatal asphyxia. Biochemical blood tests and cerebral Doppler ultrasound were measured within 6 h of age and at the 4th day old. Neurological outcomes were assessed at 1 year old., Results: Sixty-four neonates with perinatal asphyxia were enrolled. Fifty-eight (90%) had hypoxic-ischemic encephalopathy (HIE) including 20 (34%) Stage I, 21 (36%) Stage II, and 17 (29%) Stage III. In the asphyxiated infants without therapeutic hypothermia, HIE stage, PH, and base excess levels within 6 h of age were the predictors of adverse outcomes. In the asphyxiated infants receiving therapeutic hypothermia, HIE stage failed to predict outcomes. Instead, blood lactate levels and pulsatility index (PI) of medial cerebral arteries (MCA) either in 6 h of age or at the 4th day old independently predicted adverse outcomes., Conclusions: Blood lactate, which is a common accessible test at the hospital and MCA PI on cerebral ultrasound could predict adverse outcomes in asphyxiated infants receiving therapeutic hypothermia.

Background: The Orchid family is the largest families of the monocotyledons and an economically important ornamental plant worldwide. Given the pivotal role of this plant to humans, botanical researchers and breeding communities should have access to valuable genomic and transcriptomic information of this plant. Previously, we established OrchidBase, which contains expressed sequence tags (ESTs) from different tissues and developmental stages of Phalaenopsis as well as biotic and abiotic stress-treated Phalaenopsis. The database includes floral transcriptomic sequences from 10 orchid species across all the five subfamilies of Orchidaceae., Description: Recently, the whole-genome sequences of Apostasia shenzhenica, Dendrobium catenatum, and Phalaenopsis equestris were de novo assembled and analyzed. These datasets were used to develop OrchidBase 4.0, including genomic and transcriptomic data for these three orchid species. OrchidBase 4.0 offers information for gene annotation, gene expression with fragments per kilobase of transcript per millions mapped reads (FPKM), KEGG pathways and BLAST search. In addition, assembled genome sequences and location of genes and miRNAs could be visualized by the genome browser. The online resources in OrchidBase 4.0 can be accessed by browsing or using BLAST. Users can also download the assembled scaffold sequences and the predicted gene and protein sequences of these three orchid species., Conclusions: OrchidBase 4.0 is the first database that contain the whole-genome sequences and annotations of multiple orchid species. OrchidBase 4.0 is available at http://orchidbase.itps.ncku.edu.tw/., (© 2021. The Author(s).)

Orchid gynostemium, the fused organ of the androecium and gynoecium, and ovule development are unique developmental processes. Two DROOPING LEAF/CRABS CLAW (DL/CRC) genes, PeDL1 and PeDL2, were identified from the Phalaenopsis orchid genome and functionally characterized. Phylogenetic analysis indicated that the most recent common ancestor of orchids contained the duplicated DL/CRC-like genes. Temporal and spatial expression analysis indicated that PeDL genes are specifically expressed in the gynostemium and at the early stages of ovule development. Both PeDLs could partially complement an Arabidopsis crc-1 mutant. Virus-induced gene silencing (VIGS) of PeDL1 and PeDL2 affected the number of protuberant ovule initials differentiated from the placenta. Transient overexpression of PeDL1 in Phalaenopsis orchids caused abnormal development of ovule and stigmatic cavity of gynostemium. PeDL1, but not PeDL2, could form a heterodimer with Phalaenopsis equestris CINCINNATA 8 (PeCIN8). Paralogous retention and subsequent divergence of the gene sequences of PeDL1 and PeDL2 in P. equestris might result in the differentiation of function and protein behaviors. These results reveal that the ancestral duplicated DL/CRC-like genes play important roles in orchid reproductive organ innovation., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Experimental Biology. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

3-Monochloropropane-1,2-diol (3-MCPD), 2,3-epoxy-1-propanol (glycidol), and their esters are well-known food contaminants mainly formed by the heat processing of certain refined oils and coexist in various kinds of foodstuffs. However, the combined health effect and the underlying mechanism of 3-MCPD and glycidol coexposure are not well-understood. In this study, we investigated the systemic toxicity effects and the nephrotoxicity mechanisms of 3-MCPD and glycidol coexposure with in vitro and in vivo models, and next-generation sequencing (NGS) analysis. It was found that 3-MCPD and glycidol coexposure for 28 days synergistically induced toxicity in the kidney, lung, testis, and heart in C57BL/6 mice. Kidney was the most sensitive organ to coexposure, and the coexposure had a synergistic effect on inflammation and cytotoxicity through activation of the NLRP3 inflammasome, and the induction of necroptosis, and autophagic cell death in NRK-52E cells. Moreover, the NGS results revealed the genes changes associated with nephrotoxicity, inflammation and with the broad toxicity effects induced by 3-MCPD or glycidol alone or in combination, which were consistent with the results of in vitro and in vivo models. In summary, we report for the first time of the comprehensive toxicity effects and the mechanisms caused by 3-MCPD and glycidol coexposure., (Copyright © 2020 The Authors. Published by Elsevier B.V. All rights reserved.)

Background and Objective . The main purpose of this study was to develop a simple automatic diagnostic classification scheme for chemotherapy-induced peripheral neuropathy., Methods: This was a prospective cohort study that enrolled patients with colorectal or gynecologic cancer post chemotherapy for more than 1 year. The patients underwent laboratory examinations (nerve conduction studies and quantitative sensory tests), and a questionnaire about the quality of life. An unsupervised classification algorithm was used to classify the patients into groups using a small number of variables derived from the laboratory tests. A panel of five neurologists also diagnosed the types of neuropathies according to the laboratory tests. The results by the unsupervised classification algorithm and the neurologists were compared., Results: The neurologists' diagnoses showed much higher rates of entrapment syndromes (66.1%) and radiculopathies (55.1%) than polyneuropathy (motor/sensory: 33.1%/29.7%). A multivariate analysis showed that the questionnaire was not significantly correlated with the results of quantitative sensory tests ( r = 0.27) or the neurologists' diagnoses ( r = 0.27) or the neurologists' diagnoses (., Conclusion: The results of our unsupervised classification algorithm based on three variables of laboratory tests correlated well with the neurologists' diagnoses., Competing Interests: None of the authors have any disclosures of conflict of interest., (Copyright © 2020 Han-Wei Huang et al.)

Power/sample size (power) analysis estimates the likelihood of successfully finding the statistical significance in a data set. There has been a growing recognition of the importance of power analysis in the proper design of experiments. Power analysis is complex, yet necessary for the success of large studies. It is important to design a study that produces statistically accurate and reliable results. Power computation methods have been well established for both microarray-based gene expression studies and genotyping microarray-based genome-wide association studies. High-throughput sequencing (HTS) has greatly enhanced our ability to conduct biomedical studies at the highest possible resolution (per nucleotide). However, the complexity of power computations is much greater for sequencing data than for the simpler genotyping array data. Research on methods of power computations for HTS-based studies has been recently conducted but is not yet well known or widely used. In this article, we describe the power computation methods that are currently available for a range of HTS-based studies, including DNA sequencing, RNA-sequencing, microbiome sequencing and chromatin immunoprecipitation sequencing. Most importantly, we review the methods of power analysis for several types of sequencing data and guide the reader to the relevant methods for each data type.

Background: One of the most important and often neglected components of a successful RNA sequencing (RNA-Seq) experiment is sample size estimation. A few negative binomial model-based methods have been developed to estimate sample size based on the parameters of a single gene. However, thousands of genes are quantified and tested for differential expression simultaneously in RNA-Seq experiments. Thus, additional issues should be carefully addressed, including the false discovery rate for multiple statistic tests, widely distributed read counts and dispersions for different genes., Results: To solve these issues, we developed a sample size and power estimation method named RnaSeqSampleSize, based on the distributions of gene average read counts and dispersions estimated from real RNA-seq data. Datasets from previous, similar experiments such as the Cancer Genome Atlas (TCGA) can be used as a point of reference. Read counts and their dispersions were estimated from the reference's distribution; using that information, we estimated and summarized the power and sample size. RnaSeqSampleSize is implemented in R language and can be installed from Bioconductor website. A user friendly web graphic interface is provided at http://cqs.mc.vanderbilt.edu/shiny/RnaSeqSampleSize/ ., Conclusions: RnaSeqSampleSize provides a convenient and powerful way for power and sample size estimation for an RNAseq experiment. It is also equipped with several unique features, including estimation for interested genes or pathway, power curve visualization, and parameter optimization.

The genome-wide identification of microRNA transcription start sites (miRNA TSSs) is essential for understanding how miRNAs are regulated in development and disease. In this study, we developed mirSTP (mirna transcription Start sites Tracking Program), a probabilistic model for identifying active miRNA TSSs from nascent transcriptomes generated by global run-on sequencing (GRO-seq) and precision run-on sequencing (PRO-seq). MirSTP takes advantage of characteristic bidirectional transcription signatures at active TSSs in GRO/PRO-seq data, and provides accurate TSS prediction for human intergenic miRNAs at a high resolution. MirSTP performed better than existing generalized and experiment specific methods, in terms of the enrichment of various promoter-associated marks. MirSTP analysis of 27 human cell lines in 183 GRO-seq and 28 PRO-seq experiments identified TSSs for 480 intergenic miRNAs, indicating a wide usage of alternative TSSs. By integrating predicted miRNA TSSs with matched ENCODE transcription factor (TF) ChIP-seq data, we connected miRNAs into the transcriptional circuitry, which provides a valuable source for understanding the complex interplay between TF and miRNA. With mirSTP, we not only predicted TSSs for 72 miRNAs, but also identified 12 primary miRNAs with significant RNA polymerase pausing alterations after JQ1 treatment; each miRNA was further validated through BRD4 binding to its predicted promoter. MirSTP is available at http://bioinfo.vanderbilt.edu/mirSTP/., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

To assess the effect of chemotherapy on mitochondrial genome mutations in cancer survivors and their offspring, a study sequenced the full mitochondrial genome and determined the mitochondrial DNA heteroplasmic (mtDNA) mutation rate. To build a model for counts of heteroplasmic mutations in mothers and their offspring, bivariate Poisson regression was used to examine the relationship between mutation count and clinical information while accounting for the paired correlation. However, if the sequencing depth is not adequate, a limited fraction of the mtDNA will be available for variant calling. The classical bivariate Poisson regression model treats the offset term as equal within pairs; thus, it cannot be applied directly. In this research, we propose an extended bivariate Poisson regression model that has a more general offset term to adjust the length of the accessible genome for each observation. We evaluate the performance of the proposed method with comprehensive simulations, and the results show that the regression model provides unbiased parameter estimations. The use of the model is also demonstrated using the paired mtDNA dataset.

As RNA-seq rapidly develops and costs continually decrease, the quantity and frequency of samples being sequenced will grow exponentially. With proteomic investigations becoming more multivariate and quantitative, determining a study's optimal sample size is now a vital step in experimental design. Current methods for calculating a study's required sample size are mostly based on the hypothesis testing framework, which assumes each gene count can be modeled through Poisson or negative binomial distributions; however, these methods are limited when it comes to accommodating covariates. To address this limitation, we propose an estimating procedure based on the generalized linear model. This easy-to-use method constructs a representative exemplary dataset and estimates the conditional power, all without requiring complicated mathematical approximations or formulas. Even more attractive, the downstream analysis can be performed with current R/Bioconductor packages. To demonstrate the practicability and efficiency of this method, we apply it to three real-world studies, and introduce our on-line calculator developed to determine the optimal sample size for a RNA-seq study.

Triple-negative breast cancer (TNBC) and other molecularly heterogeneous malignancies present a significant clinical challenge due to a lack of high-frequency "driver" alterations amenable to therapeutic intervention. These cancers often exhibit genomic instability, resulting in chromosomal rearrangements that affect the structure and expression of protein-coding genes. However, identification of these rearrangements remains technically challenging. Using a newly developed approach that quantitatively predicts gene rearrangements in tumor-derived genetic material, we identified and characterized a novel oncogenic fusion involving the MER proto-oncogene tyrosine kinase (MERTK) and discovered a clinical occurrence and cell line model of the targetable FGFR3-TACC3 fusion in TNBC. Expanding our analysis to other malignancies, we identified a diverse array of novel and known hybrid transcripts, including rearrangements between noncoding regions and clinically relevant genes such as ALK, CSF1R, and CD274/PD-L1 The over 1,000 genetic alterations we identified highlight the importance of considering noncoding gene rearrangement partners, and the targetable gene fusions identified in TNBC demonstrate the need to advance gene fusion detection for molecularly heterogeneous cancers. Cancer Res; 76(16); 4850-60. ©2016 AACR., (©2016 American Association for Cancer Research.)

Traditionally, somatic mutations are detected by examining DNA sequence. The maturity of sequencing technology has allowed researchers to screen for somatic mutations in the whole genome. Increasingly, researchers have become interested in identifying somatic mutations through RNAseq data. With this motivation, we evaluated the practicability of detecting somatic mutations from RNAseq data. Current somatic mutation calling tools were designed for DNA sequencing data. To increase performance on RNAseq data, we developed a somatic mutation caller GLMVC based on bias reduced generalized linear model for both DNA and RNA sequencing data. Through comparison with MuTect and Varscan we showed that GLMVC performed better for somatic mutation detection using exome sequencing or RNAseq data. GLMVC is freely available for download at the following website: https://github.com/shengqh/GLMVC/wiki., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Background: Although advances in sequencing technologies have popularized the use of microRNA (miRNA) sequencing (miRNA-seq) for the quantification of miRNA expression, questions remain concerning the optimal methodologies for analysis and utilization of the data. The construction of a miRNA sequencing library selects RNA by length rather than type. However, as we have previously described, miRNAs represent only a subset of the species obtained by size selection. Consequently, the libraries obtained for miRNA sequencing also contain a variety of additional species of small RNAs. This study looks at the prevalence of these other species obtained from bone marrow aspirate specimens and explores the predictive value of these small RNAs in the determination of response to therapy in myelodysplastic syndromes (MDS)., Methods: Paired pre and post treatment bone marrow aspirate specimens were obtained from patients with MDS who were treated with either azacytidine or decitabine (24 pre-treatment specimens, 23 post-treatment specimens) with 22 additional non-MDS control specimens. Total RNA was extracted from these specimens and submitted for next generation sequencing after an additional size exclusion step to enrich for small RNAs. The species of small RNAs were enumerated, single nucleotide variants (SNVs) identified, and finally the differential expression of tRNA-derived species (tDRs) in the specimens correlated with diseasestatus and response to therapy., Results: Using miRNA sequencing data generated from bone marrow aspirate samples of patients with known MDS (N = 47) and controls (N = 23), we demonstrated that transfer RNA (tRNA) fragments (specifically tRNA halves, tRHs) are one of the most common species of small RNA isolated from size selection. Using tRNA expression values extracted from miRNA sequencing data, we identified six tRNA fragments that are differentially expressed between MDS and normal samples. Using the elastic net method, we identified four tRNAs-derived small RNAs (tDRs) that together can explain 67 % of the variation in treatment response for MDS patients. Similar analysis of specifically mitochondrial tDRs (mt-tDRs) identified 13 mt-tDRs which distinguished disease status in the samples and a single mt-tDR which predited response. Finally, 14 SNVs within the tDRs were found in at least 20 % of the MDS samples and were not observed in any of the control specimens., Discussion: This study highlights the prevalence of tDRs in RNA-seq studies focused on small RNAs. The potential etiologies of these species, both technical and biologic, are discussed as well as important challenges in the interpretation of tDR data., Conclusions: Our analysis results suggest that tRNA fragments can be accurately detected through miRNA sequencing data and that the expression of these species may be useful in the diagnosis of MDS and the prediction of response to therapy.

Recent studies have emphasized a role of adaptive immunity, and particularly T cells, in the genesis of hypertension. We sought to determine the T-cell subtypes that contribute to hypertension and renal inflammation in angiotensin II-induced hypertension. Using T-cell receptor spectratyping to examine T-cell receptor usage, we demonstrated that CD8(+) cells, but not CD4(+) cells, in the kidney exhibited altered T-cell receptor transcript lengths in Vβ3, 8.1, and 17 families in response to angiotensin II-induced hypertension. Clonality was not observed in other organs. The hypertension caused by angiotensin II in CD4(-/-) and MHCII(-/-) mice was similar to that observed in wild-type mice, whereas CD8(-/-) mice and OT1xRAG-1(-/-) mice, which have only 1 T-cell receptor, exhibited a blunted hypertensive response to angiotensin II. Adoptive transfer of pan T cells and CD8(+) T cells but not CD4(+)/CD25(-) cells conferred hypertension to RAG-1(-/-) mice. In contrast, transfer of CD4(+)/CD25(+) cells to wild-type mice receiving angiotensin II decreased blood pressure. Mice treated with angiotensin II exhibited increased numbers of kidney CD4(+) and CD8(+) T cells. In response to a sodium/volume challenge, wild-type and CD4(-/-) mice infused with angiotensin II retained water and sodium, whereas CD8(-/-) mice did not. CD8(-/-) mice were also protected against angiotensin-induced endothelial dysfunction and vascular remodeling in the kidney. These data suggest that in the development of hypertension, an oligoclonal population of CD8(+) cells accumulates in the kidney and likely contributes to hypertension by contributing to sodium and volume retention and vascular rarefaction., (© 2014 American Heart Association, Inc.)

Sample size and power determination is the first step in the experimental design of a successful study. Sample size and power calculation is required for applications for National Institutes of Health (NIH) funding. Sample size and power calculation is well established for traditional biological studies such as mouse model, genome wide association study (GWAS), and microarray studies. Recent developments in high-throughput sequencing technology have allowed RNAseq to replace microarray as the technology of choice for high-throughput gene expression profiling. However, the sample size and power analysis of RNAseq technology is an underdeveloped area. Here, we present RNAseqPS, an advanced online RNAseq power and sample size calculation tool based on the Poisson and negative binomial distributions. RNAseqPS was built using the Shiny package in R. It provides an interactive graphical user interface that allows the users to easily conduct sample size and power analysis for RNAseq experimental design. RNAseqPS can be accessed directly at http://cqs.mc.vanderbilt.edu/shiny/RNAseqPS/.

RNAseq technology is replacing microarray technology as the tool of choice for gene expression profiling. While providing much richer data than microarray, analysis of RNAseq data has been much more challenging. Among the many difficulties of RNAseq analysis, correctly adjusting for batch effect is a pivotal one for large-scale RNAseq based studies. The batch effect of RNAseq data is most obvious in microRNA (miRNA) sequencing studies. Using real miRNA sequencing (miRNAseq) data, we evaluated several batch removal techniques and discussed their effectiveness. We illustrate that by adjusting for batch effect, more reliable differentially expressed genes can be identified. Our study on batch effect in miRNAseq data can serve as a guideline for future miRNAseq studies that might contain batch effect.

Sample size determination is essential to planning clinical trials. Jung (2008) established a sample size calculation formula for paired right-censored data based on the logrank test, which has been well-studied for comparing independent survival outcomes. An alternative to rank-based methods for independent right-censored data, advocated by Pepe and Fleming (1989), tests for differences between integrated weighted Kaplan-Meier estimates and is more sensitive to the magnitude of difference in survival times between groups. In this paper, we employ the concept of the Pepe-Fleming method to determine an adequate sample size by calculating differences between Kaplan-Meier estimators considering pair-wise correlation. We specify a positive stable frailty model for the joint distribution of paired survival times. We evaluate the performance of the proposed method by simulation studies and investigate the impacts of the accrual times, follow-up times, loss to follow-up rate, and sensitivity of power under misspecification of the model. The results show that ignoring the pair-wise correlation results in overestimating the required sample size. Furthermore, the proposed method is applied to two real-world studies, and the R code for sample size calculation is made available to users.

This trial was designed to assess efficacy and safety of erlotinib with sorafenib in the treatment of patients with advanced pancreatic adenocarcinoma. An exploratory correlative study analyzing pretreatment serum samples using a multivariate protein mass spectrometry-based test (VeriStrat®), previously shown to correlate with outcomes in lung cancer patients treated with erlotinib, was performed. Patients received sorafenib 400 mg daily along with erlotinib 150 mg daily with a primary endpoint of 8-week progression free survival (PFS) rate. Pretreatment serum sample analysis by VeriStrat was done blinded to clinical and outcome data; the endpoints were PFS and overall survival (OS). Difference between groups (by VeriStrat classification) was assessed using log-rank P values; hazard ratios (HR) were obtained from Cox proportional hazards model. Thirty-six patients received study drug and were included in the survival analysis. Eight-week PFS rate of 46% (95% confidence interval (CI): 0.32-0.67) did not meet the primary endpoint of a rate ≥70%. Thirty-two patients were included in the correlative analysis, and VeriStrat "Good" patients had superior PFS (HR = 0.18, 95% CI: 0.06-0.57; P = 0.001) and OS (HR = 0.31 95% CI: 0.13-0.77, P = 0.008) compared to VeriStrat "Poor" patients. Grade 3 toxicities of this regimen included fever, anemia, diarrhea, dehydration, rash, and altered liver function. This study did not meet the primary endpoint, and this combination will not be further pursued. In this small retrospective analysis, the proteomic classification was significantly associated with clinical outcomes and is being further evaluated in ongoing studies., (© 2014 The Authors. Cancer Medicine published by John Wiley & Sons Ltd.)

Introduction: To determine the time to progression (TTP), response rate (RR), and toxicity for North American patients with relapsed small-cell lung cancer (SCLC) treated with bendamustine in the second- or third-line setting., Methods: Patients with relapsed, histologically confirmed SCLC were eligible for enrollment on study. The study population included patients with both chemotherapy-sensitive and chemotherapy-resistant disease treated with up to two prior lines of chemotherapy. Patients were treated with 120 mg/m of bendamustine on days 1 and 2 of a 21-day cycle for up to six cycles. Primary end point was TTP; secondary end points included toxicity, RR, and overall survival., Results: Fifty-nine patients were accrued, 50 patients met eligibility for enrollment. The median age of patients was 62, and 56% were men. Twenty-nine patients (58%) had chemotherapy-sensitive disease. Median TTP was 4.0 months (95% confidence interval [CI], 3.3-5.4), median overall survival was 4.8 months (95% CI, 3.8-6.3), and the RR was 26% (95% CI, 13.3%-39.5%). Patients with chemosensitive disease had a median TTP of 4.2 months (95% CI, 3.3-6.0) compared with 3.4 months (95% CI, 2.7-∞) for chemotherapy-resistant disease. The RR was 33% (95% CI, 14.2%-51.8%) in patients with chemosensitive disease and 17% (95% CI, 0%-34.4%) in those with chemoresistant disease. The most common grade 3/4 adverse events were fatigue (20%), dyspnea (12%), and anemia (12%)., Conclusion: Bendamustine has modest activity in relapsed SCLC similar to other agents evaluated in this patient population.

Background and Purpose: To determine if dose and/or dose-volume parameters to anatomic swallowing structures are predictive of gastrostomy tube (PEG) dependence from chemotherapy-intensity modulated radiotherapy (IMRT) in locally advanced head and neck cancer (LAHNC)., Methods and Materials: A retrospective study was performed on 141 consecutive patients with LAHNC (squamous cell) treated with definitive chemoIMRT with weekly concurrent carboplatin and paclitaxel. Late dysphagia was assessed by length of PEG requirement. Analysis of IMRT dose was retrospectively performed for critical swallowing structures., Results: Approximately 62% of patients required PEG, the majority placed during treatment. Mean and median time for PEG was 7.7 and 4.4 months respectively (range 1.4-43.8). Only IMRT dose to the inferior constrictor was significantly associated with length of PEG. Mean dose (of individual mean doses) was 47 Gy for prolonged PEG use versus 41 Gy for PEG ⩽ 12 months. V40 to the inferior constrictor also correlated with PEG >12 months (p = 0.02) with a mean V40 of 48% versus 41% for PEG ⩽ 12 months., Conclusions: IMRT dose to the inferior constrictor correlated with persistent dysphagia requiring prolonged PEG use. Maintaining mean inferior constrictor dose to ⩽ 41 Gy and V40 to ⩽ 41% may help minimize gastrostomy tube dependence., (Copyright © 2014. Published by Elsevier Ireland Ltd.)

Little is known about the inheritance of very low heteroplasmy mitochondria DNA (mtDNA) variations. Even with the development of new next-generation sequencing methods, the practical lower limit of measured heteroplasmy is still about 1% due to the inherent noise level of the sequencing. In this study, we sequenced the mitochondrial genome of 44 individuals using Illumina high-throughput sequencing technology and obtained high-coverage mitochondria sequencing data. Our study population contains many mother-offspring pairs. This unique study design allows us to bypass the usual heteroplasmy limitation by analyzing the correlation of mutation levels at each position in the mtDNA sequence between maternally related pairs and non-related pairs. The study showed that very low heteroplasmy variants, down to almost 0.1%, are inherited maternally and that this inheritance begins to decrease at about 0.5%, corresponding to a bottleneck of about 200 mtDNA., (Copyright © 2013. Published by Elsevier Ltd.)

Background: Sample size calculation is an important issue in the experimental design of biomedical research. For RNA-seq experiments, the sample size calculation method based on the Poisson model has been proposed; however, when there are biological replicates, RNA-seq data could exhibit variation significantly greater than the mean (i.e. over-dispersion). The Poisson model cannot appropriately model the over-dispersion, and in such cases, the negative binomial model has been used as a natural extension of the Poisson model. Because the field currently lacks a sample size calculation method based on the negative binomial model for assessing differential expression analysis of RNA-seq data, we propose a method to calculate the sample size., Results: We propose a sample size calculation method based on the exact test for assessing differential expression analysis of RNA-seq data., Conclusions: The proposed sample size calculation method is straightforward and not computationally intensive. Simulation studies to evaluate the performance of the proposed sample size method are presented; the results indicate our method works well, with achievement of desired power.

Motivation: Exome capture kits have capture efficiencies that range from 40 to 60%. A significant amount of off-target reads are from the mitochondrial genome. These unintentionally sequenced mitochondrial reads provide unique opportunities to study the mitochondria genome., Results: MitoSeek is an open-source software tool that can reliably and easily extract mitochondrial genome information from exome and whole genome sequencing data. MitoSeek evaluates mitochondrial genome alignment quality, estimates relative mitochondrial copy numbers and detects heteroplasmy, somatic mutation and structural variants of the mitochondrial genome. MitoSeek can be set up to run in parallel or serial on large exome sequencing datasets., Availability: https://github.com/riverlee/MitoSeek

Background: RNAseq technology is replacing microarray technology as the tool of choice for gene expression profiling. While providing much richer data than microarray, analysis of RNAseq data has been much more challenging. To date, there has not been a consensus on the best approach for conducting robust RNAseq analysis., Results: In this study, we designed a thorough experiment to evaluate six read count-based RNAseq analysis methods (DESeq, DEGseq, edgeR, NBPSeq, TSPM and baySeq) using both real and simulated data. We found the six methods produce similar fold changes and reasonable overlapping of differentially expressed genes based on p-values. However, all six methods suffer from over-sensitivity., Conclusions: Based on the evaluation of runtime using real data and area under the receiver operating characteristic curve (AUC-ROC) using simulated data, we found that edgeR achieves a better balance between speed and accuracy than the other methods.

Next-generation sequencing (NGS) technology has provided researchers with opportunities to study the genome in unprecedented detail. In particular, NGS is applied to disease association studies. Unlike genotyping chips, NGS is not limited to a fixed set of SNPs. Prices for NGS are now comparable to the SNP chip, although for large studies the cost can be substantial. Pooling techniques are often used to reduce the overall cost of large-scale studies. In this study, we designed a rigorous simulation model to test the practicability of estimating allele frequency from pooled sequencing data. We took crucial factors into consideration, including pool size, overall depth, average depth per sample, pooling variation, and sampling variation. We used real data to demonstrate and measure reference allele preference in DNAseq data and implemented this bias in our simulation model. We found that pooled sequencing data can introduce high levels of relative error rate (defined as error rate divided by targeted allele frequency) and that the error rate is more severe for low minor allele frequency SNPs than for high minor allele frequency SNPs. In order to overcome the error introduced by pooling, we recommend a large pool size and high average depth per sample.

Sample size determination is an important issue in the experimental design of biomedical research. Because of the complexity of RNA-seq experiments, however, the field currently lacks a sample size method widely applicable to differential expression studies utilising RNA-seq technology. In this report, we propose several methods for sample size calculation for single-gene differential expression analysis of RNA-seq data under Poisson distribution. These methods are then extended to multiple genes, with consideration for addressing the multiple testing problem by controlling false discovery rate. Moreover, most of the proposed methods allow for closed-form sample size formulas with specification of the desired minimum fold change and minimum average read count, and thus are not computationally intensive. Simulation studies to evaluate the performance of the proposed sample size formulas are presented; the results indicate that our methods work well, with achievement of desired power. Finally, our sample size calculation methods are applied to three real RNA-seq data sets.

Background: When using Illumina high throughput short read data, sometimes the genotype inferred from the positive strand and negative strand are significantly different, with one homozygous and the other heterozygous. This phenomenon is known as strand bias. In this study, we used Illumina short-read sequencing data to evaluate the effect of strand bias on genotyping quality, and to explore the possible causes of strand bias., Result: We collected 22 breast cancer samples from 22 patients and sequenced their exome using the Illumina GAIIx machine. By comparing the consistency between the genotypes inferred from this sequencing data with the genotypes inferred from SNP chip data, we found that, when using sequencing data, SNPs with extreme strand bias did not have significantly lower consistency rates compared to SNPs with low or no strand bias. However, this result may be limited by the small subset of SNPs present in both the exome sequencing and the SNP chip data. We further compared the transition and transversion ratio and the number of novel non-synonymous SNPs between the SNPs with low or no strand bias and those with extreme strand bias, and found that SNPs with low or no strand bias have better overall quality. We also discovered that the strand bias occurs randomly at genomic positions across these samples, and observed no consistent pattern of strand bias location across samples. By comparing results from two different aligners, BWA and Bowtie, we found very consistent strand bias patterns. Thus strand bias is unlikely to be caused by alignment artifacts. We successfully replicated our results using two additional independent datasets with different capturing methods and Illumina sequencers., Conclusion: Extreme strand bias indicates a potential high false-positive rate for SNPs.

Background: Exome sequencing using next-generation sequencing technologies is a cost efficient approach to selectively sequencing coding regions of human genome for detection of disease variants. A significant amount of DNA fragments from the capture process fall outside target regions, and sequence data for positions outside target regions have been mostly ignored after alignment., Result: We performed whole exome sequencing on 22 subjects using Agilent SureSelect capture reagent and 6 subjects using Illumina TrueSeq capture reagent. We also downloaded sequencing data for 6 subjects from the 1000 Genomes Project Pilot 3 study. Using these data, we examined the quality of SNPs detected outside target regions by computing consistency rate with genotypes obtained from SNP chips or the Hapmap database, transition-transversion (Ti/Tv) ratio, and percentage of SNPs inside dbSNP. For all three platforms, we obtained high-quality SNPs outside target regions, and some far from target regions. In our Agilent SureSelect data, we obtained 84,049 high-quality SNPs outside target regions compared to 65,231 SNPs inside target regions (a 129% increase). For our Illumina TrueSeq data, we obtained 222,171 high-quality SNPs outside target regions compared to 95,818 SNPs inside target regions (a 232% increase). For the data from the 1000 Genomes Project, we obtained 7,139 high-quality SNPs outside target regions compared to 1,548 SNPs inside target regions (a 461% increase)., Conclusions: These results demonstrate that a significant amount of high quality genotypes outside target regions can be obtained from exome sequencing data. These data should not be ignored in genetic epidemiology studies.

The human mitochondrial genome has an exclusively maternal mode of inheritance. Mitochondrial DNA (mtDNA) is particularly vulnerable to environmental insults due in part to an underdeveloped DNA repair system, limited to base excision and homologous recombination repair. Radiation exposure to the ovaries may cause mtDNA mutations in oocytes, which may in turn be transmitted to offspring. We hypothesized that the children of female cancer survivors who received radiation therapy may have an increased rate of mtDNA heteroplasmy mutations, which conceivably could increase their risk of developing cancer and other diseases. We evaluated 44 DNA blood samples from 17 Danish and 1 Finnish families (18 mothers and 26 children). All mothers had been treated for cancer as children and radiation doses to their ovaries were determined based on medical records and computational models. DNA samples were sequenced for the entire mitochondrial genome using the Illumina GAII system. Mother's age at sample collection was positively correlated with mtDNA heteroplasmy mutations. There was evidence of heteroplasmy inheritance in that 9 of the 18 families had at least one child who inherited at least one heteroplasmy site from his or her mother. No significant difference in single nucleotide polymorphisms between mother and offspring, however, was observed. Radiation therapy dose to ovaries also was not significantly associated with the heteroplasmy mutation rate among mothers and children. No evidence was found that radiotherapy for pediatric cancer is associated with the mitochondrial genome mutation rate in female cancer survivors and their children., (Copyright © 2012 Elsevier B.V. All rights reserved.)

Traditional wheezes detection methods are based on the frequency and durations of acoustic signal or the location of peaks from successive spectra. In these methods, the discriminative threshold used to identify peaks usually is fixed empirically. Therefore, accuracy of detected wheeze is affected by environment noise and artificial factors. The objective of this study is to classify normal and abnormal (wheezing) respiratory sounds using Cepstral analysis in Gaussian Mixture Models. The sound signal is divided in overlapped segments, which are characterized by a reduced dimension feature vectors using Mel-Frequency Cepstral Coefficients. In this study the ;speaker' is wheeze. During the test phase, an unknown sound is compared to all the GMM models and the classification decision is based on the Maximum Likelihood criterion. In these processes, identification is based on threshold value. If the threshold is bigger than zero, the sound is normal. Otherwise, the sound is wheeze. From experimental results, when the Gaussian mix number is 16, the accuracy of identification of wheeze is up to 90%.