Arbuscular mycorrhizal fungi (Glomeromycota) harbour ancient fungal tubulin genes that resemble those of the chytrids (Chytridiomycota).

The genes encoding alpha- and beta-tubulins have been widely sampled in most major fungal phyla and they are useful tools for fungal phylogeny. Here, we report the first isolation of alpha-tubulin sequences from arbuscular mycorrhizal fungi (AMF). In parallel, AMF beta-tubulins were sampled and analysed to identify the presence of paralogs of this gene. The AMF alpha-tubulin amino acid phylogeny was congruent with the results previously reported for AMF beta-tubulins and showed that AMF tubulins group together at a basal position in the fungal clade and showed high sequence similarities with members of the Chytridiomycota. This is in contrast with phylogenies for other regions of the AMF genome. The amount and nature of substitutions are consistent with an ancient divergence of both orthologs and paralogs of AMF tubulins. At the amino acid level, however, AMF tubulins have hardly evolved from those of the chytrids. This is remarkable given that these two groups are ancient and the monophyletic Glomeromycota probably diverged from basal fungal ancestors at least 500 million years ago. The specific primers we designed for the AMF tubulins, together with the high molecular variation we found among the AMF species we analysed, make AMF tubulin sequences potentially useful for AMF identification purposes.

Additive functions in boolean models of gene regulatory network modules.

Gene-on-gene regulations are key components of every living organism. Dynamical abstract models of genetic regulatory networks help explain the genome's evolvability and robustness. These properties can be attributed to the structural topology of the graph formed by genes, as vertices, and regulatory interactions, as edges. Moreover, the actual gene interaction of each gene is believed to play a key role in the stability of the structure. With advances in biology, some effort was deployed to develop update functions in Boolean models that include recent knowledge. We combine real-life gene interaction networks with novel update functions in a Boolean model. We use two sub-networks of biological organisms, the yeast cell-cycle and the mouse embryonic stem cell, as topological support for our system. On these structures, we substitute the original random update functions by a novel threshold-based dynamic function in which the promoting and repressing effect of each interaction is considered. We use a third real-life regulatory network, along with its inferred Boolean update functions to validate the proposed update function. Results of this validation hint to increased biological plausibility of the threshold-based function. To investigate the dynamical behavior of this new model, we visualized the phase transition between order and chaos into the critical regime using Derrida plots. We complement the qualitative nature of Derrida plots with an alternative measure, the criticality distance, that also allows to discriminate between regimes in a quantitative way. Simulation on both real-life genetic regulatory networks show that there exists a set of parameters that allows the systems to operate in the critical region. This new model includes experimentally derived biological information and recent discoveries, which makes it potentially useful to guide experimental research. The update function confers additional realism to the model, while reducing the complexity and solution space, thus making it easier to investigate.

Mammalian genes are transcribed with widely different bursting kinetics.

In prokaryotes and eukaryotes, most genes appear to be transcribed during short periods called transcriptional bursts, interspersed by silent intervals. We describe how such bursts generate gene-specific temporal patterns of messenger RNA (mRNA) synthesis in mammalian cells. To monitor transcription at high temporal resolution, we established various gene trap cell lines and transgenic cell lines expressing a short-lived luciferase protein from an unstable mRNA, and recorded bioluminescence in real time in single cells. Mathematical modeling identified gene-specific on- and off-switching rates in transcriptional activity and mean numbers of mRNAs produced during the bursts. Transcriptional kinetics were markedly altered by cis-regulatory DNA elements. Our analysis demonstrated that bursting kinetics are highly gene-specific, reflecting refractory periods during which genes stay inactive for a certain time before switching on again.

Whole-genome sequences of influenza A(H3N2) viruses isolated from Brazilian patients with mild illness during the 2014 season

The influenza A(H3N2) virus has circulated worldwide for almost five decades and is the dominant subtype in most seasonal influenza epidemics, as occurred in the 2014 season in South America. In this study we evaluate five whole genome sequences of influenza A(H3N2) viruses detected in patients with mild illness collected from January-March 2014. To sequence the genomes, a new generation sequencing (NGS) protocol was performed using the Ion Torrent PGM platform. In addition to analysing the common genes, haemagglutinin, neuraminidase and matrix, our work also comprised internal genes. This was the first report of a whole genome analysis with Brazilian influenza A(H3N2) samples. Considerable amino acid variability was encountered in all gene segments, demonstrating the importance of studying the internal genes. NGS of whole genomes in this study will facilitate deeper virus characterisation, contributing to the improvement of influenza strain surveillance in Brazil.

Draft genome sequences of three NDM-1-producing Enterobacteriaceae species isolated from Brazil

The emergence of multidrug-resistant Enterobacteriaceae strains producing carbapenemases, such as NDM-1, has become a major public health issue due to a high dissemination capacity and limited treatment options. Here we describe the draft genome of three NDM-1-producing isolates: Providencia rettgeri(CCBH11880), Enterobacter hormaecheisubsp. oharae(CCBH10892) and Klebsiella pneumoniae(CCBH13327), isolated in Brazil. BesidesblaNDM-1, resistance genes to aminoglycosides [aadA1, aadA2,aac(6’)-Ib-cr] and quinolones (qnrA1,qnrB4) were observed which contributed to the multidrug resistance profile. The element ISAba125 was found associated to theblaNDM-1 gene in all strains.

CATMA, a comprehensive genome-scale resource for silencing and transcript profiling of Arabidopsis genes.

BACKGROUND: The Complete Arabidopsis Transcript MicroArray (CATMA) initiative combines the efforts of laboratories in eight European countries 1 to deliver gene-specific sequence tags (GSTs) for the Arabidopsis research community. The CATMA initiative offers the power and flexibility to regularly update the GST collection according to evolving knowledge about the gene repertoire. These GST amplicons can easily be reamplified and shared, subsets can be picked at will to print dedicated arrays, and the GSTs can be cloned and used for other functional studies. This ongoing initiative has already produced approximately 24,000 GSTs that have been made publicly available for spotted microarray printing and RNA interference. RESULTS: GSTs from the CATMA version 2 repertoire (CATMAv2, created in 2002) were mapped onto the gene models from two independent Arabidopsis nuclear genome annotation efforts, TIGR5 and PSB-EuGène, to consolidate a list of genes that were targeted by previously designed CATMA tags. A total of 9,027 gene models were not tagged by any amplified CATMAv2 GST, and 2,533 amplified GSTs were no longer predicted to tag an updated gene model. To validate the efficacy of GST mapping criteria and design rules, the predicted and experimentally observed hybridization characteristics associated to GST features were correlated in transcript profiling datasets obtained with the CATMAv2 microarray, confirming the reliability of this platform. To complete the CATMA repertoire, all 9,027 gene models for which no GST had yet been designed were processed with an adjusted version of the Specific Primer and Amplicon Design Software (SPADS). A total of 5,756 novel GSTs were designed and amplified by PCR from genomic DNA. Together with the pre-existing GST collection, this new addition constitutes the CATMAv3 repertoire. It comprises 30,343 unique amplified sequences that tag 24,202 and 23,009 protein-encoding nuclear gene models in the TAIR6 and EuGène genome annotations, respectively. To cover the remaining untagged genes, we identified 543 additional GSTs using less stringent design criteria and designed 990 sequence tags matching multiple members of gene families (Gene Family Tags or GFTs) to cover any remaining untagged genes. These latter 1,533 features constitute the CATMAv4 addition. CONCLUSION: To update the CATMA GST repertoire, we designed 7,289 additional sequence tags, bringing the total number of tagged TAIR6-annotated Arabidopsis nuclear protein-coding genes to 26,173. This resource is used both for the production of spotted microarrays and the large-scale cloning of hairpin RNA silencing vectors. All information about the resulting updated CATMA repertoire is available through the CATMA database http://www.catma.org.

HCF-1 self-association via an interdigitated Fn3 structure facilitates transcriptional regulatory complex formation.

Host-cell factor 1 (HCF-1) is an unusual transcriptional regulator that undergoes a process of proteolytic maturation to generate N- (HCF-1(N)) and C- (HCF-1(C)) terminal subunits noncovalently associated via self-association sequence elements. Here, we present the crystal structure of the self-association sequence 1 (SAS1) including the adjacent C-terminal HCF-1 nuclear localization signal (NLS). SAS1 elements from each of the HCF-1(N) and HCF-1(C) subunits form an interdigitated fibronectin type 3 (Fn3) tandem repeat structure. We show that the C-terminal NLS recruited by the interdigitated SAS1 structure is required for effective formation of a transcriptional regulatory complex: the herpes simplex virus VP16-induced complex. Thus, HCF-1(N)-HCF-1(C) association via an integrated Fn3 structure permits an NLS to facilitate formation of a transcriptional regulatory complex.

Analysis of genes encoding penicillin-binding proteins in clinical isolates of Acinetobacter baumannii.

There is limited information on the role of penicillin-binding proteins (PBPs) in the resistance of Acinetobacter baumannii to β-lactams. This study presents an analysis of the allelic variations of PBP genes in A. baumannii isolates. Twenty-six A. baumannii clinical isolates (susceptible or resistant to carbapenems) from three teaching hospitals in Spain were included. The antimicrobial susceptibility profile, clonal pattern, and genomic species identification were also evaluated. Based on the six complete genomes of A. baumannii, the PBP genes were identified, and primers were designed for each gene. The nucleotide sequences of the genes identified that encode PBPs and the corresponding amino acid sequences were compared with those of ATCC 17978. Seven PBP genes and one monofunctional transglycosylase (MGT) gene were identified in the six genomes, encoding (i) four high-molecular-mass proteins (two of class A, PBP1a [ponA] and PBP1b [mrcB], and two of class B, PBP2 [pbpA or mrdA] and PBP3 [ftsI]), (ii) three low-molecular-mass proteins (two of type 5, PBP5/6 [dacC] and PBP6b [dacD], and one of type 7 (PBP7/8 [pbpG]), and (iii) a monofunctional enzyme (MtgA [mtgA]). Hot spot mutation regions were observed, although most of the allelic changes found translated into silent mutations. The amino acid consensus sequences corresponding to the PBP genes in the genomes and the clinical isolates were highly conserved. The changes found in amino acid sequences were associated with concrete clonal patterns but were not directly related to susceptibility or resistance to β-lactams. An insertion sequence disrupting the gene encoding PBP6b was identified in an endemic carbapenem-resistant clone in one of the participant hospitals.

Evolution and expression of vitellogenin genes.

A major event during the growth period of oocytes in nonmammalian animals is the accumulation of yolk. The genes coding for the yolk protein precursor, known as vitallogenin, are well characterized in a few vertebrate and invertebrate species. Studies on the evolution of these genes and on the regulatory mechanisms involved in their time, tissue- and hormone-specific expression are presented and discussed in this review.

Integrative analysis of RUNX1 downstream pathways and target genes.

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications

Extracellular protease and phospholipase C are controlled by the global regulatory gene gacA in the biocontrol strain Pseudomonas fluorescens CHA0.

Pseudomonas fluorescens strain CHA0 protects plants from various root diseases. Antibiotic metabolites synthesized by this strain play an important role in disease suppression; their production is mediated by the global activator gene gacA. Here we show by complementation that the gacA gene is also essential for the expression of two extracellular enzymes in P. fluorescens CHA0: phospholipase C and a 47-kDa metalloprotease. In contrast, the production of another exoenzyme, lipase, is not regulated by the gacA gene. Protease, phospholipase and antibiotics of P. fluorescens are all known to be optimally produced at the end of exponential growth; thus, the gacA gene appears to be a general stationary-phase regulator.

Control of pancreatic β-cell mass and function by gluco-incretin hormones : identification of novel regulatory mechanisms for the treatment of diabetes

Summary : Control of pancreatic ß-cell mass and function by gluco-incretin hormones: Identification of novel regulatory mechanisms for the treatment of diabetes The ß-cells of islets of Langerhans secrete insulin to reduce hyperglycemia. The number of pancreatic islet ß-cells and their capacity to secrete insulin is modulated in normal physiological conditions to respond to the metabolic demand of the organism. A failure of the endocrine pancreas to maintain an adequate insulin secretory capacity due to a reduced ß-cell number and function underlies the pathogenesis of both type 1 and type 2 diabetes. The molecular mechanisms controlling the glucose competence of mature ß-cells, i.e., the magnitude of their insulin secretion response to glucose, ß-cell replication, their differentiation from precursor cells and protection against apoptosis are poorly understood. To investigate these mechanisms, we studied the effects on ß-cells of the gluco-incretin hormones, glucose-dependent insulinotropic polypeptide (GIP) and glucagon-like peptide-1 (GLP-1) which are secreted by intestinal endocrine cells after food intake. Besides acutely potentiating glucose-stimulated insulin secretion, these hormones induce ß-cell differentiation from precursor cells, stimulate mature ß-cell replication, and protect them against apoptosis. Therefore, understanding the molecular basis for gluco-incretin action may lead to the uncovering of novel ß-cell regulatory events with potential application for the treatment or prevention of diabetes. Islets from mice with inactivation of both GIP and GLP-1 receptor genes (dK0) present a defect in glucose-induced insulin secretion and are more sensitive than control islets to cytokine-induced apoptosis. To search for regulatory genes, that may control both glucose competence and protection against apoptosis, we performed comparative transcriptomic analysis of islets from control and dK0 mice. We found a strong down-regulation of the IGF1 Rexpression in dK0 islets. We demonstrated in both a mouse insulin-secreting cell line and primary islets, that GLP-1 stimulated IGF-1R expression and signaling. Importantly, GLP-1induced IGF-1R-dependent Akt phosphorylation required active secretion, indicating the presence of an autocrine activation mechanism. We further showed that activation of IGF-1R signaling was dependent on the secretion of IGF-2 and IGF-2 expression was regulated by nutrients. Finally, we demonstrated that the IGF-Z/IGF-1R autocrine loop was required for GLP-1 i) to protect ß-cells against cytokine-induced apoptosis, ii) to enhance their glucose competence and iii) to increase ß-cell proliferation. Résumé : Contrôle de la masse des cellules ß pancréatiques et de leur fonction par les hormones glucoincrétines: Identification de nouveaux mécanismes régulateurs pour le traitement du diabète Les cellules ß des îlots de Langerhans sécrètent l'insuline pour diminuer l'hyperglycémie. Le nombre de cellules ß et leur capacité à sécréter l'insuline sont modulés dans les conditions physiologiques normales pour répondre à la demande métabolique de l'organisme. Un échec du pancréas endocrine à maintenir sa capacité sécrétoire d'insuline dû à une diminution du nombre et de la fonction des cellules ß conduit au diabète de type 1 et de type 2. Les mécanismes moléculaires contrôlant la compétence au glucose des cellules ß matures, tels que, l'augmentation de la sécrétion d'insuline en réponse au glucose, la réplication des cellules ß, leur différentiation à partir de cellules précurseurs et la protection contre l'apoptose sont encore peu connus. Afin d'examiner ces mécanismes, nous avons étudié les effets sur les cellules ß des hormones gluco-incrétines, glucose-dépendent insulinotropic polypeptide (G1P) et glucagon-like peptide-1 (GLP-1) qui sont sécrétées par les cellules endocrines de l'intestin après la prise alimentaire. En plus de potentialiser la sécrétion d'insuline induite par le glucose, ces hormones induisent la différentiation de cellules ß à partir de cellules précurseurs, stimulent leur prolifération et les protègent contre l'apoptose. Par conséquent, comprendre les mécanismes d'action des gluco-incrétines permettrait de découvrir de nouveaux processus régulant les cellules ß avec d'éventuelles applications dans le traitement ou la prévention du diabète. Les îlots de souris ayant une double inactivation des gènes pour les récepteurs du GIP et du GLP-1 (dK0) présentent un défaut de sécrétion d'insuline stimulée par le glucose et une sensibilité accrue à l'apoptose induite par les cytokines. Afin de déterminer les gènes régulés, qui pourraient contrôler à la fois la compétence au glucose et la protection contre l'apoptose, nous avons effectué une analyse comparative transcriptomique sur des îlots de souris contrôles et dKO. Nous avons constaté une forte diminution de l'expression d'IGF-1R dans les îlots dKO. Nous avons démontré, à la fois dans une lignée cellulaire murine sécrétant l'insuline et dans îlots primaires, que le GLP-1 stimulait l'expression d'IGF-1R et sa voie de signalisation. Par ailleurs, la phosphorylation d'Akt dépendante d'IGF1-R induite parle GLP-1 nécessite une sécrétion active, indiquant la présence d'un mécanisme d'activation autocrine. Nous avons ensuite montré que l'activation de la voie de signalisation d'IGF-1R était dépendante de la sécrétion d'IGF-2, dont l'expression est régulée par les nutriments. Finalement, nous avons démontré que la boucle autocrine IGF-2/IGF-1R est nécessaire pour le GLP-1 i) pour protéger les cellules ß contre l'apoptose induite par les cytokines, ii) pour améliorer la compétence au glucose et iii) pour augmenter la prolifération des cellules ß. Résumé tout public : Contrôle de la masse des cellules ß pancréatiques et de leur fonction par les hormones gluco-incrétines: Identification de nouveaux mécanismes régulateurs pour le traitement du diabète Chez les mammifères, la concentration de glucose sanguine (glycémie) est régulée et maintenue à une valeur relativement constante d'environ 5 mM. Cette régulation est principalement contrôlée par 2 hormones produites par les îlots pancréatiques de Langerhans: l'insuline sécrétée par les cellules ß et le glucagon sécrété par les cellules a. A la suite d'un repas, l'augmentation de la glycémie entraîne la sécrétion d'insuline ce qui permet le stockage du glucose dans le foie, les muscles et le tissu adipeux afin de diminuer le taux de glucose circulant. Lors d'un jeûne, la diminution de la glycémie permet la sécrétion de glucagon favorisant alors la production de glucose par le foie, normalisant ainsi la glycémie. Le nombre de cellules ß et leur capacité sécrétoire s'adaptent aux variations de la demande métabolique pour assurer une normoglycémie. Une destruction complète ou partielle des cellules ß conduit respectivement au diabète de type 1 et de type 2. Bien que l'augmentation de la glycémie soit le facteur stimulant de la sécrétion d'insuline, des hormones gluco-incrétines, principalement le GLP-1 (glucagon-like peptide-1) et le GIP (glucose-dependent insulinotropic polypeptide) sont libérées par l'intestin en réponse aux nutriments (glucose, acides gras) et agissent au niveau des cellules ß, potentialisant la sécrétion d'insuline induite par le glucose, stimulant leur prolifération, induisant la différentiation de cellules précurseurs en cellules ß matures et les protègent contre la mort cellulaire (apoptose). Afin d'étudier plus en détail ces mécanismes, nous avons généré des souris déficientes pour les récepteurs du GIP et du GLP-l. Les îlots pancréatiques de ces souris présentent un défaut de sécrétion d'insuline stimulée par le glucose et une sensibilité accrue à l'apoptose par rapport aux îlots de souris contrôles. Nous avons donc cherché les gènes régulés pas ces hormones contrôlant la sécrétion d'insuline et la protection contre l'apoptose. Nous avons constaté une forte diminution de l'expression du récepteur à l'IGF-1 (IGF-1R) dans les îlots de souris déficientes pour les récepteurs des gluco-incrétines. Nous avons démontré dans un model de cellules ß en culture et d'îlots que le GLP-1 augmentait l'expression d'IGF-1R et la sécrétion de son ligand (IGF-2) permettant l'activation de la voie de signalisation. Finalement, nous avons montré que l'activation de la boucle IGF-2/IGF-1R induite par le GLP-1 était nécessaire pour la protection contre l'apoptose, l'augmentation de la sécrétion et la prolifération des cellules ß.

Functional diversification of duplicate genes through subcellular adaptation of encoded proteins.

BACKGROUND: Gene duplication is the primary source of new genes with novel or altered functions. It is known that duplicates may obtain these new functional roles by evolving divergent expression patterns and/or protein functions after the duplication event. Here, using yeast (Saccharomyces cerevisiae) as a model organism, we investigate a previously little considered mode for the functional diversification of duplicate genes: subcellular adaptation of encoded proteins. RESULTS: We show that for 24-37% of duplicate gene pairs derived from the S. cerevisiae whole-genome duplication event, the two members of the pair encode proteins that localize to distinct subcellular compartments. The propensity of yeast duplicate genes to evolve new localization patterns depends to a large extent on the biological function of their progenitor genes. Proteins involved in processes with a wider subcellular distribution (for example, catabolism) frequently evolved new protein localization patterns after duplication, whereas duplicate proteins limited to a smaller number of organelles (for example, highly expressed biosynthesis/housekeeping proteins with a slow rate of evolution) rarely relocate within the cell. Paralogous proteins evolved divergent localization patterns by partitioning of ancestral localizations ('sublocalization'), but probably more frequently by relocalization to new compartments ('neolocalization'). We show that such subcellular reprogramming may occur through selectively driven substitutions in protein targeting sequences. Notably, our data also reveal that relocated proteins functionally adapted to their new subcellular environments and evolved new functional roles through changes of their physico-chemical properties, expression levels, and interaction partners. CONCLUSION: We conclude that protein subcellular adaptation represents a common mechanism for the functional diversification of duplicate genes.

Organization of lin genes and IS6100 among different strains of hexachlorocyclohexane-degrading Sphingomonas paucimobilis: evidence for horizontal gene transfer.

The organization of lin genes and IS6100 was studied in three strains of Sphingomonas paucimobilis (B90A, Sp+, and UT26) which degraded hexachlorocyclohexane (HCH) isomers but which had been isolated at different geographical locations. DNA-DNA hybridization data revealed that most of the lin genes in these strains were associated with IS6100, an insertion sequence classified in the IS6 family and initially found in Mycobacterium fortuitum. Eleven, six, and five copies of IS6100 were detected in B90A, Sp+, and UT26, respectively. IS6100 elements in B90A were sequenced from five, one, and one regions of the genomes of B90A, Sp+, and UT26, respectively, and were found to be identical. DNA-DNA hybridization and DNA sequencing of cosmid clones also revealed that S. paucimobilis B90A contains three and two copies of linX and linA, respectively, compared to only one copy of these genes in strains Sp+ and UT26. Although the copy number and the sequence of the remaining genes of the HCH degradative pathway (linB, linC, linD, and linE) were nearly the same in all strains, there were striking differences in the organization of the linA genes as a result of replacement of portions of DNA sequences by IS6100, which gave them a strange mosaic configuration. Spontaneous deletion of linD and linE from B90A and of linA from Sp+ occurred and was associated either with deletion of a copy of IS6100 or changes in IS6100 profiles. The evidence gathered in this study, coupled with the observation that the G+C contents of the linA genes are lower than that of the remaining DNA sequence of S. paucimobilis, strongly suggests that all these strains acquired the linA gene through horizontal gene transfer mediated by IS6100. The association of IS6100 with the rest of the lin genes further suggests that IS6100 played a role in shaping the current lin gene organization.

Genometric analyses of the organization of circular chromosomes: a universal pressure determines the direction of ribosomal RNA genes transcription relative to chromosome replication.

Selective pressures related to gene function and chromosomal architecture are acting on genome sequences and can be revealed, for instance, by appropriate genometric methods. Cumulative nucleotide skew analyses, i.e., GC, TA, and ORF orientation skews, predict the location of the origin of DNA replication for 88 out of 100 completely sequenced bacterial chromosomes. These methods appear fully reliable for proteobacteria, Gram-positives, and spirochetes as well as for euryarchaeotes. Based on this genome architecture information, coorientation analyses reveal that in prokaryotes, ribosomal RNA (rRNA) genes encoding the small and large ribosomal subunits are all transcribed in the same direction as DNA replication; that is, they are located along the leading strand. This result offers a simple and reliable method for circumscribing the region containing the origin of the DNA replication and reveals a strong selective pressure acting on the orientation of rRNA genes similar to the weaker one acting on the orientation of ORFs. Rate of coorientation of transfer RNA (tRNA) genes with DNA replication appears to be taxon-specific. Analyzing nucleotide biases such as GC and TA skews of genes and plotting one against the other reveals a taxonomic clusterization of species. All ribosomal RNA genes are enriched in Gs and depleted in Cs, the only so far known exception being the rRNA genes of deuterostomian mitochondria. However, this exception can be explained by the fact that in the chromosome of the human mitochondrion, the model of the deuterostomian organelle genome, DNA replication, and rRNA transcription proceed in opposite directions. A general rule is deduced from prokaryotic and mitochondrial genomes: ribosomal RNA genes that are transcribed in the same direction as the DNA replication are enriched in Gs, and those transcribed in the opposite direction are depleted in Gs.