Belowground Consequences of Intracontinental Range-Expanding Plants and Related Natives in Novel Environments

Introduced exotic plant species that originate from other continents are known to alter soil microbial community composition and nutrient cycling. Plant species that expand range to higher latitudes and altitudes as a consequence of current climate warming might as well affect the composition and functioning of native soil communities in their new range. However, the functional consequences of plant origin have been poorly studied in the case of plant range shifts. Here, we determined rhizosphere bacterial communities of four intracontinental range-expanding plant species in comparison with their four congeneric natives grown in soils collected from underneath those plant species in the field and in soils that are novel to them. We show that, when controlling for both species relatedness and soil characteristics, range-expanding plant species in higher latitude ecosystems will influence soil bacterial community composition and nutrient cycling in a manner similar to congeneric related native species. Our results highlight the importance to include phylogenetically controlled comparisons to disentangle the effect of origin from the effect of contrasting plant traits in the context of exotic plant species

Aging Uncouples Heritability and Expression-QTL in Caenorhabditis elegans

The number and distribution of gene expression QTL (eQTL) represent the genetic architecture of many complex traits, including common human diseases. We previously reported that the heritable eQTL patterns are highly dynamic with age in an N2 × CB4856 recombinant inbred population of the nematode Caenorhabditis elegans. In particular, we showed that the number of eQTL decreased with age. Here, we investigated the reason for this decrease by combining gene expression profiles at three ages in the wild types N2 and CB4856 with the reported expression profiles of the RIL population. We determined heritability and transgression (when gene expression levels in the RILs are more extreme than the parents) and investigated their relation with eQTL changes with age. Transgressive segregation was widespread but depended on physiological age. The percentage of genes with an eQTL increased with a higher heritability in young worms. However, for old worms this percentage hardly increased. Using a single marker approach, we found that almost 20% of genes with heritability >0.9 had an eQTL in developing worms. Surprisingly, only 10% was found in old worms. Using a multimarker approach, this percentage increased to almost 30% for both age groups. Comparison of the single marker to a multiple marker eQTL mapping indicated that heritable regulation of gene expression becomes more polygenic in aging worms due to multiple loci and possible epistatic interactions. We conclude that linkage studies should account for the relation between increased polygenic regulation and diminished effects at older ages

Genetical Genomics Reveals Large Scale Genotype-By-Environment Interactions in Arabidopsis thaliana

One of the major goals of quantitative genetics is to unravel the complex interactions between molecular genetic factors and the environment. The effects of these genotype-by-environment interactions also affect and cause variation in gene expression. The regulatory loci responsible for this variation can be found by genetical genomics that involves the mapping of quantitative trait loci (QTLs) for gene expression traits also called expression-QTL (eQTLs). Most genetical genomics experiments published so far, are performed in a single environment and hence do not allow investigation of the role of genotype-by-environment interactions. Furthermore, most studies have been done in a steady state environment leading to acclimated expression patterns. However a response to the environment or change therein can be highly plastic and possibly lead to more and larger differences between genotypes. Here we present a genetical genomics study on 120 Arabidopsis thaliana, Landsberg erecta × Cape Verde Islands, recombinant inbred lines (RILs) in active response to the environment by treating them with 3 h of shade. The results of this experiment are compared to a previous study on seedlings of the same RILs from a steady state environment. The combination of two highly different conditions but exactly the same RILs with a fixed genetic variation showed the large role of genotype-by-environment interactions on gene expression levels. We found environment-dependent hotspots of transcript regulation. The major hotspot was confirmed by the expression profile of a near isogenic line. Our combined analysis leads us to propose CSN5A, a COP9 signalosome component, as a candidate regulator for the gene expression response to shade

Towards increased shading potential: a combined phenotypic and genetic analysis of rice shoot architecture

Rice feeds more than half of the world’s human population. In modern rice farming, a major constraint for productivity is weed proliferation and the ecological impact of herbicide application. Increased weed competitiveness of commercial rice varieties requires enhanced shade casting to limit growth of shade-sensitive weeds and the need for herbicide. We aimed to identify traits that enhance rice shading capacity based on the canopy architecture and the underlying genetic components. We performed a phenotypic screen of a rice diversity panel comprised of 344 varieties, examining 13 canopy architecture traits linked with shading capacity in 4-week-old plants. The analysis revealed a vast range of phenotypic variation across the diversity panel. We used trait correlation and clustering to identify core traits that define shading capacity to be shoot area, number of leaves, culm and solidity (the compactness of the shoot). To simplify the complex canopy architecture, these traits were combined into a Shading Rank metric that is indicative of a plant’s ability to cast shade. Genome wide association study (GWAS) revealed genetic loci underlying canopy architecture traits, out of which five loci were substantially contributing to shading potential. Subsequent haplotype analysis further explored allelic variation and identified seven haplotypes associated with increased shading. Identification of traits contributing to shading capacity and underlying allelic variation presented in this study will serve future genomic assisted breeding programmes. The investigated diversity panel, including widely grown varieties, shows that there is big potential and genetic resources for improvement of elite breeding lines. Implementing increased shading in rice breeding will make its farming less dependent on herbicides and contribute towards more environmentally sustainable agriculture. One sentence summary Through screening a rice diversity panel for variation in shoot architecture, we identified traits corresponding to plant shading potential and their genetic constituents

Contribution of trans regulatory eQTL to cryptic genetic variation in C-elegans

BackgroundCryptic genetic variation (CGV) is the hidden genetic variation that can be unlocked by perturbing normal conditions. CGV can drive the emergence of novel complex phenotypes through changes in gene expression. Although our theoretical understanding of CGV has thoroughly increased over the past decade, insight into polymorphic gene expression regulation underlying CGV is scarce. Here we investigated the transcriptional architecture of CGV in response to rapid temperature changes in the nematode Caenorhabditis elegans. We analyzed regulatory variation in gene expression (and mapped eQTL) across the course of a heat stress and recovery response in a recombinant inbred population.ResultsWe measured gene expression over three temperature treatments: i) control, ii) heat stress, and iii) recovery from heat stress. Compared to control, exposure to heat stress affected the transcription of 3305 genes, whereas 942 were affected in recovering animals. These affected genes were mainly involved in metabolism and reproduction. The gene expression pattern in recovering animals resembled both the control and the heat-stress treatment. We mapped eQTL using the genetic variation of the recombinant inbred population and detected 2626 genes with an eQTL in the heat-stress treatment, 1797 in the control, and 1880 in the recovery. The cis-eQTL were highly shared across treatments. A considerable fraction of the trans-eQTL (40–57%) mapped to 19 treatment specific trans-bands. In contrast to cis-eQTL, trans-eQTL were highly environment specific and thus cryptic. Approximately 67% of the trans-eQTL were only induced in a single treatment, with heat-stress showing the most unique trans-eQTL.ConclusionsThese results illustrate the highly dynamic pattern of CGV across three different environmental conditions that can be evoked by a stress response over a relatively short time-span (2 h) and that CGV is mainly determined by response related trans regulatory eQTL.<br/

WormQTL-public archive and analysis web portal for natural variation data in Caenorhabditis spp

Here, we present WormQTL (http://www.wormqtl.org), an easily accessible database enabling search, comparative analysis and meta-analysis of all data on variation in Caenorhabditis spp. Over the past decade, Caenorhabditis elegans has become instrumental for molecular quantitative genetics and the systems biology of natural variation. These efforts have resulted in a valuable amount of phenotypic, high-throughput molecular and genotypic data across different developmental worm stages and environments in hundreds of C. elegans strains. WormQTL provides a workbench of analysis tools for genotype-phenotype linkage and association mapping based on but not limited to R/qtl (http://www.rqtl.org). All data can be uploaded and downloaded using simple delimited text or Excel formats and are accessible via a public web user interface for biologists and R statistic and web service interfaces for bioinformaticians, based on open source MOLGENIS and xQTL workbench software. WormQTL welcomes data submissions from other worm researcher

Field cricket genome reveals the footprint of recent, abrupt adaptation in the wild.

Evolutionary adaptation is generally thought to occur through incremental mutational steps, but large mutational leaps can occur during its early stages. These are challenging to study in nature due to the difficulty of observing new genetic variants as they arise and spread, but characterizing their genomic dynamics is important for understanding factors favoring rapid adaptation. Here, we report genomic consequences of recent, adaptive song loss in a Hawaiian population of field crickets (Teleogryllus oceanicus). A discrete genetic variant, flatwing, appeared and spread approximately 15 years ago. Flatwing erases sound-producing veins on male wings. These silent flatwing males are protected from a lethal, eavesdropping parasitoid fly. We sequenced, assembled and annotated the cricket genome, produced a linkage map, and identified a flatwing quantitative trait locus covering a large region of the X chromosome. Gene expression profiling showed that flatwing is associated with extensive genome-wide effects on embryonic gene expression. We found that flatwing male crickets express feminized chemical pheromones. This male feminizing effect, on a different sexual signaling modality, is genetically associated with the flatwing genotype. Our findings suggest that the early stages of evolutionary adaptation to extreme pressures can be accompanied by greater genomic and phenotypic disruption than previously appreciated, and highlight how abrupt adaptation might involve suites of traits that arise through pleiotropy or genomic hitchhiking

Contrasting Management Techniques and Soil Types Affect Network Connections between Soil Properties and the Tulip Microbiome

Conventional agriculture relies on high inputs of chemical fertilisers and pesticides, leading to extensive soil damage. As a result, there has been a shift towards organic cultivation to ameliorate these negative effects. However, knowledge gaps remain concerning how different agricultural management regimes can alter soil properties, root-associated biota and relationships between them. We sampled soils and bulbs from clay and sand tulip fields with conventional, organic and hybrid (i.e., intermediate between conventional and organic) management regimes to determine effects on abiotic and biotic soil properties, tulip bulb microbiomes/biota and network connectivity between them. We found that soil abiotic properties were not driven by management. However, soil microbial community composition was determined by interactions between soil type and management. Specifically, conventional management created more heterogeneous and homogenous bacterial and fungal soil communities, respectively, across soil types, while organic and hybrid soil communities differed between soil types. A partial least squares path model revealed that across all fields, soil properties determined the composition of tulip bulb biota (i.e., microorganisms, nematodes), while management techniques did not play a significant role, neither directly nor indirectly. Network analyses revealed more modular (i.e., independent) nodes of association in fields with sandy soils. Specifically, management techniques (i.e., chemical usage, ploughing depth, diminished crop rotation) formed a cluster that was independent from most soil (abiotic and biotic) and bulb biotic properties, while clay field networks showed the opposite. Our results indicate that conventional agriculture can homogenise soil microbial communities, with potential impacts on soil function and buffering capacity to stress. Increased modularity between soil properties, management techniques and tulip bulb biota in sand fields could mean increased resistance to disturbance and abiotic and biotic stress as a result of higher functional redundancy. Interactions between agricultural management regime and soil type must be taken into consideration when determining long-term influences on crop microbiome/biota associations and the potential effects on soil health and productivity

A genome-wide library of CB4856/N2 introgression lines of Caenorhabditis elegans

Recombinant inbred lines (RILs) derived from Caenorhabditis elegans wild-type N2 and CB4856 are increasingly being used for mapping genes underlying complex traits. To speed up mapping and gene discovery, introgression lines (ILs) offer a powerful tool for more efficient QTL identification. We constructed a library of 90 ILs, each carrying a single homozygous CB4856 genomic segment introgressed into the genetic background of N2. The ILs were genotyped by 123 single-nucleotide polymorphism (SNP) markers. The proportion of the CB4856 segments in most lines does not exceed 3%, and together the introgressions cover 96% of the CB4856 genome. The value of the IL library was demonstrated by identifying novel loci underlying natural variation in two ageing-related traits, i.e. lifespan and pharyngeal pumping rate. Bin mapping of lifespan resulted in six QTLs, which all have a lifespan-shortening effect on the CB4856 allele. We found five QTLs for the decrease in pumping rate, of which four colocated with QTLs found for average lifespan. This suggests pleiotropic or closely linked QTL associated with lifespan and pumping rate. Overall, the presented IL library provides a versatile resource toward easier and efficient fine mapping and functional analyses of loci and genes underlying complex traits in C. elegans

A metagenomic portrait of the microbial community responsible for two decades of bioremediation of poly-contaminated groundwater

Biodegradation of pollutants is a sustainable and cost-effective solution to groundwater pollution. Here, we investigate microbial populations involved in biodegradation of poly-contaminants in a pipeline for heavily contaminated groundwater. Groundwater moves from a polluted park to a treatment plant, where an aerated bioreactor effectively removes the contaminants. While the biomass does not settle in the reactor, sediment is collected afterwards and used to seed the new polluted groundwater via a backwash cycle. The pipeline has successfully operated since 1999, but the biological components in the reactor and the contaminated park groundwater have never been described. We sampled seven points along the pipeline, representing the entire remediation process, and characterized the changing microbial communities using genome-resolved metagenomic analysis. We assembled 297 medium- and high-quality metagenome-assembled genome sequences representing on average 46.3% of the total DNA per sample. We found that the communities cluster into two distinct groups, separating the anaerobic communities in the park groundwater from the aerobic communities inside the plant. In the park, the community is dominated by members of the genus Sulfuricurvum, while the plant is dominated by generalists from the order Burkholderiales. Known aromatic compound biodegradation pathways are four times more abundant in the plant-side communities compared to the park-side. Our findings provide a genome-resolved portrait of the microbial community in a highly effective groundwater treatment system that has treated groundwater with a complex contamination profile for two decades