Conformational solution studies of the SDS micelle-bound 11-28 fragment of two Alzheimer's beta-amyloid variants (E22K and A21G) using CD, NMR, and MD techniques

The beta-amyloid (Abeta) is the major peptide constituent of neuritic plaques in Alzheimer's disease (AD) and its aggregation is believed to play a central role in the pathogenesis of the disease. Naturally occurring mutations resulting in changes in the Abeta sequence (pos. 21-23) are associated with familial AD-like diseases with extensive cerebrovascular pathology. It was proved that the mutations alter the aggregation ability of Abeta and its neurotoxicity. Among five mutations at positions 21-23 there are two mutations with distinct clinical characteristics and potentially distinct pathogenic mechanism-the Italian (E22K) and the Flemish (A21G) mutations. In our studies we have examined the structures of the 11-28 fragment of the Italian and Flemish Abeta variants. The fragment was chosen because it has been shown to be the most important for amyloid fibril formation. The detailed structure of both variants Abeta(11-28) was determined using CD, 2D NMR, and molecular dynamics techniques under water-SDS micelle conditions. The NMR analysis revealed two distinct sets of proton resonances for the peptides. The studies of both peptides pointed out the existence of well-defined alpha-helical conformation in the Italian mutant, whereas the Flemish was found to be unstructured with the possibility of a bent structure in the central part of the peptide.     

Quality control and integration of genotypes from two calling pipelines for whole genome sequence data in the Alzheimer's disease sequencing project

The Alzheimer's Disease Sequencing Project (ADSP) performed whole genome sequencing (WGS) of 584 subjects from 111 multiplex families at three sequencing centers. Genotype calling of single nucleotide variants (SNVs) and insertion-deletion variants (indels) was performed centrally using GATK-HaplotypeCaller and Atlas V2. The ADSP Quality Control (QC) Working Group applied QC protocols to project-level variant call format files (VCFs) from each pipeline, and developed and implemented a novel protocol, termed “consensus calling,” to combine genotype calls from both pipelines into a single high-quality set. QC was applied to autosomal bi-allelic SNVs and indels, and included pipeline-recommended QC filters, variant-level QC, and sample-level QC. Low-quality variants or genotypes were excluded, and sample outliers were noted. Quality was assessed by examining Mendelian inconsistencies (MIs) among 67 parent-offspring pairs, and MIs were used to establish additional genotype-specific filters for GATK calls. After QC, 578 subjects remained. Pipeline-specific QC excluded ~12.0% of GATK and 14.5% of Atlas SNVs. Between pipelines, ~91% of SNV genotypes across all QCed variants were concordant; 4.23% and 4.56% of genotypes were exclusive to Atlas or GATK, respectively; the remaining ~0.01% of discordant genotypes were excluded. For indels, variant-level QC excluded ~36.8% of GATK and 35.3% of Atlas indels. Between pipelines, ~55.6% of indel genotypes were concordant; while 10.3% and 28.3% were exclusive to Atlas or GATK, respectively; and ~0.29% of discordant genotypes were. The final WGS consensus dataset contains 27,896,774 SNVs and 3,133,926 indels and is publicly available.     

Homozygosity Haplotype and Whole-Exome Sequencing Analysis to Identify Potentially Functional Rare Variants Involved in Multiple Sclerosis among Sardinian Families

Multiple Sclerosis (MS) is a complex multifactorial autoimmune disease, whose sex- and age-adjusted prevalence in Sardinia (Italy) is among the highest worldwide. To date, 233 loci were associated with MS and almost 20% of risk heritability is attributable to common genetic variants, but many low-frequency and rare variants remain to be discovered. Here, we aimed to contribute to the understanding of the genetic basis of MS by investigating potentially functional rare variants. To this end, we analyzed thirteen multiplex Sardinian families with Immunochip genotyping data. For five families, Whole Exome Sequencing (WES) data were also available. Firstly, we performed a non-parametric Homozygosity Haplotype analysis for identifying the Region from Common Ancestor (RCA). Then, on these potential disease-linked RCA, we searched for the presence of rare variants shared by the affected individuals by analyzing WES data. We found: (i) a variant (43181034 T > G) in the splicing region on exon 27 of CUL9; (ii) a variant (50245517 A > C) in the splicing region on exon 16 of ATP9A; (iii) a non-synonymous variant (43223539 A > C), on exon 9 of TTBK1; (iv) a non-synonymous variant (42976917 A > C) on exon 9 of PPP2R5D; and v) a variant (109859349-109859354) in 3′UTR of MYO16.     

Phylogenetic position of Rickettsia tsutsugamushi and the relationship among its antigenic variants by analyses of 16S rRNA gene sequences

Abstract The 16S rRNA gene sequences of Rickettsia tsutsugamushi and Rickettsia sibirica were determined by PCR and DNA sequencing. Phylogenetic analysis revealed that R. sibirica is positioned in a cluster of the genus Rickettsia with a similarity value of 98.1–99.6%, whereas R. tsutsugamushi is located apart from the cluster with a similarity value of 90.2–90.6%. This evidence suggests that R. tsutsugamushi should be excluded taxonomically from the genus Rickettsia . The phylogenetic classification of six antigenic variants in R. tsutsugamushi moderately reflected their antigenic relationship known in closely and distantly related strains.     

The Slowing Rate of CpG Depletion in SARS-CoV-2 Genomes Is Consistent with Adaptations to the Human Host

Depletion of CpG dinucleotides in severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) genomes has been linked to virus evolution, host-switching, virus replication, and innate immune responses. Temporal variations, if any, in the rate of CpG depletion during virus evolution in the host remain poorly understood. Here, we analyzed the CpG content of over 1.4 million full-length SARS-CoV-2 genomes representing over 170 million documented infections during the first 17 months of the pandemic. Our findings suggest that the extent of CpG depletion in SARS-CoV-2 genomes is modest. Interestingly, the rate of CpG depletion is highest during early evolution in humans and it gradually tapers off, almost reaching an equilibrium; this is consistent with adaptations to the human host. Furthermore, within the coding regions, CpG depletion occurs predominantly at codon positions 2-3 and 3-1. Loss of ZAP (Zinc-finger antiviral protein)-binding motifs in SARS-CoV-2 genomes is primarily driven by the loss of the terminal CpG within the motifs. Nonetheless, majority of the CpG depletion in SARS-CoV-2 genomes occurs outside ZAP-binding motifs. SARS-CoV-2 genomes selectively lose CpGs-motifs from a U-rich context; this may help avoid immune recognition by TLR7. SARS-CoV-2 alpha-, beta-, and delta-variants of concern have reduced CpG content compared to sequences from the beginning of the pandemic. In sum, we provide evidence that the rate of CpG depletion in virus genomes is not uniform and it greatly varies over time and during adaptations to the host. This work highlights how temporal variations in selection pressures during virus adaption may impact the rate and the extent of CpG depletion in virus genomes.     

Emergence,evolution, and vaccine production approaches of SARS-CoV-2 virus: Benefits of getting vaccinated and common questions

The emergence of coronavirus disease 2019 (COVID-19) pandemic in Wuhan city, China at the end of 2019 made it urgent to identify the origin of the causal pathogen and its molecular evolution, to appropriately design an effective vaccine. This study analyzes the evolutionary background of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2 or SARS-2) in accordance with its close relative SARS-CoV (SARS-1), which was emerged in 2002. A comparative genomic and proteomic study was conducted on SARS-2, SARS-1, and Middle East respiratory syndrome coronavirus (MERS), which was emerged in 2012. In silico analysis inferred the genetic variability among the tested viruses. The SARS-1 genome harbored 11 genes encoding 12 proteins, while SARS-2 genome contained only 10 genes encoding for 10 proteins. MERS genome contained 11 genes encoding 11 proteins. The analysis also revealed a slight variation in the whole genome size of SARS-2 comparing to its siblings resulting from sequential insertions and deletions (indels) throughout the viral genome particularly ORF1AB, spike, ORF10 and ORF8. The effective indels were observed in the gene encoding the spike protein that is responsible for viral attachment to the angiotensin-converting enzyme 2 (ACE2) cell receptor and initiating infection. These indels are responsible for the newly emerging COVID-19 variants αCoV, βCoV, γCoV and δCoV. Nowadays, few effective COVID-19 vaccines developed based on spike (S) glycoprotein were approved and become available worldwide. Currently available vaccines can relatively prevent the spread of COVID-19 and suppress the disease. The traditional (killed or attenuated virus vaccine and antibody-based vaccine) and innovated vaccine production technologies (RNA- and DNA-based vaccines and viral vectors) are summarized in this review. We finally highlight the most common questions related to COVID-19 disease and the benefits of getting vaccinated.     

基于DNA弯曲度的H2A.Z核小体定位与修饰研究

在真核生物染色质中,H2A.Z是高度保守的组蛋白变异体,与转录调控、基因组的稳定性密切相关。为了探讨组蛋白修饰、DNA弯曲度与H2A.Z核小体定位三者之间的关联,在得到实验所测的相关数据后,利用MINE算法并结合皮尔逊相关系数在酵母全基因组的转录起始位点周围探讨了三者间的线性与非线性关系。其中MIC算法可以定量的得出数据之间关联度大小的值,用于衡量数据之间是否存在着关联,而皮尔逊相关系数则用于检查是否为线性关联。结果除了发现大部分组蛋白修饰种类和核小体定位之间存在着线性关联外,还探测到有两种组蛋白修饰数据(H4ac修饰与GCN4修饰)和核小体定位数据之间存在着以往未发现的非线性关系(大致呈正余弦函数),并从数据的生物背景(组蛋白修饰与核小体位置)上探讨了出现非线性现象的原因。     

Functional implications of structural differences between variants A and B of bovine beta-lactoglobulin

The structure of the trigonal crystal form of bovine beta-lactoglobulin variant B at pH 7.1 has been determined by X-ray diffraction methods at a resolution of 2.22 A and refined to values for R and Rfree of 0.239 and 0.286, respectively. By comparison with the structure of the trigonal crystal form of bovine beta-lactoglobulin variant A at pH 7.1, which was determined previously [Qin BY et al., 1998, Biochemistry 37:14014-14023], the structural consequences of the sequence differences D64G and V118A of variants A and B, respectively, have been investigated. Only minor differences in the core calyx structure occur. In the vicinity of the mutation site D64G on loop CD (residues 61-67), there are small changes in main-chain conformation, whereas the substitution V118A on beta-strand H is unaccompanied by changes in the surrounding structure, thereby creating a void volume and weakened hydrophobic interactions with a consequent loss of thermal stability relative to variant A. A conformational difference is found for the loop EF, implicated in the pH-dependent conformational change known as the Tanford transition, but it is not clear whether this reflects differences intrinsic to the variants in solution or differences in crystallization.