Forest top canopy bacterial communities are influenced by elevation and host tree traits

Background: The phyllosphere microbiome is crucial for plant health and ecosystem functioning. While host species play a determining role in shaping the phyllosphere microbiome, host trees of the same species that are subjected to different environmental conditions can still exhibit large degrees of variation in their microbiome diversity and composition. Whether these intra-specific variations in phyllosphere microbiome diversity and composition can be observed over the broader expanse of forest landscapes remains unclear. In this study, we aim to assess the variation in the top canopy phyllosphere bacterial communities between and within host tree species in the temperate European forests, focusing on Fagus sylvatica (European beech) and Picea abies (Norway spruce).Results: We profiled the bacterial diversity, composition, driving factors, and discriminant taxa in the top canopy phyllosphere of 211 trees in two temperate forests, Veluwe National Parks, the Netherlands and Bavarian Forest National Park, Germany. We found the bacterial communities were primarily shaped by host species, and large variation existed within beech and spruce. While we showed that there was a core microbiome in all tree species examined, community composition varied with elevation, tree diameter at breast height, and leaf-specific traits (e.g., chlorophyll and P content). These driving factors of bacterial community composition also correlated with the relative abundance of specific bacterial families.Conclusions: While our results underscored the importance of host species, we demonstrated a substantial range of variation in phyllosphere bacterial diversity and composition within a host species. Drivers of these variations have implications at both the individual host tree level, where the bacterial communities differed based on tree traits, and at the broader forest landscape level, where drivers like certain highly plastic leaf traits can potentially link forest canopy bacterial community variations to forest ecosystem processes. We eventually showed close associations between forest canopy phyllosphere bacterial communities and host trees exist, and the consistent patterns emerging from these associations are critical for host plant functioning

A high-quality human reference panel reveals the complexity and distribution of genomic structural variants

Structural variation (SV) represents a major source of differences between individual human genomes and has been linked to disease phenotypes. However, the majority of studies provide neither a global view of the full spectrum of these variants nor integrate them into reference panels of genetic variation. Here, we analyse whole genome sequencing data of 769 individuals from 250 Dutch families, and provide a haplotype-resolved map of 1.9 million genome variants across 9 different variant classes, including novel forms of complex indels, and retrotransposition-mediated insertions of mobile elements and processed RNAs. A large proportion are previously under reported variants sized between 21 and 100 bp. We detect 4 megabases of novel sequence, encoding 11 new transcripts. Finally, we show 191 known, trait-associated SNPs to be in strong linkage disequilibrium with SVs and demonstrate that our panel facilitates accurate imputation of SVs in unrelated individuals

WGS-based telomere length analysis in Dutch family trios implicates stronger maternal inheritance and a role for RRM1 gene

Telomere length (TL) regulation is an important factor in ageing, reproduction and cancer development. Genetic, hereditary and environmental factors regulating TL are currently widely investigated, however, their relative contribution to TL variability is still understudied. We have used whole genome sequencing data of 250 family trios from the Genome of the Netherlands project to perform computational measurement of TL and a series of regression and genome-wide association analyses to reveal TL inheritance patterns and associated genetic factors. Our results confirm that TL is a largely heritable trait, primarily with mother’s, and, to a lesser extent, with father’s TL having the strongest influence on the offspring. In this cohort, mother’s, but not father’s age at conception was positively linked to offspring TL. Age-related TL attrition of 40 bp/year had relatively small influence on TL variability. Finally, we have identified TL-associated variations in ribonuclease reductase catalytic subunit M1 (RRM1 gene), which is known to regulate telomere maintenance in yeast. We also highlight the importance of multivariate approach and the limitations of existing tools for the analysis of TL as a polygenic heritable quantitative trait

A framework for the detection of de novo mutations in family-based sequencing data

Germline mutation detection from human DNA sequence data is challenging due to the rarity of such events relative to the intrinsic error rates of sequencing technologies and the uneven coverage across the genome. We developed PhaseByTransmission (PBT) to identify de novo single nucleotide variants and short insertions and deletions (indels) from sequence data collected in parent-offspring trios. We compute the joint probability of the data given the genotype likelihoods in the individual family members, the known familial relationships and a prior probability for the mutation rate. Candidate de novo mutations (DNMs) are reported along with their posterior probability, providing a systematic way to prioritize them for validation. Our tool is integrated in the Genome Analysis Toolkit and can be used together with the ReadBackedPhasing module to infer the parental origin of DNMs based on phase-informative reads. Using simulated data, we show that PBT outperforms existing tools, especially in low coverage data and on the X chromosome. We further show that PBT displays high validation rates on empirical parent-offspring sequencing data for whole-exome data from 104 trios and X-chromosome data from 249 parent-offspring families. Finally, we demonstrate an association between father's age at conception and the number of DNMs in female offspring's X chromosome, consistent with previous literature reports

Skewed X-inactivation is common in the general female population

X-inactivation is a well-established dosage compensation mechanism ensuring that X-chromosomal genes are expressed at comparable levels in males and females. Skewed X-inactivation is often explained by negative selection of one of the alleles. We demonstrate that imbalanced expression of the paternal and maternal X-chromosomes is common in the general population and that the random nature of the X-inactivation mechanism can be sufficient to explain the imbalance. To this end, we analyzed blood-derived RNA and whole-genome sequencing data from 79 female children and their parents from the Genome of the Netherlands project. We calculated the median ratio of the paternal over total counts at all X-chromosomal heterozygous single-nucleotide variants with coverage ≥10. We identified two individuals where the same X-chromosome was inactivated in all cells. Imbalanced expression of the two X-chromosomes (ratios ≤0.35 or ≥0.65) was observed in nearly 50% of the population. The empirically observed skewing is explained by a theoretical model where X-inactivation takes place in an embryonic stage in which eight cells give rise to the hematopoietic compartment. Genes escaping X-inactivation are expressed from both alleles and therefore demonstrate less skewing than inactivated genes. Using this characteristic, we identified three novel escapee genes (SSR4, REPS2, and SEPT6), but did not find support for many previously reported escapee genes in blood. Our collective data suggest that skewed X-inactivation is common in the general population. This may contribute to manifestation of symptoms in carriers of recessive X-linked disorders. We recommend that X-inactivation results should not be used lightly in the interpretation of X-linked variants

OTU representative sequences

Representative sequences for all OTUs for the forward (Euk20F) and Reverse (Euk302r) sequencing direction

Genetic differences among Cedrela odorata sites in Bolivia provide limited potential for fine-scale timber tracing

Illegal trade of tropical timber leads to biodiversity and economic losses worldwide. There is a need for forensic tools that allow tracing the origin of timber and verifying compliance with international and national regulations. We evaluated the potential for genetic tracing of Cedrela odorata, one of the most traded neotropical timbers, within Bolivia. Using a set of seven microsatellites (SSRs), we studied the spatial distribution and genetic diversity and tested whether populations show sufficient genetic discrimination for timber tracing at a national level. Cambium and leaves were sampled from 81 C. odorata trees from three sites, at 268–501-km distance. To explore genetic differentiation, Bayesian clustering and principal component analysis (PCA) were employed. To infer the origin of samples, we conducted kernel discriminant analysis (KDA) based on a PCA that included all alleles and a manual assessment of site-unique alleles. The PCA showed three distinct genetic clusters, but only one of them corresponded with one of the sampled sites. The KDA based on allele frequency had a 33.7% mean classification error, with a considerably lower error (8.2%) for the site which matched with one genetic cluster. The blind test on unique alleles led to a similar classification error (30%). The occurrence of multiple genetic clusters within sites suggests that Bolivian C. odorata populations contain several parental lines, resulting in limited potential for forensic tracing at a national level. Based on our findings, we recommend for additional sampling across the spatial range of C. odorata within the country to support the development of forensic techniques for this species.</p

Binary dataset of alleles from cambium and leaves of Cedrela odorata in Bolivia

Binary dataset of alleles from cambium and leaves of Cedrela odorata in Bolivi

Bioinformatic scripts

Full scripts for bioinformatic pipeline as applied for sequence data analysis for Geisen et al. 2015 Mol. Ecolog