Relaxation Of Selective Constraints Research Articles

Purifying Selection Influences the Comparison of Heterozygosities between Populations.

Heterozygosity is a fundamental measure routinely used to compare between populations to infer the level of genetic variation and their relative effective population sizes. However, such comparison is highly influenced by the magnitude of selection pressure on the genomic regions used. Using over 2 million Single Nucleotide Variants (SNVs) from chimpanzee and mouse populations, this study shows that the heterozygosities estimated using neutrally evolving sites of large populations were two times higher than those of small populations. However, this difference was only ~1.6 times for the heterozygosities estimated using nonsynonymous sites. This suggests an excess in the nonsynonymous heterozygosities due to the segregation of deleterious variants in small populations. This excess in the nonsynonymous heterozygosities of the small populations was estimated to be 23-31%. Further analysis revealed that the magnitude of the excess is modulated by effective population size (Ne) and selection intensity (s). Using chimpanzee populations, this investigation found that the excess in nonsynonymous diversity in the small population was little (6%) when the difference between the Ne values of large and small populations was small (2.4 times). Conversely, this was high (23%) when the difference in Ne was large (5.9 times). Analysis using mouse populations showed that the excess in the nonsynonymous diversity of highly constrained genes of the small population was much higher (38%) than that observed for the genes under relaxed selective constraints (21%). Similar results were observed when the expression levels of genes were used as a proxy for selection intensity. These results emphasize the use of neutral regions, less constrained genes, or lowly expressed genes when comparing the heterozygosities between populations.

Mitochondrial DNA Genomes Reveal Relaxed Purifying Selection During Human Population Expansion after the Last Glacial Maximum.

Modern humans have experienced explosive population growth in the past thousand years. We hypothesized that recent human populations have inhabited environments with relaxation of selective constraints, possibly due to the more abundant food supply after the Last Glacial Maximum. The ratio of nonsynonymous to synonymous mutations (N/S ratio) is a useful and common statistic for measuring selective constraints. In this study, we reconstructed a high-resolution phylogenetic tree using a total of 26,419 East Eurasian mitochondrial DNA genomes, which were further classified into expansion and nonexpansion groups on the basis of the frequencies of their founder lineages. We observed a much higher N/S ratio in the expansion group, especially for nonsynonymous mutations with moderately deleterious effects, indicating a weaker effect of purifying selection in the expanded clades. However, this observation on N/S ratio was unlikely in computer simulations where all individuals were under the same selective constraints. Thus, we argue that the expanded populations were subjected to weaker selective constraints than the nonexpanded populations were. The mildly deleterious mutations were retained during population expansion, which could have a profound impact on present-day disease patterns.

Functional and regulatory diversification of Period genes responsible for circadian rhythm in vertebrates.

The Period genes (Per) play essential roles in modulating the molecular circadian clock timing in a broad range of species, which regulates the physiological and cellular rhythms through the transcription-translation feedback loop. While the Period gene paralogs are widely observed among vertebrates, the evolutionary history and the functional diversification of Per genes across vertebrates are not well known. In this study, we comprehensively investigated the evolution of Per genes at the copy number and sequence levels, including de novo binding motif discovery by comparative genomics. We also determined the lineage-specific transcriptome landscape across tissues and developmental stages and phenotypic effects in public RNA-seq data sets of model species. We observed multiple lineage-specific gain and loss events Per genes, though no simple association was observed between ecological factors and Per gene numbers in each species. Among salmonid fish species, the per3 gene has been lost in the majority, whereas those retaining the per3 gene exhibit not a signature of relaxed selective constraint but rather a signature of intensified selection. We also determined the signature of adaptive diversification of the CRY-binding region in Per1 and Per3, which modulates the circadian rhythm. We also discovered putative regulatory sequences, which are lineage-specific, suggesting that these cis-regulatory elements may have evolved rapidly and divergently across different lineages. Collectively, our findings revealed the evolution of Per genes and their fine-tuned contribution to the plastic and precise regulation of circadian rhythms in various vertebrate taxa.

The evolution of seminal fluid gene expression and postmating reproductive isolation in Drosophila.

Seminal fluid protein (Sfp) genes show, in general, a higher rate of sequence divergence than genes from other categories, which is often attributed to forms of postcopulatory sexual selection or sexual conflict. Recently, the relaxation of selective constraints has been proposed as an alternative explanation for the rapid sequence evolution of Sfps and other genes with sex-limited expression. The expression of Sfp genes is a likely target of selection, but the evolution of differences in their expression levels is less understood. Here, we explore both polymorphism and divergence in Sfp gene expression between Drosophila melanogaster and Drosophila simulans, how selection might have influenced their expression, and whether changes in expression might trigger the evolution of reproductive isolating barriers. In our analysis, Sfp genes showed higher divergence, but not higher polymorphism, in expression than the average male reproductive glands gene. Sfp genes with reproductive-tissue-specific expression were enriched for both directional and stabilizing selection, while relaxed selection was the predominant mode of evolution among Sfp genes with any other nonreproductive tissue-specific or nontissue-specific expression. The knockdown of single genes known to affect intraspecific sperm competition, and with patterns of expression divergence and polymorphism suggestive of directional selection, was not enough to break down postmating reproductive isolation barriers between species. Our results identify the expression of male-specific Sfp genes as an enriched target of selection and suggest a complex molecular relationship between postcopulatory sexual selection on a single gene's expression and its effect on the onset of speciation.

Mitochondrial phylogenomics of bumblebees, Bombus (Hymenoptera: Apidae): a tale of structural variation, shifts in selection constraints, and tree discordance

Abstract The mitochondrial DNA (mtDNA) of bumblebees (Bombus) has been widely used for phylogenetic studies, but its evolution is still underexplored. Here we report a comprehensive analysis of 60 bumblebee mitogenomes, including 40 newly assembled ones, to investigate bumblebee mtDNA structure, composition, and informativeness under a phylogenetic framework. Our mtDNA dataset supports the monophyly of Bombus and its subgenera, although we found a high degree of tree discordance in deeper nodes when using different inference methods or matrix composition. Concerning mitogenome structure, our results show that tRNA genes were often rearranged, with unique rearrangements indicating shared ancestry across bumblebee subgenera, illustrating their potential for subgeneric classification. Our results also challenge the notion that faster evolving mtDNA exhibits higher gene rearrangement rates. Finally, we explicitly assessed shifts in selection constraints of mtDNA genes in obligate social parasites of subgenus Psithyrus and found that their mtDNA evolved under relaxed selective constraints. Our findings show the utility of mtDNA in providing insights into bumblebee phylogenetics, evolution, and genome trait diversification. We also highlight the potential for comparative mitogenomics to uncover previously unknown aspects of bumblebee evolution, offering exciting opportunities for future research in this field.

Comparative single-cell transcriptomic analysis of primate brains highlights human-specific regulatory evolution.

Enhanced cognitive function in humans is hypothesized to result from cortical expansion and increased cellular diversity. However, the mechanisms that drive these phenotypic innovations remain poorly understood, in part because of the lack of high-quality cellular resolution data in human and non-human primates. Here, we take advantage of single-cell expression data from the middle temporal gyrus of five primates (human, chimp, gorilla, macaque and marmoset) to identify 57 homologous cell types and generate cell type-specific gene co-expression networks for comparative analysis. Although orthologue expression patterns are generally well conserved, we find 24% of genes with extensive differences between human and non-human primates (3,383 out of 14,131), which are also associated with multiple brain disorders. To assess the functional significance of gene expression differences in an evolutionary context, we evaluate changes in network connectivity across meta-analytic co-expression networks from 19 animals. We find that a subset of these genes has deeply conserved co-expression across all non-human animals, and strongly divergent co-expression relationships in humans (139 out of 3,383, <1% of primate orthologues). Genes with human-specific cellular expression and co-expression profiles (such as NHEJ1, GTF2H2, C2 and BBS5) typically evolve under relaxed selective constraints and may drive rapid evolutionary change in brain function.

Parallel expansion and divergence of an adhesin family in pathogenic yeasts.

Opportunistic yeast pathogens arose multiple times in the Saccharomycetes class, including the recently emerged, multidrug-resistant (MDR) Candida auris. We show that homologs of a known yeast adhesin family in Candida albicans, the Hyr/Iff-like (Hil) family, are enriched in distinct clades of Candida species as a result of multiple, independent expansions. Following gene duplication, the tandem repeat-rich region in these proteins diverged extremely rapidly and generated large variations in length and β-aggregation potential, both of which are known to directly affect adhesion. The conserved N-terminal effector domain was predicted to adopt a β-helical fold followed by an α-crystallin domain, making it structurally similar to a group of unrelated bacterial adhesins. Evolutionary analyses of the effector domain in C. auris revealed relaxed selective constraint combined with signatures of positive selection, suggesting functional diversification after gene duplication. Lastly, we found the Hil family genes to be enriched at chromosomal ends, which likely contributed to their expansion via ectopic recombination and break-induced replication. Combined, these results suggest that the expansion and diversification of adhesin families generate variation in adhesion and virulence within and between species and are a key step toward the emergence of fungal pathogens.

Intronization Signatures in Coding Exons Reveal the Evolutionary Fluidity of Eukaryotic Gene Architecture.

The conventionally clear distinction between exons and introns in eukaryotic genes is actually blurred. To illustrate this point, consider sequences that are retained in mature mRNAs about 50% of the time: how should they be classified? Moreover, although it is clear that RNA splicing influences gene expression levels and is an integral part of interdependent cellular networks, introns continue to be regarded as accidental insertions; exogenous sequences whose evolutionary origin is independent of mRNA-associated processes and somewhat still elusive. Here, we present evidence that aids to resolve this disconnect between conventional views about introns and current knowledge about the role of RNA splicing in the eukaryotic cell. We first show that coding sequences flanked by cryptic splice sites are negatively selected on a genome-wide scale in Paramecium. Then, we exploit selection intensity to infer splicing-related evolutionary dynamics. Our analyses suggest that intron gain begins as a splicing error, involves a transient phase of alternative splicing, and is preferentially completed at the 5’ end of genes, which through intron gain can become highly expressed. We conclude that relaxed selective constraints may promote biological complexity in Paramecium and that the relationship between exons and introns is fluid on an evolutionary scale.

Mitochondrial Genomes of two Lycosa spiders (Araneae, Lycosidae): Genome Description and Phylogenetic Implications

We sequenced the complete mitochondrial genomes of Lycosa shansia, and Lycosa singoriensis by combining Sanger and next-generation sequencing methods and analyzed the sequenced genomes in order to explore the phylogenetic placement and the mitogenome composition and evolution of these species. The mitochondrial genome of L. shansia was 14,638 bp, whereas that of L. singoriensis was 13,686 bp. The type of genes and direction of the coding strand present in the mitogenomes were the same as those in other species of Lycosoidea, including two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 13 protein-coding genes (PCGs). The mitogenomes of the two species exhibited negative AT and positive GC skews. This indicated that the nucleotide compositions of the mitogenomes of L. singoriensis and L. shansia tended to be T and G. Both the mean and median values of Ka/Ks of ATP8 were the highest among the 13 protein-coding genes, indicating that it might have evolved more rapidly than the other protein-coding genes in both species. ATP8 may have undergone more relaxed selection constraints and accumulated more mutations. In addition, many tRNAs lacked T and D stem loops; a few had no acceptor stems. The assessed species were recovered nested within Lycosidae with high support. The present findings will be useful for future studies on the mitogenome evolution of spiders.

Decay of TRPV3 as the genomic trace of epidermal structure changes in the land-to-sea transition of mammals.

The epidermis plays an indispensable barrier function in animals. Some species have evolved unique epidermal structures to adapt to different environments. Aquatic and semi‐aquatic mammals (cetaceans, manatees, and hippopotamus) are good models to study the evolution of epidermal structures because of their exceptionally thickened stratum spinosum, the lack of stratum granulosum, and the parakeratotic stratum corneum. This study aimed to analyze an upstream regulatory gene transient receptor potential cation channel, subfamily V, member 3 (TRPV3) of epidermal differentiation so as to explore the association between TRPV3 evolution and epidermal changes in mammals. Inactivating mutations were detected in almost all the aquatic cetaceans and several terrestrial mammals. Relaxed selective pressure was examined in the cetacean lineages with inactivated TRPV3, which might contribute to its exceptionally thickened stratum spinosum as the significant thickening of stratum spinosum in TRPV3 knock‐out mouse. However, functional TRPV3 may exist in several terrestrial mammals due to their strong purifying selection, although they have “inactivating mutations.” Further, for intact sequences, relaxed selective constraints on the TRPV3 gene were also detected in aquatic cetaceans, manatees, and semi‐aquatic hippopotamus. However, they had intact TRPV3, suggesting that the accumulation of inactivating mutations might have lagged behind the relaxed selective pressure. The results of this study revealed the decay of TRPV3 being the genomic trace of epidermal development in aquatic and semi‐aquatic mammals. They provided insights into convergently evolutionary changes of epidermal structures during the transition from the terrestrial to the aquatic environment.

Noninsect-Based Diet Leads to Structural and Functional Changes of Acidic Chitinase in Carnivora.

Acidic chitinase (Chia) digests the chitin of insects in the omnivorous stomach and the chitinase activity in carnivorous Chia is significantly lower than that of the omnivorous enzyme. However, mechanistic and evolutionary insights into the functional changes in Chia remain unclear. Here we show that a noninsect-based diet has caused structural and functional changes in Chia during the course of evolution in Carnivora. By creating mouse-dog chimeric Chia proteins and modifying the amino acid sequences, we revealed that F214L and A216G substitutions led to the dog enzyme activation. In 31 Carnivora, Chia was present as a pseudogene with stop codons in the open reading frame (ORF) region. Importantly, the Chia proteins of skunk, meerkat, mongoose, and hyena, which are insect-eating species, showed high chitinolytic activity. The cat Chia pseudogene product was still inactive even after ORF restoration. However, the enzyme was activated by matching the number and position of Cys residues to an active form and by introducing five meerkat Chia residues. Mutations affecting the Chia conformation and activity after pseudogenization have accumulated in the common ancestor of Felidae due to functional constraints. Evolutionary analysis indicates that Chia genes are under relaxed selective constraint in species with noninsect-based diets except for Canidae. These results suggest that there are two types of inactivating processes in Carnivora and that dietary changes affect the structure and activity of Chia.

Multi-tissue transcriptome analysis of two Begonia species reveals dynamic patterns of evolution in the chalcone synthase gene family

Begonia is an important horticultural plant group, as well as one of the most speciose Angiosperm genera, with over 2000 described species. Genus wide studies of genome size have shown that Begonia has a highly variable genome size, and analysis of paralog pairs has previously suggested that Begonia underwent a whole genome duplication. We address the contribution of gene duplication to the generation of diversity in Begonia using a multi-tissue RNA-seq approach. We chose to focus on chalcone synthase (CHS), a gene family having been shown to be involved in biotic and abiotic stress responses in other plant species, in particular its importance in maximising the use of variable light levels in tropical plants. We used RNA-seq to sample six tissues across two closely related but ecologically and morphologically divergent species, Begonia conchifolia and B. plebeja, yielding 17,012 and 19,969 annotated unigenes respectively. We identified the chalcone synthase gene family members in our Begonia study species, as well as in Hillebrandia sandwicensis, the monotypic sister genus to Begonia, Cucumis sativus, Arabidopsis thaliana, and Zea mays. Phylogenetic analysis suggested the CHS gene family has high duplicate turnover, all members of CHS identified in Begonia arising recently, after the divergence of Begonia and Cucumis. Expression profiles were similar within orthologous pairs, but we saw high inter-ortholog expression variation. Sequence analysis showed relaxed selective constraints on some ortholog pairs, with substitutions at conserved sites. Evidence of pseudogenisation and species specific duplication indicate that lineage specific differences are already beginning to accumulate since the divergence of our study species. We conclude that there is evidence for a role of gene duplication in generating diversity through sequence and expression divergence in Begonia.

Relaxation of Selective Constraint on the Ultra-Large Mitochondrial Genomes of Arcidae (Mollusca: Bivalvia)

Relaxation of Selective Constraint on the Ultra-Large Mitochondrial Genomes of Arcidae (Mollusca: Bivalvia)

Relaxed selection underlies genome erosion in socially parasitic ant species

Inquiline ants are highly specialized and obligate social parasites that infiltrate and exploit colonies of closely related species. They have evolved many times convergently, are often evolutionarily young lineages, and are almost invariably rare. Focusing on the leaf-cutting ant genus Acromyrmex, we compared genomes of three inquiline social parasites with their free-living, closely-related hosts. The social parasite genomes show distinct signatures of erosion compared to the host lineages, as a consequence of relaxed selective constraints on traits associated with cooperative ant colony life and of inquilines having very small effective population sizes. We find parallel gene losses, particularly in olfactory receptors, consistent with inquiline species having highly reduced social behavioral repertoires. Many of the genomic changes that we uncover resemble those observed in the genomes of obligate non-social parasites and intracellular endosymbionts that branched off into highly specialized, host-dependent niches.

Transit From Autotrophism to Heterotrophism: Sequence Variation and Evolution of Chloroplast Genomes in Orobanchaceae Species.

The family Orobanchaceae including autotrophic, hemiparasitic, and holoparasitic species, is becoming a key taxa to study the evolution of chloroplast genomes in different lifestyles. But the early evolutionary trajectory in the transit from autotrophism to hemiparasitism still maintains unclear for the inadequate sampling. In this study, we compared 50 complete chloroplast genomes in Orobanchaceae, containing four newly sequenced plastomes from hemiparasitic Pedicularis, to elucidate the sequence variation patterns in the evolution of plastomes. Contrasted to the sequence and structural hypervariabilities in holoparasites, hemiparasitic plastomes exhibited high similarity to those of autotrophs in gene and GC contents. They are generally characterized with functional or physical loss of ndh/tRNA genes and the inverted small-single-copy region. Gene losses in Orobanchaceae were lineage-specific and convergent, possibly related to structural reconfiguration and expansion/contraction of the inverted region. Pseudogenization of ndh genes was unique in hemiparasites. At least in Pedicularis, the ndhF gene might be most sensitive to the environmental factors and easily pseudogenized when autotrophs transit to hemiparasites. And the changes in gene contents and structural variation potentially deeply rely on the feeding type. Selective pressure, together with mutational bias, was the dominant factor of shaping the codon usage patterns. The relaxed selective constraint, potentially with genome-based GC conversion (gBGC) and preferential codon usage, drive the fluctuation of GC contents among taxa with different lifestyles. Phylogenetic analysis in Orobanchaceae supported that parasitic species were single-originated while holoparasites were multiple-originated. Overall, the comparison of plastomes provided a good opportunity to understand the evolution process in Orobanchaceae with different lifestyles.

Disrupted architecture and fast evolution of the mitochondrial genome of Argeia pugettensis (Isopoda): implications for speciation and fitness

BackgroundArgeia pugettensis is an isopod species that parasitizes other crustaceans. Its huge native geographic range spans the Pacific from China to California, but molecular data are available only for a handful of specimens from North-American populations. We sequenced and characterised the complete mitogenome of a specimen collected in the Yellow Sea.ResultsIt exhibited a barcode (cox1) similarity level of only 87–89% with North-American populations, which is unusually low for conspecifics. Its mitogenome is among the largest in isopods (≈16.5 Kbp), mostly due to a large duplicated palindromic genomic segment (2 Kbp) comprising three genes. However, it lost a segment comprising three genes, nad4L-trnP-nad6, and many genes exhibited highly divergent sequences in comparison to isopod orthologues, including numerous mutations, deletions and insertions. Phylogenetic and selection analyses corroborated that this is one of the handful of most rapidly evolving available isopod mitogenomes, and that it evolves under highly relaxed selection constraints (as opposed to positive selection). However, its nuclear 18S gene is highly conserved, which suggests that rapid evolution is limited to its mitochondrial genome. The cox1 sequence analysis indicates that elevated mitogenomic evolutionary rates are not shared by North-American conspecifics, which suggests a breakdown of cox1 barcoding in this species.ConclusionsA highly architecturally disrupted mitogenome and decoupling of mitochondrial and nuclear rates would normally be expected to have strong negative impacts on the fitness of the organism, so the existence of this lineage is a puzzling evolutionary question. Additional studies are needed to assess the phylogenetic breadth of this disrupted mitochondrial architecture and its impact on fitness.

GBA3: a polymorphic pseudogene in humans that experienced repeated gene loss during mammalian evolution

The gene encoding the cytosolic β-glucosidase GBA3 shows pseudogenization due to a truncated allele (rs358231) that is polymorphic in humans. Since this enzyme is involved in the transformation of many plant β-glycosides, this particular case of gene loss may have been influenced by dietary adaptations during evolution. In humans, apart from the inactivating allele, we found that GBA3 accumulated additional damaging mutations, implying an extensive GBA3 loss. The allelic distribution of loss-of-function alleles revealed significant differences between human populations which can be partially related with their staple diet. The analysis of mammalian orthologs disclosed that GBA3 underwent at least nine pseudogenization events. Most events of pseudogenization occurred in carnivorous lineages, suggesting a possible link to a β-glycoside poor diet. However, GBA3 was also lost in omnivorous and herbivorous species, hinting that the physiological role of GBA3 is not fully understood and other unknown causes may underlie GBA3 pseudogenization. Such possibility relies upon a putative role in sialic acid biology, where GBA3 participates in a cellular network involving NEU2 and CMAH. Overall, our data shows that the recurrent loss of GBA3 in mammals is likely to represent an evolutionary endpoint of the relaxation of selective constraints triggered by diet-related factors.

Degradation of key photosynthetic genes in the critically endangered semi-aquatic flowering plant Saniculiphyllum guangxiense (Saxifragaceae)

BackgroundPlastid gene loss and pseudogenization has been widely documented in parasitic and mycoheterotrophic plants, which have relaxed selective constraints on photosynthetic function. More enigmatic are sporadic reports of pseudogenization and loss of important photosynthesis genes in lineages thought to be fully photosynthetic. Here we report the complete plastid genome of Saniculiphyllum guangxiense, a critically endangered and phylogenetically isolated plant lineage, along with genomic evidence of reduced chloroplast function. We also report 22 additional plastid genomes representing the diversity of its containing clade Saxifragales, characterizing gene content and placing variation in a broader phylogenetic context.ResultsWe find that the plastid genome of Saniculiphyllum has experienced pseudogenization of five genes of the ndh complex (ndhA, ndhB, ndhD, ndhF, and ndhK), previously reported in flowering plants with an aquatic habit, as well as the surprising pseudogenization of two genes more central to photosynthesis (ccsA and cemA), contrasting with strong phylogenetic conservatism of plastid gene content in all other sampled Saxifragales. These genes participate in photooxidative protection, cytochrome synthesis, and carbon uptake. Nuclear paralogs exist for all seven plastid pseudogenes, yet these are also unlikely to be functional.ConclusionsSaniculiphyllum appears to represent the greatest degree of plastid gene loss observed to date in any fully photosynthetic lineage, perhaps related to its extreme habitat specialization, yet plastid genome length, structure, and substitution rate are within the variation previously reported for photosynthetic plants. These results highlight the increasingly appreciated dynamism of plastid genomes, otherwise highly conserved across a billion years of green plant evolution, in plants with highly specialized life history traits.

Comparative chloroplast genome analysis of Artemisia (Asteraceae) in East Asia: insights into evolutionary divergence and phylogenomic implications

BackgroundArtemisia in East Asia includes a number of economically important taxa that are widely used for food, medicinal, and ornamental purposes. The identification of taxa, however, has been hampered by insufficient diagnostic morphological characteristics and frequent natural hybridization. Development of novel DNA markers or barcodes with sufficient resolution to resolve taxonomic issues of Artemisia in East Asia is significant challenge.ResultsTo establish a molecular basis for taxonomic identification and comparative phylogenomic analysis of Artemisia, we newly determined 19 chloroplast genome (plastome) sequences of 18 Artemisia taxa in East Asia, de novo-assembled and annotated the plastomes of two taxa using publicly available Illumina reads, and compared them with 11 Artemisia plastomes reported previously. The plastomes of Artemisia were 150,858–151,318 base pairs (bp) in length and harbored 87 protein-coding genes, 37 transfer RNAs, and 8 ribosomal RNA genes in conserved order and orientation. Evolutionary analyses of whole plastomes and 80 non-redundant protein-coding genes revealed that the noncoding trnH-psbA spacer was highly variable in size and nucleotide sequence both between and within taxa, whereas the coding sequences of accD and ycf1 were under weak positive selection and relaxed selective constraints, respectively. Phylogenetic analysis of the whole plastomes based on maximum likelihood and Bayesian inference analyses yielded five groups of Artemisia plastomes clustered in the monophyletic subgenus Dracunculus and paraphyletic subgenus Artemisia, suggesting that the whole plastomes can be used as molecular markers to infer the chloroplast haplotypes of Artemisia taxa. Additionally, analysis of accD and ycf1 hotspots enabled the development of novel markers potentially applicable across the family Asteraceae with high discriminatory power.ConclusionsThe complete sequences of the Artemisia plastomes are sufficiently polymorphic to be used as super-barcodes for this genus. It will facilitate the development of new molecular markers and study of the phylogenomic relationships of Artemisia species in the family Asteraceae.

Genome-wide characterization of the laccase gene family in Setaria viridis reveals members potentially involved in lignification.

Five laccase genes are potentially involved in developmental lignification in the model C4 grass Setaria viridis and their different tissue specificities suggest subfunctionalization events. Plant laccases are copper-containing glycoproteins involved in monolignol oxidation and, therefore, their activity is essential for lignin polymerization. Although these enzymes belong to large multigene families with highly redundant members, not all of them are thought to be involved in lignin metabolism. Here, we report on the genome-wide characterization of the laccase gene family in the model C4 grass Setaria viridis and further identification of the members potentially involved in monolignol oxidation. A total of 52 genes encoding laccases (SvLAC1 to SvLAC52) were found in the genome of S. viridis, and phylogenetic analyses showed that these genes were heterogeneously distributed among the characteristic six subclades of the family and are under relaxed selective constraints. The observed expansion in the total number of genes in this species was mainly caused by tandem duplications within subclade V, which accounts for 68% of the whole family. Comparative phylogenetic analyses showed that the expansion of subclade V is specifically observed for the Paniceae tribe within the Panicoideae subfamily in grasses. Five SvLAC genes (SvLAC9, SvLAC13, SvLAC15, SvLAC50, and SvLAC52) fulfilled the criteria established to identify lignin-related candidates: (1) phylogenetic proximity to previously characterized lignin-related laccases from other species, (2) similar expression pattern to that observed for lignin biosynthetic genes in the S. viridis elongating internode, and (3) high expression in S. viridis tissues undergoing active lignification. In addition, in situ hybridization experiments not only confirmed that these selected SvLAC genes were expressed in lignifying cells, but also that their expression showed different tissue specificities, suggesting subfunctionalization events within the family. These five laccase genes are strong candidates to be involved in lignin polymerization in S. viridis and might be good targets for lignin bioengineering strategies.