Two high throughput technologies to detect segmental aneuploidies identify new Williams-Beuren syndrome patients with atypical deletions.

OBJECTIVE: To develop and compare two new technologies for diagnosing a contiguous gene syndrome, the Williams-Beuren syndrome (WBS). METHODS: The first proposed method, named paralogous sequence quantification (PSQ), is based on the use of paralogous sequences located on different chromosomes and quantification of specific mismatches present at these loci using pyrosequencing technology. The second exploits quantitative real time polymerase chain reaction (QPCR) to assess the relative quantity of an analysed locus. RESULTS: A correct and unambiguous diagnosis was obtained for 100% of the analysed samples with either technique (n = 165 and n = 155, respectively). These methods allowed the identification of two patients with atypical deletions in a cohort of 182 WBS patients. Both patients presented with mild facial anomalies, mild mental retardation with impaired visuospatial cognition, supravalvar aortic stenosis, and normal growth indices. These observations are consistent with the involvement of GTF2IRD1 or GTF2I in some of the WBS facial features. CONCLUSIONS: Both PSQ and QPCR are robust, easy to interpret, and simple to set up. They represent a competitive alternative for the diagnosis of segmental aneuploidies in clinical laboratories. They have advantages over fluorescence in situ hybridisation or microsatellites/SNP genotyping for detecting short segmental aneuploidies as the former is costly and labour intensive while the latter depends on the informativeness of the polymorphisms.

Prognostic role of KRAS and BRAF in stage II and III resected colon cancer: results of the translational study on the PETACC-3, EORTC 40993, SAKK 60-00 trial.

PURPOSE: Mutations within the KRAS proto-oncogene have predictive value but are of uncertain prognostic value in the treatment of advanced colorectal cancer. We took advantage of PETACC-3, an adjuvant trial with 3,278 patients with stage II to III colon cancer, to evaluate the prognostic value of KRAS and BRAF tumor mutation status in this setting. PATIENTS AND METHODS: Formalin-fixed paraffin-embedded tissue blocks (n = 1,564) were prospectively collected and DNA was extracted from tissue sections from 1,404 cases. Planned analysis of KRAS exon 2 and BRAF exon 15 mutations was performed by allele-specific real-time polymerase chain reaction. Survival analyses were based on univariate and multivariate proportional hazard regression models. RESULTS: KRAS and BRAF tumor mutation rates were 37.0% and 7.9%, respectively, and were not significantly different according to tumor stage. In a multivariate analysis containing stage, tumor site, nodal status, sex, age, grade, and microsatellite instability (MSI) status, KRAS mutation was associated with grade (P = .0016), while BRAF mutation was significantly associated with female sex (P = .017), and highly significantly associated with right-sided tumors, older age, high grade, and MSI-high tumors (all P < 10(-4)). In univariate and multivariate analysis, KRAS mutations did not have a major prognostic value regarding relapse-free survival (RFS) or overall survival (OS). BRAF mutation was not prognostic for RFS, but was for OS, particularly in patients with MSI-low (MSI-L) and stable (MSI-S) tumors (hazard ratio, 2.2; 95% CI, 1.4 to 3.4; P = .0003). CONCLUSION: In stage II-III colon cancer, the KRAS mutation status does not have major prognostic value. BRAF is prognostic for OS in MS-L/S tumors.

Foxm1 expression in prostate epithelial cells is essential for prostate carcinogenesis.

The treatment of advanced prostate cancer (PCa) remains a challenge. Identification of new molecular mechanisms that regulate PCa initiation and progression would provide targets for the development of new cancer treatments. The Foxm1 transcription factor is highly up-regulated in tumor cells, inflammatory cells, and cells of tumor microenvironment. However, its functions in different cell populations of PCa lesions are unknown. To determine the role of Foxm1 in tumor cells during PCa development, we generated two novel transgenic mouse models, one exhibiting Foxm1 gain-of-function and one exhibiting Foxm1 loss-of-function under control of the prostate epithelial-specific Probasin promoter. In the transgenic adenocarcinoma mouse prostate (TRAMP) model of PCa that uses SV40 large T antigen to induce PCa, loss of Foxm1 decreased tumor growth and metastasis. Decreased prostate tumorigenesis was associated with a decrease in tumor cell proliferation and the down-regulation of genes critical for cell proliferation and tumor metastasis, including Cdc25b, Cyclin B1, Plk-1, Lox, and Versican. In addition, tumor-associated angiogenesis was decreased, coinciding with reduced Vegf-A expression. The mRNA and protein levels of 11β-Hsd2, an enzyme playing an important role in tumor cell proliferation, were down-regulated in Foxm1-deficient PCa tumors in vivo and in Foxm1-depleted TRAMP C2 cells in vitro. Foxm1 bound to, and increased transcriptional activity of, the mouse 11β-Hsd2 promoter through the -892/-879 region, indicating that 11β-Hsd2 was a direct transcriptional target of Foxm1. Without TRAMP, overexpression of Foxm1 either alone or in combination with inhibition of a p19(ARF) tumor suppressor caused a robust epithelial hyperplasia, but was insufficient to induce progression from hyperplasia to PCa. Foxm1 expression in prostate epithelial cells is critical for prostate carcinogenesis, suggesting that inhibition of Foxm1 is a promising therapeutic approach for prostate cancer chemotherapy.

Cytosolic sensing of extracellular self-DNA transported into monocytes by the antimicrobial peptide LL37.

The intracellular location of nucleic acid sensors prevents recognition of extracellular self-DNA released by dying cells. However, on forming a complex with the endogenous antimicrobial peptide LL37, extracellular DNA is transported into endosomal compartments of plasmacytoid dendritic cells, leading to activation of Toll-like receptor-9 and induction of type I IFNs. Whether LL37 also transports self-DNA into nonplasmacytoid dendritic cells, leading to type I IFN production via other intracellular DNA receptors is unknown. Here we found that LL37 very efficiently transports self-DNA into monocytes, leading the production of type I IFNs in a Toll-like receptor-independent manner. This type I IFN induction was mediated by double-stranded B form DNA, regardless of its sequence, CpG content, or methylation status, and required signaling through the adaptor protein STING and TBK1 kinase, indicating the involvement of cytosolic DNA sensors. Thus, our study identifies a novel link between the antimicrobial peptides and type I IFN responses involving DNA-dependent activation of cytosolic sensors in monocytes.

Crohn's disease-associated polymorphism within the PTPN2 gene affects muramyl-dipeptide-induced cytokine secretion and autophagy.

BACKGROUND: The single nucleotide polymorphism (SNP) rs2542151 within the gene locus region encoding protein tyrosine phosphatase non-receptor type 2 (PTPN2) has been associated with Crohn's disease (CD), ulcerative colitis (UC), type-I diabetes, and rheumatoid arthritis. We have previously shown that PTPN2 regulates mitogen-activated protein kinase (MAPK) signaling and cytokine secretion in human THP-1 monocytes and intestinal epithelial cells (IEC). Here, we studied whether intronic PTPN2 SNP rs1893217 regulates immune responses to the nucleotide-oligomerization domain 2 (NOD2) ligand, muramyl-dipeptide (MDP). MATERIALS AND METHODS: Genomic DNA samples from 343 CD and 663 non-IBD control patients (male and female) from a combined German, Swiss, and Polish cohort were genotyped for the presence of the PTPN2 SNPs, rs2542151, and rs1893217. PTPN2-variant rs1893217 was introduced into T(84) IEC or THP-1 cells using a lentiviral vector. RESULTS: We identified a novel association between the genetic variant, rs1893217, located in intron 7 of the PTPN2 gene and CD. Human THP-1 monocytes carrying this variant revealed increased MAPK activation as well as elevated mRNA expression of T-bet transcription factor and secretion of interferon-γ in response to the bacterial wall component, MDP. In contrast, secretion of interleukin-8 and tumor necrosis factor were reduced. In both, T(84) IEC and THP-1 monocytes, autophagosome formation was impaired. CONCLUSIONS: We identified a novel CD-associated PTPN2 variant that modulates innate immune responses to bacterial antigens. These findings not only provide key insights into the effects of a functional mutation on a clinically relevant gene, but also reveal how such a mutation could contribute to the onset of disease.

Lack of association between beta-herpesvirus infection and bronchiolitis obliterans syndrome in lung transplant recipients in the era of antiviral prophylaxis.

BACKGROUND: Cytomegalovirus (CMV), human herpesvirus-6 and -7 (HHV-6 and -7) are beta-herpesviruses that commonly reactivate and have been proposed to trigger acute rejection and chronic allograft injury. We assessed the contribution of these viruses in the development of bronchiolitis obliterans syndrome (BOS) after lung transplantation. METHODS: Quantitative real-time polymerase chain reaction of bronchoalveolar lavage samples were performed for CMV, HHV-6 and -7 in a prospective cohort of lung transplant recipients. A time-dependent Cox regression analysis was used to correlate the risk of BOS and acute rejection in patients with and without beta-herpesviruses infection. RESULTS: Ninety-three patients were included in the study over a period of 3 years. A total of 581 samples from bronchoalveolar lavage were obtained. Sixty-one patients (65.6%) had at least one positive result for one of the beta-herpesviruses: 48 patients (51.6%) for CMV and 19 patients (20.4%) for both HHV-6 and -7. Median peak viral load was 3419 copies/mL for CMV, 258 copies/mL for HHV-6, and 665 copies/mL for HHV-7. Acute rejection (>or=grade 2) occurred in 46.2% and BOS (>or=stage 1) in 19.4% of the patients. In the Cox regression model the relative risk of acute rejection or BOS was not increased in patients with any beta-herpesviruses reactivation. Acute rejection was the only independently associated risk factor for BOS. CONCLUSIONS: In lung transplant recipients receiving prolonged antiviral prophylaxis, reactivation of beta-herpesviruses within the allograft was common. However, despite high viral loads in many patients, virus replication was not associated with the development of rejection or BOS.

Infrequent replication of parvovirus B19 and erythrovirus genotypes 2 and 3 among HIV-infected patients with chronic anemia.

We investigated the role that erythroviruses (parvovirus B19 and erythrovirus genotypes 2 and 3) play in the lives of immunosuppressed HIV-infected patients with chronic anemia. We screened the serum samples of 428 patients by specific ultrasensitive real-time polymerase chain reaction assay. Sixteen patients had circulating DNA, with no apparent clinical impact. Erythrovirus-associated anemia is an extremely rare event in HIV-infected patients.

Improved detection of Rhodococcus coprophilus with a new quantitative PCR assay.

Agricultural practices, such as spreading liquid manure or the utilisation of land as animal pastures, can result in faecal contamination of water resources. Rhodococcus coprophilus is used in microbial source tracking to indicate animal faecal contamination in water. Methods previously described for detecting of R. coprophilus in water were neither sensitive nor specific. Therefore, the aim of this study was to design and validate a new quantitative polymerase chain reaction (qPCR) to improve the detection of R. coprophilus in water. The new PCR assay was based on the R. coprophilus 16S rRNA gene. The validation showed that the new approach was specific and sensitive for deoxyribunucleic acid from target host species. Compared with other PCR assays tested in this study, the detection limit of the new qPCR was between 1 and 3 log lower. The method, including a filtration step, was further validated and successfully used in a field investigation in Switzerland. Our work demonstrated that the new detection method is sensitive and robust to detect R. coprophilus in surface and spring water. Compared with PCR assays that are available in the literature or to the culture-dependent method, the new molecular approach improves the detection of R. coprophilus.

Population genetic structure in the native and introduced range of the Argentine ant/

SUMMARY : The evolution of animal societies, where some individuals forego their own reproductive opportunities to help others to reproduce, poses an evolutionary paradox that can be traced back to Darwin. Altruism may evolve through kin selection when the donor and recipient of altruistic acts are related to each other. In social insects, workers are generally highly related to the brood they rear when colonies are headed by a single queen. Yet some ants have an extraordinary social organization, called unicoloniality, whereby individuals from separate nests mix freely to form large supercolonies, which in some cases extend over hundreds of km. These supercolonies are characterised by a high number of queens (polygyny) and an absence of clear colony boundaries. This type of social organization represents an evolutionary paradox because relatedness between nestmates is effectively zero. In such conditions, kin selection cannot account for the evolution of reproductive altruism. Moreover, unicoloniality is thought to be unstable over time, because workers that can no longer aid close relatives may evolve more selfish strategies. The Argentine ant (Linepithema humile) is a highly invasive species listed among the hundred world's worst invaders by the UICN. Native from South America, L. humile has been accidentally introduced throughout the world. Native populations have been described as noninvasive with a family-based organization. In contrast, within its introduction range, they form unicolonial supercolonies that contain numerous nests without intraspecific aggression. The development of such unicolonial populations has been explained as a direct consequence of the ant's introduction into a new habitat, favouring a transition from family-based to open colonies. To determine if the social structure of the Argentine ant is fundamentally different between the native and the introduced range, we studied genetically and behaviourally native and introduced populations of L. humile over different geographic scales. Our results clearly indicated that there are no fundamental differences in the social organisation of the Argentine ant between the two ranges. Our investigations revealed that, contrary to previous claims, native populations have a unicolonial social organisation very similar to that observed in the introduced range. Consequently, the unicolonial social structure of the Argentine ant does not stem from a shift in social organization associated with introduction into new habitats but evolved in the native range and is likely a stable, evolutionarily ancient adaptation to the local environment. Our study on native populations of L. humile also gave important insight in the comprehension of the evolution of unicoloniality in the Argentine ant. Native supercolonies are relatively small compared to introduced ones and may co-habit in a same population. These supercolonies are genetically highly differentiated leading to a significant relatedness among nestmate workers when the different supercolonies of a population are taken as a reference population. This provides the necessary conditions for loin selection to operate. Furthermore, we examined a native population over time, which revealed a high supercolony extinction rate. If more competitive supercolonies are more likely to survive or replace other supercolonies, a subtle dynamical process between the spread of selfish traits within supercolony and the selective elimination of supercolonies with such traits may allow a stable equilibrium and the persistence of unicoloniality over time. Finally, a worldwide study of the Argentine ant showed that the introduced supercolonies originate from numerous independent introduction events. In conclusion, the success of the Argentine ant does not stem from a shift in social organization associated with its introduction into new habitats, but is most probably explained by the intrinsic characteristics developed in its native range. RESUME : L'altruisme de reproduction où certains individus renoncent à leur propre reproduction pour aider d'autres individus à se reproduire constitue l'un des plus grand paradoxe de l'évolution. En effet, comment expliquer l'évolution de comportements qui tendent à augmenter les chances de survie et le succès reproductif d'autres individus, alors que ces actes diminuent l'aptitude de leurs auteurs ? La théorie de la sélection de parentèle permet de résoudre ce problème. Cette théorie stipule qu'en aidant de proches parents à se reproduire, les individus peuvent promouvoir indirectement la transmission de copies de leurs propres gènes à la génération suivante. Chez les insectes sociaux, l'altruisme des ouvrières s'explique par la théorie de sélection de parentèle lorsque les colonies sont monogynes (constituées d'une seule reine) puisque les ouvrières sont fortement apparentées aux couvains dont elles s'occupent. Par contre, les espèces dites unicoloniales, dont les colonies forment des réseaux de nids appelés supercolonies, représentent toujours un paradoxe pour les théories de l'évolution puisque l'apparentement entre les différents individus d'un nid est nulle. De plus, l'unicolonialité ne devrait pas être stable sur le long terme parce que les ouvrières qui ne s'occupent plus de leur apparentés devraient développer des stratégies plus égoïstes au cours du temps. La fourmi d'Argentine (Linepithema humile) est une espèce invasive ayant un impact considérable sur son environnement. Originaire d'Amérique du Sud, elle a été introduite dans pratiquement toutes les régions du monde dont le climat est de type méditerranéen. Son incroyable succès invasif s'explique par sa structure sociale unicoloniale observée dans chacun des pays où elle a été introduite. Par contre, les rares études effectuées en Argentine ont suggéré que la fourmi d'Argentine n'était pas unicoloniale dans son aire native. L'unicolonialité chez la fourmi d'Argentine était donc considéré comme une conséquence de son introduction dans de nouveaux environnements. Durant cette thèse, nous avons vérifié si la structure sociale de cette espèce différait fondamentalement entre l'aire native et introduite. Pour cela, nous avons étudié, à différentes échelles géographiques, des populations introduites et argentines avec une approche génétique et comportementale. L'ensemble de nos résultats montrent que les différences entre les deux structure sociales ne sont pas aussi importantes que ce que l'on imaginait. Les populations natives sont aussi constituées de réseaux de nids coopérants. La taille de ses supercolonies est toutefois bien moins importante en Argentine et il n'est pas rare de trouver plusieurs supercolonies cohabitantes dans une même population. Nous avons démontré que ces réseaux de nids étaient constitués d'individus qui sont plus apparentés entre eux qu'ils ne le sont avec les individus d'autres supercolonies, ainsi l'unicolonialité dans son aire d'origine ne représente pas un réel paradoxe pour les théories de l'évolution. Finalement nous avons étudié la même population en Argentine à six ans d'intervalle et avons constaté que les supercolonies avaient un taux de survie très faible ce qui pourrait expliquer la stabilité de l'unicolonialité au cours du temps. Si les supercolonies les plus compétitives survivent mieux que les supercolonies dans lesquelles apparaissent des traits égoïstes, on devrait alors observer une dynamique entre l'apparition de traits égoïstes et l'élimination des supercolonies dans lesquelles ces traits égoïstes évolueraient. Finalement, une étude mondiale nous a montré que les supercolonies étaient originaires de nombreux événements d'introductions indépendants. En conclusion, le succès invasif de la fourmi d'Argentine n'est donc pas dû à un changement de comportement associé à son introduction mais est lié aux caractéristiques qu'elle a développées en Argentine.

Natural gene-expression variation in Down syndrome modulates the outcome of gene-dosage imbalance.

Down syndrome (DS) is characterized by extensive phenotypic variability, with most traits occurring in only a fraction of affected individuals. Substantial gene-expression variation is present among unaffected individuals, and this variation has a strong genetic component. Since DS is caused by genomic-dosage imbalance, we hypothesize that gene-expression variation of human chromosome 21 (HSA21) genes in individuals with DS has an impact on the phenotypic variability among affected individuals. We studied gene-expression variation in 14 lymphoblastoid and 17 fibroblast cell lines from individuals with DS and an equal number of controls. Gene expression was assayed using quantitative real-time polymerase chain reaction on 100 and 106 HSA21 genes and 23 and 26 non-HSA21 genes in lymphoblastoid and fibroblast cell lines, respectively. Surprisingly, only 39% and 62% of HSA21 genes in lymphoblastoid and fibroblast cells, respectively, showed a statistically significant difference between DS and normal samples, although the average up-regulation of HSA21 genes was close to the expected 1.5-fold in both cell types. Gene-expression variation in DS and normal samples was evaluated using the Kolmogorov-Smirnov test. According to the degree of overlap in expression levels, we classified all genes into 3 groups: (A) nonoverlapping, (B) partially overlapping, and (C) extensively overlapping expression distributions between normal and DS samples. We hypothesize that, in each cell type, group A genes are the most dosage sensitive and are most likely involved in the constant DS traits, group B genes might be involved in variable DS traits, and group C genes are not dosage sensitive and are least likely to participate in DS pathological phenotypes. This study provides the first extensive data set on HSA21 gene-expression variation in DS and underscores its role in modulating the outcome of gene-dosage imbalance.

Expression and functional role of the prorenin receptor in the human adrenocortical zona glomerulosa and in primary aldosteronism.

OBJECTIVES: Prorenin can be detected in plasma of hypertensive patients. If detected in patients with primary aldosteronism could implicate prorenin in the development of primary aldosteronism. To address this issue, we measured the plasma prorenin levels in primary aldosteronism patients, the expression of the prorenin receptor (PRR) in the normal human adrenocortical zona glomerulosa and aldosterone-producing adenoma (APA), and we investigated the functional effects of PRR activation in human adrenocortical cells. METHOD: Plasma renin activity, aldosterone, and active and total trypsin-activated renin were measured in primary aldosteronism patients, essential hypertensive patients, and healthy individuals, and then prorenin levels were calculated. Localization and functional role of PRR were investigated in human and rat tissues, and aldosterone-producing cells. RESULTS: Primary aldosteronism patients had detectable plasma levels of prorenin. Using digital-droplet real-time PCR, we found a high PRR-to-porphobilinogen deaminase ratio in both the normal adrenal cortex and APAs. Marked expression of the PRR gene and protein was also found in HAC15 cells. Immunoblotting, confocal, and immunogold electron microscopy demonstrated PRR at the cell membrane and intracellularly. Renin and prorenin significantly triggered both CYP11B2 expression (aldosterone synthase) and ERK1/2 phosphorylation, but only CYP11B2 transcription was prevented by aliskiren. CONCLUSION: The presence of detectable plasma prorenin in primary aldosteronism patients, and the high expression of PRR in the normal human adrenal cortex, APA tissue, CD56+ aldosterone-producing cells, along with activation of CYP11B2 synthesis and ERK1/2 phosphorylation, suggest that the circulating and locally produced prorenin may contribute to the development or maintenance of human primary aldosteronism.

Probabilistic electrical resistivity tomography of a CO2 sequestration analog

Electrical resistivity tomography (ERT) is a well-established method for geophysical characterization and has shown potential for monitoring geologic CO2 sequestration, due to its sensitivity to electrical resistivity contrasts generated by liquid/gas saturation variability. In contrast to deterministic inversion approaches, probabilistic inversion provides the full posterior probability density function of the saturation field and accounts for the uncertainties inherent in the petrophysical parameters relating the resistivity to saturation. In this study, the data are from benchtop ERT experiments conducted during gas injection into a quasi-2D brine-saturated sand chamber with a packing that mimics a simple anticlinal geological reservoir. The saturation fields are estimated by Markov chain Monte Carlo inversion of the measured data and compared to independent saturation measurements from light transmission through the chamber. Different model parameterizations are evaluated in terms of the recovered saturation and petrophysical parameter values. The saturation field is parameterized (1) in Cartesian coordinates, (2) by means of its discrete cosine transform coefficients, and (3) by fixed saturation values in structural elements whose shape and location is assumed known or represented by an arbitrary Gaussian Bell structure. Results show that the estimated saturation fields are in overall agreement with saturations measured by light transmission, but differ strongly in terms of parameter estimates, parameter uncertainties and computational intensity. Discretization in the frequency domain (as in the discrete cosine transform parameterization) provides more accurate models at a lower computational cost compared to spatially discretized (Cartesian) models. A priori knowledge about the expected geologic structures allows for non-discretized model descriptions with markedly reduced degrees of freedom. Constraining the solutions to the known injected gas volume improved estimates of saturation and parameter values of the petrophysical relationship. (C) 2014 Elsevier B.V. All rights reserved.

Two-dimensional probabilistic inversion of plane-wave electromagnetic data: Methodology, model constraints and joint inversion with electrical resistivity data

Probabilistic inversion methods based on Markov chain Monte Carlo (MCMC) simulation are well suited to quantify parameter and model uncertainty of nonlinear inverse problems. Yet, application of such methods to CPU-intensive forward models can be a daunting task, particularly if the parameter space is high dimensional. Here, we present a 2-D pixel-based MCMC inversion of plane-wave electromagnetic (EM) data. Using synthetic data, we investigate how model parameter uncertainty depends on model structure constraints using different norms of the likelihood function and the model constraints, and study the added benefits of joint inversion of EM and electrical resistivity tomography (ERT) data. Our results demonstrate that model structure constraints are necessary to stabilize the MCMC inversion results of a highly discretized model. These constraints decrease model parameter uncertainty and facilitate model interpretation. A drawback is that these constraints may lead to posterior distributions that do not fully include the true underlying model, because some of its features exhibit a low sensitivity to the EM data, and hence are difficult to resolve. This problem can be partly mitigated if the plane-wave EM data is augmented with ERT observations. The hierarchical Bayesian inverse formulation introduced and used herein is able to successfully recover the probabilistic properties of the measurement data errors and a model regularization weight. Application of the proposed inversion methodology to field data from an aquifer demonstrates that the posterior mean model realization is very similar to that derived from a deterministic inversion with similar model constraints.

Cdk5 controls lymphatic vessel development and function by phosphorylation of Foxc2.

The lymphatic system maintains tissue fluid balance, and dysfunction of lymphatic vessels and valves causes human lymphedema syndromes. Yet, our knowledge of the molecular mechanisms underlying lymphatic vessel development is still limited. Here, we show that cyclin-dependent kinase 5 (Cdk5) is an essential regulator of lymphatic vessel development. Endothelial-specific Cdk5 knockdown causes congenital lymphatic dysfunction and lymphedema due to defective lymphatic vessel patterning and valve formation. We identify the transcription factor Foxc2 as a key substrate of Cdk5 in the lymphatic vasculature, mechanistically linking Cdk5 to lymphatic development and valve morphogenesis. Collectively, our findings show that Cdk5-Foxc2 interaction represents a critical regulator of lymphatic vessel development and the transcriptional network underlying lymphatic vascular remodeling.

Error models in hydrogeology applications

Notre consommation en eau souterraine, en particulier comme eau potable ou pour l'irrigation, a considérablement augmenté au cours des années. De nombreux problèmes font alors leur apparition, allant de la prospection de nouvelles ressources à la remédiation des aquifères pollués. Indépendamment du problème hydrogéologique considéré, le principal défi reste la caractérisation des propriétés du sous-sol. Une approche stochastique est alors nécessaire afin de représenter cette incertitude en considérant de multiples scénarios géologiques et en générant un grand nombre de réalisations géostatistiques. Nous rencontrons alors la principale limitation de ces approches qui est le coût de calcul dû à la simulation des processus d'écoulements complexes pour chacune de ces réalisations. Dans la première partie de la thèse, ce problème est investigué dans le contexte de propagation de l'incertitude, oú un ensemble de réalisations est identifié comme représentant les propriétés du sous-sol. Afin de propager cette incertitude à la quantité d'intérêt tout en limitant le coût de calcul, les méthodes actuelles font appel à des modèles d'écoulement approximés. Cela permet l'identification d'un sous-ensemble de réalisations représentant la variabilité de l'ensemble initial. Le modèle complexe d'écoulement est alors évalué uniquement pour ce sousensemble, et, sur la base de ces réponses complexes, l'inférence est faite. Notre objectif est d'améliorer la performance de cette approche en utilisant toute l'information à disposition. Pour cela, le sous-ensemble de réponses approximées et exactes est utilisé afin de construire un modèle d'erreur, qui sert ensuite à corriger le reste des réponses approximées et prédire la réponse du modèle complexe. Cette méthode permet de maximiser l'utilisation de l'information à disposition sans augmentation perceptible du temps de calcul. La propagation de l'incertitude est alors plus précise et plus robuste. La stratégie explorée dans le premier chapitre consiste à apprendre d'un sous-ensemble de réalisations la relation entre les modèles d'écoulement approximé et complexe. Dans la seconde partie de la thèse, cette méthodologie est formalisée mathématiquement en introduisant un modèle de régression entre les réponses fonctionnelles. Comme ce problème est mal posé, il est nécessaire d'en réduire la dimensionnalité. Dans cette optique, l'innovation du travail présenté provient de l'utilisation de l'analyse en composantes principales fonctionnelles (ACPF), qui non seulement effectue la réduction de dimensionnalités tout en maximisant l'information retenue, mais permet aussi de diagnostiquer la qualité du modèle d'erreur dans cet espace fonctionnel. La méthodologie proposée est appliquée à un problème de pollution par une phase liquide nonaqueuse et les résultats obtenus montrent que le modèle d'erreur permet une forte réduction du temps de calcul tout en estimant correctement l'incertitude. De plus, pour chaque réponse approximée, une prédiction de la réponse complexe est fournie par le modèle d'erreur. Le concept de modèle d'erreur fonctionnel est donc pertinent pour la propagation de l'incertitude, mais aussi pour les problèmes d'inférence bayésienne. Les méthodes de Monte Carlo par chaîne de Markov (MCMC) sont les algorithmes les plus communément utilisés afin de générer des réalisations géostatistiques en accord avec les observations. Cependant, ces méthodes souffrent d'un taux d'acceptation très bas pour les problèmes de grande dimensionnalité, résultant en un grand nombre de simulations d'écoulement gaspillées. Une approche en deux temps, le "MCMC en deux étapes", a été introduite afin d'éviter les simulations du modèle complexe inutiles par une évaluation préliminaire de la réalisation. Dans la troisième partie de la thèse, le modèle d'écoulement approximé couplé à un modèle d'erreur sert d'évaluation préliminaire pour le "MCMC en deux étapes". Nous démontrons une augmentation du taux d'acceptation par un facteur de 1.5 à 3 en comparaison avec une implémentation classique de MCMC. Une question reste sans réponse : comment choisir la taille de l'ensemble d'entrainement et comment identifier les réalisations permettant d'optimiser la construction du modèle d'erreur. Cela requiert une stratégie itérative afin que, à chaque nouvelle simulation d'écoulement, le modèle d'erreur soit amélioré en incorporant les nouvelles informations. Ceci est développé dans la quatrième partie de la thèse, oú cette méthodologie est appliquée à un problème d'intrusion saline dans un aquifère côtier. -- Our consumption of groundwater, in particular as drinking water and for irrigation, has considerably increased over the years and groundwater is becoming an increasingly scarce and endangered resource. Nofadays, we are facing many problems ranging from water prospection to sustainable management and remediation of polluted aquifers. Independently of the hydrogeological problem, the main challenge remains dealing with the incomplete knofledge of the underground properties. Stochastic approaches have been developed to represent this uncertainty by considering multiple geological scenarios and generating a large number of realizations. The main limitation of this approach is the computational cost associated with performing complex of simulations in each realization. In the first part of the thesis, we explore this issue in the context of uncertainty propagation, where an ensemble of geostatistical realizations is identified as representative of the subsurface uncertainty. To propagate this lack of knofledge to the quantity of interest (e.g., the concentration of pollutant in extracted water), it is necessary to evaluate the of response of each realization. Due to computational constraints, state-of-the-art methods make use of approximate of simulation, to identify a subset of realizations that represents the variability of the ensemble. The complex and computationally heavy of model is then run for this subset based on which inference is made. Our objective is to increase the performance of this approach by using all of the available information and not solely the subset of exact responses. Two error models are proposed to correct the approximate responses follofing a machine learning approach. For the subset identified by a classical approach (here the distance kernel method) both the approximate and the exact responses are knofn. This information is used to construct an error model and correct the ensemble of approximate responses to predict the "expected" responses of the exact model. The proposed methodology makes use of all the available information without perceptible additional computational costs and leads to an increase in accuracy and robustness of the uncertainty propagation. The strategy explored in the first chapter consists in learning from a subset of realizations the relationship between proxy and exact curves. In the second part of this thesis, the strategy is formalized in a rigorous mathematical framework by defining a regression model between functions. As this problem is ill-posed, it is necessary to reduce its dimensionality. The novelty of the work comes from the use of functional principal component analysis (FPCA), which not only performs the dimensionality reduction while maximizing the retained information, but also allofs a diagnostic of the quality of the error model in the functional space. The proposed methodology is applied to a pollution problem by a non-aqueous phase-liquid. The error model allofs a strong reduction of the computational cost while providing a good estimate of the uncertainty. The individual correction of the proxy response by the error model leads to an excellent prediction of the exact response, opening the door to many applications. The concept of functional error model is useful not only in the context of uncertainty propagation, but also, and maybe even more so, to perform Bayesian inference. Monte Carlo Markov Chain (MCMC) algorithms are the most common choice to ensure that the generated realizations are sampled in accordance with the observations. Hofever, this approach suffers from lof acceptance rate in high dimensional problems, resulting in a large number of wasted of simulations. This led to the introduction of two-stage MCMC, where the computational cost is decreased by avoiding unnecessary simulation of the exact of thanks to a preliminary evaluation of the proposal. In the third part of the thesis, a proxy is coupled to an error model to provide an approximate response for the two-stage MCMC set-up. We demonstrate an increase in acceptance rate by a factor three with respect to one-stage MCMC results. An open question remains: hof do we choose the size of the learning set and identify the realizations to optimize the construction of the error model. This requires devising an iterative strategy to construct the error model, such that, as new of simulations are performed, the error model is iteratively improved by incorporating the new information. This is discussed in the fourth part of the thesis, in which we apply this methodology to a problem of saline intrusion in a coastal aquifer.