Whole-exome sequencing reveals damaging gene variants associated with hypoalphalipoproteinemia.

Low levels of high density lipoprotein-cholesterol (HDL-C) are associated with an elevated risk of arteriosclerotic coronary heart disease. Heritability of HDL-C levels is high. In this research discovery study, we used whole-exome sequencing to identify damaging gene variants that may play significant roles in determining HDL-C levels. We studied 204 individuals with a mean HDL-C level of 27.8 ± 6.4 mg/dl (range: 4-36 mg/dl). Data were analyzed by statistical gene burden testing and by filtering against candidate gene lists. We found 120 occurrences of probably damaging variants (116 heterozygous; four homozygous) among 45 of 104 recognized HDL candidate genes. Those with the highest prevalence of damaging variants were ABCA1 (n = 20), STAB1 (n = 9), OSBPL1A (n = 8), CPS1 (n = 8), CD36 (n = 7), LRP1 (n = 6), ABCA8 (n = 6), GOT2 (n = 5), AMPD3 (n = 5), WWOX (n = 4), and IRS1 (n = 4). Binomial analysis for damaging missense or loss-of-function variants identified the ABCA1 and LDLR genes at genome-wide significance. In conclusion, whole-exome sequencing of individuals with low HDL-C showed the burden of damaging rare variants in the ABCA1 and LDLR genes is particularly high and revealed numerous occurrences in HDL candidate genes, including many genes identified in genome-wide association study reports. Many of these genes are involved in cancer biology, which accords with epidemiologic findings of the association of HDL deficiency with increased risk of cancer, thus presenting a new area of interest in HDL genomics.

Evaluating the quality of the 1000 genomes project data.

BACKGROUND: Data from the 1000 Genomes project is quite often used as a reference for human genomic analysis. However, its accuracy needs to be assessed to understand the quality of predictions made using this reference. We present here an assessment of the genotyping, phasing, and imputation accuracy data in the 1000 Genomes project. We compare the phased haplotype calls from the 1000 Genomes project to experimentally phased haplotypes for 28 of the same individuals sequenced using the 10X Genomics platform. RESULTS: We observe that phasing and imputation for rare variants are unreliable, which likely reflects the limited sample size of the 1000 Genomes project data. Further, it appears that using a population specific reference panel does not improve the accuracy of imputation over using the entire 1000 Genomes data set as a reference panel. We also note that the error rates and trends depend on the choice of definition of error, and hence any error reporting needs to take these definitions into account. CONCLUSIONS: The quality of the 1000 Genomes data needs to be considered while using this database for further studies. This work presents an analysis that can be used for these assessments.

A hybrid approach for de novo human genome sequence assembly and phasing.

Despite tremendous progress in genome sequencing, the basic goal of producing a phased (haplotype-resolved) genome sequence with end-to-end contiguity for each chromosome at reasonable cost and effort is still unrealized. In this study, we describe an approach to performing de novo genome assembly and experimental phasing by integrating the data from Illumina short-read sequencing, 10X Genomics linked-read sequencing, and BioNano Genomics genome mapping to yield a high-quality, phased, de novo assembled human genome.

One-Pot Chemoenzymatic Cascade for Labeling of the Epigenetic Marker 5-Hydroxymethylcytosine.

The epigenetic DNA modification 5-hydroxymethylcytosine (5-hmC) is important for the regulation of gene expression during development and in tumorigenesis. 5-hmC can be selectively glycosylated by T4 ß-glucosyltransferase (ß-GT); introduction of an azide on the attached sugar provides a chemical handle for isolation or fluorescent tagging of 5-hmC residues by click chemistry. This approach has not been broadly adopted because of the challenging synthesis and limited commercial availability of the glycosylation substrate, 6-deoxy-6-azido-α-D-glucopyranoside. We report the enzyme-assisted synthesis of this precursor by the uridylyltransferase from Pasteurella multocida (PmGlmU). We were able to directly label 5-hmC in genomic DNA by an enzymatic cascade involving successive action of PmGlmU and ß-GT. This is a facile and cost-effective one-pot chemoenzymatic methodology for 5-hmC analysis.

A blood-brain barrier (BBB) disrupter is also a potent α-synuclein (α-syn) aggregation inhibitor: a novel dual mechanism of mannitol for the treatment of Parkinson disease (PD).

The development of disease-modifying therapy for Parkinson disease has been a main drug development challenge, including the need to deliver the therapeutic agents to the brain. Here, we examined the ability of mannitol to interfere with the aggregation process of α-synuclein in vitro and in vivo in addition to its blood-brain barrier-disrupting properties. Using in vitro studies, we demonstrated the effect of mannitol on α-synuclein aggregation. Although low concentration of mannitol inhibited the formation of fibrils, high concentration significantly decreased the formation of tetramers and high molecular weight oligomers and shifted the secondary structure of α-synuclein from α-helical to a different structure, suggesting alternative potential pathways for aggregation. When administered to a Parkinson Drosophila model, mannitol dramatically corrected its behavioral defects and reduced the amount of α-synuclein aggregates in the brains of treated flies. In the mThy1-human α-synuclein transgenic mouse model, a decrease in α-synuclein accumulation was detected in several brain regions following treatment, suggesting that mannitol promotes α-synuclein clearance in the cell bodies. It appears that mannitol has a general neuroprotective effect in the transgenic treated mice, which includes the dopaminergic system. We therefore suggest mannitol as a basis for a dual mechanism therapeutic agent for the treatment of Parkinson disease.

Quantitative structure-activity relationship analysis of ß-amyloid aggregation inhibitors.

Inhibiting the aggregation process of the ß-amyloid peptide is a promising strategy in treating Alzheimer's disease. In this work, we have collected a dataset of 80 small molecules with known inhibition levels and utilized them to develop two comprehensive quantitative structure-activity relationship models: a Bayesian model and a decision tree model. These models have exhibited high predictive accuracy: 87% of the training and test sets using the Bayesian model and 89 and 93% of the training and test sets, respectively, by the decision tree model. Subsequently these models were used to predict the activities of several new potential ß-amyloid aggregation inhibitors and these predictions were indeed validated by in vitro experiments. Key chemical features correlated with the inhibition ability were identified. These include the electro-topological state of carbonyl groups, AlogP and the number of hydrogen bond donor groups. The results demonstrate the feasibility of the developed models as tools for rapid screening, which could help in the design of novel potential drug candidates for Alzheimer's disease.

Genomic Variation and Recent Population Histories of Spotted (Strix occidentalis) and Barred (Strix varia) Owls.

Spotted owls (SOs, Strix occidentalis) are a flagship species inhabiting old-growth forests in western North America. In recent decades, their populations have declined due to ongoing reductions in suitable habitat caused by logging, wildfires, and competition with the congeneric barred owl (BO, Strix varia). The northern spotted owl (S. o. caurina) has been listed as "threatened" under the Endangered Species Act since 1990. Here, we use an updated SO genome assembly along with 51 high-coverage whole-genome sequences to examine population structure, hybridization, and recent changes in population size in SO and BO. We found that potential hybrids identified from intermediate plumage morphology were a mixture of pure BO, F1 hybrids, and F1 × BO backcrosses. Also, although SO underwent a population bottleneck around the time of the Pleistocene-Holocene transition, their population sizes rebounded and show no evidence of any historical (i.e., 100-10,000 years ago) population decline. This suggests that the current decrease in SO abundance is due to events in the past century. Finally, we estimate that western and eastern BOs have been genetically separated for thousands of years, instead of the previously assumed recent (i.e., <150 years) divergence. Although this result is surprising, it is unclear where the ancestors of western BO lived after the separation. In particular, although BO may have colonized western North America much earlier than the first recorded observations, it is also possible that the estimated divergence time reflects unsampled BO population structure within central or eastern North America.

Application of full-genome analysis to diagnose rare monogenic disorders.

Current genetic testenhancer and narrows the diagnostic intervals for rare diseases provide a diagnosis in only a modest proportion of cases. The Full-Genome Analysis method, FGA, combines long-range assembly and whole-genome sequencing to detect small variants, structural variants with breakpoint resolution, and phasing. We built a variant prioritization pipeline and tested FGA's utility for diagnosis of rare diseases in a clinical setting. FGA identified structural variants and small variants with an overall diagnostic yield of 40% (20 of 50 cases) and 35% in exome-negative cases (8 of 23 cases), 4 of these were structural variants. FGA detected and mapped structural variants that are missed by short reads, including non-coding duplication, and phased variants across long distances of more than 180 kb. With the prioritization algorithm, longer DNA technologies could replace multiple tests for monogenic disorders and expand the range of variants detected. Our study suggests that genomes produced from technologies like FGA can improve variant detection and provide higher resolution genome maps for future application.

Accurate assembly of the olive baboon (Papio anubis) genome using long-read and Hi-C data.

BACKGROUND: Baboons are a widely used nonhuman primate model for biomedical, evolutionary, and basic genetics research. Despite this importance, the genomic resources for baboons are limited. In particular, the current baboon reference genome Panu_3.0 is a highly fragmented, reference-guided (i.e., not fully de novo) assembly, and its poor quality inhibits our ability to conduct downstream genomic analyses. FINDINGS: Here we present a de novo genome assembly of the olive baboon (Papio anubis) that uses data from several recently developed single-molecule technologies. Our assembly, Panubis1.0, has an N50 contig size of ∼1.46 Mb (as opposed to 139 kb for Panu_3.0) and has single scaffolds that span each of the 20 autosomes and the X chromosome. CONCLUSIONS: We highlight multiple lines of evidence (including Bionano Genomics data, pedigree linkage information, and linkage disequilibrium data) suggesting that there are several large assembly errors in Panu_3.0, which have been corrected in Panubis1.0.

The Driver of Extreme Human-Specific Olduvai Repeat Expansion Remains Highly Active in the Human Genome.

Sequences encoding Olduvai protein domains (formerly DUF1220) show the greatest human lineage-specific increase in copy number of any coding region in the genome and have been associated, in a dosage-dependent manner, with brain size, cognitive aptitude, autism, and schizophrenia. Tandem intragenic duplications of a three-domain block, termed the Olduvai triplet, in four NBPF genes in the chromosomal 1q21.1-0.2 region, are primarily responsible for the striking human-specific copy number increase. Interestingly, most of the Olduvai triplets are adjacent to, and transcriptionally coregulated with, three human-specific NOTCH2NL genes that have been shown to promote cortical neurogenesis. Until now, the underlying genomic events that drove the Olduvai hyperamplification in humans have remained unexplained. Here, we show that the presence or absence of an alternative first exon of the Olduvai triplet perfectly discriminates between amplified (58/58) and unamplified (0/12) triplets. We provide sequence and breakpoint analyses that suggest the alternative exon was produced by an nonallelic homologous recombination-based mechanism involving the duplicative transposition of an existing Olduvai exon found in the CON3 domain, which typically occurs at the C-terminal end of NBPF genes. We also provide suggestive in vitro evidence that the alternative exon may promote instability through a putative G-quadraplex (pG4)-based mechanism. Lastly, we use single-molecule optical mapping to characterize the intragenic structural variation observed in NBPF genes in 154 unrelated individuals and 52 related individuals from 16 families and show that the presence of pG4-containing Olduvai triplets is strongly correlated with high levels of Olduvai copy number variation. These results suggest that the same driver of genomic instability that allowed the evolutionarily recent, rapid, and extreme human-specific Olduvai expansion remains highly active in the human genome.

Three patients with homozygous familial hypercholesterolemia: Genomic sequencing and kindred analysis.

BACKGROUND: Homozygous Familial Hypercholesterolemia (HoFH) is an inherited recessive condition associated with extremely high levels of low-density lipoprotein (LDL) cholesterol in affected individuals. It is usually caused by homozygous or compound heterozygous functional mutations in the LDL receptor (LDLR). A number of mutations causing FH have been reported in literature and such genetic heterogeneity presents great challenges for disease diagnosis. OBJECTIVE: We aim to determine the likely genetic defects responsible for three cases of pediatric HoFH in two kindreds. METHODS: We applied whole exome sequencing (WES) on the two probands to determine the likely functional variants among candidate FH genes. We additionally applied 10x Genomics (10xG) Linked-Reads whole genome sequencing (WGS) on one of the kindreds to identify potentially deleterious structural variants (SVs) underlying HoFH. A PCR-based screening assay was also established to detect the LDLR structural variant in a cohort of 641 patients with elevated LDL. RESULTS: In the Caucasian kindred, the FH homozygosity can be attributed to two compound heterozygous LDLR damaging variants, an exon 12 p.G592E missense mutation and a novel 3kb exon 1 deletion. By analyzing the 10xG phased data, we ascertained that this deletion allele was most likely to have originated from a Russian ancestor. In the Mexican kindred, the strikingly elevated LDL cholesterol level can be attributed to a homozygous frameshift LDLR variant p.E113fs. CONCLUSIONS: While the application of WES can provide a cost-effective way of identifying the genetic causes of FH, it often lacks sensitivity for detecting structural variants. Our finding of the LDLR exon 1 deletion highlights the broader utility of Linked-Read WGS in detecting SVs in the clinical setting, especially when HoFH patients remain undiagnosed after WES.

Genome maps across 26 human populations reveal population-specific patterns of structural variation.

Large structural variants (SVs) in the human genome are difficult to detect and study by conventional sequencing technologies. With long-range genome analysis platforms, such as optical mapping, one can identify large SVs (>2 kb) across the genome in one experiment. Analyzing optical genome maps of 154 individuals from the 26 populations sequenced in the 1000 Genomes Project, we find that phylogenetic population patterns of large SVs are similar to those of single nucleotide variations in 86% of the human genome, while ~2% of the genome has high structural complexity. We are able to characterize SVs in many intractable regions of the genome, including segmental duplications and subtelomeric, pericentromeric, and acrocentric areas. In addition, we discover ~60 Mb of non-redundant genome content missing in the reference genome sequence assembly. Our results highlight the need for a comprehensive set of alternate haplotypes from different populations to represent SV patterns in the genome.

De novo human genome assemblies reveal spectrum of alternative haplotypes in diverse populations.

The human reference genome is used extensively in modern biological research. However, a single consensus representation is inadequate to provide a universal reference structure because it is a haplotype among many in the human population. Using 10× Genomics (10×G) "Linked-Read" technology, we perform whole genome sequencing (WGS) and de novo assembly on 17 individuals across five populations. We identify 1842 breakpoint-resolved non-reference unique insertions (NUIs) that, in aggregate, add up to 2.1 Mb of so far undescribed genomic content. Among these, 64% are considered ancestral to humans since they are found in non-human primate genomes. Furthermore, 37% of the NUIs can be found in the human transcriptome and 14% likely arose from Alu-recombination-mediated deletion. Our results underline the need of a set of human reference genomes that includes a comprehensive list of alternative haplotypes to depict the complete spectrum of genetic diversity across populations.

Light-emitting self-assembled peptide nucleic acids exhibit both stacking interactions and Watson-Crick base pairing.

The two main branches of bionanotechnology involve the self-assembly of either peptides or DNA. Peptide scaffolds offer chemical versatility, architectural flexibility and structural complexity, but they lack the precise base pairing and molecular recognition available with nucleic acid assemblies. Here, inspired by the ability of aromatic dipeptides to form ordered nanostructures with unique physical properties, we explore the assembly of peptide nucleic acids (PNAs), which are short DNA mimics that have an amide backbone. All 16 combinations of the very short di-PNA building blocks were synthesized and assayed for their ability to self-associate. Only three guanine-containing di-PNAs-CG, GC and GG-could form ordered assemblies, as observed by electron microscopy, and these di-PNAs efficiently assembled into discrete architectures within a few minutes. The X-ray crystal structure of the GC di-PNA showed the occurrence of both stacking interactions and Watson-Crick base pairing. The assemblies were also found to exhibit optical properties including voltage-dependent electroluminescence and wide-range excitation-dependent fluorescence in the visible region.

The influence of chemical chaperones on enzymatic activity under thermal and chemical stresses: common features and variation among diverse chemical families.

Molecular and chemical chaperones are key components of the two main mechanisms that ensure structural stability and activity under environmental stresses. Yet, chemical chaperones are often regarded only as osmolytes and their role beyond osmotic regulation is not fully understood. Here, we systematically studied a large group of chemical chaperones, representatives of diverse chemical families, for their protective influence under either thermal or chemical stresses. Consistent with previous studies, we observed that in spite of the structural similarity between sugars and sugar alcohols, they have an apparent difference in their protective potential. Our results support the notion that the protective activity is mediated by the solvent and the presence of water is essential. In the current work we revealed that i) polyols and sugars have a completely different profile of protective activity toward trifluoroethanol and thermal stress; ii) minor changes in solvent composition that do not affect enzyme activity, yet have a great effect on the ability of osmolytes to act as protectants and iii) increasing the number of active groups of carbohydrates makes them better protectants while increasing the number of active groups of methylamines does not, as revealed by attempts to synthesize de novo designed methylamines with multiple functional groups.

Toward single-molecule optical mapping of the epigenome.

The past decade has seen an explosive growth in the utilization of single-molecule techniques for the study of complex systems. The ability to resolve phenomena otherwise masked by ensemble averaging has made these approaches especially attractive for the study of biological systems, where stochastic events lead to inherent inhomogeneity at the population level. The complex composition of the genome has made it an ideal system to study at the single-molecule level, and methods aimed at resolving genetic information from long, individual, genomic DNA molecules have been in use for the last 30 years. These methods, and particularly optical-based mapping of DNA, have been instrumental in highlighting genomic variation and contributed significantly to the assembly of many genomes including the human genome. Nanotechnology and nanoscopy have been a strong driving force for advancing genomic mapping approaches, allowing both better manipulation of DNA on the nanoscale and enhanced optical resolving power for analysis of genomic information. During the past few years, these developments have been adopted also for epigenetic studies. The common principle for these studies is the use of advanced optical microscopy for the detection of fluorescently labeled epigenetic marks on long, extended DNA molecules. Here we will discuss recent single-molecule studies for the mapping of chromatin composition and epigenetic DNA modifications, such as DNA methylation.

Beyond sequencing: optical mapping of DNA in the age of nanotechnology and nanoscopy.

Next generation sequencing (NGS) is revolutionizing all fields of biological research but it fails to extract the full range of information associated with genetic material. Optical mapping of DNA grants access to genetic and epigenetic information on individual DNA molecules up to ∼1 Mbp in length. Fluorescent labeling of specific sequence motifs, epigenetic marks and other genomic information on individual DNA molecules generates a high content optical barcode along the DNA. By stretching the DNA to a linear configuration this barcode may be directly visualized by fluorescence microscopy. We discuss the advances of these methods in light of recent developments in nano-fabrication and super-resolution optical imaging (nanoscopy) and review the latest achievements of optical mapping in the context of genomic analysis.

The effect of chemical chaperones on the assembly and stability of HIV-1 capsid protein.

Chemical chaperones are small organic molecules which accumulate in a broad range of organisms in various tissues under different stress conditions and assist in the maintenance of a correct proteostasis under denaturating environments. The effect of chemical chaperones on protein folding and aggregation has been extensively studied and is generally considered to be mediated through non-specific interactions. However, the precise mechanism of action remains elusive. Protein self-assembly is a key event in both native and pathological states, ranging from microtubules and actin filaments formation to toxic amyloids appearance in degenerative disorders, such as Alzheimer's and Parkinson's diseases. Another pathological event, in which protein assembly cascade is a fundamental process, is the formation of virus particles. In the late stage of the virus life cycle, capsid proteins self-assemble into highly-ordered cores, which encapsulate the viral genome, consequently protect genome integrity and mediate infectivity. In this study, we examined the effect of different groups of chemical chaperones on viral capsid assembly in vitro, focusing on HIV-1 capsid protein as a system model. We found that while polyols and sugars markedly inhibited capsid assembly, methylamines dramatically enhanced the assembly rate. Moreover, chemical chaperones that inhibited capsid core formation, also stabilized capsid structure under thermal denaturation. Correspondingly, trimethylamine N-oxide, which facilitated formation of high-order assemblies, clearly destabilized capsid structure under similar conditions. In contrast to the prevailing hypothesis suggesting that chemical chaperones affect proteins through preferential exclusion, the observed dual effects imply that different chaperones modify capsid assembly and stability through different mechanisms. Furthermore, our results indicate a correlation between the folding state of capsid to its tendency to assemble into highly-ordered structures.

Complete phenotypic recovery of an Alzheimer's disease model by a quinone-tryptophan hybrid aggregation inhibitor.

The rational design of amyloid oligomer inhibitors is yet an unmet drug development need. Previous studies have identified the role of tryptophan in amyloid recognition, association and inhibition. Furthermore, tryptophan was ranked as the residue with highest amyloidogenic propensity. Other studies have demonstrated that quinones, specifically anthraquinones, can serve as aggregation inhibitors probably due to the dipole interaction of the quinonic ring with aromatic recognition sites within the amyloidogenic proteins. Here, using in vitro, in vivo and in silico tools we describe the synthesis and functional characterization of a rationally designed inhibitor of the Alzheimer's disease-associated beta-amyloid. This compound, 1,4-naphthoquinon-2-yl-L-tryptophan (NQTrp), combines the recognition capacities of both quinone and tryptophan moieties and completely inhibited Abeta oligomerization and fibrillization, as well as the cytotoxic effect of Abeta oligomers towards cultured neuronal cell line. Furthermore, when fed to transgenic Alzheimer's disease Drosophila model it prolonged their life span and completely abolished their defective locomotion. Analysis of the brains of these flies showed a significant reduction in oligomeric species of Abeta while immuno-staining of the 3(rd) instar larval brains showed a significant reduction in Abeta accumulation. Computational studies, as well as NMR and CD spectroscopy provide mechanistic insight into the activity of the compound which is most likely mediated by clamping of the aromatic recognition interface in the central segment of Abeta. Our results demonstrate that interfering with the aromatic core of amyloidogenic peptides is a promising approach for inhibiting various pathogenic species associated with amyloidogenic diseases. The compound NQTrp can serve as a lead for developing a new class of disease modifying drugs for Alzheimer's disease.