Evolutionary Analyses of Gene Expression Divergence in Panicum hallii: Exploring Constitutive and Plastic Responses Using Reciprocal Transplants.

The evolution of gene expression is thought to be an important mechanism of local adaptation and ecological speciation. Gene expression divergence occurs through the evolution of cis- polymorphisms and through more widespread effects driven by trans-regulatory factors. Here, we explore expression and sequence divergence in a large sample of Panicum hallii accessions encompassing the species range using a reciprocal transplantation experiment. We observed widespread genotype and transplant site drivers of expression divergence, with a limited number of genes exhibiting genotype-by-site interactions. We used a modified FST-QST outlier approach (QPC analysis) to detect local adaptation. We identified 514 genes with constitutive expression divergence above and beyond the levels expected under neutral processes. However, no plastic expression responses met our multiple testing correction as QPC outliers. Constitutive QPC outlier genes were involved in a number of developmental processes and responses to abiotic environments. Leveraging earlier expression quantitative trait loci results, we found a strong enrichment of expression divergence, including for QPC outliers, in genes previously identified with cis and cis-environment interactions but found no patterns related to trans-factors. Population genetic analyses detected elevated sequence divergence of promoters and coding sequence of constitutive expression outliers but little evidence for positive selection on these proteins. Our results are consistent with a hypothesis of cis-regulatory divergence as a primary driver of expression divergence in P. hallii.

Transcriptome and DNA methylome divergence of inflorescence development between 2 ecotypes in Panicum hallii.

The morphological diversity of the inflorescence determines flower and seed production, which is critical for plant adaptation. Hall's panicgrass (Panicum hallii, P. hallii) is a wild perennial grass that has been developed as a model to study perennial grass biology and adaptive evolution. Highly divergent inflorescences have evolved between the 2 major ecotypes in P. hallii, the upland ecotype (P. hallii var hallii, HAL2 genotype) with compact inflorescence and large seed and the lowland ecotype (P. hallii var filipes, FIL2 genotype) with an open inflorescence and small seed. Here we conducted a comparative analysis of the transcriptome and DNA methylome, an epigenetic mark that influences gene expression regulation, across different stages of inflorescence development using genomic references for each ecotype. Global transcriptome analysis of differentially expressed genes (DEGs) and co-expression modules underlying the inflorescence divergence revealed the potential role of cytokinin signaling in heterochronic changes. Comparing DNA methylome profiles revealed a remarkable level of differential DNA methylation associated with the evolution of P. hallii inflorescence. We found that a large proportion of differentially methylated regions (DMRs) were located in the flanking regulatory regions of genes. Intriguingly, we observed a substantial bias of CHH hypermethylation in the promoters of FIL2 genes. The integration of DEGs, DMRs, and Ka/Ks ratio results characterized the evolutionary features of DMR-associated DEGs that contribute to the divergence of the P. hallii inflorescence. This study provides insights into the transcriptome and epigenetic landscape of inflorescence divergence in P. hallii and a genomic resource for perennial grass biology.

Quantitative genetic-by-soil microbiome interactions in a perennial grass affect functional traits.

Plants interact with diverse microbiomes that can impact plant growth and performance. Recent studies highlight the potential beneficial aspects of plant microbiomes, including the possibility that microbes facilitate the process of local adaptation in their host plants. Microbially mediated local adaptation in plants occurs when local host genotypes have higher fitness than foreign genotypes because of their affiliation with locally beneficial microbes. Here, plant adaptation results from genetic interactions of the host with locally beneficial microbes (e.g. host genotype-by-microbiome interactions). We used a recombinant inbred line (RIL) mapping population derived from upland and lowland ecotypes of the diploid C4 perennial bunch grass Panicum hallii to explore quantitative genetic responses to soil microbiomes focusing on functional root and shoot traits involved in ecotypic divergence. We show that the growth and development of ecotypes and their trait divergence depends on soil microbiomes. Moreover, we find that the genetic architecture is modified by soil microbiomes, revealing important plant genotype-by-microbiome interactions for quantitative traits. We detected a number of quantitative trait loci (QTL) that interact with the soil microbiome. Our results highlight the importance of microbial interactions in ecotypic divergence and trait genetic architecture in C4 perennial grasses.

Salt tolerance QTLs of an endemic rice landrace, Horkuch at seedling and reproductive stages.

Salinity has a significant negative impact on production of rice. To cope with the increased soil salinity due to climate change, we need to develop salt tolerant rice varieties that can maintain their high yield. Rice landraces indigenous to coastal Bangladesh can be a great resource to study the genetic basis of salt adaptation. In this study, we implemented a QTL analysis framework with a reciprocal mapping population developed from a salt tolerant landrace Horkuch and a high yielding rice variety IR29. Our aim was to detect genetic loci that contributes to the salt adaptive responses of the two different developmental stages of rice which are very sensitive to salinity stress. We identified 14 QTLs for 9 traits and found that most are unique to specific developmental stages. In addition, we detected a significant effect of the cytoplasmic genome on the QTL model for some traits such as leaf total potassium and filled grain weight. This underscores the importance of considering cytoplasm-nuclear interaction for breeding programs. Finally, we identified QTLs co-localization for multiple traits that highlights the possible constraint of multiple QTL selection for breeding programs due to different contributions of a donor allele for different traits.

Natural variation in growth and leaf ion homeostasis in response to salinity stress in Panicum hallii.

Soil salinity can negatively impact plants growth, development and fitness. Natural plant populations restricted to coastal environments may evolve in response to saline habitats and therefore provide insights into the process of salinity adaptation. We investigated the growth and physiological responses of coastal and inland populations of Panicum hallii to experimental salinity treatments. Coastal genotypes demonstrated less growth reduction and superior ion homeostasis compared to the inland genotypes in response to saline conditions, supporting a hypothesis of local adaptation. We identified several QTL associated with the plasticity of belowground biomass, leaf sodium and potassium content, and their ratio which underscores the genetic variation present in this species for salinity responses. Genome-wide transcriptome analysis in leaf and root tissue revealed tissue specific overexpression of genes including several cation transporters in the coastal genotype. These transporters mediate sodium ion compartmentalization and potassium ion retention and thus suggests that maintenance of ionic homeostasis of the coastal genotypes might be due to the regulation of these ion transporters. These findings contribute to our understanding of the genetics and molecular mechanisms of salinity adaptation in natural populations, and widens the scope for genetic manipulation of these candidate genes to design plants more resilient to climate change.

A Pleiotropic Flowering Time QTL Exhibits Gene-by-Environment Interaction for Fitness in a Perennial Grass.

Appropriate flowering time is a crucial adaptation impacting fitness in natural plant populations. Although the genetic basis of flowering variation has been extensively studied, its mechanisms in nonmodel organisms and its adaptive value in the field are still poorly understood. Here, we report new insights into the genetic basis of flowering time and its effect on fitness in Panicum hallii, a native perennial grass. Genetic mapping in populations derived from inland and coastal ecotypes identified flowering time quantitative trait loci (QTL) and many exhibited extensive QTL-by-environment interactions. Patterns of segregation within recombinant hybrids provide strong support for directional selection driving ecotypic divergence in flowering time. A major QTL on chromosome 5 (q-FT5) was detected in all experiments. Fine-mapping and expression studies identified a gene with orthology to a rice FLOWERING LOCUS T-like 9 (PhFTL9) as the candidate underlying q-FT5. We used a reciprocal transplant experiment to test for local adaptation and the specific impact of q-FT5 on performance. We did not observe local adaptation in terms of fitness tradeoffs when contrasting ecotypes in home versus away habitats. However, we observed that the coastal allele of q-FT5 conferred a fitness advantage only in its local habitat but not at the inland site. Sequence analyses identified an excess of low-frequency polymorphisms at the PhFTL9 promoter in the inland lineage, suggesting a role for either selection or population expansion on promoter evolution. Together, our findings demonstrate the genetic basis of flowering variation in a perennial grass and provide evidence for conditional neutrality underlying flowering time divergence.

Natural variation identifies new effectors of water-use efficiency in Arabidopsis.

Water-use efficiency (WUE) is the ratio of biomass produced per unit of water consumed; thus, it can be altered by genetic factors that affect either side of the ratio. In the present study, we exploited natural variation for WUE to discover loci affecting either biomass accumulation or water use as factors affecting WUE. Genome-wide association studies (GWAS) using integrated WUE measured through carbon isotope discrimination (δ13C) of Arabidopsis thaliana accessions identified genomic regions associated with WUE. Reverse genetic analysis of 70 candidate genes selected based on the GWAS results and transcriptome data identified 25 genes affecting WUE as measured by gravimetric and δ13C analyses. Mutants of four genes had higher WUE than wild type, while mutants of the other 21 genes had lower WUE. The differences in WUE were caused by either altered biomass or water consumption (or both). Stomatal density (SD) was not a primary cause of altered WUE in these mutants. Leaf surface temperatures indicated that transpiration differed for mutants of 16 genes, but generally biomass accumulation had a greater effect on WUE. The genes we identified are involved in diverse cellular processes, including hormone and calcium signaling, meristematic activity, photosynthesis, flowering time, leaf/vasculature development, and cell wall composition; however, none of them had been previously linked to WUE. Thus, our study successfully identified effectors of WUE that can be used to understand the genetic basis of WUE and improve crop productivity.

Geographic patterns of genomic diversity and structure in the C4 grass Panicum hallii across its natural distribution.

Geographic patterns of within-species genomic diversity are shaped by evolutionary processes, life history and historical and contemporary factors. New genomic approaches can be used to infer the influence of such factors on the current distribution of infraspecific lineages. In this study, we evaluated the genomic and morphological diversity as well as the genetic structure of the C4 grass Panicum hallii across its complex natural distribution in North America. We sampled extensively across the natural range of P. hallii in Mexico and the USA to generate double-digestion restriction-associated DNA (ddRAD) sequence data for 423 individuals from 118 localities. We used these individuals to study the divergence between the two varieties of P. hallii, P. hallii var. filipes and P. hallii var. hallii as well as the genetic diversity and structure within these groups. We also examined the possibility of admixture in the geographically sympatric zone shared by both varieties, and assessed distribution shifts related with past climatic fluctuations. There is strong genetic and morphological divergence between the varieties and consistent genetic structure defining seven genetic clusters that follow major ecoregions across the range. South Texas constitutes a hotspot of genetic diversity with the co-occurrence of all genetic clusters and admixture between the two varieties. It is likely a recolonization and convergence point of populations that previously diverged in isolation during fragmentation events following glaciation periods.

Genomic mechanisms of climate adaptation in polyploid bioenergy switchgrass.

Long-term climate change and periodic environmental extremes threaten food and fuel security1 and global crop productivity2-4. Although molecular and adaptive breeding strategies can buffer the effects of climatic stress and improve crop resilience5, these approaches require sufficient knowledge of the genes that underlie productivity and adaptation6-knowledge that has been limited to a small number of well-studied model systems. Here we present the assembly and annotation of the large and complex genome of the polyploid bioenergy crop switchgrass (Panicum virgatum). Analysis of biomass and survival among 732 resequenced genotypes, which were grown across 10 common gardens that span 1,800 km of latitude, jointly revealed extensive genomic evidence of climate adaptation. Climate-gene-biomass associations were abundant but varied considerably among deeply diverged gene pools. Furthermore, we found that gene flow accelerated climate adaptation during the postglacial colonization of northern habitats through introgression of alleles from a pre-adapted northern gene pool. The polyploid nature of switchgrass also enhanced adaptive potential through the fractionation of gene function, as there was an increased level of heritable genetic diversity on the nondominant subgenome. In addition to investigating patterns of climate adaptation, the genome resources and gene-trait associations developed here provide breeders with the necessary tools to increase switchgrass yield for the sustainable production of bioenergy.

Gene Expression analysis associated with salt stress in a reciprocally crossed rice population.

The rice landrace Horkuch, endemic to the southern saline coast of Bangladesh, is known to have salt tolerance traits and can therefore contribute to a high yielding recipient for breeding purposes. In this study, we reciprocally crossed Horkuch with high yielding but salt sensitive IR29 to detect the complement of genes that were responsible for conferring salt tolerance versus sensitivity at the seedling developmental stage. We looked at tolerant and sensitive F3 families from individual F2 segregating plants and analyzed them for differential gene expressions using RNAseq. In general, we observed higher numbers of genes differentially expressed in leaves compared to root tissues. This included both upregulation and downregulation of gene expression across our experimental factors. Gene expression decreased in sensitive leaf after stress exposure where tolerant plants showed the opposite trend. In root, tolerant plants expression decreased at higher time points of stress exposure. We also observed a strong maternal cytoplasmic effect on gene expression and this was most evident in roots where there was upregulation in functional enrichments related to phosphorylation, electron carriers, transporter and cation transmembrane activities. Stress groups (tolerant and sensitive) response in F3 families were distinctive in both cytoplasmic backgrounds and involved uniquely upregulated genes in tolerant progenies including membrane sensor proteins, enzymes involved with signaling pathways, such as those producing trehalose and G-protein coupled receptor proteins, photosynthesis-related enzymes and golgi body recycling as well as prolamin precursor proteins involved in refolding of proteins. On the other hand, sensitivity was found to be associated with differential upregulation of only a few redox proteins and higher number of apoptosis related genes compared to the tolerant response. Overall, our highly replicated experimental design was powerful and allowed the detection of relatively subtle differential expression. Our future goal is to correlate these expression differences with QTLs in this population, which would help identify the relative importance of specific genetic loci and provide a direct avenue for combining higher levels of salt tolerance with better agronomic traits in rice.

Complex interactions between day length and diurnal patterns of gene expression drive photoperiodic responses in a perennial C4 grass.

Photoperiod is a key environmental cue affecting flowering and biomass traits in plants. Key components of the photoperiodic flowering pathway have been identified in many species, but surprisingly few studies have globally examined the diurnal rhythm of gene expression with changes in day length. Using a cost-effective 3'-Tag RNA sequencing strategy, we characterize 9,010 photoperiod responsive genes with strict statistical testing across a diurnal time series in the C4 perennial grass, Panicum hallii. We show that the vast majority of photoperiod responses are driven by complex interactions between day length and sampling periods. A fine-scale contrast analysis at each sampling time revealed a detailed picture of the temporal reprogramming of cis-regulatory elements and biological processes under short- and long-day conditions. Phase shift analysis reveals quantitative variation among genes with photoperiod-dependent diurnal patterns. In addition, we identify three photoperiod enriched transcription factor families with key genes involved in photoperiod flowering regulatory networks. Finally, coexpression networks analysis of GIGANTEA homolog predicted 1,668 potential coincidence partners, including five well-known GI-interacting proteins. Our results not only provide a resource for understanding the mechanisms of photoperiod regulation in perennial grasses but also lay a foundation to increase biomass yield in biofuel crops.

Reproductive stage physiological and transcriptional responses to salinity stress in reciprocal populations derived from tolerant (Horkuch) and susceptible (IR29) rice.

Global increase in salinity levels has made it imperative to identify novel sources of genetic variation for tolerance traits, especially in rice. The rice landrace Horkuch, endemic to the saline coastal area of Bangladesh, was used in this study as the source of tolerance in reciprocal crosses with the sensitive but high-yielding IR29 variety for discovering transcriptional variation associated with salt tolerance in the resulting populations. The cytoplasmic effect of the Horkuch background in leaves under stress showed functional enrichment for signal transduction, DNA-dependent regulation and transport activities. In roots the enrichment was for cell wall organization and macromolecule biosynthesis. In contrast, the cytoplasmic effect of IR29 showed upregulation of apoptosis and downregulation of phosphorylation across tissues relative to Horkuch. Differential gene expression in leaves of the sensitive population showed downregulation of GO processes like photosynthesis, ATP biosynthesis and ion transport. Roots of the tolerant plants conversely showed upregulation of GO terms like G-protein coupled receptor pathway, membrane potential and cation transport. Furthermore, genes involved in regulating membrane potentials were constitutively expressed only in the roots of tolerant individuals. Overall our work has developed genetic resources and elucidated the likely mechanisms associated with the tolerance response of the Horkuch genotype.

Comparative genomics of two jute species and insight into fibre biogenesis.

Jute (Corchorus sp.) is one of the most important sources of natural fibre, covering ∼80% of global bast fibre production1. Only Corchorus olitorius and Corchorus capsularis are commercially cultivated, though there are more than 100 Corchorus species2 in the Malvaceae family. Here we describe high-quality draft genomes of these two species and their comparisons at the functional genomics level to support tailor-designed breeding. The assemblies cover 91.6% and 82.2% of the estimated genome sizes for C. olitorius and C. capsularis, respectively. In total, 37,031 C. olitorius and 30,096 C. capsularis genes are identified, and most of the genes are validated by cDNA and RNA-seq data. Analyses of clustered gene families and gene collinearity show that jute underwent shared whole-genome duplication ∼18.66 million years (Myr) ago prior to speciation. RNA expression analysis from isolated fibre cells reveals the key regulatory and structural genes involved in fibre formation. This work expands our understanding of the molecular basis of fibre formation laying the foundation for the genetic improvement of jute.

Enhanced Salt Tolerance Conferred by the Complete 2.3 kb cDNA of the Rice Vacuolar Na(+)/H(+) Antiporter Gene Compared to 1.9 kb Coding Region with 5' UTR in Transgenic Lines of Rice.

Soil salinity is one of the most challenging problems that restricts the normal growth and production of rice worldwide. It has therefore become very important to produce more saline tolerant rice varieties. This study shows constitutive over-expression of the vacuolar Na(+)/H(+) antiporter gene (OsNHX1) from the rice landrace (Pokkali) and attainment of enhanced level of salinity tolerance in transgenic rice plants. It also shows that inclusion of the complete un-translated regions (UTRs) of the alternatively spliced OsNHX1 gene provides a higher level of tolerance to the transgenic rice. Two separate transformation events of the OsNHX1 gene, one with 1.9 kb region containing the 5' UTR with CDS and the other of 2.3 kb, including 5' UTR, CDS, and the 3' UTR regions were performed. The transgenic plants with these two different constructs were advanced to the T3 generation and physiological and molecular screening of homozygous plants was conducted at seedling and reproductive stages under salinity (NaCl) stress. Both transgenic lines were observed to be tolerant compared to WT plants at both physiological stages. However, the transgenic lines containing the CDS with both the 5' and 3' UTR were significantly more tolerant compared to the transgenic lines containing OsNHX1 gene without the 3' UTR. At the seedling stage at 12 dS/m stress, the chlorophyll content was significantly higher (P < 0.05) and the electrolyte leakage significantly lower (P < 0.05) in the order 2.3 kb > 1.9 kb > and WT lines. Yield in g/plant in the best line from the 2.3 kb plants was significantly more (P < 0.01) compared, respectively, to the best 1.9 kb line and WT plants at stress of 6 dS/m. Transformation with the complete transcripts rather than the CDS may therefore provide more durable level of tolerance.

Unique Genotypic Differences Discovered among Indigenous Bangladeshi Rice Landraces.

Bangladesh is a reservoir of diverse rice germplasm and is home to many landraces with unique, important traits. Molecular characterization of these landraces is of value for their identification, preservation, and potential use in breeding programs. Thirty-eight rice landraces from different regions of Bangladesh including some high yielding BRRI varieties were analyzed by 34 polymorphic microsatellite markers yielding a total of 258 reproducible alleles. The analysis could locate 34 unique identifiers for 21 genotypes, making the latter potentially amenable to identity verification. An identity map for these genotypes was constructed with all the 12 chromosomes of the rice genome. Polymorphism information content (PIC) scores of the 34 SSR markers were 0.098 to 0.89 where on average 7.5 alleles were observed. A dendogram constructed using UPGMA clustered the varieties into two major groups and five subgroups. In some cases, the clustering matched with properties like aromaticity, stickiness, salt tolerance, and photoperiod insensitivity. The results will help breeders to work towards the proper utilization of these landraces for parental selection and linkage map construction for discovery of useful alleles.