Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders

Purpose

To identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs).

Methods

We coupled phenotyping with exome or genome sequencing of 467 pedigrees with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations.

Results

Analyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43/467 probands (9.2%), and prioritized variants of uncertain significance in 70/467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (e.g., MYH10, KIF21B, TGFBR2, TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (e.g., CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (e.g., TUBA4A, KIF5C, CTNNA1, KLB, FGF21), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses.

Conclusion

This study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generated single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. Seventy-five percent of elements (44 of 59) validated in an in vivo transgenic reporter assay, demonstrating that single cell accessibility is a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieved significant reduction in our variant search space and nominated candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as new candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work provides novel non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

GATK-gCNV enables the discovery of rare copy number variants from exome sequencing data

Copy number variants (CNVs) are major contributors to genetic diversity and disease. While standardized methods, such as the genome analysis toolkit (GATK), exist for detecting short variants, technical challenges have confounded uniform large-scale CNV analyses from whole-exome sequencing (WES) data. Given the profound impact of rare and de novo coding CNVs on genome organization and human disease, we developed GATK-gCNV, a flexible algorithm to discover rare CNVs from sequencing read-depth information, complete with open-source distribution via GATK. We benchmarked GATK-gCNV in 7,962 exomes from individuals in quartet families with matched genome sequencing and microarray data, finding up to 95% recall of rare coding CNVs at a resolution of more than two exons. We used GATK-gCNV to generate a reference catalog of rare coding CNVs in WES data from 197,306 individuals in the UK Biobank, and observed strong correlations between per-gene CNV rates and measures of mutational constraint, as well as rare CNV associations with multiple traits. In summary, GATK-gCNV is a tunable approach for sensitive and specific CNV discovery in WES data, with broad applications.

Convergent coexpression of autism-associated genes suggests some novel risk genes may not be detectable in large-scale genetic studies

Autism spectrum disorder (ASD) is a heritable neurodevelopmental disorder characterized by deficits in social interactions and communication. Protein-altering variants in many genes have been shown to contribute to ASD; however, understanding the convergence across many genes remains a challenge. We demonstrate that coexpression patterns from 993 human postmortem brains are significantly correlated with the transcriptional consequences of CRISPR perturbations in human neurons. Across 71 ASD risk genes, there was significant tissue-specific convergence implicating synaptic pathways. Tissue-specific convergence was further demonstrated across schizophrenia and atrial fibrillation risk genes. The degree of ASD convergence was significantly correlated with ASD association from rare variation and differential expression in ASD brains. Positively convergent genes showed intolerance to functional mutations and had shorter coding lengths than known risk genes even after removing association with ASD. These results indicate that convergent coexpression can identify potentially novel genes that are unlikely to be discovered by sequencing studies.

Statistical and functional convergence of common and rare genetic influences on autism at chromosome 16p

The canonical paradigm for converting genetic association to mechanism involves iteratively mapping individual associations to the proximal genes through which they act. In contrast, in the present study we demonstrate the feasibility of extracting biological insights from a very large region of the genome and leverage this strategy to study the genetic influences on autism. Using a new statistical approach, we identified the 33-Mb p-arm of chromosome 16 (16p) as harboring the greatest excess of autism's common polygenic influences. The region also includes the mechanistically cryptic and autism-associated 16p11.2 copy number variant. Analysis of RNA-sequencing data revealed that both the common polygenic influences within 16p and the 16p11.2 deletion were associated with decreased average gene expression across 16p. The transcriptional effects of the rare deletion and diffuse common variation were correlated at the level of individual genes and analysis of Hi-C data revealed patterns of chromatin contact that may explain this transcriptional convergence. These results reflect a new approach for extracting biological insight from genetic association data and suggest convergence of common and rare genetic influences on autism at 16p.

Developmental Variability in Autism Across 17 000 Autistic Individuals and 4000 Siblings Without an Autism Diagnosis: Comparisons by Cohort, Intellectual Disability, Genetic Etiology, and Age at Diagnosis

Importance

Presence of developmental delays in autism is well established, yet few studies have characterized variability in developmental milestone attainment in this population.

Objective

To characterize variability in the age at which autistic individuals attain key developmental milestones based on co-occurring intellectual disability (ID), presence of a rare disruptive genetic variant associated with neurodevelopmental disorders (NDD), age at autism diagnosis, and research cohort membership.

Design

The study team harmonized data from 4 cross-sectional autism cohorts: the Autism Genetics Research Exchange (n = 3284; 1997-2015), The Autism Simplex Collection (n = 694; 2008-2011), the Simons Simplex Collection (n = 2753; 2008-2011), and the Simons Foundation Powering Autism Research for Knowledge (n = 10 367; 2016-present). The last sample further included 4145 siblings without an autism diagnosis or ID.

Participants

Convenience sample of 21 243 autistic individuals or their siblings without an autism diagnosis aged 4 to 17 years.

Main Outcomes and Measures

Parents reported ages at which participants attained key milestones including smiling, sitting upright, crawling, walking, spoon-feeding self, speaking words, speaking phrases, and acquiring bladder and bowel control. A total of 5295 autistic individuals, and their biological parents, were genetically characterized to identify de novo variants in NDD-associated genes. The study team conducted time-to-event analyses to estimate and compare percentiles in time with milestone attainment across autistic individuals, subgroups of autistic individuals, and the sibling sample.

Results

Seventeen thousand ninety-eight autistic individuals (mean age, 9.15 years; 80.8% male) compared with 4145 siblings without autism or ID (mean age, 10.2 years; 50.2% female) showed delays in milestone attainment, with median (IQR) delays ranging from 0.7 (0.3-1.6) to 19.7 (11.4-32.2) months. More severe and more variable delays in autism were associated with the presence of co-occurring ID, carrying an NDD-associated rare genetic variant, and being diagnosed with autism by age 5 years. More severe and more variable delays were also associated with membership in earlier study cohorts, consistent with autism's diagnostic and ascertainment expansion over the last 30 years.

Conclusions and Relevance

As the largest summary to date of developmental milestone attainment in autism, to our knowledge, this study demonstrates substantial developmental variability across different conditions and provides important context for understanding the phenotypic and etiological heterogeneity of autism.

Rare coding variation provides insight into the genetic architecture and phenotypic context of autism

Some individuals with autism spectrum disorder (ASD) carry functional mutations rarely observed in the general population. We explored the genes disrupted by these variants from joint analysis of protein-truncating variants (PTVs), missense variants and copy number variants (CNVs) in a cohort of 63,237 individuals. We discovered 72 genes associated with ASD at false discovery rate (FDR) ≤ 0.001 (185 at FDR ≤ 0.05). De novo PTVs, damaging missense variants and CNVs represented 57.5%, 21.1% and 8.44% of association evidence, while CNVs conferred greatest relative risk. Meta-analysis with cohorts ascertained for developmental delay (DD) (n = 91,605) yielded 373 genes associated with ASD/DD at FDR ≤ 0.001 (664 at FDR ≤ 0.05), some of which differed in relative frequency of mutation between ASD and DD cohorts. The DD-associated genes were enriched in transcriptomes of progenitor and immature neuronal cells, whereas genes showing stronger evidence in ASD were more enriched in maturing neurons and overlapped with schizophrenia-associated genes, emphasizing that these neuropsychiatric disorders may share common pathways to risk.

The female protective effect against autism spectrum disorder

Autism spectrum disorder (ASD) is diagnosed three to four times more frequently in males than in females. Genetic studies of rare variants support a female protective effect (FPE) against ASD. However, sex differences in common inherited genetic risk for ASD are less studied, particularly within families. Leveraging the Danish iPSYCH resource, we found siblings of female ASD cases (n = 1,707) had higher rates of ASD than siblings of male ASD cases (n = 6,270; p < 1.0 × 10-10). In the Simons Simplex and SPARK collections, mothers of ASD cases (n = 7,436) carried more polygenic risk for ASD than fathers of ASD cases (n = 5,926; 0.08 polygenic risk score [PRS] SD; p = 7.0 × 10-7). Further, male unaffected siblings under-inherited polygenic risk (n = 1,519; p = 0.03). Using both epidemiologic and genetic approaches, our findings strongly support an FPE against ASD's common inherited influences.

Whole exome sequencing analyses reveal gene-microbiota interactions in the context of IBD

OBJECTIVE

Both the gut microbiome and host genetics are known to play significant roles in the pathogenesis of IBD. However, the interaction between these two factors and its implications in the aetiology of IBD remain underexplored. Here, we report on the influence of host genetics on the gut microbiome in IBD.

DESIGN

To evaluate the impact of host genetics on the gut microbiota of patients with IBD, we combined whole exome sequencing of the host genome and whole genome shotgun sequencing of 1464 faecal samples from 525 patients with IBD and 939 population-based controls. We followed a four-step analysis: (1) exome-wide microbial quantitative trait loci (mbQTL) analyses, (2) a targeted approach focusing on IBD-associated genomic regions and protein truncating variants (PTVs, minor allele frequency (MAF) >5%), (3) gene-based burden tests on PTVs with MAF <5% and exome copy number variations (CNVs) with site frequency <1%, (4) joint analysis of both cohorts to identify the interactions between disease and host genetics.

RESULTS

We identified 12 mbQTLs, including variants in the IBD-associated genes IL17REL, MYRF, SEC16A and WDR78. For example, the decrease of the pathway acetyl-coenzyme A biosynthesis, which is involved in short chain fatty acids production, was associated with variants in the gene MYRF (false discovery rate <0.05). Changes in functional pathways involved in the metabolic potential were also observed in participants carrying rare PTVs or CNVs in CYP2D6, GPR151 and CD160 genes. These genes are known for their function in the immune system. Moreover, interaction analyses confirmed previously known IBD disease-specific mbQTLs in TNFSF15.

CONCLUSION

This study highlights that both common and rare genetic variants affecting the immune system are key factors in shaping the gut microbiota in the context of IBD and pinpoints towards potential mechanisms for disease treatment.

Detection of de novo copy number deletions from targeted sequencing of trios

MOTIVATION

De novo copy number deletions have been implicated in many diseases, but there is no formal method to date that identifies de novo deletions in parent-offspring trios from capture-based sequencing platforms.

RESULTS

We developed Minimum Distance for Targeted Sequencing (MDTS) to fill this void. MDTS has similar sensitivity (recall), but a much lower false positive rate compared to less specific CNV callers, resulting in a much higher positive predictive value (precision). MDTS also exhibited much better scalability.

AVAILABILITY AND IMPLEMENTATION

MDTS is freely available as open source software from the Bioconductor repository.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Activity-Dependent Degradation of the Nascentome by the Neuronal Membrane Proteasome

Activity-dependent changes in neuronal function require coordinated regulation of the protein synthesis and protein degradation machinery to maintain protein homeostasis, critical for proper neuronal function. However, the biochemical evidence for this balance and coordination is largely lacking. Leveraging our recent discovery of a neuronal-specific 20S membrane proteasome complex (NMP), we began exploring how neuronal activity regulates its function. Here, we found that the NMP degrades exclusively a large fraction of ribosome-associated nascent polypeptides that are being newly synthesized during neuronal stimulation. Using deep-coverage and global mass spectrometry, we identified the nascent protein substrates of the NMP, which included products encoding immediate-early genes, such as c-Fos and Npas4. Intriguingly, we found that turnover of nascent polypeptides and not full-length proteins through the NMP occurred independent of canonical ubiquitylation pathways. We propose that these findings generally define a neuronal activity-induced protein homeostasis program of coordinated protein synthesis and degradation through the NMP.

Leaf functional trait variation associated with salt tolerance in perennial ryegrass

Salinity is one of major environmental stresses that dramatically threaten plant growth, and variations in genetic structure and functional traits have important effects on the salt tolerance of perennial ryegrass (Lolium perenne L.). The objectives of this study were to: (i) assess the inter-clonal variation of functional traits of accessions among geographic groups or between wild and commercial groups in response to salt stress; (ii) develop a mathematical model to effectively assess salt tolerance of perennial ryegrass accessions originating from different geographic populations; and (iii) determine the relation between spatial genetic structure and salt tolerance in perennial ryegrass. Wide variations were found among the accessions for seven functional traits. One regression model (F = 0.49 × F1 + 0.303 × F2 + 0.207 × F3) was established to ascertain salt tolerance of each accession. The highest variation of the traits and salt tolerance were obtained for accessions from the European group. Wild accessions exhibited more variation in functional traits and salt tolerance than commercial cultivars. Both molecular marker techniques and functional traits were used to conduct phylogenetic analysis, and the majority of accessions from the same or adjacent regions were clustered into the same group or subgroup. The perennial ryegrass accessions with similar salt tolerance had a close phylogenetic background. The patterns in functional trait variations associated with salt tolerance might allow acceleration of the process for improving salt stress resistance in perennial ryegrass.

Distribution of PBDEs in air particles from an electronic waste recycling site compared with Guangzhou and Hong Kong, South China

Air samples of total suspended particles (TSP, particles less than 30-60 microm), and particles with aerodynamic diameter smaller than 2.5 microm (PM(2.5)) were collected simultaneously at Guiyu (an electronic waste recycling site), three urban sites in Hong Kong and two urban sites in Guangzhou, South China from 16 August to 17 September 2004. Twenty-two PBDE congeners (BDE-3, -7, -15, -17, -28, -49, -71, -47, -66, -77, -100, -119, -99, -85, -126, -154, -153, -138, -156, -184, -183, -191) in TSP and PM(2.5) were measured. The results showed that the overall average concentrations of TSP and PM(2.5) collected at Guiyu were 124 and 62.1 microg m(-3), respectively. The monthly concentrations of the sum of 22 BDE congeners contained in TSP and PM(2.5) at Guiyu were 21.5 and 16.6 ng m(-3), with 74.5 and 84.3%, contributed by nine congeners (BDE-28, -47, -66, -100, -99, -154, -153, -183 and -191 respectively). This pattern was similar to Tsuen Wan site of Hong Kong. Two urban sites of Guangzhou had the same congener pattern, but were different from Yuen Long and Hok Tsui sites of Hong Kong. The results also showed that the amount of mono to penta brominated congeners, which are more toxic, accounted for 79.4-95.6% of Sigma(22)PBDEs from all sites. All congeners tested in Guiyu were up to 58-691 times higher than the other urban sites and more than 100 times higher than other studies reported elsewhere. The higher concentration in the air was due to heating or opening burning of electronic waste since PBDEs are formed when plastics containing brominated flame retardants are heated.

Seasonal variation of the particle size distribution of n-alkanes and polycyclic aromatic hydrocarbons (PAHs) in urban aerosol of Guangzhou, China

Seasonal aerosol samples have been collected by Andersen Hi-Vol pumping system equipped with a five stage cascade impactor and a backup filter (size range: 10-7.2 microm, 7.2-3.0 microm, 3.0-1.5 microm, 1.5-0.95 microm, 0.95-0.49 microm, <or=0.49 microm) in the Liwan district, Guangzhou. n-Alkanes were measured using gas chromatography and PAHs were measured using gas chromatography/mass spectrometry analysis. The bimodal log-normal distributions of n-alkanes and semi-volatile PAHs were found, while for non-volatile PAHs that was unimodal, so much as the mode of semi-volatile PAHs was similar with that of the particles. The n-alkanes and PAHs were preferably associated with fine particles. C (max) (carbon number maximum) (C(22)-C(26)), CPI (carbon preference index) (1.12-1.21), U/R (unresolved to resolved components ratio) (7.42-10.7), wax% (0.9-3.12%) and the diagnostic ratios for PAHs revealed that vehicular emission was the major source of these organic compounds during the study periods, while the contribution of epicuticular waxes emitted by terrestrial plants was minor. CPI(2) (values for petrogenic hydrocarbons), CPI(3) (values for biogenic n-alkanes) and wax% revealed that the natural preferentially accumulated in the larger aerosol while the anthropogenic in the smaller. In addition, the different MMDs (mass median diameters) for n-alkanes and PAHs were observed in different seasons. The MMDs for n-alkanes and PAHs were higher in autumn/winter than those in spring/summer. The seasonal effect was related to the hydrocarbon content in the individual particulate fractions, showing a preferential association of n-alkanes and PAHs with larger particles in the autumn/winter season.

Development of a Compound-Specific Carbon Isotope Analysis Method for Atmospheric Formaldehyde via NaHSO3 and Cysteamine Derivatization

A novel method has been developed for the compound-specific carbon isotope analysis of atmospheric formaldehyde using gas chromatography/combustion/isotope ratio mass spectrometry (GC/C/IRMS). The method allows the determination of the delta13C value for atmospheric formaldehyde at nanogram levels with higher precision and lower detection limit. In the present work, atmospheric formaldehyde was collected using NaHSO3-coated Sep-Pak silica gel cartridges, washed out by water, then derivatized by cysteamine of known delta13C value, and the delta13C value of its derivative (thiazolidine) determined by GC/C/IRMS. Finally, the delta13C value of atmospheric formaldehyde could be calculated by a simple mass balance equation between formaldehyde, cysteamine, and thiazolidine. Using three formaldehydes with different delta13C values, calibration experiments were carried out over large ranges of formaldehyde concentrations. The carbon isotope analysis method achieved excellent reproducibility and high accuracy. There was no carbon isotopic fractionation throughout the derivatization processes. The differences in the carbon isotopic compositions of thiazolidine between the measured and predicted values were always <0.5 per thousand, within the specifications of the GC/C/IRMS system. The present method was also compared with the previous 2,4-dinitrophenylhydrazine derivatization method, and this method could be performed with lower analytical error and detection limit. Using this method, four 6-h ambient atmospheric formaldehyde samples were consecutively collected from 8 to 9 March 2005. The results showed that the delta13C values of atmospheric formaldehyde were different during the daytime and nighttime. This method proved suitable for the routine operation and may provide additional insight on sources and sinks of atmospheric formaldehyde.

[The determination of porphyrin carbon isotope composition by gas chromatography-isotope ratio monitoring mass spectrometry technique]

The porphyrin carbon isotope composition can be used to explore the precursor of porphyrin, oil-oil and oil-source rock correction and calculation of paleo P CO2. The conventional method is limited because of its time consuming and large sample size (several mg of individual porphyrin) required. Therefore, it hampers the application of porphyrin carbon isotope composition into the chemistry and geoscience. The present paper describes a quantification method to prepare bis-(tert-butyldimethylsiloxy) silicon (IV) [(TBDMSO)2Si(IV)] porphyrin which is sufficiently volatile at 300 degrees C and can be used for GC-IRMS analysis. The analysis of carbon isotope composition of aetio I as the form of free base, nickel, demetalization derivative, silicon(IV) and (TBDMSO)2Si(IV) have shown that aetio I porphyrin has no obvious isotope fractionation in the whole synthesis procedure for (TBDMSO)2Si(IV) porphyrin. The carbon isotope study on the porphyrin mixtures of aetio I and OEP indicates that isotope exchange between porphyrins during the synthesis of (TBDMSO)2Si(IV) porphyrin is absent. The method can be applied to the determination of porphyrin carbon isotope compositions. The advantages of the method are time saving, less sample size and lower standard deviation.

[A study on the contents of fat and fatty acids in five shellfishes]

The contents of fat and fatty acids in Callista erycina (Linnaeus), Paphia (Paratapes) undulata (Born), Meretrix meretrix (Linnaeus), Chlamys farreri (Jones et Preston) and Patinopecten yessoensis (Jay) were studied. Fat was extracted with Bligh & Dyer (B&D) method. The lipid classes were transesterified with potassium hydroxide in methanol. Fatty acid methyl esters (FAMEs) were assayed with GC-MS and polar capillary column (HP-INNOWax 30 m x 0.25 mm i.d. x 0.25 micron). GC injector temperature was 220 degrees C. The column temperature was programmed from 150 degrees C (1 min) to 200 degrees C at 10 degrees C/min and then from 200 degrees C to 250 degrees C at 2 degrees C/min. FAMEs were identified by MS library, and part by their standards. Total identified fatty acids were over 99% for all samples. Fat contents of them were all over 1% by wet samples. And ratios between omega-3PUFA and omega-6PUFA were above 2 by and large. Patinopecten yessoensis (Jay) contains more fat and the valuable fatty acids, EPA and DHA. It is suitable to use it as the source of EPA and DHA.

Fractionation and composition of colloidal and suspended particulate materials in rivers

The association of pollutants (nutrients, heavy metals and organic compounds) with colloidal and suspended particle matter (SPM) plays a dominant role in determining their transport, fate, biogeochemistry, bioavailability and toxicity in natural waters. A scheme for the fractionation and composition of colloidal and SPM from river waters has been tested. All four separation methods, i.e. sieving, continuous flow centrifugation, tangential flow filtration, sedimentation field-flow fractionation, were for the first time used to separate five size particulate fractions from river. Significant (gram) amounts of colloidal material (<1 microm) in three size ranges, nominally 1-0.2, 0.2-0.006 and 0.006-0.003 microm were obtained. The separation scheme was able to process large samples (100 l), within reasonable times (1 day) and the apparatus was portable. The aquatic colloid size was also characterized with high resolution by using sedimentation field-flow fractionation technique. The mass-based particle size distribution for the water sample showed a broad size distribution between 0.05 and 0.4 microm with the maximum around 0.14 microm. There was a systematic increase in the content of organic carbon (estimated by loss on ignition), Mg, Ca, Na, Cu and Zn with decreasing particle size, highlighting the importance of the colloidal (<1 microm) fraction. It was concluded that the colloidal Cu and Zn concentrations in rivers might be much higher than those reported before.

[Stability of fatty acids in stored oyster]

The contents of twenty-three kinds of fatty acids in the lyophilized oyster that stored for 0, 15, 30, 45, 60, 75, 90 days were studied. Oyster fat was extracted from the powder by means of SFE. The extraction was performed for 40 min at a pressure of 37 MPa and a temperature of 50 degrees C with supercritical carbon dioxide containing 8% (V/V) ethanol at a flow-rate of 2 mL/min as liquid carbon dioxide, and the recovery of oyster fat by extraction was over 99%. After being extracted, fat was esterified with potassium hydroxide and methanol. Methyl esters of fatty acid were separated and determined by 30 m x 0.25 mm i.d. x 0.25 micron HP-INNOWax capillary column and MS detector. The injector temperature was 220 degrees C. The column temperature was programmed from 150 degrees C (1 min) to 200 degrees C at 10 degrees C/min and then from 200 degrees C to 250 degrees C at 2 degrees C/min. Twenty-three peaks were identified with gas chromatography/mass spectrometry, and quantified with area normalization method. The variations of contents of them were shown. During storage, the contents of saturated fatty acids and mono-unsaturated fatty acids were getting higher, and those of polyunsaturated fatty acids were getting lower. The decrease of them was gradual, and there was no special period of stability. And the stability of fatty acids in oyster related to the degree of unsaturation of them. The higher the unsaturation the lower the stable it was. After being stored for 90 days, the content of EPA decreased from 16.94% to 5.43% and that of DHA from 9.25% to 2.86%.

Selective A1-adenosine receptor antagonists identified using yeast Saccharomyces cerevisiae functional assays

Evaluation of a biased "library" of pyrrolo[2,3-d]pyrimidines using yeast-based functional assays expressing human A1- and A2a-adenosine receptors, led to the A1 selective antagonist 4b. A direct correlation between yeast functional activity and binding data was established. Practical compounds with polar residues at C-4 of the pyrrolopyrimidine system required H-bond donor functionality for high potency.

Design and synthesis of a pyridone-based phosphotyrosine mimetic

A novel pyridone-based tyrosine analog, 6, has been designed to mimic the binding interaction of SH2 domains with phosphotyrosine (pTyr) containing peptides. Synthesis of 6 features a key Pd catalyzed coupling of beta-iodoalanine with phosphonomethyl 4-pyridone triflate.

Structural features associated with reactive metabolite formation in clozapine analogues

Clozapine is associated with a high incidence of agranulocytosis. We had previously found that it is oxidized by granulocytes, or simply HOCl, to a reactive metabolite that irreversibly binds to the cells, and we proposed that this reactive metabolite is responsible for clozapine-induced agranulocytosis. The reactive metabolite appeared to be a nitrenium ion formed by chlorination of the nitrogen bridge between the two aromatic rings. If this is correct, analogs that contain this structural feature should also be oxidized to a reactive intermediate while those not possessing this feature would, at least, not form the same type of reactive intermediate and, therefore, may not induce agranulocytosis. We tested the first part of this hypothesis with three clozapine analogs that do contain a nitrogen bridge and three that do not. Consistent with the hypothesis, the three analogs that do contain the nitrogen bridge formed reactive intermediates that could be trapped with glutathione when oxidized by HOCl, myeloperoxidase or activated neutrophils. In contrast, we found no evidence of a reactive intermediate on oxidation of analogs that contained an oxygen or sulfur bridge rather than a nitrogen bridge. If such reactive metabolites are responsible for drug-induced agranulocytosis, it should be possible to use such a simple screening method to test drugs at an early stage in their development for the potential to induce agranulocytosis.

[Clinical study on action of artificial respiration produced by pressing in massage under anesthesia for lumbar intervertebral disc prolapse]

OBJECTIVE

To observe whether the gas exchange is sufficient for the requirement of body when the patient was prone in position and pressed in the course of massage under intravenous anesthesia with sodium pentothal (SP) for the treatment of lumbar intervertebral disc prolapse.

METHODS

The oxygen would not be given to the patients whose heart and lung functions were good during the course of anesthesia and massage, respiratory function and blood gas analysis were studied on 20 patients selected randomly.

RESULTS

Tidal volume was smaller and respiratory rate was faster than normal respiration. Saturation pulse oxygen was still normal in all patients. There were no significant changes in blood gas analysis before spinal injection and after massage (P > 0.05).

CONCLUSION

The pressing could attain the action of artificial respiration. The respiration pattern was similar to high frequency positive pressure ventilation. The action of vibration and diffusion had normal respiratory effect.

Sequence-dependent variations in the 31P NMR spectra and backbone torsional angles of wild-type and mutant Lac operator fragments

Assignment of the 31P resonances of a series of six sequenced-related tetradecamer DNA duplexes, d(TGTGAGCGCTCACA)2, d(TATGAGCGCTCATA)2, d(TCTGAGCGCTCAGA)2, d(TGTGTGCGCACACA)2, d(TGTGACGCGTCACA)2 and d(CACAGTATACTGTG)2, related to the lac operator DNA sequence was determined either by site-specific 17O labeling of the phosphoryl groups or by two-dimensional 1H-31P pure absorption phase constant time (PAC) heteronuclear correlation spectroscopy. J(H3'-P) coupling constants for each of the phosphates of the tetradecamers were obtained from 1H-31P J-resolved selective proton flip 2D spectra. By use of a modified Karplus relationship the C4'-C3'-O3'-P torsional angles (epsilon) were obtained. Comparison of the 31P chemical shifts and J(H3'-P) coupling constants of these sequences has allowed greater insight into those various factors responsible for 31P chemical shift variations in oligonucleotides and provided an important probe of the sequence-dependent structural variation of the deoxyribose phosphate backbone of DNA in solution. These sequence-specific variations in the conformation of the DNA sugar phosphate backbone of various lac operator DNA sequences can possibly explain the sequence-specific recognition of DNA by DNA binding proteins, as mediated through direct contacts between the phosphates and the protein.

Assignments of 31P NMR resonances in oligodeoxyribonucleotides: origin of sequence-specific variations in the deoxyribose phosphate backbone conformation and the 31P chemical shifts of double-helical nucleic acids

It is now possible to unambiguously assign all 31P resonances in the 31P NMR spectra of oligonucleotides by either two-dimensional NMR techniques or site-specific 17O labeling of the phosphoryl groups. Assignment of 31P signals in tetradecamer duplexes, (dTGTGAGCGCTCACA)2, (dTAT-GAGCGCTCATA)2, (dTCTGAGCGCTCAGA)2, and (dTGTGTGCGCACACA)2, and the dodecamer duplex d(CGTGAATTCGCG)2 containing one base-pair mismatch, combined with additional assignments in the literature, has allowed an analysis of the origin of the sequence-specific variation in 31P chemical shifts of DNA. The 31P chemical shifts of duplex B-DNA phosphates correlate reasonably well with some aspects of the Dickerson/Calladine sum function for variation in the helical twist of the oligonucleotides. Correlations between experimentally measured P-O and C-O torsional angles and results from molecular mechanics energy minimization calculations show that these results are consistent with the hypothesis that sequence-specific variations in 31P chemical shifts are attributable to sequence-specific changes in the deoxyribose phosphate backbone. The major structural variation responsible for these 31P shift perturbations appears to be P-O and C-O backbone torsional angles which respond to changes in the local helical structure. Furthermore, 31P chemical shifts and JH3'-P coupling constants both indicate that these backbone torsional angle variations are more permissive at the ends of the double helix than in the middle. Thus 31P NMR spectroscopy and molecular mechanics energy minimization calculations appear to be able to support sequence-specific structural variations along the backbone of the DNA in solution.

Assignment of phosphorus-31 and nonexchangeable proton resonances in a symmetrical 14 base pair lac pseudooperator DNA fragment

The 31P chemical shifts of all 13 phosphates and the chemical shifts of nearly all of the non-exchangeable protons of a symmetrical 14 base pair lac pseudooperator DNA fragment have been assigned by regiospecific labeling with oxygen-17 and two-dimensional NMR techniques. At 22 degrees C, 8 of the 13 phosphorus resonances can distinctly be resolved while the remaining 5 resonances occur in two separate overlapping regions. The 31P chemical shifts of this particular 14 base pair oligonucleotide do not follow the general observation that the more internal the phosphate is located within the oligonucleotide sequence the more upfield the 31P resonance occurs, as shown from other 31P assignment studies. Failure of this general rule is believed to be a result of helical distortions that occur along the oligonucleotide double helix, on the basis of the analysis of Callidine [Callidine, C.R. (1982) J. Mol. Biol. 161, 343-352]. Notable exceptions to the phosphate position relationship are 5'-Py-Pu-3' dinucleotide sequences, which resonate at a lower field strength than expected in agreement with similar results as reported by Ott and Eckstein [Ott, J., & Eckstein, F. (1985) Biochemistry 24, 253]. A reasonable correlation exists between 31P chemical shifts values of the 14-mer and the helical twist sum function of Calladine. The most unusual 31P resonance occurs most upfield in the 31P spectrum, which has been assigned to the second phosphate position (5'-GpT-3') from the 5' end. This unusual chemical shift may be the result of the predicted large helical twist angle that occurs at this position in the 14-mer sequence. Further, it is believed that the large helical twist represents a unique structural feature responsible for optimum binding contact between lac repressor protein and this 14-mer lac pseudooperator segment. Assignments of proton resonances were made from two-dimensional 1H-1H nuclear Overhauser effect (NOESY) connectivities in a sequential manner applicable to right-handed B-DNA, in conjunction with two-dimensional homonuclear and heteronuclear J-correlated spectroscopies (1H-1H COSY and 31P-1H HETCOR). Most nonexchangeable base proton and deoxyribose proton (except for some unresolved H4', H5', and H5" protons) resonances were assigned.