In-vivo visualization of tumor microvessel density and response to anti-angiogenic treatment by high resolution MRI in mice

Purpose

Inhibition of angiogenesis has shown clinical success in patients with cancer. Thus, imaging approaches that allow for the identification of angiogenic tumors and the detection of response to anti-angiogenic treatment are of high clinical relevance.

Experimental Design

We established an in vivo magnetic resonance imaging (MRI) approach that allows us to simultaneously image tumor microvessel density and tumor vessel size in a NSCLC model in mice.

Results

Using microvessel density imaging we demonstrated an increase in microvessel density within 8 days after tumor implantation, while tumor vessel size decreased indicating a switch from macro- to microvessels during tumor growth. Moreover, we could monitor in vivo inhibition of angiogenesis induced by the angiogenesis inhibitor PTK787, resulting in a decrease of microvessel density and a slight increase in tumor vessel size.

Conclusions

We present an in vivo imaging approach that allows us to monitor both tumor microvessel density and tumor vessel size in the tumor. Moreover, this approach enables us to assess, early-on, treatment effects on tumor microvessel density as well as on tumor vessel size. Thus, this imaging-based strategy of validating anti-angiogenic treatment effects has high potential in applications to preclinical and clinical trials.     

Informatics-guided procurement of patient samples for biomarker discovery projects in cancer research

Modern cancer research for biomarker discovery program requires solving several tasks that are directly involved with patient sample procurement. One requirement is to construct a highly efficient workflow on the clinical side for the procurement to generate a consistent supply of high quality samples for research. This undertaking needs a network of interdepartmental collaborations and participations at various levels, including physical human interactions, information technology implementations and a bioinformatics tool that is highly effective and user-friendly to busy clinicians and researchers associated with the sample procurement. Collegial participation that is sequential but continual from one department to another demands dedicated bioinformatics software coordinating between the institutional clinic and the tissue repository facility. Participants in the process include admissions, consenting process, phlebotomy, surgery center and pathology. During this multiple step procedures, clinical data are collected for detailed analytical endpoints to supplement logistics of defining and validating the discovery of biomarkers.

Serum microRNAs are promising novel biomarkers

Background

Circulating nucleic acids (CNAs) offer unique opportunities for early diagnosis of clinical conditions. Here we show that microRNAs, a family of small non-coding regulatory RNAs involved in human development and pathology, are present in bodily fluids and represent new effective biomarkers.

Methods and Results

After developing protocols for extracting and quantifying microRNAs in serum and other body fluids, the serum microRNA profiles of several healthy individuals were determined and found to be similar, validating the robustness of our methods. To address the possibility that the abundance of specific microRNAs might change during physiological or pathological conditions, serum microRNA levels in pregnant and non pregnant women were compared. In sera from pregnant women, microRNAs associated with human placenta were significantly elevated and their levels correlated with pregnancy stage.

Conclusions and Significance

Considering the central role of microRNAs in development and disease, our results highlight the medically relevant potential of determining microRNA levels in serum and other body fluids. Thus, microRNAs are a new class of CNAs that promise to serve as useful clinical biomarkers.     

Log::ProgramInfo: A Perl module to collect and log data for bioinformatics pipelines

Background

To reproduce and report a bioinformatics analysis, it is important to be able to determine the environment in which a program was run. It can also be valuable when trying to debug why different executions are giving unexpectedly different results.

Results

Log::ProgramInfo is a Perl module that writes a log file at the termination of execution of the enclosing program, to document useful execution characteristics. This log file can be used to re-create the environment in order to reproduce an earlier execution. It can also be used to compare the environments of two executions to determine whether there were any differences that might affect (or explain) their operation.

Availability

The source is available on CPAN (Macdonald and Boutros, Log-ProgramInfo. http://search.cpan.org/~boutroslb/Log-ProgramInfo/).

Conclusion

Using Log::ProgramInfo in programs creating result data for publishable research, and including the Log::ProgramInfo output log as part of the publication of that research is a valuable method to assist others to duplicate the programming environment as a precursor to validating and/or extending that research.

     

Field-based artificial crossings indicate partial compatibility of reciprocal crosses between <Emphasis Type="Italic">Pinus sylvestris</Emphasis> and <Emphasis Type="Italic">Pinus mugo</Emphasis> and unexpected chloroplast DNA inheritance

Crossability relationships between Scots pine (Pinus sylvestris L.) and mountain dwarf pine (Pinus mugo Turra) was studied, using artificial pollination approach. Partial compatibility of the reciprocal crossings of these species was proved experimentally, validating the idea of a spontaneous formation of their hybrid swarms under natural conditions. The hybrids were validated using organellar DNA markers and nuclear DNA microsatellites. Based on the percentage of filled seeds, the interspecific crossings were less efficient than the intraspecific cross-pollinations of P. sylvestris and P. mugo individuals. Both species were found to intercross readily with individuals of their putative hybrid swarm, P. mugo exhibiting a higher hybridological affinity towards putatively hybrid individuals than P. sylvestris. Validation of the hybrids confirmed the paternal inheritance of chloroplast DNA (cpDNA) in the combination P. sylvestris × P. mugo only. Surprisingly, in the reciprocal crossing P. mugo × P. sylvestris, maternal inheritance of cpDNA was revealed. Obtained results offer a new insight into the direction and intensity of gene flow within the hybrid swarms of Scots pine and mountain dwarf pine.     

Investigation of major genetic alterations in neuroblastoma

Neuroblastoma (NB) is the most common extracranial solid tumor in childhood. This malignancy shows a wide spectrum of clinical outcome and its prognosis is conditioned by manifold biological and genetic factors. We investigated the tumor genetic profile and clinical data of 29 patients with NB by multiplex ligation-dependent probe amplification (MLPA) to assess therapeutic risk. In 18 of these tumors, MYCN status was assessed by fluorescence in situ hybridization (FISH). Copy number variation was also determined for confirming MLPA findings in two 6p loci. We found 2p, 7q and 17q gains, and 1p and 11q losses as the most frequent chromosome alterations in this cohort. FISH confirmed all cases of MYCN amplification detected by MLPA. In view of unexpected 6p imbalance, copy number variation of two 6p loci was assessed for validating MLPA findings. Based on clinical data and genetic profiles, patients were stratified in pretreatment risk groups according to international consensus. MLPA proved to be effective for detecting multiple genetic alterations in all chromosome regions as requested by the International Neuroblastoma Risk Group (INRG) for therapeutic stratification. Moreover, this technique proved to be cost effective, reliable, only requiring standard PCR equipment, and attractive for routine analysis. However, the observed 6p imbalances made PKHD1 and DCDC2 inadequate for control loci. This must be considered when designing commercial MLPA kits for NB. Finally, four patients showed a normal MLPA profile, suggesting that NB might have a more complex genetic pattern than the one assessed by presently available MLPA kits.     

Development,Calibration, and Validation of a U.S. White Male Population-Based Simulation Model of Esophageal Adenocarcinoma

Background

The incidence of esophageal adenocarcinoma (EAC) has risen rapidly in the U.S. and western world. The aim of the study was to begin the investigation of this rapid rise by developing, calibrating, and validating a mathematical disease simulation model of EAC using available epidemiologic data.

Methods

The model represents the natural history of EAC, including the essential biologic health states from normal mucosa to detected cancer. Progression rates between health states were estimated via calibration, which identified distinct parameter sets producing model outputs that fit epidemiologic data; specifically, the prevalence of pre-cancerous lesions and EAC cancer incidence from the published literature and Surveillance, Epidemiology, and End Results (SEER) data. As an illustrative example of a clinical and policy application, the calibrated and validated model retrospectively analyzed the potential benefit of an aspirin chemoprevention program.

Results

Model outcomes approximated calibration targets; results of the model''s fit and validation are presented. Approximately 7,000 cases of EAC could have been prevented over a 30-year period if all white males started aspirin chemoprevention at age 40 in 1965.

Conclusions

The model serves as the foundation for future analyses to determine a cost-effective screening and management strategy to prevent EAC morbidity and mortality.     

The human telomeric proteome during telomere replication

Telomere shortening can cause detrimental diseases and contribute to aging. It occurs due to the end replication problem in cells lacking telomerase. Furthermore, recent studies revealed that telomere shortening can be attributed to difficulties of the semi-conservative DNA replication machinery to replicate the bulk of telomeric DNA repeats. To investigate telomere replication in a comprehensive manner, we develop QTIP-iPOND - Quantitative Telomeric chromatin Isolation Protocol followed by isolation of Proteins On Nascent DNA - which enables purification of proteins that associate with telomeres specifically during replication. In addition to the core replisome, we identify a large number of proteins that specifically associate with telomere replication forks. Depletion of several of these proteins induces telomere fragility validating their importance for telomere replication. We also find that at telomere replication forks the single strand telomere binding protein POT1 is depleted, whereas histone H1 is enriched. Our work reveals the dynamic changes of the telomeric proteome during replication, providing a valuable resource of telomere replication proteins. To our knowledge, this is the first study that examines the replisome at a specific region of the genome.     

Design Characteristics Influence Performance of Clinical Prediction Rules in Validation: A Meta-Epidemiological Study

Background

Many new clinical prediction rules are derived and validated. But the design and reporting quality of clinical prediction research has been less than optimal. We aimed to assess whether design characteristics of validation studies were associated with the overestimation of clinical prediction rules’ performance. We also aimed to evaluate whether validation studies clearly reported important methodological characteristics.

Methods

Electronic databases were searched for systematic reviews of clinical prediction rule studies published between 2006 and 2010. Data were extracted from the eligible validation studies included in the systematic reviews. A meta-analytic meta-epidemiological approach was used to assess the influence of design characteristics on predictive performance. From each validation study, it was assessed whether 7 design and 7 reporting characteristics were properly described.

Results

A total of 287 validation studies of clinical prediction rule were collected from 15 systematic reviews (31 meta-analyses). Validation studies using case-control design produced a summary diagnostic odds ratio (DOR) 2.2 times (95% CI: 1.2–4.3) larger than validation studies using cohort design and unclear design. When differential verification was used, the summary DOR was overestimated by twofold (95% CI: 1.2 -3.1) compared to complete, partial and unclear verification. The summary RDOR of validation studies with inadequate sample size was 1.9 (95% CI: 1.2 -3.1) compared to studies with adequate sample size. Study site, reliability, and clinical prediction rule was adequately described in 10.1%, 9.4%, and 7.0% of validation studies respectively.

Conclusion

Validation studies with design shortcomings may overestimate the performance of clinical prediction rules. The quality of reporting among studies validating clinical prediction rules needs to be improved.     

Map-based cloning and functional analysis of the chromogen gene <Emphasis Type="Italic">C</Emphasis> in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

The chromogen gene C is critical for anthocyanin regulation in rice, and apiculus color is an important agronomic trait in selective breeding and variety purification. Mapbased cloning and in-depth functional analysis of the C gene will be useful for understanding the molecular mechanism of anthocyanin biosynthesis and for rice breeding. Japonica landrace Lijiangxintuanheigu (LTH) has red apiculi and purple stigmas. Genetic analysis showed that red apiculus and purple stigma in LTH co-segregated indicating control by a single dominant gene, or by two completely linked genes. Using 1,851 recessive individuals from two F₂ populations, the target gene OsC was delimited to a 70.8 kb interval on chromosome 6 that contains the rice homologue of the maize anthocyanin regulatory gene C1. When the entire OsC gene and its full-length cDNA cloned from LTH were transformed into japonica cultivar Kitaake with colorless apiculi and stigmas all positive transformants had red apiculi but non-colored stigmas, validating that OsC alone was responsible for the apiculus color and represented the functional C gene. OsC was constitutively expressed in all tissues examined, with strongest expression in leaf blades. These results set a foundation to clarify the regulatory mechanisms of OsC in the anthocyanin biosynthetic pathway.     

Genetic identification of SNP markers linked to a new grape phylloxera resistant locus in <Emphasis Type="Italic">Vitis cinerea</Emphasis> for marker-assisted selection

Background

Grape phylloxera (Daktulosphaira vitifoliae Fitch) is a major insect pest that negatively impacts commercial grapevine performance worldwide. Consequently, the use of phylloxera resistant rootstocks is an essential component of vineyard management. However, the majority of commercially available rootstocks used in viticulture production provide limited levels of grape phylloxera resistance, in part due to the adaptation of phylloxera biotypes to different Vitis species. Therefore, there is pressing need to develop new rootstocks better adapted to specific grape growing regions with complete resistance to grape phylloxera biotypes.

Results

Grapevine rootstock breeding material, including an accession of Vitis cinerea and V. aestivalis, DRX55 ([M. rotundifolia x V. vinifera] x open pollinated) and MS27-31 (M. rotundifolia specific hybrid), provided complete resistance to grape phylloxera in potted plant assays. To map the genetic factor(s) of grape phylloxera resistance, a F₁ V. cinerea x V. vinifera Riesling population was screened for resistance. Heritability analysis indicates that the V. cinerea accession contained a single allele referred as RESISTANCE TO DAKTULOSPHAIRA VITIFOLIAE 2 (RDV2) that confers grape phylloxera resistance. Using genetic maps constructed with pseudo-testcross markers for V. cinerea and Riesling, a single phylloxera resistance locus was identified in V. cinerea. After validating SNPs at the RDV2 locus, interval and linkage mapping showed that grape phylloxera resistance mapped to linkage group 14 at position 16.7 cM.

Conclusion

The mapping of RDV2 and the validation of markers linked to grape phylloxera resistance provides the basis to breed new rootstocks via marker-assisted selection that improve vineyard performance.

     

Growth impairment after TBI of leukemia survivors children: a model- based investigation

Background

Children receiving Total Body Irradiation (TBI) in preparation for Hematopoietic Stem Cell Transplantation (HSCT) are at risk for Growth Hormone Deficiency (GHD), which sometimes severely compromises their Final Height (FH). To better represent the impact of such therapies on growth we apply a mathematical model, which accounts both for the gompertzian-like growth trend and the hormone-related ‘spurts’, and evaluate how the parameter values estimated on the children undergoing TBI differ from those of the matched normal population.

Methods

25 patients long-term childhood lymphoblastic and myeloid acute leukaemia survivors followed at Pediatric Onco-Hematology, Stem Cell Transplantation and Cellular Therapy Division, Regina Margherita Children’s Hospital (Turin, Italy) were retrospectively analysed for assessing the influence of TBI on their longitudinal growth and for validating a new method to estimate the GH therapy effects. Six were treated with GH therapy after a GHD diagnosis.

Results

We show that when TBI was performed before puberty overall growth and pubertal duration were significantly impaired, but such growth limitations were completely reverted in the small sample (6 over 25) of children who underwent GH replacement therapies.

Conclusion

Since in principle the model could account for any additional growth ‘spurt’ induced by therapy, it may become a useful ‘simulation’ tool for paediatricians for comparing the predicted therapy effectiveness depending on its timing and dosage.

     

Generation of comprehensive transposon insertion mutant library for the model archaeon, <Emphasis Type="Italic">Haloferax volcanii</Emphasis>, and its use for gene discovery

Background

Archaea share fundamental properties with bacteria and eukaryotes. Yet, they also possess unique attributes, which largely remain poorly characterized. Haloferax volcanii is an aerobic, moderately halophilic archaeon that can be grown in defined media. It serves as an excellent archaeal model organism to study the molecular mechanisms of biological processes and cellular responses to changes in the environment. Studies on haloarchaea have been impeded by the lack of efficient genetic screens that would facilitate the identification of protein functions and respective metabolic pathways.

Results

Here, we devised an insertion mutagenesis strategy that combined Mu in vitro DNA transposition and homologous-recombination-based gene targeting in H. volcanii. We generated an insertion mutant library, in which the clones contained a single genomic insertion. From the library, we isolated pigmentation-defective and auxotrophic mutants, and the respective insertions pinpointed a number of genes previously known to be involved in carotenoid and amino acid biosynthesis pathways, thus validating the performance of the methodologies used. We also identified mutants that had a transposon insertion in a gene encoding a protein of unknown or putative function, demonstrating that novel roles for non-annotated genes could be assigned.

Conclusions

We have generated, for the first time, a random genomic insertion mutant library for a halophilic archaeon and used it for efficient gene discovery. The library will facilitate the identification of non-essential genes behind any specific biochemical pathway. It represents a significant step towards achieving a more complete understanding of the unique characteristics of halophilic archaea.

     

Viral Cre-LoxP tools aid genome engineering in mammalian cells

Background

Targeted nucleases have transformed genome editing technology, providing more efficient methods to make targeted changes in mammalian genome. In parallel, there is an increasing demand of Cre-LoxP technology for complex genome manipulation such as large deletion, addition, gene fusion and conditional removal of gene sequences at the target site. However, an efficient and easy-to-use Cre-recombinase delivery system remains lacking.

Results

We designed and constructed two sets of expression vectors for Cre-recombinase using two highly efficient viral systems, the integrative lentivirus and non-integrative adeno associated virus. We demonstrate the effectiveness of those methods in Cre-delivery into stably-engineered HEK293 cells harboring LoxP-floxed red fluorescent protein (RFP) and puromycin (Puro) resistant reporters. The delivered Cre recombinase effectively excised the floxed RFP-Puro either directly or conditionally, therefore validating the function of these molecular tools. Given the convenient options of two selections markers, these viral-based systems offer a robust and easy-to-use tool for advanced genome editing, expanding complicated genome engineering to a variety of cell types and conditions.

Conclusions

We have developed and functionally validated two viral-based Cre-recombinase delivery systems for efficient genome manipulation in various mammalian cells. The ease of gene delivery with the built-in reporters and inducible element enables live cell monitoring, drug selection and temporal knockout, broadening applications of genome editing.

     

Cloning of fatty acid desaturase-coding sequence (<Emphasis Type="Italic">Lufad3</Emphasis>) from flax and its functional validation in rice

Flax contains very high levels of alpha-linolenic acid (57%) and a fatty acid desaturase 3-coding sequence (Lufad3) of flax has been amplified from the RNA isolated from developing seeds. The deduced amino acid sequence of LuFAD3 showed the presence of three histidine motifs, six membrane spanning domains and an endoplasmic reticulum (ER) retention signal KSK, indicating its plausible localization into the ER. Flax is not amenable for genetic transformation and not suitable for functional validation of Lufad3 gene. Hence, rice with well-developed genetic transformation has been selected as heterologous host system. Coding sequence of Lufad3 driven by maize Ubi1 promoter has been introduced into indica rice by Agrobacterium tumefaciens-mediated genetic transformation. Southern analysis of putative transformants (T ₀) revealed signals at variable lengths at >5 kbp indicating random integration of transgene into the genomes of different transformants. The Mendelian segregation observed for selectable marker gene hyg in both T₁ and T₂ generations confirmed stable inheritance and single-site integration of transgenes. As compared to untransformed control (UC), homozygous transgenic rice expressing Lufad3 showed higher levels of essential α-linolenic acid in leaves and seeds validating its functionality.     

Selection of reference genes for quantitative real-time PCR in <Emphasis Type="Italic">Casuarina equisetifolia</Emphasis> under salt stress

Real time quantitative PCR (qPCR) is widely used in gene expression analysis for its accuracy and sensitivity. Reference genes serving as endogenous controls are necessary for gene normalization. In order to select an appropriate reference gene to normalize gene expression in Casuarina equisetifolia under salt stress, 10 potential reference genes were evaluated using real time qPCR in the leaves and roots of plants grown under different NaCl concentrations and treatment durations. GeNorm, NormFinder, and BestKeeper analyses reveal that elongation factor 1-alpha (EF1α) and ubiquitin-conjugating enzyme E2 (UBC) were the most appropriate reference genes for real time qPCR under salt stress. However, β-tubulin (βTUB) and actin 7, which were widely used as reference genes in other plant species, were not always stably expressed. The combination of EF1α, UBC, uncharacterized protein 2, DNAJ homolog subfamily A member 2, and glyceraldehyde-3-phosphate dehydrogenase should be ideal reference genes for normalizing gene expression data in all samples under salt stress. It indicates the need for reference gene selection for normalizing gene expression in C. equisetifolia. In addition, the suitability of reference genes selected was confirmed by validating the expression of WRKY29-like and expansin-like B1. The results enable analysis of salt response mechanism and gene expression in C. equisetifolia.     

Development of cost-effective single nucleotide polymorphism marker assays for genetic diversity analysis in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Competitive allele-specific PCR (KASPar) assay is a user-friendly system that provides flexibility in the numbers of single nucleotide polymorphisms (SNPs) and genotypes. Based on Illumina-GA-IIx genomic data from 10 genotypes with a broad genetic background, 3183 SNPs were selected for KASPar assays development, and 568 were finally converted and selected for Brassica rapa germplasm characterization (17.8%) on the basis of reproducibility, missing data rate, and uniform genetic distribution. High levels of polymorphism of these markers across 231 B. rapa genotypes were verified, illustrating by high polymorphic information content (averaged 0.34), minor allele frequency (0.37), genetic diversity (0.45), and the low observed heterozygosity (0.10). Based on the SNP dataset, structure and principal coordinates analysis, and neighbor-joining phylogenetic methods were used to examine the population structure and gave highly consistent results. The 231 accessions were divided into the four primary subspecies, representing 99 accessions from B. rapa ssp. pekinensis, 85 from B. rapa ssp. chinensis, 30 from B. rapa ssp. rapifera, and 17 from B. rapa ssp. oleifera and were further subdivided into 12 lower-order clusters according to different morphotypes. The genetic variability and pairwise fixation index analysis revealed that the ssp. pekinensis accessions possess the most extensive genetic variation among the four subspecies. The KASPar system is highly useful for validating SNPs and will be valuable for genetics research and breeding applications in B. rapa.     

Weighted set enrichment of gene expression data

Background

Sets of genes that are known to be associated with each other can be used to interpret microarray data. This gene set approach to microarray data analysis can illustrate patterns of gene expression which may be more informative than analyzing the expression of individual genes. Various statistical approaches exist for the analysis of gene sets. There are three main classes of these methods: over-representation analysis, functional class scoring, and pathway topology based methods.

Methods

We propose weighted hypergeometric and weighted chi-squared methods in order to assign a rank to the degree to which each gene participates in the enrichment. Each gene is assigned a weight determined by the absolute value of its log fold change, which is then raised to a certain power. The power value can be adjusted as needed. Datasets from the Gene Expression Omnibus are used to test the method. The significantly enriched pathways are validated through searching the literature in order to determine their relevance to the dataset.

Results

Although these methods detect fewer significantly enriched pathways, they can potentially produce more relevant results. Furthermore, we compare the results of different enrichment methods on a set of microarray studies all containing data from various rodent neuropathic pain models.

Discussion

Our method is able to produce more consistent results than other methods when evaluated on similar datasets. It can also potentially detect relevant pathways that are not identified by the standard methods. However, the lack of biological ground truth makes validating the method difficult.