Homology Mapping Research Articles

Lichen Biosynthetic Gene Clusters Part II: Homology Mapping Suggests a Functional Diversity

Lichens are renowned for their diverse natural products though little is known of the genetic programming dictating lichen natural product biosynthesis. We sequenced the genome of Cladonia uncialis and profiled its secondary metabolite biosynthetic gene clusters. Through a homology searching approach, we can now propose specific functions for gene products as well as the biosynthetic pathways that are encoded in several of these gene clusters. This analysis revealed that the lichen genome encodes the required enzymes for patulin and betaenones A-C biosynthesis, fungal toxins not known to be produced by lichens. Within several gene clusters, some (but not all) genes are genetically similar to genes devoted to secondary metabolite biosynthesis in Fungi. These lichen clusters also contain accessory tailoring genes without such genetic similarity, suggesting that the encoded tailoring enzymes perform distinct chemical transformations. We hypothesize that C. uncialis gene clusters have evolved by shuffling components of ancestral fungal clusters to create new series of chemical steps, leading to the production of hitherto undiscovered derivatives of fungal secondary metabolites.

Using genomic data to study insecticide resistance in the house fly, Musca domestica

The house fly, Musca domestica, is a major pest at livestock facilities throughout the world. Insecticides have been the most common control strategy for flies, but many populations have evolved resistance. The speed by which we are able to identify the mutations responsible for resistance has been a major challenge for the development of high throughput resistance monitoring assays as new insecticides are introduced for control. This is particularly true for mutations that cause trans regulation of a gene, which then results in resistance. In this paper we take advantage of the conserved homology of dipteran chromosomes to assign 3069 genes to chromosomes. Of these, 234 were of toxicological interest (CYPs, esterases/hydrolases, glutathione S-transferases (GSTs) and target sites). The chromosomal location of genes known from linkage analysis studies matched the location predicted by homology mapping in ten out of ten cases, indicating a high reliability of our approach. The CYPs, esterases/hydrolases and GSTs were not randomly distributed throughout the genome. They clustered on chromosomes, but the pattern was different between the CYPs, esterases/hydrolases and GSTs. Examples are provided for how the availability of the house fly genome, combined with an ability to assign genes to chromosomes, will help to accelerate research in house flies.

Genome-Wide Prediction of Potential Vaccine Candidates for Campylobacter jejuni Using Reverse Vaccinology.

Campylobacteriosis is a deadly disease which has developed resistance to most of the available chemotherapeutic agents. Although various studies provide evidence of acquired immunity following exposure to Campylobacter jejuni, no effective vaccine has been developed, still. Hence, there is an urgent need to identify potential vaccine candidates for Campylobacter species. In the proposed study, Campylobacter jejuni subsp. jejuni serotype O:2 (strain NCTC 11168) was taken and computational approach was employed to screen C. jejuni genome for promising vaccine candidates. From 1623 protein-coding sequences, 37 potential antigens were screened for epitope prediction based on surface association, consensus antigenicity predictions, solubility, transmembrane domain, and ortholog analysis. Comprehensive immunogenic analysis of these 37 antigens revealed that antigen Q0PA22 shows the greatest potential for experimental immunogenicity analysis. It has several potential CD4+and CD8+T-cell epitopes, as well as high probability of B-cell epitope regions as compared to well-characterized antigen Omp18 (Uniprot ID:Q0PC24). Among the highest scoring predicted epitopes, an optimal set of epitopes with respect to overall immunogenicity in target populations for campylobacteriosis viz. Europe, North America and Southwest Asia was determined. An epitope AMLTYMQWL from antigen no. 6(Q0PA22) binds to the most prevalent allele HLA-A*0201, and this epitope has most immunogenicity for all the target populations. In addition, this epitope exhibited highly significant TCR-pMHC interactions having a joint Z value of 4.87. Homology mapping studies of the predicted epitope show best homology to a well-studied antigenic peptide from influenza virus H5N1. Therefore, the predicted epitope might be a suitable vaccine candidate.

Quasistabilization and basepoint moving maps in link Floer homology

We analyze the effect of adding, removing, and moving basepoints on link Floer homology. We prove that adding or removing basepoints via a procedure called quasi-stabilization is a natural operation on a certain version of link Floer homology, which we call $CFL_{UV}^\infty$. We consider the effect on the full link Floer complex of moving basepoints, and develop a simple calculus for moving basepoints on the link Floer complexes. We apply it to compute the effect of several diffeomorphisms corresponding to moving basepoints. Using these techniques we prove a conjecture of Sarkar about the map on the full link Floer complex induced by a finger move along a link component.

Saliency detection based on salient edges and remarkable discriminating for superpixel pairs

Saliency detection is an essential pre-processing step for intensifying image objects in many computer vision fields. Other than most bottom-up methods, in this paper, we propose a novel saliency model based on high-level image edges and low-level feature contrast. Edges are inherent features for an image, and it can locate the object via a salient selection. With this theory, an accurate object contour is generated after two salient discriminations for edge-contiguous superpixel couples. After position is confirmed, foreground and background of the image can be divided by the contour. With an ingrowth model, we then obtain a foreground (or background) seeds referencing spatial adjacent relation. According to the seed set, a homologous saliency map is computed via a seed-based saliency approach, which is proposed on the basis of affinity matrix. Compared with some pre-existing algorithms, low-level information is utilized purposefully by salient edges in the presented method, which makes the extracted fore regions more precise.

A Comparative Study of Pygmy Hippopotamus (Choeropsis liberiensis) Karyotype by Cross-Species Chromosome Painting

The family Hippopotamidae is comprised of two genera with two living species, the common hippo (Hippopotamus amphibius) and the pygmy hippo (Choeropsis liberiensis). Unlike the common hippo, the karyotype of C. liberiensis has not yet been investigated via cross-species chromosome painting methods. We established chromosomal homologies between the pygmy hippo, pig, and cattle by fluorescence in situ hybridization using whole chromosome, arm-specific, region specific, and bacterial artificial chromosome (BAC) probes. Probes from the 18 pig autosomes painted 45 conserved chromosomal segments in the pygmy hippo genome. The pygmy hippo and cattle homology map was deduced from our hybridization results of painting probes to pygmy hippo chromosomes with a combination of previously published dromedary hybridization data. On the pygmy hippo and cattle homology map, 29 cattle autosomes revealed 39 conservative segments on pygmy hippo chromosomes. For a more detailed structural analysis of genome rearrangements and X chromosome structure, we used cattle region specific and BAC probes. Our report demonstrates that cattle probes are useful not only in comparative studies within Ruminantia, but also in more phylogenetically distant Artiodactyla species.

XGSA: A statistical method for cross-species gene set analysis.

Gene set analysis is a powerful tool for determining whether an experimentally derived set of genes is statistically significantly enriched for genes in other pre-defined gene sets, such as known pathways, gene ontology terms, or other experimentally derived gene sets. Current gene set analysis methods do not facilitate comparing gene sets across different organisms as they do not explicitly deal with homology mapping between species. There lacks a systematic investigation about the effect of complex gene homology on cross-species gene set analysis. In this study, we show that not accounting for the complex homology structure when comparing gene sets in two species can lead to false positive discoveries, especially when comparing gene sets that have complex gene homology relationships. To overcome this bias, we propose a straightforward statistical approach, called XGSA, that explicitly takes the cross-species homology mapping into consideration when doing gene set analysis. Simulation experiments confirm that XGSA can avoid false positive discoveries, while maintaining good statistical power compared to other ad hoc approaches for cross-species gene set analysis. We further demonstrate the effectiveness of XGSA with two real-life case studies that aim to discover conserved or species-specific molecular pathways involved in social challenge and vertebrate appendage regeneration. The R source code for XGSA is available under a GNU General Public License at http://github.com/VCCRI/XGSA CONTACT: jho@victorchang.edu.au.

A Common Signal Patch Drives AP-1 Protein-dependent Golgi Export of Inwardly Rectifying Potassium Channels

Nearly all members of the inwardly rectifying potassium (Kir) channel family share a cytoplasmic domain structure that serves as an unusual AP-1 clathrin adaptor-dependent Golgi export signal in one Kir channel, Kir2.1 (KCNJ2), raising the question whether Kir channels share a common Golgi export mechanism. Here we explore this idea, focusing on two structurally and functionally divergent Kir family members, Kir2.3 (KCNJ4) and Kir4.1/5.1 (KCNJ10/16), which have ∼50% amino identity. We found that Golgi export of both channels is blocked upon siRNA-mediated knockdown of the AP-1 γ subunit, as predicted for the common AP-1-dependent trafficking process. A comprehensive mutagenic analysis, guided by homology mapping in atomic resolution models of Kir2.1, Kir2.3, and Kir4.1/5.1, identified a common structure that serves as a recognition site for AP-1 binding and governs Golgi export. Larger than realized from previous studies with Kir2.1, the signal is created by a patch of residues distributed at the confluence of cytoplasmic N and C termini. The signal involves a stretch of hydrophobic residues from the C-terminal region that form a hydrophobic cleft, an adjacent cluster of basic residues within the N terminus, and a potential network of salt bridges that join the N- and C-terminal poles together. Because patch formation and AP-1 binding are dependent on proper folding of the cytoplasmic domains, the signal provides a common quality control mechanism at the Golgi for Kir channels. These findings identify a new proteostatic mechanism that couples protein folding of channels to forward trafficking in the secretory pathway.

Stein fillings of homology 3-spheres and mapping class groups

In this article, using combinatorial techniques of mapping class groups, we show that a Stein fillable integral homology 3-sphere supported by an open book with page a 4-holed sphere admits a unique Stein filling up to symplectic deformation. Furthermore, according to a property of deforming symplectic fillings of rational homology 3-spheres into strong symplectic fillings, we also show that a symplectically fillable integral homology 3-sphere supported by an open book with page a 4-holed sphere admits a unique symplectic filling up to symplectic deformation and blow-up.

Comparison of Stochastic and Deterministic Methods for Mapping Environmental Noise from Opencast Quarries

An accurate mapping of the acoustic spatial variability generated by quarrying plants constitutes an important step in the assessment of the environmental compatibility of the extraction site with the surrounding area. This paper refers to the results of a study comparing performances obtained by stochastic and deterministic methods to determine the acoustic climate of the areas surrounding opencast quarries. To this aim noise levels are monitored in the areas surrounding a limestone quarry site in central Italy. In the first phase an analytical prediction model, based on ISO 9613 and ISO 3746, is developed using the sampled sound power level of each source as input data. In the second phase, in contrast, the stochastic method of Ordinary Kriging is used to plot an acoustic map. This is drawn up by interpolating the scattered sound pressure level data from sampled points which are identified as nodes of a rectangular grid covering the chosen area. The analytical prediction model output, which consists in the sound pressure level at various target points, is then compared with the homologous map obtained from the Ordinary Kriging method. Finally the resulting maps are evaluated with respect to statistical measures and cross validation procedures in order to determine the most suitable method.

An Approach for Predicting Essential Genes Using Multiple Homology Mapping and Machine Learning Algorithms.

Investigation of essential genes is significant to comprehend the minimal gene sets of cell and discover potential drug targets. In this study, a novel approach based on multiple homology mapping and machine learning method was introduced to predict essential genes. We focused on 25 bacteria which have characterized essential genes. The predictions yielded the highest area under receiver operating characteristic (ROC) curve (AUC) of 0.9716 through tenfold cross-validation test. Proper features were utilized to construct models to make predictions in distantly related bacteria. The accuracy of predictions was evaluated via the consistency of predictions and known essential genes of target species. The highest AUC of 0.9552 and average AUC of 0.8314 were achieved when making predictions across organisms. An independent dataset from Synechococcus elongatus, which was released recently, was obtained for further assessment of the performance of our model. The AUC score of predictions is 0.7855, which is higher than other methods. This research presents that features obtained by homology mapping uniquely can achieve quite great or even better results than those integrated features. Meanwhile, the work indicates that machine learning-based method can assign more efficient weight coefficients than using empirical formula based on biological knowledge.

Recognition of ZnT8, Proinsulin, and Homologous MAP Peptides in Sardinian Children at Risk of T1D Precedes Detection of Classical Islet Antibodies.

As numerous studies put in evidence the increasing incidence of type 1 diabetes (T1D) in children, an early diagnosis is of great importance to define correct treatment and diet. Currently, the identification of classical islet autoantibodies is the primary biomarker for diagnosis in subjects at risk, especially in pediatric patients. Recent studies suggest that detection of antibodies against ZnT8 protein in preclinical phase can predict the development of T1D. We previously demonstrated a significant association of Mycobacterium avium subspecies paratuberculosis (MAP) with T1D in adult Sardinian patients. To enforce this finding, we investigated the presence of antibodies against ZnT8 and proinsulin (PI) with respective homologous epitopes: MAP3865c133–141/ZnT8186–194, MAP3865c125–133/ZnT8178–186, MAP2404c70–85/PI46–61, and MAP1,4αgbp157–173/PI64–80, in 23 children at risk for T1D, formerly involved in the TRIGR study, and 22 healthy controls (HCs). Positivity to anti-MAP and homologous human peptides was detected in 48% of at-risk subjects compared to 5,85% HCs, preceding appearance of islet autoantibodies. Being MAP easily transmitted to humans with infected cow's milk and detected in retail infant formulas, MAP epitopes could be present in extensively hydrolyzed formula and act as antigens stimulating β-cell autoimmunity.

Comparison of Expression and Epigenetic Profiles in Human and Mouse Erythropoiesis and Megakaryopoiesis Using a Systems Biology Model

To date numerous datasets of gene expression and epigenetic profiles for mouse and human hematopoietic cells have been generated. While individual data sets for a particular cell type have been correlated, no approach exists to harness all expression and epigenetic profiles for the different types of hematopoietic cells. Our goal is to develop a systems biology platform to compare epigenetic profiles of hematopoietic cells towards a better understanding of epigenetic mechanisms governing hematopoiesis. To provide the necessary foundation to support systematic studies of hematopoiesis, we have developed the Systems Biology Repository (SBR, http://sbrblood.nhgri.nih.gov), a data "ranch" for organizing and analyzing transcriptome and epigenome data cells throughout differentiation.To populate SBR, we extracted, curated, annotated, and integrated all human and mouse hematopoietic datasets available through the Encyclopedia of DNA Elements (ENCODE), the Gene Expression Omnibus (GEO) and the Short Read Repository (SRR). These include genome-wide profiles of DNA methylation, histone methylation and acetylation, transcription factor occupancy (ChIPSeq), chromatin accessibility (DNaseISeq, ATACSeq, FAIRESeq), and coding as well as non-coding transcriptional profiles (RNASeq).To demonstrate the utility of SBR, we conducted three different analyses. The first was a vertical study of HistoneSeq (H3K4me1, H3K4me2, H3K4me3, and H3K27ac), DNA methylation and RNASeq profiles during mouse erythroid differentiation. We found a global decrease in DNA methylation from hematopoietic stem and progenitor cells (HSC) through common myeloid progenitors (CMP), erythroid progenitor cells (MEP) and erythroblasts (ERY; 92936 peaks in HSC to 14422 in ERY). The number of expressed genes (using a tags per million cutoff of 10) increased in erythroid progenitors (8901 in HSC to 10778 in CMP and 10670 in MEP) before decreasing in ERY (8654). 62% of histone marks delineating active enhancers (H3K27ac, H3K4me1) are present in both HSC and ERY, while 48% arise de novo during differentiation. In contrast, only 16% of active promoter specific histone marks (H3K4me2, H3K4me3) are present in both HSC and ERY.For a horizontal analysis we compared the DNA methylation, RNASeq, histone modification (H3K4me1, H3K4me2, H3K4me3, and H3K27ac) and transcription factor binding (GATA1 and NFE2) profiles of erythroblasts (ERY) and megakaryocytes (MEG). We found a similar relationship between gene expression and the histone and DNA methylation profiles in each cell type but differences between expression and in transcription factor occupancy. DNA methylation and H3K4me3 was enriched in the gene body of expressed genes (>36%) for both ERY (p ≤ 0.001) and MEG (p ≤ 0.01). In contrast DNA methylation was enriched in the upstream and downstream regions of non-coding RNA genes (p ≤ 0.001). Transcription factor occupancy was cell type specific: 79% of GATA1 sites are in ERY and 72% of NFE2 sites are in MEG. In erythroblasts, DNA methylation and GATA1 binding in the gene body are associated with gene silencing (4 fold difference, p ≤ 0.001), while in megakaryocytes, DNA methylation and NFE2 binding in the gene body are associated with gene activation (8 fold difference, p ≤ 0.001).We used the Mouse Genome Informatics homology map data to perform a cross-species comparison of the expression profiles of mouse and human multipotent progenitors (MPP), proerythroblasts and orthochromatic erythroblasts. We found a total of 5247 genes expressed at significantly different levels (p ≤ 0.001) between human and mouse MPP, while only 2010 genes were expressed at significantly similar levels (p ≤ 0.001). At the proerythroblast and orthochromatic erythroblast stages 7696 genes and 6571 genes were expressed at significantly different levels (p ≤ 0.001) between human and mouse respectively, while 2024 and 2560 genes were expressed at significantly similar levels (p ≤ 0.001). These data are consistent with previous studies showing differences in the transcriptional profiles of mouse and human hematopoietic cells.In summary, SBR provides a foundation to model the genetic and epigenetic landscape in both the mouse and human hematopoietic system, and enables functional correlations to be made between the species. As SBR is expanded to include data from patient cells, it will be possible to model epigenetic changes associated with disease. DisclosuresNo relevant conflicts of interest to declare.

Genome-wide analysis of alternative splicing events in Hordeum vulgare: Highlighting retention of intron-based splicing and its possible function through network analysis

In this study, using homology mapping of assembled expressed sequence tags against the genomic data, we identified alternative splicing events in barley. Results demonstrated that intron retention is frequently associated with specific abiotic stresses. Network analysis resulted in discovery of some specific sub-networks between miRNAs and transcription factors in genes with high number of alternative splicing, such as cross talk between SPL2, SPL10 and SPL11 regulated by miR156 and miR157 families. To confirm the alternative splicing events, elongation factor protein (MLOC_3412) was selected followed by experimental verification of the predicted splice variants by Semi quantitative Reverse Transcription PCR (qRT-PCR). Our novel integrative approach opens a new avenue for functional annotation of alternative splicing through regulatory-based network discovery.

MitoCarta2.0: an updated inventory of mammalian mitochondrial proteins

Mitochondria are complex organelles that house essential pathways involved in energy metabolism, ion homeostasis, signalling and apoptosis. To understand mitochondrial pathways in health and disease, it is crucial to have an accurate inventory of the organelle's protein components. In 2008, we made substantial progress toward this goal by performing in-depth mass spectrometry of mitochondria from 14 organs, epitope tagging/microscopy and Bayesian integration to assemble MitoCarta (www.broadinstitute.org/pubs/MitoCarta): an inventory of genes encoding mitochondrial-localized proteins and their expression across 14 mouse tissues. Using the same strategy we have now reconstructed this inventory separately for human and for mouse based on (i) improved gene transcript models, (ii) updated literature curation, including results from proteomic analyses of mitochondrial sub-compartments, (iii) improved homology mapping and (iv) updated versions of all seven original data sets. The updated human MitoCarta2.0 consists of 1158 human genes, including 918 genes in the original inventory as well as 240 additional genes. The updated mouse MitoCarta2.0 consists of 1158 genes, including 967 genes in the original inventory plus 191 additional genes. The improved MitoCarta 2.0 inventory provides a molecular framework for system-level analysis of mammalian mitochondria.

On a Hodge theoretic property of the Künneth map in periodic cyclic homology

We prove that the standard Künneth map in periodic cyclic homology of differential Z/2-graded algebras is compatible with a generalization of the Hodge filtration and explain how this result is related to various Thom–Sebastiani type theorems in singularity theory.

MitoMiner v3.1, an update on the mitochondrial proteomics database.

Mitochondrial proteins remain the subject of intense research interest due to their implication in an increasing number of different conditions including mitochondrial and metabolic disease, cancer, and neuromuscular degenerative and age-related disorders. However, the mitochondrial proteome has yet to be accurately and comprehensively defined, despite many studies. To support mitochondrial research, we developed MitoMiner (http://mitominer.mrc-mbu.cam.ac.uk), a freely accessible mitochondrial proteomics database. MitoMiner integrates different types of subcellular localisation evidence with protein information from public resources, and so provides a comprehensive central resource for data on mitochondrial protein localisation. Here we report important updates to the database including the addition of subcellular immunofluorescent staining results from the Human Protein Atlas, computational predictions of mitochondrial targeting sequences, and additional large-scale mass-spectrometry and GFP tagging data sets. This evidence is shared across the 12 species in MitoMiner (now including Schizosaccharomyces pombe) by homology mapping. MitoMiner provides multiple ways of querying the data including simple text searches, predefined queries and custom queries created using the interactive QueryBuilder. For remote programmatic access, API's are available for several programming languages. This combination of data and flexible querying makes MitoMiner a unique platform to investigate mitochondrial proteins, with application in mitochondrial research and prioritising candidate mitochondrial disease genes.

Multidirectional chromosome painting substantiates the occurrence of extensive genomic reshuffling within Accipitriformes.

BackgroundPrevious cross-species painting studies with probes from chicken (Gallus gallus) chromosomes 1–10 and a paint pool of nineteen microchromosomes have revealed that the drastic karyotypic reorganization in Accipitridae is due to extensive synteny disruptions and associations. However, the number of synteny association events and identities of microchromosomes involved in such synteny associations remain undefined, due to the lack of paint probes derived from individual chicken microchromosomes. Moreover, no genome-wide homology map between Accipitridae species and other avian species with atypical karyotype organization has been reported till now, and the karyotype evolution within Accipitriformes remains unclear.ResultsTo delineate the synteny-conserved segments in Accipitridae, a set of painting probes for the griffon vulture, Gyps fulvus (2n = 66) was generated from flow-sorted chromosomes. Together with previous generated probes from the stone curlew, Burhinus oedicnemus (2n = 42), a Charadriiformes species with atypical karyotype organization, we conducted multidirectional chromosome painting, including reciprocal chromosome painting between B. oedicnemus and G. fulvus and cross-species chromosome painting between B. oedicnemus and two accipitrid species (the Himalayan griffon, G. himalayensis 2n = 66, and the common buzzard, Buteo buteo, 2n = 68). In doing so, genome-wide homology maps between B. oedicnemus and three Accipitridae species were established. From there, a cladistic analysis using chromosomal characters and mapping of chromosomal changes on a consensus molecular phylogeny were conducted in order to search for cytogenetic signatures for different lineages within Accipitriformes.ConclusionOur study confirmed that the genomes of the diurnal birds of prey, especially the genomes of species in Accipitriformes excluding Cathartidae, have been extensively reshuffled when compared to other bird lineages. The chromosomal rearrangements involved include both fusions and fissions. Our chromosome painting data indicated that the Palearctic common buzzard (BBU) shared several common chromosomal rearrangements with some Old World vultures, and was found to be more closely related to other Accipitridae than to Neotropical buteonine raptors from the karyotypic perspective. Using both a chromosome-based cladistic analysis as well as by mapping of chromosomal differences onto a molecular-based phylogenetic tree, we revealed a number of potential cytogenetic signatures that support the clade of Pandionidae (PHA) + Accipitridae. In addition, our cladistic analysis using chromosomal characters appears to support the placement of osprey (PHA) in Accipitridae.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-015-0484-0) contains supplementary material, which is available to authorized users.

Genome-Scale Identification of Cell-Wall-Related Genes in Switchgrass through Comparative Genomics and Computational Analyses of Transcriptomic Data

Large numbers of plant cell-wall (CW)-related genes have been identified or predicted in several plant genomes such as Arabidopsis thaliana, Oryza sativa (rice), and Zea mays (maize), as results of intensive studies of these organisms in the past 2 decades. However, no such gene list has been identified in switchgrass (Panicum virgatum), a key bioenergy crop. Here, we present a computational study for prediction of CW genes in switchgrass using a two-step procedure: (i) homology mapping of all annotated CW genes in the fore-mentioned species to switchgrass, giving rise to a total of 991 genes, and (ii) candidate prediction of CW genes based on switchgrass genes co-expressed with the 991 genes under a large number of experimental conditions. Specifically, our co-expression analyses using the 991 genes as seeds led to the identification of 104 large clusters of co-expressed genes, each referred to as a co-expression module (CEM), covering 830 of the 991 genes plus 823 additional genes that are strongly co-expressed with some of the 104 CEMs. These 1653 genes represent our prediction of CW genes in switchgrass, 112 of which are homologous to predicted CW genes in Arabidopsis. Functional inference of these genes is conducted to derive the possible functional relations among these predicted CW genes. Overall, these data may offer a highly useful information source for cell-wall biologists of switchgrass as well as plants in general.

Upscaling the porosity of the Callovo-Oxfordian mudstone from the pore scale to the formation scale; insights from the 3H-PMMA autoradiography technique and SEM BSE imaging

The Callovo-Oxfordian mudstone (Meuse/Haute-Marne, France) is currently considered as the host rock barrier for a deep geological repository. The intimate relationships between the porosity and mineralogy of this host rock were investigated at the small scale (μm–mm) and large scale (m–hm). At the small scale, we have adapted the 3H-PMMA autoradiographic method to map the porosity of the Callovo-Oxfordian mudstone. The 3H-PMMA autoradiographic method was improved in terms of its spatial resolution. 3H-PMMA porosity maps were then compared to homologous mineral maps (clay minerals, carbonates and tectosilicates) built from scanning electron microscopy images (using back-scattered electron imaging). Based on an inversion procedure, the specific porosity of each mineral group was estimated from the mineral and porosity maps. We found that the spatial distribution of porosity at the small scale is mainly controlled by the spatial distribution of the clay matrix (the average porosity of the clay matrix is 40–45%), whereas quartz and carbonate mineral grains have low porosities (0–4%). At the geological formation scale, the porosity and mineralogy distributions were determined by logging tool techniques (nuclear magnetic resonance and spectral gamma-ray). The coupled evolution of clay content and porosity with depth was analyzed according to the porosity/mineralogy relationship defined at the small scale. Finally, we modeled the evolution of the porosity of the Callovo-Oxfordian mudstone with depth by considering the clay content and the effect of physical compaction during burial.