Coordination of gene expression between organellar and nuclear genomes

Following the acquisition of chloroplasts and mitochondria by eukaryotic cells during endosymbiotic evolution, most of the genes in these organelles were either lost or transferred to the nucleus. Encoding organelle-destined proteins in the nucleus allows for host control of the organelle. In return, organelles send signals to the nucleus to coordinate nuclear and organellar activities. In photosynthetic eukaryotes, additional interactions exist between mitochondria and chloroplasts. Here we review recent advances in elucidating the intracellular signalling pathways that coordinate gene expression between organelles and the nucleus, with a focus on photosynthetic plants.

Transposable elements and the evolution of regulatory networks

The control and coordination of eukaryotic gene expression rely on transcriptional and post-transcriptional regulatory networks. Although progress has been made in mapping the components and deciphering the function of these networks, the mechanisms by which such intricate circuits originate and evolve remain poorly understood. Here I revisit and expand earlier models and propose that genomic repeats, and in particular transposable elements, have been a rich source of material for the assembly and tinkering of eukaryotic gene regulatory systems.

Population genetic inference using a fixed number of segregating sites: a reassessment

Coalescent theory is commonly used to perform population genetic inference at the nucleotide level. Here, we examine the procedure that fixes the number of segregating sites (henceforth the FS procedure). In this approach a fixed number of segregating sites (S) are placed on a coalescent tree (independently of the total and internode lengths of the tree). Thus, although widely used, the FS procedure does not strictly follow the assumptions of coalescent theory and must be considered an approximation of (i) the standard procedure that uses a fixed population mutation parameter theta, and (ii) procedures that condition on the number of segregating sites. We study the differences in the false positive rate for nine statistics by comparing the FS procedure with the procedures (i) and (ii), using several evolutionary models with single-locus and multilocus data. Our results indicate that for single-locus data the FS procedure is accurate for the equilibrium neutral model, but problems arise under the alternative models studied; furthermore, for multilocus data, the FS procedure becomes inaccurate even for the standard neutral model. Therefore, we recommend a procedure that fixes the theta value (or alternatively, procedures that condition on S and take into account the uncertainty of theta) for analysing evolutionary models with multilocus data. With single-locus data, the FS procedure should not be employed for models other than the standard neutral model.

Microbiology in the post-genomic era

Genomics has revolutionized every aspect of microbiology. Now, 13 years after the first bacterial genome was sequenced, it is important to pause and consider what has changed in microbiology research as a consequence of genomics. In this article, we review the evolving field of bacterial typing and the genomic technologies that enable comparative analysis of multiple genomes and the metagenomes of complex microbial environments, and address the implications of the genomic era for the future of microbiology.

Streptococcus pyogenes pSM19035 requires dynamic assembly of ATP-bound ParA and ParB on parS DNA during plasmid segregation

The accurate partitioning of Firmicute plasmid pSM19035 at cell division depends on ATP binding and hydrolysis by homodimeric ATPase δ₂ (ParA) and binding of ω₂ (ParB) to its cognate parS DNA. The 1.83 Å resolution crystal structure of δ₂ in a complex with non-hydrolyzable ATPγS reveals a unique ParA dimer assembly that permits nucleotide exchange without requiring dissociation into monomers. In vitro, δ₂ had minimal ATPase activity in the absence of ω₂ and parS DNA. However, stoichiometric amounts of ω₂ and parS DNA stimulated the δ₂ ATPase activity and mediated plasmid pairing, whereas at high (4:1) ω₂ : δ₂ ratios, stimulation of the ATPase activity was reduced and δ₂ polymerized onto DNA. Stimulation of the δ₂ ATPase activity and its polymerization on DNA required ability of ω₂ to bind parS DNA and its N-terminus. In vivo experiments showed that δ₂ alone associated with the nucleoid, and in the presence of ω₂ and parS DNA, δ₂ oscillated between the nucleoid and the cell poles and formed spiral-like structures. Our studies indicate that the molar ω₂ : δ₂ ratio regulates the polymerization properties of (δ•ATP•Mg²⁺)₂ on and depolymerization from parS DNA, thereby controlling the temporal and spatial segregation of pSM19035 before cell division.

Global analysis of aberrant pre-mRNA splicing in glioblastoma using exon expression arrays

Background

Tumor-predominant splice isoforms were identified during comparative in silico sequence analysis of EST clones, suggesting that global aberrant alternative pre-mRNA splicing may be an epigenetic phenomenon in cancer. We used an exon expression array to perform an objective, genome-wide survey of glioma-specific splicing in 24 GBM and 12 nontumor brain samples. Validation studies were performed using RT-PCR on glioma cell lines, patient tumor and nontumor brain samples.

Results

In total, we confirmed 14 genes with glioma-specific splicing; seven were novel events identified by the exon expression array (A2BP1, BCAS1, CACNA1G, CLTA, KCNC2, SNCB, and TPD52L2). Our data indicate that large changes (> 5-fold) in alternative splicing are infrequent in gliomagenesis (< 3% of interrogated RefSeq entries). The lack of splicing changes may derive from the small number of splicing factors observed to be aberrantly expressed.

Conclusion

While we observed some tumor-specific alternative splicing, the number of genes showing exclusive tumor-specific isoforms was on the order of tens, rather than the hundreds suggested previously by in silico mining. Given the important role of alternative splicing in neural differentiation, there may be selective pressure to maintain a majority of splicing events in order to retain glial-like characteristics of the tumor cells.

Cuticular hydrocarbon analysis of an awake behaving fly using direct analysis in real-time time-of-flight mass spectrometry

In mammals and insects, pheromones strongly influence social behaviors such as aggression and mate recognition. In Drosophila melanogaster, pheromones in the form of cuticular hydrocarbons play prominent roles in courtship. GC/MS is the primary analytical tool currently used to study Drosophila cuticular hydrocarbons. Although GC/MS is highly reproducible and sensitive, it requires that the fly be placed in a lethal solution of organic solvent, thereby impeding further behavioral studies. We present a technique for the analysis of hydrocarbons and other surface molecules from live animals by using direct analysis in real-time (DART) MS. Cuticular hydrocarbons were sampled from the surface of a restrained, awake behaving fly by using several brief, carefully controlled depressions of the abdomen with a small steel probe. DART mass spectral analysis of the probe detected ions with mass-to-charge ratio (m/z) of the protonated molecule corresponding to many of the previously identified unsaturated hydrocarbons. Six additional cuticular hydrocarbons also were identified. Consistent with previous GC/MS studies, male and female differences in chemical composition were evident. Spatial differences in the expression profile also were observed on males. Sampling from an individual female first as a virgin and then 45 and 90 min after successful copulation showed that mass signals likely to correspond to cis-vaccenyl acetate, tricosene, and pentacosene increased in relative intensity after courtship. This method provides near-instantaneous analysis of an individual animal's chemical profile in parallel with behavioral studies and could be extended to other models of pheromone-mediated behavior.

Malin: maximum likelihood analysis of intron evolution in eukaryotes

Summary: Malin is a software package for the analysis of eukaryotic gene structure evolution. It provides a graphical user interface for various tasks commonly used to infer the evolution of exon–intron structure in protein-coding orthologs. Implemented tasks include the identification of conserved homologous intron sites in protein alignments, as well as the estimation of ancestral intron content, lineage-specific intron losses and gains. Estimates are computed either with parsimony, or with a probabilistic model that incorporates rate variation across lineages and intron sites.

Availability: Malin is available as a stand-alone Java application, as well as an application bundle for MacOS X, at the website http://www.iro.umontreal.ca/~csuros/introns/malin/. The software is distributed under a BSD-style license.

Contact:csuros@iro.umontreal.ca

Apical loop-internal loop RNA pseudoknots: a new type of stimulator of -1 translational frameshifting in bacteria

Nearly all members of a widespread family of bacterial transposable elements related to insertion sequence 3 (IS3), therefore called the IS3 family, very likely use programmed -1 ribosomal frameshifting to produce their transposase, a protein required for mobility. Comparative analysis of the potential frameshift signals in this family suggested that most of the insertion sequences from the IS51 group contain in their mRNA an elaborate pseudoknot that could act as a recoding stimulator. It results from a specific intramolecular interaction between an apical loop and an internal loop from two stem-loop structures. Directed mutagenesis, chemical probing, and gel mobility assays of the frameshift region of one element from the IS51 group, IS3411, provided clear evidences of the existence of the predicted structure. Modeling was used to generate a three-dimensional molecular representation of the apical loop-internal loop complex. We could demonstrate that mutations affecting the stability of the structure reduce both frameshifting and transposition, thus establishing the biological importance of this new type of RNA structure for the control of transposition level.

Histone modifications and nuclear architecture: a review

Epigenetic modifications, such as acetylation, phosphorylation, methylation, ubiquitination, and ADP ribosylation, of the highly conserved core histones, H2A, H2B, H3, and H4, influence the genetic potential of DNA. The enormous regulatory potential of histone modification is illustrated in the vast array of epigenetic markers found throughout the genome. More than the other types of histone modification, acetylation and methylation of specific lysine residues on N-terminal histone tails are fundamental for the formation of chromatin domains, such as euchromatin, and facultative and constitutive heterochromatin. In addition, the modification of histones can cause a region of chromatin to undergo nuclear compartmentalization and, as such, specific epigenetic markers are non-randomly distributed within interphase nuclei. In this review, we summarize the principles behind epigenetic compartmentalization and the functional consequences of chromatin arrangement within interphase nuclei.

Dimorphism in fungal plant pathogens

Fungi are mostly sessile organisms, and thus have evolved ways to cope with environmental changes. Many fungi produce 'dormant' structures, which allow them to survive periods of unfavorable conditions. Another ingenious active approach to a changing environment has been adopted by the 'dimorphic fungi', which simply shift their thallic organization as a way to adapt and thrive in the new conditions. Dimorphism is extensively exploited by both plant and animal pathogenic fungi, where the encounter with the host prompts a shift in the mode of growth. In this review, we focus on the phenomenon of dimorphism among plant pathogenic fungi through discussion of several relatively well-studied exemplar species.

Review. Meiotic drive and sex determination: molecular and cytological mechanisms of sex ratio adjustment in birds

Differences in relative fitness of male and female offspring across ecological and social environments should favour the evolution of sex-determining mechanisms that enable adjustment of brood sex ratio to the context of breeding. Despite the expectation that genetic sex determination should not produce consistent bias in primary sex ratios, extensive and adaptive modifications of offspring sex ratio in relation to social and physiological conditions during reproduction are often documented. Such discordance emphasizes the need for empirical investigation of the proximate mechanisms for modifying primary sex ratios, and suggests epigenetic effects on sex-determining mechanisms as the most likely candidates. Birds, in particular, are thought to have an unusually direct opportunity to modify offspring sex ratio because avian females are heterogametic and because the sex-determining division in avian meiosis occurs prior to ovulation and fertilization. However, despite evidence of strong epigenetic effects on sex determination in pre-ovulatory avian oocytes, the mechanisms behind such effects remain elusive. Our review of molecular and cytological mechanisms of avian meiosis uncovers a multitude of potential targets for selection on biased segregation of sex chromosomes, which may reflect the diversity of mechanisms and levels on which such selection operates in birds. Our findings indicate that pronounced differences between sex chromosomes in size, shape, size of protein bodies, alignment at the meiotic plate, microtubule attachment and epigenetic markings should commonly produce biased segregation of sex chromosomes as the default state, with secondary evolution of compensatory mechanisms necessary to maintain unbiased meiosis. We suggest that it is the epigenetic effects that modify such compensatory mechanisms that enable context-dependent and precise adjustment of primary sex ratio in birds. Furthermore, we highlight the features of avian meiosis that can be influenced by maternal hormones in response to environmental stimuli and may account for the precise and adaptive patterns of offspring sex ratio adjustment observed in some species.

Convergent mechanisms of genome evolution of large and giant DNA viruses

We have taken advantage of the availability of the genome sequences of a collection of large and giant viruses infecting bacteria (T4 family) and eukaryotes (NCLDV group) to assess some of the evolutionary forces which might have shaped their genomes. Despite having apparently different ancestors, these two groups of viruses are affected by convergent evolutionary forces. Both types of virus probably originated from a simple and ancient viral ancestor with a small subset of 30-35 genes encoding replication and structural proteins. The genome size and diversity of the descendants most likely grew progressively by (i) lineage-specific gene duplications, (ii) lateral gene transfers of cellular genes and (iii) accretion of diverse families of mobile genetic elements. These results argue against the hypothesis that giant viruses derive from a regressive cell.

Changes of heavy metal, metallothionein and heat shock proteins in Sertoli cells induced by cadmium exposure

In this study, we examined the levels of Cadmium (Cd), iron (Fe) and zinc (Zn), which were considered to be involved in Sertoli cell damage caused by Cd exposure. We also examined metallothionein (MT), heat shock protein 70 (Hsp70) and heme oxygenase-1 (HO-1) expressions in Sertoli cells induced by Cd exposure. Evaluation by the in-air micro-particle induced X-ray emission (PIXE) method revealed that Cd and Fe distribution was increased in the cytoplasm of Sertoli cells after Cd exposure. By contrast, Zn was decreased in the cytoplasm of Sertoli cells after Cd exposure. It was suggested that the target of Cd toxicity was the cytoplasm of Sertoli cells, Fe was considered to enhance damage to Sertoli cells caused by Cd exposure. The DNA fragmentation rate was determined by ELISA after Cd exposure to Sertoli cells. It remained essentially unchanged with 2.5 microM Cd exposure of Sertoli cells; however, MT, Hsp70 and HO-1 were significantly increased by Cd exposure. As a result, Cd-induced MT was protected Sertoli cells against apoptosis, and Cd-induced HO-1 was involved in protection against oxidative stress. Incidentally, MT, Hsp70 and HO-1 showed similar responses to Cd exposure.

The primary open-angle glaucoma gene WDR36 functions in ribosomal RNA processing and interacts with the p53 stress-response pathway

Primary open-angle glaucoma (POAG) is a genetically complex neuropathy that affects retinal ganglion cells and is a leading cause of blindness worldwide. WDR36, a gene of unknown function, was recently identified as causative for POAG at locus GLC1G. Subsequent studies found disease-associated variants in control populations, leaving the role of WDR36 in this disease unclear. To address this issue, we determined the function of WDR36. We studied Wdr36 in zebrafish and found it is the functional homolog of yeast Utp21. Utp21 is cell essential and functions in the nucleolar processing of 18S rRNA, which is required for ribosome biogenesis. Evidence for functional homology comes from sequence alignment, ubiquitous expression, sub-cellular localization to the nucleolus and loss-of-function phenotypes that include defects in 18S rRNA processing and abnormal nucleolar morphology. Additionally, we show that loss of Wdr36 function leads to an activation of the p53 stress-response pathway, suggesting that co-inheritance of defects in p53 pathway genes may influence the impact of WDR36 variants on POAG. Although these results overall do not provide evidence for or against a role of WDR36 in POAG, they do provide important baseline information for future studies.

Intraocular pressure elevation induces mitochondrial fission and triggers OPA1 release in glaucomatous optic nerve

PURPOSE

To determine whether elevation of intraocular pressure (IOP) triggers mitochondrial fission and ultrastructural changes and alters optic atrophy type 1 (OPA1) expression and distribution in the optic nerve (ON) of glaucomatous DBA/2J mice.

METHODS

IOP in the eyes of DBA/2J mice was measured, and mitochondrial structural changes were assessed by conventional electron microscopy (EM) and EM tomography. Cytochrome c oxidase IV subunit 1 (COX), OPA1, and Dnm1, a rat homologue of dynamin-related protein-1, mRNA were measured by quantitative (q)PCR. COX and OPA1 protein distribution was assessed by immunocytochemistry and Western blot.

RESULTS

Excavation of the optic nerve head (ONH), axon loss, and COX reduction were evident in 10-month-old glaucomatous ONHs of eyes with >20 mm Hg IOP elevation. EM analysis showed mitochondrial fission, matrix swelling, substantially reduced cristae volume, and abnormal cristae depletion in 10-month-old glaucomatous ONH axons. The mean length of mitochondrial cross section in these axons decreased from 858.2 +/- 515.3 nm in 3-month-old mice to 583.3 +/- 298.6 nm in 10-month-old glaucomatous mice (P < 0.001). Moderate reductions of COX mRNA were observed in the 10-month-old DBA/2J mice's ONHs. Larger reductions of OPA1 immunoreactivity and gene expression were coupled with larger increases of Dnm1 gene expression in 10-month-old glaucomatous ONH. Subcellular fractionation analysis indicates increased release of both OPA1 and cytochrome c from mitochondria in 10-month-old glaucomatous ONs.

CONCLUSIONS

IOP elevation may directly damage mitochondria in the ONH axons by promoting reduction of COX, mitochondrial fission and cristae depletion, alterations of OPA1 and Dnm1 expression, and induction of OPA1 release. Thus, interventions to preserve mitochondria may be useful for protecting against ON degeneration in glaucoma.

Triad pattern algorithm for predicting strong promoter candidates in bacterial genomes

Background

Bacterial promoters, which increase the efficiency of gene expression, differ from other promoters by several characteristics. This difference, not yet widely exploited in bioinformatics, looks promising for the development of relevant computational tools to search for strong promoters in bacterial genomes.

Results

We describe a new triad pattern algorithm that predicts strong promoter candidates in annotated bacterial genomes by matching specific patterns for the group I σ⁷⁰ factors of Escherichia coli RNA polymerase. It detects promoter-specific motifs by consecutively matching three patterns, consisting of an UP-element, required for interaction with the α subunit, and then optimally-separated patterns of -35 and -10 boxes, required for interaction with the σ⁷⁰ subunit of RNA polymerase. Analysis of 43 bacterial genomes revealed that the frequency of candidate sequences depends on the A+T content of the DNA under examination. The accuracy of in silico prediction was experimentally validated for the genome of a hyperthermophilic bacterium, Thermotoga maritima, by applying a cell-free expression assay using the predicted strong promoters. In this organism, the strong promoters govern genes for translation, energy metabolism, transport, cell movement, and other as-yet unidentified functions.

Conclusion

The triad pattern algorithm developed for predicting strong bacterial promoters is well suited for analyzing bacterial genomes with an A+T content of less than 62%. This computational tool opens new prospects for investigating global gene expression, and individual strong promoters in bacteria of medical and/or economic significance.

Regulation of the cellular heat shock response in Caenorhabditis elegans by thermosensory neurons

Temperature pervasively affects all cellular processes. In response to a rapid increase in temperature, all cells undergo a heat shock response, an ancient and highly conserved program of stress-inducible gene expression, to reestablish cellular homeostasis. In isolated cells, the heat shock response is initiated by the presence of misfolded proteins and therefore thought to be cell-autonomous. In contrast, we show that within the metazoan Caenorhabditis elegans, the heat shock response of somatic cells is not cell-autonomous but rather depends on the thermosensory neuron, AFD, which senses ambient temperature and regulates temperature-dependent behavior. We propose a model whereby this loss of cell autonomy serves to integrate behavioral, metabolic, and stress-related responses to establish an organismal response to environmental change.

Circadian clock function in Arabidopsis thaliana: time beyond transcription

The past decade has seen a remarkable advance in our understanding of the plant circadian system, mostly in Arabidopsis thaliana. It is now well established that Arabidopsis clock genes and their protein products operate through autoregulatory feedback loops that promote rhythmic oscillations in cellular, metabolic and physiological activities. This article reviews recent studies that have provided evidence for new mechanisms of clock organization and function. These mechanisms include protein-protein interactions and the regulation of protein stability, which, together, directly connect light signalling to the Arabidopsis circadian system. Evidence of rhythmic changes in chromatin structure has also opened new and exciting ways for regulation of clock gene expression. All of these mechanisms ensure an appropriate synchronization with the environment, which is crucial for successful plant growth and development.

An inducible and reversible mouse genetic rescue system

Inducible and reversible regulation of gene expression is a powerful approach for uncovering gene function. We have established a general method to efficiently produce reversible and inducible gene knockout and rescue in mice. In this system, which we named iKO, the target gene can be turned on and off at will by treating the mice with doxycycline. This method combines two genetically modified mouse lines: a) a KO line with a tetracycline-dependent transactivator replacing the endogenous target gene, and b) a line with a tetracycline-inducible cDNA of the target gene inserted into a tightly regulated (TIGRE) genomic locus, which provides for low basal expression and high inducibility. Such a locus occurs infrequently in the genome and we have developed a method to easily introduce genes into the TIGRE site of mouse embryonic stem (ES) cells by recombinase-mediated insertion. Both KO and TIGRE lines have been engineered for high-throughput, large-scale and cost-effective production of iKO mice. As a proof of concept, we have created iKO mice in the apolipoprotein E (ApoE) gene, which allows for sensitive and quantitative phenotypic analyses. The results demonstrated reversible switching of ApoE transcription, plasma cholesterol levels, and atherosclerosis progression and regression. The iKO system shows stringent regulation and is a versatile genetic system that can easily incorporate other techniques and adapt to a wide range of applications.

We describe a technology for the creation of inducible and reversible gene inactivation in mice. It combines two genetically modified mouse lines: a knock-out line with a tetracycline transactivator replacing the endogenous target gene, and a line in which a tetracycline-inducible cDNA of the target gene has been inserted into a specific genomic locus. A critical component of this system is the unique chromosomal loci we have identified and engineered that offer a platform for easy insertion of any gene of interest for tightly controlled expression. Because of its simple binary nature, allowing independent modification of each of the two components and possibility of use in a high-throughput mode, we believe that our system will be useful for multiple applications, such as introducing mutant or humanized form of the target gene as well as functional manipulating tools. We have applied this technology to the Apolipoprotein E (ApoE) gene and have demonstrated that: a) the expression of ApoE is strictly dependent on the presence of doxycycline, a tetracycline group antibiotic, in the mouse diet, b) in the absence of doxycycline (ApoE repressed) atherosclerotic plaques are formed, confirming the importance of ApoE in the process, and c) upon re-induction of ApoE in the animals with doxicyclin, atherosclerosis regressed.

Developmental plasticity of lymphocytes

Experimental perturbation of signaling or transcription factor networks has been used to study the developmental potential of lymphoid progenitors, lineage-committed precursors and mature lymphocytes. Common lymphoid progenitors and uncommitted pro-T cells can be efficiently diverted into myeloid or erythroid lineages by ectopic cytokine signaling or retroviral expression of the myeloid C/EBPalpha or erythroid GATA1 transcription factor. Forced C/EBPalpha expression furthermore induces direct transdifferentiation of immature thymocytes or B cells into macrophages. Notably, conditional inactivation of the B cell commitment factor Pax5 is sufficient to convert mature B cells into functional T cells via dedifferentiation to uncommitted progenitors. Together these experiments have uncovered an unanticipated developmental plasticity of lymphocytes, which may account for lineage switches observed in human malignancies.

Retromer sorting: a pathogenic pathway in late-onset Alzheimer disease

During the tail end of the 20th century, a "golden period" in Alzheimer disease (AD) research, many of the pathogenic molecules of the autosomal dominant form of the disease were isolated. These molecular defects, however, do not exist in "sporadic" late-onset AD, the form of the disease that accounts for more than 95% of all cases. Pinpointing the pathogenic molecules of late-onset AD has, therefore, become an urgent goal, both for understanding disease mechanisms and for opening up novel therapeutic avenues. The retromer sorting pathway transports cargo along the endosome-trans-Golgi network, and retromer defects were first implicated in late-onset AD by a study that combined brain imaging with microarray. A range of studies have confirmed that defects in this pathway can play a pathogenic role in the disease. Herein, these findings will be reviewed, the details of the retromer sorting pathway will be discussed, and a biological model that can account for the disease's regional selectivity will be elaborated.

Reactive nitrogen species-dependent DNA damage in EBV-associated nasopharyngeal carcinoma: the relation to STAT3 activation and EGFR expression

Nasopharyngeal carcinoma (NPC) is strongly associated with Epstein-Barr virus (EBV) infection. Recently, reactive nitrogen and oxygen species are considered to participate in inflammation-related carcinogenesis through DNA damage. In our study, we obtained biopsy and surgical specimens of nasopharyngeal tissues from NPC patients in southern China, and performed double immunofluorescent staining to examine the formation of 8-nitroguanine, a nitrative DNA lesion and 8-oxo-7,8-dihydro-2'-deoxyguanosine, an oxidative DNA lesion, in these specimens. Strong DNA lesions were observed in cancer cells and inflammatory cells in stroma of NPC patients. Intensive immunoreactivity of iNOS was detected in the cytoplasm of 8-nitroguanine-positive cancer cells. DNA lesions and iNOS expression were also observed in epithelial cells of EBV-positive patients with chronic nasopharyngitis, although their intensities were significantly weaker than those in NPC patients. In EBV-negative subjects, no or little DNA lesions and iNOS expression were observed. EGFR and phosphorylated STAT3 were strongly expressed in cancer cells of NPC patients, but NF-kappaB was not expressed, suggesting that STAT3-dependent mechanism is important for NPC carcinogenesis. IL-6 was expressed mainly in inflammatory cells of nasopharyngeal tissues of EBV-infected patients. EBV-encoded RNAs (EBERs) and latent membrane protein 1 (LMP1) were detected in cancer cells from all EBV-infected patients. In vitro cell system, nuclear accumulation of EGFR was observed in LMP1-expressing cells, and IL-6 induced phosphorylated STAT3 and iNOS. These data suggest that nuclear accumulation of EGFR and STAT3 activation by IL-6 play the key role in iNOS expression and resultant DNA damage, leading to EBV-mediated NPC.

Epigenetic Basis for the Differentiation Potential of Mesenchymal and Embryonic Stem Cells

SUMMARY: Stem cells have the ability to self-renew, and give rise to one or more differentiated cell types. Embryonic stem cells can differentiate into all cell types of the body and have unlimited self-renewal capacity. Somatic stem cells are found in many adult tissues. They have an extensive but finite lifespan and can differentiate into a more restricted range of cell types. Increasing evidence indicates that the multilineage differentiation ability of stem cells is defined by the potential for expression of developmentally regulated transcription factors and of lineage specification genes. Gene expression, or as emphasized here, the potential for gene expression, is largely controlled by epigenetic modifications of DNA (DNA methylation) and chromatin (such as post-translational histone modifications) in the regulatory regions of specific genes. Epigenetic modifications can also influence the timing of DNA replication. We highlight here how mechanisms by which genes are poised for transcription in undifferentiated stem cells are being uncovered through the mapping of DNA methylation profiles on differentiation-regulated promoters and at the genome-wide level, histone modifications, and transcription factor binding. Epigenetic marks on developmentally regulated and lineage specification genes in stem cells seem to define a state of pluripotency.

Chloroplast genome sequencing analysis of Heterosigma akashiwo CCMP452 (West Atlantic) and NIES293 (West Pacific) strains

BACKGROUND

Heterokont algae form a monophyletic group within the stramenopile branch of the tree of life. These organisms display wide morphological diversity, ranging from minute unicells to massive, bladed forms. Surprisingly, chloroplast genome sequences are available only for diatoms, representing two (Coscinodiscophyceae and Bacillariophyceae) of approximately 18 classes of algae that comprise this taxonomic cluster. A universal challenge to chloroplast genome sequencing studies is the retrieval of highly purified DNA in quantities sufficient for analytical processing. To circumvent this problem, we have developed a simplified method for sequencing chloroplast genomes, using fosmids selected from a total cellular DNA library. The technique has been used to sequence chloroplast DNA of two Heterosigma akashiwo strains. This raphidophyte has served as a model system for studies of stramenopile chloroplast biogenesis and evolution.

RESULTS

H. akashiwo strain CCMP452 (West Atlantic) chloroplast DNA is 160,149 bp in size with a 21,822-bp inverted repeat, whereas NIES293 (West Pacific) chloroplast DNA is 159,370 bp in size and has an inverted repeat of 21,665 bp. The fosmid cloning technique reveals that both strains contain an isomeric chloroplast DNA population resulting from an inversion of their single copy domains. Both strains contain multiple small inverted and tandem repeats, non-randomly distributed within the genomes. Although both CCMP452 and NIES293 chloroplast DNAs contains 197 genes, multiple nucleotide polymorphisms are present in both coding and intergenic regions. Several protein-coding genes contain large, in-frame inserts relative to orthologous genes in other plastids. These inserts are maintained in mRNA products. Two genes of interest in H. akashiwo, not previously reported in any chloroplast genome, include tyrC, a tyrosine recombinase, which we hypothesize may be a result of a lateral gene transfer event, and an unidentified 456 amino acid protein, which we hypothesize serves as a G-protein-coupled receptor. The H. akashiwo chloroplast genomes share little synteny with other algal chloroplast genomes sequenced to date.

CONCLUSION

The fosmid cloning technique eliminates chloroplast isolation, does not require chloroplast DNA purification, and reduces sequencing processing time. Application of this method has provided new insights into chloroplast genome architecture, gene content and evolution within the stramenopile cluster.

Oxidative stress causes telomere damage in Fanconi anaemia cells - a possible predisposition for malignant transformation

Fanconi anaemia (FA) is an autosomal recessive and X-linked disease characterized by severe genetic instability and increased incidence of cancer. One explanation for this instability may be the cellular hypersensitivity to oxidative stress leading to chromosomal breaks. This study explored the possible oxidative damage to telomeres of FA lymphocyte cell line, HSC536/N, and its possible effect on telomere function. We postulated that combination of oxidative damage with overexpression of telomerase may provide a possible model for malignant transformation in FA. The cells were grown in the presence of telomerase inhibitor and exposed for 1 month to H(2)O(2) combined with various antioxidants. This exposure caused shortening of telomere length and damage to the telomere single stranded overhang, which was prevented by several oxidants. This shortening was associated with development of severe telomere dysfunction. Control cells did not exhibit this sensitivity to H(2)O(2). Telomere dysfunction did not evoke damage response in FA cells, in contrast to normal P53 upregulation in control cells. Reconstitution of telomerase activity protected FA telomeres from further oxidative damage. These results suggest a scenario in which oxidative stress causes telomere shortening and ensuing telomere dysfunction may form the basis for malignant transformation in FA cells. Upregulation of telomerase activity in sporadic FA cells may perpetuate that process, thus explaining the malignant character of FA cells in vivo.

Strategies for reliable exploitation of evolutionary concepts in high throughput biology

The recent availability of the complete genome sequences of a large number of model organisms, together with the immense amount of data being produced by the new high-throughput technologies, means that we can now begin comparative analyses to understand the mechanisms involved in the evolution of the genome and their consequences in the study of biological systems. Phylogenetic approaches provide a unique conceptual framework for performing comparative analyses of all this data, for propagating information between different systems and for predicting or inferring new knowledge. As a result, phylogeny-based inference systems are now playing an increasingly important role in most areas of high throughput genomics, including studies of promoters (phylogenetic footprinting), interactomes (based on the presence and degree of conservation of interacting proteins), and in comparisons of transcriptomes or proteomes (phylogenetic proximity and co-regulation/co-expression). Here we review the recent developments aimed at making automatic, reliable phylogeny-based inference feasible in large-scale projects. We also discuss how evolutionary concepts and phylogeny-based inference strategies are now being exploited in order to understand the evolution and function of biological systems. Such advances will be fundamental for the success of the emerging disciplines of systems biology and synthetic biology, and will have wide-reaching effects in applied fields such as biotechnology, medicine and pharmacology.

Pairing and synapsis in wild type Arabidopsis thaliana

A spreading technique was used to perform a structural analysis of prophase I nuclei in pollen mother cells (PMCs) of wild-type Arabidopsis thaliana. In leptotene, all chromosomes developed fully axial elements before a presynaptic alignment was observed. Pairing and synapsis start in regions close to the telomeres at early zygotene. Interstitial synaptonemal complex (SC) stretches were found to occur at several sites per bivalent at mid zygotene. Within individual bivalents, extensive regions of SC formation often existed at the same time as other extensive regions that were unsynapsed. Also in the same nucleus, one bivalent might have several SC segments, while other bivalents have only a few. The classical bouquet was not so evident as in other plant species. Length measurements of the five pachytene bivalents have allowed the elaboration of a pachytene karyotype. Pachytene chromatin compaction in Arabidopsis was significantly less than that observed in the other species analysed and this is paralleled with a higher recombination rate (centimorgans per megabase).

A high through-put reverse genetic screen identifies two genes involved in remote memory in mice

Previous studies have revealed that the initial stages of memory formation require several genes involved in synaptic, transcriptional and translational mechanisms. In contrast, very little is known about the molecular and cellular mechanisms underlying later stages of memory, including remote memory (i.e. 7-day memory). To identify genes required for remote memory, we screened randomly selected mouse strains harboring known mutations. In our primary reverse genetic screen, we identified 4 putative remote memory mutant strains out of a total of 54 lines analyzed. Additionally, we found 11 other mutant strains with other abnormal profiles. Secondary screens confirmed that mutations of integrin β2 (Itgβ2) and steryl-O-acyl transferase 1 (Soat1) specifically disrupted remote memory. This study identifies some of the first genes required for remote memory, and suggests that screens of targeted mutants may be an efficient strategy to identify molecular requirements for this process.

Primate chromosome evolution: ancestral karyotypes, marker order and neocentromeres

In 1992 the Japanese macaque was the first species for which the homology of the entire karyotype was established by cross-species chromosome painting. Today, there are chromosome painting data on more than 50 species of primates. Although chromosome painting is a rapid and economical method for tracking translocations, it has limited utility for revealing intrachromosomal rearrangements. Fortunately, the use of BAC-FISH in the last few years has allowed remarkable progress in determining marker order along primate chromosomes and there are now marker order data on an array of primate species for a good number of chromosomes. These data reveal inversions, but also show that centromeres of many orthologous chromosomes are embedded in different genomic contexts. Even if the mechanisms of neocentromere formation and progression are just beginning to be understood, it is clear that these phenomena had a significant impact on shaping the primate genome and are fundamental to our understanding of genome evolution. In this report we complete and integrate the dataset of BAC-FISH marker order for human syntenies 1, 2, 4, 5, 8, 12, 17, 18, 19, 21, 22 and the X. These results allowed us to develop hypotheses about the content, marker order and centromere position in ancestral karyotypes at five major branching points on the primate evolutionary tree: ancestral primate, ancestral anthropoid, ancestral platyrrhine, ancestral catarrhine and ancestral hominoid. Current models suggest that between-species structural rearrangements are often intimately related to speciation. Comparative primate cytogenetics has become an important tool for elucidating the phylogeny and the taxonomy of primates. It has become increasingly apparent that molecular cytogenetic data in the future can be fruitfully combined with whole-genome assemblies to advance our understanding of primate genome evolution as well as the mechanisms and processes that have led to the origin of the human genome.

The surprisingly diverse ways that prokaryotes move

Prokaryotic cells move through liquids or over moist surfaces by swimming, swarming, gliding, twitching or floating. An impressive diversity of motility mechanisms has evolved in prokaryotes. Movement can involve surface appendages, such as flagella that spin, pili that pull and Mycoplasma 'legs' that walk. Internal structures, such as the cytoskeleton and gas vesicles, are involved in some types of motility, whereas the mechanisms of some other types of movement remain mysterious. Regardless of the type of motility machinery that is employed, most motile microorganisms use complex sensory systems to control their movements in response to stimuli, which allows them to migrate to optimal environments.

Transcriptome Analysis of ESTs from a Chaetognath Reveals a Deep-Branching Clade of Retrovirus-Like Retrotransposons

Chaetognaths constitute a small marine phylum exhibiting several characteristic which are highly unusual in animal genomes, including two classes of both rRNA and protein ribosomal genes. As in this phylum presence of retrovirus-like elements has never been documented, analysis of a published expressed sequence tag (EST) collection of the chaetognath Spadella cephaloptera has been made. Twelve sequences representing transcript sections of reverse transcriptase domain of active retrotransposons were isolated from~11,000 ESTs. Five of them are originated from Gypsy retrovirus-like elements, whereas the other are transcripts from a Bel-Pao LTR-retrotransposon, a Penelope-like element and LINE retrotransposons. Moreover, a part of a putative integrase has also been found. Phylogenetic analyses suggest a deep-branching clade of the retrovirus-like elements, which is in agreement with the probably Cambrian origin of the phylum. Moreover, retrotransposons have not been found in telomeric-like transcripts which are probably constituted by both vertebrate and arthropod canonical repeats.

The mosaic genome of Anaeromyxobacter dehalogenans strain 2CP-C suggests an aerobic common ancestor to the delta-proteobacteria

Anaeromyxobacter dehalogenans strain 2CP-C is a versaphilic delta-Proteobacterium distributed throughout many diverse soil and sediment environments. 16S rRNA gene phylogenetic analysis groups A. dehalogenans together with the myxobacteria, which have distinguishing characteristics including strictly aerobic metabolism, sporulation, fruiting body formation, and surface motility. Analysis of the 5.01 Mb strain 2CP-C genome substantiated that this organism is a myxobacterium but shares genotypic traits with the anaerobic majority of the delta-Proteobacteria (i.e., the Desulfuromonadales). Reflective of its respiratory versatility, strain 2CP-C possesses 68 genes coding for putative c-type cytochromes, including one gene with 40 heme binding motifs. Consistent with its relatedness to the myxobacteria, surface motility was observed in strain 2CP-C and multiple types of motility genes are present, including 28 genes for gliding, adventurous (A-) motility and 17 genes for type IV pilus-based motility (i.e., social (S-) motility) that all have homologs in Myxococcus xanthus. Although A. dehalogenans shares many metabolic traits with the anaerobic majority of the delta-Proteobacteria, strain 2CP-C grows under microaerophilic conditions and possesses detoxification systems for reactive oxygen species. Accordingly, two gene clusters coding for NADH dehydrogenase subunits and two cytochrome oxidase gene clusters in strain 2CP-C are similar to those in M. xanthus. Remarkably, strain 2CP-C possesses a third NADH dehydrogenase gene cluster and a cytochrome cbb₃ oxidase gene cluster, apparently acquired through ancient horizontal gene transfer from a strictly anaerobic green sulfur bacterium. The mosaic nature of the A. dehalogenans strain 2CP-C genome suggests that the metabolically versatile, anaerobic members of the delta-Proteobacteria may have descended from aerobic ancestors with complex lifestyles.

Origin of INSL3-mediated testicular descent in therian mammals

Testicular descent is a unique physiological adaptation found in therian mammals allowing optimal spermatogenesis below core body temperature. Recent studies show that INSL3, produced by Leydig cells, and its receptor LGR8 (RXFP2) are essential for mediating the transabdominal phase of testicular descent during early development. However, the origin and genetic basis for this physiological adaptation is not clear. Using syntenic mapping and the functional characterization of contemporary and resurrected relaxin family hormones, we show that derivation of INSL3-mediated testicular descent involved the duplication of an ancestral RLN3-like gene that encodes an indiscriminate ligand for LGR7 (RXFP1) and LGR8. This event was followed by acquisition of the LGR7-selective characteristics by a daughter gene (RLN3) prior to the evolution of the common ancestor of monotremes, marsupials, and placentals. A subsequent mutation of the other daughter gene (INSL3) occurred before the emergence of therian mammals, which then led to the derivation of the reciprocal LGR8-specific characteristics of INSL3. The stepwise evolution of these independent signaling pathways through gene duplication and subsequent divergence is consistent with Darwinian theory of selection and adaptation, and the temporal proximity suggests an association between these genetic events and the concurrent evolution of testicular descent in ancestral therian mammals.

SalK/SalR, a two-component signal transduction system, is essential for full virulence of highly invasive Streptococcus suis serotype 2

Background

Streptococcus suis serotype 2 (S. suis 2, SS2) has evolved into a highly infectious entity, which caused the two recent large-scale outbreaks of human SS2 epidemic in China, and is characterized by a toxic shock-like syndrome. However, the molecular pathogenesis of this new emerging pathogen is still poorly understood.

Methodology/Principal Findings

89K is a newly predicted pathogenicity island (PAI) which is specific to Chinese epidemic strains isolated from these two SS2 outbreaks. Further bioinformatics analysis revealed a unique two-component signal transduction system (TCSTS) located in the candidate 89K PAI, which is orthologous to the SalK/SalR regulatory system of Streptococcus salivarius. Knockout of salKR eliminated the lethality of SS2 in experimental infection of piglets. Functional complementation of salKR into the isogenic mutant ΔsalKR restored its soaring pathogenicity. Colonization experiments showed that the ΔsalKR mutant could not colonize any susceptible tissue of piglets when administered alone. Bactericidal assays demonstrated that resistance of the mutant to polymorphonuclear leukocyte (PMN)-mediated killing was greatly decreased. Expression microarray analysis exhibited a transcription profile alteration of 26 various genes down-regulated in the ΔsalKR mutant.

Conclusions/Significance

These findings suggest that SalK/SalR is requisite for the full virulence of ethnic Chinese isolates of highly pathogenic SS2, thus providing experimental evidence for the validity of this bioinformatically predicted PAI.

The function of the clock-associated transcriptional regulator CCA1 (CIRCADIAN CLOCK-ASSOCIATED 1) in Arabidopsis thaliana

In the model higher plant Arabidopsis thaliana, the CCA1 (CIRCADIAN CLOCK-ASSOCIATED 1) gene plays important circadian clock-associated roles. The CCA1 protein is a member of a small subfamily of single MYB-related transcription factors. This family consists of several homologous CCA1-like transcription factors, including the closest homolog LHY (LATE ELONGATED HYPOCOTYL). To gain insight into the molecular function of CCA1 and its homologs, here we took a unique genetic approach that was recently developed for Arabidopsis thaliana. Through this strategy, referred to as CRES-T (Chimeric REpressor Silencing Technology), a transgenic plant was constructed to produce a dominant negative transcriptional repressor (designated CCA1-SRDX). By employing the resulting transgenic lines, together with previously established cca1 lhy double mutant and CCA1-ox (over-expressing) plants, their circadian clock-associated phenotypes were examined and compared with each other. The observed clock-associated phenotypes of the CCA1-SRDX plants were very similar to those of CCA1-ox, but not to those of cca1 lhy, suggesting that CCA1 acts predominantly as a transcriptional repressor in nature. However, the developmental morphology (or architecture) of adult CCA1-SRDX plants were quite different from that of CCA1-ox, suggesting that CCA1 might also be implicated, directly or indirectly, in an as yet unknown circadian-associated output pathway at a late developmental stage.

Recombinant adeno-associated viral vectors are deficient in provoking a DNA damage response

Adeno-associated virus type 2 (AAV2) provokes a DNA damage response that mimics a stalled replication fork. We have previously shown that this response is dependent on ataxia telangiectasia-mutated and Rad3-related kinase and involves recruitment of DNA repair proteins into foci associated with AAV2 DNA. Here, we investigated whether recombinant AAV2 (rAAV2) vectors are able to produce a similar response. Surprisingly, the results show that both single-stranded and double-stranded green fluorescent protein-expressing rAAV2 vectors are defective in producing such a response. We show that the DNA damage signaling initiated by AAV2 was not due to the virus-encoded Rep or viral capsid proteins. UV-inactivated AAV2 induced a response similar to that of untreated AAV2. This type of DNA damage response was not provoked by other DNA molecules, such as single-stranded bacteriophage M13 or plasmid DNAs. Rather, the results indicate that the ability of AAV2 to produce a DNA damage response can be attributed to the presence of cis-acting AAV2 DNA sequences, which are absent in rAAV2 vectors and could function as origins of replication creating stalled replication complexes. This hypothesis was tested by using a single-stranded rAAV2 vector containing the p5 AAV2 sequence that has previously been shown to enhance AAV2 replication. This vector was indeed able to trigger DNA damage signaling. These findings support the conclusion that efficient formation of AAV2 replication complexes is required for this AAV2-induced DNA damage response and provide an explanation for the poor response in rAAV2-infected cells.

Transcription induces strand-specific mutations at the 5' end of human genes

A regional analysis of nucleotide substitution rates along human genes and their flanking regions allows us to quantify the effect of mutational mechanisms associated with transcription in germ line cells. Our analysis reveals three distinct patterns of substitution rates. First, a sharp decline in the deamination rate of methylated CpG dinucleotides, which is observed in the vicinity of the 5' end of genes. Second, a strand asymmetry in complementary substitution rates, which extends from the 5' end to 1 kbp downstream from the 3' end, associated with transcription-coupled repair. Finally, a localized strand asymmetry, an excess of C-->T over G-->A substitution in the nontemplate strand confined to the first 1-2 kbp downstream of the 5' end of genes. We hypothesize that higher exposure of the nontemplate strand near the 5' end of genes leads to a higher cytosine deamination rate. Up to now, only the somatic hypermutation (SHM) pathway has been known to mediate localized and strand-specific mutagenic processes associated with transcription in mammalia. The mutational patterns in SHM are induced by cytosine deaminase, which just targets single-stranded DNA. This DNA conformation is induced by R-loops, which preferentially occur at the 5' ends of genes. We predict that R-loops are extensively formed in the beginning of transcribed regions in germ line cells.

Alternative splicing: a missing piece in the puzzle of intron gain

Spliceosomal introns, a hallmark of eukaryotic gene organization, were an unexpected discovery. After three decades, crucial issues such as when and how introns first appeared in evolution remain unsettled. An issue yet to be answered is how intron positions arise de novo. Phylogenetic investigations concur that intron positions continue to emerge, at least in some lineages. Yet genomic scans for the sources of introns occupying new positions have been fruitless. Two alternative solutions to this paradox are: (i) formation of new intron positions halted before the recent past and (ii) it continues to occur, but through processes different from those generally assumed. One process generally dismissed is intron sliding--the relocation of a preexisting intron over short distances--because of supposed associated deleterious effects. The puzzle of intron gain arises owing to a pervasive operational definition of introns, which sees them as precisely demarcated segments of the genome separated from the neighboring nonintronic DNA by unmovable limits. Intron homology is defined as position homology. Recent studies of pre-mRNA processing indicate that this assumption needs to be revised. We incorporate recent advances on the evolutionarily frequent process of alternative splicing, by which exons of primary transcripts are spliced in different patterns, into a new model of intron sliding that accounts for the diversity of intron positions. We posit that intron positional diversity is driven by two overlapping processes: (i) background process of continuous relocation of preexisting introns by sliding and (ii) spurts of extensive gain/loss of new intron sequences.

Understanding and re-engineering nucleoprotein machines to cure human disease

The mammalian nucleus is filled with self-organizing, nanometer-scale nucleoprotein machines that carry out DNA replication, RNA biogenesis and DNA repair. We discuss, as a model, the nonhomologous end-joining (NHEJ) machine, which repairs DNA double-strand breaks. The NHEJ machine consists of six core polypeptides and 10-20 ancillary polypeptides. A full understanding of its design principles will require measuring the behavior of single NHEJ complexes in living cells, using a Nano Toolbox that includes bright, stable, biocompatible fluorophores, efficient protein and nucleic acid-tagging strategies, and sensitive, high-resolution imaging methods. Taking inspiration from natural examples, it might be possible to adapt and redesign the NHEJ machine to precisely correct mutations responsible for common human diseases, such as K-ras in lung cancer or human papillomavirus E6 and E7 genes in cervical and oral cancers.

Transcriptional regulatory network for sexual differentiation in fission yeast

Microarray analysis of the transcriptome of fission yeast after genetic perturbation of 6 genes known to have a role in sexual differentiation reveals insights into the regulatory principles controlling the gene expression program driving this process.

Background

Changes in gene expression are hallmarks of cellular differentiation. Sexual differentiation in fission yeast (Schizosaccharomyces pombe) provides a model system for gene expression programs accompanying and driving cellular specialization. The expression of hundreds of genes is modulated in successive waves during meiosis and sporulation in S. pombe, and several known transcription factors are critical for these processes.

Results

We used DNA microarrays to investigate meiotic gene regulation by examining transcriptomes after genetic perturbations (gene deletion and/or overexpression) of rep1, mei4, atf21 and atf31, which encode known transcription factors controlling sexual differentiation. This analysis reveals target genes at a genome-wide scale and uncovers combinatorial control by Atf21p and Atf31p. We also studied two transcription factors not previously implicated in sexual differentiation whose meiotic induction depended on Mei4p: Rsv2p induces stress-related genes during spore formation, while Rsv1p represses glucose-metabolism genes. Our data further reveal negative feedback interactions: both Rep1p and Mei4p not only activate specific gene expression waves (early and middle genes, respectively) but are also required for repression of genes induced in the previous waves (Ste11p-dependent and early genes, respectively).

Conclusion

These data give insight into regulatory principles controlling the extensive gene expression program driving sexual differentiation and highlight sophisticated interactions and combinatorial control among transcription factors. Besides triggering simultaneous expression of gene waves, transcription factors also repress genes in the previous wave and induce other factors that in turn regulate a subsequent wave. These dependencies ensure an ordered and timely succession of transcriptional waves during cellular differentiation.

Genetic subtraction profiling identifies genes essential for Arabidopsis reproduction and reveals interaction between the female gametophyte and the maternal sporophyte

Genetic subtraction and expression profiling of wild-type Arabidopsis and a sporophytic mutant lacking an embryo sac identified 1,260 genes expressed in the embryo sac; a total of 527 genes were identified for their expression in ovules of mutants lacking an embryo sac.

Background

The embryo sac contains the haploid maternal cell types necessary for double fertilization and subsequent seed development in plants. Large-scale identification of genes expressed in the embryo sac remains cumbersome because of its inherent microscopic and inaccessible nature. We used genetic subtraction and comparative profiling by microarray between the Arabidopsis thaliana wild-type and a sporophytic mutant lacking an embryo sac in order to identify embryo sac expressed genes in this model organism. The influences of the embryo sac on the surrounding sporophytic tissues were previously thought to be negligible or nonexistent; we investigated the extent of these interactions by transcriptome analysis.

Results

We identified 1,260 genes as embryo sac expressed by analyzing both our dataset and a recently reported dataset, obtained by a similar approach, using three statistical procedures. Spatial expression of nine genes (for instance a central cell expressed trithorax-like gene, an egg cell expressed gene encoding a kinase, and a synergid expressed gene encoding a permease) validated our approach. We analyzed mutants in five of the newly identified genes that exhibited developmental anomalies during reproductive development. A total of 527 genes were identified for their expression in ovules of mutants lacking an embryo sac, at levels that were twofold higher than in the wild type.

Conclusion

Identification of embryo sac expressed genes establishes a basis for the functional dissection of embryo sac development and function. Sporophytic gain of expression in mutants lacking an embryo sac suggests that a substantial portion of the sporophytic transcriptome involved in carpel and ovule development is, unexpectedly, under the indirect influence of the embryo sac.