Physical Peculiarity of Two Sites in Human Promoters: Universality and Diverse Usage in Gene Function

Since the discovery of physical peculiarities around transcription start sites (TSSs) and a site corresponding to the TATA box, research has revealed only the average features of these sites. Unsettled enigmas include the individual genes with these features and whether they relate to gene function. Herein, using 10 physical properties of DNA, including duplex DNA free energy, base stacking energy, protein-induced deformability, and stabilizing energy of Z-DNA, we clarified for the first time that approximately 97% of the promoters of 21,056 human protein-coding genes have distinctive physical properties around the TSS and/or position -27; of these, nearly 65% exhibited such properties at both sites. Furthermore, about 55% of the 21,056 genes had a minimum value of regional duplex DNA free energy within TSS-centered ±300 bp regions. Notably, distinctive physical properties within the promoters and free energies of the surrounding regions separated human protein-coding genes into five groups; each contained specific gene ontology (GO) terms. The group represented by immune response genes differed distinctly from the other four regarding the parameter of the free energies of the surrounding regions. A vital suggestion from this study is that physical-feature-based analyses of genomes may reveal new aspects of the organization and regulation of genes.

The Ribosomal DNA Loci of the Ancient Monocot Pistia stratiotes L. (Araceae) Contain Different Variants of the 35S and 5S Ribosomal RNA Gene Units

The freshwater plant water lettuce (Pistia stratiotes L.) grows in warm climatic zones and is used for phytoremediation and biomass production. P. stratiotes belongs to the Araceae, an ecologically and structurally diverse early monocot family, but the phylogenetic relationships among Araceae members are poorly understood. Ribosomal DNAs (rDNAs), including the 35S and 5S rDNA, encode the RNA components of ribosomes and are widely used in phylogenetic and evolutionary studies of various plant taxa. Here, we comprehensively characterized the chromosomal locations and molecular organization of 35S and 5S rDNA genes in water lettuce using karyological and molecular methods. Fluorescence in situ hybridization revealed a single location for the 35S and 5S rDNA loci, each on a different pair of the species' 28 chromosomes. Molecular cloning and nucleotide sequencing of 35S rDNA of P. stratiotes, the first representative Araceae sensu stricto in which such a study was performed, displayed typical structural characteristics. The full-length repeat showed high sequence conservation of the regions producing the 18S, 5.8S, and 25S rRNAs and divergence of the internal transcribed spacers ITS1 and ITS2 as well as the large intergenic spacer (IGS). Alignments of the deduced sequence of 18S rDNA with the sequences available for other Araceae and representatives of other clades were used for phylogenetic analysis. Examination of 11 IGS sequences revealed significant intra-genomic length variability due to variation in subrepeat number, with four types of units detected within the 35S rDNA locus of the P. stratiotes genome (estimated size 407 Mb/1C). Similarly, the 5S rDNA locus harbors gene units comprising a conserved 119-bp sequence encoding 5S rRNA and two types of non-transcribed spacer (NTS) sequences. Type I was classified into four subtypes, which apparently originated via progressive loss of subrepeats within the duplicated NTS region containing the 3' part of the 5S rRNA gene. The minor Type II NTS is shorter than Type I and differs in nucleotide composition. Some DNA clones containing two or three consecutive 5S rDNA repeats harbored 5S rDNA genes with different types of NTSs, confirming the mosaic composition of the 5S rDNA locus.

The complete mitochondrial genome of the white seabream Diplodus sargus (Perciformes: Sparidae) from the Tyrrhenian sea

The white seabream Diplodus sargus (Linnaeus, 1758) is a species of interest for commercial fisheries throughout its range of distribution and it is also reared using aquaculture techniques. Herein, we present the first complete sequence and annotation of the mitochondrial genome of this species. The D. sargus mitogenome is 16,515 base pairs in length and contains 13 protein-coding genes, 2 rRNA, 22 tRNA, and 2 non-coding regions (D-loop and L-origin). The overall nucleotide composition is: 27.3% A, 28.9% C, 26.8% T, and 17.0% G. Maximum likelihood analyses placed D. sargus as a sister species of Diplodus puntazzo. This study provides valuable information for further studying identification methods and evolutionary relationships of Sparidae species.

Mosaic Arrangement of the 5S rDNA in the Aquatic Plant Landoltia punctata (Lemnaceae)

Duckweeds are a group of monocotyledonous aquatic plants in the Araceae superfamily, represented by 37 species divided into five genera. Duckweeds are the fastest growing flowering plants and are distributed around the globe; moreover, these plants have multiple applications, including biomass production, wastewater remediation, and making pharmaceutical proteins. Dotted duckweed (Landoltia punctata), the sole species in genus Landoltia, is one of the most resilient duckweed species. The ribosomal DNA (rDNA) encodes the RNA components of ribosomes and represents a significant part of plant genomes but has not been comprehensively studied in duckweeds. Here, we characterized the 5S rDNA genes in L. punctata by cloning and sequencing 25 PCR fragments containing the 5S rDNA repeats. No length variation was detected in the 5S rDNA gene sequence, whereas the nontranscribed spacer (NTS) varied from 151 to 524 bp. The NTS variants were grouped into two major classes, which differed both in nucleotide sequence and the type and arrangement of the spacer subrepeats. The dominant class I NTS, with a characteristic 12-bp TC-rich sequence present in 3-18 copies, was classified into four subclasses, whereas the minor class II NTS, with shorter, 9-bp nucleotide repeats, was represented by two identical sequences. In addition to these diverse subrepeats, class I and class II NTSs differed in their representation of cis-elements and the patterns of predicted G-quadruplex structures, which may influence the transcription of the 5S rDNA. Similar to related duckweed species in the genus Spirodela, L. punctata has a relatively low rDNA copy number, but in contrast to Spirodela and the majority of other plants, the arrangement of the 5S rDNA units demonstrated an unusual, heterogeneous pattern in L. punctata, as revealed by analyzing clones containing double 5S rDNA neighboring units. Our findings may further stimulate the research on the evolution of the plant rDNA and discussion of the molecular forces driving homogenization of rDNA repeats in concerted evolution.

Genome Evolutionary Dynamics Meets Functional Genomics: A Case Story on the Identification of SLC25A44

Gene clusters are becoming promising tools for gene identification. The study reveals the purposive genomic distribution of genes toward higher inheritance rates of intact metabolic pathways/phenotypes and, thereby, higher fitness. The co-localization of co-expressed, co-interacting, and functionally related genes was found as genome-wide trends in humans, mouse, golden eagle, rice fish, Drosophila, peanut, and Arabidopsis. As anticipated, the analyses verified the co-segregation of co-localized events. A negative correlation was notable between the likelihood of co-localization events and the inter-loci distances. The evolution of genomic blocks was also found convergent and uniform along the chromosomal arms. Calling a genomic block responsible for adjacent metabolic reactions is therefore recommended for identification of candidate genes and interpretation of cellular functions. As a case story, a function in the metabolism of energy and secondary metabolites was proposed for Slc25A44, based on its genomic local information. Slc25A44 was further characterized as an essential housekeeping gene which has been under evolutionary purifying pressure and belongs to the phylogenetic ETC-clade of SLC25s. Pathway enrichment mapped the Slc25A44s to the energy metabolism. The expression of peanut and human Slc25A44s in oocytes and Saccharomyces cerevisiae strains confirmed the transport of common precursors for secondary metabolites and ubiquinone. These results suggest that SLC25A44 is a mitochondrion-ER-nucleus zone transporter with biotechnological applications. Finally, a conserved three-amino acid signature on the cytosolic face of transport cavity was found important for rational engineering of SLC25s.

Insights Into the Molecular Evolution of AT-Hook Motif Nuclear Localization Genes in Brassica napus

AT-hook motif nuclear localization (AHL) proteins belong to a family of transcription factors, and play important roles in plant growth and development and response to various stresses through protein-DNA and protein-protein interactions. To better understand the Brassica napus AHL gene family, AHL genes in B. napus and related species were analyzed. Using Arabidopsis as a reference, 122 AHL gene family members were first identified in B. napus. According to the phylogenetic tree and gene organization, the BnaAHLs were classified into two clades (Clade-A and Clade-B) and three types (Type-I, Type-II, and Type-III). Gene organization and motif distribution analysis suggested that the AHL gene family is relatively conserved during evolution. These BnaAHLs are unevenly distributed on 38 chromosomes and expanded by whole-genome duplication (WGD) or segmental duplication. And large-scale loss events have also occurred in evolution. All types of BnaAHLs are subject to purification or neutral selection, while some positive selection sites are also identified in Type-II and Type-III groups. At the same time, the purification effect of Type-I members are stronger than that of the others. In addition, RNA-seq data and cis-acting element analysis also suggested that the BnaAHLs play important roles in B. napus growth and development, as well as in response to some abiotic and biotic stresses. Protein-protein interaction analysis identified some important BnaAHL-binding proteins, which also play key roles in plant growth and development. This study is helpful to fully understand the origin and evolution of the AHL gene in B. napus, and lays the foundation for their functional studies.

The complete mitochondrial genome of the sharpsnout seabream Diplodus puntazzo (Perciformes: Sparidae)

The sharpsnout seabream Diplodus puntazzo Walbaum, 1792 is a target species of small-scale fishery activities and is cage-cultured for human consumption. Nonetheless, genetic information on this species is limited. We here first sequence its complete mitochondrial genome. The sequence is composed of 16,638 base pairs, accounting for 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and 2 non-coding regions (D-loop and L-origin). The overall nucleotide composition is: 27.4% A, 28.9% C, 26.9% T, and 16.8% G. Maximum likelihood analyses placed D. puntazzo close to Acanthopagrus and some Pagellus species.

Spatiotemporal coupling and decoupling of gene transcription with DNA replication origins during embryogenesis in C. elegans

The primary task of developing embryos is genome replication, yet how DNA replication is integrated with the profound cellular changes that occur through development is largely unknown. Using an approach to map DNA replication at high resolution in C. elegans, we show that replication origins are marked with specific histone modifications that define gene enhancers. We demonstrate that the level of enhancer associated modifications scale with the efficiency at which the origin is utilized. By mapping replication origins at different developmental stages, we show that the positions and activity of origins is largely invariant through embryogenesis. Contrary to expectation, we find that replication origins are specified prior to the broad onset of zygotic transcription, yet when transcription initiates it does so in close proximity to the pre-defined replication origins. Transcription and DNA replication origins are correlated, but the association breaks down when embryonic cell division ceases. Collectively, our data indicate that replication origins are fundamental organizers and regulators of gene activity through embryonic development.

Genome-Wide Identification and Expression Profiling of Tomato Hsp20 Gene Family in Response to Biotic and Abiotic Stresses

The Hsp20 genes are involved in the response of plants to environment stresses including heat shock and also play a vital role in plant growth and development. They represent the most abundant small heat shock proteins (sHsps) in plants, but little is known about this family in tomato (Solanum lycopersicum), an important vegetable crop in the world. Here, we characterized heat shock protein 20 (SlHsp20) gene family in tomato through integration of gene structure, chromosome location, phylogenetic relationship, and expression profile. Using bioinformatics-based methods, we identified at least 42 putative SlHsp20 genes in tomato. Sequence analysis revealed that most of SlHsp20 genes possessed no intron or a relatively short intron in length. Chromosome mapping indicated that inter-arm and intra-chromosome duplication events contributed remarkably to the expansion of SlHsp20 genes. Phylogentic tree of Hsp20 genes from tomato and other plant species revealed that SlHsp20 genes were grouped into 13 subfamilies, indicating that these genes may have a common ancestor that generated diverse subfamilies prior to the mono-dicot split. In addition, expression analysis using RNA-seq in various tissues and developmental stages of cultivated tomato and the wild relative Solanum pimpinellifolium revealed that most of these genes (83%) were expressed in at least one stage from at least one genotype. Out of 42 genes, 4 genes were expressed constitutively in almost all the tissues analyzed, implying that these genes might have specific housekeeping function in tomato cell under normal growth conditions. Two SlHsp20 genes displayed differential expression levels between cultivated tomato and S. pimpinellifolium in vegetative (leaf and root) and reproductive organs (floral bud and flower), suggesting inter-species diversification for functional specialization during the process of domestication. Based on genome-wide microarray analysis, we showed that the transcript levels of SlHsp20 genes could be induced profusely by abiotic and biotic stresses such as heat, drought, salt, Botrytis cinerea, and Tomato Spotted Wilt Virus (TSWV), indicating their potential roles in mediating the response of tomato plants to environment stresses. In conclusion, these results provide valuable information for elucidating the evolutionary relationship of Hsp20 gene family and functional characterization of the SlHsp20 gene family in the future.

Nano-dissection and sequencing of DNA at single sub-nuclear structures

The relative positioning of gene loci within a mammalian nucleus is non-random and plays a role in gene regulation. Some sub-nuclear structures may represent "hubs" that bring specific genetic loci into close proximity where co-regulatory mechanisms can operate. The identification of loci in proximity to a shared sub-nuclear structure can provide insights into the function of the associated structure, and reveal relationships between the loci sharing a common association. A technique is introduced based on the nano-dissection of DNA from thin sections of cells by high-precision nano-tools operated inside a scanning electron microscope. The ability to dissect and identify gene loci occupying a shared site at a single sub-nuclear structure is demonstrated here for the first time. The technique is applied to the nano-dissection of DNA in vicinity of a single promyelocytic leukemia nuclear body (PML NB), and reveals novel loci from several chromosomes that are confirmed to associate at PML NBs with statistical significance in a cell population. Furthermore, it is demonstrated that pairs of loci from different chromosomes congregate at the same nuclear body. It is proposed that this technique is the first that allows the de novo determination of gene loci associations with single nuclear sub-structures.

The evolution and functional significance of nested gene structures in Drosophila melanogaster

Nearly 10% of the genes in the genome of Drosophila melanogaster are in nested structures, in which one gene is completely nested within the intron of another gene (nested and including gene, respectively). Even though the coding sequences and untranslated regions of these nested/including gene pairs do not overlap, their intimate structures and the possibility of shared regulatory sequences raise questions about the evolutionary forces governing the origination and subsequent functional and evolutionary impacts of these structures. In this study, we show that nested genes experience weaker evolutionary constraint, have faster rates of protein evolution, and are expressed in fewer tissues than other genes, while including genes show the opposite patterns. Surprisingly, despite completely overlapping with each other, nested and including genes are less likely to display correlated gene expression and biological function than the nearby yet nonoverlapping genes. Interestingly, significantly fewer nested genes are transcribed from the same strand as the including gene. We found that same-strand nested genes are more likely to be single-exon genes. In addition, same-strand including genes are less likely to have known lethal or sterile phenotypes than opposite-strand including genes only when the corresponding nested genes have introns. These results support our hypothesis that selection against potential erroneous mRNA splicing when nested and including genes are on the same strand plays an important role in the evolution of nested gene structures.

Molecular, pharmacological, and signaling properties of octopamine receptors from honeybee (Apis mellifera) brain

G protein-coupled receptors are important regulators of cellular signaling processes. Within the large family of rhodopsin-like receptors, those binding to biogenic amines form a discrete subgroup. Activation of biogenic amine receptors leads to transient changes of intracellular Ca²⁺-([Ca²⁺](i)) or 3',5'-cyclic adenosine monophosphate ([cAMP](i)) concentrations. Both second messengers modulate cellular signaling processes and thereby contribute to long-lasting behavioral effects in an organism. In vivo pharmacology has helped to reveal the functional effects of different biogenic amines in honeybees. The phenolamine octopamine is an important modulator of behavior. Binding of octopamine to its receptors causes elevation of [Ca²⁺](i) or [cAMP](i). To date, only one honeybee octopamine receptor that induces Ca²⁺ signals has been molecularly and pharmacologically characterized. Here, we examined the pharmacological properties of four additional honeybee octopamine receptors. When heterologously expressed, all receptors induced cAMP production after binding to octopamine with EC₅₀(s) in the nanomolar range. Receptor activity was most efficiently blocked by mianserin, a substance with antidepressant activity in vertebrates. The rank order of inhibitory potency for potential receptor antagonists was very similar on all four honeybee receptors with mianserin >> cyproheptadine > metoclopramide > chlorpromazine > phentolamine. The subroot of octopamine receptors activating adenylyl cyclases is the largest that has so far been characterized in arthropods, and it should now be possible to unravel the contribution of individual receptors to the physiology and behavior of honeybees.

Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp nasonia (hymenoptera: pteromalidae)

We sequenced the nearly complete mtDNA of 3 species of parasitic wasps, Nasonia vitripennis (2 strains), Nasonia giraulti, and Nasonia longicornis, including all 13 protein-coding genes and the 2 rRNAs, and found unusual patterns of mitochondrial evolution. The Nasonia mtDNA has a unique gene order compared with other insect mtDNAs due to multiple rearrangements. The mtDNAs of these wasps also show nucleotide substitution rates over 30 times faster than nuclear protein-coding genes, indicating among the highest substitution rates found in animal mitochondria (normally <10 times faster). A McDonald and Kreitman test shows that the between-species frequency of fixed replacement sites relative to silent sites is significantly higher compared with within-species polymorphisms in 2 mitochondrial genes of Nasonia, atp6 and atp8, indicating directional selection. Consistent with this interpretation, the Ka/Ks (nonsynonymous/synonymous substitution rates) ratios are higher between species than within species. In contrast, cox1 shows a signature of purifying selection for amino acid sequence conservation, although rates of amino acid substitutions are still higher than for comparable insects. The mitochondrial-encoded polypeptides atp6 and atp8 both occur in F0F1ATP synthase of the electron transport chain. Because malfunction in this fundamental protein severely affects fitness, we suggest that the accelerated accumulation of replacements is due to beneficial mutations necessary to compensate mild-deleterious mutations fixed by random genetic drift or Wolbachia sweeps in the fast evolving mitochondria of Nasonia. We further propose that relatively high rates of amino acid substitution in some mitochondrial genes can be driven by a "Compensation-Draft Feedback"; increased fixation of mildly deleterious mutations results in selection for compensatory mutations, which lead to fixation of additional deleterious mutations in nonrecombining mitochondrial genomes, thus accelerating the process of amino acid substitutions.

Genome organization in three dimensions: thinking outside the line

Gene organization on nuclear chromosomes is usually depicted as a linear array, but at least some regions of the genome are localized to specific subnuclear positions in interphase nuclei. Studies in yeast have found that centromeres and telomeres are found around the nuclear periphery, and that tRNA genes are gathered at the nucleolus, along with the ribosomal RNA gene cluster. These 325 loci alone impose significant constraints on the three dimensional organization of chromosomes in the nucleus, and there is mounting experimental evidence that transcription by RNA polymerase II is strongly affected by proximity to these regions. Given these observations, one consideration in understanding nuclear gene regulation might be the degree to which spatial positioning affects at least a subset of gene families.

Intragenic DNA spacers interrupt the ovalbumin gene

We have performed restriction endonuclease mapping to examine the linear organization of the ovalbumin gene in chromosomal DNA. Treatment of genomic DNA with restriction endonucleases that do not cleave the ovalbumin mRNA sequence results in the generation of multiple DNA fragments capable of annealing with ovalbumin-specific probes in molecular hybridization reactions. These data strongly suggest that the linear order of DNA sequences coding for ovalbumin is interrupted by at least two intragenic DNA spacers absent from the corresponding RNA. At least one of these spacer sequences interrupts the coding sequence; therefore the chromosomal ovalbumin gene is not colinear with its translational product. We can discern no difference in the sequence organization about this gene in producer and nonproducer somatic cells, suggesting that the presence of intragenic spacers does not reflect that transcriptional activity of this gene. Furthermore, profiles obtained for gametes are identical to those observed for somatic cells, indicating that significant translocation during development is not responsible for the generation of this split sequence. Divergence of the intragenic spacer, however, is observed between individual chickens; thus multiple alleles may exist for this gene that are identifiable by differences in the organization of spacer with no apparent phenotypic evidence for their presence.