An ancient whole-genome duplication event and its contribution to flavor compounds in the tea plant (Camellia sinensis)

Tea, coffee, and cocoa are the three most popular nonalcoholic beverages in the world and have extremely high economic and cultural value. The genomes of four tea plant varieties have recently been sequenced, but there is some debate regarding the characterization of a whole-genome duplication (WGD) event in tea plants. Whether the WGD in the tea plant is shared with other plants in order Ericales and how it contributed to tea plant evolution remained unanswered. Here we re-analyzed the tea plant genome and provided evidence that tea experienced only WGD event after the core-eudicot whole-genome triplication (WGT) event. This WGD was shared by the Polemonioids-Primuloids-Core Ericales (PPC) sections, encompassing at least 17 families in the order Ericales. In addition, our study identified eight pairs of duplicated genes in the catechins biosynthesis pathway, four pairs of duplicated genes in the theanine biosynthesis pathway, and one pair of genes in the caffeine biosynthesis pathway, which were expanded and retained following this WGD. Nearly all these gene pairs were expressed in tea plants, implying the contribution of the WGD. This study shows that in addition to the role of the recent tandem gene duplication in the accumulation of tea flavor-related genes, the WGD may have been another main factor driving the evolution of tea flavor.

Evolution of an assembly factor-based subunit contributed to a novel NDH-PSI supercomplex formation in chloroplasts

Chloroplast NADH dehydrogenase-like (NDH) complex is structurally related to mitochondrial Complex I and forms a supercomplex with two copies of Photosystem I (the NDH-PSI supercomplex) via linker proteins Lhca5 and Lhca6. The latter was acquired relatively recently in a common ancestor of angiosperms. Here we show that NDH-dependent Cyclic Electron Flow 5 (NDF5) is an NDH assembly factor in Arabidopsis. NDF5 initiates the assembly of NDH subunits (PnsB2 and PnsB3) and Lhca6, suggesting that they form a contact site with Lhca6. Our analysis of the NDF5 ortholog in Physcomitrella and angiosperm genomes reveals the subunit PnsB2 to be newly acquired via tandem gene duplication of NDF5 at some point in the evolution of angiosperms. Another Lhca6 contact subunit, PnsB3, has evolved from a protein unrelated to NDH. The structure of the largest photosynthetic electron transport chain complex has become more complicated by acquiring novel subunits and supercomplex formation with PSI.

Reference Transcriptome Data in Silkworm Bombyx mori

Herein, we performed RNA-seq analysis of ten major tissues/subparts of silkworm larvae. The sequences were mapped onto the reference genome assembly and the reference transcriptome data were successfully constructed. The reference data provided a nearly complete sequence for sericin-1, a major silk gene with a complex structure. We also markedly improved the gene model for other genes. The transcriptomic expression was investigated in each tissue and a number of transcripts were identified that were exclusively expressed in tissues such as the testis. Transcripts strongly expressed in the midgut formed tight genomic clusters, suggesting that they originated from tandem gene duplication. Transcriptional factor genes expressed in specific tissues or the silk gland subparts were also identified. We successfully constructed reference transcriptome data in the silkworm and found that a number of transcripts showed unique expression profiles. These results will facilitate basic studies on the silkworm and accelerate its applications, which will contribute to further advances in lepidopteran and entomological research as well as the practical use of these insects.

Extensive structural variation in the Bowman-Birk inhibitor family in common wheat (Triticum aestivum L.)

Background Bowman-Birk inhibitors (BBI) are a family of serine-type protease inhibitors that modulate endogenous plant proteolytic activities during different phases of development. They also inhibit exogenous proteases as a component of plant defense mechanisms, and their overexpression can confer resistance to phytophagous herbivores and multiple fungal and bacterial pathogens. Dicot BBIs are multifunctional, with a “double-headed” structure containing two separate inhibitory loops that can bind and inhibit trypsin and chymotrypsin proteases simultaneously. By contrast, monocot BBIs have a non-functional chymotrypsin inhibitory loop, although they have undergone internal duplication events giving rise to proteins with multiple BBI domains. Results We used a Hidden Markov Model (HMM) profile-based search to identify 57 BBI genes in the common wheat (Triticum aestivum L.) genome. The BBI genes are unevenly distributed, with large gene clusters in the telomeric regions of homoeologous group 1 and 3 chromosomes that likely arose through a series of tandem gene duplication events. The genomes of wheat progenitors also contain contiguous clusters of BBI genes, suggesting this family underwent expansion before the domestication of common wheat. However, the BBI gene family varied in size among different cultivars, showing this family remains dynamic. Because of these expansions, the BBI gene family is larger in wheat than other monocots such as maize, rice and Brachypodium. We found BBI proteins in common wheat with intragenic homologous duplications of cysteine-rich functional domains, including one protein with four functional BBI domains. This diversification may expand the spectrum of target substrates. Expression profiling suggests that some wheat BBI proteins may be involved in regulating endogenous proteases during grain development, while others were induced in response to biotic and abiotic stresses, suggesting a role in plant defense. Conclusions Genome-wide characterization reveals that the BBI gene family in wheat is subject to a high rate of homologous tandem duplication and deletion events, giving rise to a diverse set of encoded proteins. This information will facilitate the functional characterization of individual wheat BBI genes to determine their role in wheat development and stress responses, and their potential application in breeding.

Specific Evolution and Gene Family Expansion of Complement 3 and Regulatory Factor H in Fish

The complement system comprises a large family of plasma proteins that play a central role in innate and adaptive immunity. To better understand the evolution of the complement system in vertebrates and the contribution of complement to fish immunity comprehensive in silico and expression analysis of the gene repertoire was made. Particular attention was given to C3 and the evolutionary related proteins C4 and C5 and to one of the main regulatory factors of C3b, factor H (Cfh). Phylogenetic and gene linkage analysis confirmed the standing hypothesis that the ancestral c3/c4/c5 gene duplicated early. The duplication of C3 (C3.1 and C3.2) and C4 (C4.1 and C4.2) was likely a consequence of the (1R and 2R) genome tetraploidization events at the origin of the vertebrates. In fish, gene number was not conserved and multiple c3 and cfh sequence related genes were encountered, and phylogenetic analysis of each gene generated two main clusters. Duplication of c3 and cfh genes occurred across the teleosts in a species-specific manner. In common, with other immune gene families the c3 gene expansion in fish emerged through a process of tandem gene duplication. Gilthead sea bream (Sparus aurata), had nine c3 gene transcripts highly expressed in liver although as reported in other fish, extra-hepatic expression also occurs. Differences in the sequence and protein domains of the nine deduced C3 proteins in the gilthead sea bream and the presence of specific cysteine and N-glycosylation residues within each isoform was indicative of functional diversity associated with structure. The diversity of C3 and other complement proteins as well as Cfh in teleosts suggests they may have an enhanced capacity to activate complement through direct interaction of C3 isoforms with pathogenic agents.

Extensive structural variation in the Bowman-Birk inhibitor family in common wheat (Triticum aestivum L.)

Background Bowman-Birk inhibitors (BBI) are a family of serine-type protease inhibitors that modulate endogenous plant proteolytic activities during different phases of development. They also inhibit exogenous proteases as a component of plant defense mechanisms, and their overexpression can confer resistance to phytophagous herbivores and multiple fungal and bacterial pathogens. Dicot BBIs are multifunctional, with a “double-headed” structure containing two separate inhibitory loops that can bind and inhibit trypsin and chymotrypsin proteases simultaneously. By contrast, monocot BBIs have a non-functional chymotrypsin inhibitory loop, although they have undergone internal duplication events giving rise to proteins with multiple BBI domains.Results We used a Hidden Markov Model (HMM) profile-based search to identify 57 BBI genes in the common wheat (Triticum aestivum L.) genome. The BBI genes are unevenly distributed, with large gene clusters in the telomeric regions of homeologous group 1 and 3 chromosomes that likely arose through a series of tandem gene duplication events. The genomes of wheat progenitors also contain contiguous clusters of BBI genes, suggesting this family underwent expansion before the domestication of common wheat. However, the BBI gene family varied in size among different cultivars, showing this family remains dynamic. Because of these expansions, the BBI gene family is larger in wheat than other monocots such as maize, rice and Brachypodium.We found BBI proteins in common wheat with intragenic homologous duplications of cysteine-rich functional domains, including one protein with four functional BBI domains. This diversification may expand the spectrum of target substrates. Expression profiling suggests that some wheat BBI proteins may be involved in regulating endogenous proteases during grain development, while others were induced in response to biotic and abiotic stresses, suggesting a role in plant defense.Conclusions Genome-wide characterization reveals that the BBI gene family in wheat is subject to a high rate of homologous tandem duplication and deletion events, giving rise to a diverse set of encoded proteins. This information will facilitate the functional characterization of individual wheat BBI genes to determine their role in wheat development and stress responses and their potential application in breeding.

Comparative Genomic and Proteomic Analyses of Three Widespread Phytophthora Species: Phytophthora chlamydospora, Phytophthora gonapodyides and Phytophthora pseudosyringae

The Phytophthora genus includes some of the most devastating plant pathogens. Here we report draft genome sequences for three ubiquitous Phytophthora species—Phytophthora chlamydospora, Phytophthora gonapodyides, and Phytophthora pseudosyringae. Phytophthora pseudosyringae is an important forest pathogen that is abundant in Europe and North America. Phytophthora chlamydospora and Ph. gonapodyides are globally widespread species often associated with aquatic habitats. They are both regarded as opportunistic plant pathogens. The three sequenced genomes range in size from 45 Mb to 61 Mb. Similar to other oomycete species, tandem gene duplication appears to have played an important role in the expansion of effector arsenals. Comparative analysis of carbohydrate-active enzymes (CAZymes) across 44 oomycete genomes indicates that oomycete lifestyles may be linked to CAZyme repertoires. The mitochondrial genome sequence of each species was also determined, and their gene content and genome structure were compared. Using mass spectrometry, we characterised the extracellular proteome of each species and identified large numbers of proteins putatively involved in pathogenicity and osmotrophy. The mycelial proteome of each species was also characterised using mass spectrometry. In total, the expression of approximately 3000 genes per species was validated at the protein level. These genome resources will be valuable for future studies to understand the behaviour of these three widespread Phytophthora species.

Whole genome comparisons of ergot fungi reveals the divergence and evolution of species within the genus Claviceps are the result of varying mechanisms driving genome evolution and host range expansion

The genus Claviceps has been known for centuries as an economically important fungal genera for pharmacology and agricultural research. Only recently have researchers begun to unravel the evolutionary history of the genus, with origins in South America and classification of four distinct sections through ecological, morphological, and metabolic features (Claviceps sects. Citrinae, Paspalorum, Pusillae, and Claviceps). The first three sections are additionally characterized by narrow host range, while sect. Claviceps is considered evolutionarily more successful and adaptable as it has the largest host range and biogeographical distribution. However, the reasons for this success and adaptability remain unclear. Our study elucidates factors influencing adaptability by sequencing and annotating 50 Claviceps genomes, representing 21 species, for a comprehensive comparison of genome architecture and plasticity in relation to host range potential. Our results show the trajectory from specialized one-speed genomes (sects. Citrinae and Paspalorum) towards adaptive two-speed genomes (sects. Pusillae and Claviceps) through co-localization of transposable elements around predicted effectors and a putative loss of repeat-induced point mutation resulting in unconstrained tandem gene duplication coinciding with increased host range potential and speciation. Alterations of genomic architecture and plasticity can substantially influence and shape the evolutionary trajectory of fungal pathogens and their adaptability. Furthermore, our study provides a large increase in available genomic resources to propel future studies of Claviceps in pharmacology and agricultural research, as well as, research into deeper understanding of the evolution of adaptable plant pathogens.

Integrative Analysis of the Core Fruit Lignification Toolbox in Pear Reveals Targets for Fruit Quality Bioengineering

Stone cell content is an important factor affecting pear fruit flavor. Lignin, a major component of pear stone cells, hinders the quality and value of commercial fruit. The completion of the Chinese white pear (Pyrus bretschneideri) genome sequence provides an opportunity to perform integrative analysis of the genes encoding the eleven protein families (i.e., PAL, C4H, 4CL, HCT, C3H, CSE, CCoAOMT, CCR, F5H, COMT, and CAD) in the phenylpropanoid pathway. Here, a systematic study based on expression patterns and phylogenetic analyses was performed to identify the members of each gene family potentially involved in the lignification in the Chinese white pear. The phylogenetic analysis suggested that 35 P. bretschneideri genes belong to bona fide lignification clade members. Compared to other plants, some multigene families are expanded by tandem gene duplication, such as HCT, C3H, COMT, and CCR. RNA sequencing was used to study the expression patterns of the genes in different tissues, including leaf, petal, bud, sepal, ovary, stem, and fruit. Eighteen genes presented a high expression in fruit, indicating that these genes may be involved in the biosynthesis of lignin in pear fruit. Similarly to what has been observed for Populus trichocarpa, a bimolecular fluorescence complementation (BiFC) experiment indicated that P. bretschneideri C3H and C4H might also interact with each other to regulate monolignol biosynthesis in P. bretschneideri, ultimately affecting the stone cell content in pear fruits. The identification of the major genes involved in lignin biosynthesis in pear fruits provides the basis for the development of strategies to improve fruit quality.

The Draft Genome of Kochia scoparia and the Mechanism of Glyphosate Resistance via Transposon-Mediated EPSPS Tandem Gene Duplication

Increased copy number of the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) gene confers resistance to glyphosate, the world’s most-used herbicide. There are typically three to eight EPSPS copies arranged in tandem in glyphosate-resistant populations of the weed kochia (Kochia scoparia). Here, we report a draft genome assembly from a glyphosate-susceptible kochia individual. Additionally, we assembled the EPSPS locus from a glyphosate-resistant kochia plant by sequencing select bacterial artificial chromosomes from a kochia bacterial artificial chromosome library. Comparing the resistant and susceptible EPSPS locus allowed us to reconstruct the history of duplication in the structurally complex EPSPS locus and uncover the genes that are coduplicated with EPSPS, several of which have a corresponding change in transcription. The comparison between the susceptible and resistant assemblies revealed two dominant repeat types. Additionally, we discovered a mobile genetic element with a FHY3/FAR1-like gene predicted in its sequence that is associated with the duplicated EPSPS gene copies in the resistant line. We present a hypothetical model based on unequal crossing over that implicates this mobile element as responsible for the origin of the EPSPS gene duplication event and the evolution of herbicide resistance in this system. These findings add to our understanding of stress resistance evolution and provide an example of rapid resistance evolution to high levels of environmental stress.