Evolutionary significance of gene expression divergence

Recent large-scale studies of evolutionary changes in gene expression among mammalian species have led to the proposal that gene expression divergence may be neutral with respect to organismic fitness. Here, we employ a comparative analysis of mammalian gene sequence divergence and gene expression divergence to test the hypothesis that the evolution of gene expression is predominantly neutral. Two models of neutral gene expression evolution are considered: 1-purely neutral evolution (i.e., no selective constraint) of gene expression levels and patterns and 2-neutral evolution accompanied by selective constraint. With respect to purely neutral evolution, levels of change in gene expression between human-mouse orthologs are correlated with levels of gene sequence divergence that are determined largely by purifying selection. In contrast, evolutionary changes of tissue-specific gene expression profiles do not show such a correlation with sequence divergence. However, divergence of both gene expression levels and profiles are significantly lower for orthologous human-mouse gene pairs than for pairs of randomly chosen human and mouse genes. These data clearly point to the action of selective constraint on gene expression divergence and are inconsistent with the purely neutral model; however, there is likely to be a neutral component in evolution of gene expression, particularly, in tissues where the expression of a given gene is low and functionally irrelevant. The model of neutral evolution with selective constraint predicts a regular, clock-like accumulation of gene expression divergence. However, relative rate tests of the divergence among human-mouse-rat orthologous gene sets reveal clock-like evolution for gene sequence divergence, and to a lesser extent for gene expression level divergence, but not for the divergence of tissue-specific gene expression profiles. Taken together, these results indicate that gene expression divergence is subject to the effects of purifying selective constraint and suggest that it might also be substantially influenced by positive Darwinian selection.     

Evolutionary analysis of PHLPP1 gene in humans and non-human primates

The chromosome 18q22-23 region has been shown to be implicated in bipolar disorder (BPAD) by several studies. PHLPP1 gene, in the locus (chromosome 18q22-23), is involved in circadian pathways and bears modules like ‘PH domain and leucine rich repeat protein phosphatase’. This gene also contains a polyglutamine (CAG or PolyQ) repeat motif at the carboxyl terminal end. A comparative analysis of the PolyQ repeats of the PHLPP1 gene in humans, non-human primates and other species has been attempted in order to investigate the possible significance of repeat length as seen in other triplet-repeat associated diseases. Sequencing of the CAG repeat in humans and in non-human primates revealed that the CAG repeat is not polymorphic in humans; whereas, in other species it shows an area of high variability, both in length and sequence composition. Despite the conservation of circadian clock components in different species, there is remarkable diversity in the protein structure, regulation and biochemical functions of the circadian orthologs. These can be due to specific adaptations in accordance with the physiology of the particular species providing a species-specific biological advantage.     

Genome-wide associations of gene expression variation in humans

The exploration of quantitative variation in human populations has become one of the major priorities for medical genetics. The successful identification of variants that contribute to complex traits is highly dependent on reliable assays and genetic maps. We have performed a genome-wide quantitative trait analysis of 630 genes in 60 unrelated Utah residents with ancestry from Northern and Western Europe using the publicly available phase I data of the International HapMap project. The genes are located in regions of the human genome with elevated functional annotation and disease interest including the ENCODE regions spanning 1% of the genome, Chromosome 21 and Chromosome 20q12–13.2. We apply three different methods of multiple test correction, including Bonferroni, false discovery rate, and permutations. For the 374 expressed genes, we find many regions with statistically significant association of single nucleotide polymorphisms (SNPs) with expression variation in lymphoblastoid cell lines after correcting for multiple tests. Based on our analyses, the signal proximal (cis-) to the genes of interest is more abundant and more stable than distal and trans across statistical methodologies. Our results suggest that regulatory polymorphism is widespread in the human genome and show that the 5-kb (phase I) HapMap has sufficient density to enable linkage disequilibrium mapping in humans. Such studies will significantly enhance our ability to annotate the non-coding part of the genome and interpret functional variation. In addition, we demonstrate that the HapMap cell lines themselves may serve as a useful resource for quantitative measurements at the cellular level.     

Towards safe, non-viral therapeutic gene expression in humans

The potential dangers of using viruses to deliver and integrate DNA into host cells in gene therapy have been poignantly highlighted in recent clinical trials. Safer, non-viral gene delivery approaches have been largely ignored in the past because of their inefficient delivery and the resulting transient transgene expression. However, recent advances indicate that efficient, long-term gene expression can be achieved by non-viral means. In particular, integration of DNA can be targeted to specific genomic sites without deleterious consequences and it is possible to maintain transgenes as small episomal plasmids or artificial chromosomes. The application of these approaches to human gene therapy is gradually becoming a reality.     

Evolutionary rate and gene expression across different brain regions

Background

The evolutionary rate of a protein is a basic measure of evolution at the molecular level. Previous studies have shown that genes expressed in the brain have significantly lower evolutionary rates than those expressed in somatic tissues.

Results

We study the evolutionary rates of genes expressed in 21 different human brain regions. We find that genes highly expressed in the more recent cortical regions of the brain have lower evolutionary rates than genes highly expressed in subcortical regions. This may partially result from the observation that genes that are highly expressed in cortical regions tend to be highly expressed in subcortical regions, and thus their evolution faces a richer set of functional constraints. The frequency of mammal-specific and primate-specific genes is higher in the highly expressed gene sets of subcortical brain regions than in those of cortical brain regions. The basic inverse correlation between evolutionary rate and gene expression is significantly stronger in brain versus nonbrain tissues, and in cortical versus subcortical regions. Extending upon this cortical/subcortical trend, this inverse correlation is generally more marked for tissues that are located higher along the cranial vertical axis during development, giving rise to the possibility that these tissues are also more evolutionarily recent.

Conclusions

We find that cortically expressed genes are more conserved than subcortical ones, and that gene expression levels exert stronger constraints on sequence evolution in cortical versus subcortical regions. Taken together, these findings suggest that cortically expressed genes are under stronger selective pressure than subcortically expressed genes.     

Evolutionary diversification of SPANX-N sperm protein gene structure and expression

The sperm protein associated with nucleus in the X chromosome (SPANX) genes cluster at Xq27 in two subfamilies, SPANX-A/D and SPANX-N. SPANX-A/D is specific for hominoids and is fairly well characterized. The SPANX-N gave rise to SPANX-A/D in the hominoid lineage approximately 7 MYA. Given the proposed role of SPANX genes in spermatogenesis, we have extended studies to SPANX-N gene evolution, variation, regulation of expression, and intra-sperm localization. By immunofluorescence analysis, SPANX-N proteins are localized in post-meiotic spermatids exclusively, like SPANX-A/D. But in contrast to SPANX-A/D, SPANX-N are found in all ejaculated spermatozoa rather than only in a subpopulation, are localized in the acrosome rather than in the nuclear envelope, and are expressed at a low level in several nongametogenic adult tissues as well as many cancers. Presence of a binding site for CTCF and its testis-specific paralogue BORIS in the SPANX promoters suggests, by analogy to MAGE-A1 and NY-ESO-1, that their activation in spermatogenesis is mediated by the programmed replacement of CTCF by BORIS. Based on the relative density of CpG, the more extended expression of SPANX-N compared to SPANX-A/D in nongametogenic tissues is likely attributed to differences in promoter methylation. Our findings suggest that the recent duplication of SPANX genes in hominoids was accompanied by different localization of SPANX-N proteins in post-meiotic sperm and additional expression in several nongonadal tissues. This suggests a corresponding functional diversification of SPANX gene families in hominoids. SPANX proteins thus provide unique targets to investigate their roles in the function of spermatozoa, selected malignancies, and for SPANX-N, in other tissues as well.     

Evolutionary ecology of opsin gene sequence,expression and repertoire

Linking molecular evolution to biological function is a long‐standing challenge in evolutionary biology. Some of the best examples of this involve opsins, the genes that encode the molecular basis of light reception. In this issue of Molecular Ecology, three studies examine opsin gene sequence, expression and repertoire to determine how natural selection has shaped the visual system. First, Escobar‐Camacho et al. ( 2017 ) use opsin repertoire and expression in three Amazonian cichlid species to show that a shift in sensitivity towards longer wavelengths is coincident with the long‐wavelength‐dominated Amazon basin. Second, Stieb et al. ( 2017 ) explore opsin sequence and expression in reef‐dwelling damselfish and find that UV‐ and long‐wavelength vision are both important, but likely for different ecological functions. Lastly, Suvorov et al. ( 2017 ) study an expansive opsin repertoire in the insect order Odonata and find evidence that copy number expansion is consistent with the permanent heterozygote model of gene duplication. Together these studies emphasize the utility of opsin genes for studying both the local adaptation of sensory systems and, more generally, gene family evolution.     

Regulatory implications of translational frameshifting in cellular gene expression

The genetic code, once thought to be rigid, hag been found to be quite fiexible, permitting several different reading alternatives. One of these is translatlonal frameshifting, a process programmed in the mRNA sequence and which enables a +1 or -1 shift from the reading frame of the initiation codon. So far, the Involvement of translatlonal frameahifting in gene expression has been described mainly in viruses (particularly retroviruses), retrotransposons, and bacterial insertion elements, in this MicroReview., we present a survey of the cellular genes, mostly in Escherichia coil, which have been found to be expressed through a transiational frameshifting process, as well as a discussion of the regulatory implications of this process.     

Evolutionary constraints and expression analysis of gene duplications in Rhodobacter Sphaeroides 2.4.1

ABSTRACT: BACKGROUND: Gene duplication is a major force that contributes to the evolution of new metabolic functions in all organisms. Rhodobacter sphaeroides 2.4.1 is a bacterium that displays a wide degree of metabolic versatility and genome complexity and therefore is a fitting model for the study of gene duplications in bacteria. A comprehensive analysis of 234 duplicate gene-pairs in R. sphaeroides was performed using structural constraint and expression analysis. RESULTS: The results revealed that most gene-pairs in in-paralogs are maintained under negative selection (omega [less than or equal to] 0.3), but the strength of selection differed among in-paralog gene-pairs. Although in-paralogs located on different replicons are maintained under purifying selection, the duplicated genes distributed between the primary chromosome (CI) and the second chromosome (CII) are relatively less selectively constrained than the gene-pairs located within each chromosome. The mRNA expression patterns of duplicate gene-pairs were examined through microarray analysis of this organism grown under seven different growth conditions. Results revealed that that ~62% of paralogs have similar expression patterns (cosine [greater than or equal to] 0.90) over all of these growth conditions, while only ~7% of paralogs are very different in their expression patterns (cosine < 0.50). CONCLUSIONS: The overall findings of the study suggest that only a small proportion of paralogs contribute to the metabolic diversity and the evolution of novel metabolic functions in R. sphaeroides. In addition, the lack of relationships between structural constraints and gene-pair expression suggests that patterns of gene-pair expression are likely associated with conservation or divergence of gene-pair promoter regions and other coregulation mechanisms.     

Regulation of Fto/Ftm gene expression in mice and humans

Two recent, large whole-genome association studies (GWAS) in European populations have associated a approximately 47-kb region that contains part of the FTO gene with high body mass index (BMI). The functions of FTO and adjacent FTM in human biology are not clear. We examined expression of these genes in organs of mice segregating for monogenic obesity mutations, exposed to underfeeding/overfeeding, and to 4 degrees C. Fto/Ftm expression was reduced in mesenteric adipose tissue of mice segregating for the Ay, Lep ob, Lepr db, Cpe fat, or tub mutations, and there was a similar trend in other tissues. These effects were not due to adiposity per se. Hypothalamic Fto and Ftm expression were decreased by fasting in lean and obese animals and by cold exposure in lean mice. The fact that responses of Fto and Ftm expression to these manipulations were almost indistinguishable suggested that the genes might be coregulated. The putative overlapping regulatory region contains at least two canonical CUTL1 binding sites. One of these nominal CUTL1 sites includes rs8050136, a SNP associated with high body mass. The A allele of rs8050136 preferentially bound CUTL1[corrected] in human fibroblast DNA. 70% knockdown of CUTL1 expression in human fibroblasts decreased FTO and FTM expression by 90 and 65%, respectively. Animals and humans with various genetic interruptions of FTO or FTM have phenotypes reminiscent of aspects of the Bardet-Biedl obesity syndrome, a confirmed "ciliopathy." FTM has recently been shown to be a ciliary basal body protein.     

Evolutionary implications of matK indels in Poaceae

Insertion/deletion events (indels) and nucleotide substitutions at the extreme 3' end of the chloroplast gene matK have been identified that distinguish certain major lineages of grasses. A 1-bp (base pair) deletion creating a shift in the open reading frame (ORF) and a point mutation support the positions of Streptochaeta and Anomochloa as the two most basal lineages in Poaceae. Another 1-bp deletion resulting in early termination of the ORF is unique to Ehrharta, a member of the taxonomically disputable tribe Ehrharteae. A 6-bp insertion supports monophyly of subfamilies Panicoideae, Arundinoideae, Centothecoideae, and Chloridoideae (PACC). This marker appears useful in defining PACC clade members and may have potential in providing insight into the sister-group relationship between PACC and other lineages. Alignment of deduced amino acid sequences from bryophytes, gymnosperms, and angiosperms shows that this region is relatively conserved, but variation is notably higher in Poaceae. The evolutionary implications of these changes in grasses and other plant families are addressed.     

Epigene conversion: a proposal with implications for gene mapping in humans.

Epigenetic modification of DNA is now recognized as a potentially important factor in the inheritance and expression of some mutations; its ability to complicate human genetic analysis is concurrently becoming apparent. One unusual form of epigenetic modification, dominant position-effect variegation (PEV), has been used as a model for Huntington disease. In dominant PEV, a fully dominant mutant phenotype results from stable epigenetic inactivation of an allele adjacent to the structural alteration (cis-inactivation) combined with a complementary inactivation of the homologous normal allele (trans-inactivation). We now propose that trans-inactivation of the normal allele may occasionally persist through meiosis. Such "epigene conversion" occurring at the Huntington disease locus in a few percent of meioses would largely account for the published anomalies in that region's genetic map. This concept could also explain anomalous linkage map data for other disease-causing alleles in humans.     

Evolutionary and diagnostic implications of intragenomic heterogeneity in the 16S rRNA gene in Aeromonas strains

Sequencing 16S rRNA genes (SSU) cloned from Aeromonas strains revealed that strains contained up to six copies differing by < or = 1.5%. The SSU copies from Aeromonas veronii LMG13695 clustered with sequences from four Aeromonas species. These results demonstrate intragenomic heterogeneity of SSU and suggest caution when using SSU to identify aeromonads.     

Evolutionary implications of bacterial polyketide synthases

Polyketide synthases (PKS) perform a stepwise biosynthesis of diverse carbon skeletons from simple activated carboxylic acid units. The products of the complex pathways possess a wide range of pharmaceutical properties, including antibiotic, antitumor, antifungal, and immunosuppressive activities. We have performed a comprehensive phylogenetic analysis of multimodular and iterative PKS of bacteria and fungi and of the distinct types of fatty acid synthases (FAS) from different groups of organisms based on the highly conserved ketoacyl synthase (KS) domains. Apart from enzymes that meet the classification standards we have included enzymes involved in the biosynthesis of mycolic acids, polyunsaturated fatty acids (PUFA), and glycolipids in bacteria. This study has revealed that PKS and FAS have passed through a long joint evolution process, in which modular PKS have a central position. They appear to have derived from bacterial FAS and primary iterative PKS and, in addition, share a common ancestor with animal FAS and secondary iterative PKS. Furthermore, we have carried out a phylogenomic analysis of all modular PKS that are encoded by the complete eubacterial genomes currently available in the database. The phylogenetic distribution of acyltransferase and KS domain sequences revealed that multiple gene duplications, gene losses, as well as horizontal gene transfer (HGT) have contributed to the evolution of PKS I in bacteria. The impact of these factors seems to vary considerably between the bacterial groups. Whereas in actinobacteria and cyanobacteria the majority of PKS I genes may have evolved from a common ancestor, several lines of evidence indicate that HGT has strongly contributed to the evolution of PKS I in proteobacteria. Discovery of new evolutionary links between PKS and FAS and between the different PKS pathways in bacteria may help us in understanding the selective advantage that has led to the evolution of multiple secondary metabolite biosyntheses within individual bacteria.