Midlatitude daytime D region ionosphere variations measured from radio atmospherics

We measured the midlatitude daytime ionospheric D region electron density profile height variations in July and August 2005 near Duke University by using radio atmospherics (or sferics for short), which are the high-power, broadband very low frequency (VLF) signals launched by lightning discharges. As expected, the measured daytime D region electron density profile heights showed temporal variations quantitatively correlated with solar zenith angle changes. In the midlatitude geographical regions near Duke University, the observed quiet time heights decreased from ∼80 km near sunrise to ∼71 km near noon when the solar zenith angle was minimum. The measured height quantitative dependence on the solar zenith angle was slightly different from the low-latitude measurement given in a previous work. We also observed unexpected spatial variations not linked to the solar zenith angle on some days, with 15% of days exhibiting regional differences larger than 0.5 km. In these 2 months, 14 days had sudden height drops caused by solar flare X-rays, with a minimum height of 63.4 km observed. The induced height change during a solar flare event was approximately proportional to the logarithm of the X-ray flux. In the long waveband (wavelength, 1-8 Å), an increase in flux by a factor of 10 resulted in 6.3 km decrease of the height at the flux peak time, nearly a perfect agreement with the previous measurement. During the rising and decaying phases of the solar flare, the height changes correlated more consistently with the short, rather than the long, wavelength X-ray flux changes. © 2010 by the American Geophysical Union.

The characterization of twenty sequenced human genomes.

We present the analysis of twenty human genomes to evaluate the prospects for identifying rare functional variants that contribute to a phenotype of interest. We sequenced at high coverage ten "case" genomes from individuals with severe hemophilia A and ten "control" genomes. We summarize the number of genetic variants emerging from a study of this magnitude, and provide a proof of concept for the identification of rare and highly-penetrant functional variants by confirming that the cause of hemophilia A is easily recognizable in this data set. We also show that the number of novel single nucleotide variants (SNVs) discovered per genome seems to stabilize at about 144,000 new variants per genome, after the first 15 individuals have been sequenced. Finally, we find that, on average, each genome carries 165 homozygous protein-truncating or stop loss variants in genes representing a diverse set of pathways.

Intergenic and genic sequence lengths have opposite relationships with respect to gene expression.

Eukaryotic genomes are mostly composed of noncoding DNA whose role is still poorly understood. Studies in several organisms have shown correlations between the length of the intergenic and genic sequences of a gene and the expression of its corresponding mRNA transcript. Some studies have found a positive relationship between intergenic sequence length and expression diversity between tissues, and concluded that genes under greater regulatory control require more regulatory information in their intergenic sequences. Other reports found a negative relationship between expression level and gene length and the interpretation was that there is selection pressure for highly expressed genes to remain small. However, a correlation between gene sequence length and expression diversity, opposite to that observed for intergenic sequences, has also been reported, and to date there is no testable explanation for this observation. To shed light on these varied and sometimes conflicting results, we performed a thorough study of the relationships between sequence length and gene expression using cell-type (tissue) specific microarray data in Arabidopsis thaliana. We measured median gene expression across tissues (expression level), expression variability between tissues (expression pattern uniformity), and expression variability between replicates (expression noise). We found that intergenic (upstream and downstream) and genic (coding and noncoding) sequences have generally opposite relationships with respect to expression, whether it is tissue variability, median, or expression noise. To explain these results we propose a model, in which the lengths of the intergenic and genic sequences have opposite effects on the ability of the transcribed region of the gene to be epigenetically regulated for differential expression. These findings could shed light on the role and influence of noncoding sequences on gene expression.

Structural basis for receptor subtype-specific regulation revealed by a chimeric beta 3/beta 2-adrenergic receptor.

The physiological significance of multiple G-protein-coupled receptor subtypes, such as the beta-adrenergic receptors (beta ARs), remains obscure, since in many cases several subtypes activate the same effector and utilize the same physiological agonists. We inspected the deduced amino acid sequences of the beta AR subtypes for variations in the determinants for agonist regulation as a potential basis for subtype differentiation. Whereas the beta 2AR has a C terminus containing 11 serine and threonine residues representing potential sites for beta AR kinase phosphorylation, which mediates rapid agonist-promoted desensitization, only 3 serines are present in the comparable region of the beta 3AR, and they are in a nonfavorable context. The beta 3AR also lacks sequence homology in regions which are important for agonist-mediated sequestration and down-regulation of the beta 2AR, although such determinants are less well defined. We therefore tested the idea that the agonist-induced regulatory properties of the two receptors might differ by expressing both subtypes in CHW cells and exposing them to the agonist isoproterenol. The beta 3AR did not display short-term agonist-promoted functional desensitization or sequestration, or long-term down-regulation. To assign a structural basis for these subtype-specific differences in agonist regulation, we constructed a chimeric beta 3/beta 2AR which comprised the beta 3AR up to proline-365 of the cytoplasmic tail and the C terminus of the beta 2AR. When cells expressing this chimeric beta 3/beta 2AR were exposed to isoproterenol, functional desensitization was observed. Whole-cell phosphorylation studies showed that the beta 2AR displayed agonist-dependent phosphorylation, but no such phosphorylation could be demonstrated with the beta 3AR, even when beta AR kinase was overexpressed. In contrast, the chimeric beta 3/beta 2AR did display agonist-dependent phosphorylation, consistent with its functional desensitization. In addition to conferring functional desensitization and phosphorylation to the beta 3AR, the C-terminal tail of the beta 2AR also conferred agonist-promoted sequestration and long-term receptor down-regulation.

Comparative analyses of seven algorithms for copy number variant identification from single nucleotide polymorphism arrays.

Determination of copy number variants (CNVs) inferred in genome wide single nucleotide polymorphism arrays has shown increasing utility in genetic variant disease associations. Several CNV detection methods are available, but differences in CNV call thresholds and characteristics exist. We evaluated the relative performance of seven methods: circular binary segmentation, CNVFinder, cnvPartition, gain and loss of DNA, Nexus algorithms, PennCNV and QuantiSNP. Tested data included real and simulated Illumina HumHap 550 data from the Singapore cohort study of the risk factors for Myopia (SCORM) and simulated data from Affymetrix 6.0 and platform-independent distributions. The normalized singleton ratio (NSR) is proposed as a metric for parameter optimization before enacting full analysis. We used 10 SCORM samples for optimizing parameter settings for each method and then evaluated method performance at optimal parameters using 100 SCORM samples. The statistical power, false positive rates, and receiver operating characteristic (ROC) curve residuals were evaluated by simulation studies. Optimal parameters, as determined by NSR and ROC curve residuals, were consistent across datasets. QuantiSNP outperformed other methods based on ROC curve residuals over most datasets. Nexus Rank and SNPRank have low specificity and high power. Nexus Rank calls oversized CNVs. PennCNV detects one of the fewest numbers of CNVs.

Interactions of chromatin context, binding site sequence content, and sequence evolution in stress-induced p53 occupancy and transactivation.

Cellular stresses activate the tumor suppressor p53 protein leading to selective binding to DNA response elements (REs) and gene transactivation from a large pool of potential p53 REs (p53REs). To elucidate how p53RE sequences and local chromatin context interact to affect p53 binding and gene transactivation, we mapped genome-wide binding localizations of p53 and H3K4me3 in untreated and doxorubicin (DXR)-treated human lymphoblastoid cells. We examined the relationships among p53 occupancy, gene expression, H3K4me3, chromatin accessibility (DNase 1 hypersensitivity, DHS), ENCODE chromatin states, p53RE sequence, and evolutionary conservation. We observed that the inducible expression of p53-regulated genes was associated with the steady-state chromatin status of the cell. Most highly inducible p53-regulated genes were suppressed at baseline and marked by repressive histone modifications or displayed CTCF binding. Comparison of p53RE sequences residing in different chromatin contexts demonstrated that weaker p53REs resided in open promoters, while stronger p53REs were located within enhancers and repressed chromatin. p53 occupancy was strongly correlated with similarity of the target DNA sequences to the p53RE consensus, but surprisingly, inversely correlated with pre-existing nucleosome accessibility (DHS) and evolutionary conservation at the p53RE. Occupancy by p53 of REs that overlapped transposable element (TE) repeats was significantly higher (p<10-7) and correlated with stronger p53RE sequences (p<10-110) relative to nonTE-associated p53REs, particularly for MLT1H, LTR10B, and Mer61 TEs. However, binding at these elements was generally not associated with transactivation of adjacent genes. Occupied p53REs located in L2-like TEs were unique in displaying highly negative PhyloP scores (predicted fast-evolving) and being associated with altered H3K4me3 and DHS levels. These results underscore the systematic interaction between chromatin status and p53RE context in the induced transactivation response. This p53 regulated response appears to have been tuned via evolutionary processes that may have led to repression and/or utilization of p53REs originating from primate-specific transposon elements.

Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection.

DNaseI footprinting is an established assay for identifying transcription factor (TF)-DNA interactions with single base pair resolution. High-throughput DNase-seq assays have recently been used to detect in vivo DNase footprints across the genome. Multiple computational approaches have been developed to identify DNase-seq footprints as predictors of TF binding. However, recent studies have pointed to a substantial cleavage bias of DNase and its negative impact on predictive performance of footprinting. To assess the potential for using DNase-seq to identify individual binding sites, we performed DNase-seq on deproteinized genomic DNA and determined sequence cleavage bias. This allowed us to build bias corrected and TF-specific footprint models. The predictive performance of these models demonstrated that predicted footprints corresponded to high-confidence TF-DNA interactions. DNase-seq footprints were absent under a fraction of ChIP-seq peaks, which we show to be indicative of weaker binding, indirect TF-DNA interactions or possible ChIP artifacts. The modeling approach was also able to detect variation in the consensus motifs that TFs bind to. Finally, cell type specific footprints were detected within DNase hypersensitive sites that are present in multiple cell types, further supporting that footprints can identify changes in TF binding that are not detectable using other strategies.

Quantitative genetics of CTCF binding reveal local sequence effects and different modes of X-chromosome association.

Associating genetic variation with quantitative measures of gene regulation offers a way to bridge the gap between genotype and complex phenotypes. In order to identify quantitative trait loci (QTLs) that influence the binding of a transcription factor in humans, we measured binding of the multifunctional transcription and chromatin factor CTCF in 51 HapMap cell lines. We identified thousands of QTLs in which genotype differences were associated with differences in CTCF binding strength, hundreds of them confirmed by directly observable allele-specific binding bias. The majority of QTLs were either within 1 kb of the CTCF binding motif, or in linkage disequilibrium with a variant within 1 kb of the motif. On the X chromosome we observed three classes of binding sites: a minority class bound only to the active copy of the X chromosome, the majority class bound to both the active and inactive X, and a small set of female-specific CTCF sites associated with two non-coding RNA genes. In sum, our data reveal extensive genetic effects on CTCF binding, both direct and indirect, and identify a diversity of patterns of CTCF binding on the X chromosome.

Impacts of land use change and climate variations on annual inflow into Miyun Reservoir,Beijing, China

The Miyun Reservoir, the only surface water source for Beijing city, has experienced water supply decline in recent decades. Previous studies suggest that both land use change and climate contribute to the changes of water supply in this critical watershed. However, the specific causes of the decline in the Miyun Reservoir are debatable under a non-stationary climate in the past 4 decades. The central objective of this study was to quantify the separate and collective contributions of land use change and climate variability to the decreasing inflow into the Miyun Reservoir during 1961–2008. Different from previous studies on this watershed, we used a comprehensive approach to quantify the timing of changes in hydrology and associated environmental variables using the long-term historical hydrometeorology and remote-sensing-based land use records. To effectively quantify the different impacts of the climate variation and land use change on streamflow during different sub-periods, an annual water balance model (AWB), the climate elasticity model (CEM), and a rainfall–runoff model (RRM) were employed to conduct attribution analysis synthetically. We found a significant (p  <  0.01) decrease in annual streamflow, a significant positive trend in annual potential evapotranspiration (p  <  0.01), and an insignificant (p  >  0.1) negative trend in annual precipitation during 1961–2008. We identified two streamflow breakpoints, 1983 and 1999, by the sequential Mann–Kendall test and double-mass curve. Climate variability alone did not explain the decrease in inflow to the Miyun Reservoir. Reduction of water yield was closely related to increase in actual evapotranspiration due to the expansion of forestland and reduction in cropland and grassland, and was likely exacerbated by increased water consumption for domestic and industrial uses in the basin. The contribution to the observed streamflow decline from land use change fell from 64–92 % during 1984–1999 to 36–58 % during 2000–2008, whereas the contribution from climate variation climbed from 8–36 % during the 1984–1999 to 42–64 % during 2000–2008. Model uncertainty analysis further demonstrated that climate warming played a dominant role in streamflow reduction in the most recent decade (i.e., 2000s). We conclude that future climate change and variability will further challenge the water supply capacity of the Miyun Reservoir to meet water demand. A comprehensive watershed management strategy needs to consider the climate variations besides vegetation management in the study basin.