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

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.Polyploidy or whole-genome duplication is thought to be a prominent evolutionary force in eukaryotes (Wolfe, 2001; Udall and Wendel, 2006), especially for flowering plants (Blanc and Wolfe, 2004; Van de Peer et al., 2009). Almost 95% of angiosperms show evidence of having undergone at least one round of whole-genome duplication in their evolutionary history, suggesting that most extant diploid flowering plants have evolved from ancient polyploids (Cui et al., 2006; Soltis et al., 2009). Gene duplications in the event of polyploidization provide sources for evolutionary novelties that could benefit plants (Lukens et al., 2004; Chen, 2007). Divergence after gene duplication could result in three primary evolutionary fates of duplicated genes: pseudogenization, neofunctionalization, and subfunctionalization (Force et al., 1999; Conant and Wolfe, 2008; Liu and Adams, 2010). Pseudogenization implies that duplicated genes with redundant functions lose their function by accumulating negative mutations; neofunctionalization denotes that the redundant gene evolves a new adaptive function; while subfunctionalization causes the duplicated genes to adopt a different part of the function of an ancestral gene (Rodríguez-Trelles et al., 2003; Flagel and Wendel, 2009; Liu and Adams, 2010).The Brassica genus consists of three basic diploid species: Brassica rapa (AA; n = 10), Brassica nigra (BB; n = 8), and Brassica oleracea (CC; n = 9), and their derivative allotetraploid species: Brassica juncea (AABB; n = 18), Brassica napus (AACC; n = 19), and Brassica carinata (BBCC; n = 17; Beilstein et al., 2006). Comparative genetic mapping demonstrated that these Brassica spp., which diverged from Arabidopsis thaliana approximately 20 to 40 million years ago (Lagercrantz and Lydiate, 1996; Blanc et al., 2003; Town et al., 2006), descended from a common ancestor after whole-genome triplication (Parkin et al., 2002). Collinear comparison showed that for each of 24 ancestral genomic blocks defined in the ancestral karyotype in A. thaliana, three syntenic copies were identified in each of the diploid Brassica spp. genomes, with only one exception (Schranz et al., 2006; Wang et al., 2011; Cheng et al., 2013). Fractionation (gene loss from homologous genomic regions) and chromosomal rearrangements were prevalent in the diploidization process of the hexaploid Brassica spp. common ancestor (Lagercrantz, 1998; Town et al., 2006; Ziolkowski et al., 2006; Mun et al., 2009). Based on gene density differences caused by varying gene loss rates in the three collinear genomic block copies, these genomic blocks were classified into three subgenomes in the diploid Brassica spp. genomes: the least fractionated (LF), the medium fractionated (MF1), and the most fractionated (MF2) subgenomes (Wang et al., 2011; Tang and Lyons, 2012; Cheng et al., 2013). The different rates of gene loss of the three subgenomes support a two-step origin of the Brassiceae ancestral genome involving a tetraploidization process followed by substantial fractionation of the subgenomes MF1 and MF2 and more recently hybridization with the third subgenome LF to form a hexaploid (Wang et al., 2011; Tang et al., 2012). Although collinearity and changes in genomic structure, including duplications, deletions, and rearrangements of the Brassica spp. and A. thaliana are well studied, there is limited knowledge of the molecular and functional divergence of duplicated or homologous genes in the Brassica spp.To utilize heterosis in B. napus breeding, hybrid production is based mainly on male sterility. Currently, the recessive epistatic genic male-sterile three-type line system 7365ABC is widely used for oilseed heterosis due to its advantages, producing 100% sterile offspring for realizing the triple-cross hybrid (Huang et al., 2007; Xia et al., 2012). The male sterility in this system is controlled by two genes, a recessive male sterile gene Bnms3 and a epistatic gene BnRf (Zhou et al., 2012). Bnms3 was recently reported to be a homolog to A. thaliana Tic40 (Dun et al., 2011). Tic40 was identified as a member of the TIC (for translocon at the inner envelope membrane of chloroplasts) complex that functions as a cochaperone to coordinate Tic110 and the stromal chaperone heat shock protein93 (Hsp93) (Chou et al., 2003, 2006). The A. thaliana tic40 mutant displayed a chlorotic phenotype throughout development (Chou et al., 2003) but a male-fertile phenotype with mature pollen grains (Dun et al., 2011). Interestingly, one allele of Bnms3, BnaC.Tic40, can rescue the fertility of the B. napus male-sterile line 7365A (Dun et al., 2011).In this study, we identified and characterized Tic40 homologs in B. napus and three basic diploid Brassica spp. We suggested that neofunctionalization of BnaC9.Tic40 after the divergence of the Brassica spp. A, B, and C genomes was caused by amino acid replacements in the C terminus of BnaC9.Tic40. In addition, we validated that the allelic genes BnaC9.Tic40 (equivalent to BnaC.Tic40 as described by Dun et al. [2011]) and bnac9.tic40 originate from the C and A genomes, respectively, and became allelic due to a homologous chromosomal rearrangement in B. napus. These results provide further knowledge for the effective utilization of the restoring gene of 7365A and a better insight into the functional divergence of homologous duplicated genes in paleoploid Brassica spp.