Population transcriptomic characterization of the genetic and expression variation of a candidate progenitor of Miscanthus energy crops

Population transcriptogenomics highlights impaired metabolism and small population sizes in tree frogs living in the Chernobyl Exclusion Zone

BACKGROUND

Individual functional modifications shape the ability of wildlife populations to cope with anthropogenic environmental changes. But instead of adaptive response, human-altered environments can generate a succession of deleterious functional changes leading to the extinction of the population. To study how persistent anthropogenic changes impacted local species' population status, we characterised population structure, genetic diversity and individual response of gene expression in the tree frog Hyla orientalis along a gradient of radioactive contamination around the Chernobyl nuclear power plant.

RESULTS

We detected lower effective population size in populations most exposed to ionizing radiation in the Chernobyl Exclusion Zone that is not compensated by migrations from surrounding areas. We also highlighted a decreased body condition of frogs living in the most contaminated area, a distinctive transcriptomics signature and stop-gained mutations in genes involved in energy metabolism. While the association with dose will remain correlational until further experiments, a body of evidence suggests the direct or indirect involvement of radiation exposure in these changes.

CONCLUSIONS

Despite ongoing migration and lower total dose rates absorbed than at the time of the accident, our results demonstrate that Hyla orientalis specimens living in the Chernobyl Exclusion Zone are still undergoing deleterious changes, emphasizing the long-term impacts of the nuclear disaster.

Mitochondrial Heteroplasmy and PCR Amplification Bias Lead to Wrong Species Delimitation with High Confidence in the South American and Antarctic Marine Bivalve Aequiyoldia eightsii Species Complex

The species delimitation of the marine bivalve species complex Aequiyoldia eightsii in South America and Antarctica is complicated by mitochondrial heteroplasmy and amplification bias in molecular barcoding. In this study, we compare different data sources (mitochondrial cytochrome c oxidase subunit I (COI) sequences; nuclear and mitochondrial SNPs). Whilst all the data suggest that populations on either side of the Drake Passage belong to different species, the picture is less clear within Antarctic populations, which harbor three distinct mitochondrial lineages (p-dist ≈ 6%) that coexist in populations and in a subset of individuals with heteroplasmy. Standard barcoding procedures lead to amplification bias favoring either haplotype unpredictably and thus overestimate the species richness with high confidence. However, nuclear SNPs show no differentiation akin to the trans-Drake comparison, suggesting that the Antarctic populations represent a single species. Their distinct haplotypes likely evolved during periods of temporary allopatry, whereas recombination eroded similar differentiation patterns in the nuclear genome after secondary contact. Our study highlights the importance of using multiple data sources and careful quality control measures to avoid bias and increase the accuracy of molecular species delimitation. We recommend an active search for mitochondrial heteroplasmy and haplotype-specific primers for amplification in DNA-barcoding studies.

Genetic diversity and gene expression diversity shape the adaptive pattern of the aquatic plant Batrachium bungei along an altitudinal gradient on the Qinghai-Tibet plateau

It is an intriguing issue of evolutionary biology how genetic diversity and gene expression diversity shape the adaptive patterns. Comparative transcriptomic studies of wild populations in extreme environments provide critical insights into the relative contribution of genetic and expressive components. In this study, we analyzed the genetic diversity and gene expression diversity of 20 populations of the aquatic plant Batrachium bungei along elevations ranging from 2690 to 4896 m on the Qinghai-Tibet plateau (QTP). Based on single nucleotide polymorphisms (SNPs) and gene expression data from 100 individuals of B. bungei, we found that variation in genetic sequence was more sensitive to detect weak differentiation than gene expression. Using 292,613 high-quality SNPs, we documented a significant phylogeographical structure, a low within-population genetic diversity, and a high inter-population genetic differentiation in B. bungei populations. Analysis of relationship between geographic distance, genetic distance, and gene expression similarity showed that geographic isolation shaped gene flow patterns but not gene expression patterns. We observed a negative relationship between genetic diversity and gene expression diversity within and among B. bungei populations, and we demonstrated that as environmental conditions worsen with increasing altitude, genetic diversity played an increased role in maintaining the stability of populations, while the corresponding role of gene expression diversity decreased. These results suggested that genetic diversity and gene expression diversity might act as a complementary mechanism contributing to the long-term survival of B. bungei in extreme environments.

Population transcriptomic sequencing reveals allopatric divergence and local adaptation in Pseudotaxus chienii (Taxaceae)

BACKGROUND

Elucidating the effects of geography and selection on genetic variation is critical for understanding the relative importance of adaptation in driving differentiation and identifying the environmental factors underlying its occurrence. Adaptive genetic variation is common in tree species, especially widely distributed long-lived species. Pseudotaxus chienii can occupy diverse habitats with environmental heterogeneity and thus provides an ideal material for investigating the process of population adaptive evolution. Here, we characterize genetic and expression variation patterns and investigate adaptive genetic variation in P. chienii populations.

RESULTS

We generated population transcriptome data and identified 13,545 single nucleotide polymorphisms (SNPs) in 5037 unigenes across 108 individuals from 10 populations. We observed lower nucleotide diversity (π = 0.000701) among the 10 populations than observed in other gymnosperms. Significant negative correlations between expression diversity and nucleotide diversity in eight populations suggest that when the species adapts to the surrounding environment, gene expression and nucleotide diversity have a reciprocal relationship. Genetic structure analyses indicated that each distribution region contains a distinct genetic group, with high genetic differentiation among them due to geographical isolation and local adaptation. We used FST outlier, redundancy analysis, and latent factor mixed model methods to detect molecular signatures of local adaptation. We identified 244 associations between 164 outlier SNPs and 17 environmental variables. The mean temperature of the coldest quarter, soil Fe and Cu contents, precipitation of the driest month, and altitude were identified as the most important determinants of adaptive genetic variation. Most candidate unigenes with outlier signatures were related to abiotic and biotic stress responses, and the monoterpenoid biosynthesis and ubiquitin-mediated proteolysis KEGG pathways were significantly enriched in certain populations and deserve further attention in other long-lived trees.

CONCLUSIONS

Despite the strong population structure in P. chienii, genomic data revealed signatures of divergent selection associated with environmental variables. Our research provides SNPs, candidate unigenes, and biological pathways related to environmental variables to facilitate elucidation of the genetic variation in P. chienii in relation to environmental adaptation. Our study provides a promising tool for population genomic analyses and insights into the molecular basis of local adaptation.

Chromosome-scale assembly and analysis of biomass crop Miscanthus lutarioriparius genome

Miscanthus, a rhizomatous perennial plant, has great potential for bioenergy production for its high biomass and stress tolerance. We report a chromosome-scale assembly of Miscanthus lutarioriparius genome by combining Oxford Nanopore sequencing and Hi-C technologies. The 2.07-Gb assembly covers 96.64% of the genome, with contig N50 of 1.71 Mb. The centromere and telomere sequences are assembled for all 19 chromosomes and chromosome 10, respectively. Allotetraploid origin of the M. lutarioriparius is confirmed using centromeric satellite repeats. The tetraploid genome structure and several chromosomal rearrangements relative to sorghum are clearly demonstrated. Tandem duplicate genes of M. lutarioriparius are functional enriched not only in terms related to stress response, but cell wall biosynthesis. Gene families related to disease resistance, cell wall biosynthesis and metal ion transport are greatly expanded and evolved. The expansion of these families may be an important genomic basis for the enhancement of remarkable traits of M. lutarioriparius.

Water Use Efficiency and Stress Tolerance of the Potential Energy Crop Miscanthus lutarioriparius Grown on the Loess Plateau of China

As a potential energy crop with high biomass yield, Miscanthus lutarioriparius (M. lutarioriparius), endemic to the Long River Range in central China, needs to be investigated for its acclimation to stressful climatic and soil conditions often found on the marginal land. In this study, traits related to acclimation and yield, including survival rates, plant height (PH), stem diameter (SD), tiller number (TN), water use efficiency (WUE), and photosynthetic rates (A), were examined for 41 M. lutarioriparius populations that transplanted to the arid and cold Loess Plateau of China. The results showed that the average survival rate of M. lutarioriparius populations was only 4.16% over the first winter but the overwinter rate increased to 35.03% after the second winter, suggesting that plants having survived the first winter could have acclaimed to the low temperature. The strikingly high survival rates over the second winter were found to be 95.83% and 80.85%, respectively, for HG18 and HG39 populations. These populations might be especially valuable for the selection of energy crops for such an area. Those individuals surviving for the two consecutive winters showed significantly higher WUE than those measured after the first winter. The high WUE and low stomatal conductance (g_s) observed in survived individuals could have been responsible for their acclimation to this new and harsh environment. A total of 61 individuals with productive growth traits and strong resistance to cold and drought were identified for further energy crop development. This study showed that the variation of M. lutarioriparius held great potential for developing energy crops following continuous field selection.

Genomic signals found using RNA sequencing show signatures of selection and subtle population differentiation in walleye (Sander vitreus) in a large freshwater ecosystem

RNA sequencing is an effective approach for studying aquatic species yielding both physiological and genomic data. However, its population genetic applications are not well-characterized. We investigate this possible role for RNA sequencing for population genomics in Lake Winnipeg, Manitoba, Canada, walleye (Sander vitreus). Lake Winnipeg walleye represent the largest component of the second-largest freshwater fishery in Canada. In the present study, large female walleye were sampled via nonlethal gill biopsy over two years at three spawning sites representing a latitudinal gradient in the lake. Genetic variation from sequenced mRNA was analyzed for neutral and adaptive markers to investigate population structure and possible adaptive variation. We find low population divergence (F ST = 0.0095), possible northward gene flow, and outlier loci that vary latitudinally in transcripts associated with cell membrane proteins and cytoskeletal function. These results indicate that Lake Winnipeg walleye may be effectively managed as a single demographically connected metapopulation with contributing subpopulations and suggest genomic differences possibly underlying observed phenotypic differences. Despite its high cost relative to other genotyping methods, RNA sequencing data can yield physiological in addition to genetic information discussed here. We therefore argue that it is useful for addressing diverse molecular questions in the conservation of freshwater species.

Utilization of Tissue Ploidy Level Variation in de Novo Transcriptome Assembly of Pinus sylvestris

Compared to angiosperms, gymnosperms lag behind in the availability of assembled and annotated genomes. Most genomic analyses in gymnosperms, especially conifer tree species, rely on the use of de novo assembled transcriptomes. However, the level of allelic redundancy and transcript fragmentation in these assembled transcriptomes, and their effect on downstream applications have not been fully investigated. Here, we assessed three assembly strategies for short-reads data, including the utility of haploid megagametophyte tissue during de novo assembly as single-allele guides, for six individuals and five different tissues in Pinus sylvestris. We then contrasted haploid and diploid tissue genotype calls obtained from the assembled transcriptomes to evaluate the extent of paralog mapping. The use of the haploid tissue during assembly increased its completeness without reducing the number of assembled transcripts. Our results suggest that current strategies that rely on available genomic resources as guidance to minimize allelic redundancy are less effective than the application of strategies that cluster redundant assembled transcripts. The strategy yielding the lowest levels of allelic redundancy among the assembled transcriptomes assessed here was the generation of SuperTranscripts with Lace followed by CD-HIT clustering. However, we still observed some levels of heterozygosity (multiple gene fragments per transcript reflecting allelic redundancy) in this assembled transcriptome on the haploid tissue, indicating that further filtering is required before using these assemblies for downstream applications. We discuss the influence of allelic redundancy when these reference transcriptomes are used to select regions for probe design of exome capture baits and for estimation of population genetic diversity.

Evolutionary and ecological functional genomics, from lab to the wild

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.

Coexpression network revealing the plasticity and robustness of population transcriptome during the initial stage of domesticating energy crop Miscanthus lutarioriparius

Coexpression network revealing genes with Co-variation Expression pattern (CE) and those with Top rank of Expression fold change (TE) played different roles in responding to new environment of Miscanthus lutarioriparius. Variation in gene expression level, the product of genetic and/or environmental perturbation, determines the robustness-to-plasticity spectrum of a phenotype in plants. Understanding how expression variation of plant population response to a new field is crucial to domesticate energy crops. Weighted Gene Coexpression Network Analysis (WGCNA) was used to explore the patterns of expression variation based on 72 Miscanthus lutarioriparius transcriptomes from two contrasting environments, one near the native habitat and the other in one harsh domesticating region. The 932 genes with Co-variation Expression pattern (CE) and other 932 genes with Top rank of Expression fold change (TE) were identified and the former were strongly associated with the water use efficiency (r ≥ 0.55, P ≤ 10^-7). Functional enrichment of CE genes were related to three organelles, which well matched the annotation of twelve motifs identified from their conserved noncoding sequence; while TE genes were mostly related to biotic and/or abiotic stress. The expression robustness of CE genes with high genetic diversity kept relatively stable between environments while the harsh environment reduced the expression robustness of TE genes with low genetic diversity. The expression plasticity of CE genes was increased less than that of TE genes. These results suggested that expression variation of CE genes and TE genes could account for the robustness and plasticity of acclimation ability of Miscanthus, respectively. The patterns of expression variation revealed by transcriptomic network would shed new light on breeding and domestication of energy crops.