CpG dinucleotides and the mutation rate of non-CpG DNA

Locus-level L1 DNA methylation profiling reveals the epigenetic and transcriptional interplay between L1s and their integration sites

Long interspersed element 1 (L1) retrotransposons are implicated in human disease and evolution. Their global activity is repressed by DNA methylation, but deciphering the regulation of individual copies has been challenging. Here, we combine short- and long-read sequencing to unveil L1 methylation heterogeneity across cell types, families, and individual loci and elucidate key principles involved. We find that the youngest primate L1 families are specifically hypomethylated in pluripotent stem cells and the placenta but not in most tumors. Locally, intronic L1 methylation is intimately associated with gene transcription. Conversely, the L1 methylation state can propagate to the proximal region up to 300 bp. This phenomenon is accompanied by the binding of specific transcription factors, which drive the expression of L1 and chimeric transcripts. Finally, L1 hypomethylation alone is typically insufficient to trigger L1 expression due to redundant silencing pathways. Our results illuminate the epigenetic and transcriptional interplay between retrotransposons and their host genome.

Mutagenesis of human genomes by endogenous mobile elements on a population scale

Several large-scale Illumina whole-genome sequencing (WGS) and whole-exome sequencing (WES) projects have emerged recently that have provided exceptional opportunities to discover mobile element insertions (MEIs) and study the impact of these MEIs on human genomes. However, these projects also have presented major challenges with respect to the scalability and computational costs associated with performing MEI discovery on tens or even hundreds of thousands of samples. To meet these challenges, we have developed a more efficient and scalable version of our mobile element locator tool (MELT) called CloudMELT. We then used MELT and CloudMELT to perform MEI discovery in 57,919 human genomes and exomes, leading to the discovery of 104,350 nonredundant MEIs. We leveraged this collection (1) to examine potentially active L1 source elements that drive the mobilization of new Alu, L1, and SVA MEIs in humans; (2) to examine the population distributions and subfamilies of these MEIs; and (3) to examine the mutagenesis of GENCODE genes, ENCODE-annotated features, and disease genes by these MEIs. Our study provides new insights on the L1 source elements that drive MEI mutagenesis and brings forth a better understanding of how this mutagenesis impacts human genomes.

NOREC4DNA: using near-optimal rateless erasure codes for DNA storage

BACKGROUND

DNA is a promising storage medium for high-density long-term digital data storage. Since DNA synthesis and sequencing are still relatively expensive tasks, the coding methods used to store digital data in DNA should correct errors and avoid unstable or error-prone DNA sequences. Near-optimal rateless erasure codes, also called fountain codes, are particularly interesting codes to realize high-capacity and low-error DNA storage systems, as shown by Erlich and Zielinski in their approach based on the Luby transform (LT) code. Since LT is the most basic fountain code, there is a large untapped potential for improvement in using near-optimal erasure codes for DNA storage.

RESULTS

We present NOREC4DNA, a software framework to use, test, compare, and improve near-optimal rateless erasure codes (NORECs) for DNA storage systems. These codes can effectively be used to store digital information in DNA and cope with the restrictions of the DNA medium. Additionally, they can adapt to possible variable lengths of DNA strands and have nearly zero overhead. We describe the design and implementation of NOREC4DNA. Furthermore, we present experimental results demonstrating that NOREC4DNA can flexibly be used to evaluate the use of NORECs in DNA storage systems. In particular, we show that NORECs that apparently have not yet been used for DNA storage, such as Raptor and Online codes, can achieve significant improvements over LT codes that were used in previous work. NOREC4DNA is available on https://github.com/umr-ds/NOREC4DNA .

CONCLUSION

NOREC4DNA is a flexible and extensible software framework for using, evaluating, and comparing NORECs for DNA storage systems.

Analysis and annotation of DNA methylation in two nonhuman primate species using the Infinium Human Methylation 450K and EPIC BeadChips

Aim: Nonhuman primates are essential for research on many human diseases. The Infinium Human Methylation450/EPIC BeadChips are popular tools for the study of the methylation state across the human genome at affordable cost. Methods: We performed a precise evaluation and re-annotation of the BeadChip probes for the analysis of genome-wide DNA methylation patterns in rhesus macaques and African green monkeys through in silico analyses combined with functional validation by pyrosequencing. Results: Up to 165,847 of the 450K and 261,545 probes of the EPIC BeadChip can be reliably used. The annotation files are provided in a format compatible with a variety of standard bioinformatic pipelines. Conclusion: Our study will facilitate high-throughput DNA methylation analyses in Macaca mulatta and Chlorocebus sabaeus.

Distinct nucleotide patterns among three subgenomes of bread wheat and their potential origins during domestication after allopolyploidization

BACKGROUND

The speciation and fast global domestication of bread wheat have made a great impact on three subgenomes of bread wheat. DNA base composition is an essential genome feature, which follows the individual-strand base equality rule and [AT]-increase pattern at the genome, chromosome, and polymorphic site levels among thousands of species. Systematic analyses on base compositions of bread wheat and its wild progenitors could facilitate further understanding of the evolutionary pattern of genome/subgenome-wide base composition of allopolyploid species and its potential causes.

RESULTS

Genome/subgenome-wide base-composition patterns were investigated by using the data of polymorphic site in 93 accessions from worldwide populations of bread wheat, its diploid and tetraploid progenitors, and their corresponding reference genome sequences. Individual-strand base equality rule and [AT]-increase pattern remain in recently formed hexaploid species bread wheat at the genome, subgenome, chromosome, and polymorphic site levels. However, D subgenome showed the fastest [AT]-increase across polymorphic site from Aegilops tauschii to bread wheat than that on A and B subgenomes from wild emmer to bread wheat. The fastest [AT]-increase could be detected almost all chromosome windows on D subgenome, suggesting different mechanisms between D and other two subgenomes. Interestingly, the [AT]-increase is mainly contributed by intergenic regions at non-selective sweeps, especially the fastest [AT]-increase of D subgenome. Further transition frequency and sequence context analysis indicated that three subgenomes shared same mutation type, but D subgenome owns the highest mutation rate on high-frequency mutation type. The highest mutation rate on D subgenome was further confirmed by using a bread-wheat-private SNP set. The exploration of loci/genes related to the [AT] value of D subgenome suggests the fastest [AT]-increase of D subgenome could be involved in DNA repair systems distributed on three subgenomes of bread wheat.

CONCLUSIONS

The highest mutation rate is detected on D subgenome of bread wheat during domestication after allopolyploidization, leading to the fastest [AT]-increase pattern of D subgenome. The phenomenon may come from the joint action of multiple repair systems inherited from its wild progenitors.

Cryptic genetic variation enhances primate L1 retrotransposon survival by enlarging the functional coiled coil sequence space of ORF1p

Accounting for continual evolution of deleterious L1 retrotransposon families, which can contain hundreds to thousands of members remains a major issue in mammalian biology. L1 activity generated upwards of 40% of some mammalian genomes, including humans where they remain active, causing genetic defects and rearrangements. L1 encodes a coiled coil-containing protein that is essential for retrotransposition, and the emergence of novel primate L1 families has been correlated with episodes of extensive amino acid substitutions in the coiled coil. These results were interpreted as an adaptive response to maintain L1 activity, however its mechanism remained unknown. Although an adventitious mutation can inactivate coiled coil function, its effect could be buffered by epistatic interactions within the coiled coil, made more likely if the family contains a diverse set of coiled coil sequences-collectively referred to as the coiled coil sequence space. Amino acid substitutions that do not affect coiled coil function (i.e., its phenotype) could be "hidden" from (not subject to) purifying selection. The accumulation of such substitutions, often referred to as cryptic genetic variation, has been documented in various proteins. Here we report that this phenomenon was in effect during the latest episode of primate coiled coil evolution, which occurred 30-10 MYA during the emergence of primate L1Pa7-L1Pa3 families. First, we experimentally demonstrated that while coiled coil function (measured by retrotransposition) can be eliminated by single epistatic mutations, it nonetheless can also withstand extensive amino acid substitutions. Second, principal component and cluster analysis showed that the coiled coil sequence space of each of the L1Pa7-3 families was notably increased by the presence of distinct, coexisting coiled coil sequences. Thus, sampling related networks of functional sequences rather than traversing discrete adaptive states characterized the persistence L1 activity during this evolutionary event.

Variation in base composition underlies functional and evolutionary divergence in non-LTR retrotransposons

Background

Non-LTR retrotransposons often exhibit base composition that is markedly different from the nucleotide content of their host’s gene. For instance, the mammalian L1 element is AT-rich with a strong A bias on the positive strand, which results in a reduced transcription. It is plausible that the A-richness of mammalian L1 is a self-regulatory mechanism reflecting a trade-off between transposition efficiency and the deleterious effect of L1 on its host. We examined if the A-richness of L1 is a general feature of non-LTR retrotransposons or if different clades of elements have evolved different nucleotide content. We also investigated if elements belonging to the same clade evolved towards different base composition in different genomes or if elements from different clades evolved towards similar base composition in the same genome.

Results

We found that non-LTR retrotransposons differ in base composition among clades within the same host but also that elements belonging to the same clade differ in base composition among hosts. We showed that nucleotide content remains constant within the same host over extended period of evolutionary time, despite mutational patterns that should drive nucleotide content away from the observed base composition.

Conclusions

Our results suggest that base composition is evolving under selection and may be reflective of the long-term co-evolution between non-LTR retrotransposons and their host. Finally, the coexistence of elements with drastically different base composition suggests that these elements may be using different strategies to persist and multiply in the genome of their host.

Unparalleled Rapid Evolution of KIR Genes in Rhesus and Cynomolgus Macaque Populations

The killer cell Ig-like receptors (KIR) modulate immune responses through interactions with MHC class I molecules. The KIR region in large cohorts of rhesus and cynomolgus macaque populations were characterized, and the experimental design enabled the definition of a considerable number of alleles (n = 576) and haplotypes, which are highly variable with regard to architecture. Although high levels of polymorphism were recorded, only a few alleles are shared between species and populations. The rapid evolution of allelic polymorphism, accumulated by point mutations, was further confirmed by the emergence of a novel KIR allele in a rhesus macaque family. In addition to allelic variation, abundant orthologous and species-specific KIR genes were identified, the latter of which are frequently generated by fusion events. The concerted action of both genetic mechanisms, in combination with differential selective pressures at the population level, resulted in the unparalleled rapid evolution of the KIR gene region in two closely related macaque species. The variation of the KIR gene repertoire at the species and population level might have an impact on the outcome of preclinical studies with macaque models.

Genome-wide nucleotide patterns and potential mechanisms of genome divergence following domestication in maize and soybean

BACKGROUND

Plant domestication provides a unique model to study genome evolution. Many studies have been conducted to examine genes, genetic diversity, genome structure, and epigenome changes associated with domestication. Interestingly, domesticated accessions have significantly higher [A] and [T] values across genome-wide polymorphic sites than accessions sampled from the corresponding progenitor species. However, the relative contributions of different genomic regions to this genome divergence pattern and underlying mechanisms have not been well characterized.

RESULTS

Here, we investigate the genome-wide base-composition patterns by analyzing millions of SNPs segregating among 100 accessions from a teosinte-maize comparison set and among 302 accessions from a wild-domesticated soybean comparison set. We show that non-genic part of the genome has a greater contribution than genic SNPs to the [AT]-increase observed between wild and domesticated accessions in maize and soybean. The separation between wild and domesticated accessions in [AT] values is significantly enlarged in non-genic and pericentromeric regions. Motif frequency and sequence context analyses show the motifs (PyCG) related to solar-UV signature are enriched in these regions, particularly when they are methylated. Additional analysis using population-private SNPs also implicates the role of these motifs in relatively recent mutations. With base-composition across polymorphic sites as a genome phenotype, genome scans identify a set of putative candidate genes involved in UV damage repair pathways.

CONCLUSIONS

The [AT]-increase is more pronounced in genomic regions that are non-genic, pericentromeric, transposable elements; methylated; and with low recombination. Our findings establish important links among UV radiation, mutation, DNA repair, methylation, and genome evolution.

Identification of recurrent USP48 and BRAF mutations in Cushing's disease

Cushing's disease results from corticotroph adenomas of the pituitary that hypersecrete adrenocorticotropin (ACTH), leading to excess glucocorticoid and hypercortisolism. Mutations of the deubiquitinase gene USP8 occur in 35-62% of corticotroph adenomas. However, the major driver mutations in USP8 wild-type tumors remain elusive. Here, we report recurrent mutations in the deubiquitinase gene USP48 (predominantly encoding p.M415I or p.M415V; 21/91 subjects) and BRAF (encoding p.V600E; 15/91 subjects) in corticotroph adenomas with wild-type USP8. Similar to USP8 mutants, both USP48 and BRAF mutants enhance the promoter activity and transcription of the gene encoding proopiomelanocortin (POMC), which is the precursor of ACTH, providing a potential mechanism for ACTH overproduction in corticotroph adenomas. Moreover, primary corticotroph tumor cells harboring BRAF V600E are sensitive to the BRAF inhibitor vemurafenib. Our study thus contributes to the understanding of the molecular mechanism of the pathogenesis of corticotroph adenoma and informs therapeutic targets for this disease.

APOBEC: From mutator to editor

APOBECs (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like) are a family of cytidine deaminases that prefer single-stranded nucleic acids as substrates. Besides their physiological functions, APOBEC family members have been found to cause hypermutations of cancer genomes, which could be correlated with cancer development and poor prognosis. Recently, APOBEC family members have been combined with the versatile CRISPR/Cas9 system to perform targeted base editing or induce hypermutagenesis. This combination improved the CRISPR/Cas9-mediated gene editing at single-base precision, greatly enhancing its usefulness. Here, we review the physiological functions and structural characteristics of APOBEC family members and their roles as endogenous mutators that contribute to hypermutations during carcinogenesis. We also review the various iterations of the APOBEC-CRISPR/Cas9 gene-editing tools, pointing out their features and limitations as well as the possibilities for future developments.

A defect in myoblast fusion underlies Carey-Fineman-Ziter syndrome

Multinucleate cellular syncytial formation is a hallmark of skeletal muscle differentiation. Myomaker, encoded by Mymk (Tmem8c), is a well-conserved plasma membrane protein required for myoblast fusion to form multinucleated myotubes in mouse, chick, and zebrafish. Here, we report that autosomal recessive mutations in MYMK (OMIM 615345) cause Carey-Fineman-Ziter syndrome in humans (CFZS; OMIM 254940) by reducing but not eliminating MYMK function. We characterize MYMK-CFZS as a congenital myopathy with marked facial weakness and additional clinical and pathologic features that distinguish it from other congenital neuromuscular syndromes. We show that a heterologous cell fusion assay in vitro and allelic complementation experiments in mymk knockdown and mymkinsT/insT zebrafish in vivo can differentiate between MYMK wild type, hypomorphic and null alleles. Collectively, these data establish that MYMK activity is necessary for normal muscle development and maintenance in humans, and expand the spectrum of congenital myopathies to include cell-cell fusion deficits.

Intrinsic Molecular Processes: Impact on Mutagenesis

Mutations provide resources for genome evolution by generating genetic variability. In addition, mutations act as a driving force leading to disease pathogenesis, and thus have important implications for disease diagnosis, prognosis, and treatment. Understanding the mechanisms underlying how mutations occur is therefore of prime importance for elucidating evolutionary and pathogenic processes. Recent genomics studies have revealed that mutations occur non-randomly across the human genome. In particular, the distribution of mutations is highly associated with intrinsic molecular processes including transcription, chromatin organization, DNA replication timing, and DNA repair. Interplay between intrinsic processes and extrinsic mutagenic exposure may thus imprint a characteristic mutational landscape on tumors. We discuss the impact of intrinsic molecular processes on mutation acquisition in cancer.

Simplification, Innateness, and the Absorption of Meaning from Context: How Novelty Arises from Gradual Network Evolution

How does new genetic information arise? Traditional thinking holds that mutation happens by accident and then spreads in the population by either natural selection or random genetic drift. There have been at least two fundamental conceptual problems with imagining an alternative. First, it seemed that the only alternative is a mutation that responds "smartly" to the immediate environment; but in complex multicellulars, it is hard to imagine how this could be implemented. Second, if there were mechanisms of mutation that "knew" what genetic changes would be favored in a given environment, this would have only begged the question of how they acquired that particular knowledge to begin with. This paper offers an alternative that avoids these problems. It holds that mutational mechanisms act on information that is in the genome, based on considerations of simplicity, parsimony, elegance, etc. (which are different than fitness considerations). This simplification process, under the performance pressure exerted by selection, not only leads to the improvement of adaptations but also creates elements that have the capacity to serve in new contexts they were not originally selected for. Novelty, then, arises at the system level from emergent interactions between such elements. Thus, mechanistically driven mutation neither requires Lamarckian transmission nor closes the door on novelty, because the changes it implements interact with one another globally in surprising and beneficial ways. Finally, I argue, for example, that genes used together are fused together; that simplification leads to complexity; and that evolution and learning are conceptually linked.

Single genome retrieval of context-dependent variability in mutation rates for human germline

BACKGROUND

Accurate knowledge of the core components of substitution rates is of vital importance to understand genome evolution and dynamics. By performing a single-genome and direct analysis of 39,894 retrotransposon remnants, we reveal sequence context-dependent germline nucleotide substitution rates for the human genome.

RESULTS

The rates are characterised through rate constants in a time-domain, and are made available through a dedicated program (Trek) and a stand-alone database. Due to the nature of the method design and the imposed stringency criteria, we expect our rate constants to be good estimates for the rates of spontaneous mutations. Benefiting from such data, we study the short-range nucleotide (up to 7-mer) organisation and the germline basal substitution propensity (BSP) profile of the human genome; characterise novel, CpG-independent, substitution prone and resistant motifs; confirm a decreased tendency of moieties with low BSP to undergo somatic mutations in a number of cancer types; and, produce a Trek-based estimate of the overall mutation rate in human.

CONCLUSIONS

The extended set of rate constants we report may enrich our resources and help advance our understanding of genome dynamics and evolution, with possible implications for the role of spontaneous mutations in the emergence of pathological genotypes and neutral evolution of proteomes.

Prediction of the gene expression in normal lung tissue by the gene expression in blood

BACKGROUND

Comparative analysis of gene expression in human tissues is important for understanding the molecular mechanisms underlying tissue-specific control of gene expression. It can also open an avenue for using gene expression in blood (which is the most easily accessible human tissue) to predict gene expression in other (less accessible) tissues, which would facilitate the development of novel gene expression based models for assessing disease risk and progression. Until recently, direct comparative analysis across different tissues was not possible due to the scarcity of paired tissue samples from the same individuals.

METHODS

In this study we used paired whole blood/lung gene expression data from the Genotype-Tissue Expression (GTEx) project. We built a generalized linear regression model for each gene using gene expression in lung as the outcome and gene expression in blood, age and gender as predictors.

RESULTS

For ~18 % of the genes, gene expression in blood was a significant predictor of gene expression in lung. We found that the number of single nucleotide polymorphisms (SNPs) influencing expression of a given gene in either blood or lung, also known as the number of quantitative trait loci (eQTLs), was positively associated with efficacy of blood-based prediction of that gene's expression in lung. This association was strongest for shared eQTLs: those influencing gene expression in both blood and lung.

CONCLUSIONS

In conclusion, for a considerable number of human genes, their expression levels in lung can be predicted using observable gene expression in blood. An abundance of shared eQTLs may explain the strong blood/lung correlations in the gene expression.

DNA Demethylation Dynamics in the Human Prenatal Germline

Global DNA demethylation in humans is a fundamental process that occurs in pre-implantation embryos and reversion to naive ground state pluripotent stem cells (PSCs). However, the extent of DNA methylation reprogramming in human germline cells is unknown. Here, we performed whole-genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq) of human prenatal germline cells from 53 to 137 days of development. We discovered that the transcriptome and methylome of human germline is distinct from both human PSCs and the inner cell mass (ICM) of human blastocysts. Using this resource to monitor the outcome of global DNA demethylation with reversion of primed PSCs to the naive ground state, we uncovered hotspots of ultralow methylation at transposons that are protected from demethylation in the germline and ICM. Taken together, the human germline serves as a valuable in vivo tool for monitoring the epigenome of cells that have emerged from a global DNA demethylation event.

CXCL8 histone H3 acetylation is dysfunctional in airway smooth muscle in asthma: regulation by BET

Asthma is characterized by airway inflammation and remodeling and CXCL8 is a CXC chemokine that drives steroid-resistant neutrophilic airway inflammation. We have shown that airway smooth muscle (ASM) cells isolated from asthmatic individuals secrete more CXCL8 than cells from nonasthmatic individuals. Here we investigated chromatin modifications at the CXCL8 promoter in ASM cells from nonasthmatic and asthmatic donors to further understand how CXCL8 is dysregulated in asthma. ASM cells from asthmatic donors had increased histone H3 acetylation, specifically histone H3K18 acetylation, and increased binding of histone acetyltransferase p300 compared with nonasthmatic donors but no differences in CXCL8 DNA methylation. The acetylation reader proteins Brd3 and Brd4 were bound to the CXCL8 promoter and Brd inhibitors inhibited CXCL8 secretion from ASM cells by disrupting Brd4 and RNA polymerase II binding to the CXCL8 promoter. Our results show a novel dysregulation of CXCL8 transcriptional regulation in asthma characterized by a promoter complex that is abnormal in ASM cells isolated from asthmatic donors and can be modulated by Brd inhibitors. Brd inhibitors may provide a new therapeutic strategy for steroid-resistant inflammation.

Breaking bad: The mutagenic effect of DNA repair

Species survival depends on the faithful replication of genetic information, which is continually monitored and maintained by DNA repair pathways that correct replication errors and the thousands of lesions that arise daily from the inherent chemical lability of DNA and the effects of genotoxic agents. Nonetheless, neutrally evolving DNA (not under purifying selection) accumulates base substitutions with time (the neutral mutation rate). Thus, repair processes are not 100% efficient. The neutral mutation rate varies both between and within chromosomes. For example it is 10-50 fold higher at CpGs than at non-CpG positions. Interestingly, the neutral mutation rate at non-CpG sites is positively correlated with CpG content. Although the basis of this correlation was not immediately apparent, some bioinformatic results were consistent with the induction of non-CpG mutations by DNA repair at flanking CpG sites. Recent studies with a model system showed that in vivo repair of preformed lesions (mismatches, abasic sites, single stranded nicks) can in fact induce mutations in flanking DNA. Mismatch repair (MMR) is an essential component for repair-induced mutations, which can occur as distant as 5 kb from the introduced lesions. Most, but not all, mutations involved the C of TpCpN (G of NpGpA) which is the target sequence of the C-preferring single-stranded DNA specific APOBEC deaminases. APOBEC-mediated mutations are not limited to our model system: Recent studies by others showed that some tumors harbor mutations with the same signature, as can intermediates in RNA-guided endonuclease-mediated genome editing. APOBEC deaminases participate in normal physiological functions such as generating mutations that inactivate viruses or endogenous retrotransposons, or that enhance immunoglobulin diversity in B cells. The recruitment of normally physiological error-prone processes during DNA repair would have important implications for disease, aging and evolution. This perspective briefly reviews both the bioinformatic and biochemical literature relevant to repair-induced mutagenesis and discusses future directions required to understand the mechanistic basis of this process.

Identification of human-specific AluS elements through comparative genomics

Mobile elements are responsible for ~45% of the human genome. Among them is the Alu element, accounting for 10% of the human genome (>1.1million copies). Several studies of Alu elements have reported that they are frequently involved in human genetic diseases and genomic rearrangements. In this study, we investigated the AluS subfamily, which is a relatively old Alu subfamily and has the highest copy number in primate genomes. Previously, a set of 263 human-specific AluS insertions was identified in the human genome. To validate these, we compared each of the human-specific AluS loci with its pre-insertion site in other primate genomes, including chimpanzee, gorilla, and orangutan. We obtained 24 putative human-specific AluS candidates via the in silico analysis and manual inspection, and then tried to verify them using PCR amplification and DNA sequencing. Through the PCR product sequencing, we were able to detect two instances of near-parallel Alu insertions in nearby sites that led to computational false negatives. Finally, we computationally and experimentally verified 23 human-specific AluS elements. We reported three alternative Alu insertion events, which are accompanied by filler DNA and/or Alu retrotransposition mediated-deletion. Bisulfite sequencing was carried out to examine DNA methylation levels of human-specific AluS elements. The results showed that fixed AluS elements are hypermethylated compared with polymorphic elements, indicating a possible relation between DNA methylation and Alu fixation in the human genome.

Early developmental conditioning of later health and disease: physiology or pathophysiology?

Extensive experimental animal studies and epidemiological observations have shown that environmental influences during early development affect the risk of later pathophysiological processes associated with chronic, especially noncommunicable, disease (NCD). This field is recognized as the developmental origins of health and disease (DOHaD). We discuss the extent to which DOHaD represents the result of the physiological processes of developmental plasticity, which may have potential adverse consequences in terms of NCD risk later, or whether it is the manifestation of pathophysiological processes acting in early life but only becoming apparent as disease later. We argue that the evidence suggests the former, through the operation of conditioning processes induced across the normal range of developmental environments, and we summarize current knowledge of the physiological processes involved. The adaptive pathway to later risk accords with current concepts in evolutionary developmental biology, especially those concerning parental effects. Outside the normal range, effects on development can result in nonadaptive processes, and we review their underlying mechanisms and consequences. New concepts concerning the underlying epigenetic and other mechanisms involved in both disruptive and nondisruptive pathways to disease are reviewed, including the evidence for transgenerational passage of risk from both maternal and paternal lines. These concepts have wider implications for understanding the causes and possible prevention of NCDs such as type 2 diabetes and cardiovascular disease, for broader social policy and for the increasing attention paid in public health to the lifecourse approach to NCD prevention.

Interaction-based evolution: how natural selection and nonrandom mutation work together

BACKGROUND

The modern evolutionary synthesis leaves unresolved some of the most fundamental, long-standing questions in evolutionary biology: What is the role of sex in evolution? How does complex adaptation evolve? How can selection operate effectively on genetic interactions? More recently, the molecular biology and genomics revolutions have raised a host of critical new questions, through empirical findings that the modern synthesis fails to explain: for example, the discovery of de novo genes; the immense constructive role of transposable elements in evolution; genetic variance and biochemical activity that go far beyond what traditional natural selection can maintain; perplexing cases of molecular parallelism; and more.

PRESENTATION OF THE HYPOTHESIS

Here I address these questions from a unified perspective, by means of a new mechanistic view of evolution that offers a novel connection between selection on the phenotype and genetic evolutionary change (while relying, like the traditional theory, on natural selection as the only source of feedback on the fit between an organism and its environment). I hypothesize that the mutation that is of relevance for the evolution of complex adaptation-while not Lamarckian, or "directed" to increase fitness-is not random, but is instead the outcome of a complex and continually evolving biological process that combines information from multiple loci into one. This allows selection on a fleeting combination of interacting alleles at different loci to have a hereditary effect according to the combination's fitness.

TESTING AND IMPLICATIONS OF THE HYPOTHESIS

This proposed mechanism addresses the problem of how beneficial genetic interactions can evolve under selection, and also offers an intuitive explanation for the role of sex in evolution, which focuses on sex as the generator of genetic combinations. Importantly, it also implies that genetic variation that has appeared neutral through the lens of traditional theory can actually experience selection on interactions and thus has a much greater adaptive potential than previously considered. Empirical evidence for the proposed mechanism from both molecular evolution and evolution at the organismal level is discussed, and multiple predictions are offered by which it may be tested.

REVIEWERS

This article was reviewed by Nigel Goldenfeld (nominated by Eugene V. Koonin), Jürgen Brosius and W. Ford Doolittle.

Comparative analysis of context-dependent mutagenesis using human and mouse models

Substitution rates strongly depend on their nucleotide context. One of the most studied examples is the excess of C > T mutations in the CG context in various groups of organisms, including vertebrates. Studies on the molecular mechanisms underlying this mutation regularity have provided insights into evolution, mutagenesis, and cancer development. Recently several other hypermutable motifs were identified in the human genome. There is an increased frequency of T > C mutations in the second position of the words ATTG and ATAG and an increased frequency of A > C mutations in the first position of the word ACAA. For a better understanding of evolution, it is of interest whether these mutation regularities are human specific or present in other vertebrates, as their presence might affect the validity of currently used substitution models and molecular clocks. A comprehensive analysis of mutagenesis in 4 bp mutation contexts requires a vast amount of mutation data. Such data may be derived from the comparisons of individual genomes or from single nucleotide polymorphism (SNP) databases. Using this approach, we performed a systematical comparison of mutation regularities within 2-4 bp contexts in Mus musculus and Homo sapiens and uncovered that even closely related organisms may have notable differences in context-dependent mutation regularities.

Segmenting the human genome based on states of neutral genetic divergence

Many studies have demonstrated that divergence levels generated by different mutation types vary and covary across the human genome. To improve our still-incomplete understanding of the mechanistic basis of this phenomenon, we analyze several mutation types simultaneously, anchoring their variation to specific regions of the genome. Using hidden Markov models on insertion, deletion, nucleotide substitution, and microsatellite divergence estimates inferred from human-orangutan alignments of neutrally evolving genomic sequences, we segment the human genome into regions corresponding to different divergence states--each uniquely characterized by specific combinations of divergence levels. We then parsed the mutagenic contributions of various biochemical processes associating divergence states with a broad range of genomic landscape features. We find that high divergence states inhabit guanine- and cytosine (GC)-rich, highly recombining subtelomeric regions; low divergence states cover inner parts of autosomes; chromosome X forms its own state with lowest divergence; and a state of elevated microsatellite mutability is interspersed across the genome. These general trends are mirrored in human diversity data from the 1000 Genomes Project, and departures from them highlight the evolutionary history of primate chromosomes. We also find that genes and noncoding functional marks [annotations from the Encyclopedia of DNA Elements (ENCODE)] are concentrated in high divergence states. Our results provide a powerful tool for biomedical data analysis: segmentations can be used to screen personal genome variants--including those associated with cancer and other diseases--and to improve computational predictions of noncoding functional elements.

Epigenetics and miRNA emerge as key regulators of smooth muscle cell phenotype and function

Regulation of phenotypic plasticity in smooth muscle requires an understanding of the mechanisms regulating phenotype-specific genes and the processes dysregulated during pathogenesis. Decades of study in airway smooth muscle has provided extensive knowledge of the gene expression patterns and signaling pathways necessary to maintain and alter smooth muscle cell phenotype. With this solid foundation, the importance and complexity of inheritable epigenetic modifications and mechanisms silencing gene expression have now emerged as fundamental components regulating aspects of inflammation, proliferation and remodeling.

Point mutation instability (PIN) mutator phenotype as model for true back mutations seen in hereditary tyrosinemia type 1 - a hypothesis

Hereditary tyrosinemia type 1 (HT1) is an autosomal recessive disorder affecting fumarylacetoacetate hydrolase (FAH), the last enzyme in the tyrosine catabolism pathway. The liver mosaicism observed in HT1 patients is due to the reversion to the wild type of one allele of the original point mutation in fah. It is generally accepted that these reversions are true back mutations; however, the mechanism is still unresolved. Previous reports excluded intragenic recombination, mitotic recombination, or homologous recombination with a pseudogene as possible mechanisms of mutation reversion in HT1. Sequence analysis did not reveal DNA motifs, tandem repeats or other sequence peculiarities that may be involved in mutation reversion. We propose the hypothesis that a point mutation instability mutator (PIN) phenotype brought about by the sustained stress environment created by the accumulating metabolites in the cell is the driver of the true back mutations in HT1. The metabolites accumulating in HT1 create a sustained stress environment by activating the extracellular signal-regulated kinase (ERK) and AKT survival pathways, inducing aberrant mitosis and development of death resistant cells, depleting glutathione, and impairing DNA ligase IV and possibly DNA polymerases δ and ε. This continual production of proliferative and stress-related survival signals in the cellular environment coupled with the mutagenicity of FAA, may instigate a mutator phenotype and could end in tumorigenesis and/or mutation reversion. The establishment of a PIN-mutator phenotype therefore not only seems to be a possible mechanism underlying the true back mutations, but also contributes to explaining the clinical heterogeneity seen in hereditary tyrosinemia type 1.

The impact of the organism on its descendants

Historically, evolutionary biologists have taken the view that an understanding of development is irrelevant to theories of evolution. However, the integration of several disciplines in recent years suggests that this position is wrong. The capacity of the organism to adapt to challenges from the environment can set up conditions that affect the subsequent evolution of its descendants. Moreover, molecular events arising from epigenetic processes can be transmitted from one generation to the next and influence genetic mutation. This in turn can facilitate evolution in the conditions in which epigenetic change was first initiated.

Variation in the mutation rate across mammalian genomes

It has been known for many years that the mutation rate varies across the genome. However, only with the advent of large genomic data sets is the full extent of this variation becoming apparent. The mutation rate varies over many different scales, from adjacent sites to whole chromosomes, with the strongest variation seen at the smallest scales. Some of these patterns have clear mechanistic bases, but much of the rate variation remains unexplained, and some of it is deeply perplexing. Variation in the mutation rate has important implications in evolutionary biology and underexplored implications for our understanding of hereditary disease and cancer.

On the sequence-directed nature of human gene mutation: the role of genomic architecture and the local DNA sequence environment in mediating gene mutations underlying human inherited disease

Different types of human gene mutation may vary in size, from structural variants (SVs) to single base-pair substitutions, but what they all have in common is that their nature, size and location are often determined either by specific characteristics of the local DNA sequence environment or by higher order features of the genomic architecture. The human genome is now recognized to contain "pervasive architectural flaws" in that certain DNA sequences are inherently mutation prone by virtue of their base composition, sequence repetitivity and/or epigenetic modification. Here, we explore how the nature, location and frequency of different types of mutation causing inherited disease are shaped in large part, and often in remarkably predictable ways, by the local DNA sequence environment. The mutability of a given gene or genomic region may also be influenced indirectly by a variety of noncanonical (non-B) secondary structures whose formation is facilitated by the underlying DNA sequence. Since these non-B DNA structures can interfere with subsequent DNA replication and repair and may serve to increase mutation frequencies in generalized fashion (i.e., both in the context of subtle mutations and SVs), they have the potential to serve as a unifying concept in studies of mutational mechanisms underlying human inherited disease.

Placenta-specific expression of the interleukin-2 (IL-2) receptor β subunit from an endogenous retroviral promoter

The long terminal repeat (LTR) sequences of endogenous retroviruses and retroelements contain promoter elements and are known to form chimeric transcripts with nearby cellular genes. Here we show that an LTR of the THE1D retroelement family has been domesticated as an alternative promoter of human IL2RB, the gene encoding the β subunit of the IL-2 receptor. The LTR promoter confers expression specifically in the placental trophoblast as opposed to its native transcription in the hematopoietic system. Rather than sequence-specific determinants, DNA methylation was found to regulate transcription initiation and splicing efficiency in a tissue-specific manner. Furthermore, we detected the cytoplasmic signaling domain of the IL-2Rβ protein in the placenta, suggesting that IL-2Rβ undergoes preferential proteolytic cleavage in this tissue. These findings implicate novel functions for this cytokine receptor subunit in the villous trophoblast and reveal an intriguing example of ancient LTR exaptation to drive tissue-specific gene expression.

A codon substitution model that incorporates the effect of the GC contents, the gene density and the density of CpG islands of human chromosomes

Background

Developing a model for codon substitutions is essential for the analyses of protein sequences. Recent studies on the mutation rates in the non-coding regions have shown that CpG mutation rates in the human genome are negatively correlated to the local GC content and to the densities of functional elements. This study aimed at understanding the effect of genomic features, namely, GC content, gene density, and frequency of CpG islands, on the rates of codon substitution in human chromosomes.

Results

Codon substitution rates of CpG to TpG mutations, TpG to CpG mutations, and non-CpG transitions and transversions in humans were estimated by comparing the coding regions of thousands of human and chimpanzee genes and inferring their ancestral sequences by using macaque genes as the outgroup. Since the genomic features are depending on each other, partial regression coefficients of these features were obtained.

Conclusion

The substitution rates of codons depend on gene densities of the chromosomes. Transcription-associated mutation is one such pressure. On the basis of these results, a model of codon substitutions that incorporates the effect of genomic features on codon substitution in human chromosomes was developed.

Substitution patterns are under different influences in primates and rodents

There are large-scale variations of the GC-content along mammalian chromosomes that have been called isochore structures. Primates and rodents have different isochore structures, which suggests that these lineages exhibit different modes of GC-content evolution. It has been shown that, in the human lineage, GC-biased gene conversion (gBGC), a neutral process associated with meiotic recombination, acts on GC-content evolution by influencing A or T to G or C substitution rates. We computed genome-wide substitution patterns in the mouse lineage from multiple alignments and compared them with substitution patterns in the human lineage. We found that in the mouse lineage, gBGC is active but weaker than in the human lineage and that male-specific recombination better predicts GC-content evolution than female-specific recombination. Furthermore, we were able to show that G or C to A or T substitution rates are predicted by a combination of different factors in both lineages. A or T to G or C substitution rates are most strongly predicted by meiotic recombination in the human lineage but by CpG odds ratio (the observed CpG frequency normalized by the expected CpG frequency) in the mouse lineage, suggesting that substitution patterns are under different influences in primates and rodents.

Relationship between amino acid composition and gene expression in the mouse genome

BACKGROUND

Codon bias is a phenomenon that refers to the differences in the frequencies of synonymous codons among different genes. In many organisms, natural selection is considered to be a cause of codon bias because codon usage in highly expressed genes is biased toward optimal codons. Methods have previously been developed to predict the expression level of genes from their nucleotide sequences, which is based on the observation that synonymous codon usage shows an overall bias toward a few codons called major codons. However, the relationship between codon bias and gene expression level, as proposed by the translation-selection model, is less evident in mammals.

FINDINGS

We investigated the correlations between the expression levels of 1,182 mouse genes and amino acid composition, as well as between gene expression and codon preference. We found that a weak but significant correlation exists between gene expression levels and amino acid composition in mouse. In total, less than 10% of variation of expression levels is explained by amino acid components. We found the effect of codon preference on gene expression was weaker than the effect of amino acid composition, because no significant correlations were observed with respect to codon preference.

CONCLUSION

These results suggest that it is difficult to predict expression level from amino acid components or from codon bias in mouse.

Integrated genetic and epigenetic analysis identifies haplotype-specific methylation in the FTO type 2 diabetes and obesity susceptibility locus

Recent multi-dimensional approaches to the study of complex disease have revealed powerful insights into how genetic and epigenetic factors may underlie their aetiopathogenesis. We examined genotype-epigenotype interactions in the context of Type 2 Diabetes (T2D), focussing on known regions of genomic susceptibility. We assayed DNA methylation in 60 females, stratified according to disease susceptibility haplotype using previously identified association loci. CpG methylation was assessed using methylated DNA immunoprecipitation on a targeted array (MeDIP-chip) and absolute methylation values were estimated using a Bayesian algorithm (BATMAN). Absolute methylation levels were quantified across LD blocks, and we identified increased DNA methylation on the FTO obesity susceptibility haplotype, tagged by the rs8050136 risk allele A (p = 9.40×10⁻⁴, permutation p = 1.0×10⁻³). Further analysis across the 46 kb LD block using sliding windows localised the most significant difference to be within a 7.7 kb region (p = 1.13×10⁻⁷). Sequence level analysis, followed by pyrosequencing validation, revealed that the methylation difference was driven by the co-ordinated phase of CpG-creating SNPs across the risk haplotype. This 7.7 kb region of haplotype-specific methylation (HSM), encapsulates a Highly Conserved Non-Coding Element (HCNE) that has previously been validated as a long-range enhancer, supported by the histone H3K4me1 enhancer signature. This study demonstrates that integration of Genome-Wide Association (GWA) SNP and epigenomic DNA methylation data can identify potential novel genotype-epigenotype interactions within disease-associated loci, thus providing a novel route to aid unravelling common complex diseases.

Origin and evolution of a placental-specific microRNA family in the human genome

Background

MicroRNAs (miRNAs) are a class of short regulatory RNAs encoded in the genome of DNA viruses, some single cell organisms, plants and animals. With the rapid development of technology, more and more miRNAs are being discovered. However, the origin and evolution of most miRNAs remain obscure. Here we report the origin and evolution dynamics of a human miRNA family.

Results

We have shown that all members of the miR-1302 family are derived from MER53 elements. Although the conservation scores of the MER53-derived pre-miRNA sequences are low, we have identified 36 potential paralogs of MER53-derived miR-1302 genes in the human genome and 58 potential orthologs of the human miR-1302 family in placental mammals. We suggest that in placental species, this miRNA family has evolved following the birth-and-death model of evolution. Three possible mechanisms that can mediate miRNA duplication in evolutionary history have been proposed: the transposition of the MER53 element, segmental duplications and Alu-mediated recombination. Finally, we have found that the target genes of miR-1302 are over-represented in transportation, localization, and system development processes and in the positive regulation of cellular processes. Many of them are predicted to function in binding and transcription regulation.

Conclusions

The members of miR-1302 family that are derived from MER53 elements are placental-specific miRNAs. They emerged at the early stage of the recent 180 million years since eutherian mammals diverged from marsupials. Under the birth-and-death model, the miR-1302 genes have experienced a complex expansion with some members evolving by segmental duplications and some by Alu-mediated recombination events.

Sequence composition similarities with the 7SL RNA are highly predictive of functional genomic features

Transposable elements derived from the 7SL RNA gene, such as Alu elements in primates, have had remarkable success in several mammalian lineages. The results presented here show a broad spectrum of functions for genomic segments that display sequence composition similarities with the 7SL RNA gene. Using thoroughly documented loci, we report that DNaseI-hypersensitive sites can be singled out in large genomic sequences by an assessment of sequence composition similarities with the 7SL RNA gene. We apply a root word frequency approach to illustrate a distinctive relationship between the sequence of the 7SL RNA gene and several classes of functional genomic features that are not presumed to be of transposable origin. Transposable elements that show noticeable similarities with the 7SL sequence include Alu sequences, as expected, but also long terminal repeats and the 5′-untranslated regions of long interspersed repetitive elements. In sequences masked for repeated elements, we find, when using the 7SL RNA gene as query sequence, distinctive similarities with promoters, exons and distal gene regulatory regions. The latter being the most notoriously difficult to detect, this approach may be useful for finding genomic segments that have regulatory functions and that may have escaped detection by existing methods.

Chromosomal distribution of disease genes in the human genome

Genes are nonrandomly distributed in the human genome, both within and between chromosomes. Thus, genes of similar function and common evolutionary origin are often clustered, as are genes with similar expression profiles. We now report that the >2400 genes known to underlie human monogenic inherited disease are non-randomly distributed in the genome over and above the general nonrandomness evident in the distribution of human genes. Further, a subset of 315 inherited disease genes subject to gross deletion was found to exhibit a degree of clustering that was twice that manifested by disease genes in general. The clustering of human disease genes is likely to have important implications for understanding the genotype-phenotype relationship in contiguous gene syndromes as well as those conditions characterized by multigene deletions or complex chromosomal rearrangements.

The impact of CpG island on defining transcriptional activation of the mouse L1 retrotransposable elements

Background

L1 retrotransposable elements are potent insertional mutagens responsible for the generation of genomic variation and diversification of mammalian genomes, but reliable estimates of the numbers of actively transposing L1 elements are mostly nonexistent. While the human and mouse genomes contain comparable numbers of L1 elements, several phylogenetic and L1Xplore analyses in the mouse genome suggest that 1,500–3,000 active L1 elements currently exist and that they are still expanding in the genome. Conversely, the human genome contains only 150 active L1 elements. In addition, there is a discrepancy among the nature and number of mouse L1 elements in L1Xplore and the mouse genome browser at the UCSC and in the literature. To date, the reason why a high copy number of active L1 elements exist in the mouse genome but not in the human genome is unknown, as are the potential mechanisms that are responsible for transcriptional activation of mouse L1 elements.

Methodology/Principal Findings

We analyzed the promoter sequences of the 1,501 potentially active mouse L1 elements retrieved from the GenBank and L1Xplore databases and evaluated their transcription factors binding sites and CpG content. To this end, we found that a substantial number of mouse L1 elements contain altered transcription factor YY1 binding sites on their promoter sequences that are required for transcriptional initiation, suggesting that only a half of L1 elements are capable of being transcriptionally active. Furthermore, we present experimental evidence that previously unreported CpG islands exist in the promoters of the most active T_F family of mouse L1 elements. The presence of sequence variations and polymorphisms in CpG islands of L1 promoters that arise from transition mutations indicates that CpG methylation could play a significant role in determining the activity of L1 elements in the mouse genome.

Conclusions

A comprehensive analysis of mouse L1 promoters suggests that the number of transcriptionally active elements is significantly lower than the total number of full-length copies from the three active mouse L1 families. Like human L1 elements, the CpG islands and potentially the transcription factor YY1 binding sites are likely to be required for transcriptional initiation of mouse L1 elements.

The mutational spectrum of non-CpG DNA varies with CpG content

The accumulation of base substitutions (mutations) not subject to natural selection is the neutral mutation rate. Because this rate reflects the in vivo processes involved in maintaining the integrity of genetic information, the factors that affect the neutral mutation rate are of considerable interest. Mammals exhibit two dramatically different neutral mutation rates: the CpG mutation rate, wherein the C of most CpGs (i.e., methyl-CpG) mutate at 10-50 times that of C in any other context or of any other base. The latter mutations constitute the non-CpG rate. The high CpG rate results from the spontaneous deamination of methyl-C to T and incomplete restoration of the ensuing T:G mismatches to C:Gs. Here, we determined the neutral non-CpG mutation rate as a function of CpG content by comparing sequence divergence of thousands of pairs of neutrally evolving chimpanzee and human orthologs that differ primarily in CpG content. Both the mutation rate and the mutational spectrum (transition/transversion ratio) of non-CpG residues change in parallel as sigmoidal (logistic) functions of CpG content. As different mechanisms generate transitions and transversions, these results indicate that both mutation rate and mutational processes are contingent on the local CpG content. We consider several possible mechanisms that might explain how CpG exerts these effects.

The (r)evolution of SINE versus LINE distributions in primate genomes: sex chromosomes are important

The densities of transposable elements (TEs) in the human genome display substantial variation both within individual chromosomes and among chromosome types (autosomes and the two sex chromosomes). Finding an explanation for this variability has been challenging, especially in light of genome landscapes unique to the sex chromosomes. Here, using a multiple regression framework, we investigate primate Alu and L1 densities shaped by regional genome features and location on a particular chromosome type. As a result of our analysis, first, we build statistical models explaining up to 79% and 44% of variation in Alu and L1 element density, respectively. Second, we analyze sex chromosome versus autosome TE densities corrected for regional genomic effects. We discover that sex-chromosome bias in Alu and L1 distributions not only persists after accounting for these effects, but even presents differences in patterns, confirming preferential Alu integration in the male germline, yet likely integration of L1s in both male and female germlines or in early embryogenesis. Additionally, our models reveal that local base composition (measured by GC content and density of L1 target sites) and natural selection (inferred via density of most conserved elements) are significant to predicting densities of L1s. Interestingly, measurements of local double-stranded breaks (a 13-mer associated with genome instability) strongly correlate with densities of Alu elements; little evidence was found for the role of recombination-driven deletion in driving TE distributions over evolutionary time. Thus, Alu and L1 densities have been influenced by the combination of distinct local genome landscapes and the unique evolutionary dynamics of sex chromosomes.

Characteristics of the CArG-SRF binding context in mammalian genomes

Serum response factor (SRF), a member of the MADS family, binds a 10-bp cis element known as the CArG box. However, despite our extensive knowledge of SRF and the CArG box, limited information is available on the CArG-SRF binding context or how CArG flanking sequences are defined. We used statistical tests and computer simulation to find characteristics of the CArG-SRF binding context. Based on the combination of published literature and a search of DBTSS, 150 and 136 functional CArG boxes together with 10 bp flanking on each side were found in mouse and human genomes, respectively. By statistical analysis of the 30 positions we found some new conserved positions of interest (P < 0.005) such as -15, -8, and +8, in addition to the ten highly conserved positions of the CArG box. Intriguingly, studies comparing the flanking positions between consensus CArG boxes and CArG-like boxes showed that there are more conserved positions in the latter. Moreover, CpG content within the CArG-SRF binding context is much higher than that within introns. Collectively, these results suggest that there are some special pre-existing features in the flanking sequences of functional CArG boxes probably contributing to SRF selectively recognizing and binding to the functional CArG from millions of functionless CArG boxes in mammalian genomes. This is a significant step toward understanding the mechanism of transcriptional regulation of SRF-dependent genes.

Understanding what determines the frequency and pattern of human germline mutations

Surprising findings about human germline mutation have come from applying new technologies to detect rare mutations in germline DNA, from analysing DNA sequence divergence between humans and closely related species, and from investigating human polymorphic variation. In this Review we discuss how these approaches affect our current understanding of the roles of sex, age, mutation hot spots, germline selection and genomic factors in determining human nucleotide substitution mutation patterns and frequencies. To enhance our understanding of mutation and disease, more extensive molecular data on the human germ line with regard to mutation origin, DNA repair, epigenetic status and the effect of newly arisen mutations on gamete development are needed.

Mutation patterns in the human genome: more variable than expected

Why are some genomic positions more mutable than others? The identification of cryptic mutation hotspots in the human genome indicates that the determinants of mutation rates are more complex than anticipated.