Likelihood and approximate likelihood analyses of genetic structure in a linear habitat: performance and robustness to model mis-specification

We evaluate the performance of maximum likelihood (ML) analysis of allele frequency data in a linear array of populations. The parameters are a mutation rate and either the dispersal rate in a stepping stone model or a dispersal rate and a scale parameter in a geometric dispersal model. An approximate procedure known as maximum product of approximate conditional (PAC) likelihood is found to perform as well as ML. Mis-specification biases may occur because the importance sampling algorithm is formally defined in term of mutation and migration rates scaled by the total size of the population, and this size may differ widely in the statistical model and in reality. As could be expected, ML generally performs well when the statistical model is correctly specified. Otherwise, mutation rate estimates are much closer to mutation probability scaled by number of demes in the statistical model than scaled by number of demes in reality when mutation probability is high and dispersal is most limited. This mis-specification bias actually has practical benefits. However, opposite results are found in opposite conditions. Migration rate estimates show roughly similar trends, but they may not always be easily interpreted as low-bias estimates of dispersal rate under any scaling. Estimation of the dispersal scale parameter is also affected by mis-specification of the number of demes, and the different biases compensate each other in such a way that good estimation of the so-called neighborhood size (or more precisely the product of population density and mean-squared parent-offspring dispersal distance) is achieved. Results congruent with these findings are found in an application to a damselfly data set.     

Optimization of rotamers prior to template minimization improves stability predictions made by computational protein design

Computational protein design (CPD) predictions are highly dependent on the structure of the input template used. However, it is unclear how small differences in template geometry translate to large differences in stability prediction accuracy. Herein, we explored how structural changes to the input template affect the outcome of stability predictions by CPD. To do this, we prepared alternate templates by Rotamer Optimization followed by energy Minimization (ROM) and used them to recapitulate the stability of 84 protein G domain β1 mutant sequences. In the ROM process, side-chain rotamers for wild-type (WT) or mutant sequences are optimized on crystal or nuclear magnetic resonance (NMR) structures prior to template minimization, resulting in alternate structures termed ROM templates. We show that use of ROM templates prepared from sequences known to be stable results predominantly in improved prediction accuracy compared to using the minimized crystal or NMR structures. Conversely, ROM templates prepared from sequences that are less stable than the WT reduce prediction accuracy by increasing the number of false positives. These observed changes in prediction outcomes are attributed to differences in side-chain contacts made by rotamers in ROM templates. Finally, we show that ROM templates prepared from sequences that are unfolded or that adopt a nonnative fold result in the selective enrichment of sequences that are also unfolded or that adopt a nonnative fold, respectively. Our results demonstrate the existence of a rotamer bias caused by the input template that can be harnessed to skew predictions toward sequences displaying desired characteristics.     

西南三江(金沙江、澜沧江和怒江)流域化学风化过程

选择中国西南三江流域作为研究对象,结合水化学及溶解性无机碳碳同位素,借助正演模型,分析河水溶解性组分来源及混入比例,说明硫酸参与区域碳酸盐矿物风化过程及CO2净释放量。结果表明:(1)西南三江流域河水受复杂岩性控制,金沙江河水以Na-Ca-Cl-HCO3为主,与流域上游蒸发盐矿物溶解有关,澜沧江和怒江则以Ca-HCO3为主,显示碳酸盐和硅酸盐矿物溶解的影响;(2)蒸发盐矿物溶解是金沙江干流河水阳离子主要来源,其贡献均值为52%,对澜沧江和怒江河水贡献较小,平均值分别为11%和2%。碳酸盐矿物溶解是澜沧江和怒江干流河水阳离子主要来源,其贡献均值分别为70%和78%,对金沙江干流河水的贡献较小,平均为38%,但对金沙江支流河水的贡献较大,平均为74%。硅酸盐风化产物对金沙江、澜沧江和怒江干流的平均贡献分别为8%、16%和15%,对金沙江支流河水的贡献为15%;(3)金沙江、澜沧江和怒江流域硅酸盐矿物风化速率分别为1.39、3.27和4.27 t·km-2·a-1,CO2消耗速率分别为0.34×105、1.18×105和1.40×105mol·km-2·a-1,碳酸盐矿物风化速率分别为16.93、33.13和33.54 t·km-2·a-1,CO2消耗速率分别为1.40×105、2.47×105和2.22×105mol·km-2·a-1;(4)西南三江流域河水硫酸盐主要来源于硫化物矿物氧化,硫酸参与碳酸盐矿物风化净释放CO2量分别为0.73×105、3.01×105和4.27×105mol·km-2·a-1,高于中国贵州地区,应当在研究区域碳循环以及全球气候变化过程中引起重视。     

A radioisotope-nondependent high-sensitivity method for measuring the activity of glioblastoma-related O-methylguanine DNA methyltransferase

O⁶-Methylguanine DNA methyltransferase (MGMT) cancels the anticancer effect of temozolomide (drug for glioblastoma), which introduces methylation to DNA. Therefore, developing an MGMT inhibitor is a promising strategy for the treatment of this cancer. For this purpose, a sensitive detection method that does not depend on the conventional radioisotope (RI) method was developed. This was realized by a fluorescence-based method that measured the amount of cleavable restriction sites demethylated by the action of MGMT; this method was enhanced by introducing a polymerase chain reaction (PCR) amplification step. As an assay of enzyme activity, 20-fold higher sensitivity (subnanomolar) was attained compared with our and others’ fluorescence-based approaches.     

Evaluation and elimination of inhibitory effects of salts and heavy metal ions on biodegradation of Congo red by Pseudomonas sp. mutant

In this study, it was attempted to evaluate the influences and also recommended some elimination methods for inhibitory effects offered by salts and heavy metal ions. Congo red dye solution treated with mutant Pseudomonas sp. was taken as a model system for study. The salts used in this study are NaCl, CaCl₂ and MgSO₄·7H₂O. Though the growth was inhibited at concentrations above 4 g/l, toleration was achieved by acclimatization process. In case of heavy metal ions, Cr (VI) showed low inhibition up to 500 mg/l of concentration, compared to Zn (II) and Cu (II). It was due to the presence of chromium reductase enzyme which was confirmed by SDS-PAGE. Zn (II) and Cu (II) ion inhibitions were eliminated by chelation with EDTA. The critical ion concentrations obtained as per Han-Levenspiel model for Cr (VI), Zn (II) and Cu (II) were 0.8958, 0.3028 and 0.204 g/l respectively.     

Total mercury in stranded marine mammals from the Oregon and southern Washington coasts

Muscle samples from 105 marine mammals stranded along the Oregon and Washington coasts (2002–2009) were tested for levels of total mercury (THg) by Cold Vapor Atomic Fluorescence Spectrometry. The THg present is in the form of the highly toxic methylmercury. After normalizing tissue to 75% water weight, Steller sea lions and northern elephant seals exhibited the highest mean concentrations of THg followed by harbor seals, harbor porpoises, and California sea lions, 0.34 ± 0.278, 0.34 ± 0.485, 0.21 ± 0.216, 0.17 ± 0.169, and 0.15 ± 0.126 mg/kg normalized wet weight (ww), respectively. The mean normalized values demonstrate limited muscle methylmercury accumulation in these species in the Pacific Northwest. However, actual ww concentrations in some of the stranded carcasses may pose a risk to scavengers. Normalizing muscle mercury concentrations eliminated the variability from desiccation, and allowed for a clearer indication of the amount of mercury the animal accumulated before stranding.     

The Cus efflux system removes toxic ions via a methionine shuttle

Gram-negative bacteria, such as Escherichia coli, frequently utilize tripartite efflux complexes in the resistance-nodulation-cell division (RND) family to expel diverse toxic compounds from the cell. These efflux systems span the entire cell envelope to mediate the phenomenon of bacterial multidrug resistance. The three parts of the efflux complexes are: (1) a membrane fusion protein (MFP) connecting (2) a substrate-binding inner membrane transporter to (3) an outer membrane-anchored channel in the periplasmic space. One such efflux system CusCBA is responsible for extruding biocidal Cu(I) and Ag(I) ions. We recently determined the crystal structures of both the inner membrane transporter CusA and MFP CusB of the CusCBA tripartite efflux system from E. coli. These are the first structures of the heavy-metal efflux (HME) subfamily of the RND efflux pumps. Here, we summarize the structural information of these two efflux proteins and present the accumulated evidence that this efflux system utilizes methionine residues to bind and export Cu(I)/Ag(I). Genetic and structural analyses suggest that the CusA pump is capable of picking up the metal ions from both the periplasm and cytoplasm. We propose a stepwise shuttle mechanism for this pump to extrude metal ions from the cell.     

Historical introgression and the role of selective vs. neutral processes in structuring nuclear genetic variation (AFLP) in a circumpolar marine fish, the capelin (Mallotus villosus)

The capelin (Mallotus villosus) is a widespread marine fish species for which previous work has identified geographically distinct mtDNA clades, the frontiers of which are well within adult and larval dispersal capabilities. Here, we use AFLPs to test for the presence of nuclear gene flow among clades. In addition, we evaluate genetic structuring within one clade, the Northwest Atlantic (NWA). We found that each of the mtDNA clades corresponds with a unique nuclear DNA genetic cluster. Within the NWA clade, we detected individuals with small but significant amounts of genetic ancestry from other clades, likely due to historical introgression. Further support for historical introgression comes from analyses of variance in locus-specific differentiation, which support introgression between some clades and divergence without gene flow between others. Within the NWA, we identified two genetic clusters that correspond to sites in geographically adjacent areas. However, these clusters differ primarily at 'outlier' loci, and a genetic subdivision (K=2) was not supported by genetic clustering programs using neutral loci. Significant neutral F(ST) differentiation was found only between sites that otherwise differed at outlier loci. Thus, these populations may be in the initial stages of 'isolation by adaptation'. These results suggest strong between-clade reproductive isolation despite opportunities for gene flow and support the hypothesis that selection can contribute to divergence in otherwise 'open' systems.     

Total Chemical Synthesis,Assembly of Human Torque Teno Virus Genome

Torque teno virus(TTV) is a nonenveloped virus containing a single-stranded,circular DNA genome of approximately 3.8kb.We completely synthesized the 3 808 nucleotides of the TTV(SANBAN isolate) genome,which contains a hairpin structure and a GC-rich region.More than 100 overlapping oligonucleotides were chemically synthesized and assembled by polymerase chain assembly reaction(PCA),and the synthesis was completed with splicing by overlap extension(SOEing).This study establishes the methodological basis of t...