Evolutionary genetics: making the most of redundancy

Telomerase RNA gene paralogs in plants - the usual pathway to unusual telomeres

Telomerase, telomeric DNA and associated proteins together represent a complex, finely tuned and functionally conserved mechanism that ensures genome integrity by protecting and maintaining chromosome ends. Changes in its components can threaten an organism's viability. Nevertheless, molecular innovation in telomere maintenance has occurred multiple times during eukaryote evolution, giving rise to species/taxa with unusual telomeric DNA sequences, telomerase components or telomerase-independent telomere maintenance. The central component of telomere maintenance machinery is telomerase RNA (TR) as it templates telomere DNA synthesis, its mutation can change telomere DNA and disrupt its recognition by telomere proteins, thereby leading to collapse of their end-protective and telomerase recruitment functions. Using a combination of bioinformatic and experimental approaches, we examine a plausible scenario of evolutionary changes in TR underlying telomere transitions. We identified plants harbouring multiple TR paralogs whose template regions could support the synthesis of diverse telomeres. In our hypothesis, formation of unusual telomeres is associated with the occurrence of TR paralogs that can accumulate mutations, and through their functional redundancy, allow for the adaptive evolution of the other telomere components. Experimental analyses of telomeres in the examined plants demonstrate evolutionary telomere transitions corresponding to TR paralogs with diverse template regions.

Variation in the genomic basis of parallel phenotypic and ecological divergence in benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Sympatric adaptive phenotypic divergence should be underlain by genomic differentiation between sub-populations. When divergence drives similar patterns of phenotypic and ecological variation within species we expect evolution to draw on common allelic variation. We investigated divergence histories and genomic signatures of adaptive divergence between benthic and pelagic morphs of Icelandic Arctic charr. Divergence histories for each of four populations were reconstructed using coalescent modelling and 14,187 single nucleotide polymorphisms. Sympatric divergence with continuous gene flow was supported in two populations while allopatric divergence with secondary contact was supported in one population; we could not differentiate between demographic models in the fourth population. We detected parallel patterns of phenotypic divergence along benthic-pelagic evolutionary trajectories among populations. Patterns of genomic differentiation between benthic and pelagic morphs were characterized by outlier loci in many narrow peaks of differentiation throughout the genome, which may reflect the eroding effects of gene flow on nearby neutral loci. We then used genome-wide association analyses to relate both phenotypic (body shape and size) and ecological (carbon and nitrogen stable isotopes) variation to patterns of genomic differentiation. Many peaks of genomic differentiation were associated with phenotypic and ecological variation in the three highly divergent populations, suggesting a genomic basis for adaptive divergence. We detected little evidence for a parallel genomic basis of differentiation as most regions and outlier loci were not shared among populations. Our results show that adaptive divergence can have varied genomic consequences in populations with relatively recent common origins, similar divergence histories, and parallel phenotypic divergence.

Retention of fatty acyl desaturase 1 (fads1) in Elopomorpha and Cyclostomata provides novel insights into the evolution of long-chain polyunsaturated fatty acid biosynthesis in vertebrates

Background

Provision of long-chain polyunsaturated fatty acids (LC-PUFA) in vertebrates occurs through the diet or via endogenous production from C₁₈ precursors through consecutive elongations and desaturations. It has been postulated that the abundance of LC-PUFA in the marine environment has remarkably modulated the gene complement and function of Fads in marine teleosts. In vertebrates two fatty acyl desaturases, namely Fads1 and Fads2, encode ∆5 and ∆6 desaturases, respectively. To fully clarify the evolutionary history of LC-PUFA biosynthesis in vertebrates, we investigated the gene repertoire and function of Fads from species placed at key evolutionary nodes.

Results

We demonstrate that functional Fads1Δ5 and Fads2∆6 arose from a tandem gene duplication in the ancestor of vertebrates, since they are present in the Arctic lamprey. Additionally, we show that a similar condition was retained in ray-finned fish such as the Senegal bichir and spotted gar, with the identification of fads1 genes in these lineages. Functional characterisation of the isolated desaturases reveals the first case of a Fads1 enzyme with ∆5 desaturase activity in the Teleostei lineage, the Elopomorpha. In contrast, in Osteoglossomorpha genomes, while no fads1 was identified, two separate fads2 duplicates with ∆6 and ∆5 desaturase activities respectively were uncovered.

Conclusions

We conclude that, while the essential genetic components involved LC-PUFA biosynthesis evolved in the vertebrate ancestor, the full completion of the LC-PUFA biosynthesis pathway arose uniquely in gnathostomes.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1271-5) contains supplementary material, which is available to authorized users.

Evolution of the ability to modulate host chemokine networks via gene duplication in human cytomegalovirus (HCMV)

Human cytomegalovirus (HCMV) is a widespread pathogen that is particularly skillful at evading immune detection and defense mechanisms, largely due to extensive co-evolution with its host. One aspect of this co-evolution involves the acquisition of virally encoded G protein-coupled receptors (GPCRs) with homology to the chemokine receptor family. GPCRs are the largest family of cell surface proteins, found in organisms from yeast to humans, and they regulate a variety of cellular processes including development, sensory perception, and immune cell trafficking. The US27 and US28 genes are encoded by human and primate CMVs, but homologs are not found in the genomes of viruses infecting rodents or other species. Phylogenetic analysis was used to investigate the US27 and US28 genes, which are adjacent in the unique short (US) region of the HCMV genome, and their relationship to one another and to human chemokine receptor genes. The results indicate that both US27 and US28 share the same common ancestor with human chemokine receptor CX3CR1, suggesting that a single host gene was captured and a subsequent viral gene duplication event occurred. The US28 gene product (pUS28) has maintained the function of the ancestral gene and has the ability to bind and signal in response to CX3CL1/fractalkine, the natural ligand for CX3CR1. In contrast, pUS27 does not bind to any known chemokine ligand, and the sequence has diverged significantly, highlighted by the fact that pUS27 currently exhibits greater sequence similarity to human CCR1. While the evolutionary advantage of the gene duplication and neofunctionalization event remains unclear, the US27 and US28 genes are highly conserved among different HCMV strains and retained even in laboratory strains that have lost many virulence genes, suggesting that US27 and US28 have each evolved distinct, important functions during virus infection.

Analyses of Old "Prokaryotic" Proteins Indicate Functional Diversification in Arabidopsis and Oryza sativa

During evolution, various processes such as duplication, divergence, recombination, and many other events leads to the evolution of new genes with novel functions. These evolutionary events, thus significantly impact the evolution of cellular, physiological, morphological, and other phenotypic trait of organisms. While evolving, eukaryotes have acquired large number of genes from the earlier prokaryotes. This work is focused upon identification of old "prokaryotic" proteins in Arabidopsis and Oryza sativa genome, further highlighting their possible role(s) in the two genomes. Our results suggest that with respect to their genome size, the fraction of old "prokaryotic" proteins is higher in Arabidopsis than in Oryza sativa. The large fractions of such proteins encoding genes were found to be localized in various endo-symbiotic organelles. The domain architecture of the old "prokaryotic" proteins revealed similar distribution in both Arabidopsis and Oryza sativa genomes showing their conserved evolution. In Oryza sativa, the old "prokaryotic" proteins were more involved in developmental processes, might be due to constant man-made selection pressure for better agronomic traits/productivity. While in Arabidopsis, these proteins were involved in metabolic functions. Overall, the analysis indicates the distinct pattern of evolution of old "prokaryotic" proteins in Arabidopsis and Oryza sativa.

Molecular Evolution and Functional Diversification of Replication Protein A1 in Plants

Replication protein A (RPA) is a heterotrimeric, single-stranded DNA binding complex required for eukaryotic DNA replication, repair, and recombination. RPA is composed of three subunits, RPA1, RPA2, and RPA3. In contrast to single RPA subunit genes generally found in animals and yeast, plants encode multiple paralogs of RPA subunits, suggesting subfunctionalization. Genetic analysis demonstrates that five Arabidopsis thaliana RPA1 paralogs (RPA1A to RPA1E) have unique and overlapping functions in DNA replication, repair, and meiosis. We hypothesize here that RPA1 subfunctionalities will be reflected in major structural and sequence differences among the paralogs. To address this, we analyzed amino acid and nucleotide sequences of RPA1 paralogs from 25 complete genomes representing a wide spectrum of plants and unicellular green algae. We find here that the plant RPA1 gene family is divided into three general groups termed RPA1A, RPA1B, and RPA1C, which likely arose from two progenitor groups in unicellular green algae. In the family Brassicaceae the RPA1B and RPA1C groups have further expanded to include two unique sub-functional paralogs RPA1D and RPA1E, respectively. In addition, RPA1 groups have unique domains, motifs, cis-elements, gene expression profiles, and pattern of conservation that are consistent with proposed functions in monocot and dicot species, including a novel C-terminal zinc-finger domain found only in plant RPA1C-like sequences. These results allow for improved prediction of RPA1 subunit functions in newly sequenced plant genomes, and potentially provide a unique molecular tool to improve classification of Brassicaceae species.

Bridging sequence diversity and tissue-specific expression by DNA methylation in genes of the mouse prolactin superfamily

Much of the DNA in genomes is organized within gene families and hierarchies of gene superfamilies. DNA methylation is the main epigenetic event involved in gene silencing and genome stability. In the present study, we analyzed the DNA methylation status of the prolactin (PRL) superfamily to obtain insight into its tissue-specific expression and the evolution of its sequence diversity. The PRL superfamily in mice consists of two dozen members, which are expressed in a tissue-specific manner. The genes in this family have CpG-less sequences, and they are located within a 1-Mb region as a gene cluster on chromosome 13. We tentatively grouped the family into several gene clusters, depending on location and gene orientation. We found that all the members had tissue-dependent differentially methylated regions (T-DMRs) around the transcription start site. The T-DMRs are hypermethylated in nonexpressing tissues and hypomethylated in expressing cells, supporting the idea that the expression of the PRL superfamily genes is subject to epigenetic regulation. Interestingly, the DNA methylation patterns of T-DMRs are shared within a cluster, while the patterns are different among the clusters. Finally, we reconstituted the nucleotide sequences of T-DMRs by converting TpG to CpG based on the consideration of a possible conversion of 5-methylcytosine to thymine by spontaneous deamination during the evolutionary process. On the phylogenic tree, the reconstituted sequences were well matched with the DNA methylation pattern of T-DMR and orientation. Our study suggests that DNA methylation is involved in tissue-specific expression and sequence diversity during evolution.

Molecular evolution of miraculin-like proteins in soybean Kunitz super-family

Miraculin-like proteins (MLPs) belong to soybean Kunitz super-family and have been characterized from many plant families like Rutaceae, Solanaceae, Rubiaceae, etc. Many of them possess trypsin inhibitory activity and are involved in plant defense. MLPs exhibit significant sequence identity (~30-95%) to native miraculin protein, also belonging to Kunitz super-family compared with a typical Kunitz family member (~30%). The sequence and structure-function comparison of MLPs with that of a classical Kunitz inhibitor have demonstrated that MLPs have evolved to form a distinct group within Kunitz super-family. Sequence analysis of new genes along with available MLP sequences in the literature revealed three major groups for these proteins. A significant feature of Rutaceae MLP type 2 sequences is the presence of phosphorylation motif. Subtle changes are seen in putative reactive loop residues among different MLPs suggesting altered specificities to specific proteases. In phylogenetic analysis, Rutaceae MLP type 1 and type 2 proteins clustered together on separate branches, whereas native miraculin along with other MLPs formed distinct clusters. Site-specific positive Darwinian selection was observed at many sites in both the groups of Rutaceae MLP sequences with most of the residues undergoing positive selection located in loop regions. The results demonstrate the sequence and thereby the structure-function divergence of MLPs as a distinct group within soybean Kunitz super-family due to biotic and abiotic stresses of local environment.

Functions of the duplicated hik31 operons in central metabolism and responses to light, dark, and carbon sources in Synechocystis sp. strain PCC 6803

There are two closely related hik31 operons involved in signal transduction on the chromosome and the pSYSX plasmid in the cyanobacterium Synechocystis sp. strain PCC 6803. We studied the growth, cell morphology, and gene expression in operon and hik mutants for both copies, under different growth conditions, to examine whether the duplicated copies have the same or different functions and gene targets and whether they are similarly regulated. Phenotype analysis suggested that both operons regulated common and separate targets in the light and the dark. The chromosomal operon was involved in the negative control of autotrophic events, whereas the plasmid operon was involved in the positive control of heterotrophic events. Both the plasmid and double operon mutant cells were larger and had division defects. The growth data also showed a regulatory role for the chromosomal hik gene under high-CO(2) conditions and the plasmid operon under low-O(2) conditions. Metal stress experiments indicated a role for the chromosomal hik gene and operon in mediating Zn and Cd tolerance, the plasmid operon in Co tolerance, and the chromosomal operon and plasmid hik gene in Ni tolerance. We conclude that both operons are differentially and temporally regulated. We suggest that the chromosomal operon is the primarily expressed copy and the plasmid operon acts as a backup to maintain appropriate gene dosages. Both operons share an integrated regulatory relationship and are induced in high light, in glucose, and in active cell growth. Additionally, the plasmid operon is induced in the dark with or without glucose.

Systematic analysis of potato acid invertase genes reveals that a cold-responsive member, StvacINV1, regulates cold-induced sweetening of tubers

Acid invertase is believed to play a regulatory role during plant developmental processes and to respond to environmental stimuli. The expression profiles of the entire acid invertase family are not yet available for potato. By searching existing databases, it was determined that there are at least six acid invertase genes in potato, including four cell-wall invertase genes and two vacuolar invertase genes. They were subjected to comparative expression profiling in various organs of potato plants and in stored tubers to exploit their potential functions. The results revealed that each gene exhibited a unique expression pattern, which differed in transcript abundance or showed organ-specific features, pointing to the possible involvement of individual genes in plant development. The vacuolar invertase gene StvacINV1 had the highest expression level among three genes detected in the potato tubers. Further storage experiments showed that StvacINV1 was strongly induced by low temperatures, which is consistent with glucose accumulation in cold-stored tubers. Suppression of StvacINV1 by the antisense transformation in potato confirmed that lower StvacINV1 transcript abundance in transgenic tubers is related to lower reducing sugar content and lighter chip color in comparison with the wild type. The evidence strongly suggests that StvacINV1 is a gene involved in regulation of cold-induced sweetening of potato tubers. This provides an avenue for studying the mechanism involved in the regulation of the cold-induced sweetening trait and for agronomic enhancement.

Evolutionary origins of Brassicaceae specific genes in Arabidopsis thaliana

Background

All sequenced genomes contain a proportion of lineage-specific genes, which exhibit no sequence similarity to any genes outside the lineage. Despite their prevalence, the origins and functions of most lineage-specific genes remain largely unknown. As more genomes are sequenced opportunities for understanding evolutionary origins and functions of lineage-specific genes are increasing.

Results

This study provides a comprehensive analysis of the origins of lineage-specific genes (LSGs) in Arabidopsis thaliana that are restricted to the Brassicaceae family. In this study, lineage-specific genes within the nuclear (1761 genes) and mitochondrial (28 genes) genomes are identified. The evolutionary origins of two thirds of the lineage-specific genes within the Arabidopsis thaliana genome are also identified. Almost a quarter of lineage-specific genes originate from non-lineage-specific paralogs, while the origins of ~10% of lineage-specific genes are partly derived from DNA exapted from transposable elements (twice the proportion observed for non-lineage-specific genes). Lineage-specific genes are also enriched in genes that have overlapping CDS, which is consistent with such novel genes arising from overprinting. Over half of the subset of the 958 lineage-specific genes found only in Arabidopsis thaliana have alignments to intergenic regions in Arabidopsis lyrata, consistent with either de novo origination or differential gene loss and retention, with both evolutionary scenarios explaining the lineage-specific status of these genes. A smaller number of lineage-specific genes with an incomplete open reading frame across different Arabidopsis thaliana accessions are further identified as accession-specific genes, most likely of recent origin in Arabidopsis thaliana. Putative de novo origination for two of the Arabidopsis thaliana-only genes is identified via additional sequencing across accessions of Arabidopsis thaliana and closely related sister species lineages. We demonstrate that lineage-specific genes have high tissue specificity and low expression levels across multiple tissues and developmental stages. Finally, stress responsiveness is identified as a distinct feature of Brassicaceae-specific genes; where these LSGs are enriched for genes responsive to a wide range of abiotic stresses.

Conclusion

Improving our understanding of the origins of lineage-specific genes is key to gaining insights regarding how novel genes can arise and acquire functionality in different lineages. This study comprehensively identifies all of the Brassicaceae-specific genes in Arabidopsis thaliana and identifies how the majority of such lineage-specific genes have arisen. The analysis allows the relative importance (and prevalence) of different evolutionary routes to the genesis of novel ORFs within lineages to be assessed. Insights regarding the functional roles of lineage-specific genes are further advanced through identification of enrichment for stress responsiveness in lineage-specific genes, highlighting their likely importance for environmental adaptation strategies.

Evolution and function of the uterine serpins (SERPINA14)

Uterine serpins (recently designated as SERPINA14) are hormonally induced proteins secreted in large quantities by the endometrial epithelium during pregnancy. The SERPINA14 proteins belong to the serine proteinase inhibitor (serpin) superfamily, but their apparent lack of inhibitory activity toward serine proteinases suggests that these proteins evolved a different function from the anti-proteinase activity typically found in most members of the serpin superfamily. The gene is present in a limited group of mammals in the Laurasiatheria superorder (ruminants, horses, pigs, dolphins and some carnivores) while being absent in primates, rodents, lagomorphs and marsupials. Thus, the gene is likely to have evolved by gene duplication after divergence of Laurasiatheria and to play an important role in pregnancy. That role may vary between species. In sheep, SERPINA14 probably serves an immunoregulatory role to prevent rejection of the fetal allograft. It is inhibitory to lymphocyte proliferation and natural killer cell function. In the pig, SERPINA14 is involved in iron transport to the fetus by binding to and stabilizing the iron-binding protein uteroferrin. It is possible that SERPINA14 has undergone divergence in function since the original emergence of the gene in a common ancestor of species possessing SERPINA14.

The evolution of skeletal muscle performance: gene duplication and divergence of human sarcomeric alpha-actinins

In humans, there are two skeletal muscle alpha-actinins, encoded by ACTN2 and ACTN3, and the ACTN3 genotype is associated with human athletic performance. Remarkably, approximately 1 billion people worldwide are deficient in alpha-actinin-3 due to the common ACTN3 R577X polymorphism. The alpha-actinins are an ancient family of actin-binding proteins with structural, signalling and metabolic functions. The skeletal muscle alpha-actinins diverged approximately 250-300 million years ago, and ACTN3 has since developed restricted expression in fast muscle fibres. Despite ACTN2 and ACTN3 retaining considerable sequence similarity, it is likely that following duplication there was a divergence in function explaining why alpha-actinin-2 cannot completely compensate for the absence of alpha-actinin-3. This paper focuses on the role of skeletal muscle alpha-actinins, and how possible changes in functions between these duplicates fit in the context of gene duplication paradigms.

Aspects on evolution of fungal beta-lactam biosynthesis gene clusters and recruitment of trans-acting factors

Penicillins and cephalosporins are beta-lactam antibiotics. The formation of hydrophobic penicillins has been reported in fungi only, notably Penicillium chrysogenum and Aspergillus (Emericella) nidulans, whereas the hydrophilic cephalosporins are produced by both fungi, e.g., Acremonium chrysogenum (cephalosporin C), and bacteria. The producing bacteria include Gram-negatives and Gram-positives, e.g., Streptomyces clavuligerus (cephamycin C) and Lysobacter lactamgenus (cephabacins), respectively. The evolutionary origin of beta-lactam biosynthesis genes has been the subject of discussion for many years, and two main hypotheses have been proposed: (i) horizontal gene transfer (HGT) from bacteria to fungi or (ii) vertical decent. There are strong arguments in favour of HGT, e.g., unlike most other fungal genes, beta-lactam biosynthesis genes are clustered and some of these genes lack introns. In contrast to S. clavuligerus, all regulators of fungal beta-lactam biosynthesis genes represent wide-domain regulators that are not part of the gene cluster. If bacterial regulators were co-transferred with the gene cluster from bacteria to fungi, most likely they would have been non-functional in eukaryotes and lost during evolution. Recently, the penicillin biosynthesis gene aatB was discovered, which is not part of the penicillin biosynthesis gene cluster and is even located on a different chromosome. The aatB gene is regulated by the same regulators AnCF and AnBH1 as the penicillin biosynthesis gene aatA (penDE). Data suggest that aatA and aatB are paralogues derived by duplication of a common ancestor gene. This data supports a model in which part of the beta-lactam biosynthesis gene cluster was transferred to some fungi, i.e., the acvA and ipnA gene without a regulatory gene. We propose that during the assembly of aatA and acvA-ipnA into a single gene cluster, recruitment of transcriptional regulators occurred along with acquisition of the duplicated aatA ancestor gene and its cis-acting sites.

Characterization of SCaMC-3-like/slc25a41, a novel calcium-independent mitochondrial ATP-Mg/Pi carrier

The SCaMCs (small calcium-binding mitochondrial carriers) constitute a subfamily of mitochondrial carriers responsible for the ATP-Mg/P(i) exchange with at least three paralogues in vertebrates. SCaMC members are proteins with two functional domains, the C-terminal transporter domain and the N-terminal domain which harbours calcium-binding EF-hands and faces the intermembrane space. In the present study, we have characterized a shortened fourth paralogue, SCaMC-3L (SCaMC-3-like; also named slc25a41), which lacks the calcium-binding N-terminal extension. SCaMC-3L orthologues are found exclusively in mammals, showing approx. 60% identity to the C-terminal half of SCaMC-3, its closest paralogue. In mammalian genomes, SCaMC-3 and SCaMC-3L genes are adjacent on the same chromosome, forming a head-to-tail tandem array, and show identical exon-intron boundaries, indicating that SCaMC-3L could have arisen from an SCaMC-3 ancestor by a partial duplication event which occurred prior to mammalian radiation. Expression and functional data suggest that, following the duplication event, SCaMC-3L has acquired more restrictive functions. Unlike the broadly expressed longer SCaMCs, mouse SCaMC-3L shows a limited expression pattern; it is preferentially expressed in testis and, at lower levels, in brain. SCaMC-3L transport activity was studied in yeast deficient in Sal1p, the calcium-dependent mitochondrial ATP-Mg/P(i) carrier, co-expressing SCaMC-3L and mitochondrial-targeted luciferase, and it was found to perform ATP-Mg/P(i) exchange, in a similar manner to Sal1p or other ATP-Mg/P(i) carriers. However, metabolite transport through SCaMC-3L is calcium-independent, representing a novel mechanism involved in adenine nucleotide transport across the inner mitochondrial membrane, different to ADP/ATP translocases or long SCaMC paralogues.

Identification of the novel penicillin biosynthesis gene aatB of Aspergillus nidulans and its putative evolutionary relationship to this fungal secondary metabolism gene cluster

The final step of penicillin biosynthesis in the filamentous fungus Aspergillus nidulans is catalysed by isopenicillin N acyltransferase encoded by the aatA gene. Because there is no bacterial homologue, its evolutionary origin remained obscure. As shown here,disruption of aatA still enabled penicillin production. Genome mining led to the discovery of the aatB gene(AN6775.3) which has a similar structure and expression pattern as aatA. Disruption of aatB resulted in a reduced penicillin titre. Surface plasmon resonance analysis and Northern blot analysis indicated that the promoters of both aatA and aatB are bound and regulated by the same transcription factors AnCF and AnBH1f. In contrast to aatA, aatB does not encode a peroxisomal targeting signal (PTS1). Overexpression of a mutated aatB(PTS1) gene in an aatA-disruption strain(leading to peroxisomal localization of AatB)increased the penicillin titre more than overexpression of the wild-type aatB. Homologues of aatA are exclusively part of the penicillin biosynthesis gene cluster,whereas aatB homologues also exist in non-producing fungi. Our findings suggest that aatB is a paralogue of aatA. They extend the model of evolution of the penicillin biosynthesis gene cluster by recruitment of a biosynthesis gene and its cis-regulatory sites upon gene duplication.

Asymmetric histone modifications between the original and derived loci of human segmental duplications

A systematic analysis of histone modifications between human segmental duplications shows that two seemingly identical genomic copies have distinct epigenomic properties.

Background

Sequencing and annotation of several mammalian genomes have revealed that segmental duplications are a common architectural feature of primate genomes; in fact, about 5% of the human genome is composed of large blocks of interspersed segmental duplications. These segmental duplications have been implicated in genomic copy-number variation, gene novelty, and various genomic disorders. However, the molecular processes involved in the evolution and regulation of duplicated sequences remain largely unexplored.

Results

In this study, the profile of about 20 histone modifications within human segmental duplications was characterized using high-resolution, genome-wide data derived from a ChIP-Seq study. The analysis demonstrates that derivative loci of segmental duplications often differ significantly from the original with respect to many histone methylations. Further investigation showed that genes are present three times more frequently in the original than in the derivative, whereas pseudogenes exhibit the opposite trend. These asymmetries tend to increase with the age of segmental duplications. The uneven distribution of genes and pseudogenes does not, however, fully account for the asymmetry in the profile of histone modifications.

Conclusion

The first systematic analysis of histone modifications between segmental duplications demonstrates that two seemingly 'identical' genomic copies are distinct in their epigenomic properties. Results here suggest that local chromatin environments may be implicated in the discrimination of derived copies of segmental duplications from their originals, leading to a biased pseudogenization of the new duplicates. The data also indicate that further exploration of the interactions between histone modification and sequence degeneration is necessary in order to understand the divergence of duplicated sequences.

Vitamin B12 partners the CarH repressor to downregulate a photoinducible promoter in Myxococcus xanthus

A light-inducible promoter, P(B), drives expression of the carB operon in Myxococcus xanthus. Repressed by CarA in the dark, P(B) is activated when CarS, produced in the light, sequesters CarA to prevent operator-CarA binding. The MerR-type, N-terminal domain of CarA, which mediates interactions with both operator and CarS, is linked to a C-terminal oligomerization module with a predicted cobalamin-binding motif. Here, we show that although CarA does bind vitamin B12, mutating the motif involved has no effect on its ability to repress P(B). Intriguingly, P(B) could be repressed in the dark even with no CarA, so long as B12 and an intact CarA operator were present. We have discovered that this effect of B12 depends on the gene immediately downstream of carA. Its product, CarH, also consists of a MerR-type, N-terminal domain that specifically recognizes the CarA operator and CarS, linked to a predicted B12-binding C-terminal oligomerization module. The B12-mediated repression of P(B) in the dark is relieved by deleting carH, by mutating the DNA- or B12-binding residues of CarH, or by illumination. Our findings unveil parallel regulatory circuits that control a light-inducible promoter using a transcriptional factor repertoire that includes a paralogous gene pair and vitamin B12.