Transposases are mobile genetic elements suggested to have an important role in bacterial genome plasticity and host adaptation but their transcriptional activity in natural bacterial communities is largely unexplored. Here we analyzed metagenomes and -transcriptomes of size fractionated (0.1-0.8, 0.8-3.0 and 3.0-200 μm) bacterial communities from the brackish Baltic Sea, and adjacent marine waters. The Baltic Sea transposase levels, up to 1.7% of bacterial genes and 2% of bacterial transcripts, were considerably higher than in marine waters and similar to levels reported for extreme environments. Large variations in expression were found between transposase families and groups of bacteria, with a two-fold higher transcription in Cyanobacteria than in any other phylum. The community-level results were corroborated at the genus level by Synechococcus transposases reaching up to 5.2% of genes and 6.9% of transcripts, which is in contrast to marine Synechococcus that largely lack these genes. Levels peaked in Synechococcus from the largest size fraction, suggesting high frequencies of lateral gene transfer and high genome plasticity in colony-forming picocyanobacteria. Together, the results support an elevated rate of transposition-based genome change and adaptation in bacterial populations of the Baltic Sea, and possibly also of other highly dynamic estuarine waters.

In 2011, a nationwide outbreak of tularaemia occurred in Norway with 180 recorded cases. It was associated with the largest peak in lemming density seen in 40 years. Francisella tularensis was isolated from 18 patients. To study the geographical distribution of F.tularensis genotypes in Norway and correlate genotype with epidemiology and clinical presentation,we performed whole genome sequencing of patient isolates. All 18 genomes from the outbreak carried genetic signatures of F. tularensis subsp. holarctica and were assigned to genetic clades using canonical single nucleotide polymorphisms. Ten isolates were assigned to major genetic clade B.6 (subclade B.7),seven to clade B.12, and one to clade B.4. The B.6 subclade B.7 was most common in southern and central Norway, while clade B.12 was evenly distributed between the southern, central and northern parts of the country. There was no association between genotype and clinical presentation of tularaemia, time of year or specimen type. We found extensive sequence similarity with F. tularensis subsp. holarctica genomes from high-endemic tularaemia areas in Sweden.Finding nearly identical genomes across large geographical distances in Norway and Sweden imply a life cycle of the bacterium without replication between the outbreaks and raise new questions about long-range migration mechanisms.

Detailed linkage and recombination rate maps are necessary to use the full potential of genome sequencing and population genomic analyses. We used a custom collared flycatcher 50 K SNP array to develop a high-density linkage map with 37 262 markers assigned to 34 linkage groups in 33 autosomes and the Z chromosome. The best-order map contained 4215 markers, with a total distance of 3132 cM and a mean genetic distance between markers of 0.12 cM. Facilitated by the array being designed to include markers from most scaffolds, we obtained a second-generation assembly of the flycatcher genome that approaches full chromosome sequences (N50 super-scaffold size 20.2 Mb and with 1.042 Gb (of 1.116 Gb) anchored to and mostly ordered and oriented along chromosomes). We found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5-2.0 event (6.6-7.5 Mb) per My have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution. The mean recombination rate was 3.1 cM/Mb and correlated closely with chromosome size, from 2 cM/Mb for chromosomes >100 Mb to >10 cM/Mb for chromosomes <10 Mb. This size dependence seemed entirely due to an obligate recombination event per chromosome; if 50 cM was subtracted from the genetic lengths of chromosomes, the rate per physical unit DNA was constant across chromosomes. Flycatcher recombination rate showed similar variation along chromosomes as chicken but lacked the large interior recombination deserts characteristic of zebra finch chromosomes.

Variants in the growth factor receptor-bound protein 10 (GRB10) gene were in a GWAS meta-analysis associated with reduced glucose-stimulated insulin secretion and increased risk of type 2 diabetes (T2D) if inherited from the father, but inexplicably reduced fasting glucose when inherited from the mother. GRB10 is a negative regulator of insulin signaling and imprinted in a parent-of-origin fashion in different tissues. GRB10 knock-down in human pancreatic islets showed reduced insulin and glucagon secretion, which together with changes in insulin sensitivity may explain the paradoxical reduction of glucose despite a decrease in insulin secretion. Together, these findings suggest that tissue-specific methylation and possibly imprinting of GRB10 can influence glucose metabolism and contribute to T2D pathogenesis. The data also emphasize the need in genetic studies to consider whether risk alleles are inherited from the mother or the father.

Aging is a major risk factor for many chronic diseases. This study aimed to examine the effects of antihypertensive, lipid-lowering, and antidiabetic drugs on biological aging. We included 672 participants and 2746 repeated measurements from the Swedish Adoption/Twin Study of Aging. Self-reported medicine uses were categorized into antidiabetic, antihypertensive, and lipid-lowering drugs. A total of 12 biomarkers for biological aging (BA biomarkers) were included as outcomes. Conditional generalized estimating equations were applied conditioning on individuals to estimate the drug effect on BA biomarker level within the same person when using or not using the drug. Chronological age, body mass index, smoking status, number of multiple medication uses, blood pressure, blood glucose level, and apoB/apoA ratio were adjusted for as covariates in the model. Overall, using antihypertensive drugs was associated with a decrease in one DNA-methylation age (PCGrimAge: beta = - 0.39, 95%CI = - 0.67 to - 0.12). When looking into drug subcategories, calcium channel blockers (CCBs) were associated with a decrease in several DNA-methylation ages (PCHorvathAge beta = - 1.28, 95%CI = - 2.34 to - 0.21; PCSkin&bloodAge beta = - 1.34, 95%CI = - 2.61 to - 0.07; PCPhenoAge beta = - 1.74, 95%CI = - 2.58 to - 0.89; PCGrimAge beta = - 0.57, 95%CI = - 0.96 to - 0.17) and in functional biological ages (functional age index beta = - 2.18, 95%CI = - 3.65 to - 0.71; frailty index beta = - 1.31, 95%CI = - 2.43 to - 0.18). However, the results within other drug subcategories were inconsistent. Calcium channel blockers may decrease biological aging captured by the BA biomarkers measured at epigenetic and functional level. Future studies are warranted to confirm these effects and understand the underlying biological mechanisms.

The genetic complexity of polygenic traits represents a captivating and intricate facet of biological inheritance. Unlike Mendelian traits controlled by a single gene, polygenic traits are influenced by multiple genetic loci, each exerting a modest effect on the trait. This cumulative impact of numerous genes, interactions among them, environmental factors, and epigenetic modifications results in a multifaceted architecture of genetic contributions to complex traits. Given the well-characterized genome, diverse traits, and range of genetic resources, chicken (Gallus gallus) was employed as a model organism to dissect the intricate genetic makeup of a previously identified major Quantitative Trait Loci (QTL) for body weight on chromosome 1. A multigenerational advanced intercross line (AIL) of 3215 chickens whose genomes had been sequenced to an average of 0.4x was analyzed using genome-wide association study (GWAS) and variance-heterogeneity GWAS (vGWAS) to identify markers associated with 8-week body weight. Additionally, epistatic interactions were studied using the natural and orthogonal interaction (NOIA) model. Six genetic modules, two from GWAS and four from vGWAS, were strongly associated with the studied trait. We found evidence of both additive- and non-additive interactions between these modules and constructed a putative local epistasis network for the region. Our screens for functional alleles revealed a missense variant in the gene ribonuclease H2 subunit B (RNASEH2B), which has previously been associated with growth-related traits in chickens and Darwin's finches. In addition, one of the most strongly associated SNPs identified is located in a non-coding region upstream of the long non-coding RNA, ENSGALG00000053256, previously suggested as a candidate gene for regulating chicken body weight. By studying large numbers of individuals from a family material using approaches to capture both additive and non-additive effects, this study advances our understanding of genetic complexities in a highly polygenic trait and has practical implications for poultry breeding and agriculture.

Novel sequences (NSs), not present in the human reference genome, are abundant and remain largely unexplored. Here, we utilize de novo assembly to study NS in 1,000 Swedish individuals first sequenced as part of the SweGen project revealing a total of 46 Mb in 61,044 distinct contigs of sequences not present in GRCh38. The contigs were aligned to recently published catalogs of Icelandic and Pan-African NSs, as well as the chimpanzee genome, revealing a great diversity of shared sequences. Analyzing the positioning of NS across the chimpanzee genome, we find that 2,807 NS align confidently within 143 chimpanzee orthologs of human genes. Aligning the whole genome sequencing data to the chimpanzee genome, we discover ancestral NS common throughout the Swedish population. The NSs were searched for repeats and repeat elements: revealing a majority of repetitive sequence (56%), and enrichment of simple repeats (28%) and satellites (15%). Lastly, we align the unmappable reads of a subset of the thousand genomes data to our collection of NS, as well as the previously published Pan-African NS: revealing that both the Swedish and Pan-African NS are widespread, and that the Swedish NSs are largely a subset of the Pan-African NS. Overall, these results highlight the importance of creating a more diverse reference genome and illustrate that significant amounts of the NS may be of ancestral origin.