Chromosomal rearrangement and genetic structure in the shrews of the Sorex araneus group

ABSTRACT The role of chromosomal rearrangements in the speciation process is much debated and many theoretical models have been developed. The shrews of the Sorex araneus group offer extraordinary opportunities to study the relationship between chromosomal variation and speciation. Indeed, this group of morphologically very similar species received a great deal of attention due to its karyotypic variability, which is mainly attributed to Robertsonian fusions. To explore the impact of karyotypic changes on genetic differentiation, we first studied the relationship between genetic and karyotypic structure among Alpine species and among chromosome races of the S. araneus group using Bayesian admixture analyses. The results of these analyses confirmed the taxonomic status of the studied species even though introgression can still be detected between species. Moreover, the strong spatial sub-structure highlighted the role of historical factors (e.g. geographical isolation) on genetic structure. Next, we studied gene flow at the chromosome level to address the question of the impact of chromosomal rearrangements on genetic differentiation. We used flow sorted chromosomes from three different karyotypic taxa of the S. araneus group to map microsatellite markers at the chromosóme arm level. We have been able to map 24 markers and to show that the karyotypic organisation of these taxa is well conserved, which suggests that these markers can be used for further inter-taxa studies. A general prediction of chromosomal speciation models is that genetic differentiation between two taxa should be larger across rearranged chromosomes than across chromosomes common to both taxa. We combined two approaches using mapped microsatellites to test this prediction. First, we studied the genetic differentiation among five shrew taxa placed at different evolutionary levels (i.e. within and among species). In this large scale study, we detected an overall significant difference in genetic structure between rearranged vs. common chromosomes. Moreover, this effect varied among pairwise comparisons, which allowed us to differentiate the role of the karyotypic complexity of hybrids and of the evolutionary divergence between taxa. Secondly, we compared the levels of gene flow measured across common vs. rearranged chromosomes in two karyotypically different hybrid zones (strong vs. low complexity of hybrids), which show similar levels of genetic structure. We detected a significantly stronger genetic structure across rearranged chromosomes in the hybrid zone showing the highest level of hybrid complexity. The large variance observed among loci suggested that other factors, such as the position of markers within the chromosome, also certainly affects genetic structure. In conclusion, our results strongly support the role of chromosomal rearrangements in the reproductive barrier and suggest their importance in speciation process of the S. araneus group. RESUME Le rôle des réarrangements chromosomiques dans les processus de spéciation est fortement débattu et de nombreux modèles théoriques ont été développés sur le sujet. Les musaraignes du groupe Sorex araneus présentent de nombreuses opportunités pour étudier les relations entre les variations chromosomiques et la spéciation. En effet, ce groupe d'espèces morphologiquement très proches a attiré l'attention des chercheurs en raison de sa variabilité caryotypique principalement attribuée à des fusions Robertsoniennes. Pour explorer l'impact des changements caryotypiques sur la différenciation génétique, nous avons tout d'abord étudié les relations entre la structure génétique et caryotypique de races chromosomiques et d'espèces alpine du groupe S. araneus en utilisant des analyses Bayesiennes d' « admixture ». Les résultats de ces analyses ont confirmé le statut taxonomique des espèces étudiées bien que nous ayons détecté de l'introgression entre espèces. L'observation d'une sous structure spatiale relativement forte souligne l'importance des facteurs historiques (telle que l'isolation géographique) sur la structure génétique de ce groupe. Ensuite, nous avons étudié le flux de gène au niveau des chromosomes pour aborder de manière directe la question de l'impact des réarrangements chromosomiques sur la différenciation génétique. En conséquence, nous avons utilisé des tris de chromosomes de trois taxons du groupe S. araneus pour localiser des marqueurs microsatellites au niveau du bras chromosomique. Au cours de cette étude, nous avons pu localiser 24 marqueurs et montrer une forte conservation dans l'organisation du caryotype de ces taxa. Ce résultat suggère que leur utilisation est appropriée pour des études entre taxa. Une prédiction générale à tous les modèles de spéciation chromosomique correspond à la plus grande différenciation génétique des chromosomes réarrangés que des chromosomes communs. Nous avons combiné deux approches utilisant des microsatellites localisés au niveau du bras chromosomique pour tester cette prédiction. Premièrement, nous avons étudié la différenciation génétique entre cinq taxa du groupe S. araneus se trouvant à des niveaux évolutifs différents (i.e. à l'intérieur et entre espèce). Au cours de cette étude, nous avons détecté une différenciation globale significativement plus élevée sur les chromosomes réarrangés. Cet effet varie entre les comparaisons, ce qui nous a permis de souligner le rôle de la complexité caryotypique des hybrides et du niveau de divergence évolutive entre taxa. Deuxièmement, nous avons comparé le flux de gènes des chromosomes communs et réarrangés dans deux zones d'hybridation caryotypiquement différentes (forte vs. Faible complexité des hybrides) mais présentant un niveau de différenciation génétique similaire. Ceci nous a permis de détecter une structure génétique significativement plus élevée sur les chromosomes réarrangés au centre de la zone d'hybridation présentant la plus grande complexité caryotypic. La forte variance observée entre loci souligne en outre le fait que d'autres facteurs, tel que la position du marqueur sur le chromosome, affectent probablement aussi la structure génétique mesurée. En conclusion, nos résultats supportent fortement le rôle des réarrangements chromosomiques dans la barrière reproductive entre espèces ainsi que leur importance dans les processus de spéciation des musaraignes du groupe S. araneus.

Uncertainties in the prediction of spatial variability of soil CO2 emissions and related properties

The soil CO2 emission has high spatial variability because it depends strongly on soil properties. The purpose of this study was to (i) characterize the spatial variability of soil respiration and related properties, (ii) evaluate the accuracy of results of the ordinary kriging method and sequential Gaussian simulation, and (iii) evaluate the uncertainty in predicting the spatial variability of soil CO2 emission and other properties using sequential Gaussian simulations. The study was conducted in a sugarcane area, using a regular sampling grid with 141 points, where soil CO2 emission, soil temperature, air-filled pore space, soil organic matter and soil bulk density were evaluated. All variables showed spatial dependence structure. The soil CO2 emission was positively correlated with organic matter (r = 0.25, p < 0.05) and air-filled pore space (r = 0.27, p < 0.01) and negatively with soil bulk density (r = -0.41, p < 0.01). However, when the estimated spatial values were considered, the air-filled pore space was the variable mainly responsible for the spatial characteristics of soil respiration, with a correlation of 0.26 (p < 0.01). For all variables, individual simulations represented the cumulative distribution functions and variograms better than ordinary kriging and E-type estimates. The greatest uncertainties in predicting soil CO2 emission were associated with areas with the highest estimated values, which produced estimates from 0.18 to 1.85 t CO2 ha-1, according to the different scenarios considered. The knowledge of the uncertainties generated by the different scenarios can be used in inventories of greenhouse gases, to provide conservative estimates of the potential emission of these gases.

Automatic Prediction of Facial Trait Judgments: Appearance vs. Structural Models

Evaluating other individuals with respect to personality characteristics plays a crucial role in human relations and it is the focus of attention for research in diverse fields such as psychology and interactive computer systems. In psychology, face perception has been recognized as a key component of this evaluation system. Multiple studies suggest that observers use face information to infer personality characteristics. Interactive computer systems are trying to take advantage of these findings and apply them to increase the natural aspect of interaction and to improve the performance of interactive computer systems. Here, we experimentally test whether the automatic prediction of facial trait judgments (e.g. dominance) can be made by using the full appearance information of the face and whether a reduced representation of its structure is sufficient. We evaluate two separate approaches: a holistic representation model using the facial appearance information and a structural model constructed from the relations among facial salient points. State of the art machine learning methods are applied to a) derive a facial trait judgment model from training data and b) predict a facial trait value for any face. Furthermore, we address the issue of whether there are specific structural relations among facial points that predict perception of facial traits. Experimental results over a set of labeled data (9 different trait evaluations) and classification rules (4 rules) suggest that a) prediction of perception of facial traits is learnable by both holistic and structural approaches; b) the most reliable prediction of facial trait judgments is obtained by certain type of holistic descriptions of the face appearance; and c) for some traits such as attractiveness and extroversion, there are relationships between specific structural features and social perceptions.

An evaluation of volume-based morphometry for prediction of mild cognitive impairment and Alzheimer's disease.

Voxel-based morphometry from conventional T1-weighted images has proved effective to quantify Alzheimer's disease (AD) related brain atrophy and to enable fairly accurate automated classification of AD patients, mild cognitive impaired patients (MCI) and elderly controls. Little is known, however, about the classification power of volume-based morphometry, where features of interest consist of a few brain structure volumes (e.g. hippocampi, lobes, ventricles) as opposed to hundreds of thousands of voxel-wise gray matter concentrations. In this work, we experimentally evaluate two distinct volume-based morphometry algorithms (FreeSurfer and an in-house algorithm called MorphoBox) for automatic disease classification on a standardized data set from the Alzheimer's Disease Neuroimaging Initiative. Results indicate that both algorithms achieve classification accuracy comparable to the conventional whole-brain voxel-based morphometry pipeline using SPM for AD vs elderly controls and MCI vs controls, and higher accuracy for classification of AD vs MCI and early vs late AD converters, thereby demonstrating the potential of volume-based morphometry to assist diagnosis of mild cognitive impairment and Alzheimer's disease.

The use of quantile trees in the prediction of the diameter distribution of a stand

Abstract

Genome-wide prediction of matrix attachment regions that increase gene expression in mammalian cells.

Gene transfer in eukaryotic cells and organisms suffers from epigenetic effects that result in low or unstable transgene expression and high clonal variability. Use of epigenetic regulators such as matrix attachment regions (MARs) is a promising approach to alleviate such unwanted effects. Dissection of a known MAR allowed the identification of sequence motifs that mediate elevated transgene expression. Bioinformatics analysis implied that these motifs adopt a curved DNA structure that positions nucleosomes and binds specific transcription factors. From these observations, we computed putative MARs from the human genome. Cloning of several predicted MARs indicated that they are much more potent than the previously known element, boosting the expression of recombinant proteins from cultured cells as well as mediating high and sustained expression in mice. Thus we computationally identified potent epigenetic regulators, opening new strategies toward high and stable transgene expression for research, therapeutic production or gene-based therapies.

How a haemosporidian parasite of bats gets around: the genetic structure of a parasite, vector and host compared.

Parasite population structure is often thought to be largely shaped by that of its host. In the case of a parasite with a complex life cycle, two host species, each with their own patterns of demography and migration, spread the parasite. However, the population structure of the parasite is predicted to resemble only that of the most vagile host species. In this study, we tested this prediction in the context of a vector-transmitted parasite. We sampled the haemosporidian parasite Polychromophilus melanipherus across its European range, together with its bat fly vector Nycteribia schmidlii and its host, the bent-winged bat Miniopterus schreibersii. Based on microsatellite analyses, the wingless vector, and not the bat host, was identified as the least structured population and should therefore be considered the most vagile host. Genetic distance matrices were compared for all three species based on a mitochondrial DNA fragment. Both host and vector populations followed an isolation-by-distance pattern across the Mediterranean, but not the parasite. Mantel tests found no correlation between the parasite and either the host or vector populations. We therefore found no support for our hypothesis; the parasite population structure matched neither vector nor host. Instead, we propose a model where the parasite's gene flow is represented by the added effects of host and vector dispersal patterns.

Nonlinear prediction based on score function

The linear prediction coding of speech is based in the assumption that the generation model is autoregresive. In this paper we propose a structure to cope with the nonlinear effects presents in the generation of the speech signal. This structure will consist of two stages, the first one will be a classical linear prediction filter, and the second one will model the residual signal by means of two nonlinearities between a linear filter. The coefficients of this filter are computed by means of a gradient search on the score function. This is done in order to deal with the fact that the probability distribution of the residual signal still is not gaussian. This fact is taken into account when the coefficients are computed by a ML estimate. The algorithm based on the minimization of a high-order statistics criterion, uses on-line estimation of the residue statistics and is based on blind deconvolution of Wiener systems [1]. Improvements in the experimental results with speech signals emphasize on the interest of this approach.

Structure-Based, Rational Design of T Cell Receptors.

Adoptive cell transfer using engineered T cells is emerging as a promising treatment for metastatic melanoma. Such an approach allows one to introduce T cell receptor (TCR) modifications that, while maintaining the specificity for the targeted antigen, can enhance the binding and kinetic parameters for the interaction with peptides (p) bound to major histocompatibility complexes (MHC). Using the well-characterized 2C TCR/SIYR/H-2K(b) structure as a model system, we demonstrated that a binding free energy decomposition based on the MM-GBSA approach provides a detailed and reliable description of the TCR/pMHC interactions at the structural and thermodynamic levels. Starting from this result, we developed a new structure-based approach, to rationally design new TCR sequences, and applied it to the BC1 TCR targeting the HLA-A2 restricted NY-ESO-1157-165 cancer-testis epitope. Fifty-four percent of the designed sequence replacements exhibited improved pMHC binding as compared to the native TCR, with up to 150-fold increase in affinity, while preserving specificity. Genetically engineered CD8(+) T cells expressing these modified TCRs showed an improved functional activity compared to those expressing BC1 TCR. We measured maximum levels of activities for TCRs within the upper limit of natural affinity, K D = ∼1 - 5 μM. Beyond the affinity threshold at K D < 1 μM we observed an attenuation in cellular function, in line with the "half-life" model of T cell activation. Our computer-aided protein-engineering approach requires the 3D-structure of the TCR-pMHC complex of interest, which can be obtained from X-ray crystallography. We have also developed a homology modeling-based approach, TCRep 3D, to obtain accurate structural models of any TCR-pMHC complexes when experimental data is not available. Since the accuracy of the models depends on the prediction of the TCR orientation over pMHC, we have complemented the approach with a simplified rigid method to predict this orientation and successfully assessed it using all non-redundant TCR-pMHC crystal structures available. These methods potentially extend the use of our TCR engineering method to entire TCR repertoires for which no X-ray structure is available. We have also performed a steered molecular dynamics study of the unbinding of the TCR-pMHC complex to get a better understanding of how TCRs interact with pMHCs. This entire rational TCR design pipeline is now being used to produce rationally optimized TCRs for adoptive cell therapies of stage IV melanoma.

The usefulness of phonological information for the acquisition of grammatical structure in child language development

During the process of language development, one of the most important tasks that children must face is that of identifying the grammatical category to which words in their language belong. This is essential in order to be able to form grammatically correct utterances. How do children proceed in order to classify words in their language and assign them to their corresponding grammatical category? The present study investigates the usefulness of phonological information for the categorization of nouns in English, given the fact that it is phonology the first source of information that might be available to prelinguistic infants who lack access to semantic information or complex morphosyntactic information. We analyse four different corpora containing linguistic samples of English speaking mothers addressing their children in order to explore the reliability with which words are represented in mothers’ speech based on several phonological criteria. The results of the analysis confirm the prediction that most of the words to which English learning infants are exposed during the first two years of life can be accounted for in terms of their phonological resemblance

Application of GA-PLS and GA-KPLS calculations for the prediction of the retention indices of essential oils

Genetic algorithm and partial least square (GA-PLS) and kernel PLS (GA-KPLS) techniques were used to investigate the correlation between retention indices (RI) and descriptors for 117 diverse compounds in essential oils from 5 Pimpinella species gathered from central Turkey which were obtained by gas chromatography and gas chromatography-mass spectrometry. The square correlation coefficient leave-group-out cross validation (LGO-CV) (Q²) between experimental and predicted RI for training set by GA-PLS and GA-KPLS was 0.940 and 0.963, respectively. This indicates that GA-KPLS can be used as an alternative modeling tool for quantitative structure-retention relationship (QSRR) studies.

Modeling of the Fluid-Structure Interaction in Inelastic Piping Systems

In many engineering applications, compliant piping systems conveying liquids are subjected to inelastic deformations due to severe pressure surges such as plastic tubes in modern water supply transmission lines and metallic pipings in nuclear power plants. In these cases the design of such systems may require an adequate modeling of the interactions between the fluid dynamics and the inelastic structural pipe motions. The reliability of the prediction of fluid-pipe behavior depends mainly on the adequacy of the constitutive equations employed in the analysis. In this paper it is proposed a systematic and general approach to consistently incorporate different kinds of inelastic behaviors of the pipe material in a fluid-structure interaction analysis. The main feature of the constitutive equations considered in this work is that a very simple numerical technique can be used for solving the coupled equations describing the dynamics of the fluid and pipe wall. Numerical examples concerning the analysis of polyethylene and stainless steel pipe networks are presented to illustrate the versatility of the proposed approach.

Prediction of exposed domains of envelope glycoprotein in Indian HIV-1 isolates and experimental confirmation of their immunogenicity in humans

We describe the impact of subtype differences on the seroreactivity of linear antigenic epitopes in envelope glycoprotein of HIV-1 isolates from different geographical locations. By computer analysis, we predicted potential antigenic sites of envelope glycoprotein (gp120 and gp4l) of this virus. For this purpose, after fetching sequences of proteins of interest from data banks, values of hydrophilicity, flexibility, accessibility, inverted hydrophobicity, and secondary structure were considered. We identified several potential antigenic epitopes in a B subtype strain of envelope glycoprotein of HIV-1 (IIIB). Solid- phase peptide synthesis methods of Merrifield and Fmoc chemistry were used for synthesizing peptides. These synthetic peptides corresponded mainly to the C2, V3 and CD4 binding sites of gp120 and some parts of the ectodomain of gp41. The reactivity of these peptides was tested by ELISA against different HIV-1-positive sera from different locations in India. For two of these predicted epitopes, the corresponding Indian consensus sequences (LAIERYLKQQLLGWG and DIIGDIRQAHCNISEDKWNET) (subtype C) were also synthesized and their reactivity was tested by ELISA. These peptides also distinguished HIV-1-positive sera of Indians with C subtype infections from sera from HIV-negative subjects.

Griscelli syndrome-type 2 in twin siblings: case report and update on RAB27A human mutations and gene structure

Griscelli syndrome (GS) is a rare autosomal recessive disorder caused by mutation in the MYO5A (GS1, Elejalde), RAB27A (GS2) or MLPH (GS3) genes. Typical features of all three subtypes of this disease include pigmentary dilution of the hair and skin and silvery-gray hair. Whereas the GS3 phenotype is restricted to the pigmentation dysfunction, GS1 patients also show primary neurological impairment and GS2 patients have severe immunological deficiencies that lead to recurrent infections and hemophagocytic syndrome. We report here the diagnosis of GS2 in 3-year-old twin siblings, with silvery-gray hair, immunodeficiency, hepatosplenomegaly and secondary severe neurological symptoms that culminated in multiple organ failure and death. Light microscopy examination of the hair showed large, irregular clumps of pigments characteristic of GS. A homozygous nonsense mutation, C-T transition (c.550C>T), in the coding region of the RAB27A gene, which leads to a premature stop codon and prediction of a truncated protein (R184X), was found. In patient mononuclear cells, RAB27A mRNA levels were the same as in cells from the parents, but no protein was detected. In addition to the case report, we also present an updated summary on the exon/intron organization of the human RAB27A gene, a literature review of GS2 cases, and a complete list of the human mutations currently reported in this gene. Finally, we propose a flow chart to guide the early diagnosis of the GS subtypes and Chédiak-Higashi syndrome.

The relationship between stem-form, stand-closure and site -conditions : the influence of environmental conditions on tree allometry and forest structure in west-central Alberta

Changes in the configuration of a tree stern result insignificant differences in its total volume and in the proportion of that volume that is merchantable timber. Tree allometry, as represented by stem-fo~, is the result of the vertical force of gravity and the horizontal force of wind. The effect of wind force is demonstrated in the relationship between stem-form, standclosure and site-conditions. An increase in wind force on the individual tree due to a decrease in stand density should produce a more tapered tree. The density of the stand is determined by the conditions that the trees are growing under. The ability of the tree to respond to increased wind force may also be a function of these conditions . This stem-form/stand-closure/site-conditions relationship was examined using a pre-existing database from westcentral Alberta. This database consisted of environmental, vegetation, soils and timber data covering a wide range of sites. There were 653 sample trees with 82 variables that formed the basis of the analysis. There were eight tree species consisting of Pinus contorta, Picea mariana, Picea engelmannii x glauca, Abies lasiocarpa, Larix laricina, Populus tremuloides, Betula papyrifera and Populus balsamifera plus a comprehensive all-species data set. As the actual conformation of the stern is very individual, stem-fo~was represented by the diameter at breast height to total height r~tio. The four stand-closure variables, crown closure, total basal area, total volume and total number of stems were reduced to total basal area and total number of stems utilizing a bivariate correlation matrix by species. Site-conditions were subdivided into macro, meso and micro variables and reduced in number 3 using cross-tabulations, bivariate correlation and principal components analysis as screening tools. The stem-fo~/stand-closure relationship was examined using bivariate correlation coefficients for stem-fo~ with total number of stems and stem-fo~ with total basal area. The stem-fo~/site-conditions and the stand-closure/site- conditions relationships were examined using multiple correlation coefficients. The stem-form/stand-closure/site-conditions relationship was examined using multiple correlation coefficients in separate analyses for both total number of stems and total basal area. An increase in stand-closure produced a decrease in stem-form for both total number of stems and total basal area for most species. There was a significant relationship between stem-form and site-conditions and between stand-closure and site-conditions for both total number of stems and total basal area for most species. There was a significant relationship between the stemform and site-conditions, including the stand-closure, for most species; total number of stems was involved independently of the site-conditions in the prediction of stem-form and total basal area was not. Larix laricina and Betula papyrifera were the exceptions to the trends observed with most species. The influence of both stand-closure (total number of stems in particular) and site-conditions (elevation in particular) suggest that forest management practices should include these- ecological parameters in determining appropriate restocking levels.