Modeling the evolution of regulatory elements by simultaneous detection and alignment with phylogenetic pair HMMs.

The computational detection of regulatory elements in DNA is a difficult but important problem impacting our progress in understanding the complex nature of eukaryotic gene regulation. Attempts to utilize cross-species conservation for this task have been hampered both by evolutionary changes of functional sites and poor performance of general-purpose alignment programs when applied to non-coding sequence. We describe a new and flexible framework for modeling binding site evolution in multiple related genomes, based on phylogenetic pair hidden Markov models which explicitly model the gain and loss of binding sites along a phylogeny. We demonstrate the value of this framework for both the alignment of regulatory regions and the inference of precise binding-site locations within those regions. As the underlying formalism is a stochastic, generative model, it can also be used to simulate the evolution of regulatory elements. Our implementation is scalable in terms of numbers of species and sequence lengths and can produce alignments and binding-site predictions with accuracy rivaling or exceeding current systems that specialize in only alignment or only binding-site prediction. We demonstrate the validity and power of various model components on extensive simulations of realistic sequence data and apply a specific model to study Drosophila enhancers in as many as ten related genomes and in the presence of gain and loss of binding sites. Different models and modeling assumptions can be easily specified, thus providing an invaluable tool for the exploration of biological hypotheses that can drive improvements in our understanding of the mechanisms and evolution of gene regulation.

Four-decade responses of soil trace elements to an aggrading old-field forest: B, Mn, Zn, Cu, and Fe.

In the ancient and acidic Ultisol soils of the Southern Piedmont, USA, we studied changes in trace element biogeochemistry over four decades, a period during which formerly cultivated cotton fields were planted with pine seedlings that grew into mature forest stands. In 16 permanent plots, we estimated 40-year accumulations of trace elements in forest biomass and O horizons (between 1957 and 1997), and changes in bioavailable soil fractions indexed by extractions of 0.05 mol/L HCl and 0.2 mol/L acid ammonium oxalate (AAO). Element accumulations in 40-year tree biomass plus O horizons totaled 0.9, 2.9, 4.8, 49.6, and 501.3 kg/ha for Cu, B, Zn, Mn, and Fe, respectively. In response to this forest development, samples of the upper 0.6-m of mineral soil archived in 1962 and 1997 followed one of three patterns. (1) Extractable B and Mn were significantly depleted, by -4.1 and -57.7 kg/ha with AAO, depletions comparable to accumulations in biomass plus O horizons, 2.9 and 49.6 kg/ha, respectively. Tree uptake of B and Mn from mineral soil greatly outpaced resupplies from atmospheric deposition, mineral weathering, and deep-root uptake. (2) Extractable Zn and Cu changed little during forest growth, indicating that nutrient resupplies kept pace with accumulations by the aggrading forest. (3) Oxalate-extractable Fe increased substantially during forest growth, by 275.8 kg/ha, about 10-fold more than accumulations in tree biomass (28.7 kg/ha). The large increases in AAO-extractable Fe in surficial 0.35-m mineral soils were accompanied by substantial accretions of Fe in the forest's O horizon, by 473 kg/ha, amounts that dwarfed inputs via litterfall and canopy throughfall, indicating that forest Fe cycling is qualitatively different from that of other macro- and micronutrients. Bioturbation of surficial forest soil layers cannot account for these fractions and transformations of Fe, and we hypothesize that the secondary forest's large inputs of organic additions over four decades has fundamentally altered soil Fe oxides, potentially altering the bioavailability and retention of macro- and micronutrients, contaminants, and organic matter itself. The wide range of responses among the ecosystem's trace elements illustrates the great dynamics of the soil system over time scales of decades.

The upstream enhancer elements of the G6PC promoter are critical for optimal G6PC expression in murine glycogen storage disease type Ia.

Glycogen storage disease type-Ia (GSD-Ia) patients deficient in glucose-6-phosphatase-α (G6Pase-α or G6PC) manifest impaired glucose homeostasis characterized by fasting hypoglycemia, growth retardation, hepatomegaly, nephromegaly, hyperlipidemia, hyperuricemia, and lactic acidemia. Two efficacious recombinant adeno-associated virus pseudotype 2/8 (rAAV8) vectors expressing human G6Pase-α have been independently developed. One is a single-stranded vector containing a 2864-bp of the G6PC promoter/enhancer (rAAV8-GPE) and the other is a double-stranded vector containing a shorter 382-bp minimal G6PC promoter/enhancer (rAAV8-miGPE). To identify the best construct, a direct comparison of the rAAV8-GPE and the rAAV8-miGPE vectors was initiated to determine the best vector to take forward into clinical trials. We show that the rAAV8-GPE vector directed significantly higher levels of hepatic G6Pase-α expression, achieved greater reduction in hepatic glycogen accumulation, and led to a better toleration of fasting in GSD-Ia mice than the rAAV8-miGPE vector. Our results indicated that additional control elements in the rAAV8-GPE vector outweigh the gains from the double-stranded rAAV8-miGPE transduction efficiency, and that the rAAV8-GPE vector is the current choice for clinical translation in human GSD-Ia.