Kaining Hu

Genome-wide association study reveals the genetic architecture of flowering time in rapeseed (Brassica napus L.)

Flowering time adaptation is a major breeding goal in the allopolyploid species Brassica napus. To investigate the genetic architecture of flowering time, a genome-wide association study (GWAS) of flowering time was conducted with a diversity panel comprising 523 B. napus cultivars and inbred lines grown in eight different environments. Genotyping was performed with a Brassica 60K Illumina Infinium SNP array. A total of 41 single-nucleotide polymorphisms (SNPs) distributed on 14 chromosomes were found to be associated with flowering time, and 12 SNPs located in the confidence intervals of quantitative trait loci (QTL) identified in previous researches based on linkage analyses. Twenty-five candidate genes were orthologous to Arabidopsis thaliana flowering genes. To further our understanding of the genetic factors influencing flowering time in different environments, GWAS was performed on two derived traits, environment sensitivity and temperature sensitivity. The most significant SNPs were found near Bn-scaff_16362_1-p380982, just 13 kb away from BnaC09g41990D, which is orthologous to A. thaliana CONSTANS (CO), an important gene in the photoperiod flowering pathway. These results provide new insights into the genetic control of flowering time in B. napus and indicate that GWAS is an effective method by which to reveal natural variations of complex traits in B. napus.

Genome-Wide Association Study Dissecting the Genetic Architecture Underlying the Branch Angle Trait in Rapeseed ( Brassica napus L.)

The rapeseed branch angle is an important morphological trait because an adequate branch angle enables more efficient light capture under high planting densities. Here, we report that the average angle of the five top branches provides a reliable representation of the average angle of all branches. Statistical analyses revealed a significantly positive correlation between the branch angle and multiple plant-type and yield-related traits. The 60 K Brassica Infinium® single nucleotide polymorphism (SNP) array was utilized to genotype an association panel with 520 diverse accessions. A genome-wide association study was performed to determine the genetic architecture of branch angle, and 56 loci were identified as being significantly associated with the branch angle trait via three models, including a robust, novel, nonparametric Anderson-Darling (A-D) test. Moreover, these loci explained 51.1% of the phenotypic variation when a simple additive model was applied. Within the linkage disequilibrium (LD) decay ranges of 53 loci, we observed plausible candidates orthologous to documented Arabidopsis genes, such as LAZY1, SGR2, SGR4, SGR8, SGR9, PIN3, PIN7, CRK5, TIR1, and APD7. These results provide insight into the genetic basis of the branch angle trait in rapeseed and might facilitate marker-based breeding for improvements in plant architecture.

Genome-Wide Association Study Provides Insight into the Genetic Control of Plant Height in Rapeseed (Brassica napus L.)

Plant height is a key morphological trait of rapeseed. In this study, we measured plant height of a rapeseed population across six environments. This population contains 476 inbred lines representing the major Chinese rapeseed genepool and 44 lines from other countries. The 60K Brassica Infinium® SNP array was utilized to genotype the association panel. A genome-wide association study (GWAS) was performed via three methods, including a robust, novel, nonparametric Anderson–Darling (A–D) test. Consequently, 68 loci were identified as significantly associated with plant height (P < 5.22 × 10−5), and more than 70% of the loci (48) overlapped the confidence intervals of reported QTLs from nine mapping populations. Moreover, 24 GWAS loci were detected with selective sweep signals, which reflected the signatures of historical semi-dwarf breeding. In the linkage disequilibrium (LD) decay range up—and downstream of 65 loci (r2 > 0.1), we found plausible candidates orthologous to the documented Arabidopsis genes involved in height regulation. One significant association found by GWAS colocalized with the established height locus BnRGA in rapeseed. Our results provide insights into the genetic basis of plant height in rapeseed and may facilitate marker-based breeding.

Comparative Analysis of the Brassica napus Root and Leaf Transcript Profiling in Response to Drought Stress

Drought stress is one of the major abiotic factors affecting Brassica napus (B. napus) productivity. In order to identify genes of potential importance to drought stress and obtain a deeper understanding of the molecular mechanisms regarding the responses of B. napus to dehydration stress, we performed large-scale transcriptome sequencing of B. napus plants under dehydration stress using the Illumina sequencing technology. In this work, a relatively drought tolerant B. napus line, Q2, identified in our previous study, was used. Four cDNA libraries constructed from mRNAs of control and dehydration-treated root and leaf were sequenced by Illumina technology. A total of 6018 and 5377 differentially expressed genes (DEGs) were identified in root and leaf. In addition, 1745 genes exhibited a coordinated expression profile between the two tissues under drought stress, 1289 (approximately 74%) of which showed an inverse relationship, demonstrating different regulation patterns between the root and leaf. The gene ontology (GO) enrichment test indicated that up-regulated genes in root were mostly involved in “stimulus” “stress” biological process, and activated genes in leaf mainly functioned in “cell” “cell part” components. Furthermore, a comparative network related to plant hormone signal transduction and AREB/ABF, AP2/EREBP, NAC, WRKY and MYC/MYB transcription factors (TFs) provided a view of different stress tolerance mechanisms between root and leaf. Some of the DEGs identified may be candidates for future research aimed at detecting drought-responsive genes and will be useful for understanding the molecular mechanisms of drought tolerance in root and leaf of B. napus.

Neofunctionalization of Duplicated Tic40 Genes Caused a Gain-of-Function Variation Related to Male Fertility in Brassica oleracea Lineages

Divergent duplicated genes in the Brassica oleracea gain a novel function that could restore the fertility of a Brassica napus male-sterile mutant. Gene duplication followed by functional divergence in the event of polyploidization is a major contributor to evolutionary novelties. The Brassica genus evolved from a common ancestor after whole-genome triplication. Here, we studied the evolutionary and functional features of Brassica spp. homologs to Tic40 (for translocon at the inner membrane of chloroplasts with 40 kDa). Four Tic40 loci were identified in allotetraploid Brassica napus and two loci in each of three basic diploid Brassica spp. Although these Tic40 homologs share high sequence identities and similar expression patterns, they exhibit altered functional features. Complementation assays conducted on Arabidopsis thaliana tic40 and the B. napus male-sterile line 7365A suggested that all Brassica spp. Tic40 homologs retain an ancestral function similar to that of AtTic40, whereas BolC9.Tic40 in Brassica oleracea and its ortholog in B. napus, BnaC9.Tic40, in addition, evolved a novel function that can rescue the fertility of 7365A. A homologous chromosomal rearrangement placed bnac9.tic40 originating from the A genome (BraA10.Tic40) as an allele of BnaC9.Tic40 in the C genome, resulting in phenotypic variation for male sterility in the B. napus near-isogenic two-type line 7365AB. Assessment of the complementation activity of chimeric B. napus Tic40 domain-swapping constructs in 7365A suggested that amino acid replacements in the carboxyl terminus of BnaC9.Tic40 cause this functional divergence. The distribution of these amino acid replacements in 59 diverse Brassica spp. accessions demonstrated that the neofunctionalization of Tic40 is restricted to B. oleracea and its derivatives and thus occurred after the divergence of the Brassica spp. A, B, and C genomes.

Altered Transcription and Neofunctionalization of Duplicated Genes Rescue the Harmful Effects of a Chimeric Gene in Brassica napus.

Brassica plants have altered transcript levels and coevolved neofunctionalization of duplicated genes to block developmental defects resulting from detrimental chimeric genes. Chimeric genes contribute to the evolution of diverse functions in plants and animals. However, new chimeric genes also increase the risk of developmental defects. Here, we show that the chimeric gene Brassica napus male sterile 4 (Bnams4b) is responsible for genic male sterility in the widely used canola line 7365A (Bnams3 ms3ms4bms4b). Bnams4b originated via exon shuffling ∼4.6 million years ago. It causes defects in the normal functions of plastids and induces aborted anther formation and/or albino leaves and buds. Evidence of the age of the mutation, its tissue expression pattern, and its sublocalization indicated that it coevolved with BnaC9.Tic40 (BnaMs3). In Arabidopsis thaliana, Bnams4b results in complete male sterility that can be rescued by BnaC9.Tic40, suggesting that BnaC9.Tic40 might restore fertility through effects on protein level. Another suppressor gene, Bnams4a, rescues sterility by reducing the level of transcription of Bnams4b. Our results suggest that Brassica plants have coevolved altered transcription patterns and neofunctionalization of duplicated genes that can block developmental defects resulting from detrimental chimeric genes.

Overexpression of the Novel Arabidopsis Gene At5g02890 Alters Inflorescence Stem Wax Composition and Affects Phytohormone Homeostasis

The cuticle is composed of cutin and cuticular wax. It covers the surfaces of land plants and protects them against environmental damage. At5g02890 encodes a novel protein in Arabidopsis thaliana. In the current study, protein sequence analysis showed that At5g02890 is highly conserved in the Brassicaceae. Arabidopsis lines overexpressing At5g02890 (OE-At5g02890 lines) and an At5g02890 orthologous gene from Brassica napus (OE-Bn1 lines) exhibited glossy stems. Chemical analysis revealed that overexpression of At5g02890 caused significant reductions in the levels of wax components longer than 28 carbons (C28) in inflorescence stems, whereas the levels of wax molecules of chain length C28 or shorter were significantly increased. Transcriptome analysis indicated that nine of 11 cuticular wax synthesis-related genes with different expression levels in OE-At5g02890 plants are involved in very-long-chain fatty acid (VLCFA) elongation. At5g02890 is localized to the endoplasmic reticulum (ER), which is consistent with its function in cuticular wax biosynthesis. These results demonstrate that the overexpression of At5g02890 alters cuticular wax composition by partially blocking VLCFA elongation of C28 and higher. In addition, detailed analysis of differentially expressed genes associated with plant hormones and endogenous phytohormone levels in wild-type and OE-At5g02890 plants indicated that abscisic acid (ABA), jasmonic acid (JA), and jasmonoyl-isoleucine (JA-Ile) biosynthesis, as well as polar auxin transport, were also affected by overexpression of At5g02890. Taken together, these findings indicate that overexpression of At5g02890 affects both cuticular wax biosynthesis and phytohormone homeostasis in Arabidopsis.

Trilocular phenotype in Brassica juncea L. resulted from interruption of CLAVATA1 gene homologue (BjMc1) transcription

As a desirable agricultural trait, multilocular trait of rapeseed (Brassica rapa; Brassica napus; Brassica juncea), always represents higher yield per plant compared with bilocular plants. We previously isolated a trilocular gene locus, Bjmc1, and identified a set of molecular markers linked to the trilocular gene. With a map-based cloning, we identified that the BjMc1 was located in B genome of Brassica juncea, and it was a CLAVATA1 (CLV1) gene homologue. The insertion of a copia-LTR retrotransposable element 1 (RTE1) into the coding region of BjMc1 interrupted its transcription in rapeseed, leading to the trilocular phenotype. Phylogenetic analysis showed that Mc1 genes were conserved and widespread in land plants. Two amino acid sites had undergone positive selection in the ancestor of Mc1 genes, and then purifying selection was the dominant force after the divergence of dicots and monocots from their common ancestor in the evolutionary process, indicating that Mc1 genes are conserved in modern land plants. Our results provided new insights in molecular regulatory mechanism of multilocularity in rapeseed, and better understanding of molecular mechanism in crop yield improvement.

Transcriptomic Analysis of Seed Coats in Yellow-Seeded Brassica napus Reveals Novel Genes That Influence Proanthocyanidin Biosynthesis

Yellow seeds are a favorable trait for Brassica crops breeding due to better quality than their black-seeded counterparts. Here, we compared the Brassica napus seed coat transcriptomes between yellow- and brown-seeded near-isogenic lines (Y-NIL and B-NIL) that were developed from the resynthesized yellow-seeded line No. 2127-17. A total of 4,974 differentially expressed genes (DEG) were identified during seed development, involving 3,128 up-regulated and 1,835 down-regulated genes in yellow seed coats. Phenylpropanoid and flavonoid biosynthesis pathways were enriched in down-regulated genes, whereas the top two pathways for up-regulated genes were plant–pathogen interaction and plant hormone signal transduction. Twelve biosynthetic genes and three regulatory genes involved in the flavonoid pathway exhibited similar expression patterns in seed coats during seed development, of which the down-regulation mainly contributed to the reduction of proanthocyanidins (PAs) in yellow seed coats, indicating that these genes associated with PA biosynthesis may be regulated by an unreported common regulator, possibly corresponding to the candidate for the dominant black-seeded gene D in the NILs. Three transcription factor (TF) genes, including one bHLH gene and two MYB-related genes that are located within the previous seed coat color quantitative trait locus (QTL) region on chromosome A09, also showed similar developmental expression patterns to the key PA biosynthetic genes and they might thus potentially involved participate in flavonoid biosynthesis regulation. Our study identified novel potential TFs involved in PAs accumulation and will provide pivotal information for identifying the candidate genes for seed coat color in B. napus.

Breeding signature of combining ability improvement revealed by a genomic variation map from recurrent selection population in Brassica napus

Combining ability is crucial for parent selection in crop hybrid breeding. The present investigation and results had revealed the underlying genetic factors which might contribute in adequate combining ability, further assisting in enhancing heterosis and stability. Here, we conducted a large-scale analysis of genomic variation in order to define genomic regions affecting the combining ability in recurrent selection population of rapeseed. A population of 175 individuals was genotyped with the Brassica60K SNP chip. 525 hybrids were assembled with three different testers and used to evaluate the general combining ability (GCA) in three environments. By detecting the changes of the genomic variation, we identified 376 potential genome regions, spanning 3.03% of rapeseed genome which provided QTL-level resolution on potentially selected variants. More than 96% of these regions were located in the C subgenome, indicating that C subgenome had sustained stronger selection pressure in the breeding program than the A subgenome. In addition, a high level of linkage disequilibrium in rapeseed genome was detected, suggesting that marker-assisted selection for the population improvement might be easily implemented. This study outlines the evidence for high GCA on a genomic level and provided underlying molecular mechanism for recurrent selection improvement in B. napus.