Histone modifications within the human X centromere region.

Human centromeres are multi-megabase regions of highly ordered arrays of alpha satellite DNA that are separated from chromosome arms by unordered alpha satellite monomers and other repetitive elements. Complexities in assembling such large repetitive regions have limited detailed studies of centromeric chromatin organization. However, a genomic map of the human X centromere has provided new opportunities to explore genomic architecture of a complex locus. We used ChIP to examine the distribution of modified histones within centromere regions of multiple X chromosomes. Methylation of H3 at lysine 4 coincided with DXZ1 higher order alpha satellite, the site of CENP-A localization. Heterochromatic histone modifications were distributed across the 400-500 kb pericentromeric regions. The large arrays of alpha satellite and gamma satellite DNA were enriched for both euchromatic and heterochromatic modifications, implying that some pericentromeric repeats have multiple chromatin characteristics. Partial truncation of the X centromere resulted in reduction in the size of the CENP-A/Cenp-A domain and increased heterochromatic modifications in the flanking pericentromere. Although the deletion removed approximately 1/3 of centromeric DNA, the ratio of CENP-A to alpha satellite array size was maintained in the same proportion, suggesting that a limited, but defined linear region of the centromeric DNA is necessary for kinetochore assembly. Our results indicate that the human X centromere contains multiple types of chromatin, is organized similarly to smaller eukaryotic centromeres, and responds to structural changes by expanding or contracting domains.

Comparative genomic data of the Avian Phylogenomics Project.

BACKGROUND: The evolutionary relationships of modern birds are among the most challenging to understand in systematic biology and have been debated for centuries. To address this challenge, we assembled or collected the genomes of 48 avian species spanning most orders of birds, including all Neognathae and two of the five Palaeognathae orders, and used the genomes to construct a genome-scale avian phylogenetic tree and perform comparative genomics analyses (Jarvis et al. in press; Zhang et al. in press). Here we release assemblies and datasets associated with the comparative genome analyses, which include 38 newly sequenced avian genomes plus previously released or simultaneously released genomes of Chicken, Zebra finch, Turkey, Pigeon, Peregrine falcon, Duck, Budgerigar, Adelie penguin, Emperor penguin and the Medium Ground Finch. We hope that this resource will serve future efforts in phylogenomics and comparative genomics. FINDINGS: The 38 bird genomes were sequenced using the Illumina HiSeq 2000 platform and assembled using a whole genome shotgun strategy. The 48 genomes were categorized into two groups according to the N50 scaffold size of the assemblies: a high depth group comprising 23 species sequenced at high coverage (>50X) with multiple insert size libraries resulting in N50 scaffold sizes greater than 1 Mb (except the White-throated Tinamou and Bald Eagle); and a low depth group comprising 25 species sequenced at a low coverage (~30X) with two insert size libraries resulting in an average N50 scaffold size of about 50 kb. Repetitive elements comprised 4%-22% of the bird genomes. The assembled scaffolds allowed the homology-based annotation of 13,000 ~ 17000 protein coding genes in each avian genome relative to chicken, zebra finch and human, as well as comparative and sequence conservation analyses. CONCLUSIONS: Here we release full genome assemblies of 38 newly sequenced avian species, link genome assembly downloads for the 7 of the remaining 10 species, and provide a guideline of genomic data that has been generated and used in our Avian Phylogenomics Project. To the best of our knowledge, the Avian Phylogenomics Project is the biggest vertebrate comparative genomics project to date. The genomic data presented here is expected to accelerate further analyses in many fields, including phylogenetics, comparative genomics, evolution, neurobiology, development biology, and other related areas.

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.