Conditioning of stochastic 3-D fracture networks to hydrological and geophysical data

The geometry and connectivity of fractures exert a strong influence on the flow and transport properties of fracture networks. We present a novel approach to stochastically generate three-dimensional discrete networks of connected fractures that are conditioned to hydrological and geophysical data. A hierarchical rejection sampling algorithm is used to draw realizations from the posterior probability density function at different conditioning levels. The method is applied to a well-studied granitic formation using data acquired within two boreholes located 6 m apart. The prior models include 27 fractures with their geometry (position and orientation) bounded by information derived from single-hole ground-penetrating radar (GPR) data acquired during saline tracer tests and optical televiewer logs. Eleven cross-hole hydraulic connections between fractures in neighboring boreholes and the order in which the tracer arrives at different fractures are used for conditioning. Furthermore, the networks are conditioned to the observed relative hydraulic importance of the different hydraulic connections by numerically simulating the flow response. Among the conditioning data considered, constraints on the relative flow contributions were the most effective in determining the variability among the network realizations. Nevertheless, we find that the posterior model space is strongly determined by the imposed prior bounds. Strong prior bounds were derived from GPR measurements and helped to make the approach computationally feasible. We analyze a set of 230 posterior realizations that reproduce all data given their uncertainties assuming the same uniform transmissivity in all fractures. The posterior models provide valuable statistics on length scales and density of connected fractures, as well as their connectivity. In an additional analysis, effective transmissivity estimates of the posterior realizations indicate a strong influence of the DFN structure, in that it induces large variations of equivalent transmissivities between realizations. The transmissivity estimates agree well with previous estimates at the site based on pumping, flowmeter and temperature data.

Properties of functional brain networks correlate with frequency of psychogenic non-epileptic seizures.

Abnormalities in the topology of brain networks may be an important feature and etiological factor for psychogenic non-epileptic seizures (PNES). To explore this possibility, we applied a graph theoretical approach to functional networks based on resting state EEGs from 13 PNES patients and 13 age- and gender-matched controls. The networks were extracted from Laplacian-transformed time-series by a cross-correlation method. PNES patients showed close to normal local and global connectivity and small-world structure, estimated with clustering coefficient, modularity, global efficiency, and small-worldness (SW) metrics, respectively. Yet the number of PNES attacks per month correlated with a weakness of local connectedness and a skewed balance between local and global connectedness quantified with SW, all in EEG alpha band. In beta band, patients demonstrated above-normal resiliency, measured with assortativity coefficient, which also correlated with the frequency of PNES attacks. This interictal EEG phenotype may help improve differentiation between PNES and epilepsy. The results also suggest that local connectivity could be a target for therapeutic interventions in PNES. Selective modulation (strengthening) of local connectivity might improve the skewed balance between local and global connectivity and so prevent PNES events.

Nuclear receptor Rev-erb alpha (Nr1d1) functions in concert with Nr2e3 to regulate transcriptional networks in the retina.

The majority of diseases in the retina are caused by genetic mutations affecting the development and function of photoreceptor cells. The transcriptional networks directing these processes are regulated by genes such as nuclear hormone receptors. The nuclear hormone receptor gene Rev-erb alpha/Nr1d1 has been widely studied for its role in the circadian cycle and cell metabolism, however its role in the retina is unknown. In order to understand the role of Rev-erb alpha/Nr1d1 in the retina, we evaluated the effects of loss of Nr1d1 to the developing retina and its co-regulation with the photoreceptor-specific nuclear receptor gene Nr2e3 in the developing and mature retina. Knock-down of Nr1d1 expression in the developing retina results in pan-retinal spotting and reduced retinal function by electroretinogram. Our studies show that NR1D1 protein is co-expressed with NR2E3 in the outer neuroblastic layer of the developing mouse retina. In the adult retina, NR1D1 is expressed in the ganglion cell layer and is co-expressed with NR2E3 in the outer nuclear layer, within rods and cones. Several genes co-targeted by NR2E3 and NR1D1 were identified that include: Nr2c1, Recoverin, Rgr, Rarres2, Pde8a, and Nupr1. We examined the cyclic expression of Nr1d1 and Nr2e3 over a twenty-four hour period and observed that both nuclear receptors cycle in a similar manner. Taken together, these studies reveal a novel role for Nr1d1, in conjunction with its cofactor Nr2e3, in regulating transcriptional networks critical for photoreceptor development and function.

Using active learning for monitoring networks design: The example of wind power plants sites evaluation

Locating new wind farms is of crucial importance for energy policies of the next decade. To select the new location, an accurate picture of the wind fields is necessary. However, characterizing wind fields is a difficult task, since the phenomenon is highly nonlinear and related to complex topographical features. In this paper, we propose both a nonparametric model to estimate wind speed at different time instants and a procedure to discover underrepresented topographic conditions, where new measuring stations could be added. Compared to space filling techniques, this last approach privileges optimization of the output space, thus locating new potential measuring sites through the uncertainty of the model itself.

Modeling resting-state functional networks when the cortex falls asleep: local and global changes.

The transition from wakefulness to sleep represents the most conspicuous change in behavior and the level of consciousness occurring in the healthy brain. It is accompanied by similarly conspicuous changes in neural dynamics, traditionally exemplified by the change from "desynchronized" electroencephalogram activity in wake to globally synchronized slow wave activity of early sleep. However, unit and local field recordings indicate that the transition is more gradual than it might appear: On one hand, local slow waves already appear during wake; on the other hand, slow sleep waves are only rarely global. Studies with functional magnetic resonance imaging also reveal changes in resting-state functional connectivity (FC) between wake and slow wave sleep. However, it remains unclear how resting-state networks may change during this transition period. Here, we employ large-scale modeling of the human cortico-cortical anatomical connectivity to evaluate changes in resting-state FC when the model "falls asleep" due to the progressive decrease in arousal-promoting neuromodulation. When cholinergic neuromodulation is parametrically decreased, local slow waves appear, while the overall organization of resting-state networks does not change. Furthermore, we show that these local slow waves are structured macroscopically in networks that resemble the resting-state networks. In contrast, when the neuromodulator decrease further to very low levels, slow waves become global and resting-state networks merge into a single undifferentiated, broadly synchronized network.

HIV-1 drug resistance transmission networks in southwest Switzerland.

To determine viral subtypes and resistance mutations to antiretroviral treatment (ART) in untreated HIV-1 acutely infected subjects from Southwest Switzerland. Clinical samples were obtained from the HIV primary infection cohort from Lausanne. Briefly, pol gene was amplified by nested PCR and sequenced to generate a 1?kb sequence spanning protease and reverse transcriptase key protein regions. Nucleotide sequences were used to assess viral genotype and ART resistance mutations. Blood specimens and medical information were obtained from 30 patients. Main viral subtypes corresponded to clade B, CRF02_AG, and F1. Resistant mutations to PIs consisted of L10V and accessory mutations 16E and 60E present in all F1 clades. The NNRTI major resistant mutation 103N was detected in all F1 viruses and in other 2 clades. Additionally, we identified F1 sequences from other 6 HIV infected and untreated individuals from Southwest Switzerland, harboring nucleotide motifs and resistance mutations to ART as observed in the F1 strains from the cohort. These data reveal a high transmission rate (16.6%) for NNRTI resistant mutation 103N in a cohort of HIV acute infection. Three of the 5 resistant strains were F1 clades closely related to other F1 isolates from HIV-1 infection untreated patients also coming from Southwest Switzerland. Overall, we provide strong evidence towards an HIV-1 resistant transmission network in Southwest Switzerland. These findings have relevant implications for the local molecular mapping of HIV-1 and future ART surveillance studies in the region.

Evidence for transcript networks composed of chimeric RNAs in human cells.

The classic organization of a gene structure has followed the Jacob and Monod bacterial gene model proposed more than 50 years ago. Since then, empirical determinations of the complexity of the transcriptomes found in yeast to human has blurred the definition and physical boundaries of genes. Using multiple analysis approaches we have characterized individual gene boundaries mapping on human chromosomes 21 and 22. Analyses of the locations of the 5' and 3' transcriptional termini of 492 protein coding genes revealed that for 85% of these genes the boundaries extend beyond the current annotated termini, most often connecting with exons of transcripts from other well annotated genes. The biological and evolutionary importance of these chimeric transcripts is underscored by (1) the non-random interconnections of genes involved, (2) the greater phylogenetic depth of the genes involved in many chimeric interactions, (3) the coordination of the expression of connected genes and (4) the close in vivo and three dimensional proximity of the genomic regions being transcribed and contributing to parts of the chimeric RNAs. The non-random nature of the connection of the genes involved suggest that chimeric transcripts should not be studied in isolation, but together, as an RNA network.

EEG-based functional brain networks: does the network size matter?

Functional connectivity in human brain can be represented as a network using electroencephalography (EEG) signals. These networks--whose nodes can vary from tens to hundreds--are characterized by neurobiologically meaningful graph theory metrics. This study investigates the degree to which various graph metrics depend upon the network size. To this end, EEGs from 32 normal subjects were recorded and functional networks of three different sizes were extracted. A state-space based method was used to calculate cross-correlation matrices between different brain regions. These correlation matrices were used to construct binary adjacency connectomes, which were assessed with regards to a number of graph metrics such as clustering coefficient, modularity, efficiency, economic efficiency, and assortativity. We showed that the estimates of these metrics significantly differ depending on the network size. Larger networks had higher efficiency, higher assortativity and lower modularity compared to those with smaller size and the same density. These findings indicate that the network size should be considered in any comparison of networks across studies.