Re-analysis of the Levanluhta skeletal material : Sex and stature estimation of individuals in an Iron Age water burial in Finland

Levanluhta, an Iron Age water burial site in Finland, and its material consisting of commingled skeletal remains and artifacts, has been studied by several researchers over the past 100 years, resulting in multiple interpretations of the people and the site. Previous skeletal analyses have concluded that the majority of the individuals represented in the remains were females and children and were of relatively short stature, so possibly nutritionally deprived. This study re-analyzed the commingled adult human remains with updated methods. The methods applied in this study to estimate sex and stature were based on more representative European reference samples than the previously applied methods. The methods included morphology, osteometrics, and computed tomography (CT) scans. Our results indicated that depending on the reference data, the majority of the individual adult bones including os coxae (73%, n = 45) and long bones (humerus 83%-89%, n = 52; radius 72%-89%, n = 47; ulna 50%-65%, n = 58; femur 92%-100%, n = 25; tibia 77%-85%, n = 26) were classified as females based on their size and morphology. The cross-sectional bone properties of humerii, femora, and tibiae visualized using CT scanning also supported these findings. However, the cranial morphology did not show as clear female-biased sex ratio as other methods (42% females, 33% males, 24% undetermined, n = 33). In females, the mean stature based on the tibia (155.3 cm, n = 10) was within the range of the coeval European females and did not necessarily indicate nutritional deprivation, which is in line with previously published stable isotope findings from the site. The mean stature based on the tibia suggested that the Levanluhta males were short (164.0 cm, n = 3), but final interpretations were limited due to the small number of male individuals. The current study affirmed that the Levanluhta skeletal assemblage was female biased and gave new insights into interpretation of the stature.Peer reviewe

Landscape of multi-nucleotide variants in 125,748 human exomes and 15,708 genomes

Multi-nucleotide variants (MNVs), defined as two or more nearby variants existing on the same haplotype in an individual, are a clinically and biologically important class of genetic variation. However, existing tools typically do not accurately classify MNVs, and understanding of their mutational origins remains limited. Here, we systematically survey MNVs in 125,748 whole exomes and 15,708 whole genomes from the Genome Aggregation Database (gnomAD). We identify 1,792,248 MNVs across the genome with constituent variants falling within 2bp distance of one another, including 18,756 variants with a novel combined effect on protein sequence. Finally, we estimate the relative impact of known mutational mechanisms - CpG deamination, replication error by polymerase zeta, and polymerase slippage at repeat junctions - on the generation of MNVs. Our results demonstrate the value of haplotype-aware variant annotation, and refine our understanding of genome-wide mutational mechanisms of MNVs. Multi-nucleotide variants (MNV) are genetic variants in close proximity of each other on the same haplotype whose functional impact is difficult to predict if they reside in the same codon. Here, Wang et al. use the gnomAD dataset to assemble a catalogue of MNVs and estimate their global mutation rate.Peer reviewe

Free energy barrier for melittin reorientation from a membrane-bound state to a transmembrane state

An important step in a phospholipid membrane pore formation by melittin antimicrobial peptide is a reorientation of the peptide from a surface into a transmembrane conformation. In this work we perform umbrella sampling simulations to calculate the potential of mean force (PMF) for the reorientation of melittin from a surface-bound state to a transmembrane state and provide a molecular level insight into understanding peptide and lipid properties that influence the existence of the free energy barrier. The PMFs were calculated for a peptide to lipid (P/L) ratio of 1/128 and 4/128. We observe that the free energy barrier is reduced when the P/L ratio increased. In addition, we study the cooperative effect; specifically we investigate if the barrier is smaller for a second melittin reorientation, given that another neighboring melittin was already in the transmembrane state. We observe that indeed the barrier of the PMF curve is reduced in this case, thus confirming the presence of a cooperative effect

Transcript expression-aware annotation improves rare variant interpretation

The acceleration of DNA sequencing in samples from patients and population studies has resulted in extensive catalogues of human genetic variation, but the interpretation of rare genetic variants remains problematic. A notable example of this challenge is the existence of disruptive variants in dosage-sensitive disease genes, even in apparently healthy individuals. Here, by manual curation of putative loss-of-function (pLoF) variants in haploinsufficient disease genes in the Genome Aggregation Database (gnomAD)(1), we show that one explanation for this paradox involves alternative splicing of mRNA, which allows exons of a gene to be expressed at varying levels across different cell types. Currently, no existing annotation tool systematically incorporates information about exon expression into the interpretation of variants. We develop a transcript-level annotation metric known as the 'proportion expressed across transcripts', which quantifies isoform expression for variants. We calculate this metric using 11,706 tissue samples from the Genotype Tissue Expression (GTEx) project(2) and show that it can differentiate between weakly and highly evolutionarily conserved exons, a proxy for functional importance. We demonstrate that expression-based annotation selectively filters 22.8% of falsely annotated pLoF variants found in haploinsufficient disease genes in gnomAD, while removing less than 4% of high-confidence pathogenic variants in the same genes. Finally, we apply our expression filter to the analysis of de novo variants in patients with autism spectrum disorder and intellectual disability or developmental disorders to show that pLoF variants in weakly expressed regions have similar effect sizes to those of synonymous variants, whereas pLoF variants in highly expressed exons are most strongly enriched among cases. Our annotation is fast, flexible and generalizable, making it possible for any variant file to be annotated with any isoform expression dataset, and will be valuable for the genetic diagnosis of rare diseases, the analysis of rare variant burden in complex disorders, and the curation and prioritization of variants in recall-by-genotype studies.Peer reviewe

Novel subgroups of adult-onset diabetes and their association with outcomes : a data-driven cluster analysis of six variables

Background Diabetes is presently classified into two main forms, type 1 and type 2 diabetes, but type 2 diabetes in particular is highly heterogeneous. A refined classification could provide a powerful tool to individualise treatment regimens and identify individuals with increased risk of complications at diagnosis. Methods We did data-driven cluster analysis (k-means and hierarchical clustering) in patients with newly diagnosed diabetes (n=8980) from the Swedish All New Diabetics in Scania cohort. Clusters were based on six variables (glutamate decarboxylase antibodies, age at diagnosis, BMI, HbA(1c), and homoeostatic model assessment 2 estimates of beta-cell function and insulin resistance), and were related to prospective data from patient records on development of complications and prescription of medication. Replication was done in three independent cohorts: the Scania Diabetes Registry (n=1466), All New Diabetics in Uppsala (n=844), and Diabetes Registry Vaasa (n=3485). Cox regression and logistic regression were used to compare time to medication, time to reaching the treatment goal, and risk of diabetic complications and genetic associations. Findings We identified five replicable clusters of patients with diabetes, which had significantly different patient characteristics and risk of diabetic complications. In particular, individuals in cluster 3 (most resistant to insulin) had significantly higher risk of diabetic kidney disease than individuals in clusters 4 and 5, but had been prescribed similar diabetes treatment. Cluster 2 (insulin deficient) had the highest risk of retinopathy. In support of the clustering, genetic associations in the clusters differed from those seen in traditional type 2 diabetes. Interpretation We stratified patients into five subgroups with differing disease progression and risk of diabetic complications. This new substratification might eventually help to tailor and target early treatment to patients who would benefit most, thereby representing a first step towards precision medicine in diabetes.Peer reviewe

Genetics of migraine in the age of genome-wide association studies

Genetic factors importantly contribute to migraine. However, unlike for rare monogenic forms of migraine, approaches to identify genes for common forms of migraine have been of limited success. Candidate gene association studies were often negative and positive results were often not replicated or replication failed. Further, the significance of positive results from linkage studies remains unclear owing to the inability to pinpoint the genes under the peaks that may be involved in migraine. Problems hampering these studies include limited sample sizes, methods of migraine ascertainment, and the heterogeneous clinical phenotype. Three genome-wide association studies are available now and have successfully identified four new genetic variants associated with migraine. One new variant (rs1835740) modulates glutamate homeostasis, thus integrates well with current concepts of neurotransmitter disturbances. This variant may be more specific for severe forms of migraine such as migraine with aura than migraine without aura. Another variant (rs11172113) implicates the lipoprotein receptor LRP1, which may interact with neuronal glutamate receptors, thus also providing a link to the glutamate pathway. In contrast, rs10166942 is in close proximity to TRPM8, which codes for a cold and pain sensor. For the first time this links a gene explicitly implicated in pain related pathways to migraine. The potential function of the fourth variant rs2651899 (PRDM16) in migraine is unclear. All these variants only confer a small to moderate change in risk for migraine, which concurs with migraine being a heterogeneous disorder. Ongoing large international collaborations will likely identify additional gene variants for migraine