Spatial transcriptomics uncover sucrose post-phloem transport during maize kernel development.

Maize kernels are complex biological systems composed of three genetic sources, namely maternal tissues, progeny embryos, and progeny endosperms. The lack of gene expression profiles with spatial information has limited the understanding of the specific functions of each cell population, and hindered the exploration of superior genes in kernels. In our study, we conduct microscopic sectioning and spatial transcriptomics analysis during the grain filling stage of maize kernels. This enables us to visualize the expression patterns of all genes through electronical RNA in situ hybridization, and identify 11 cell populations and 332 molecular marker genes. Furthermore, we systematically elucidate the spatial storage mechanisms of the three major substances in maize kernels: starch, protein, and oil. These findings provide valuable insights into the functional genes that control agronomic traits in maize kernels.

New insight into the genetic basis of oil content based on noninvasive three-dimensional phenotyping and tissue-specific transcriptome in Brassica napus.

BACKGROUND: Increasing seed oil content is the most important breeding goal in Brassica napus, and phenotyping is crucial to dissect its genetic basis in crops. To date, QTL mapping for oil content has been based on whole seeds, and the lipid distribution is far from uniform in different tissues of seeds in B. napus. In this case, the phenotype based on whole seeds was unable to sufficiently reveal the complex genetic characteristics of seed oil content. RESULTS: Here, the three-dimensional (3D) distribution of lipid was determined for B. napus seeds by magnetic resonance imaging (MRI) and 3D quantitative analysis, and ten novel oil content-related traits were obtained by subdividing the seeds. Based on a high-density genetic linkage map, 35 QTLs were identified for 4 tissues, the outer cotyledon (OC), inner cotyledon (IC), radicle (R) and seed coat (SC), which explained up to 13.76% of the phenotypic variation. Notably, 14 tissue-specific QTLs were reported for the first time, 7 of which were novel. Moreover, haplotype analysis showed that the favorable alleles for different seed tissues exhibited cumulative effects on oil content. Furthermore, tissue-specific transcriptomes revealed that more active energy and pyruvate metabolism influenced carbon flow in the IC, OC and R than in the SC at the early and middle seed development stages, thus affecting the distribution difference in oil content. Combining tissue-specific QTL mapping and transcriptomics, 86 important candidate genes associated with lipid metabolism were identified that underlie 19 unique QTLs, including the fatty acid synthesis rate-limiting enzyme-related gene CAC2, in the QTLs for OC and IC. CONCLUSIONS: The present study provides further insight into the genetic basis of seed oil content at the tissue-specific level.

A genome-wide association study identifies a transporter for zinc uploading to maize kernels.

The Zn content in cereal seeds is an important trait for crop production as well as for human health. However, little is known about how Zn is loaded to plant seeds. Here, through a genome-wide association study (GWAS), we identify the Zn-NA (nicotianamine) transporter gene ZmYSL2 that is responsible for loading Zn to maize kernels. High promoter sequence variation in ZmYSL2 most likely drives the natural variation in Zn concentrations in maize kernels. ZmYSL2 is specifically localized on the plasma membrane facing the maternal tissue of the basal endosperm transfer cell layer (BETL) and functions in loading Zn-NA into the BETL. Overexpression of ZmYSL2 increases the Zn concentration in the kernels by 31.6%, which achieves the goal of Zn biofortification of maize. These findings resolve the mystery underlying the loading of Zn into plant seeds, providing an efficient strategy for breeding or engineering maize varieties with enriched Zn nutrition.

Transcriptome analysis of thymic tissues from Chinese Partridge Shank chickens with or without Newcastle disease virus LaSota vaccine injection via high-throughput RNA sequencing.

The LaSota strain of Newcastle disease virus (NDV) is a commonly used vaccine to control Newcastle disease. However, improper immunization is a common reason for vaccine failure. Hence, it is imperative to thoroughly explore innate immunity-related molecular regulatory responses to the LaSota vaccine. In this text, 140 long non-coding RNAs (lncRNAs), 8 microRNAs (miRNAs), and 1514 mRNAs were identified to be differentially expressed by RNA sequencing analysis in the thymic tissues of Chinese Partridge Shank chickens after LaSota vaccine inoculation. Moreover, 70 dysregulated genes related to innate immunity were identified based on GO, Reactome pathway, and InnateDB annotations and differential expression analysis. Additionally, dysregulated lncRNAs and innate immunity-related mRNAs that could interact with dysregulated miRNAs were identified based on bioinformatics prediction analysis via the miRanda software and differential expression analysis. Among these transcripts, expression patterns of five lncRNAs, seven miRNAs, and six mRNAs were further examined by RT-qPCR assay. Both RNA-seq and RT-qPCR outcomes showed that 10 transcripts (MSTRG.22689.1, ENSGALT00000065826, ENSGALT00000059336, ENSGALT00000060887, gga-miR-6575-5p, gga-miR-6631-5p, gga-miR-1727, paraoxonase 2 (PON2), mitogen-activated protein kinase 10, and cystic fibrosis transmembrane conductance regulator (CFTR) were highly expressed, and 4 transcripts (MSTRG.188121.10, gga-miR-6655-5p, gga-miR-6548-3p, and matrix metallopeptidase 9 (MMP9) were low expressed after NDV infection. Additionally, two potential competing endogenous RNA networks (ENSGALT00000060887/gga-miR-6575-5p/PON2 or MSTRG.188121.10/gga-miR-6631-5p/MMP9) and some co-expression axes (ENSGALT00000065826/gga-miR-6631-5p, MSTRG.188121.10/gga-miR-6575-5p, MSTRG.188121.10/CFTR, ENSGALT00000060887/MMP9) were identified based on RT-qPCR and co-expression analyses. In conclusion, we identified multiple dysregulated lncRNAs, miRNAs, and mRNAs after LaSota infection and some potential regulatory networks for these dysregulated transcripts.

ABA-induced phosphorylation of basic leucine zipper 29, ABSCISIC ACID INSENSITIVE 19, and Opaque2 by SnRK2.2 enhances gene transactivation for endosperm filling in maize.

Opaque2 (O2) functions as a central regulator of the synthesis of starch and storage proteins and the O2 gene is transcriptionally regulated by a hub coordinator of seed development and grain filling, ABSCISIC ACID INSENSITIVE 19 (ZmABI19), in maize (Zea mays). Here, we identified a second hub coordinator, basic Leucine Zipper 29 (ZmbZIP29) that interacts with ZmABI19 to regulate O2 expression. Like zmabi19, zmbzip29 mutations resulted in a dramatic decrease of transcript and protein levels of O2 and thus a significant reduction of starch and storage proteins. zmbzip29 seeds developed slower and had a smaller size at maturity than those of the wild type. The zmbzip29;zmabi19 double mutant displayed more severe seed phenotypes and a greater reduction of storage reserves compared to the single mutants, whereas overexpression of the two transcription factors enhanced O2 expression, storage-reserve accumulation, and kernel weight. ZmbZIP29, ZmABI19, and O2 expression was induced by abscisic acid (ABA). With ABA treatment, ZmbZIP29 and ZmABI19 synergistically transactivated the O2 promoter. Through liquid chromatography tandem-mass spectrometry analysis, we established that the residues threonine(T) 57 in ZmABI19, T75 in ZmbZIP29, and T387 in O2 were phosphorylated, and that SnRK2.2 was responsible for the phosphorylation. The ABA-induced phosphorylation at these sites was essential for maximum transactivation of downstream target genes for endosperm filling in maize.

A major yellow-seed QTL on chromosome A09 significantly increases the oil content and reduces the fiber content of seed in Brassica napus.

KEY MESSAGE: A major yellow-seed QTL on chromosome A09 significantly increases the oil content and reduces the fiber content of seed in Brassica napus. The yellow-seed trait (YST) has always been a main breeding objective for rapeseed because yellow-seeded B. napus generally contains higher oil contents, fewer pigments and polyphenols and lower fiber content than black-seeded B. napus, although the mechanism controlling this correlation remains unclear. In this study, QTL mapping was implemented for YST based on a KN double haploid population derived from the hybridization of yellow-seeded B. napus N53-2 with a high oil content and black-seeded Ken-C8 with a relatively low oil content. Ten QTLs were identified, including four stable QTLs that could be detected in multiple environments. A major QTL, cqSC-A09, on chromosome A09 was identified by both QTL mapping and BSR-Seq technology, and explained more than 41% of the phenotypic variance. The major QTL cqSC-A09 for YST not only controls the seed color but also affects the oil and fiber contents in seeds. More importantly, the advantageous allele could increase the oil content and reduce the pigment and fiber content at the same time. This is the first QTL reported to control seed color, oil content and fiber content simultaneously with a large effect and has great application value for breeding high oil varieties with high seed quality. Important candidate genes, including BnaA09. JAZ1, BnaA09. GH3.3 and BnaA09. LOX3, were identified for cqSC-A09 by combining sequence variation annotation, expression differences and an interaction network, which lays a foundation for further cloning and breeding applications in the future.

Thymic transcriptome analysis after Newcastle disease virus inoculation in chickens and the influence of host small RNAs on NDV replication.

Newcastle disease (ND), caused by virulent Newcastle disease virus (NDV), is a contagious viral disease affecting various birds and poultry worldwide. In this project, differentially expressed (DE) circRNAs, miRNAs and mRNAs were identified by high-throughput RNA sequencing (RNA-Seq) in chicken thymus at 24, 48, 72 or 96 h post LaSota NDV vaccine injection versus pre-inoculation group. The vital terms or pathways enriched by vaccine-influenced genes were tested through KEGG and GO analysis. DE genes implicated in innate immunity were preliminarily screened out through GO, InnateDB and Reactome Pathway databases. The interaction networks of DE innate immune genes were established by STRING website. Considering the high expression of gga-miR-6631-5p across all the four time points, DE circRNAs or mRNAs with the possibility to bind to gga-miR-6631-5p were screened out. Among DE genes that had the probability to interact with gga-miR-6631-5p, 7 genes were found to be related to innate immunity. Furthermore, gga-miR-6631-5p promoted LaSota NDV replication by targeting insulin induced gene 1 (INSIG1) in DF-1 chicken fibroblast cells. Taken together, our data provided the comprehensive information about molecular responses to NDV LaSota vaccine in Chinese Partridge Shank Chickens and elucidated the vital roles of gga-miR-6631-5p/INSIG1 axis in LaSota NDV replication.

The B3 domain-containing transcription factor ZmABI19 coordinates expression of key factors required for maize seed development and grain filling.

Grain filling in maize (Zea mays) is regulated by a group of spatiotemporally synchronized transcription factors (TFs), but the factors that coordinate their expression remain unknown. We used the promoter of the grain filling-specific TF gene Opaque2 (O2) to screen upstream regulatory factors and identified a B3 domain TF, ZmABI19, that directly binds to the O2 promoter for transactivation. zmabi19 mutants displayed developmental defects in the endosperm and embryo, and mature kernels were opaque and reduced in size. The accumulation of zeins, starch and lipids dramatically decreased in zmabi19 mutants. RNA sequencing revealed an alteration of the nutrient reservoir activity and starch and sucrose metabolism in zmabi19 endosperms, and plant phytohormone signal transduction and lipid metabolism in zmabi19 embryos. Chromatin immunoprecipitation followed by sequencing coupled with differential expression analysis identified 106 high-confidence direct ZmABI19 targets. ZmABI19 directly regulates multiple key grain filling TFs including O2, Prolamine-box binding factor 1, ZmbZIP22, NAC130, and Opaque11 in the endosperm and Viviparous1 in the embryo. A number of phytohormone-related genes were also bound and regulated by ZmABI19. Our results demonstrate that ZmABI19 functions as a grain filling initiation regulator. ZmABI19 roles in coupling early endosperm and embryo development are also discussed.

Transcriptomic analysis of chicken immune response to infection of different doses of Newcastle disease vaccine.

Newcastle disease virus (NDV) is a contagious poultry paramyxovirus, leading to substantial economic losses to the poultry industry. Here, RNA-seq was carried out to investigate the altered expression of immune-related genes in chicken thymus within 96 h in response to NDV infection. In NDV-infected chicken thymus tissues, comparative transcriptome analysis revealed 1386 differentially expressed genes (DEGs) at 24 h with 989 up- and 397 down-regulated genes, 728 DEGs at 48 h with 567 up- and 161 down-regulated genes, 1514 DEGs at 72 h with 1016 up- and 498 down-regulated genes, and 1196 DEGs at 96 h with 522 up- and 674 down-regulated genes, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that these candidate targets mainly participate in biological processes or biochemical, metabolic and signal transduction processes. Notably, there is large enrichment in biological processes, cell components and metabolic processes, which may be related to NDV pathogenicity. In addition, the expression of five immune-related DEGs identified by RNA-seq was validated by quantitative real-time polymerase chain reaction (qRT-PCR). Our results indicated that the expression levels of AvBD5, IL16, IL22 and IL18R1 were obviously up-regulated, and Il-18 expression was also changed, but not significantly, which play key roles in the defense against NDV. Overall, we identified several candidate targets that may be involved in the regulation of NDV infection, which provide new insights into the complicated regulatory mechanisms of virus-host interactions, and explore new strategies for protecting chickens against the virus.

Genetic dissection of harvest index and related traits through genome-wide quantitative trait locus mapping in Brassica napus L.

The harvest index (HI) is the ratio of grain yield to the total biomass and represents the harvestable yield of crops. In Brassica napus, the HI is lower than that of other economically important crops, and limited relevant studies have been carried out regarding this issue. In this study, phenotypic analyses of 11 related traits showed the complexity of HI and the feasibility of cultivating desirable varieties with high HI. Quantitative trait loci (QTL) mapping based on a high-density genetic map identified 160 QTL, 163 epistatic loci pairs for HI and three closely related traits: seed yield (SY), biomass yield (BY) and plant height (PH), including two, five and three major QTL for HI, SY and PH, respectively. The related candidate genes underlying the QTL and epistatic loci with coding region variation were identified and investigated, including BnaA02g14010D, homologous to OsTB1, which functions as a negative regulator for lateral branching, and BnaA02g18890D, homologous to OsGW2, which controls grain width and weight. The complex correlation of HI with related traits, numerous QTL and epistatic loci and the candidate genes identified here provide new insights into the genetic architecture of HI, which might further enhance effective breeding strategies for yield improvement in rapeseed.

Genetic dissection of the mechanism of flowering time based on an environmentally stable and specific QTL in Brassica napus.

Flowering time is an important agronomic trait that is highly influenced by the environment. To elucidate the genetic mechanism of flowering time in rapeseed (Brassica napus L.), a genome-wide QTL analysis was performed in a doubled haploid population grown in winter, semi-winter and spring ecological conditions. Fifty-five consensus QTLs were identified after combining phenotype and genomic data, including 12 environment-stable QTLs and 43 environment-specific QTLs. Importantly, six major QTLs for flowering time were identified, of which two were considered environment-specific QTLs in spring ecological condition and four were considered environment-stable QTLs in winter and semi-winter ecological conditions. Through QTL comparison, 18 QTLs were colocalized with QTLs from six other published studies. Combining the candidate genes with their functional annotation, in 49 of 55 consensus QTLs, 151 candidate genes in B. napus corresponding to 95 homologous genes in Arabidopsis thaliana related to flowering were identified, including BnaC03g32910D (CO), BnaA02g12130D (FT) and BnaA03g13630D (FLC). Most of the candidate genes were involved in different flowering regulatory pathways. Based on re-sequencing and differences in sequence annotation between the two parents, we found that regions containing some candidate genes have numerous non-frameshift InDels and many non- synonymous mutations, which might directly lead to gene functional variation. Flowering time was negativly correlated with seed yield and thousand seed weight based on a QTL comparison of flowering time and seed yield traits, which has implications in breeding new early-maturing varieties of B. napus. Moreover, a putative flowering regulatory network was constructed, including the photoperiod, circadian clock, vernalization, autonomous and gibberellin pathways. Multiple copies of genes led to functional difference among the different copies of homologous genes, which also increased the complexity of the flowering regulatory networks. Taken together, the present results not only provide new insights into the genetic regulatory network underlying the control of flowering time but also improve our understanding of flowering time regulatory pathways in rapeseed.

3D Reconstruction of Lipid Droplets in the Seed of Brassica napus.

Rapeseed is one of the most important and widely cultured oilseed crops for food and nonfood purposes worldwide. Neutral lipids are stored in lipid droplets (LDs) as fuel for germination and subsequent seedling growth. Most of the LD detection in seeds was still in 2D levels, and some of the details might have been lost in previous studies. In the present work, the configuration of LDs in seeds was obtained by confocal imaging combined with 3D reconstruction technology in Brassica napus. The size and shape of LDs, LD numbers, cell interval spaces and cell size were observed and compared at 3D levels in the seeds of different materials with high and low oil content. It was also revealed that different cells located in the same tissue exhibited various oil contents according to the construction at the 3D level, which was not previously reported in B. napus. The present work provides a new way to understand the differential in cell populations and enhance the seed oil content at the single cell level within seeds.

Synteny analysis of genes and distribution of loci controlling oil content and fatty acid profile based on QTL alignment map in Brassica napus.

BACKGROUND: Deciphering the genetic architecture of a species is a good way to understand its evolutionary history, but also to tailor its profile for breeding elite cultivars with desirable traits. Aligning QTLs from diverse population in one map and utilizing it for comparison, but also as a basis for multiple analyses assure a stronger evidence to understand the genetic system related to a given phenotype. RESULTS: In this study, 439 genes involved in fatty acid (FA) and triacylglycerol (TAG) biosyntheses were identified in Brassica napus. B. napus genome showed mixed gene loss and insertion compared to B. rapa and B. oleracea, and C genome had more inserted genes. Identified QTLs for oil (OC-QTLs) and fatty acids (FA-QTLs) from nine reported populations were projected on the physical map of the reference genome "Darmor-bzh" to generate a map. Thus, 335 FA-QTLs and OC-QTLs could be highlighted and 82 QTLs were overlapping. Chromosome C3 contained 22 overlapping QTLs with all trait studied except for C18:3. In total, 218 candidate genes which were potentially involved in FA and TAG were identified in 162 QTLs confidence intervals and some of them might affect many traits. Also, 76 among these candidate genes were found inside 57 overlapping QTLs, and candidate genes for oil content were in majority (61/76 genes). Then, sixteen genes were found in overlapping QTLs involving three populations, and the remaining 60 genes were found in overlapping QTLs of two populations. Interaction network and pathway analysis of these candidate genes indicated ten genes that might have strong influence over the other genes that control fatty acids and oil formation. CONCLUSION: The present results provided new information for genetic basis of FA and TAG formation in B. napus. A map including QTLs from numerous populations was built, which could serve as reference to study the genome profile of B. napus, and new potential genes emerged which might affect seed oil. New useful tracks were showed for the selection of population or/and selection of interesting genes for breeding improvement purpose.

Genetic dissection of seed oil and protein content and identification of networks associated with oil content in Brassica napus.

High-density linkage maps can improve the precision of QTL localization. A high-density SNP-based linkage map containing 3207 markers covering 3072.7 cM of the Brassica napus genome was constructed in the KenC-8 × N53-2 (KNDH) population. A total of 67 and 38 QTLs for seed oil and protein content were identified with an average confidence interval of 5.26 and 4.38 cM, which could explain up to 22.24% and 27.48% of the phenotypic variation, respectively. Thirty-eight associated genomic regions from BSA overlapped with and/or narrowed the SOC-QTLs, further confirming the QTL mapping results based on the high-density linkage map. Potential candidates related to acyl-lipid and seed storage underlying SOC and SPC, respectively, were identified and analyzed, among which six were checked and showed expression differences between the two parents during different embryonic developmental periods. A large primary carbohydrate pathway based on potential candidates underlying SOC- and SPC-QTLs, and interaction networks based on potential candidates underlying SOC-QTLs, was constructed to dissect the complex mechanism based on metabolic and gene regulatory features, respectively. Accurate QTL mapping and potential candidates identified based on high-density linkage map and BSA analyses provide new insights into the complex genetic mechanism of oil and protein accumulation in the seeds of rapeseed.

Proteomic Dissection of Seed Germination and Seedling Establishment in Brassica napus.

The success of seed germination and establishment of a normal seedling are key determinants of plant species propagation. At present, only a few studies have focused on the genetic control of seed germination by using a proteomic approach in Brassica napus. In the present study, the protein expression pattern of seed germination was investigated using differential fluorescence two-dimensional gel electrophoresis in B. napus. One hundred and thirteen differentially expressed proteins (DEPs) that were mainly involved in storage (23.4%), energy metabolism (18.9%), protein metabolism (16.2%), defense/disease (12.6%), seed maturation (11.7%), carbohydrate metabolism (4.5%), lipid metabolism (4.5%), amino acids metabolism (3.6%), cell growth/division (3.6%), and some unclear functions (2.7%) were observed by proteomic analysis. Seventeen genes corresponding to 11 DEPs were identified within or near the associated linkage disequilibrium regions related to seed germination and vigor quantitative traits reported in B. napus in previous studies. The expression pattern of proteins showed that heterotrophic metabolism could be activated in the process of seed germination and that the onset of defense mechanisms might start during seed germination. These findings will help generate a more in-depth understanding of the mobilization of seed storage reserves and regulation mechanisms of the germination process in B. napus.