Genomic imbalance determines positive and negative modulation of gene expression in diploid maize

Genomic imbalance caused by changing the dosage of individual chromosomes (aneuploidy) has a more detrimental effect than varying the dosage of complete sets of chromosomes (ploidy). We examined the impact of both increased and decreased dosage of 15 distal and 1 interstitial chromosomal regions via RNA-seq of maize (Zea mays) mature leaf tissue to reveal new aspects of genomic imbalance. The results indicate that significant changes in gene expression in aneuploids occur both on the varied chromosome (cis) and the remainder of the genome (trans), with a wider spread of modulation compared with the whole-ploidy series of haploid to tetraploid. In general, cis genes in aneuploids range from a gene-dosage effect to dosage compensation, whereas for trans genes the most common effect is an inverse correlation in that expression is modulated toward the opposite direction of the varied chromosomal dosage, although positive modulations also occur. Furthermore, this analysis revealed the existence of increased and decreased effects in which the expression of many genes under genome imbalance are modulated toward the same direction regardless of increased or decreased chromosomal dosage, which is predicted from kinetic considerations of multicomponent molecular interactions. The findings provide novel insights into understanding mechanistic aspects of gene regulation.

Genomic imbalance caused by the addition or subtraction of chromosomal segments leads to modulations of gene expression inversely or directly related to chromosomal dosage but also to nonlinear responses.     

High-resolution array copy number analyses for detection of deletion,gain, amplification and copy-neutral LOH in primary neuroblastoma tumors: Four cases of homozygous deletions of the CDKN2A gene

Background

Neuroblastoma is a very heterogeneous pediatric tumor of the sympathetic nervous system showing clinically significant patterns of genetic alterations. Favorable tumors usually have near-triploid karyotypes with few structural rearrangements. Aggressive stage 4 tumors often have near-diploid or near-tetraploid karyotypes and structural rearrangements. Whole genome approaches for analysis of genome-wide copy number have been used to analyze chromosomal abnormalities in tumor samples. We have used array-based copy number analysis using oligonucleotide single nucleotide polymorphisms (SNP) arrays to analyze the chromosomal structure of a large number of neuroblastoma tumors of different clinical and biological subsets.

Results

Ninety-two neuroblastoma tumors were analyzed with 50 K and/or 250 K SNP arrays from Affymetrix, using CNAG3.0 software. Thirty percent of the tumors harbored 1p deletion, 22% deletion of 11q, 26% had MYCN amplification and 45% 17q gain. Most of the tumors with 1p deletion were found among those with MYCN amplification. Loss of 11q was most commonly seen in tumors without MYCN amplification. In the case of MYCN amplification, two types were identified. One type displayed simple continuous amplicons; the other type harbored more complex rearrangements. MYCN was the only common gene in all cases with amplification. Complex amplification on chromosome 12 was detected in two tumors and three different overlapping regions of amplification were identified. Two regions with homozygous deletions, four cases with CDKN2A deletions in 9p and one case with deletion on 3p (the gene RBMS3) were also detected in the tumors.

Conclusion

SNP arrays provide useful tools for high-resolution characterization of significant chromosomal rearrangements in neuroblastoma tumors. The mapping arrays from Affymetrix provide both copy number and allele-specific information at a resolution of 10–12 kb. Chromosome 9p, especially the gene CDKN2A, is subject to homozygous (four cases) and heterozygous deletions (five cases) in neuroblastoma tumors.

     

Coexistence of diploid, triploid and tetraploid crucian carp (Carassius auratus) in natural waters

Background

Crucian carp (abbreviated CC) belongs to the genus of Carassius within the family of Cyprinidae. It has been one of the most important freshwater species for Chinese aquaculture and is especially abundant in the Dongting water system of Hunan province. CC used to be considered as all diploid forms. However, coexistence of diploid (abbreviated 2nCC), triploid (abbreviated 3nCC) and tetraploid crucian carp (abbreviated 4nCC) population of the Dongting water system was first found by our recently researches.

Results

We examined the ploidy level and compared biological characteristics in different ploidy CC. In reproductive mode, 2nCC was bisexual generative and 4nCC generated all-female offspring by gynogenesis. However, 3nCC generated progenies in two different ways. 3nCC produced bisexual triploid offspring fertilized with 3nCC spermatozoa, while it produced all-female triploid offspring by gynogenesis when its ova were activated by heterogenous spermatozoa. The complete mitochondrial DNA of three different ploidy fishes was sequenced and analyzed, suggesting no significant differences. Interestingly, microchromosomes were found only in 3nCC, which were concluded to be the result of hybridization. Allogenetic DNA fragments of Sox genes were obtained in 3nCC and 4nCC, which were absent in 2nCC. Phylogenetics analysis based on Sox4 gene indicated 3nCC and 4nCC formed a separate group from 2nCC.

Conclusions

In summary, this is the first report of the co-existence of three types of different ploidy crucian carps in natural waters in China. It was proved that the coexistence of different ploidy CC was reproductively maintained. We further hypothesized that 3nCC and 4nCC were allopolyploids that resulted from hybridization. The different ploidy CC population we obtained in this study possesses great significance for the study of polyploidization and the evolution of vertebrates.     

A proof of the DBRF-MEGN method, an algorithm for deducing minimum equivalent gene networks

Background

We previously developed the DBRF-MEGN (difference-based regulation finding-minimum equivalent gene network) method, which deduces the most parsimonious signed directed graphs (SDGs) consistent with expression profiles of single-gene deletion mutants. However, until the present study, we have not presented the details of the method's algorithm or a proof of the algorithm.

Results

We describe in detail the algorithm of the DBRF-MEGN method and prove that the algorithm deduces all of the exact solutions of the most parsimonious SDGs consistent with expression profiles of gene deletion mutants.

Conclusions

The DBRF-MEGN method provides all of the exact solutions of the most parsimonious SDGs consistent with expression profiles of gene deletion mutants.     

Meiotic chromosome pairing behaviour of natural tetraploids and induced autotetraploids of Actinidia chinensis

Key message

Non-preferential chromosome pairing was identified in tetraploid Actinidia chinensis and a higher mean multivalent frequency in pollen mother cells was found in colchine-induced tetraploids of A. chinensis compared with naturally occurring tetraploids.

Abstract

Diploid and tetraploid Actinidia chinensis are used for the development of kiwifruit cultivars. Diploid germplasm can be exploited in a tetraploid breeding programme via unreduced (2n) gametes and chemical-induced chromosome doubling of diploid cultivars and selections. Meiotic chromosome behaviour in diploid A. chinensis ‘Hort16A’ and colchicine-induced tetraploids from ‘Hort16A’ was analysed and compared with that in a diploid male and tetraploid males of A. chinensis raised from seeds sourced from the wild in China. Both naturally occurring and induced tetraploids formed multivalents, but colchicine-induced tetraploids showed a higher mean multivalent frequency in the pollen mother cells. Lagging chromosomes at anaphase I and II were observed at low frequencies in the colchicine-induced tetraploids. To investigate whether preferential or non-preferential chromosome pairing occurs in tetraploid A. chinensis, the inheritance of microsatellite alleles was analysed in the tetraploid progeny of crosses between A. chinensis (4x) and A. arguta (4x). The frequencies of inherited microsatellite allelic combinations in the hybrids suggested that non-preferential chromosome pairing had occurred in the tetraploid A. chinensis parent.     

The Vaccinium corymbosum FLOWERING LOCUS T-like gene (VcFT): a flowering activator reverses photoperiodic and chilling requirements in blueberry

Key message

The blueberry FLOWERING LOCUS T ( FT )-like gene ( VcFT ) cloned from the cDNA of a tetraploid, northern highbush blueberry ( Vaccinium corymbosum L.) is able to reverse the photoperiodic and chilling requirements and drive early and continuous flowering.

Abstract

Blueberry is a woody perennial bush with a longer juvenile period than annual crops, requiring vernalization to flower normally. Few studies have been reported on the molecular mechanism of flowering in blueberry or other woody plants. Because FLOWERING LOCUS T (FT) from Arabidopsis thaliana plays a multifaceted role in generating mobile molecular signals to regulate plant flowering time, isolation and functional analysis of the blueberry (Vaccinium corymbosum L.) FT-like gene (VcFT) will facilitate the elucidation of molecular mechanisms of flowering in woody plants. Based on EST sequences, a 525-bpVcFT was identified and cloned from the cDNA of a tetraploid, northern highbush blueberry cultivar, Bluecrop. Ectopic expression of 35S:VcFT in tobacco induced flowering an average of 28 days earlier than wild-type plants. Expression of the 35S:VcFT in the blueberry cultivar Aurora resulted in an extremely early flowering phenotype, which flowered not only during in vitro culture, a growth stage when nontransgenic shoots had not yet flowered, but also in 6–10-week old, soil-grown transgenic plants, in contrast to the fact that at least 1 year and 800 chilling hours are required for the appearance of the first flower of both nontransgenic ‘Aurora’ and transgenic controls with the gusA. These results demonstrate that the VcFT is a functional floral activator and overexpression of the VcFT is able to reverse the photoperiodic and chilling requirements and drive early and continuous flowering.     

Chromosomal distribution of H3K4me2, H3K9me2 and 5-methylcytosine: variations associated with polyploidy and hybridization in <Emphasis Type="Italic">Brachiaria</Emphasis> (Poaceae)

Key message

Assessment of chromosomal distribution of modified histones and 5-methylcytosine shown that there are diversification of chromosomal types among species of Brachiaria and its interspecific hybrids.

Abstract

Histone post-translational modifications and DNA methylation are epigenetic processes that are involved in structural and functional organization of the genome. This study compared the chromosomal distribution of modified histones and 5-methylcytosine (5-mCyt) in species and interspecific hybrids of Brachiaria with different ploidy levels and reproduction modes. The relation between H3K9me2 and 5-mCyt was observed in the nucleolus organizer region, centromeric central domain and pericentromeric region. H3K4me2 was detected in euchromatic domains, mainly in the terminal chromosomal regions. Comparison of chromosomal distribution among species and hybrids showed greater variation of chromosomal types for the H3K9me2 in B. decumbens (tetraploid and apomictic species) and the 963 hybrid, while, for the H3K4me2, the variation was higher in B. brizantha and B. decumbens (tetraploid and apomictic species) and 963 hybrid. The chromosome distribution of 5-mCyt was similar between B. brizantha and B. decumbens, which differ from the distribution observed in B. ruziziensis (diploid and sexual species). Significant alterations in DNA methylation were observed in the artificially tetraploidized B. ruziziensis and in the interspecific hybrids, possibly as result of hybridization and polyploidization processes. The monitoring of histone modifications and DNA methylation allowed categorizing nuclear and chromosomal distribution of these epigenetic marks, thus contributing to the knowledge of composition and structure of the genome/epigenome of Brachiaria species and hybrids. These data can be useful for speciation and genome evolution studies in genus Brachiaria, and represent important markers to explore relationships between genomes.

     

Fitting hidden Markov models of protein domains to a target species: application to Plasmodium falciparum

Background

The identification of gene sets that are significantly impacted in a given condition based on microarray data is a crucial step in current life science research. Most gene set analysis methods treat genes equally, regardless how specific they are to a given gene set.

Results

In this work we propose a new gene set analysis method that computes a gene set score as the mean of absolute values of weighted moderated gene t-scores. The gene weights are designed to emphasize the genes appearing in few gene sets, versus genes that appear in many gene sets. We demonstrate the usefulness of the method when analyzing gene sets that correspond to the KEGG pathways, and hence we called our method P athway A nalysis with D own-weighting of O verlapping G enes (PADOG). Unlike most gene set analysis methods which are validated through the analysis of 2-3 data sets followed by a human interpretation of the results, the validation employed here uses 24 different data sets and a completely objective assessment scheme that makes minimal assumptions and eliminates the need for possibly biased human assessments of the analysis results.

Conclusions

PADOG significantly improves gene set ranking and boosts sensitivity of analysis using information already available in the gene expression profiles and the collection of gene sets to be analyzed. The advantages of PADOG over other existing approaches are shown to be stable to changes in the database of gene sets to be analyzed. PADOG was implemented as an R package available at: http://bioinformaticsprb.med.wayne.edu/PADOG/or http://www.bioconductor.org.     

Hybridization success is largely limited to homoploid <Emphasis Type="Italic">Prunus</Emphasis> hybrids: a multidisciplinary approach

Prunus fruticosa is a rare shrub occurring in Eurasian thermophilous forest-steppe alliances. The species frequently hybridizes with cultivated Prunus species in Europe (allochthonous tetraploid P. cerasus and partly indigenous diploid P. avium). Propidium iodide flow cytometry, distance-based morphometrics, elliptic Fourier analysis and embryology were employed to evaluate the extent of hybridization in six Slovak populations. Flow cytometric analyses revealed three ploidy levels: diploid (P. avium), triploid (P. × mohacsyana) and tetraploid (P. fruticosa, P. × eminens and P. cerasus). In addition, P. fruticosa and P. cerasus, at the tetraploid level, were found to differ in absolute genome size. An embryological evaluation suggested the existence of a triploid block in P. × mohacsyana and significant potential for hybridization among tetraploid taxa (indicated also by a continuous distribution of genome size data and further mirrored by morphometrics). Although hybrids significantly differ in ploidy level and embryological characteristics, they are almost indistinguishable using morphological characters. Hybridization with P. cerasus thus turns out to be a significant threat to wild populations of P. fruticosa compared to the relatively weak influence of P. avium.     

Rate of red blood cell destruction varies in different strains of mice infected with Plasmodium berghei-ANKA after chronic exposure

Background

The ATP-binding cassette (ABC) superfamily is one of the largest evolutionarily conserved families of proteins. ABC proteins play key roles in cellular detoxification of endobiotics and xenobiotics. Overexpression of certain ABC proteins, among them the multidrug resistance associated protein (MRP), contributes to drug resistance in organisms ranging from human neoplastic cells to parasitic protozoa. In the present study, the Plasmodium berghei mrp gene (pbmrp) was partially characterized and the predicted protein was classified using bioinformatics in order to explore its putative involvement in drug resistance.

Methods

The pbmrp gene from the P. berghei drug sensitive, N clone, was sequenced using a PCR strategy. Classification and domain organization of pbMRP were determined with bioinformatics. The Plasmodium spp. MRPs were aligned and analysed to study their conserved motifs and organization. Gene copy number and organization were determined via Southern blot analysis in both N clone and the chloroquine selected line, RC. Chromosomal Southern blots and RNase protection assays were employed to determine the chromosomal location and expression levels of pbmrp in blood stages.

Results

The pbmrp gene is a single copy, intronless gene with a predicted open reading frame spanning 5820 nucleotides. Bioinformatic analyses show that this protein has distinctive features characteristic of the ABCC sub-family. Multiple sequence alignments reveal a high degree of conservation in the nucleotide binding and transmembrane domains within the MRPs from the Plasmodium spp. analysed. Expression of pbmrp was detected in asexual blood stages. Gene organization, copy number and mRNA expression was similar in both lines studied. A chromosomal translocation was observed in the chloroquine selected RC line, from chromosome 13/14 to chromosome 8, when compared to the drug sensitive N clone.

Conclusion

In this study, the pbmrp gene was sequenced and classified as a member of the ABCC sub-family. Multiple sequence alignments reveal that this gene is homologous to the Plasmodium y. yoelii and Plasmodium knowlesi mrp, and the Plasmodium vivax and Plasmodium falciparum mrp2 genes. There were no differences in gene organization, copy number, or mRNA expression between N clone and the RC line, but a chromosomal translocation of pbmrp from chromosome 13/14 to chromosome 8 was detected in RC.     

The CIN4 Chromosomal Instability qPCR Classifier Defines Tumor Aneuploidy and Stratifies Outcome in Grade 2 Breast Cancer

Purpose

Quantifying chromosomal instability (CIN) has both prognostic and predictive clinical utility in breast cancer. In order to establish a robust and clinically applicable gene expression-based measure of CIN, we assessed the ability of four qPCR quantified genes selected from the 70-gene Chromosomal Instability (CIN70) expression signature to stratify outcome in patients with grade 2 breast cancer.

Methods

AURKA, FOXM1, TOP2A and TPX2 (CIN4), were selected from the CIN70 signature due to their high level of correlation with histological grade and mean CIN70 signature expression in silico. We assessed the ability of CIN4 to stratify outcome in an independent cohort of patients diagnosed between 1999 and 2002. 185 formalin-fixed, paraffin-embedded (FFPE) samples were included in the qPCR measurement of CIN4 expression. In parallel, ploidy status of tumors was assessed by flow cytometry. We investigated whether the categorical CIN4 score derived from the CIN4 signature was correlated with recurrence-free survival (RFS) and ploidy status in this cohort.

Results

We observed a significant association of tumor proliferation, defined by Ki67 and mitotic index (MI), with both CIN4 expression and aneuploidy. The CIN4 score stratified grade 2 carcinomas into good and poor prognostic cohorts (mean RFS: 83.8±4.9 and 69.4±8.2 months, respectively, p = 0.016) and its predictive power was confirmed by multivariate analysis outperforming MI and Ki67 expression.

Conclusions

The first clinically applicable qPCR derived measure of tumor aneuploidy from FFPE tissue, stratifies grade 2 tumors into good and poor prognosis groups.     

Assessing unmodified 70-mer oligonucleotide probe performance on glass-slide microarrays

Background

Long oligonucleotide microarrays are potentially more cost- and management-efficient than cDNA microarrays, but there is little information on the relative performance of these two probe types. The feasibility of using unmodified oligonucleotides to accurately measure changes in gene expression is also unclear.

Results

Unmodified sense and antisense 70-mer oligonucleotides representing 75 known rat genes and 10 Arabidopsis control genes were synthesized, printed and UV cross-linked onto glass slides. Printed alongside were PCR-amplified cDNA clones corresponding to the same genes, enabling us to compare the two probe types simultaneously. Our study was designed to evaluate the mRNA profiles of heart and brain, along with Arabidopsis cRNA spiked into the labeling reaction at different relative copy number. Hybridization signal intensity did not correlate with probe type but depended on the extent of UV irradiation. To determine the effect of oligonucleotide concentration on hybridization signal, 70-mers were serially diluted. No significant change in gene-expression ratio or loss in hybridization signal was detected, even at the lowest concentration tested (6.25 μm). In many instances, signal intensity actually increased with decreasing concentration. The correlation coefficient between oligonucleotide and cDNA probes for identifying differentially expressed genes was 0.80, with an average coefficient of variation of 13.4%. Approximately 8% of the genes showed discordant results with the two probe types, and in each case the cDNA results were more accurate, as determined by real-time PCR.

Conclusions

Microarrays of UV cross-linked unmodified oligonucleotides provided sensitive and specific measurements for most of the genes studied.