The Road to Sorghum Domestication: Evidence From Nucleotide Diversity and Gene Expression Patterns.

Native African cereals (sorghum, millets) ensure food security to millions of low-income people from low fertility and drought-prone regions of Africa and Asia. In spite of their agronomic importance, the genetic bases of their phenotype and adaptations are still not well-understood. Here we focus on Sorghum bicolor, which is the fifth cereal worldwide for grain production and constitutes the staple food for around 500 million people. We leverage transcriptomic resources to address the adaptive consequences of the domestication process. Gene expression and nucleotide variability were analyzed in 11 domesticated and nine wild accessions. We documented a downregulation of expression and a reduction of diversity both in nucleotide polymorphism (30%) and gene expression levels (18%) in domesticated sorghum. These findings at the genome-wide level support the occurrence of a global reduction of diversity during the domestication process, although several genes also showed patterns consistent with the action of selection. Nine hundred and forty-nine genes were significantly differentially expressed between wild and domesticated gene pools. Their functional annotation points to metabolic pathways most likely contributing to the sorghum domestication syndrome, such as photosynthesis and auxin metabolism. Coexpression network analyzes revealed 21 clusters of genes sharing similar expression patterns. Four clusters (totaling 2,449 genes) were significantly enriched in differentially expressed genes between the wild and domesticated pools and two were also enriched in domestication and improvement genes previously identified in sorghum. These findings reinforce the evidence that the combined and intricated effects of the domestication and improvement processes do not only affect the behaviors of a few genes but led to a large rewiring of the transcriptome. Overall, these analyzes pave the way toward the identification of key domestication genes valuable for genetic resources characterization and breeding purposes.

Transcriptional Regulation of Sorghum Stem Composition: Key Players Identified Through Co-expression Gene Network and Comparative Genomics Analyses.

Most sorghum biomass accumulates in stem secondary cell walls (SCW). As sorghum stems are used as raw materials for various purposes such as feed, energy and fiber reinforced polymers, identifying the genes responsible for SCW establishment is highly important. Taking advantage of studies performed in model species, most of the structural genes contributing at the molecular level to the SCW biosynthesis in sorghum have been proposed while their regulatory factors have mostly not been determined. Validation of the role of several MYB and NAC transcription factors in SCW regulation in Arabidopsis and a few other species has been provided. In this study, we contributed to the recent efforts made in grasses to uncover the mechanisms underlying SCW establishment. We reported updated phylogenies of NAC and MYB in 9 different species and exploited findings from other species to highlight candidate regulators of SCW in sorghum. We acquired expression data during sorghum internode development and used co-expression analyses to determine groups of co-expressed genes that are likely to be involved in SCW establishment. We were able to identify two groups of co-expressed genes presenting multiple evidences of involvement in SCW building. Gene enrichment analysis of MYB and NAC genes provided evidence that while NAC SECONDARY WALL THICKENING PROMOTING FACTOR NST genes and SECONDARY WALL-ASSOCIATED NAC DOMAIN PROTEIN gene functions appear to be conserved in sorghum, NAC master regulators of SCW in sorghum may not be as tissue compartmentalized as in Arabidopsis. We showed that for every homolog of the key SCW MYB in Arabidopsis, a similar role is expected for sorghum. In addition, we unveiled sorghum MYB and NAC that have not been identified to date as being involved in cell wall regulation. Although specific validation of the MYB and NAC genes uncovered in this study is needed, we provide a network of sorghum genes involved in SCW both at the structural and regulatory levels.

Transcriptome data from three endemic Myrtaceae species from New Caledonia displaying contrasting responses to myrtle rust (Austropuccinia psidii).

The myrtle rust disease, caused by the fungus Austropuccinia psidii, infects a wide range of host species within the Myrtaceae family worldwide. Since its first report in 2013 in New Caledonia, it was found on various types of native environments where Myrtaceae are the dominant or codominant species, as well as in several commercial nurseries. It is now considered as a significant threat to ecosystems biodiversity and Myrtaceae-related economy. The use of predictive molecular markers for resistance against myrtle rust is currently the most cost-effective and ecological approach to control the disease. Such an approach for neo Caledonian endemic Myrtaceae species was not possible because of the lack of genomic resources. The recent advancement in new generation sequencing technologies accompanied with relevant bioinformatics tools now provide new research opportunity for work in non-model organism at the transcriptomic level. The present study focuses on transcriptome analysis on three Myrtaceae species endemic to New Caledonia (Arillastrum gummiferum, Syzygium longifolium and Tristaniopsis glauca) that display contrasting responses to the pathogen (non-infected vs infected). Differential gene expression (DGE) and variant calling analysis were conducted on each species. We combined a dual approach by using 1) the annotated reference genome of a related Myrtaceae species (Eucalyptus grandis) and 2) a de novo transcriptomes of each species.

Evolutionary forces affecting synonymous variations in plant genomes.

Base composition is highly variable among and within plant genomes, especially at third codon positions, ranging from GC-poor and homogeneous species to GC-rich and highly heterogeneous ones (particularly Monocots). Consequently, synonymous codon usage is biased in most species, even when base composition is relatively homogeneous. The causes of these variations are still under debate, with three main forces being possibly involved: mutational bias, selection and GC-biased gene conversion (gBGC). So far, both selection and gBGC have been detected in some species but how their relative strength varies among and within species remains unclear. Population genetics approaches allow to jointly estimating the intensity of selection, gBGC and mutational bias. We extended a recently developed method and applied it to a large population genomic dataset based on transcriptome sequencing of 11 angiosperm species spread across the phylogeny. We found that at synonymous positions, base composition is far from mutation-drift equilibrium in most genomes and that gBGC is a widespread and stronger process than selection. gBGC could strongly contribute to base composition variation among plant species, implying that it should be taken into account in plant genome analyses, especially for GC-rich ones.

Small RNA profiling reveals regulation of Arabidopsis miR168 and heterochromatic siRNA415 in response to fungal elicitors.

BACKGROUND: Small RNAs (sRNAs), including small interfering RNAs (siRNAs) and microRNAs (miRNAs), have emerged as important regulators of eukaryotic gene expression. In plants, miRNAs play critical roles in development, nutrient homeostasis and abiotic stress responses. Accumulating evidence also reveals that sRNAs are involved in plant immunity. Most studies on pathogen-regulated sRNAs have been conducted in Arabidopsis plants infected with the bacterial pathogen Pseudomonas syringae, or treated with the flagelin-derived elicitor peptide flg22 from P. syringae. This work investigates sRNAs that are regulated by elicitors from the fungus Fusarium oxysporum in Arabidopsis. RESULTS: Microarray analysis revealed alterations on the accumulation of a set of sRNAs in response to elicitor treatment, including miRNAs and small RNA sequences derived from massively parallel signature sequencing. Among the elicitor-regulated miRNAs was miR168 which regulates ARGONAUTE1, the core component of the RNA-induced silencing complex involved in miRNA functioning. Promoter analysis in transgenic Arabidopsis plants revealed transcriptional activation of MIR168 by fungal elicitors. Furthermore, transgenic plants expressing a GFP-miR168 sensor gene confirmed that the elicitor-induced miR168 is active. MiR823, targeting Chromomethylase3 (CMT3) involved in RNA-directed DNA methylation (RdDM) was also found to be regulated by fungal elicitors. In addition to known miRNAs, microarray analysis allowed the identification of an elicitor-inducible small RNA that was incorrectly annotated as a miRNA. Studies on Arabidopsis mutants impaired in small RNA biogenesis demonstrated that this sRNA, is a heterochromatic-siRNA (hc-siRNA) named as siRNA415. Hc-siRNAs are known to be involved in RNA-directed DNA methylation (RdDM). SiRNA415 is detected in several plant species. CONCLUSION: Results here presented support a transcriptional regulatory mechanism underlying MIR168 expression. This finding highlights the importance of miRNA functioning in adaptive processes of Arabidopsis plants to fungal infection. The results of this study also lay a foundation for the involvement of RdDM processes through the activity of siRNA415 and miR823 in mediating regulation of immune responses in Arabidopsis plants.

Characterization of a cassiicolin-encoding gene from Corynespora cassiicola, pathogen of rubber tree (Hevea brasiliensis).

Corynespora Leaf Fall (CLF) is a major disease of rubber tree (Hevea brasiliensis) caused by the Ascomycota Corynespora cassiicola. Here we describe the cloning and characterization of a gene encoding cassiicolin (Cas), a glycosylated cystein-rich small secreted protein (SSP) identified as a potential CLF disease effector in rubber tree. Three isolates with contrasted levels of aggressiveness were analyzed comparatively. The cassiicolin gene was detected - and the toxin successfully purified - from the isolates with high and medium aggressiveness (CCP and CCAM3 respectively) but not from the isolate with the lowest aggressiveness (CCAM1), suggesting the existence of a different disease effector in the later. CCP and CCAM3 carried strictly identical cassiicolin genes and produced toxins of identical mass, as evidence by mass spectrometry analysis, thus suggesting conserved post-translational modifications in addition to sequence identity. The differences in aggressiveness between CCP and CCAM3 may be attributed to differences in cassiicolin transcript levels rather than qualitative variations in cassiicolin structure. Cassiicolin may play an important role in the early phase of infection since a peak of cassiicolin transcripts occurred in 1 or 2 days after inoculation (before the occurrence of the first symptoms), in both the tolerant and the susceptible cultivars.

Characterization of a defensin gene expressed in oil palm inflorescences: induction during tissue culture and possible association with epigenetic somaclonal variation events.

From differential display studies performed on oil palm (Elaeis guineensis Jacq.) tissue cultures bearing or lacking an epigenetic homeotic flowering abnormality, known as mantled, EGAD1, a gene coding for a putative plant defensin, has been identified and characterized. In whole plants, transcripts of the EGAD1 gene were detected only in inflorescences. The closest characterized relative of the oil palm EGAD1 gene is the Petunia PPT gene, which is expressed principally in the pistil of the flower. The 77 amino acid polypeptide encoded by the EGAD1 gene displays strong similarities with a number of plant defensin proteins, which are thought to play a protective role and which have been shown in some cases to possess antifungal properties. Oil palm tissue cultures exhibit a generally strong induction of accumulation of EGAD1 transcripts, which were detected to differing extents at all stages of the tissue culture regeneration process. The 5' flanking region of the EGAD1 gene was found to contain two different types of potential cis-acting DNA element previously identified in the promoters of plant defence-related genes, which may explain the observed expression in tissue cultures. At the callus stage of the in vitro regeneration procedure, a differential accumulation of EGAD1 transcripts was observed which correlated with the presence or absence of the mantled flowering abnormality. EGAD1 gene expression may therefore be a marker of epigenetic somaclonal variation events.