Genomic basis of multiphase evolution driving divergent selection of zinc-finger homeodomain genes.

Gene families divergently evolve and become adapted as different genes with specific structures and functions in living organisms. We performed comprehensive structural and functional analyses of Zinc-finger homeodomain genes (ZF-HDs), including Mini zinc-finger genes (MIFs) and Zinc-finger with homeodomain genes (ZHDs), displaying competitive functions each other. Intensive annotation updates for 90 plant genomes verified that most MIFs (MIF-Is) exhibited distinct motif compositions from ZHDs, although some MIFs (MIF-Zs) contained ZHD-specific motifs. Phylogenetic analyses suggested that MIF-Zs and ZHDs originated from the same ancestral gene, whereas MIF-Is emerged from a distinct progenitor. We used a gene-editing system to identify a novel function of MIF-Is in rice: regulating the surface material patterns in anthers and pollen through transcriptional regulation by interacting ZHDs. Kingdom-wide investigations determined that (i) ancestral MIFs diverged into MIF-Is and MIF-Zs in the last universal common ancestor, (ii) integration of HD into the C-terminal of MIF-Zs created ZHDs after emergence of green plants and (iii) MIF-Is and ZHDs subsequently expanded independently into specific plant lineages, with additional formation of MIF-Zs from ZHDs. Our comprehensive analysis provides genomic evidence for multiphase evolution driving divergent selection of ZF-HDs.

Identification and expression analyses of B3 genes reveal lineage-specific evolution and potential roles of REM genes in pepper.

BACKGROUND: The B3 gene family, one of the largest plant-specific transcription factors, plays important roles in plant growth, seed development, and hormones. However, the B3 gene family, especially the REM subfamily, has not been systematically and functionally studied. RESULTS: In this study, we performed genome-wide re-annotation of B3 genes in five Solanaceae plants, Arabidopsis thaliana, and Oryza sativa, and finally predicted 1,039 B3 genes, including 231 (22.2%) newly annotated genes. We found a striking abundance of REM genes in pepper species (Capsicum annuum, Capsicum baccatum, and Capsicum chinense). Comparative motif analysis revealed that REM and other subfamilies (ABI3/VP1, ARF, RAV, and HSI) consist of different amino acids. We verified that the large number of REM genes in pepper were included in the specific subgroup (G8) through the phylogenetic analysis. Chromosome location and evolutionary analyses suggested that the G8 subgroup genes evolved mainly via a pepper-specific recent tandem duplication on chromosomes 1 and 3 after speciation between pepper and other Solanaceae. RNA-seq analyses suggested the potential functions of REM genes under salt, heat, cold, and mannitol stress conditions in pepper (C. annuum). CONCLUSIONS: Our study provides evolutionary and functional insights into the REM gene family in pepper.

Comparative and expression analyses of AP2/ERF genes reveal copy number expansion and potential functions of ERF genes in Solanaceae.

BACKGROUND: The AP2/ERF gene family is a superfamily of transcription factors that are important in the response of plants to abiotic stress and development. However, comprehensive research of the AP2/ERF genes in the Solanaceae family is lacking. RESULTS: Here, we updated the annotation of AP2/ERF genes in the genomes of eight Solanaceae species, as well as Arabidopsis thaliana and Oryza sativa. We identified 2,195 AP2/ERF genes, of which 368 (17%) were newly identified. Based on phylogenetic analyses, we observed expansion of the copy number of these genes, especially those belonging to specific Ethylene-Responsive Factor (ERF) subgroups of the Solanaceae. From the results of chromosomal location and synteny analyses, we identified that the AP2/ERF genes of the pepper (Capsicum annuum), the tomato (Solanum lycopersicum), and the potato (Solanum tuberosum) belonging to ERF subgroups form a tandem array and most of them are species-specific without orthologs in other species, which has led to differentiation of AP2/ERF gene repertory among Solanaceae. We suggest that these genes mainly emerged through recent gene duplication after the divergence of these species. Transcriptome analyses showed that the genes have a putative function in the response of the pepper and tomato to abiotic stress, especially those in ERF subgroups. CONCLUSIONS: Our findings will provide comprehensive information on AP2/ERF genes and insights into the structural, evolutionary, and functional understanding of the role of these genes in the Solanaceae.

Genetic diversity assessment and genome-wide association study reveal candidate genes associated with component traits in sweet potato (Ipomoea batatas (L.) Lam).

The Korean sweet potatoes were bred by various cultivars introduced from Japanese, American, Porto Rico, China, and Burundi. This issue enriched their genetic diversity but also resulted in a mixture of cultivars. For genotyping, we collected and sequenced 66 sweet potato germplasms from different localities around Korea, including 36 modern cultivars, 5 local cultivars, and 25 foreign cultivars. This identified 447.6 million trimmed reads and 324.8 million mapping reads and provided 39,424 single nucleotide polymorphisms (SNPs) markers. Phylogenetic clustering and population structure analysis distinctly classified these germplasms into 5 genetic groups, group 1, group 2, group 3, group 4, and group 5, containing 20, 15, 10, 7, and 14 accessions, respectively. Sixty-three significant SNPs were selected by genome-wide association for sugar composition-related traits (fructose, glucose, and total sugars), total starch, amylose content, and total carotenoid of the storage root. A total of 37 candidate genes encompassing these significant SNPs were identified, among which, 7 genes were annotated to involve in sugar and starch metabolism, including galactose metabolism (itf04g30630), starch and sucrose metabolism (itf03g13270, itf15g09320), carbohydrate metabolism (itf14g10250), carbohydrate and amino acid metabolism (itf12g19270), and amino sugar and nucleotide sugar metabolism (itf03g21950, itf15g04880). This results indicated that sugar and starch are important characteristics to determine the genetic diversity of sweet potatoes. These findings not only illustrate the importance of component traits to genotyping sweet potatoes but also explain an important reason resulting in genetic diversity of sweet potato.

Recurrent mutations promote widespread structural and functional divergence of MULE-derived genes in plants.

Transposable element (TE)-derived genes are increasingly recognized as major sources conferring essential traits in agriculturally important crops but underlying evolutionary mechanisms remain obscure. We updated previous annotations and constructed 18,744 FAR-RED IMPAIRED RESPONSE1 (FAR1) genes, a transcription factor family derived from Mutator-like elements (MULEs), from 80 plant species, including 15,546 genes omitted in previous annotations. In-depth sequence comparison of the updated gene repertoire revealed that FAR1 genes underwent continuous structural divergence via frameshift and nonsense mutations that caused premature translation termination or specific domain truncations. CRISPR/Cas9-based genome editing and transcriptome analysis determined a novel gene involved in fertility-regulating transcription of rice pollen, denoting the functional capacity of our re-annotated gene models especially in monocots which had the highest copy numbers. Genomic evidence showed that the functional gene adapted by obtaining a shortened form through a frameshift mutation caused by a tandem duplication of a 79-bp sequence resulting in premature translation termination. Our findings provide improved resources for comprehensive studies of FAR1 genes with beneficial agricultural traits and unveil novel evolutionary mechanisms generating structural divergence and subsequent adaptation of TE-derived genes in plants.

Identification of novel PHD-finger genes in pepper by genomic re-annotation and comparative analyses.

BACKGROUND: The plant homeodomain (PHD)-finger gene family that belongs to zinc-finger genes, plays an important role in epigenetics by regulating gene expression in eukaryotes. However, inaccurate annotation of PHD-finger genes hinders further downstream comparative, evolutionary, and functional studies. RESULTS: We performed genome-wide re-annotation in Arabidopsis thaliana (Arabidopsis), Oryza sativa (rice), Capsicum annuum (pepper), Solanum tuberosum (potato), and Solanum lycopersicum (tomato) to better understand the role of PHD-finger genes in these species. Our investigation identified 875 PHD-finger genes, of which 225 (26% of total) were newly identified, including 57 (54%) novel PHD-finger genes in pepper. The PHD-finger genes of the five plant species have various integrated domains that may be responsible for the diversification of structures and functions of these genes. Evolutionary analyses suggest that PHD-finger genes were expanded recently by lineage-specific duplication, especially in pepper and potato, resulting in diverse repertoires of PHD-finger genes among the species. We validated the expression of six newly identified PHD-finger genes in pepper with qRT-PCR. Transcriptome analyses suggest potential functions of PHD-finger genes in response to various abiotic stresses in pepper. CONCLUSIONS: Our data, including the updated annotation of PHD-finger genes, provide useful information for further evolutionary and functional analyses to better understand the roles of the PHD-finger gene family in pepper.

De novo phasing resolves haplotype sequences in complex plant genomes.

Genome phasing is a recently developed assembly method that separates heterozygous eukaryotic genomic regions and builds haplotype-resolved assemblies. Because differences between haplotypes are ignored in most published de novo genomes, assemblies are available as consensus genomes consisting of haplotype mixtures, thus increasing the need for genome phasing. Here, we review the operating principles and characteristics of several freely available and widely used phasing tools (TrioCanu, FALCON-Phase, and ALLHiC). An examination of downstream analyses using haplotype-resolved genome assemblies in plants indicated significant differences among haplotypes regarding chromosomal rearrangements, sequence insertions, and expression of specific alleles that contribute to the acquisition of the biological characteristics of plant species. Finally, we suggest directions to solve addressing limitations of current genome-phasing methods. This review provides insights into the current progress, limitations, and future directions of de novo genome phasing, which will enable researchers to easily access and utilize genome-phasing in studies involving highly heterozygous complex plant genomes.

Admixture of divergent genomes facilitates hybridization across species in the family Brassicaceae.

Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species are rare in nature due to reproductive barriers but how such hurdles can be overcome is largely unknown. Here we report the hybrid genome structure of xBrassicoraphanus, a synthetic allotetraploid of Brassica rapa and Raphanus sativus. We performed cytogenetic analysis and de novo genome assembly to examine chromosome behaviors and genome integrity in the hybrid. Transcriptome analysis was conducted to investigate expression of duplicated genes in conjunction with epigenome analysis to address whether genome admixture entails epigenetic reconfiguration. Allotetraploid xBrassicoraphanus retains both parental chromosomes without genome rearrangement. Meiotic synapsis formation and chromosome exchange are avoided between nonhomologous progenitor chromosomes. Reconfiguration of transcription network occurs, and less divergent cis-elements of duplicated genes are associated with convergent expression. Genome-wide DNA methylation asymmetry between progenitors is largely maintained but, notably, B. rapa-originated transposable elements are transcriptionally silenced in xBrassicoraphanus through gain of DNA methylation. Our results demonstrate that hybrid genome stabilization and transcription compatibility necessitate epigenome landscape adjustment and rewiring of cis-trans interactions. Overall, this study suggests that a certain extent of genome divergence facilitates hybridization across species, which may explain the great diversification and expansion of angiosperms during evolution.

Comparative analysis of de novo genomes reveals dynamic intra-species divergence of NLRs in pepper.

BACKGROUND: Peppers (Capsicum annuum L.) containing distinct capsaicinoids are the most widely cultivated spices in the world. However, extreme genomic diversity among species represents an obstacle to breeding pepper. RESULTS: Here, we report de novo genome assemblies of Capsicum annuum 'Early Calwonder (non-pungent, ECW)' and 'Small Fruit (pungent, SF)' along with their annotations. In total, we assembled 2.9 Gb of ECW and SF genome sequences, representing over 91% of the estimated genome sizes. Structural and functional annotation of the two pepper genomes generated about 35,000 protein-coding genes each, of which 93% were assigned putative functions. Comparison between newly and publicly available pepper gene annotations revealed both shared and specific gene content. In addition, a comprehensive analysis of nucleotide-binding and leucine-rich repeat (NLR) genes through whole-genome alignment identified five significant regions of NLR copy number variation (CNV). Detailed comparisons of those regions revealed that these CNVs were generated by intra-specific genomic variations that accelerated diversification of NLRs among peppers. CONCLUSIONS: Our analyses unveil an evolutionary mechanism responsible for generating CNVs of NLRs among pepper accessions, and provide novel genomic resources for functional genomics and molecular breeding of disease resistance in Capsicum species.

TGFam-Finder: a novel solution for target-gene family annotation in plants.

Whole-genome annotation error that omits essential protein-coding genes hinders further research. We developed Target Gene Family Finder (TGFam-Finder), an alternative tool for the structural annotation of protein-coding genes containing target domain(s) of interest in plant genomes. TGFam-Finder took considerably reduced annotation run-time and improved accuracy compared to conventional annotation tools. Large-scale re-annotation of 50 plant genomes identified an average of 150, 166 and 86 additional far-red-impaired response 1, nucleotide-binding and leucine-rich-repeat, and cytochrome P450 genes, respectively, that were missed in previous annotations. We detected significantly higher number of translated genes in the new annotations using mass spectrometry data from seven plant species compared to previous annotations. TGFam-Finder along with the new gene models can provide an optimized platform for comprehensive functional, comparative, and evolutionary studies in plants.

Genome and evolution of the shade-requiring medicinal herb Panax ginseng.

Panax ginseng C. A. Meyer, reputed as the king of medicinal herbs, has slow growth, long generation time, low seed production and complicated genome structure that hamper its study. Here, we unveil the genomic architecture of tetraploid P. ginseng by de novo genome assembly, representing 2.98 Gbp with 59 352 annotated genes. Resequencing data indicated that diploid Panax species diverged in association with global warming in Southern Asia, and two North American species evolved via two intercontinental migrations. Two whole genome duplications (WGD) occurred in the family Araliaceae (including Panax) after divergence with the Apiaceae, the more recent one contributing to the ability of P. ginseng to overwinter, enabling it to spread broadly through the Northern Hemisphere. Functional and evolutionary analyses suggest that production of pharmacologically important dammarane-type ginsenosides originated in Panax and are produced largely in shoot tissues and transported to roots; that newly evolved P. ginseng fatty acid desaturases increase freezing tolerance; and that unprecedented retention of chlorophyll a/b binding protein genes enables efficient photosynthesis under low light. A genome-scale metabolic network provides a holistic view of Panax ginsenoside biosynthesis. This study provides valuable resources for improving medicinal values of ginseng either through genomics-assisted breeding or metabolic engineering.

Expansion of sesquiterpene biosynthetic gene clusters in pepper confers nonhost resistance to the Irish potato famine pathogen.

Chemical barriers contribute to nonhost resistance, which is defined as the resistance of an entire plant species to nonadapted pathogen species. However, the molecular basis of metabolic defense in nonhost resistance remains elusive. Here, we report genetic evidence for the essential role of phytoalexin capsidiol in nonhost resistance of pepper (Capsicum spp.) to potato late blight Phytophthora infestans using transcriptome and genome analyses. Two different genes for capsidiol biosynthesis, 5-epi-aristolochene synthase (EAS) and 5-epi-aristolochene-1,3-dihydroxylase (EAH), belong to multigene families. However, only a subset of EAS/EAH gene family members were highly induced upon P. infestans infection, which was associated with parallel accumulation of capsidiol in P. infestans-infected pepper. Silencing of EAS homologs in pepper resulted in a significant decrease in capsidiol accumulation and allowed the growth of nonadapted P. infestans that is highly sensitive to capsidiol. Phylogenetic and genomic analyses of EAS/EAH multigene families revealed that the emergence of pathogen-inducible EAS/EAH genes in Capsicum-specific genomic regions rendered pepper a nonhost of P. infestans. This study provides insights into evolutionary aspects of nonhost resistance based on the combination of a species-specific phytoalexin and sensitivity of nonadapted pathogens.

Divergent evolution of multiple virus-resistance genes from a progenitor in Capsicum spp.

Plants have evolved hundreds of nucleotide-binding and leucine-rich domain proteins (NLRs) as potential intracellular immune receptors, but the evolutionary mechanism leading to the ability to recognize specific pathogen effectors is elusive. Here, we cloned Pvr4 (a Potyvirus resistance gene in Capsicum annuum) and Tsw (a Tomato spotted wilt virus resistance gene in Capsicum chinense) via a genome-based approach using independent segregating populations. The genes both encode typical NLRs and are located at the same locus on pepper chromosome 10. Despite the fact that these two genes recognize completely different viral effectors, the genomic structures and coding sequences of the two genes are strikingly similar. Phylogenetic studies revealed that these two immune receptors diverged from a progenitor gene of a common ancestor. Our results suggest that sequence variations caused by gene duplication and neofunctionalization may underlie the evolution of the ability to specifically recognize different effectors. These findings thereby provide insight into the divergent evolution of plant immune receptors.

RNA-seq pinpoints a Xanthomonas TAL-effector activated resistance gene in a large-crop genome.

Transcription activator-like effector (TALE) proteins of the plant pathogenic bacterial genus Xanthomonas bind to and transcriptionally activate host susceptibility genes, promoting disease. Plant immune systems have taken advantage of this mechanism by evolving TALE binding sites upstream of resistance (R) genes. For example, the pepper Bs3 and rice Xa27 genes are hypersensitive reaction plant R genes that are transcriptionally activated by corresponding TALEs. Both R genes have a hallmark expression pattern in which their transcripts are detectable only in the presence and not the absence of the corresponding TALE. By transcriptome profiling using next-generation sequencing (RNA-seq), we tested whether we could avoid laborious positional cloning for the isolation of TALE-induced R genes. In a proof-of-principle experiment, RNA-seq was used to identify a candidate for Bs4C, an R gene from pepper that mediates recognition of the Xanthomonas TALE protein AvrBs4. We identified one major Bs4C candidate transcript by RNA-seq that was expressed exclusively in the presence of AvrBs4. Complementation studies confirmed that the candidate corresponds to the Bs4C gene and that an AvrBs4 binding site in the Bs4C promoter directs its transcriptional activation. Comparison of Bs4C with a nonfunctional allele that is unable to recognize AvrBs4 revealed a 2-bp polymorphism within the TALE binding site of the Bs4C promoter. Bs4C encodes a structurally unique R protein and Bs4C-like genes that are present in many solanaceous genomes seem to be as tightly regulated as pepper Bs4C. These findings demonstrate that TALE-specific R genes can be cloned from large-genome crops with a highly efficient RNA-seq approach.

Haplotype-resolved genome assembly and resequencing analysis provide insights into genome evolution and allelic imbalance in Pinus densiflora.

Haplotype-level allelic characterization facilitates research on the functional, evolutionary and breeding-related features of extremely large and complex plant genomes. We report a 21.7-Gb chromosome-level haplotype-resolved assembly in Pinus densiflora. We found genome rearrangements involving translocations and inversions between chromosomes 1 and 3 of Pinus species and a proliferation of specific long terminal repeat (LTR) retrotransposons (LTR-RTs) in P. densiflora. Evolutionary analyses illustrated that tandem and LTR-RT-mediated duplications led to an increment of transcription factor (TF) genes in P. densiflora. The haplotype sequence comparison showed allelic imbalances, including presence-absence variations of genes (PAV genes) and their functional contributions to flowering and abiotic stress-related traits in P. densiflora. Allele-aware resequencing analysis revealed PAV gene diversity across P. densiflora accessions. Our study provides insights into key mechanisms underlying the evolution of genome structure, LTR-RTs and TFs within the Pinus lineage as well as allelic imbalances and diversity across P. densiflora.

Evolution of the large genome in Capsicum annuum occurred through accumulation of single-type long terminal repeat retrotransposons and their derivatives.

Although plant genome sizes are extremely diverse, the mechanism underlying the expansion of huge genomes that did not experience whole-genome duplication has not been elucidated. The pepper, Capsicum annuum, is an excellent model for studies of genome expansion due to its large genome size (2700 Mb) and the absence of whole genome duplication. As most of the pepper genome structure has been identified as constitutive heterochromatin, we investigated the evolution of this region in detail. Our findings show that the constitutive heterochromatin in pepper was actively expanded 20.0-7.5 million years ago through a massive accumulation of single-type Ty3/Gypsy-like elements that belong to the Del subgroup. Interestingly, derivatives of the Del elements, such as non-autonomous long terminal repeat retrotransposons and long-unit tandem repeats, played important roles in the expansion of constitutive heterochromatic regions. This expansion occurred not only in the existing heterochromatic regions but also into the euchromatic regions. Furthermore, our results revealed a repeat of unit length 18-24 kb. This repeat was found not only in the pepper genome but also in the other solanaceous species, such as potato and tomato. These results represent a characteristic mechanism for large genome evolution in plants.

Two different domain architectures generate structural and functional diversity among bZIP genes in the Solanaceae family.

The bZIP gene family is one of the largest transcription factor families and has important roles in plant growth, development, and stress responses. However, bZIP genes in the Solanaceae family have not been extensively investigated. Here, we conducted genome-wide re-annotation in nine Solanaceae species and Arabidopsis thaliana. We annotated 935 bZIP genes, including 107 (11%) that were newly identified. Structural analyses of bZIP genes in the Solanaceae family revealed that the bZIP domain displayed two types of architectures depending on the presence of an additional domain, suggesting that these architectures generate diversified structures and functions. Motif analyses indicated that the two types of bZIP genes had distinct sequences adjacent to the bZIP domain. Phylogenetic analyses suggested that the two types of bZIP genes distinctly evolved and ultimately adapted in different lineages. Transcriptome analyses in pepper (Capsicum annuum) and tomato (Solanum lycopersicum) revealed putative functional diversity between the two types of bZIP genes in response to various abiotic stresses. This study extensively updated bZIP gene family annotations and provided novel evolutionary and functional evidence for the role of bZIP genes in Solanaceae plants. Our findings provide evolutionary and functional characteristics of bZIP genes for a better understanding of their roles in Solanaceae plants.

Comparative analysis of pepper and tomato reveals euchromatin expansion of pepper genome caused by differential accumulation of Ty3/Gypsy-like elements.

BACKGROUND: Among the Solanaceae plants, the pepper genome is three times larger than that of tomato. Although the gene repertoire and gene order of both species are well conserved, the cause of the genome-size difference is not known. To determine the causes for the expansion of pepper euchromatic regions, we compared the pepper genome to that of tomato. RESULTS: For sequence-level analysis, we generated 35.6 Mb of pepper genomic sequences from euchromatin enriched 1,245 pepper BAC clones. The comparative analysis of orthologous gene-rich regions between both species revealed insertion of transposons exclusively in the pepper sequences, maintaining the gene order and content. The most common type of the transposon found was the LTR retrotransposon. Phylogenetic comparison of the LTR retrotransposons revealed that two groups of Ty3/Gypsy-like elements (Tat and Athila) were overly accumulated in the pepper genome. The FISH analysis of the pepper Tat elements showed a random distribution in heterochromatic and euchromatic regions, whereas the tomato Tat elements showed heterochromatin-preferential accumulation. CONCLUSIONS: Compared to tomato pepper euchromatin doubled its size by differential accumulation of a specific group of Ty3/Gypsy-like elements. Our results could provide an insight on the mechanism of genome evolution in the Solanaceae family.

Trichome density of main stem is tightly linked to PepMoV resistance in chili pepper (Capsicum annuum L.).

A relationship between pepper trichome and pepper mottle virus (PepMoV) resistance was examined. In an intraspecific F(2) mapping population from the cross between Capsicum annuum CM334 (trichome-bearing and PepMoV resistant) and Chilsungcho (glabrous and PepMoV susceptible), major QTLs for both traits were identified by composite interval mapping in linkage group (LG) 24 corresponding a telomere region on pepper chromosome 10. Ptel1 of putative trichome enhancing locus was a common major QTL for trichome density on the main stem and calyx. Ptel1 apart from HpmsE031 at a 1.03 cM interval was specifically associated to the trichome density on the main stem, whereas Ptel2 near m104 marker on LG2 was specific for the calyx trichome. Epistatic analysis indicated that Ptel1 engaged in controlling the trichome density by mutual interactions with the organ-specific QTLs. For PepMoV resistance, two QTLs (Pep1 and Pep2) were identified on the LG 24. Pep1 was located with Ptel1 in the R-gene cluster (RGC) for potyvirus resistance including Pvr4 with broad spectrum resistance to potyviruses. Pep1 flanking TG420 marker seemed to be the major factors determining correlation with PepMoV resistance. These results indicate that the level of trichome density on pepper main stem can be used as a morphological marker for Pvr4 in pepper breeding.

Comparative Analysis of Re-Annotated Genes Provides Insight into Evolutionary Divergence and Expressions of Aquaporin Family in Pepper.

Aquaporins (AQPs) are known to have a vital role in water transport in all living organisms including agriculturally important crops, but a comprehensive genomic study of AQPs in pepper has not been implemented. Here, we updated previous gene annotations and generated a total of 259 AQP genes from five plants, including pepper. Phylogenetic and motif analyses revealed that a large proportion of pepper AQP genes belong to the specific subgroup of tonoplast intrinsic protein (TIP) subfamily, TIP4. Chromosomal localization and estimated duplication times illustrated that genes in TIP4 formed a tandem array on the short arm of chromosome 1, resulting from pepper-specific expansion after its divergence with Solanaceae species. Transcriptome analyses under various abiotic stress conditions revealed that transport-, photosystem-, and thylakoid-related genes were generally enriched in expression clusters containing AQP genes in pepper. These results provide valuable genomic resources and insight into the evolutionary mechanism that generate genomic diversity of the AQP gene family in pepper.