Wavelet-based parametric functional mapping of developmental trajectories with high-dimensional data

The biological and statistical advantages of functional mapping result from joint modeling of the mean-covariance structures for developmental trajectories of a complex trait measured at a series of time points. While an increased number of time points can better describe the dynamic pattern of trait development, significant difficulties in performing functional mapping arise from prohibitive computational times required as well as from modeling the structure of a high-dimensional covariance matrix. In this article, we develop a statistical model for functional mapping of quantitative trait loci (QTL) that govern the developmental process of a quantitative trait on the basis of wavelet dimension reduction. By breaking an original signal down into a spectrum by taking its averages (smooth coefficients) and differences (detail coefficients), we used the discrete Haar wavelet shrinkage technique to transform an inherently high-dimensional biological problem into its tractable low-dimensional representation within the framework of functional mapping constructed by a Gaussian mixture model. Unlike conventional nonparametric modeling of wavelet shrinkage, we incorporate mathematical aspects of developmental trajectories into the smooth coefficients used for QTL mapping, thus preserving the biological relevance of functional mapping in formulating a number of hypothesis tests at the interplay between gene actions/interactions and developmental patterns for complex phenotypes. This wavelet-based parametric functional mapping has been statistically examined and compared with full-dimensional functional mapping through simulation studies. It holds great promise as a powerful statistical tool to unravel the genetic machinery of developmental trajectories with large-scale high-dimensional data.     

Phosphorylation of synaptic vesicle protein 2 modulates binding to synaptotagmin

Synaptic vesicle protein 2 (SV2) is a component of all synaptic vesicles that is required for normal neurotransmission. Here we report that in intact synaptic terminals SV2 is a phosphoprotein. Phosphopeptide mapping studies indicate that a major site of phosphorylation is located on the cytoplasmic amino terminus. SV2 is phosphorylated on serine and threonine but not on tyrosine residues, indicating that it is a substrate for serine/threonine kinases. Phosphopeptide mapping, in gel kinase assays, and surveys of kinase inhibitors suggest that casein kinase I is a primary SV2 kinase. The amino terminus of SV2 was previously shown to mediate its interaction with synaptotagmin, a calcium-binding protein also required for normal neurotransmission. Comparison of synaptotagmin binding with phosphorylated and unphosphorylated SV2 amino-terminal peptides reveals an increase in binding with phosphorylation. These results suggest that the affinity of SV2 for synaptotagmin is modulated by phosphorylation of SV2.     

The application of next-generation sequencing in the autozygosity mapping of human recessive diseases

Autozygosity, or the inheritance of two copies of an ancestral allele, has the potential to not only reveal phenotypes caused by biallelic mutations in autosomal recessive genes, but to also facilitate the mapping of such mutations by flagging the surrounding haplotypes as tractable runs of homozygosity (ROH), a process known as autozygosity mapping. Since SNPs replaced microsatellites as markers for the purpose of genomewide identification of ROH, autozygosity mapping of Mendelian genes has witnessed a significant acceleration. Historically, successful mapping traditionally required favorable family structure that permits the identification of an autozygous interval that is amenable to candidate gene selection and confirmation by Sanger sequencing. This requirement presented a major bottleneck that hindered the utilization of simplex cases and many multiplex families with autosomal recessive phenotypes. However, the advent of next-generation sequencing that enables massively parallel sequencing of DNA has largely bypassed this bottleneck and thus ushered in an era of unprecedented pace of Mendelian disease gene discovery. The ability to identify a single causal mutation among a massive number of variants that are uncovered by next-generation sequencing can be challenging, but applying autozygosity as a filter can greatly enhance the enrichment process and its throughput. This review will discuss the power of combining the best of both techniques in the mapping of recessive disease genes and offer some tips to troubleshoot potential limitations.     

Computational solvent mapping reveals the importance of local conformational changes for broad substrate specificity in mammalian cytochromes P450

Computational solvent mapping moves small organic molecules as probes around a protein surface, finds favorable binding positions, clusters the conformations, and ranks the clusters on the basis of their average free energy. Prior mapping studies of enzymes, crystallized in either substrate-free or substrate-bound form, have shown that the largest number of solvent probe clusters invariably overlaps in the active site. We have applied this method to five cytochromes P450. As expected, the mapping of two bacterial P450s, P450 cam (CYP101) and P450 BM-3 (CYP102), identified the substrate-binding sites in both ligand-bound and ligand-free P450 structures. However, the mapping finds the active site only in the ligand-bound structures of the three mammalian P450s, 2C5, 2C9, and 2B4. Thus, despite the large cavities seen in the unbound structures of these enzymes, the features required for binding small molecules are formed only in the process of substrate binding. The ability of adjusting their binding sites to substrates that differ in size, shape, and polarity is likely to be responsible for the broad substrate specificity of these mammalian P450s. Similar behavior was seen at "hot spots" of protein-protein interfaces that can also bind small molecules in grooves created by induced fit. In addition, the binding of S-warfarin to P450 2C9 creates a high-affinity site for a second ligand, which may help to explain the prevalence of drug-drug interactions involving this and other mammalian P450s.     

The peptide antibiotic microcin 25 is imported through the TonB pathway and the SbmA protein.

Selection of spontaneous mutants for insensitivity to the peptide antibiotic microcin 25 led to the isolation of five categories of mutants. Phenotypic and mapping studies showed the mutations to be located in the fhuA, exb, tonB, and sbmA genes. The latter encodes a cytoplasmic membrane protein which is also required for the penetration of microcin B17.     

Mapping genes for polygenic disorders: considerations for study design in the complex trait of inflammatory bowel disease

While the methodology for the mapping of Mendelian disorders is well established, the practical and theoretical steps required for successful gene identification in a complex trait are still difficult to predict. A number of analytical models and simulations based on repetitive drawings from predefined statistical distributions are available. To supplement these analytical models, we developed an integrated simulation approach by directly simulating entire populations under a disease model based on epidemiological data. Random mating, nonoverlapping populations and the absence of differential fitness were assumed. Samples were drawn from these homogeneous and heterogeneous populations and analyzed with established analysis tools. We investigated the properties of linkage and association studies in inflammatory bowel disease - modeled as a six-locus polygenic disorder - as an example of this approach. In nonparametric linkage studies, lod scores varied widely, with the median required sample size depending on the locus-specific relative sibling risk. A fine mapping resolution <4 cM was found to require nonparametric lod scores >10. Family-based association studies (TDT test) and case-control studies showed a similar sensitivity and can identify risk loci in populations with moderate levels of linkage disequilibrium in sample sizes of 500-800 triplets. Case-control association studies were prone to false-positive results if applied in heterogeneous populations, with the false-positive rate increasing with sample size because population heterogeneity is detected with increasing power.     

Alumorphs--human DNA polymorphisms detected by polymerase chain reaction using Alu-specific primers

The simultaneous analysis of multiple loci could substantially increase the efficiency of mapping studies. Toward this goal, we used the polymerase chain reaction to amplify multiple DNA fragments originating from dispersed genomic segments that are flanked by Alu repeats. Analysis of different human DNA samples revealed numerous amplification products distinguishable by size, some of which vary between individuals. A family study demonstrated that these polymorphic fragments are inherited in a Mendelian fashion. Because of the ubiquitous distribution of Alu repeats, these markers, called "alumorphs," could be useful for linkage mapping of the human genome. A major advantage of alumorphs is that no prior knowledge of DNA sequence of marker loci is required. This approach may find general application for any genome where interspersed repetitive sequences are found.     

A Functional and Structural Analysis of the Sex Combs Reduced Locus of Drosophila Melanogaster

We have undertaken a developmental genetic analysis of the homeotic gene Sex combs reduced (Scr) of Drosophila melanogaster by examining embryonic and adult phenotypes of mutations affecting Scr gene function. Molecular mapping of Scr breakpoint lesions has defined a segment of greater than 70 kb of DNA necessary for proper Scr gene function. This region is split by the fushi tarazu (ftz) gene, with lesions affecting embryonic Scr function molecularly mapping to the region proximal (5') to ftz and those exhibiting polyphasic semilethality predominantly mapping distal (3') to ftz. Gain-of-function mutations are associated with genomic rearrangements and map throughout the Scr locus. Our analysis has revealed that the Scr locus encompasses genetic elements that are responsible for functions in both the embryonic and larval to adult periods of development. From these studies, we conclude that Scr is a complex genetic locus with an extensive regulatory region that directs functions required for normal head and thoracic development in both the embryo and the adult and that the regulation of Scr during these two periods is distinct.     

A composite method for mapping quantitative trait loci without interference of female achiasmatic and gender effects in silkworm,Bombyx mori

The silkworm, Bombyx mori, is an economically important insect that was domesticated more than 5000 years ago. Its major economic traits focused on by breeders are quantitative traits, and an accurate and efficient QTL mapping method is necessary to explore their genetic architecture. However, current widely used QTL mapping models are not well suited for silkworm because they ignore female achiasmate and gender effects. In this study, we propose a composite method combining rational population selection and special mapping methods to map QTL in silkworm. By determining QTL for cocoon shell weight (CSW), we demonstrated the effectiveness of this method. In the CSW mapping process, only 56 markers were used and five loci or chromosomes were detected, more than in previous studies. Additionally, loci on chromosomes 1 and 11 dominated and accounted for 35.10% and 15.03% of the phenotypic variance respectively. Unlike previous studies, epistasis was detected between loci on chromosomes 11 and 22. These mapping results demonstrate the power and convenience of this method for QTL mapping in silkworm, and this method may inspire the development of similar approaches for other species with special genetic characteristics.     

Quality-controlled R-loop meta-analysis reveals the characteristics of R-loop consensus regions

R-loops are three-stranded nucleic acid structures formed from the hybridization of RNA and DNA. While the pathological consequences of R-loops have been well-studied to date, the locations, classes, and dynamics of physiological R-loops remain poorly understood. R-loop mapping studies provide insight into R-loop dynamics, but their findings are challenging to generalize. This is due to the narrow biological scope of individual studies, the limitations of each mapping modality, and, in some cases, poor data quality. In this study, we reprocessed 810 R-loop mapping datasets from a wide array of biological conditions and mapping modalities. From this data resource, we developed an accurate R-loop data quality control method, and we reveal the extent of poor-quality data within previously published studies. We then identified a set of high-confidence R-loop mapping samples and used them to define consensus R-loop sites called ‘R-loop regions’ (RL regions). In the process, we identified a stark divergence between RL regions detected by S9.6 and dRNH-based mapping methods, particularly with respect to R-loop size, location, and colocalization with RNA binding factors. Taken together, this work provides a much-needed method to assess R-loop data quality and offers novel context regarding the differences between dRNH- and S9.6-based R-loop mapping approaches.     

<Emphasis Type="Italic">O</Emphasis>-GlcNAcylation of the <Emphasis Type="Italic">Plum pox virus</Emphasis> capsid protein catalyzed by SECRET AGENT: characterization of <Emphasis Type="Italic">O</Emphasis>-GlcNAc sites by electron transfer dissociation mass spectrometry

The capsid protein of Plum pox virus (PPV-CP) is modified with O-linked β-N-acetylglucosamine (O-GlcNAc). In Arabidopsis thaliana this modification is made by an O-GlcNAc transferase named SECRET AGENT (SEC). Modification of PPV-CP by SEC is hypothesized to have a direct role in the infection process, because virus titer and rate of spread are reduced in SEC mutants. Previous studies used deletion mapping and site-directed mutagenesis to identify four O-GlcNAc sites on the capsid protein that are modified by Escherichia coli-expressed SEC. The infection process was not affected when two of these sites were mutated suggesting that O-GlcNAcylation of these sites does not have a significant role in the infection process or that a subset of the modifications is sufficient. Since it is possible that the mutational mapping approach missed or incorrectly identified O-GlcNAc sites, the modifications produced by E. coli-expressed SEC were characterized using mass spectrometry. O-GlcNAcylated peptides were enzymatically tagged with galactose, the products were enriched on immobilized Ricinus communis agglutinin I and sequenced by electron transfer dissociation (ETD) mass spectrometry. Five O-GlcNAc sites on PPV-CP were identified. Two of these sites were not identified in by the previous mutational mapping. In addition, one site previously predicted by mutation mapping was not detected, but modification of this site was not supported when the mutation mapping was repeated. This study suggests that mapping modification sites by ETD mass spectrometry is more comprehensive and accurate than mutational mapping.     

Identification of novel, less toxic PTP-LAR inhibitors using in silico strategies: pharmacophore modeling, SADMET-based virtual screening and docking

Human leukocyte antigen-related (PTP-LAR) is a receptor-like transmembrane phosphatase and a potential target for diabetes, obesity and cancer. In the present study, a sequence of in silico strategies (pharmacophore mapping, a 3D database searching, SADMET screening, and docking and toxicity studies) was performed to identify eight novel nontoxic PTP-LAR inhibitors. Twenty different pharmacophore hypotheses were generated using two methods; the best (hypothesis 2) consisted of three hydrogen-bond acceptor (A), one ring aromatic (R), and one hydrophobic aliphatic (Z) features. This hypothesis was used to screen molecules from several databases, such as Specs, IBS, MiniMaybridge, NCI, and an in-house PTP inhibitor database. In order to overcome the general bioavailability problem associated with phosphatases, the hits obtained were filtered by Lipinski’s rule of five and SADMET properties and validated by molecular docking studies using the available crystal structure 1LAR. These docking studies suggested the ligand binding pattern and interactions required for LAR inhibition. The docking analysis also revealed that sulfonylurea derivatives with an isoquinoline or naphthalene scaffold represent potential LAR drugs. The screening protocol was further validated using ligand pharmacophore mapping studies, which showed that the abovementioned interactions are indeed crucial and that the screened molecules can be presumed to possess potent inhibitory activities.     

Clustering and variance maps for cryo-electron tomography using wedge-masked differences

Cryo-electron tomography provides 3D imaging of frozen hydrated biological samples with nanometer resolution. Reconstructed volumes suffer from low signal-to-noise-ratio (SNR)(1) and artifacts caused by systematically missing tomographic data. Both problems can be overcome by combining multiple subvolumes with varying orientations, assuming they contain identical structures. Clustering (unsupervised classification) is required to ensure or verify population homogeneity, but this process is complicated by the problems of poor SNR and missing data, the factors that led to consideration of multiple subvolumes in the first place. Here, we describe a new approach to clustering and variance mapping in the face of these difficulties. The combined subvolume is taken as an estimate of the true subvolume, and the effect of missing data is computed for individual subvolumes. Clustering and variance mapping then proceed based on differences between expected and observed subvolumes. We show that this new method is faster and more accurate than two current, widely used techniques.     

High-throughput RFLP genotyping method for large genomes based on a chemiluminescent detection system

Genetic linkage mapping based on RFLPs is a valuable genomics tool for studying organisms with no genome sequence information. However, the generally used Southern hybridization method based on the radioisotope³²P is not ideal for genotyping large mapping populations. We have overcome limitations of the alternative chemiluminescent detection system and developed a high-throughput RFLP genotyping method suitable for large-scale mapping studies of large genomes. Important elements in our process are PCR labeling of probes, complete removal of post-PCR unincorporated nucleotides via column-based purification methods, use of a 1:4 DIG-[11]-dUTP:dTTP ratio, and using a rocker instead of an orbital shaker during hybridization and post-hybridization processing of membranes. Using this method, we mapped the large genome of the homosporous fern speciesCeratopteris richardii by genotyping a mapping population of 513 doubled haploid line (DHL) progeny of a cross between two completely homozygous parental lines. Our genotyping method can robustly detect sub-picogram quantities of DNA fragments from a large number of samples and can be applied to linkage mapping studies of other organisms with large genomes.     

Rice Yellow Mottle Virus stress responsive genes from susceptible and tolerant rice genotypes

Background  

The effects of viral infection involve concomitant plant gene variations and cellular changes. A simple system is required to assess the complexity of host responses to viral infection. The genome of the Rice yellow mottle virus (RYMV) is a single-stranded RNA with a simple organisation. It is the most well-known monocotyledon virus model. Several studies on its biology, structure and phylogeography have provided a suitable background for further genetic studies. 12 rice chromosome sequences are now available and provide strong support for genomic studies, particularly physical mapping and gene identification.