Maize (Zea mays L.) genome diversity as revealed by RNA-sequencing

Efficient mutagenesis and genotyping of maize inbreds using biolistics, multiplex CRISPR/Cas9 editing, and Indel-Selective PCR

CRISPR/Cas9 based genome editing has advanced our understanding of a myriad of important biological phenomena. Important challenges to multiplex genome editing in maize include assembly of large complex DNA constructs, few genotypes with efficient transformation systems, and costly/labor-intensive genotyping methods. Here we present an approach for multiplex CRISPR/Cas9 genome editing system that delivers a single compact DNA construct via biolistics to Type I embryogenic calli, followed by a novel efficient genotyping assay to identify desirable editing outcomes. We first demonstrate the creation of heritable mutations at multiple target sites within the same gene. Next, we successfully created individual and stacked mutations for multiple members of a gene family. Genome sequencing found off-target mutations are rare. Multiplex genome editing was achieved for both the highly transformable inbred line H99 and Illinois Low Protein1 (ILP1), a genotype where transformation has not previously been reported. In addition to screening transformation events for deletion alleles by PCR, we also designed PCR assays that selectively amplify deletion or insertion of a single nucleotide, the most common outcome from DNA repair of CRISPR/Cas9 breaks by non-homologous end-joining. The Indel-Selective PCR (IS-PCR) method enabled rapid tracking of multiple edited alleles in progeny populations. The 'end to end' pipeline presented here for multiplexed CRISPR/Cas9 mutagenesis can be applied to accelerate maize functional genomics in a broader diversity of genetic backgrounds.

Identification and characterisation of 'No apical meristem; Arabidopsis transcription activation factor; Cup-shape cotyledon' (NAC) family transcription factors involved in sugar accumulation and abscisic acid signalling in grape (Vitis vinifera)

The 'No apical meristem; Arabidopsis transcription activation factor; Cup-shape cotyledon' (NAC) transcription factors are pivotal in plant development and stress response. Sucrose-non-fermenting-related protein kinase 1.2 (SnRK1) is a key enzyme in glucose metabolism and ABA signalling. In this study, we used grape (Vitis vinifera ) calli to explore NAC's roles in sugar and ABA pathways and its relationship with VvSnRK1.2 . We identified 19 VvNACs highly expressed at 90days after blooming, coinciding with grape maturity and high sugar accumulation, and 11 VvNACs randomly selected from 19 were demonstrated in response to sugar and ABA treatments. VvNAC26 showed significant response to sugar and ABA treatments, and its protein, as a nucleus protein, had transcriptional activation in yeast. We obtained the overexpression (OE-VvNAC26 ) and RNA-inhibition (RNAi-VvNAC26 ) of VvNAC26 in transgenic calli by Agrobacterium tumefaciens -mediated transformation. We found that VvNAC26 negatively influenced fructose content. Under sugar and ABA treatments, VvNAC26 negatively influenced the expression of most sugar-related genes, while positively influencing the expression of most ABA pathway-related genes. Dual-luciferase reporter experiments demonstrated that VvNAC26 significantly upregulates VvSnRK1.2 promoter expression in tobacco (Nicotiana benthamiana ) leaves, although this process in grape calli requires ABA. The levels of sugar content, sugar-related genes, and ABA-related genes fluctuated significantly in OE-VvNAC26 +RNAi-VvSnRK1.2 and OE-VvSnRK1.2 +RNAi-VvNAC26 transgenic calli. These findings indicated that VvNAC26 regulates sugar metabolism and ABA pathway, displaying synergistic interactions with VvSnRK1.2 .

A complete telomere-to-telomere assembly of the maize genome

A complete telomere-to-telomere (T2T) finished genome has been the long pursuit of genomic research. Through generating deep coverage ultralong Oxford Nanopore Technology (ONT) and PacBio HiFi reads, we report here a complete genome assembly of maize with each chromosome entirely traversed in a single contig. The 2,178.6 Mb T2T Mo17 genome with a base accuracy of over 99.99% unveiled the structural features of all repetitive regions of the genome. There were several super-long simple-sequence-repeat arrays having consecutive thymine-adenine-guanine (TAG) tri-nucleotide repeats up to 235 kb. The assembly of the entire nucleolar organizer region of the 26.8 Mb array with 2,974 45S rDNA copies revealed the enormously complex patterns of rDNA duplications and transposon insertions. Additionally, complete assemblies of all ten centromeres enabled us to precisely dissect the repeat compositions of both CentC-rich and CentC-poor centromeres. The complete Mo17 genome represents a major step forward in understanding the complexity of the highly recalcitrant repetitive regions of higher plant genomes.

Transcriptome analysis revealed the expression levels of genes related to abscisic acid and auxin biosynthesis in grapevine (Vitis vinifera L.) under root restriction

The root system is essential for the stable growth of plants. Roots help anchor plants in the soil and play a crucial role in water uptake, mineral nutrient absorption and endogenous phytohormone formation. Root-restriction (RR) cultivation, a powerful technique, confines plant roots to a specific soil space. In the present study, roots of one-year-old "Muscat Hamburg" grapevine under RR and control (nR) treatments harvested at 70 and 125 days after planting were used for transcriptome sequencing, and in total, 2031 (nR7 vs. nR12), 1445 (RR7 vs. RR12), 1532 (nR7 vs. RR7), and 2799 (nR12 vs. RR12) differentially expressed genes (DEGs) were identified. Gene Ontology (GO) enrichment analysis demonstrated that there were several genes involved in the response to different phytohormones, including abscisic acid (ABA), auxin (IAA), ethylene (ETH), gibberellins (GAs), and cytokinins (CTKs). Among them, multiple genes, such as PIN2 and ERF113, are involved in regulating vital plant movements by various phytohormone pathways. Moreover, following RR cultivation, DEGs were enriched in the biological processes of plant-type secondary cell wall biosynthesis, the defense response, programmed cell death involved in cell development, and the oxalate metabolic process. Furthermore, through a combined analysis of the transcriptome and previously published microRNA (miRNA) sequencing results, we found that multiple differentially expressed miRNAs (DEMs) and DEG combinations in different comparison groups exhibited opposite trends, indicating that the expression levels of miRNAs and their target genes were negatively correlated. Furthermore, RR treatment indeed significantly increased the ABA content at 125 days after planting and significantly decreased the IAA content at 70 days after planting. Under RR cultivation, most ABA biosynthesis-related genes were upregulated, while most IAA biosynthesis-related genes were downregulated. These findings lay a solid foundation for further establishing the network through which miRNAs regulate grapevine root development through target genes and for further exploring the molecular mechanism through which endogenous ABA and IAA regulate root architecture development in grapevine.

Single-parent expression complementation contributes to phenotypic heterosis in maize hybrids

The dominance model of heterosis explains the superior performance of F1-hybrids via the complementation of deleterious alleles by beneficial alleles in many genes. Genes active in one parent but inactive in the second lead to single-parent expression (SPE) complementation in maize (Zea mays L.) hybrids. In this study, SPE complementation resulted in approximately 700 additionally active genes in different tissues of genetically diverse maize hybrids on average. We established that the number of SPE genes is significantly associated with mid-parent heterosis (MPH) for all surveyed phenotypic traits. In addition, we highlighted that maternally (SPE_B) and paternally (SPE_X) active SPE genes enriched in gene co-expression modules are highly correlated within each SPE type but separated between these two SPE types. While SPE_B-enriched co-expression modules are positively correlated with phenotypic traits, SPE_X-enriched modules displayed a negative correlation. Gene ontology term enrichment analyses indicated that SPE_B patterns are associated with growth and development, whereas SPE_X patterns are enriched in defense and stress response. In summary, these results link the degree of phenotypic MPH to the prevalence of gene expression complementation observed by SPE, supporting the notion that hybrids benefit from SPE complementation via its role in coordinating maize development in fluctuating environments.

Standardized genome-wide function prediction enables comparative functional genomics: a new application area for Gene Ontologies in plants

Background

Genome-wide gene function annotations are useful for hypothesis generation and for prioritizing candidate genes potentially responsible for phenotypes of interest. We functionally annotated the genes of 18 crop plant genomes across 14 species using the GOMAP pipeline.

Results

By comparison to existing GO annotation datasets, GOMAP-generated datasets cover more genes, contain more GO terms, and are similar in quality (based on precision and recall metrics using existing gold standards as the basis for comparison). From there, we sought to determine whether the datasets across multiple species could be used together to carry out comparative functional genomics analyses in plants. To test the idea and as a proof of concept, we created dendrograms of functional relatedness based on terms assigned for all 18 genomes. These dendrograms were compared to well-established species-level evolutionary phylogenies to determine whether trees derived were in agreement with known evolutionary relationships, which they largely are. Where discrepancies were observed, we determined branch support based on jackknifing then removed individual annotation sets by genome to identify the annotation sets causing unexpected relationships.

Conclusions

GOMAP-derived functional annotations used together across multiple species generally retain sufficient biological signal to recover known phylogenetic relationships based on genome-wide functional similarities, indicating that comparative functional genomics across species based on GO data holds promise for generating novel hypotheses about comparative gene function and traits.

Primary Root and Mesocotyl Elongation in Maize Seedlings: Two Organs with Antagonistic Growth below the Soil Surface

Maize illustrates one of the most complex cases of embryogenesis in higher plants that results in the development of early embryo with distinctive organs such as the mesocotyl, seminal and primary roots, coleoptile, and plumule. After seed germination, the elongation of root and mesocotyl follows opposite directions in response to specific tropisms (positive and negative gravitropism and hydrotropism). Tropisms represent the differential growth of an organ directed toward several stimuli. Although the life cycle of roots and mesocotyl takes place in darkness, their growth and functions are controlled by different mechanisms. Roots ramify through the soil following the direction of the gravity vector, spreading their tips into new territories looking for water; when water availability is low, the root hydrotropic response is triggered toward the zone with higher moisture. Nonetheless, there is a high range of hydrotropic curvatures (angles) in maize. The processes that control root hydrotropism and mesocotyl elongation remain unclear; however, they are influenced by genetic and environmental cues to guide their growth for optimizing early seedling vigor. Roots and mesocotyls are crucial for the establishment, growth, and development of the plant since both help to forage water in the soil. Mesocotyl elongation is associated with an ancient agriculture practice known as deep planting. This tradition takes advantage of residual soil humidity and continues to be used in semiarid regions of Mexico and USA. Due to the genetic diversity of maize, some lines have developed long mesocotyls capable of deep planting while others are unable to do it. Hence, the genetic and phenetic interaction of maize lines with a robust hydrotropic response and higher mesocotyl elongation in response to water scarcity in time of global heating might be used for developing more resilient maize plants.

Transcriptome analysis reveals the mechanism associated with dynamic changes in fatty acid and phytosterol content in foxtail millet (Setaria italica) during seed development

Foxtail millet (Setaria italica) is an excellent source of beneficial natural fatty acids and phytosterols. However, the mechanisms underlying the dynamic changes of fatty acids and phytosterols during seed development are unknown. In this study, a comprehensive dynamic change analysis of the bioactive compounds during seed development was conducted in two cultivars with different crude fat content (high-fat, JG 35 [5.40%]; and low-fat, JG 39 [2.90%]). GC-FID/MS analysis showed that the proportion of unsaturated fatty acids (UFAs) were higher than the saturated fatty acids (SFAs). UFA content first increased, then decreased during seed development, while SFA content showed the opposite trend. Oil contents continuously increased with seed development, especially at the S2 stage. Phytosterol contents initially increased, then decreased with seed development. Transcriptome analysis revealed that 152 genes were associated with fatty acid metabolism and phytosterol biosynthesis, of which 46 and 62 were related to UFA and phytosterol biosynthesis, respectively. Furthermore, the key genes involved in fatty acid synthesis (ACCase and FATA/B), triacylglycerol biosynthesis (LACS, GPAT, and DGAT), and phytosterols synthesis (CAS1, STM1, EGR6, and DWF1) were overexpressed. This led to maximum UFA, oil, and phytosterol accumulation in JG 35 at the S2 stage. This study reveals the mechanism behind the dynamic changes of fatty acid and phytosterol contents in foxtail millet during seed development.

Single-parent expression drives dynamic gene expression complementation in maize hybrids

Single-parent expression (SPE) is defined as gene expression in only one of the two parents. SPE can arise from differential expression between parental alleles, termed non-presence/absence (non-PAV) SPE, or from the physical absence of a gene in one parent, termed PAV SPE. We used transcriptome data of diverse Zea mays (maize) inbreds and hybrids, including 401 samples from five different tissues, to test for differences between these types of SPE genes. Although commonly observed, SPE is highly genotype and tissue specific. A positive correlation was observed between the genetic distance of the two inbred parents and the number of SPE genes identified. Regulatory analysis showed that PAV SPE and non-PAV SPE genes are mainly regulated by cis effects, with a small fraction under trans regulation. Polymorphic transposable element insertions in promoter sequences contributed to the high level of cis regulation for PAV SPE and non-PAV SPE genes. PAV SPE genes were more frequently expressed in hybrids than non-PAV SPE genes. The expression of parentally silent alleles in hybrids of non-PAV SPE genes was relatively rare but occurred in most hybrids. Non-PAV SPE genes with expression of the silent allele in hybrids are more likely to exhibit above high parent expression level than hybrids that do not express the silent allele, leading to non-additive expression. This study provides a comprehensive understanding of the nature of non-PAV SPE and PAV SPE genes and their roles in gene expression complementation in maize hybrids.

Multiple-genotypes transcriptional analysis revealed candidates genes and nucleotide variants for improvement of quality characteristics in tea (Camellia sinensis (L.) O. Kuntze)

Tea quality is a polygenic trait that exhibits tremendous genetic variability due to accumulation of array of secondary metabolites. To elucidate global molecular insights controlling quality attributes, metabolite profiling and transcriptome sequencing of twelve diverse tea cultivars was performed in tea shoots harvested during quality season. RP-HPLC-DAD analysis of quality parameters revealed significant difference in catechins, theanine and caffeine contents. Transcriptome sequencing resulted into 50,107 non-redundant transcripts with functional annotations of 81.6% (40,847) of the transcripts. Interestingly, 2872 differentially expressed transcripts exhibited significant enrichment in 38 pathways (FDR ≤ 0.05) including secondary metabolism, amino acid and carbon metabolism. Thirty-eight key candidates reportedly involved in biosynthesis of fatty acid derived volatiles, volatile terpenes, glycoside hydrolysis and key quality related pathways (flavonoid, caffeine and theanine-biosynthesis) were highly expressed in catechins-rich tea cultivars. Furthermore, enrichment of candidates involved in flavonoid biosynthesis, transcriptional regulation, volatile terpene and biosynthesis of fatty acid derived volatile in Protein-Protein Interactome network revealed well-coordinated regulation of quality characteristics in tea. Additionally, ascertainment of 23,649 non-synonymous SNPs and validation of candidate SNPs present in quality related genes suggests their potential utility in genome-wide mapping and marker development for expediting breeding of elite compound-rich tea cultivars.

Intricate genetic variation networks control the adventitious root growth angle in apple

Background

The root growth angle (RGA) typically determines plant rooting depth, which is significant for plant anchorage and abiotic stress tolerance. Several quantitative trait loci (QTLs) for RGA have been identified in crops. However, the underlying mechanisms of the RGA remain poorly understood, especially in apple rootstocks. The objective of this study was to identify QTLs, validate genetic variation networks, and develop molecular markers for the RGA in apple rootstock.

Results

Bulked segregant analysis by sequencing (BSA-seq) identified 25 QTLs for RGA using 1955 hybrids of the apple rootstock cultivars ‘Baleng Crab’ (Malus robusta Rehd., large RGA) and ‘M9’ (M. pumila Mill., small RGA). With RNA sequencing (RNA-seq) and parental resequencing, six major functional genes were identified and constituted two genetic variation networks for the RGA. Two single nucleotide polymorphisms (SNPs) of the MdLAZY1 promoter damaged the binding sites of MdDREB2A and MdHSFB3, while one SNP of MdDREB2A and MdIAA1 affected the interactions of MdDREB2A/MdHSFB3 and MdIAA1/MdLAZY1, respectively. A SNP within the MdNPR5 promoter damaged the interaction between MdNPR5 and MdLBD41, while one SNP of MdLBD41 interrupted the MdLBD41/MdbHLH48 interaction that affected the binding ability of MdLBD41 on the MdNPR5 promoter. Twenty six SNP markers were designed on candidate genes in each QTL interval, and the marker effects varied from 0.22°-26.11°.

Conclusions

Six diagnostic markers, SNP592, G122, b13, Z312, S1272, and S1288, were used to identify two intricate genetic variation networks that control the RGA and may provide new insights into the accuracy of the molecular markers. The QTLs and SNP markers can potentially be used to select deep-rooted apple rootstocks.

The neonicotinoid insecticide thiamethoxam enhances expression of stress-response genes in Zea mays in an environmentally specific pattern

Recent studies indicate that thiamethoxam (TMX), a neonicotinoid insecticide, can affect plant responses to environmental stressors, such as neighboring weeds. The molecular mechanisms behind both stable and environmentally specific responses to TMX likely involve genes related to defense and stress responses. We investigated the effect of a TMX seed treatment on global gene expression in maize coleoptiles both under normal conditions and under low ratio red to far-red (R:FR) light stress induced by the presence of neighboring plants. The neighboring plant treatment upregulated genes involved in biotic and abiotic stress responses and affected specific photosynthesis and cell-growth related genes. Low R:FR light may enhance maize resistance to herbivores and pathogens. TMX appears to compromise resistance. The TMX treatment stably repressed many genes that encode proteins involved in biotic stress responses, as well as cell-growth genes. Notably, TMX effects on many genes' expression were conditional on the environment. In response to low R:FR, plants treated with TMX engage genes in the JA pathway, as well as other stress-related response pathways. Neighboring weeds may condition TMX-treated plants to become more stress tolerant.

Genome-wide transcriptome and proteome profiles indicate an active role of alternative splicing during de-etiolation of maize seedlings

AS events affect genes encoding protein domain composition and make the single gene produce more proteins with a certain number of genes to satisfy the establishment of photosynthesis during de-etiolation. The drastic switch from skotomorphogenic to photomorphogenic development is an excellent system to elucidate rapid developmental responses to environmental stimuli in plants. To decipher the effects of different light wavelengths on de-etiolation, we illuminated etiolated maize seedlings with blue, red, blue-red mixed and white light, respectively. We found that blue light alone has the strongest effect on photomorphogenesis and that this effect can be attributed to the higher number and expression levels of photosynthesis and chlorosynthesis proteins. Deep sequencing-based transcriptome analysis revealed gene expression changes under different light treatments and a genome-wide alteration in alternative splicing (AS) profiles. We discovered 41,188 novel transcript isoforms for annotated genes, which increases the percentage of multi-exon genes with AS to 63% in maize. We provide peptide support for all defined types of AS, especially retained introns. Further in silico prediction revealed that 58.2% of retained introns have changes in domains compared with their most similar annotated protein isoform. This suggests that AS acts as a protein function switch allowing rapid light response through the addition or removal of functional domains. The richness of novel transcripts and protein isoforms also demonstrates the potential and importance of integrating proteomics into genome annotation in maize.

Characterization of Genetic and Allelic Diversity Amongst Cultivated and Wild Lentil Accessions for Germplasm Enhancement

Intensive breeding of cultivated lentil has resulted in a relatively narrow genetic base, which limits the options to increase crop productivity through selection. Assessment of genetic diversity in the wild gene pool of lentil, as well as characterization of useful and novel alleles/genes that can be introgressed into elite germplasm, presents new opportunities and pathways for germplasm enhancement, followed by successful crop improvement. In the current study, a lentil collection consisting of 467 wild and cultivated accessions that originated from 10 diverse geographical regions was assessed, to understand genetic relationships among different lentil species/subspecies. A total of 422,101 high-confidence SNP markers were identified against the reference lentil genome (cv. CDC Redberry). Phylogenetic analysis clustered the germplasm collection into four groups, namely, Lens culinaris/Lens orientalis, Lens lamottei/Lens odemensis, Lens ervoides, and Lens nigricans. A weak correlation was observed between geographical origin and genetic relationship, except for some accessions of L. culinaris and L. ervoides. Genetic distance matrices revealed a comparable level of variation within the gene pools of L. culinaris (Nei’s coefficient 0.01468–0.71163), L. ervoides (Nei’s coefficient 0.01807–0.71877), and L. nigricans (Nei’s coefficient 0.02188–1.2219). In order to understand any genic differences at species/subspecies level, allele frequencies were calculated from a subset of 263 lentil accessions. Among all cultivated and wild lentil species, L. nigricans exhibited the greatest allelic differentiation across the genome compared to all other species/subspecies. Major differences were observed on six genomic regions with the largest being on Chromosome 1 (c. 1 Mbp). These results indicate that L. nigricans is the most distantly related to L. culinaris and additional structural variations are likely to be identified from genome sequencing studies. This would provide further insights into evolutionary relationships between cultivated and wild lentil germplasm, for germplasm improvement and introgression.

Transcriptome Analysis by RNA-Seq Reveals Genes Related to Plant Height in Two Sets of Parent-hybrid Combinations in Easter lily (Lilium longiflorum)

In this study, two different hybrids of Easter lily (Lilium longiflorum), obtained from two cross combinations, along with their four parents were sequenced by high–throughput RNA–sequencing (RNA–Seq) to find out differentially expressed gene in parent-hybrid combinations. The leaf mRNA profiles of two hybrids and their four parents were RNA–sequenced with a view to identify the potential candidate genes related to plant height heterosis. In both cross combinations, based to morphological traits mid–parent heterosis (MPH) was higher than high–parent heterosis (HPH) for plant height, leaf length, and number of flowers whereas HPH was higher than MPH for flowering time. A total of 4,327 differentially expressed genes (DEGs) were identified through RNA–Seq between the hybrids and their parents based on fold changes (FC) ≥ 2 for up– and ≤ –2 for down–regulation. Venn diagram analysis revealed that there were 703 common DEGs in two hybrid combinations, those were either up– or down–regulated. Most of the commonly expressed DEGs exhibited higher non–additive effects especially overdominance (75.9%) rather than additive (19.4%) and dominance (4.76%) effects. Among the 384 functionally annotated DEGs identified through Blast2GO tool, 12 DEGs were up–regulated and 16 of them were down–regulated in a similar fashion in both hybrids as revealed by heat map analysis. These 28 universally expressed DEGs were found to encode different types of proteins and enzymes those might regulate heterosis by modulating growth, development and stress–related functions in lily. In addition, gene ontology (GO) analysis of 260 annotated DEGs revealed that biological process might play dominant role in heterotic expression. In this first report of transcriptome sequencing in Easter lily, the notable universally up-regulated DEGs annotated ABC transporter A family member–like, B3 domain–containing, disease resistance RPP13/1, auxin–responsive SAUR68–like, and vicilin–like antimicrobial peptides 2–2 proteins those were perhaps associated with plant height heterosis. The genes expressed universally due to their overdominace function perhaps influenced MPH for greater plant height― largely by modulating biological processes involved therein. The genes identified in this study might be exploited in heterosis breeding for plant height of L. longiflorum.

The seasonal development dynamics of the yak hair cycle transcriptome

Background

Mammalian hair play an important role in mammals’ ability to adapt to changing climatic environments. The seasonal circulation of yak hair helps them adapt to high altitude but the regulation mechanisms of the proliferation and differentiation of hair follicles (HFs) cells during development are still unknown. Here, using time series data for transcriptome and hormone contents, we systematically analyzed the mechanism regulating the periodic expression of hair development in the yak and reviewed how different combinations of genetic pathways regulate HFs development and cycling.

Results

This study used high-throughput RNA sequencing to provide a detailed description of global gene expression in 15 samples from five developmental time points during the yak hair cycle. According to clustering analysis, we found that these 15 samples could be significantly grouped into three phases, which represent different developmental periods in the hair cycle. A total of 2316 genes were identified in these three consecutive developmental periods and their expression patterns could be divided into 9 clusters. In the anagen, genes involved in activating hair follicle growth are highly expressed, such as the WNT pathway, FGF pathway, and some genes related to hair follicle differentiation. In the catagen, genes that inhibit differentiation and promote hair follicle cell apoptosis are highly expressed, such as BMP4, and Wise. In the telogen, genes that inhibit hair follicle activity are highly expressed, such as DKK1 and BMP1. Through co-expression analysis, we revealed a number of modular hub genes highly associated with hormones, such as SLF2, BOP1 and DPP8. They may play unique roles in hormonal regulation of events associated with the hair cycle.

Conclusions

Our results revealed the expression pattern and molecular mechanisms of the seasonal hair cycle in the yak. The findings will be valuable in further understanding the alpine adaptation mechanism in the yak, which is important in order to make full use of yak hair resources and promote the economic development of pastoral plateau areas.

Comparative analysis of sugarcane root transcriptome in response to the plant growth-promoting Burkholderia anthina MYSP113

The diazotrophic Burkholderia anthina MYSP113 is a vital plant growth-promoting bacteria and sugarcane root association. The present study based on a detailed analysis of sugarcane root transcriptome by using the HiSeq-Illumina platform in response to the strain MYSP113. The bacterium was initially isolated from the rhizosphere of sugarcane. To better understand biological, cellular, and molecular mechanisms, a de novo transcriptomic assembly of sugarcane root was performed. HiSeq-Illumina platformwas employed for the sequencing of an overall of 16 libraries at a 2×100 bp configuration. Differentially expressed genes (DEGs) analysis identified altered gene expression in 370 genes (total of 199 up-regulated genes and 171 down-regulated genes). Deciphering the huge datasets, concerning the functioning and production of biological systems, a high throughput genome sequencing analysis was attempted with Gene ontology functional analyses and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The report revealed a total of 148930 unigenes. 70414 (47.28%) of them were annotated successfully to Gene Ontology (GO) terms. 774 at 45 days, 4985 of 30 days and 15 days of 6846 terms were significantly regulated. GO analysis revealed that many genes involved in the metabolic, oxidation-reduction process and biological regulatory processes in response to strain MYSP113 and significantly enriched as compare to the control. Moreover, KEGG enriched results show that differentially expressed genes were classified into different pathway categories involved in various processes, such as nitrogen metabolism, plant hormone signal transduction, etc. The sample correlation analyses could help examine the similarity at the gene expression level. The reliability of the observed differential gene expression patterns was validated with quantitative real-time PCR (qRT-PCR). Additionally, plant enzymes activities such as peroxidase and superoxide dismutase were significantly increased in plant roots after the inoculation of strain MYSP113. The results of the research may help in understanding the plant growth-promoting rhizobacteria and plant interaction.

Pan-genomics in the human genome era

Since the early days of the genome era, the scientific community has relied on a single 'reference' genome for each species, which is used as the basis for a wide range of genetic analyses, including studies of variation within and across species. As sequencing costs have dropped, thousands of new genomes have been sequenced, and scientists have come to realize that a single reference genome is inadequate for many purposes. By sampling a diverse set of individuals, one can begin to assemble a pan-genome: a collection of all the DNA sequences that occur in a species. Here we review efforts to create pan-genomes for a range of species, from bacteria to humans, and we further consider the computational methods that have been proposed in order to capture, interpret and compare pan-genome data. As scientists continue to survey and catalogue the genomic variation across human populations and begin to assemble a human pan-genome, these efforts will increase our power to connect variation to human diversity, disease and beyond.

Robust non-syntenic gene expression patterns in diverse maize hybrids during root development

Distantly related maize (Zea mays L.) inbred lines exhibit an exceptional degree of structural genomic diversity, which is probably unique among plants. This study systematically investigated the developmental and genotype-dependent regulation of the primary root transcriptomes of a genetically diverse panel of maize F1-hybrids and their parental inbred lines. While we observed substantial transcriptomic changes during primary root development, we demonstrated that hybrid-associated gene expression patterns, including differential, non-additive, and allele-specific transcriptome profiles, are particularly robust to these developmental fluctuations. For instance, differentially expressed genes with preferential expression in hybrids were highly conserved during development in comparison to their parental counterparts. Similarly, in hybrids a major proportion of non-additively expressed genes with expression levels between the parental values were particularly conserved during development. Importantly, in these expression patterns non-syntenic genes that evolved after the separation of the maize and sorghum lineages were systemically enriched. Furthermore, non-syntenic genes were substantially linked to the conservation of all surveyed gene expression patterns during primary root development. Among all F1-hybrids, between ~40% of the non-syntenic genes with unexpected allelic expression ratios and ~60% of the non-syntenic differentially and non-additively expressed genes were conserved and therefore robust to developmental changes. Hence, the enrichment of non-syntenic genes during primary root development might be involved in the developmental adaptation of maize roots and thus the superior performance of hybrids.

De novo Assembly of Transcriptomes From a B73 Maize Line Introgressed With a QTL for Resistance to Gray Leaf Spot Disease Reveals a Candidate Allele of a Lectin Receptor-Like Kinase

Gray leaf spot (GLS) disease in maize, caused by the fungus Cercospora zeina, is a threat to maize production globally. Understanding the molecular basis for quantitative resistance to GLS is therefore important for food security. We developed a de novo assembly pipeline to identify candidate maize resistance genes. Near-isogenic maize lines with and without a QTL for GLS resistance on chromosome 10 from inbred CML444 were produced in the inbred B73 background. The B73-QTL line showed a 20% reduction in GLS disease symptoms compared to B73 in the field (p = 0.01). B73-QTL leaf samples from this field experiment conducted under GLS disease pressure were RNA sequenced. The reads that did not map to the B73 or C. zeina genomes were expected to contain novel defense genes and were de novo assembled. A total of 141 protein-coding sequences with B73-like or plant annotations were identified from the B73-QTL plants exposed to C. zeina. To determine whether candidate gene expression was induced by C. zeina, the RNAseq reads from C. zeina-challenged and control leaves were mapped to a master assembly of all of the B73-QTL reads, and differential gene expression analysis was conducted. Combining results from both bioinformatics approaches led to the identification of a likely candidate gene, which was a novel allele of a lectin receptor-like kinase named L-RLK-CML that (i) was induced by C. zeina, (ii) was positioned in the QTL region, and (iii) had functional domains for pathogen perception and defense signal transduction. The 817AA L-RLK-CML protein had 53 amino acid differences from its 818AA counterpart in B73. A second "B73-like" allele of L-RLK was expressed at a low level in B73-QTL. Gene copy-specific RT-qPCR confirmed that the l-rlk-cml transcript was the major product induced four-fold by C. zeina. Several other expressed defense-related candidates were identified, including a wall-associated kinase, two glutathione s-transferases, a chitinase, a glucan beta-glucosidase, a plasmodesmata callose-binding protein, several other receptor-like kinases, and components of calcium signaling, vesicular trafficking, and ethylene biosynthesis. This work presents a bioinformatics protocol for gene discovery from de novo assembled transcriptomes and identifies candidate quantitative resistance genes.

A genetic characterization of Korean waxy maize (Zea mays L.) landraces having flowering time variation by RNA sequencing

Maize is the second-most produced crop in the Korean peninsula and has been continuously cultivated since the middle of the 16th century, when it was originally introduced from China. Even with this extensive cultivation history, the diversity and properties of Korean landraces have not been investigated at the nucleotide sequence level. We collected 12 landraces with various flowering times and performed RNA-seq in the early vegetative stage. The transcriptomes of 12 Korean landraces have been analyzed for their genetic variations in coding sequence and genetic relationships to other maize germplasm. The Korean landraces showed specific genetic characteristics and were closely related to a Chinese inbred line. Flowering-time related gene profiles pointed to multiple causes for the variation of flowering time within Korean landraces; the profiles revealed significant positive and negative correlations among genes, allowing us to infer possible mechanisms for flowering time variation in maize. Our results demonstrate the value of transcriptome-based genetic and gene expression profiles for information on possible breeding resources, which is particularly needed in Korean waxy landraces.

Advances in Molecular Genetics and Genomics of African Rice (Oryza glaberrima Steud)

African rice (Oryza glaberrima) has a pool of genes for resistance to diverse biotic and abiotic stresses, making it an important genetic resource for rice improvement. African rice has potential for breeding for climate resilience and adapting rice cultivation to climate change. Over the last decade, there have been tremendous technological and analytical advances in genomics that have dramatically altered the landscape of rice research. Here we review the remarkable advances in knowledge that have been witnessed in the last few years in the area of genetics and genomics of African rice. Advances in cheap DNA sequencing technologies have fuelled development of numerous genomic and transcriptomic resources. Genomics has been pivotal in elucidating the genetic architecture of important traits thereby providing a basis for unlocking important trait variation. Whole genome re-sequencing studies have provided great insights on the domestication process, though key studies continue giving conflicting conclusions and theories. However, the genomic resources of African rice appear to be under-utilized as there seems to be little evidence that these vast resources are being productively exploited for example in practical rice improvement programmes. Challenges in deploying African rice genetic resources in rice improvement and the genomics efforts made in addressing them are highlighted.

Gene Networks Underlying Cannabinoid and Terpenoid Accumulation in Cannabis

Glandular trichomes are specialized anatomical structures that accumulate secretions with important biological roles in plant-environment interactions. These secretions also have commercial uses in the flavor, fragrance, and pharmaceutical industries. The capitate-stalked glandular trichomes of Cannabis sativa (cannabis), situated on the surfaces of the bracts of the female flowers, are the primary site for the biosynthesis and storage of resins rich in cannabinoids and terpenoids. In this study, we profiled nine commercial cannabis strains with purportedly different attributes, such as taste, color, smell, and genetic origin. Glandular trichomes were isolated from each of these strains, and cell type-specific transcriptome data sets were acquired. Cannabinoids and terpenoids were quantified in flower buds. Statistical analyses indicated that these data sets enable the high-resolution differentiation of strains by providing complementary information. Integrative analyses revealed a coexpression network of genes involved in the biosynthesis of both cannabinoids and terpenoids from imported precursors. Terpene synthase genes involved in the biosynthesis of the major monoterpenes and sesquiterpenes routinely assayed by cannabis testing laboratories were identified and functionally evaluated. In addition to cloning variants of previously characterized genes, specifically CsTPS14CT [(-)-limonene synthase] and CsTPS15CT (β-myrcene synthase), we functionally evaluated genes that encode enzymes with activities not previously described in cannabis, namely CsTPS18VF and CsTPS19BL (nerolidol/linalool synthases), CsTPS16CC (germacrene B synthase), and CsTPS20CT (hedycaryol synthase). This study lays the groundwork for developing a better understanding of the complex chemistry and biochemistry underlying resin accumulation across commercial cannabis strains.

RNA-Seq analysis of soft rush (Juncus effusus): transcriptome sequencing, de novo assembly, annotation, and polymorphism identification

BACKGROUND

Juncus effusus L. (family: Juncaceae; order: Poales) is a helophytic rush growing in temperate damp or wet terrestrial habitats and is of almost cosmopolitan distribution. The species has been studied intensively with respect to its interaction with co-occurring plants as well as microbes being involved in major biogeochemical cycles. J. effusus has biotechnological value as component of Constructed Wetlands where the plant has been employed in phytoremediation of contaminated water. Its genome has not been sequenced.

RESULTS

In this study we carried out functional annotation and polymorphism analysis of de novo assembled RNA-Seq data from 18 genotypes using 249 million paired-end Illumina HiSeq reads and 2.8 million 454 Titanium reads. The assembly comprised 158,591 contigs with a mean contig length of 780 bp. The assembly was annotated using the dammit! annotation pipeline, which queries the databases OrthoDB, Pfam-A, Rfam, and runs BUSCO (Benchmarking Single-Copy Ortholog genes). In total, 111,567 contigs (70.3%) were annotated with functional descriptions, assigned gene ontology terms, and conserved protein domains, which resulted in 30,932 non-redundant gene sequences. Results of BUSCO and KEGG pathway analyses were similar for J. effusus as for the well-studied members of the Poales, Oryza sativa and Sorghum bicolor. A total of 566,433 polymorphisms were identified in transcribed regions with an average frequency of 1 polymorphism in every 171 bases.

CONCLUSIONS

The transcriptome assembly was of high quality and genome coverage was sufficient for global analyses. This annotated knowledge resource can be utilized for future gene expression analysis, genomic feature comparisons, genotyping, primer design, and functional genomics in J. effusus.

Multiple Maize Reference Genomes Impact the Identification of Variants by Genome-Wide Association Study in a Diverse Inbred Panel

Use of a single reference genome for genome-wide association studies (GWAS) limits the gene space represented to that of a single accession. This limitation can complicate identification and characterization of genes located within presence-absence variations (PAVs). In this study, we present the draft de novo genome assembly of 'PHJ89', an 'Oh43'-type inbred line of maize ( L.). From three separate reference genome assemblies ('B73', 'PH207', and PHJ89) that represent the predominant germplasm groups of maize, we generated three separate whole-seedling gene expression profiles and single nucleotide polymorphism (SNP) matrices from a panel of 942 diverse inbred lines. We identified 34,447 (B73), 39,672 (PH207), and 37,436 (PHJ89) transcripts that are not present in the respective reference genome assemblies. Genome-wide association studies were conducted in the 942 inbred panel with both the SNP and expression data values to map (SCMV) resistance. Highlighting the impact of alternative reference genomes in gene discovery, the GWAS results for SCMV resistance with expression values as a surrogate measure of PAV resulted in robust detection of the physical location of a known resistance gene when the B73 reference that contains the gene was used, but not the PH207 reference. This study provides the valuable resource of the Oh43-type PHJ89 genome assembly as well as SNP and expression data for 942 individuals generated from three different reference genomes.

Dynamic Patterns of Gene Expression Additivity and Regulatory Variation throughout Maize Development

Gene expression variation is a key component underlying phenotypic variation and heterosis. Transcriptome profiling was performed on 23 different tissues or developmental stages of two maize inbreds, B73 and Mo17, as well as their F₁ hybrid. The obtained large-scale datasets provided opportunities to monitor the developmental dynamics of differential expression, additivity for gene expression, and regulatory variation. The transcriptome can be divided into ∼30 000 genes that are expressed in at least one tissue of one inbred and an additional ∼10 000 ″silent" genes that are not expressed in any tissue of any genotype, 90% of which are non-syntenic relative to other grasses. Many (∼74%) of the expressed genes exhibit differential expression in at least one tissue. However, the majority of genes with differential expression do not exhibit consistent differential expression in different tissues. These genes often exhibit tissue-specific differential expression with equivalent expression in other tissues, and in many cases they switch the directionality of differential expression in different tissues. This suggests widespread variation for tissue-specific regulation of gene expression between the two maize inbreds B73 and Mo17. Nearly 5000 genes are expressed in only one parent in at least one tissue (single parent expression) and 97% of these genes are expressed at mid-parent levels or higher in the hybrid, providing extensive opportunities for hybrid complementation in heterosis. In general, additive expression patterns are much more common than non-additive patterns, and this trend is more pronounced for genes with strong differential expression or single parent expression. There is relatively little evidence for non-additive expression patterns that are maintained in multiple tissues. The analysis of allele-specific expression allowed classification of cis- and trans-regulatory variation. Genes with cis-regulatory variation often exhibit additive expression and tend to have more consistent regulatory variation throughout development. In contrast, genes with trans-regulatory variation are enriched for non-additive patterns and often show tissue-specific differential expression. Taken together, this study provides a deeper understanding of regulatory variation and the degree of additive gene expression throughout maize development. The dynamic nature of differential expression, additivity, and regulatory variation imply abundant variability for tissue-specific regulatory mechanisms and suggest that connections between transcriptome and phenome will require expression data from multiple tissues.

The Tannat genome: Unravelling its unique characteristics

Tannat (Vitis vinifera) is the most cultivated grapevine variety in Uruguay for the production of high quality red wines. Its berries have unusually high levels of polyphenolic compounds (anthocyanins and tannins), producing wines with intense purple colour and high antioxidant properties. Remarkably, more than 40% of its tannins are galloylated, which determines a greater antioxidant power. Technologies of massive sequencing allow the characterization of genomic variants between different clutivars. The Tannat genome was sequenced with a 134X coverage using the Illumina technology, and was annotated using transcriptomes (RNA-Seq) of different berry tissues. When comparing the genomes of Tannat with Pinot Noir PN40024 (reference genome) we found an expansion of the gene families related to the biosynthesis of polyphenols. A search base on the recently reported epicatechin galloyl transferase (ECGT) from tea leaves determined five putative genes encoding the ECGT in Tannat. Genetic analysis of one of the transcription factor that regulates the anthocyanin synthesis during berry ripening, VvMYBA1, shows the presence of Gret1 retrotransposon in one of the VvMYBA1 alleles in the Tannat clones analysed. This work makes original contributions about the molecular bases of the biosynthesis of anthocyanins and tannins during the development of the flagship grape of Uruguay.

Whole-Genome Resequencing and Pan-Transcriptome Reconstruction Highlight the Impact of Genomic Structural Variation on Secondary Metabolite Gene Clusters in the Grapevine Esca Pathogen Phaeoacremonium minimum

The Ascomycete fungus Phaeoacremonium minimum is one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence among Pm. minimum isolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference of Pm. minimum comprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included the Pm. minimum pan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.

Genomic prediction applied to high-biomass sorghum for bioenergy production

The increasing cost of energy and finite oil and gas reserves have created a need to develop alternative fuels from renewable sources. Due to its abiotic stress tolerance and annual cultivation, high-biomass sorghum (Sorghum bicolor L. Moench) shows potential as a bioenergy crop. Genomic selection is a useful tool for accelerating genetic gains and could restructure plant breeding programs by enabling early selection and reducing breeding cycle duration. This work aimed at predicting breeding values via genomic selection models for 200 sorghum genotypes comprising landrace accessions and breeding lines from biomass and saccharine groups. These genotypes were divided into two sub-panels, according to breeding purpose. We evaluated the following phenotypic biomass traits: days to flowering, plant height, fresh and dry matter yield, and fiber, cellulose, hemicellulose, and lignin proportions. Genotyping by sequencing yielded more than 258,000 single-nucleotide polymorphism markers, which revealed population structure between subpanels. We then fitted and compared genomic selection models BayesA, BayesB, BayesCπ, BayesLasso, Bayes Ridge Regression and random regression best linear unbiased predictor. The resulting predictive abilities varied little between the different models, but substantially between traits. Different scenarios of prediction showed the potential of using genomic selection results between sub-panels and years, although the genotype by environment interaction negatively affected accuracies. Functional enrichment analyses performed with the marker-predicted effects suggested several interesting associations, with potential for revealing biological processes relevant to the studied quantitative traits. This work shows that genomic selection can be successfully applied in biomass sorghum breeding programs.

Genomic and epigenomic immunity in common bean: the unusual features of NB-LRR gene family

In plants, a key class of genes comprising most of disease resistance (R) genes encodes Nucleotide-binding leucine-rich repeat (NL) proteins. Access to common bean (Phaseolus vulgaris) genome sequence provides unparalleled insight into the organization and evolution of this large gene family (∼400 NL) in this important crop. As observed in other plant species, most common bean NL are organized in cluster of genes. However, a particularity of common bean is that these clusters are often located in subtelomeric regions close to terminal knobs containing the satellite DNA khipu. Phylogenetically related NL are spread between different chromosome ends, suggesting frequent exchanges between non-homologous chromosomes. NL peculiar location, in proximity to heterochromatic regions, led us to study their DNA methylation status using a whole-genome cytosine methylation map. In common bean, NL genes displayed an unusual body methylation pattern since half of them are methylated in the three contexts, reminiscent of the DNA methylation pattern of repeated sequences. Moreover, 90 NL were also abundantly targeted by 24 nt siRNA, with 90% corresponding to methylated NL genes. This suggests the existence of a transcriptional gene silencing mechanism of NL through the RdDM (RNA-directed DNA methylation) pathway in common bean that has not been described in other plant species.

A transcriptome-SNP-derived linkage map of Apios americana (potato bean) provides insights about genome re-organization and synteny conservation in the phaseoloid legumes

KEY MESSAGE

We report a linkage map for Apios americana and describe synteny with selected warm-season legumes. A translocation event in common bean and soybean is confirmed against Apios and Vigna species. Apios (Apios americana; "apios"), a tuberous perennial legume in the Phaseoleae tribe, was widely used as a food by Native Americans. Work in the last 40 years has led to several improved breeding lines. Aspects of the pollination biology (complex floral structure and tripping mechanism) have made controlled crosses difficult, and the previous reports indicated that the plant is likely primarily an outcrosser. We used a pseudo-testcross strategy to construct a genetic map specific to the maternal parent. The map was built using single-nucleotide polymorphism markers identified by comparing the expressed sequences of individuals in the mapping population against a de novo maternal reference transcriptome assembly. The apios map consists of 11 linkage groups and 1121 recombinationally distinct loci, covering ~ 938.6 cM. By sequencing the transcriptomes of all potential pollen parents, we were able to identify the probable pollen donors and to discover new aspects of the pollination biology in apios. No selfing was observed, but multiple pollen parents were seen within individual pods. Comparisons with genome sequences in other species in the Phaseoleae showed extended synteny for most apios linkage groups. This synteny supports the robustness of the map, and also sheds light on the history of the Phaseoleae, as apios is relatively early diverging in this tribe. We detected a translocation event that separates apios and two Vigna species from Phaseolus vulgaris and Glycine max. This apios mapping work provides a general protocol for sequencing-based construction of high-density linkage maps in outcrossing species with heterogeneous pollen parents.

Circular RNAs mediated by transposons are associated with transcriptomic and phenotypic variation in maize

Circular RNAs (circRNAs) are covalently closed RNA molecules. Recent studies have shown that circRNAs can arise from the transcripts of transposons. Given the prevalence of transposons in the maize genome and dramatic genomic variation driven by transposons, we hypothesize that transposons in maize may be involved in the formation of circRNAs and further modulate phenotypic variation. We performed circRNA-Seq on B73 seedling leaves and uncovered 2804 high-confidence maize circRNAs, which show distinct genomic features. Comprehensive analyses demonstrated that sequences related to LINE1-like elements (LLEs) and their Reverse Complementary Pairs (LLERCPs) are significantly enriched in the flanking regions of circRNAs. Interestingly, as the number of LLERCPs increase, the accumulation of circRNAs varies, whereas that of linear transcripts decreases. Furthermore, genes with LLERCP-mediated circRNAs are enriched among loci that are associated with phenotypic variation. These results suggest that circRNAs are likely to be involved in the modulation of phenotypic variation by LLERCPs. Further, we showed that the presence/absence variation of LLERCPs was associated with expression variation of circRNA-circ1690 and was related to ear height, potentially through the interplay between circRNAs and functional linear transcripts. Our first study of maize circRNAs uncovers a potential new way for transposons to modulate transcriptomic and phenotypic variations.

Transcriptomic characterization and potential marker development of contrasting sugarcane cultivars

Sugarcane (Saccharum officinarum L.) is an important crop for sugar production and bioenergy worldwide. In this study, we performed transcriptome sequencing for six contrasting sugarcane genotypes involved in leaf abscission, tolerance to pokkah boeng disease and drought stress. More than 465 million high-quality reads were generated, which were de novo assembled into 93,115 unigenes. Based on a similarity search, 43,526 (46.74%) unigenes were annotated against at least one of the public databases. Functional classification analyses showed that these unigenes are involved in a wide range of metabolic pathways. Comparative transcriptome analysis revealed that many unigenes involved in response to abscisic acid and ethylene were up-regulated in the easy leaf abscission genotype, and unigenes associated with response to jasmonic acid and salicylic acid were up-regulated in response to the pokkah boeng disease in the tolerance genotype. Moreover, unigenes related to peroxidase, antioxidant activity and signal transduction were up-regulated in response to drought stress in the tolerant genotype. Finally, we identified a number of putative markers, including 8,630 simple sequence repeats (SSRs) and 442,152 single-nucleotide polymorphisms (SNPs). Our data will be important resources for future gene discovery, molecular marker development, and genome studies in sugarcane.

The limited role of differential fractionation in genome content variation and function in maize (Zea mays L.) inbred lines

Maize is a diverse paleotetraploid species with considerable presence/absence variation and copy number variation. One mechanism through which presence/absence variation can arise is differential fractionation. Fractionation refers to the loss of duplicate gene pairs from one of the maize subgenomes during diploidization. Differential fractionation refers to non-shared gene loss events between individuals following a whole-genome duplication event. We investigated the prevalence of presence/absence variation resulting from differential fractionation in the syntenic portion of the genome using two whole-genome de novo assemblies of the inbred lines B73 and PH207. Between these two genomes, syntenic genes were highly conserved with less than 1% of syntenic genes being subject to differential fractionation. The few variably fractionated syntenic genes that were identified are unlikely to contribute to functional phenotypic variation, as there is a significant depletion of these genes in annotated gene sets. In further comparisons of 60 diverse inbred lines, non-syntenic genes were six times more likely to be variable than syntenic genes, suggesting that comparisons among additional genome assemblies are not likely to result in the discovery of large-scale presence/absence variation among syntenic genes.

A High-Density Genetic Map of an Allohexaploid Brassica Doubled Haploid Population Reveals Quantitative Trait Loci for Pollen Viability and Fertility

A doubled haploid (DH) mapping population was obtained from microspore culture of an allohexaploid F₁ from the cross between two recently-synthesized allohexaploid Brassica lines. We used single nucleotide polymorphism (SNP) genetic variation based on restriction-site associated DNA (RAD) sequencing to construct a high density genetic linkage map of the population. RAD libraries were constructed from the genomic DNA of both parents and 146 DH progenies. A total of 2.87 G reads with an average sequencing depth of 2.59 × were obtained in the parents and of 1.41 × in the progeny. A total of 290,422 SNPs were identified from clustering of RAD reads, from which we developed 7,950 high quality SNP markers that segregated normally (1:1) in the population. The linkage map contained all 27 chromosomes from the parental A, B and C genomes with a total genetic distance of 5725.19 cM and an average of 0.75 cM between adjacent markers. Genetic distance on non-integrated linkage groups was 1534.23 cM, or 21% of total genetic distance. Out of 146 DH progenies, 91 had a complete set of 27 chromosomes as expected of a hexaploid species, and 21 out of 27 chromosomes showed high collinearity between the physical and linkage maps. The loss of chromosome(s) or chromosome segment(s) in the DH population was associated with a reduction in pollen viability. Twenty-five additive QTL were associated with pollen viability and fertility-related traits (seed number, seed yield, pod length, plant height, 1000-seed weight). In addition, 44 intra-genomic and 18 inter-genomic epistatic QTL pairs were detected for 4 phenotypic traits. This provides confidence that the DH population may be selected for improved pollen viability and fertility in a future allohexaploid Brassica species.

Brassinosteroid and gibberellin control of seedling traits in maize (Zea mays L.)

In this study, we established two doubled haploid (DH) libraries with a total of 207 DH lines. We applied BR and GA inhibitors to all DH lines at seedling stage and measured seedling BR and GA inhibitor responses. Moreover, we evaluated field traits for each DH line (untreated). We conducted genome-wide association studies (GWAS) with 62,049 genome wide SNPs to explore the genetic control of seedling traits by BR and GA. In addition, we correlate seedling stage hormone inhibitor response with field traits. Large variation for BR and GA inhibitor response and field traits was observed across these DH lines. Seedling stage BR and GA inhibitor response was significantly correlate with yield and flowering time. Using three different GWAS approaches to balance false positive/negatives, multiple SNPs were discovered to be significantly associated with BR/GA inhibitor responses with some localized within gene models. SNPs from gene model GRMZM2G013391 were associated with GA inhibitor response across all three GWAS models. This gene is expressed in roots and shoots and was shown to regulate GA signaling. These results show that BRs and GAs have a great impact for controlling seedling growth. Gene models from GWAS results could be targets for seeding traits improvement.

The Genetics of Plant Metabolism

Plant metabolic studies have traditionally focused on the role and regulation of the enzymes catalyzing key reactions within specific pathways. Within the past 20 years, reverse genetic approaches have allowed direct determination of the effects of the deficiency, or surplus, of a given protein on the biochemistry of a plant. In parallel, top-down approaches have also been taken, which rely on screening broad, natural genetic diversity for metabolic diversity. Here, we compare and contrast the various strategies that have been adopted to enhance our understanding of the natural diversity of metabolism. We also detail how these approaches have enhanced our understanding of both specific and global aspects of the genetic regulation of metabolism. Finally, we discuss how such approaches are providing important insights into the evolution of plant secondary metabolism.

Analysis of the first Taraxacum kok-saghyz transcriptome reveals potential rubber yield related SNPs

Taraxacum kok-saghyz (TK) is a potential alternative crop for natural rubber (NR) production, due to its high molecular weight rubber, short breeding cycle, and diverse environmental adaptation. However, improvements in rubber yield and agronomically relevant traits are still required before it can become a commercially-viable crop. An RNA-Seq based transcriptome was developed from a pool of roots from genotypes with high and low rubber yield. A total of 55,532 transcripts with lengths over 200 bp were de novo assembled. As many as 472 transcripts were significantly homologous to 49 out of 50 known plant putative rubber biosynthesis related genes. 158 transcripts were significantly differentially expressed between high rubber and low rubber genotypes. 21,036 SNPs were different in high and low rubber TK genotypes. Among these, 50 SNPs were found within 39 transcripts highly homologous to 49 publically-searched rubber biosynthesis related genes. 117 SNPs were located within 36 of the differentially expressed gene sequences. This comprehensive TK transcriptomic reference, and large set of SNPs including putative exonic markers associated with rubber related gene homologues and differentially expressed genes, provides a solid foundation for further genetic dissection of rubber related traits, comparative genomics and marker-assisted selection for the breeding of TK.

Enhancing genetic gain in the era of molecular breeding

As one of the important concepts in conventional quantitative genetics and breeding, genetic gain can be defined as the amount of increase in performance that is achieved annually through artificial selection. To develop pro ducts that meet the increasing demand of mankind, especially for food and feed, in addition to various industrial uses, breeders are challenged to enhance the potential of genetic gain continuously, at ever higher rates, while they close the gaps that remain between the yield potential in breeders' demonstration trials and the actual yield in farmers' fields. Factors affecting genetic gain include genetic variation available in breeding materials, heritability for traits of interest, selection intensity, and the time required to complete a breeding cycle. Genetic gain can be improved through enhancing the potential and closing the gaps, which has been evolving and complemented with modern breeding techniques and platforms, mainly driven by molecular and genomic tools, combined with improved agronomic practice. Several key strategies are reviewed in this article. Favorable genetic variation can be unlocked and created through molecular and genomic approaches including mutation, gene mapping and discovery, and transgene and genome editing. Estimation of heritability can be improved by refining field experiments through well-controlled and precisely assayed environmental factors or envirotyping, particularly for understanding and controlling spatial heterogeneity at the field level. Selection intensity can be significantly heightened through improvements in the scale and precision of genotyping and phenotyping. The breeding cycle time can be shortened by accelerating breeding procedures through integrated breeding approaches such as marker-assisted selection and doubled haploid development. All the strategies can be integrated with other widely used conventional approaches in breeding programs to enhance genetic gain. More transdisciplinary approaches, team breeding, will be required to address the challenge of maintaining a plentiful and safe food supply for future generations. New opportunities for enhancing genetic gain, a high efficiency breeding pipeline, and broad-sense genetic gain are also discussed prospectively.

Combining ability and testcross performance of drought-tolerant maize inbred lines under stress and non-stress environments in Kenya

Drought and poor soil fertility are among the major abiotic stresses affecting maize productivity in sub-Saharan Africa. Maize breeding efforts at the International Maize and Wheat Improvement Center (CIMMYT) have focused on incorporating drought stress tolerance and nitrogen-use efficiency (NUE) into tropical maize germplasm. The objectives of this study were to estimate the general combining ability (GCA) and specific combining ability (SCA) of selected maize inbred lines under drought stress (DS), low-nitrogen (LN) and optimum moisture and nitrogen (optimum) conditions, and to assess the yield potential and stability of experimental hybrids under these management conditions. Forty-nine experimental three-way cross hybrids, generated from a 7 × 7 line by tester crosses, and six commercial checks were evaluated across 11 optimum, DS and LN sites in Kenya in 2014 using an alpha lattice design with two replicates per entry at each site. DS reduced both grain yield (GY) and plant height (PH), while anthesis-silking interval (ASI) increased under both DS and LN. Hybrids 'L4/T2' and 'L4/T1' were found to be superior and stable, while inbreds 'L4' and 'L6' were good combiners for GY and other secondary traits across sites. Additive variance played a greater role for most traits under the three management conditions, suggesting that further progress in the improvement of these traits should be possible. GY under optimum conditions was positively correlated with GY under both DS and LN conditions, but GY under DS and LN was not correlated. Our results suggest the feasibility for simultaneous improvement in grain yield performance of genotypes under optimum, DS and LN conditions.

Genotyping-by-sequencing approaches to characterize crop genomes: choosing the right tool for the right application

In the last decade, the revolution in sequencing technologies has deeply impacted crop genotyping practice. New methods allowing rapid, high-throughput genotyping of entire crop populations have proliferated and opened the door to wider use of molecular tools in plant breeding. These new genotyping-by-sequencing (GBS) methods include over a dozen reduced-representation sequencing (RRS) approaches and at least four whole-genome resequencing (WGR) approaches. The diversity of methods available, each often producing different types of data at different cost, can make selection of the best-suited method seem a daunting task. We review the most common genotyping methods used today and compare their suitability for linkage mapping, genomewide association studies (GWAS), marker-assisted and genomic selection and genome assembly and improvement in crops with various genome sizes and complexity. Furthermore, we give an outline of bioinformatics tools for analysis of genotyping data. WGR is well suited to genotyping biparental cross populations with complex, small- to moderate-sized genomes and provides the lowest cost per marker data point. RRS approaches differ in their suitability for various tasks, but demonstrate similar costs per marker data point. These approaches are generally better suited for de novo applications and more cost-effective when genotyping populations with large genomes or high heterozygosity. We expect that although RRS approaches will remain the most cost-effective for some time, WGR will become more widespread for crop genotyping as sequencing costs continue to decrease.

Stability of Single-Parent Gene Expression Complementation in Maize Hybrids upon Water Deficit Stress

Heterosis is the superior performance of F1 hybrids compared with their homozygous, genetically distinct parents. In this study, we monitored the transcriptomic divergence of the maize (Zea mays) inbred lines B73 and Mo17 and their reciprocal F1 hybrid progeny in primary roots under control and water deficit conditions simulated by polyethylene glycol treatment. Single-parent expression (SPE) of genes is an extreme instance of gene expression complementation, in which genes are active in only one of two parents but are expressed in both reciprocal hybrids. In this study, 1,997 genes only expressed in B73 and 2,024 genes only expressed in Mo17 displayed SPE complementation under control and water deficit conditions. As a consequence, the number of active genes in hybrids exceeded the number of active genes in the parental inbred lines significantly independent of treatment. SPE patterns were substantially more stable to expression changes by water deficit treatment than other genotype-specific expression profiles. While, on average, 75% of all SPE patterns were not altered in response to polyethylene glycol treatment, only 17% of the remaining genotype-specific expression patterns were not changed by water deficit. Nonsyntenic genes that lack syntenic orthologs in other grass species, and thus evolved late in the grass lineage, were significantly overrepresented among SPE genes. Hence, the significant overrepresentation of nonsyntenic genes among SPE patterns and their stability under water limitation might suggest a function of these genes during the early developmental manifestation of heterosis under fluctuating environmental conditions in hybrid progeny of the inbred lines B73 and Mo17.

Transcriptome Analysis of Flowering Time Genes under Drought Stress in Maize Leaves

Flowering time is an important factor determining yield and seed quality in maize. A change in flowering time is a strategy used to survive abiotic stresses. Among abiotic stresses, drought can increase anthesis-silking intervals (ASI), resulting in negative effects on maize yield. We have analyzed the correlation between flowering time and drought stress using RNA-seq and bioinformatics tools. Our results identified a total of 619 genes and 126 transcripts whose expression was altered by drought stress in the maize B73 leaves under short-day condition. Among drought responsive genes, we also identified 20 genes involved in flowering times. Gene Ontology (GO) enrichment analysis was used to predict the functions of the drought-responsive genes and transcripts. GO categories related to flowering time included reproduction, flower development, pollen-pistil interaction, and post-embryonic development. Transcript levels of several genes that have previously been shown to affect flowering time, such as PRR37, transcription factor HY5, and CONSTANS, were significantly altered by drought conditions. Furthermore, we also identified several drought-responsive transcripts containing C₂H₂ zinc finger, CCCH, and NAC domains, which are frequently involved in transcriptional regulation and may thus have potential to alter gene expression programs to change maize flowering time. Overall, our results provide a genome-wide analysis of differentially expressed genes (DEGs), novel transcripts, and isoform variants expressed during the reproductive stage of maize plants subjected to drought stress and short-day condition. Further characterization of the drought-responsive transcripts identified in this study has the potential to advance our understanding of the mechanisms that regulate flowering time under drought stress.

Analysis of Plant Pan-Genomes and Transcriptomes with GET_HOMOLOGUES-EST, a Clustering Solution for Sequences of the Same Species

The pan-genome of a species is defined as the union of all the genes and non-coding sequences found in all its individuals. However, constructing a pan-genome for plants with large genomes is daunting both in sequencing cost and the scale of the required computational analysis. A more affordable alternative is to focus on the genic repertoire by using transcriptomic data. Here, the software GET_HOMOLOGUES-EST was benchmarked with genomic and RNA-seq data of 19 Arabidopsis thaliana ecotypes and then applied to the analysis of transcripts from 16 Hordeum vulgare genotypes. The goal was to sample their pan-genomes and classify sequences as core, if detected in all accessions, or accessory, when absent in some of them. The resulting sequence clusters were used to simulate pan-genome growth, and to compile Average Nucleotide Identity matrices that summarize intra-species variation. Although transcripts were found to under-estimate pan-genome size by at least 10%, we concluded that clusters of expressed sequences can recapitulate phylogeny and reproduce two properties observed in A. thaliana gene models: accessory loci show lower expression and higher non-synonymous substitution rates than core genes. Finally, accessory sequences were observed to preferentially encode transposon components in both species, plus disease resistance genes in cultivated barleys, and a variety of protein domains from other families that appear frequently associated with presence/absence variation in the literature. These results demonstrate that pan-genome analyses are useful to explore germplasm diversity.

Pan- and core- network analysis of co-expression genes in a model plant

Genome-wide gene expression experiments have been performed using the model plant Arabidopsis during the last decade. Some studies involved construction of coexpression networks, a popular technique used to identify groups of co-regulated genes, to infer unknown gene functions. One approach is to construct a single coexpression network by combining multiple expression datasets generated in different labs. We advocate a complementary approach in which we construct a large collection of 134 coexpression networks based on expression datasets reported in individual publications. To this end we reanalyzed public expression data. To describe this collection of networks we introduced concepts of 'pan-network' and 'core-network' representing union and intersection between a sizeable fractions of individual networks, respectively. We showed that these two types of networks are different both in terms of their topology and biological function of interacting genes. For example, the modules of the pan-network are enriched in regulatory and signaling functions, while the modules of the core-network tend to include components of large macromolecular complexes such as ribosomes and photosynthetic machinery. Our analysis is aimed to help the plant research community to better explore the information contained within the existing vast collection of gene expression data in Arabidopsis.

Transcriptome analysis in switchgrass discloses ecotype difference in photosynthetic efficiency

BACKGROUND

Switchgrass, a warm-season perennial grass studied as a potential dedicated biofuel feedstock, is classified into two main taxa - lowland and upland ecotypes - that differ in morphology and habitat of adaptation. But there is limited information on their inherent molecular variations.

RESULTS

Transcriptome analysis by RNA-sequencing (RNA-Seq) was conducted for lowland and upland ecotypes to document their gene expression variations. Mapping of transcriptome to the reference genome (Panicum virgatum v1.1) revealed that the lowland and upland ecotypes differ substantially in sets of genes transcribed as well as levels of expression. Differential gene expression analysis exhibited that transcripts related to photosynthesis efficiency and development and photosystem reaction center subunits were upregulated in lowlands compared to upland genotype. On the other hand, catalase isozymes, helix-loop-helix, late embryogenesis abundant group I, photosulfokinases, and S-adenosyl methionine synthase gene transcripts were upregulated in the upland compared to the lowlands. At ≥100x coverage and ≥5% minor allele frequency, a total of 25,894 and 16,979 single nucleotide polymorphism (SNP) markers were discovered for VS16 (upland ecotype) and K5 (lowland ecotype) against the reference genome. The allele combination of the SNPs revealed that the transition mutations are more prevalent than the transversion mutations.

CONCLUSIONS

The gene ontology (GO) analysis of the transcriptome indicated lowland ecotype had significantly higher representation for cellular components associated with photosynthesis machinery controlling carbon fixation. In addition, using the transcriptome data, SNP markers were detected, which were distributed throughout the genome. The differentially expressed genes and SNP markers detected in this study would be useful resources for traits mapping and gene transfer across ecotypes in switchgrass breeding for increased biomass yield for biofuel conversion.

NtPHYB1K326, a homologous gene of Arabidopsis PHYB, positively regulates the content of phenolic compounds in tobacco

Polyphenols are important secondary metabolites and bioactive compounds in plants. Light is a vital abiotic factor that greatly impacts the content of polyphenols in plants. In spite of their importance the mechanism of polyphenol regulation still remains unknown in tobacco. A phytochrome B homolog, NtPHYB1K326, was isolated from Nicotiana tabacum cv. K326 to investigate the role of light receptors in the regulation of polyphenol metabolism in tobacco leaves. Furthermore, role of NtPHYB1K326 in polyphenol metabolism was analyzed by over-expression and RNAi-silencing approaches. Consistent and complemented results indicated involvement of NtPHYB1K326 in the regulation of polyphenol metabolism in tobacco leaves. Moreover, high levels of NtPHYB1K326 transcripts favor the accumulation of chlorogenic acid and its isomers, the key polyphenol component in tobacco leaves. Transcriptome analysis was also carried out for exploring the regulation mechanism of NtPHYB1K326 in the polyphenol metabolism. Compared with WT, 1665 and 1421 differentially-expressed genes were found in NtPHYB1K326-GFP and NtPHYB1K326-RNAi transgenic lines, respectively. Among these, about 30 genes were related to phenylpropanoid pathway, which is predominantly involved in synthesis of polyphenols. Further evidences from quantitative RT-PCR confirmed that NtPHYB1K326 may control phenylpropanoid pathway through regulating the transcription of PAL4 (phenylalanine ammonialyase 4), 4CL1 (4-coumarate:coenzyme A ligase 1) and COMT (caffeic acid 3-O-methyltransferase) genes.

Draft Assembly of Elite Inbred Line PH207 Provides Insights into Genomic and Transcriptome Diversity in Maize

Intense artificial selection over the last 100 years has produced elite maize (Zea mays) inbred lines that combine to produce high-yielding hybrids. To further our understanding of how genome and transcriptome variation contribute to the production of high-yielding hybrids, we generated a draft genome assembly of the inbred line PH207 to complement and compare with the existing B73 reference sequence. B73 is a founder of the Stiff Stalk germplasm pool, while PH207 is a founder of Iodent germplasm, both of which have contributed substantially to the production of temperate commercial maize and are combined to make heterotic hybrids. Comparison of these two assemblies revealed over 2500 genes present in only one of the two genotypes and 136 gene families that have undergone extensive expansion or contraction. Transcriptome profiling revealed extensive expression variation, with as many as 10,564 differentially expressed transcripts and 7128 transcripts expressed in only one of the two genotypes in a single tissue. Genotype-specific genes were more likely to have tissue/condition-specific expression and lower transcript abundance. The availability of a high-quality genome assembly for the elite maize inbred PH207 expands our knowledge of the breadth of natural genome and transcriptome variation in elite maize inbred lines across heterotic pools.