Mass conservative three-dimensional water tracer distribution from Markov chain Monte Carlo inversion of time-lapse ground-penetrating radar data

Time-lapse geophysical measurements are widely used to monitor the movement of water and solutes through the subsurface. Yet commonly used deterministic least squares inversions typically suffer from relatively poor mass recovery, spread overestimation, and limited ability to appropriately estimate nonlinear model uncertainty. We describe herein a novel inversion methodology designed to reconstruct the three-dimensional distribution of a tracer anomaly from geophysical data and provide consistent uncertainty estimates using Markov chain Monte Carlo simulation. Posterior sampling is made tractable by using a lower-dimensional model space related both to the Legendre moments of the plume and to predefined morphological constraints. Benchmark results using cross-hole ground-penetrating radar travel times measurements during two synthetic water tracer application experiments involving increasingly complex plume geometries show that the proposed method not only conserves mass but also provides better estimates of plume morphology and posterior model uncertainty than deterministic inversion results.

Relationships between grain yield and accumulation of biomass, nitrogen and phosphorus in common bean cultivars

Shoot biomass is considered a relevant component for crop yield, but relationships between biological productivity and grain yield in legume crops are usually difficult to establish. Two field experiments were carried out to investigate the relationships between grain yield, biomass production and N and P accumulation at reproductive stages of common bean (Phaseolus vulgaris) cultivars. Nine and 18 cultivars were grown on 16 m² plots in 1998 and 1999, respectively, with four replications. Crop biomass was sampled at four growth stages (flowering R6, pod setting R7, beginning of pod filling R8, and mid-pod filling R8.5), grain yield was measured at maturity, and N and P concentrations were determined in plant tissues. In both years, bean cultivars differed in grain yield, in root mass at R6 and R7 stages, and in shoot mass at R6 and R8.5, whereas at R7 and R8 differences in shoot mass were significant in 1998 only. In both years, grain yield did not correlate with shoot mass at R6 and R7 and with root mass at R6. Grain yield correlated with shoot mass at R8 in 1999 but not in 1998, with shoot mass at R8.5 and with root mass at R7 in both years. Path coefficient analysis indicated that shoot mass at R8.5 had a direct effect on grain yield in both years, that root mass at R7 had a direct effect on grain yield in 1998, and that in 1999 the amounts of N and P in shoots at R8.5 had indirect effects on grain yield via shoot mass at R8.5. A combined analysis of both experiments revealed that biomass accumulation, N and P in shoots at R6 and R7 as well as root mass at R6 were similar in both years. In 1998 however bean accumulated more root mass at R7 and more biomass and N and P in shoots at R8 and R8.5, resulting in a 57 % higher grain yield in 1998. This indicates that grain yield of different common bean cultivars is not intrinsically associated with vegetative vigor at flowering and that mechanisms during pod filling can strongly influence the final crop yield. The establishment of a profuse root system during pod setting, associated with the continuous N and P acquisition during early pod filling, seems to be relevant for higher grain yields of common bean.

Common variants in the ATP2B1 gene are associated with susceptibility to hypertension: the Japanese Millennium Genome Project.

Hypertension is one of the most common complex genetic disorders. We have described previously 38 single nucleotide polymorphisms (SNPs) with suggestive association with hypertension in Japanese individuals. In this study we extend our previous findings by analyzing a large sample of Japanese individuals (n=14 105) for the most associated SNPs. We also conducted replication analyses in Japanese of susceptibility loci for hypertension identified recently from genome-wide association studies of European ancestries. Association analysis revealed significant association of the ATP2B1 rs2070759 polymorphism with hypertension (P=5.3×10(-5); allelic odds ratio: 1.17 [95% CI: 1.09 to 1.26]). Additional SNPs in ATP2B1 were subsequently genotyped, and the most significant association was with rs11105378 (odds ratio: 1.31 [95% CI: 1.21 to 1.42]; P=4.1×10(-11)). Association of rs11105378 with hypertension was cross-validated by replication analysis with the Global Blood Pressure Genetics consortium data set (odds ratio: 1.13 [95% CI: 1.05 to 1.21]; P=5.9×10(-4)). Mean adjusted systolic blood pressure was highly significantly associated with the same SNP in a meta-analysis with individuals of European descent (P=1.4×10(-18)). ATP2B1 mRNA expression levels in umbilical artery smooth muscle cells were found to be significantly different among rs11105378 genotypes. Seven SNPs discovered in published genome-wide association studies were also genotyped in the Japanese population. In the combined analysis with replicated 3 genes, FGF5 rs1458038, CYP17A1, rs1004467, and CSK rs1378942, odds ratio of the highest risk group was 2.27 (95% CI: 1.65 to 3.12; P=4.6×10(-7)) compared with the lower risk group. In summary, this study confirmed common genetic variation in ATP2B1, as well as FGF5, CYP17A1, and CSK, to be associated with blood pressure levels and risk of hypertension.

Testing for multicointegration in panel data with common factors

The paper addresses the concept of multicointegration in panel data frame- work. The proposal builds upon the panel data cointegration procedures developed in Pedroni (2004), for which we compute the moments of the parametric statistics. When individuals are either cross-section independent or cross-section dependence can be re- moved by cross-section demeaning, our approach can be applied to the wider framework of mixed I(2) and I(1) stochastic processes analysis. The paper also deals with the issue of cross-section dependence using approximate common factor models. Finite sample performance is investigated through Monte Carlo simulations. Finally, we illustrate the use of the procedure investigating inventories, sales and production relationship for a panel of US industries.

HCFC1 is a common component of active human CpG-island promoters and coincides with ZNF143, THAP11, YY1, and GABP transcription factor occupancy.

In human transcriptional regulation, DNA-sequence-specific factors can associate with intermediaries that orchestrate interactions with a diverse set of chromatin-modifying enzymes. One such intermediary is HCFC1 (also known as HCF-1). HCFC1, first identified in herpes simplex virus transcription, has a poorly defined role in cellular transcriptional regulation. We show here that, in HeLa cells, HCFC1 is observed bound to 5400 generally active CpG-island promoters. Examination of the DNA sequences underlying the HCFC1-binding sites revealed three sequence motifs associated with the binding of (1) ZNF143 and THAP11 (also known as Ronin), (2) GABP, and (3) YY1 sequence-specific transcription factors. Subsequent analysis revealed colocalization of HCFC1 with these four transcription factors at ∼90% of the 5400 HCFC1-bound promoters. These studies suggest that a relatively small number of transcription factors play a major role in HeLa-cell transcriptional regulation in association with HCFC1.

Common Variable Immunodeficiency: molecular pathways and clinical manifestations

Common variable immunodeficiency (CVID), is a disease that is characterized by hypogammaglobulinemia as well as a defect in T, B and dendritic cells. This leads to recurrent bacterial infection mainly caused by Streptococcus pneumoniae, Klebsiella pneumoniae and Haemophilus influenzae, as well as inflammatory manifestations, i.e. granulomateous disease, gastro-intestinal disorders and chronic lung disease. Intravenous Immunoglobulin (IVIg) therapy reduces CVID susceptibility to bacterial infections to some extend. We analyzed clinical aspects of patients from our database. We recently showed that bacteria-specific CD4 T cells of CVID patients were impaired. We therefor postulated that CVID patients may harbor an acquired T-cell deficiency also called exhaustion. To test this hypothesis, we performed a comprehensive investigation of the functional profiles of bacteria-specific CD4 T cells isolated from 31 healthy individuals and 30 CVID patients. In the present study, we demonstrated that bacteria-specific but not virus-specific CD4 T cells in CVID patients harbored reduced proliferation capacity and expressed high level of PD-1. Interestingly, the blockade of PD-1/PD-1 ligands interactions restored partially bacteria but not virus-specific CD4 T-cell proliferation. Finally, we showed that 1) the level of endotoxins inversely correlates with IgG concentration, 2) IVIG treated CVID patients harbored reduced endotoxemia and 3) IgG concentration exceeding 7 mg/mL strongly reduces both the proportion of CVID patients with detectable endotoxemia and the concentration of endotoxins in plasma. Taken together our observations, suggest that primary B-cell defect(s) in CVID patients leads to recurrent bacterial infections that are associated to an acquired (secondary) impairment of CD4 T cells which may in turn exacerbate the lack of protection against extracellular bacteria.

Integration of ground-penetrating radar and microgravimetric methods to map shallow caves

Ground-penetrating radar (GPR) and microgravimetric surveys have been conducted in the southern Jura mountains of western Switzerland in order to map subsurface karstic features. The study site, La Grande Rolaz cave, is an extensive system in which many portions have been mapped. By using small station spacing and careful processing for the geophysical data, and by modeling these data with topographic information from within the cave, accurate interpretations have been achieved. The constraints on the interpreted geologic models are better when combining the geophysical methods than when using only one of the methods, despite the general limitations of two-dimensional (2D) profiling. For example, microgravimetry can complement GPR methods for accurately delineating a shallow cave section approximately 10 X 10 mt in size. Conversely, GPR methods can be complementary in determining cavity depths and in verifying the presence of off-line features and numerous areas of small cavities and fractures, which may be difficult to resolve in microgravimetric data.

Domain Adaptation in remote sensing: increasing the portability of land-cover classifiers

Among the types of remote sensing acquisitions, optical images are certainly one of the most widely relied upon data sources for Earth observation. They provide detailed measurements of the electromagnetic radiation reflected or emitted by each pixel in the scene. Through a process termed supervised land-cover classification, this allows to automatically yet accurately distinguish objects at the surface of our planet. In this respect, when producing a land-cover map of the surveyed area, the availability of training examples representative of each thematic class is crucial for the success of the classification procedure. However, in real applications, due to several constraints on the sample collection process, labeled pixels are usually scarce. When analyzing an image for which those key samples are unavailable, a viable solution consists in resorting to the ground truth data of other previously acquired images. This option is attractive but several factors such as atmospheric, ground and acquisition conditions can cause radiometric differences between the images, hindering therefore the transfer of knowledge from one image to another. The goal of this Thesis is to supply remote sensing image analysts with suitable processing techniques to ensure a robust portability of the classification models across different images. The ultimate purpose is to map the land-cover classes over large spatial and temporal extents with minimal ground information. To overcome, or simply quantify, the observed shifts in the statistical distribution of the spectra of the materials, we study four approaches issued from the field of machine learning. First, we propose a strategy to intelligently sample the image of interest to collect the labels only in correspondence of the most useful pixels. This iterative routine is based on a constant evaluation of the pertinence to the new image of the initial training data actually belonging to a different image. Second, an approach to reduce the radiometric differences among the images by projecting the respective pixels in a common new data space is presented. We analyze a kernel-based feature extraction framework suited for such problems, showing that, after this relative normalization, the cross-image generalization abilities of a classifier are highly increased. Third, we test a new data-driven measure of distance between probability distributions to assess the distortions caused by differences in the acquisition geometry affecting series of multi-angle images. Also, we gauge the portability of classification models through the sequences. In both exercises, the efficacy of classic physically- and statistically-based normalization methods is discussed. Finally, we explore a new family of approaches based on sparse representations of the samples to reciprocally convert the data space of two images. The projection function bridging the images allows a synthesis of new pixels with more similar characteristics ultimately facilitating the land-cover mapping across images.

Drift compensation of gas sensor array data by common principal component analysis

A new drift compensation method based on Common Principal Component Analysis (CPCA) is proposed. The drift variance in data is found as the principal components computed by CPCA. This method finds components that are common for all gasses in feature space. The method is compared in classification task with respect to the other approaches published where the drift direction is estimated through a Principal Component Analysis (PCA) of a reference gas. The proposed new method ¿ employing no specific reference gas, but information from all gases ¿has shown the same performance as the traditional approach with the best-fitted reference gas. Results are shown with data lasting 7-months including three gases at different concentrations for an array of 17 polymeric sensors.

Testing for multicointegration in panel data with common factors

The paper addresses the concept of multicointegration in panel data frame- work. The proposal builds upon the panel data cointegration procedures developed in Pedroni (2004), for which we compute the moments of the parametric statistics. When individuals are either cross-section independent or cross-section dependence can be re- moved by cross-section demeaning, our approach can be applied to the wider framework of mixed I(2) and I(1) stochastic processes analysis. The paper also deals with the issue of cross-section dependence using approximate common factor models. Finite sample performance is investigated through Monte Carlo simulations. Finally, we illustrate the use of the procedure investigating inventories, sales and production relationship for a panel of US industries.

A common ancestor DNA motif for invertebrate and vertebrate hormone response elements.

The ecdysone-responsive DNA sequence of the Drosophila hsp27 gene promoter contains four direct and inverted repeats reminiscent of those that compose the vertebrate palindromic estrogen response element (ERE) and the thyroid hormone/retinoic acid response element (TRE/RRE). Interestingly, a 3 bp substitution in the wild-type Hsp27 ecdysone response element (EcdRE) increases both its similarity with the vertebrate ERE and TRE/RRE and its capacity to confer ecdysone responsiveness to a heterologous promoter. Remarkably, increasing the spacing between the inverted repeats of this strong EcdRE by two nucleotides converts it into an ERE. Inversely, decreasing the spacing between the two inverted repeats of the vertebrate consensus palindromic ERE, from three to one nucleotide, converts it into a functional EcdRE. Thus, the only difference between an invertebrate EcdRE and a vertebrate palindromic ERE or TRE/RRE is in the spacing between the conserved inverted repeated motifs forming these palindromic HREs. The finding that the sequence motif 5'-GGTCA-3' present in the vertebrate ERE and TRE/RRE is also a functionally important characteristic of an invertebrate HRE, suggests that a common ancestor regulatory DNA sequence gave rise to all HREs known so far. We discuss the possibility that this progenitor motif is the GGTCA sequence.

Time consistent Pareto solutions in common access resource gameswith asymmetric players

In the analysis of equilibrium policies in a di erential game, if agents have different time preference rates, the cooperative (Pareto optimum) solution obtained by applying the Pontryagin's Maximum Principle becomes time inconsistent. In this work we derive a set of dynamic programming equations (in discrete and continuous time) whose solutions are time consistent equilibrium rules for N-player cooperative di erential games in which agents di er in their instantaneous utility functions and also in their discount rates of time preference. The results are applied to the study of a cake-eating problem describing the management of a common property exhaustible natural resource. The extension of the results to a simple common property renewable natural resource model in in nite horizon is also discussed.

Leaf area of common bean genotypes during early pod filling as related to plant adaptation to limited phosphorus supply

Low phosphorus supply markedly limits leaf growth and genotypes able to maintain adequate leaf area at low P could adapt better to limited-P conditions. This work aimed to investigate the relationship between leaf area production of common bean (Phaseolus vulgaris) genotypes during early pod filling and plant adaptation to limited P supply. Twenty-four genotypes, comprised of the four growth habits in the species and two weedy accessions, were grown at two P level applied to the soil (20 and 80 mg kg-1) in 4 kg pots and harvested at two growth stages (pod setting and early pod filling). High P level markedly increased the leaf number and leaf size (leaf area per leaf), slightly increased specific leaf area but did not affect the net assimilation rate. At low P level most genotypic variation for plant dry mass was associated with leaf size, whereas at high P level this variation was associated primarily with the number of leaves and secondarily with leaf size, specific leaf area playing a minor role at both P level. Determinate bush genotypes presented a smaller leaf area, fewer but larger leaves with higher specific leaf area and lower net assimilation rate. Climbing genotypes showed numerous leaves, smaller and thicker leaves with a higher net assimilation rate. Indeterminate bush and indeterminate prostrate genotypes presented the highest leaf area, achieved through intermediate leaf number, leaf size and specific leaf area. The latter groups were better adapted to limited P. It is concluded that improved growth at low P during early pod filling was associated with common bean genotypes able to maintain leaf expansion through leaves with greater individual leaf area.

Uso de sistema GPR (Ground Penetring Radar) na avaliação de atributos de um solo sob plantio de cana-de-açúcar

O GPR utiliza técnica indireta para investigação das estruturas de subsuperfície, sendo caracterizado como método não-invasivo, o qual permite extrair informações ao longo do perfil de solo sem perfurar, sondar ou escavar. Este trabalho teve o objetivo de avaliar a utilização do GPR no estabelecimento de atributos físicos do solo em área de cultivo com cana-de-açúcar, além de estabelecer sua correlação com o método do cilindro volumétrico. A área de estudo está localizada no município de Bom Sucesso - PR, sendo utilizada com cultura da cana há mais de 20 anos. Foram obtidas amostras de solos indeformadas com cilindro volumétrico nas profundidades de 0-0,1, 0,1-0,2, 0,2-0,3, 0,3-0,4 e 0,4-0,5 m, ao longo de três transeptos paralelos, sendo o central referente à linha de plantio (transepto 2) e os outros dois correspondentes às entrelinhas laterais (transeptos 1 e 3). Os valores de densidade e porosidade total, obtidos por meio das amostras físicas, foram inicialmente interpretados para verificar o estado do solo nas linhas e nas entrelinhas de cana, nas cinco profundidades. Posteriormente, foram correlacionados com os valores das constantes dielétricas (K) adquiridos pelo georradar, porém apenas para a camada de 0,4-0,5 m, uma vez que as camadas superficiais até a profundidade de 0,4 m não puderam ser avaliadas devido às limitações técnicas do equipamento, aliadas às condições físicas e químicas do solo, que acabaram por prejudicar as reflexões. A densidade mostrou-se menor nas linhas, quando comparada às entrelinhas, para as camadas de 0-0,1 e 0,2-0,3 m. A porosidade total foi menor nas entrelinhas do que nas linhas na camada de 0-0,1 m. Os valores de K foram maiores nas linhas, onde a umidade e a porosidade total foram maiores, e menores nas entrelinhas, onde as densidades encontravam se próximas ao limite crítico. A correlação de K com a densidade e porosidade total não foi significativa. Pequenas modificações na porosidade do solo em função do aumento da densidade foram suficientes para diminuir o valor de K nas entrelinhas, devido à redução na quantidade de água nos poros. Esse fato pode ser confirmado pela correlação de K com a umidade. Assim, pode-se dizer que o georradar foi extremamente sensível à variação de umidade no solo, apresentando melhores resultados para esta variável do que para a densidade e porosidade total determinada por meio do cilindro volumétrico.