A numerical simulation of neural fields on curved geometries.

Despite the highly convoluted nature of the human brain, neural field models typically treat the cortex as a planar two-dimensional sheet of ne;urons. Here, we present an approach for solving neural field equations on surfaces more akin to the cortical geometries typically obtained from neuroimaging data. Our approach involves solving the integral form of the partial integro-differential equation directly using collocation techniques alongside efficient numerical procedures for determining geodesic distances between neural units. To illustrate our methods, we study localised activity patterns in a two-dimensional neural field equation posed on a periodic square domain, the curved surface of a torus, and the cortical surface of a rat brain, the latter of which is constructed using neuroimaging data. Our results are twofold: Firstly, we find that collocation techniques are able to replicate solutions obtained using more standard Fourier based methods on a flat, periodic domain, independent of the underlying mesh. This result is particularly significant given the highly irregular nature of the type of meshes derived from modern neuroimaging data. And secondly, by deploying efficient numerical schemes to compute geodesics, our approach is not only capable of modelling macroscopic pattern formation on realistic cortical geometries, but can also be extended to include cortical architectures of more physiological relevance. Importantly, such an approach provides a means by which to investigate the influence of cortical geometry upon the nucleation and propagation of spatially localised neural activity and beyond. It thus promises to provide model-based insights into disorders like epilepsy, or spreading depression, as well as healthy cognitive processes like working memory or attention.

Ornithine carbamoyltransferase deficiency: molecular characterization of 29 families.

Ornithine carbamoyltransferase deficiency is the most common inherited defect of the urea cycle. We examined 28 male and 9 female patients from 29 families and identified 25 distinct mutations in OTC, 14 of which were novel. Three novel missense mutations (p.Ala102Pro, p.Pro158Ser, p.Lys210Glu) and a novel deletion of the Leu43 are not directly involved either in the enzyme active site or in the intersubunit interactions; however, the mutations include conserved residues involved in intramolecular interaction network essential for the function of the enzyme. Three novel large deletions - a 444 kb deletion affecting RPGR, OTC and TSPAN7, a 10 kb-deletion encompassing OTC exons 5 and 6 and a 24.5 kb-deletion encompassing OTC exons 9 and 10 - have probably been initiated by double strand breaks at recombination-promoting motifs with subsequent non-homologous end joining repair. Finally, we present a manifesting heterozygote carrying a hypomorphic mutation p.Arg129His in combination with unfavorably skewed X-inactivation in three peripheral tissues.

SPRED1 mutations (Legius syndrome): another clinically useful genotype for dissecting the neurofibromatosis type 1 phenotype.

OBJECTIVE: Mutations of the SPRED1 gene, one of a family of Sprouty (Spry)/Spred proteins known to "downregulate" mitogen activated protein kinase (MAPK) signalling, have been identified in patients with a mild neurofibromatosis type 1 (NF1) phenotype with pigmentary changes but no neurofibromas (Legius syndrome).To ascertain the frequency of SPRED1 mutations as a cause of this phenotype and to investigate whether other SPRED/SPRY genes may be causal, a panel of unrelated mild NF1 patients were screened for mutations of the SPRED1-3 and the SPRY1-4 genes. METHODS: 85 patients with a mild NF1 phenotype were screened for SPRED1 mutations. 44 patients negative for both NF1 and SPRED1 mutations were then screened for SPRED2-3 and SPRY1-4 mutations. Complexity analysis was applied to analyse the flanking sequences surrounding the identified SPRED1 mutations for the presence of direct and inverted repeats or symmetric sequence elements in order to infer probable mutational mechanism. RESULTS: SPRED1 mutations were identified in 6 cases; 5 were novel and included 3 nonsense (R16X, E73X, R262X), 2 frameshift (c.1048_c1049 delGG, c.149_1152del 4 bp), and a single missense mutation (V44D). Short direct or inverted repeats detected immediately adjacent to some SPRED1 mutations may have led to the formation of the microdeletions and base pair substitutions. DISCUSSION: The identification of SPRED1 gene mutation in NF1-like patients has major implications for counselling NF1 families.

Germline and somatic NF1 gene mutation spectrum in NF1-associated malignant peripheral nerve sheath tumors (MPNSTs).

About 10% of neurofibromatosis type 1 (NF1) patients develop malignant peripheral nerve sheath tumors (MPNSTs) and represent considerable patient morbidity and mortality. Elucidation of the genetic mechanisms by which inherited and acquired NF1 disease gene variants lead to MPNST development is important. A study was undertaken to identify the constitutional and somatic NF1 mutations in 34 MPNSTs from 27 NF1 patients. The NF1 germline mutations identified in 22 lymphocytes DNA from these patients included seven novel mutations and a large 1.4-Mb deletion. The NF1 germline mutation spectrum was similar to that previously identified in adult NF1 patients without MPNST. Somatic NF1 mutations were identified in tumor DNA from 31 out of 34 MPNSTs, of which 28 were large genomic deletions. The high prevalence (>90%) of such deletions in MPNST contrast with the =or<20% found in benign neurofibromas and is indicative of the involvement of different mutational mechanisms in these tumors. Coinactivation of the TP53 gene by deletion, or by point mutation along with NF1 gene inactivation, is known to exacerbate disease symptoms in NF1, therefore TP53 gene inactivation was screened. DNA from 20 tumors showed evidence for loss of heterozygosity (LOH) across the TP53 region in 11 samples, with novel TP53 point mutations in four tumors.

An absence of cutaneous neurofibromas associated with a 3-bp inframe deletion in exon 17 of the NF1 gene (c.2970-2972 delAAT): evidence of a clinically significant NF1 genotype-phenotype correlation.

Neurofibromatosis type 1 (NF1) is characterized by cafe-au-lait spots, skinfold freckling, and cutaneous neurofibromas. No obvious relationships between small mutations (<20 bp) of the NF1 gene and a specific phenotype have previously been demonstrated, which suggests that interaction with either unlinked modifying genes and/or the normal NF1 allele may be involved in the development of the particular clinical features associated with NF1. We identified 21 unrelated probands with NF1 (14 familial and 7 sporadic cases) who were all found to have the same c.2970-2972 delAAT (p.990delM) mutation but no cutaneous neurofibromas or clinically obvious plexiform neurofibromas. Molecular analysis identified the same 3-bp inframe deletion (c.2970-2972 delAAT) in exon 17 of the NF1 gene in all affected subjects. The Delta AAT mutation is predicted to result in the loss of one of two adjacent methionines (codon 991 or 992) ( Delta Met991), in conjunction with silent ACA-->ACG change of codon 990. These two methionine residues are located in a highly conserved region of neurofibromin and are expected, therefore, to have a functional role in the protein. Our data represent results from the first study to correlate a specific small mutation of the NF1 gene to the expression of a particular clinical phenotype. The biological mechanism that relates this specific mutation to the suppression of cutaneous neurofibroma development is unknown.

Gross deletions and translocations in human genetic disease.

Translocations and gross deletions constitute an important cause of both cancer and inherited disease. Such gene rearrangements are non-randomly distributed in the human genome as a consequence of selection for growth advantage and/or the inherent potential of some DNA sequences to be frequently involved in breakage and recombination. Chromosomal rearrangements are generated by a variety of recombinational processes, each characterised by mechanism-specific DNA sequence features. Various types of recombinogenic motifs have been shown to promote non-homologous end joining whilst direct repeats may mediate homologous recombination. In addition, repetitive sequence elements can facilitate the formation of secondary structure between DNA ends at translocation or gross deletion breakpoints, and in so doing, may play a role in illegitimate recombination. Although results from DNA breakpoint studies are broadly consistent with a role for homologous unequal recombination in deletion mutagenesis and a role for non-homologous recombination in the generation of translocations, homologous recombination and non-homologous end joining are unlikely to be mutually exclusive mechanisms. Thus, chromosomal rearrangements will often represent the net result of multiple highly complex molecular interactions that are not always readily explicable.

Two novel severe mutations in the pancreatic secretory trypsin inhibitor gene (SPINK1) cause familial and/or hereditary pancreatitis.

Mutations in the serine protease inhibitor Kazal type 1 gene (SPINK1) encoding pancreatic secretory trypsin inhibitor (PSTI) have recently been found to be associated with chronic pancreatitis. Nevertheless, knowledge of severe mutations is particularly scarce, both in terms of number and in the extent of clinical information. The aim of this study was to expand the known spectrum of such mutations. 46 unrelated families, each including at least two pancreatitis patients and carrying neither cationic trypsinogen (PRSS1) mutations nor the frequent SPINK1 N34S mutation, participated in this study. The four exons and their flanking sequences of the SPINK1 gene were screened by denaturing high performance liquid chromatography analysis (DHPLC); and mutations were identified by direct sequencing. A heterozygous microdeletion mutation (c.27delC), which occurs within a symmetric element, was identified in two families. In one family, c.27delC showed segregation with the disease across two generations, with a penetrance of up to 75%. But in the other family, however, the same mutation manifested as a low-penetrance susceptibility factor. In addition, a novel heterozygous splicing mutation, c.87+1G>A (G>A substitution at nucleotide +1 of intron 2) was found in one family with familial pancreatitis. Our results also helped to resolve the sharply differing views about PSTI's role in pancreatitis.

[A rapid method for detecting interconnections between functionally and/or evolutionary close biological sequences]. / Bystryi metod vyiavleniia vzaimosviazei v podborkakh funktsional'no i/ili évoliutsionno blizkikh biologicheskikh tekstov.

It is now common practice to retrieve, by key words, highly specialized selections of sequences from general-purpose databases such as EMBL, GenBank, etc. The sequences included in a selection are often interconnected, which means that there are duplications, embeddings, intersections, homology, common structural elements. Knowledge of these interconnections is necessary for further processing of the sequences. We propose a rapid (single scan) method for identification of such interconnections by means of complexity analysis that generalizes the Lempel-Ziv approach. Analysis of a selection of 5'-flanking regions of vertebrate growth hormone genes from EMBL is presented as an example.

Changes in primary DNA sequence complexity influence the phenotypic consequences of mutations in human gene regulatory regions.

No general rules have been proposed to account for the functional consequences of gene regulatory mutations. In a first attempt to establish the nature of such rules, an analysis was performed of the DNA sequence context of 153 different single base-pair substitutions in the regulatory regions of 65 different human genes underlying inherited disease. Use of a recently proposed measure of DNA sequence complexity (taking into account the level of structural repetitiveness of a DNA sequence, rather than simply the oligonucleotide composition) has served to demonstrate that the concomitant change in local DNA sequence complexity surrounding a substituted nucleotide is related to the likelihood of a regulatory mutation coming to clinical attention. Mutations that led to an increase in complexity exhibited higher odds ratios in favour of pathological consequences than mutations that led to a decrease or left complexity unchanged. This relationship, however, was discernible only for pyrimidine-to-purine transversions. Odds ratios for other types of substitution were not found to be significantly associated with local changes in sequence complexity, even though a trend similar to that observed for Y-->R transversions was also apparent for transitions. These findings suggest that the maintenance of a defined level of DNA sequence complexity, or at least the avoidance of an increase in sequence complexity, is a critical prerequisite for the function of gene regulatory regions.

Promoter shuffling has occurred during the evolution of the vertebrate growth hormone gene.

Comparative studies of vertebrate gene promoter regions seldom detect gross rearrangements ('promoter shuffling') since such analyses usually employ relatively similar DNA sequences. Conversely, attempts to compare evolutionarily more divergent promoter sequences have been largely unsuccessful owing to the inability of conventional alignment procedures to deal with gross rearrangements. These limitations have been circumvented in the present study by using the novel technique of complexity analysis to identify modular components ('blocks') in the growth hormone (GH) gene promoter sequences of some 22 vertebrate species, from salmon to human. Significant rearrangement of blocks was found to have occurred, indicating that they have evolved as independent units. Some blocks appear to be ubiquitous, whereas others are restricted to a specific taxon. Considerable variation between orthologous GH gene promoters was apparent in terms of block length, copy number and relative location. It may be inferred that a wide variety of different mutational mechanisms have operated upon the GH gene promoter over evolutionary time. These include gross changes such as deletion, duplication, amplification, elongation, contraction, transposition, inversion and fusion, as well as the slow, steady accumulation of single base-pair substitutions. Thus the patchwork structure of the modular GH promoter region, and those of its paralogous GH2 and prolactin (PRL) counterparts, have continually been shuffled into new combinations through the rearrangement of pre-existing blocks. Although some of these changes may have had no influence on promoter function, others could have served to alter either the level of gene expression or the responsiveness of the promoter to external stimuli.

Identification of an intronic regulatory element in the human protein C (PROC) gene.

Regulatory DNA elements responsible for human protein C (PROC) gene expression have previously been identified in the upstream promoter region and first (untranslated) exon of the gene. Here we show that an additional sequence element located more than 500 bp downstream of the core promoter within intron 1 further enhances PROC promoter-driven reporter gene expression in human hepatoma cells. In common with core promoter constructs used in previous studies, the activity of this 3'-extended regulatory region is diminished by a naturally occurring promoter mutation. However, in contrast to constructs lacking intronic sequence, the promoter/intron regulatory region is repressed rather than activated by the transcription factor HNF-1. Using both conventional alignment procedures and complexity analysis to study the human and canine PROC sequences, we identified two conserved intronic regions, which were tested for their involvement in gene regulation. High-level gene expression from the intron-coupled promoter was dependent upon the integrity of a 142 bp sequence element, a duplicate copy of which is located in an upstream region of the PROC gene that possesses enhancer activity. These findings emphasise the potential importance of intragenic sequences for gene regulation and serve to illustrate that the results of PROC promoter/reporter gene experiments are critically dependent upon the sequence context. The identification of such intragenic elements is relevant to the analysis of human genetic disease since it will facilitate the detection and functional evaluation of regulatory mutations and polymorphisms.

Evolution of the proximal promoter region of the mammalian growth hormone gene.

The evolutionary relationship between the proximal growth hormone (GH) gene promoter sequences of 12 mammalian species was explored by comparison of their trinucleotide composition and by multiple sequence alignment. Both approaches yielded results that were consistent with the known fossil record-based phylogeny of the analysed sequences, suggesting that the two methods of tree reconstruction might be equally efficient and reliable. The pattern of evolution inferred for the mammalian GH gene promoters was found to vary both temporally and spatially. Thus, two distinct regions devoid of any evolutionary changes exist in primates, but only one of these 'gaps' is also observed in rodents, and neither is seen in ruminants. Furthermore, different evolutionary rates must have prevailed during different periods of evolutionary time and in different lineages, with a dramatic increase in evolutionary rate apparent in primates. Since a similar pattern of discontinuity has been previously noted for the evolution of the GH-coding regions, it may reflect the action of positive selection operating upon the GH gene as a single cohesive unit. Strong evidence for the action of gene conversion between primate GH gene promoters is provided by the fact that the human GH1 and GH2 sequences, which are thought to have diverged before the divergence of Old World monkeys from great apes, are more similar to one another than either is to the rhesus monkey GH2 promoter. Finally, it was noted that a number of nucleotide positions in the GH1 gene promoter that are polymorphic in humans appear to be highly conserved in mammals. This apparent conundrum, which could represent a caveat for the interpretation of phylogenetic footprinting studies, is potentially explicable in terms either of reduced genetic diversity in highly inbred animal species or insufficient population data from non-human species.

On the complexity measures of genetic sequences.

MOTIVATION: It is well known that the regulatory regions of genomes are highly repetitive. They are rich in direct, symmetric and complemented repeats, and there is no doubt about the functional significance of these repeats. Among known measures of complexity, the Ziv-Lempel complexity measure reflects most adequately repeats occurring in the text. But this measure does not take into account isomorphic repeats. By isomorphic repeats we mean fragments that are identical (or symmetric) modulo some permutation of the alphabet letters. RESULTS: In this paper, two complexity measures of symbolic sequences are proposed that generalize the Ziv-Lempel complexity measure by taking into account any isomorphic repeats in the text (rather than just direct repeats as in Ziv-Lempel). The first of them, the complexity vector, is designed for small alphabets such as the alphabet of nucleotides. The second is based on a search for the longest isomorphic fragment in the history of sequence synthesis and can be used for alphabets of arbitrary cardinality. These measures have been used for recognition of structural regularities in DNA sequences. Some interesting structures related to the regulatory region of the human growth hormone are reported.

Feature selection for genetic sequence classification.

MOTIVATION: Most of the existing methods for genetic sequence classification are based on a computer search for homologies in nucleotide or amino acid sequences. The standard sequence alignment programs scale very poorly as the number of sequences increases or the degree of sequence identity is <30%. Some new computationally inexpensive methods based on nucleotide or amino acid compositional analysis have been proposed, but prediction results are still unsatisfactory and depend on the features chosen to represent the sequences. RESULTS: In this paper, a feature selection method based on the Gamma (or near-neighbour) test is proposed. If there is a continuous or smooth map from feature space to the classification target values, the Gamma test gives an estimate for the mean-squared error of the classification, despite the fact that one has no a priori knowledge of the smooth mapping. We can search a large space of possible feature combinations for a combination which gives a smallest estimated mean-squared error using a genetic algorithm. The method was used for feature selection and classification of the large subunits of rRNA according to RDP (Ribosomal Database Project) phylogenetic classes. The sequences were represented by dinucleotide frequency distribution. The nearest-neighbour criterion has been used to estimate the predictive accuracy of the classification based on the selected features. For examples discussed, we found that the classification according to the first nearest neighbour is correct for 80% of the test samples. If we consider the set of the 10 nearest neighbours, then 94% of the test samples are classified correctly. AVAILABILITY: The principal novel component of this method is the Gamma test and this can be downloaded compiled for Unix Sun 4, Windows 95 and MS-DOS from http://www.cs.cf.ac.uk/ec/ CONTACT: s.margetts@cs.cf.ac.uk