Jiaqin Shi

Genome-wide microsatellite characterization and marker development in the sequenced Brassica crop species.

Although much research has been conducted, the pattern of microsatellite distribution has remained ambiguous, and the development/utilization of microsatellite markers has still been limited/inefficient in Brassica, due to the lack of genome sequences. In view of this, we conducted genome-wide microsatellite characterization and marker development in three recently sequenced Brassica crops: Brassica rapa, Brassica oleracea and Brassica napus. The analysed microsatellite characteristics of these Brassica species were highly similar or almost identical, which suggests that the pattern of microsatellite distribution is likely conservative in Brassica. The genomic distribution of microsatellites was highly non-uniform and positively or negatively correlated with genes or transposable elements, respectively. Of the total of 115 869, 185 662 and 356 522 simple sequence repeat (SSR) markers developed with high frequencies (408.2, 343.8 and 356.2 per Mb or one every 2.45, 2.91 and 2.81 kb, respectively), most represented new SSR markers, the majority had determined physical positions, and a large number were genic or putative single-locus SSR markers. We also constructed a comprehensive database for the newly developed SSR markers, which was integrated with public Brassica SSR markers and annotated genome components. The genome-wide SSR markers developed in this study provide a useful tool to extend the annotated genome resources of sequenced Brassica species to genetic study/breeding in different Brassica species.

Evolutionary Dynamics of Microsatellite Distribution in Plants: Insight from the Comparison of Sequenced Brassica, Arabidopsis and Other Angiosperm Species

Despite their ubiquity and functional importance, microsatellites have been largely ignored in comparative genomics, mostly due to the lack of genomic information. In the current study, microsatellite distribution was characterized and compared in the whole genomes and both the coding and non-coding DNA sequences of the sequenced Brassica, Arabidopsis and other angiosperm species to investigate their evolutionary dynamics in plants. The variation in the microsatellite frequencies of these angiosperm species was much smaller than those for their microsatellite numbers and genome sizes, suggesting that microsatellite frequency may be relatively stable in plants. The microsatellite frequencies of these angiosperm species were significantly negatively correlated with both their genome sizes and transposable elements contents. The pattern of microsatellite distribution may differ according to the different genomic regions (such as coding and non-coding sequences). The observed differences in many important microsatellite characteristics (especially the distribution with respect to motif length, type and repeat number) of these angiosperm species were generally accordant with their phylogenetic distance, which suggested that the evolutionary dynamics of microsatellite distribution may be generally consistent with plant divergence/evolution. Importantly, by comparing these microsatellite characteristics (especially the distribution with respect to motif type) the angiosperm species (aside from a few species) all clustered into two obviously different groups that were largely represented by monocots and dicots, suggesting a complex and generally dichotomous evolutionary pattern of microsatellite distribution in angiosperms. Polyploidy may lead to a slight increase in microsatellite frequency in the coding sequences and a significant decrease in microsatellite frequency in the whole genome/non-coding sequences, but have little effect on the microsatellite distribution with respect to motif length, type and repeat number. Interestingly, several microsatellite characteristics seemed to be constant in plant evolution, which can be well explained by the general biological rules.

Linkage and regional association analysis reveal two new tightly-linked major-QTLs for pod number and seed number per pod in rapeseed ( Brassica napus L. )

To facilitate the pseudochromosomes assembly and gene cloning in rapeseed, we developed a reference genetic population/map (named BnaZNF2) from two sequenced cultivars, Zhongshuang11 and No.73290, those exhibit significant differences in many traits, particularly yield components. The BnaZNF2 genetic map exhibited perfect collinearity with the physical map of B. napus, indicating its high quality. Comparative mapping revealed several genomic rearrangements between B. napus and B. rapa or B. oleracea. A total of eight and 16 QTLs were identified for pod number and seed number per pod, respectively, and of which three and five QTLs are identical to previously identified ones, whereas the other five and 11 are novel. Two new major QTL respectively for pod number and seed number per pod, qPN.A06-1 and qSN.A06-1 (R2 = 22.8% and 32.1%), were colocalised with opposite effects, and only qPN.A06-1 was confirmed and narrowed by regional association analysis to 180 kb including only 33 annotated genes. Conditional QTL analysis and subsequent NILs test indicated that tight linkage, rather than pleiotropy, was the genetic causation of their colocalisation. Our study demonstrates potential of this reference genetic population/map for precise QTL mapping and as a base for positional gene cloning in rapeseed.

Genetic architecture and mechanism of seed number per pod in rapeseed: elucidated through linkage and near-isogenic line analysis.

Seed number per pod (SNPP) is one of the major yield components and breeding targets in rapeseed that shows great variation and is invaluable for genetic improvement. To elucidate the genetic architecture and uncover the mechanism of SNPP, we identified five quantitative trait loci (QTLs) using the BnaZNRIL population, which were integrated with those of previous studies by physical map to demonstrate a complex and relatively complete genetic architecture of SNPP. A major QTL, qSN.A6, was successfully fine-mapped from 1910 to 267 kb using near-isogenic line (NIL). In addition, qSN.A6 exhibited an antagonistic pleiotropy on seed weight (SW), which is caused by a physiological interaction in which SNPP acts “upstream” of SW. Because the negative effect of qSN.A6 on SW cannot fully counteract its positive effect on SNPP, it also enhanced the final yield (17.4%), indicating its great potential for utilization in breeding. The following genetic and cytological experiments further confirmed that the different rate of ovule abortion was responsible for the ~5 seed difference between Zhongshuang11 and NIL-qSN.A6. This systematic approach to dissecting the comprehensive genetic architecture of SNPP and characterizing the underlying mechanism has advanced the understanding of SNPP and will facilitate the development of high-yield cultivars.

The Natural Variation of Seed Weight Is Mainly Controlled by Maternal Genotype in Rapeseed (Brassica napus L.).

Seed weight is a very important and complex trait in rapeseed (Brassica napus L.). The seed weight of rapeseed shows great variation in its natural germplasm resources; however, the morphological, cytological and genetic causes of this variation have remained unclear. In the present study, nine highly pure inbred rapeseed lines with large seed weight variation and different genetic backgrounds were selected for morphological, cytological and genetic studies on seed weight. The results showed the following: (1) Seed weight showed an extremely significant correlation and coordinated variation with seed size (including seed diameter, seed surface area and seed volume), but it showed no significant correlation with bulk density, which suggests that seed weight is determined by size rather than bulk density. (2) Seed weight showed a higher correlation with the cell numbers of seed coats and cotyledons than the cell sizes of seed coats and cotyledons, which suggests that cell number is more tightly correlated with final seed weight. (3) Seed weight was mainly controlled by the maternal genotype, with little or no xenia and cytoplasmic effects. This is the first report on the morphological and cytological causes of seed weight natural variation in rapeseed. We concluded that the natural variation of seed weight is mainly controlled by maternal genotype. This finding lays a foundation for genetic and breeding studies of seed weight in rapeseed and opens a new field of research on the regulation of seed traits in plants.

An integrated analysis of QTL mapping and RNA sequencing provides further insights and promising candidates for pod number variation in rapeseed (Brassica napus L.)

AbstractBackgroundAs the most important yield component in rapeseed (Brassica napus L.), pod number is determined by a series of successive growth and development processes. Pod number shows extensive variation in rapeseed natural germplasm, which is valuable for genetic improvement. However, the genetic and especially the molecular mechanism for this kind of variation are poorly understood. In this study, we conducted QTL mapping and RNA sequencing, respectively, using the BnaZNRIL population and its two parental cultivars Zhongshuang11 and No.73290 which showed significant difference in pod number, primarily due to the difference in floral organ number.ResultA total of eight QTLs for pod number were identified using BnaZNRIL population with a high-density SNP linkage map, each was distributed on seven linkage groups and explained 5.8–11.9% of phenotypic variance. Then, they were integrated with those previously detected in BnaZNF2 population (deriving from same parents) and resulted in 15 consensus-QTLs. Of which, seven QTLs were identical to other studies, whereas the other eight should be novel.RNA sequencing of the shoot apical meristem (SAM) at the formation stage of floral bud primordia identified 9135 genes that were differentially expressed between the two parents. Gene ontology (GO) analysis showed that the top two enriched groups were S-assimilation, providing an essential nutrient for the synthesis of diverse metabolites, and polyamine metabolism, serving as second messengers that play an essential role in flowering genes initiation. KEGG analysis showed that the top three overrepresented pathways were carbohydrate (707 genes), amino acid (390 genes) and lipid metabolisms (322 genes). In silico mapping showed that 647 DEGs were located within the confidence intervals of 15 consensus QTLs. Based on annotations of Arabidopsis homologs corresponding to DEGs, nine genes related to meristem growth and development were considered as promising candidates for six QTLs.ConclusionIn this study, we discovered the first repeatable major QTL for pod number in rapeseed. In addition, RNA sequencing was performed for SAM in rapeseed, which provides new insights into the determination of floral organ number. Furthermore, the integration of DEGs and QTLs identified promising candidates for further gene cloning and mechanism study.

Microarray Expression Analysis of the Main Inflorescence in Brassica napus

The effect of the number of pods on the main inflorescence (NPMI) on seed yield in Brassica napus plants grown at high density is a topic of great economic and scientific interest. Here, we sought to identify patterns of gene expression that determine the NPMI during inflorescence differentiation. We monitored gene expression profiles in the main inflorescence of two B. napus F6 RIL pools, each composed of nine lines with a low or high NPMI, and their parental lines, Zhongshuang 11 (ZS11) and 73290, using a Brassica 90K elements oligonucleotide array. We identified 4,805 genes that were differentially expressed (≥1.5 fold-change) between the low- and high-NPMI samples. Of these, 82.8% had been annotated and 17.2% shared no significant homology with any known genes. About 31 enriched GO clusters were identified amongst the differentially expressed genes (DEGs), including those involved in hormone responses, development regulation, carbohydrate metabolism, signal transduction, and transcription regulation. Furthermore, 92.8% of the DEGs mapped to chromosomes that originated from B. rapa and B. oleracea, and 1.6% of the DEGs co-localized with two QTL intervals (PMI10 and PMI11) known to be associated with the NPMI. Overexpression of BnTPI, which co-localized with PMI10, in Arabidopsis suggested that this gene increases the NPMI. This study provides insight into the molecular factors underlying inflorescence architecture, NPMI determination and, consequently, seed yield in B. napus.

Genetic Dissection of Root Morphological Traits Related to Nitrogen Use Efficiency in Brassica napus L. under Two Contrasting Nitrogen Conditions

As the major determinant for nutrient uptake, root system architecture (RSA) has a massive impact on nitrogen use efficiency (NUE). However, little is known the molecular control of RSA as related to NUE in rapeseed. Here, a rapeseed recombinant inbred line population (BnaZNRIL) was used to investigate root morphology (RM, an important component for RSA) and NUE-related traits under high-nitrogen (HN) and low-nitrogen (LN) conditions by hydroponics. Data analysis suggested that RM-related traits, particularly root size had significantly phenotypic correlations with plant dry biomass and N uptake irrespective of N levels, but no or little correlation with N utilization efficiency (NUtE), providing the potential to identify QTLs with pleiotropy or specificity for RM- and NUE-related traits. A total of 129 QTLs (including 23 stable QTLs, which were repeatedly detected at least two environments or different N levels) were identified and 83 of them were integrated into 22 pleiotropic QTL clusters. Five RM-NUE, ten RM-specific and three NUE-specific QTL clusters with same directions of additive-effect implied two NUE-improving approaches (RM-based and N utilization-based directly) and provided valuable genomic regions for NUE improvement in rapeseed. Importantly, all of four major QTLs and most of stable QTLs (20 out of 23) detected here were related to RM traits under HN and/or LN levels, suggested that regulating RM to improve NUE would be more feasible than regulating N efficiency directly. These results provided the promising genomic regions for marker-assisted selection on RM-based NUE improvement in rapeseed.

Maternal control of seed weight in rapeseed (Brassica napus L.): the causal link between the size of pod (mother, source) and seed (offspring, sink)

Summary Seed size/weight is one of the key traits related to plant domestication and crop improvement. In rapeseed (Brassica napus L.) germplasm, seed weight shows extensive variation, but its regulatory mechanism is poorly understood. To identify the key mechanism of seed weight regulation, a systematic comparative study was performed. Genetic, morphological and cytological evidence showed that seed weight was controlled by maternal genotype, through the regulation of seed size mainly via cell number. The physiological evidence indicated that differences in the pod length might result in differences in pod wall photosynthetic area, carbohydrates and the final seed weight. We also identified two pleiotropic major quantitative trait loci that acted indirectly on seed weight via their effects on pod length. RNA‐seq results showed that genes related to pod development and hormones were significantly differentially expressed in the pod wall; genes related to development, cell division, nutrient reservoir and ribosomal proteins were all up‐regulated in the seeds of the large‐seed pool. Finally, we proposed a potential seed weight regulatory mechanism that is specific to rapeseed and novel in plants. The results demonstrate a causal link between the size of the pod (mother, source) and the seed (offspring, sink) in rapeseed, which provides novel insight into the maternal control of seed weight and will open a new research field in plants.

Genetic and Cytological Analyses of the Natural Variation of Seed Number per Pod in Rapeseed (Brassica napus L.)

Seed number is one of the key traits related to plant evolution/domestication and crop improvement/breeding. In rapeseed germplasm, the seed number per pod (SNPP) shows a very wide variation from several to nearly 30; however, the underlying causations/mechanisms for this variation are poorly known. In the current study, the genetic and cytological bases for the natural variation of SNPP in rapeseed was firstly and systematically investigated using the representative four high-SNPP and five low-SNPP lines. The results of self- or cross-pollination experiment between the high- and low-SNPP lines showed that the natural variation of SNPP was mainly controlled by maternal effect (mean = 0.79), followed by paternal effect (mean = 0.21). Analysis of the data using diploid seed embryo–cytoplasmic–maternal model further showed that the maternal genotype, embryo, and cytoplasm effects, respectively, explained 47.6, 35.2, and 7.5% of the genetic variance. In addition, the analysis of combining ability showed that for the SNPP of hybrid F1 was mainly determined by the general combining ability of parents (63.0%), followed by special combining ability of parental combination (37.0%). More importantly, the cytological observation showed that the SNPP difference between the high- and low-SNPP lines was attributable to the accumulative differences in its components. Of which, the number of ovules, the proportion of fertile ovules, the proportion of fertile ovules to be fertilized, and the proportion of fertilized ovules to develop into seeds accounted for 30.7, 18.2, 7.1, and 43.9%, respectively. The accordant results of both genetic and cytological analyses provide solid evidences and systematic insights to further understand the mechanisms underlying the natural variation of SNPP, which will facilitate the development of high-yield cultivars in rapeseed.