Signature of diversifying selection on members of the pentatricopeptide repeat protein family in Arabidopsis lyrata

Adaptive evolution driving the young duplications in six Rosaceae species

Background

In plant genomes, high proportions of duplicate copies reveals that gene duplications play an important role in the evolutionary processes of plant species. A series of gene families under positive selection after recent duplication events in plant genomes indicated the evolution of duplicates driven by adaptive evolution. However, the genome-wide evolutionary features of young duplicate genes among closely related species are rarely reported.

Results

In this study, we conducted a systematic survey of young duplicate genes at genome-wide levels among six Rosaceae species, whose whole-genome sequencing data were successively released in recent years. A total of 35,936 gene families were detected among the six species, in which 60.25% were generated by young duplications. The 21,650 young duplicate gene families could be divided into two expansion types based on their duplication patterns, species-specific and lineage-specific expansions. Our results showed the species-specific expansions advantaging over the lineage-specific expansions. In the two types of expansions, high-frequency duplicate domains exhibited functional preference in response to environmental stresses.

Conclusions

The functional preference of the young duplicate genes in both the expansion types showed that they were inclined to respond to abiotic or biotic stimuli. Moreover, young duplicate genes under positive selection in both species-specific and lineage-specific expansions suggested that they were generated to adapt to the environmental factors in Rosaceae species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07422-7.

Association Mapping of Verticillium Wilt Disease in a Worldwide Collection of Cotton (Gossypium hirsutum L.)

Cotton (Gossypium spp.) is the best plant fiber source in the world and provides the raw material for industry. Verticillium wilt caused by Verticillium dahliae Kleb. is accepted as a major disease of cotton production. The most practical way to deal with verticillium wilt is to develop resistant/tolerant varieties after cultural practices. One of the effective selections in plant breeding is the use of marker-assisted selection (MAS) via quantitative trait loci (QTL). Therefore, in this study, we aimed to discover the genetic markers associated with the disease. Through the association mapping analysis, common single nucleotide polymorphism (SNP) markers were obtained using 4730 SNP alleles. As a result, twenty-three markers were associated with defoliating (PYDV6 isolate) pathotype, twenty-one markers with non-defoliating (Vd11 isolate) pathotype, ten QTL with Disease Severity Index (DSI) of the leaves at the 50–60% boll opening period and eight markers were associated with DSI in the stem section. Some of the markers that show significant associations are located on protein coding genes such as protein Mpv17-like, 21 kDa protein-like, transcription factor MYB113-like, protein dehydration-induced 19 homolog 3-like, F-box protein CPR30-like, extracellular ribonuclease LE-like, putative E3 ubiquitin-protein ligase LIN, pentatricopeptide repeat-containing protein At3g62890-like, fructose-1,6-bisphosphatase, tubby-like F-box protein 8, endoglucanase 16-like, glucose-6-phosphate/phosphate translocator 2, metal tolerance protein 11-like, VAN3-binding protein-like, transformation/transcription domain-associated protein-like, pyruvate kinase isozyme A, ethylene-responsive transcription factor CRF2-like, molybdate transporter 2-like, IRK-interacting protein-like, glycosylphosphatidylinositol anchor attachment 1 protein, U3 small nucleolar RNA-associated protein 4-like, microtubule-associated protein futsch-like, transport and Golgi organization 2 homolog, splicing factor 3B subunit 3-like, mediator of RNA polymerase II transcription subunit 15a-like, putative ankyrin repeat protein, and protein networked 1D-like. It has been reported in previous studies that most of these genes are associated with biotic and abiotic stress factors. As a result, once validated, it would be possible to use the markers obtained in the study in Marker Assisted Selection (MAS) breeding.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

Mitochondrial-nuclear co-evolution leads to hybrid incompatibility through pentatricopeptide repeat proteins

Mitochondrial-nuclear incompatibility has a major role in reproductive isolation between species. However, the underlying mechanism and driving force of mitochondrial-nuclear incompatibility remain elusive. Here, we report a pentatricopeptide repeat-containing (PPR) protein, Ccm1, and its interacting partner, 15S rRNA, to be involved in hybrid incompatibility between two yeast species, Saccharomyces cerevisiae and Saccharomyces bayanus S. bayanus-Ccm1 has reduced binding affinity for S. cerevisiae-15S rRNA, leading to respiratory defects in hybrid cells. This incompatibility can be rescued by single mutations on several individual PPR motifs, demonstrating the highly evolvable nature of PPR proteins. When we examined other PPR proteins in the closely related Saccharomyces sensu stricto yeasts, about two-thirds of them showed detectable incompatibility. Our results suggest that fast co-evolution between flexible PPR proteins and their mitochondrial RNA substrates may be a common driving force in the development of mitochondrial-nuclear hybrid incompatibility.

Comparative genomic analysis of the compound Brassica napus Rf locus

Background

The plant trait of cytoplasmically-inherited male sterility (CMS) and its suppression by nuclear restorer-of-fertility (Rf) genes can be viewed as a genetic arms race between the mitochondrial and nuclear genomes. Most nuclear Rf genes have been shown to encode P-type pentatricopeptide repeat proteins (PPRs). Phylogenetic analysis of P-class PPRs from sequenced plants genomes has shown that Rf-proteins cluster in a distinct clade of P-class PPRs, RFL-PPRs, that display hallmarks of positive evolutionary selection. Genes encoding RFL-PPRs (RFLs) within a given plant genome tend to be closely related both in sequence and position, but a detailed understanding of how such species-specific expansion occurs is lacking. In the canola, (oilseed rape) species Brassica napus, previous work has indicated the nuclear restorer genes for the two native forms of CMS, Rfn (for nap CMS) and Rfp (pol CMS), represent alternate haplotypes, or alleles, of a single nuclear locus.

Results

Fine genetic mapping indicates that Rfn does indeed localize to the same genomic region as Rfp. We find this region is enriched in RFL genes, three of which, based on their position and expression, represent potential candidates for Rfn; one of these genes, designated PPR4, is a preferred candidate in that it is not expressed in the nap CMS line. Comparison of the corresponding regions of the genomes of B. rapa, B. oleracea, Arabidopsis thaliana and A. lyrata provides insight into the expansion of this group of RFL genes in different lines of evolutionary descent.

Conclusions

Unlike other nuclear restorer loci containing multiple RFL genes, the RFL genes in the Rf region of B. napus are not present in tandem arrays but rather are dispersed in genomic location. The genes do not share similar flanking non-coding regions and do not contain introns, indicating that they have duplicated primarily through a retrotransposition-mediated process. In contrast, segmental duplication has been responsible for the distribution of the 10 sequences we annotated as RFL genes in the corresponding region of the A. lyrata genome. Our observations define the Brassica Rf locus and indicate that different mechanisms may be responsible for the proliferation of RFL genes even among closely related genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3117-0) contains supplementary material, which is available to authorized users.

Genomic Conflicts that Cause Pollen Mortality and Raise Reproductive Barriers in Arabidopsis thaliana

Species differentiation and the underlying genetics of reproductive isolation are central topics in evolutionary biology. Hybrid sterility is one kind of reproductive barrier that can lead to differentiation between species. Here, we analyze the complex genetic basis of the intraspecific hybrid male sterility that occurs in the offspring of two distant natural strains of Arabidopsis thaliana, Shahdara and Mr-0, with Shahdara as the female parent. Using both classical and quantitative genetic approaches as well as cytological observation of pollen viability, we demonstrate that this particular hybrid sterility results from two causes of pollen mortality. First, the Shahdara cytoplasm induces gametophytic cytoplasmic male sterility (CMS) controlled by several nuclear loci. Second, several segregation distorters leading to allele-specific pollen abortion (pollen killers) operate in hybrids with either cytoplasm. The complete sterility of the hybrid with the Shahdara cytoplasm results from the genetic linkage of the two causes of pollen mortality, i.e., CMS nuclear determinants and pollen killers. Furthermore, natural variation at these loci in A. thaliana is associated with different male-sterility phenotypes in intraspecific hybrids. Our results suggest that the genomic conflicts that underlie segregation distorters and CMS can concurrently lead to reproductive barriers between distant strains within a species. This study provides a new framework for identifying molecular mechanisms and the evolutionary history of loci that contribute to reproductive isolation, and possibly to speciation. It also suggests that two types of genomic conflicts, CMS and segregation distorters, may coevolve in natural populations.

Genome-Wide Association Study on Resistance to Stalk Rot Diseases in Grain Sorghum

Stalk rots are important biotic constraints to sorghum production worldwide. Several pathogens may be associated with the disease, but Macrophomina phaseolina and Fusarium thapsinum are recognized as the major causal organisms. The diseases become more aggressive when drought and high-temperature stress occur during grain filling. Progress in genetic improvement efforts has been slow due to lack of effective phenotyping protocol and the strong environmental effect on disease incidence and severity. Deployment of modern molecular tools is expected to accelerate efforts to develop resistant hybrids. This study was aimed at identifying genomic regions associated with resistance to both causal organisms. A sorghum diversity panel consisting of 300 genotypes assembled from different parts of the world was evaluated for response to infection by both pathogens. Community resources of 79,132 single nucleotide polymorphic (SNP) markers developed on the panel were used in association studies using a multi-locus mixed model to map loci associated with stalk rot resistance. Adequate genetic variation was observed for resistance to both pathogens. Structure analysis grouped the genotypes into five subpopulations primarily based on the racial category of the genotypes. Fourteen loci and a set of candidate genes appear to be involved in connected functions controlling plant defense response. However, each associated SNP had relatively small effect on the traits, accounting for 19-30% of phenotypic variation. Linkage disequilibrium analyses suggest that significant SNPs are genetically independent. Estimation of frequencies of associated alleles revealed that durra and caudatum subpopulations were enriched for resistant alleles, but the results suggest complex molecular mechanisms underlying resistance to both pathogens.

Detecting selection in the blue crab, Callinectes sapidus, using DNA sequence data from multiple nuclear protein-coding genes

The identification of genes involved in the adaptive evolution of non-model organisms with uncharacterized genomes constitutes a major challenge. This study employed a rigorous and targeted candidate gene approach to test for positive selection on protein-coding genes of the blue crab, Callinectes sapidus. Four genes with putative roles in physiological adaptation to environmental stress were chosen as candidates. A fifth gene not expected to play a role in environmental adaptation was used as a control. Large samples (n>800) of DNA sequences from C. sapidus were used in tests of selective neutrality based on sequence polymorphisms. In combination with these, sequences from the congener C. similis were used in neutrality tests based on interspecific divergence. In multiple tests, significant departures from neutral expectations and indicative of positive selection were found for the candidate gene trehalose 6-phosphate synthase (tps). These departures could not be explained by any of the historical population expansion or bottleneck scenarios that were evaluated in coalescent simulations. Evidence was also found for balancing selection at ATP-synthase subunit 9 (atps) using a maximum likelihood version of the Hudson, Kreitmen, and Aguadé test, and positive selection favoring amino acid replacements within ATP/ADP translocase (ant) was detected using the McDonald-Kreitman test. In contrast, test statistics for the control gene, ribosomal protein L12 (rpl), which presumably has experienced the same demographic effects as the candidate loci, were not significantly different from neutral expectations and could readily be explained by demographic effects. Together, these findings demonstrate the utility of the candidate gene approach for investigating adaptation at the molecular level in a marine invertebrate for which extensive genomic resources are not available.

Cytoplasmic male sterility contributes to hybrid incompatibility between subspecies of Arabidopsis lyrata

In crosses between evolutionarily diverged populations, genomic incompatibilities may result in sterile hybrids, indicating evolution of reproductive isolation. In several plant families, crosses within a population can also lead to male sterile progeny because of conflict between the maternally and biparentally inherited genomes. We examined hybrid fertility between subspecies of the perennial outcrossing self-incompatible Lyrate rockcress (Arabidopsis lyrata) in large reciprocal F2 progenies and three generations of backcrosses. In one of the reciprocal F2 progenies, almost one-fourth of the plants were male-sterile. Correspondingly, almost one-half of the plants in one of the four reciprocal backcross progenies expressed male sterility. In an additional four independent F2 and backcross families, three segregated male sterility. The observed asymmetrical hybrid incompatibility is attributable to male sterility factors in one cytoplasm, for which the other population lacks effective fertility restorers. Genotyping of 96 molecular markers and quantitative trait locus mapping revealed that only 60% of the plants having the male sterile cytoplasm and lacking the corresponding restorers were phenotypically male-sterile. Genotyping data showed that there is only one restorer locus, which mapped to a 600-kb interval at the top of chromosome 2 in a region containing a cluster of pentatricopeptide repeat genes. Male fertility showed no trade-off with seed production. We discuss the role of cytoplasm and genomic conflict in incipient speciation and conclude that cytoplasmic male sterility–lowering hybrid fitness is a transient effect with limited potential to form permanent reproductive barriers between diverged populations of hermaphrodite self-incompatible species.

The evolutionary ecology of cytonuclear interactions in angiosperms

Interactions between cytoplasmic and nuclear genomes have significant evolutionary consequences. In angiosperms, the most common cytonuclear interaction is between mitochondrial genes that disrupt pollen production (cytoplasmic male sterility, CMS) and nuclear genes that restore it (nuclear male fertility restorers, Rf). The outcome of CMS/Rf interactions can depend on whether Rf alleles have negative pleiotropic effects on fitness. Although these fitness costs are often considered to be independent of the ecological context, we argue that the effects of Rf alleles on fitness should be context dependent. Thus, measuring the cost of restoration across a range of environments could help explain geographic and phylogenetic variation in the distribution of Rf alleles and the outcome of CMS/Rf interactions.

Genome evolution in outcrossing versus selfing versus asexual species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes, and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. However, while asexual species essentially conform to theoretical predictions, the situation is less simple in selfing species. We discuss the possible reasons to potentially explain this paradox. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

Selection patterns on restorer-like genes reveal a conflict between nuclear and mitochondrial genomes throughout angiosperm evolution

Eukaryotic cells have harbored mitochondria for at least 1.5 billion years in an apparently mutually beneficial symbiosis. Studies on the agronomically important crop trait cytoplasmic male sterility (CMS) have suggested the semblance of a host-parasite relationship between the nuclear and mitochondrial genomes, but molecular evidence for this is lacking. Key players in CMS systems are the fertility restorer (Rf) genes required for the development of a functional male gametophyte in plants carrying a mitochondrial CMS gene. In the majority of cases, Rf genes encode pentatricopeptide repeat (PPR) proteins. We show that most angiosperms for which extensive genomic sequence data exist contain multiple PPR genes related to Rf genes. These Rf-like genes show a number of characteristic features compared with other PPR genes, including chromosomal clustering and unique patterns of evolution, notably high rates of nonsynonymous to synonymous substitutions, suggesting diversifying selection. The highest probabilities of diversifying selection were seen for amino acid residues 1, 3, and 6 within the PPR motif. PPR proteins are involved in RNA processing, and mapping the selection data to a predicted consensus structure of an array of PPR motifs suggests that these residues are likely to form base-specific contacts to the RNA ligand. We suggest that the selection patterns on Rf-like genes reveal a molecular "arms-race" between the nuclear and mitochondrial genomes that has persisted throughout most of the evolutionary history of angiosperms.

Speciation genes in plants

BACKGROUND

Analyses of speciation genes--genes that contribute to the cessation of gene flow between populations--can offer clues regarding the ecological settings, evolutionary forces and molecular mechanisms that drive the divergence of populations and species. This review discusses the identities and attributes of genes that contribute to reproductive isolation (RI) in plants, compares them with animal speciation genes and investigates what these genes can tell us about speciation.

SCOPE

Forty-one candidate speciation genes were identified in the plant literature. Of these, seven contributed to pre-pollination RI, one to post-pollination, prezygotic RI, eight to hybrid inviability, and 25 to hybrid sterility. Genes, gene families and genetic pathways that were frequently found to underlie the evolution of RI in different plant groups include the anthocyanin pathway and its regulators (pollinator isolation), S RNase-SI genes (unilateral incompatibility), disease resistance genes (hybrid necrosis), chimeric mitochondrial genes (cytoplasmic male sterility), and pentatricopeptide repeat family genes (cytoplasmic male sterility).

CONCLUSIONS

The most surprising conclusion from this review is that identities of genes underlying both prezygotic and postzygotic RI are often predictable in a broad sense from the phenotype of the reproductive barrier. Regulatory changes (both cis and trans) dominate the evolution of pre-pollination RI in plants, whereas a mix of regulatory mutations and changes in protein-coding genes underlie intrinsic postzygotic barriers. Also, loss-of-function mutations and copy number variation frequently contribute to RI. Although direct evidence of positive selection on speciation genes is surprisingly scarce in plants, analyses of gene family evolution, along with theoretical considerations, imply an important role for diversifying selection and genetic conflict in the evolution of RI. Unlike in animals, however, most candidate speciation genes in plants exhibit intraspecific polymorphism, consistent with an important role for stochastic forces and/or balancing selection in development of RI in plants.