Gene duplication as a mechanism of genomic adaptation to a changing environment

Genomic resources for energy cane breeding in the post genomics era

Sugarcane is one of the most sustainable energy crops among cultivated crops presenting the highest tonnage of cultivated plants. Its high productivity of sugar, bioethanol and bioelectricity make it a promising green alternative to petroleum. Furthermore, the myriad of products that can be derived from sugarcane biomass has been driving breeding programs towards varieties with a higher yield of fiber and a more vigorous and sustainable performance: the energy cane. Here we provide an overview of the energy cane including plant description, breeding efforts, types, and end-uses. In addition, we describe recently published genomic resources for the development of this crop, discuss current knowledge of cell wall metabolism, bioinformatic tools and databases available for the community.

Insights into the genome structure of four acetogenic bacteria with specific reference to the Wood-Ljungdahl pathway

Acetogenic bacteria are obligate anaerobes with the ability of converting carbon dioxide and other one‐carbon substrates into acetate through the Wood–Ljungdahl (WL) pathway. These substrates are becoming increasingly important feedstock in industrial microbiology. The main potential industrial application of acetogenic bacteria is the production of metabolites that constitute renewable energy sources (biofuel); such bacteria are of particular interest for this purpose thanks to their low energy requirements for large‐scale cultivation. Here, we report new genome sequences for four species, three of them are reported for the first time, namely Acetobacterium paludosum DSM 8237, Acetobacterium tundrae DSM 917, Acetobacterium bakii DSM 8239, and Alkalibaculum bacchi DSM 221123. We performed a comparative genomic analysis focused on the WL pathway's genes and their encoded proteins, using Acetobacterium woodii as a reference genome. The Average Nucleotide Identity (ANI) values ranged from 70% to 95% over an alignment length of 5.4–6.5 Mbp. The core genome consisted of 363 genes, whereas the number of unique genes in a single genome ranged from 486 in A. tundrae to 2360 in A.bacchi. No significant rearrangements were detected in the gene order for the Wood–Ljungdahl pathway however, two species showed variations in genes involved in formate metabolism: A. paludosum harbor two copies of fhs1, and A. bakii a truncated fdhF1. The analysis of protein networks highlighted the expansion of protein orthologues in A. woodii compared to A. bacchi, whereas protein networks involved in the WL pathway were more conserved. This study has increased our understanding on the evolution of the WL pathway in acetogenic bacteria.

Current and Promising Approaches to Identify Horizontal Gene Transfer Events in Metagenomes

High-throughput shotgun metagenomics sequencing has enabled the profiling of myriad natural communities. These data are commonly used to identify gene families and pathways that were potentially gained or lost in an environment and which may be involved in microbial adaptation. Despite the widespread interest in these events, there are no established best practices for identifying gene gain and loss in metagenomics data. Horizontal gene transfer (HGT) represents several mechanisms of gene gain that are especially of interest in clinical microbiology due to the rapid spread of antibiotic resistance genes in natural communities. Several additional mechanisms of gene gain and loss, including gene duplication, gene loss-of-function events, and de novo gene birth are also important to consider in the context of metagenomes but have been less studied. This review is largely focused on detecting HGT in prokaryotic metagenomes, but methods for detecting these other mechanisms are first discussed. For this article to be self-contained, we provide a general background on HGT and the different possible signatures of this process. Lastly, we discuss how improved assembly of genomes from metagenomes would be the most straight-forward approach for improving the inference of gene gain and loss events. Several recent technological advances could help improve metagenome assemblies: long-read sequencing, determining the physical proximity of contigs, optical mapping of short sequences along chromosomes, and single-cell metagenomics. The benefits and limitations of these advances are discussed and open questions in this area are highlighted.

Divergent Copies of a Cryptosporidium parvum-Specific Subtelomeric Gene

Subtype families of Cryptosporidium parvum differ in host range, with IIa and IId being found in a broad range of animals, IIc in humans, and IIo and IIp in some rodents. Previous studies indicated that the subtelomeric cgd6_5520-5510 gene in C. parvum is lost in many Cryptosporidium species, and could potentially contribute to the broad host range of the former. In this study, we identified the presence of a second copy of the gene in some C. parvum subtype families with a broad host range, and showed sequence differences among them. The sequence differences in the cgd6_5520-5510 gene were not segregated by the sequence type of the 60 kDa glycoprotein gene. Genetic recombination appeared to have played a role in generating divergent nucleotide sequences between copies and among subtype families. These data support the previous conclusion on the potential involvement of the insulinase-like protease encoded by the subtelomeric cgd6_5520-5510 gene in the broad host range of C. parvum IIa and IId subtypes.

Plant evolution and environmental adaptation unveiled by long-read whole-genome sequencing of Spirodela

Aquatic plants have to adapt to the environments distinct from where land plants grow. A critical aspect of adaptation is the dynamics of sequence repeats, not resolved in older sequencing platforms due to incomplete and fragmented genome assemblies from short reads. Therefore, we used PacBio long-read sequencing of the Spirodela polyrhiza genome, reaching a 44-fold increase of contiguity with an N50 (a median of contig lengths) of 831 kb and filling 95.4% of gaps left from the previous version. Reconstruction of repeat regions indicates that sequentially nested long terminal repeat (LTR) retrotranspositions occur early in monocot evolution, featured with both prokaryote-like gene-rich regions and eukaryotic repeat islands. Protein-coding genes are reduced to 18,708 gene models supported by 492,435 high-quality full-length PacBio complementary DNA (cDNA) sequences. Different from land plants, the primitive architecture of Spirodela's adventitious roots and lack of lateral roots and root hairs are consistent with dispensable functions of nutrient absorption. Disease-resistant genes encoding antimicrobial peptides and dirigent proteins are expanded by tandem duplications. Remarkably, disease-resistant genes are not only amplified, but also highly expressed, consistent with low levels of 24-nucleotide (nt) small interfering RNA (siRNA) that silence the immune system of land plants, thereby protecting Spirodela against a wide spectrum of pathogens and pests. The long-read sequence information not only sheds light on plant evolution and adaptation to the environment, but also facilitates applications in bioenergy and phytoremediation.

Fluorescence in situ hybridization in plants: recent developments and future applications

Fluorescence in situ hybridization (FISH) was developed more than 30 years ago and has been the most paradigm-changing technique in cytogenetic research. FISH has been used to answer questions related to structure, mutation, and evolution of not only individual chromosomes but also entire genomes. FISH has served as an important tool for chromosome identification in many plant species. This review intends to summarize and discuss key technical development and applications of FISH in plants since 2006. The most significant recent advance of FISH is the development and application of probes based on synthetic oligonucleotides (oligos). Oligos specific to a repetitive DNA sequence, to a specific chromosomal region, or to an entire chromosome can be computationally identified, synthesized in parallel, and fluorescently labeled. Oligo probes designed from conserved DNA sequences from one species can be used among genetically related species, allowing comparative cytogenetic mapping of these species. The advances with synthetic oligo probes will significantly expand the applications of FISH especially in non-model plant species. Recent achievements and future applications of FISH and oligo-FISH are discussed.

Complex DNA structures trigger copy number variation across the Plasmodium falciparum genome

Antimalarial resistance is a major obstacle in the eradication of the human malaria parasite, Plasmodium falciparum. Genome amplifications, a type of DNA copy number variation (CNV), facilitate overexpression of drug targets and contribute to parasite survival. Long monomeric A/T tracks are found at the breakpoints of many Plasmodium resistance-conferring CNVs. We hypothesize that other proximal sequence features, such as DNA hairpins, act with A/T tracks to trigger CNV formation. By adapting a sequence analysis pipeline to investigate previously reported CNVs, we identified breakpoints in 35 parasite clones with near single base-pair resolution. Using parental genome sequence, we predicted the formation of stable hairpins within close proximity to all future breakpoint locations. Especially stable hairpins were predicted to form near five shared breakpoints, establishing that the initiating event could have occurred at these sites. Further in-depth analyses defined characteristics of these 'trigger sites' across the genome and detected signatures of error-prone repair pathways at the breakpoints. We propose that these two genomic signals form the initial lesion (hairpins) and facilitate microhomology-mediated repair (A/T tracks) that lead to CNV formation across this highly repetitive genome. Targeting these repair pathways in P. falciparum may be used to block adaptation to antimalarial drugs.

Depth distributions of nitrite reductase (nirK) gene variants reveal spatial dynamics of thaumarchaeal ecotype populations in coastal Monterey Bay

Ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota are key players in nutrient cycling, yet large gaps remain in our understanding of their ecology and metabolism. Despite multiple lines of evidence pointing to a central role for copper-containing nitrite reductase (NirK) in AOA metabolism, the thaumarchaeal nirK gene is rarely studied in the environment. In this study, we examine the diversity of nirK in the marine pelagic environment, in light of previously described ecological patterns of pelagic thaumarchaeal populations. Phylogenetic analyses show that nirK better resolves diversification patterns of marine Thaumarchaeota, compared to the conventionally used marker gene amoA. Specifically, we demonstrate that the three major phylogenetic clusters of marine nirK correspond to the three 'ecotype' populations of pelagic Thaumarchaeota. In this context, we further examine the relative distributions of the three variant groups in metagenomes and metatranscriptomes representing two depth profiles in coastal Monterey Bay. Our results reveal that nirK effectively tracks the dynamics of thaumarchaeal ecotype populations, particularly finer-scale diversification patterns within major lineages. We also find evidence for multiple copies of nirK per genome in a fraction of thaumarchaeal cells in the water column, which must be taken into account when using it as a molecular marker.

Genomic dissection and transcriptional profiling of Cysteine-rich receptor-like kinases in five cereals and functional characterization of TaCRK68-A

Cysteine-rich receptor-like kinases (CRK) constitute one of the largest subfamily of receptor-like kinases, which play crucial roles in plant development and stress response. In total, 43, 37, 36, 38 and 170 CRK genes including duplicated genes were identified in the genome of Brachypodium distachyon, Hordeum vulgare, Oryza sativa, Sorghum bicolor and Triticum aestivum, respectively. These CRK proteins were tightly clustered into four phylogenetic groups and exhibited close syntenic relationship among orthologous genes. Majority of CRK proteins contain a transmembrane domain for plasma membrane localization. The organization of exon/intron, domains and motifs were variably conserved. Tissue-specific expression suggested the involvement of certain CRK genes in plant development. Modulated expression revealed their specific stress-responsive functions. Co-expression and interaction analysis indicated their role in signaling. Ks value and divergence time analysis suggested duplication of TaCRK genes before the hybridization of T. aestivum sub-genomes. Expression comparison of duplicated TaCRK genes revealed functional retention, neofunctionalization or pseudo-functionalization. Recombinant expression of a stress-responsive gene TaCRK68-A in Escherichia coli and Saccharomyces cerevisiae displayed enhanced tolerance against heat, drought, cold and salinity stresses. The study suggested vital functions of CRKs during development and stresses, and provides the basis for functional characterization of each gene in future studies.

Establishment of an induced memory response in Pseudomonas aeruginosa during infection of a eukaryotic host

In a given habitat, bacterial cells often experience recurrent exposures to the same environmental stimulus. The ability to memorize the past event and to adjust current behaviors can lead to efficient adaptation to the recurring stimulus. Here we demonstrate that the versatile bacterium Pseudomonas aeruginosa adopts a virulence phenotype after serial passage in the invertebrate model host Galleria mellonella. The virulence phenotype was not linked to the acquisition of genetic variations and was sustained for several generations, despite cultivation of the ex vivo virulence-adapted P. aeruginosa cells under rich medium conditions in vitro. Transcriptional reprogramming seemed to be induced by a host-specific food source, as reprogramming was also observed upon cultivation of P. aeruginosa in rich medium supplemented with polyunsaturated long-chain fatty acids. The establishment of induced memory responses adds a time dimension and seems to fill the gap between long-term evolutionary genotypic adaptation and short-term induced individual responses. Efforts to unravel the fundamental mechanisms that underlie the carry-over effect to induce such memory responses will continue to be of importance as hysteretic behavior can serve survival of bacterial populations in changing and challenging habitats.

Contraction of the ROS Scavenging Enzyme Glutathione S-Transferase Gene Family in Cetaceans

Cetaceans are a group of marine mammals whose ancestors were adaptated for life on land. Life in an aquatic environment poses many challenges for air-breathing mammals. Diving marine mammals have adapted to rapid reoxygenation and reactive oxygen species (ROS)-mediated reperfusion injury. Here, we considered the evolution of the glutathione transferase (GST) gene family which has important roles in the detoxification of endogenously-derived ROS and environmental pollutants. We characterized the cytosolic GST gene family in 21 mammalian species; cetaceans, sirenians, pinnipeds, and their terrestrial relatives. All seven GST classes were identified, showing that GSTs are ubiquitous in mammals. Some GST genes are the product of lineage-specific duplications and losses, in line with a birth-and-death evolutionary model. We detected sites with signatures of positive selection that possibly influence GST structure and function, suggesting that adaptive evolution of GST genes is important for defending mammals from various types of noxious environmental compounds. We also found evidence for loss of alpha and mu GST subclass genes in cetacean lineages. Notably, cetaceans have retained a homolog of at least one of the genes GSTA1, GSTA4, and GSTM1; GSTs that are present in both the cytosol and mitochondria. The observed variation in number and selection pressure on GST genes suggest that the gene family structure is dynamic within cetaceans.

Neofunctionalization of Mitochondrial Proteins and Incorporation into Signaling Networks in Plants

Because of their symbiotic origin, many mitochondrial proteins are well conserved across eukaryotic kingdoms. It is however less obvious how specific lineages have obtained novel nuclear-encoded mitochondrial proteins. Here, we report a case of mitochondrial neofunctionalization in plants. Phylogenetic analysis of genes containing the Domain of Unknown Function 295 (DUF295) revealed that the domain likely originated in Angiosperms. The C-terminal DUF295 domain is usually accompanied by an N-terminal F-box domain, involved in ubiquitin ligation via binding with ASK1/SKP1-type proteins. Due to gene duplication, the gene family has expanded rapidly, with 94 DUF295-related genes in Arabidopsis thaliana alone. Two DUF295 family subgroups have uniquely evolved and quickly expanded within Brassicaceae. One of these subgroups has completely lost the F-box, but instead obtained strongly predicted mitochondrial targeting peptides. We show that several representatives of this DUF295 Organellar group are effectively targeted to plant mitochondria and chloroplasts. Furthermore, many DUF295 Organellar genes are induced by mitochondrial dysfunction, whereas F-Box DUF295 genes are not. In agreement, several Brassicaceae-specific DUF295 Organellar genes were incorporated in the evolutionary much older ANAC017-dependent mitochondrial retrograde signaling pathway. Finally, a representative set of DUF295 T-DNA insertion mutants was created. No obvious aberrant phenotypes during normal growth and mitochondrial dysfunction were observed, most likely due to the large extent of gene duplication and redundancy. Overall, this study provides insight into how novel mitochondrial proteins can be created via “intercompartmental” gene duplication events. Moreover, our analysis shows that these newly evolved genes can then be specifically integrated into relevant, pre-existing coexpression networks.

Adaptation of S. cerevisiae to Fermented Food Environments Reveals Remarkable Genome Plasticity and the Footprints of Domestication

The budding yeast Saccharomyces cerevisiae can be found in the wild and is also frequently associated with human activities. Despite recent insights into the phylogeny of this species, much is still unknown about how evolutionary processes related to anthropogenic niches have shaped the genomes and phenotypes of S. cerevisiae. To address this question, we performed population-level sequencing of 82 S. cerevisiae strains from wine, flor, rum, dairy products, bakeries, and the natural environment (oak trees). These genomic data enabled us to delineate specific genetic groups corresponding to the different ecological niches and revealed high genome content variation across the groups. Most of these strains, compared with the reference genome, possessed additional genetic elements acquired by introgression or horizontal transfer, several of which were population-specific. In addition, several genomic regions in each population showed evidence of nonneutral evolution, as shown by high differentiation, or of selective sweeps including genes with key functions in these environments (e.g., amino acid transport for wine yeast). Linking genetics to lifestyle differences and metabolite traits has enabled us to elucidate the genetic basis of several niche-specific population traits, such as growth on galactose for cheese strains. These data indicate that yeast has been subjected to various divergent selective pressures depending on its niche, requiring the development of customized genomes for better survival in these environments. These striking genome dynamics associated with local adaptation and domestication reveal the remarkable plasticity of the S. cerevisiae genome, revealing this species to be an amazing complex of specialized populations.

A major role for noncoding regulatory mutations in the evolution of enzyme activity

This study investigates how evolutionary changes in enzyme activity occur. Multiple species of Drosophila flies have adapted to food with different levels of alcohol. This study uncovers genetic changes responsible for these repeated adaptive events, focusing on the main enzyme responsible for alcohol metabolism, Alcohol dehydrogenase. Better alcohol metabolism could be achieved either through changes to the enzyme itself or through changes in DNA regulatory sequences that affect how many enzyme molecules are produced. In four different cases, it was found that regulatory changes were the most frequent contributors to evolution. These findings have important implications because most studies of enzyme evolution focus exclusively on changes to protein sequence, and thus a significant source of adaptive changes may be overlooked.

The quantitative evolution of protein activity is a common phenomenon, yet we know little about any general mechanistic tendencies that underlie it. For example, an increase (or decrease) in enzyme activity may evolve from changes in protein sequence that alter specific activity, or from changes in gene expression that alter the amount of protein produced. The latter in turn could arise via mutations that affect gene transcription, posttranscriptional processes, or copy number. Here, to determine the types of genetic changes underlying the quantitative evolution of protein activity, we dissected the basis of ecologically relevant differences in Alcohol dehydrogenase (Adh) enzyme activity between and within several Drosophila species. By using recombinant Adh transgenes to map the functional divergence of ADH enzyme activity in vivo, we find that amino acid substitutions explain only a minority (0 to 25%) of between- and within-species differences in enzyme activity. Instead, noncoding substitutions that occur across many parts of the gene (enhancer, promoter, and 5′ and 3′ untranslated regions) account for the majority of activity differences. Surprisingly, one substitution in a transcriptional Initiator element has occurred in parallel in two species, indicating that core promoters can be an important natural source of the tuning of gene activity. Furthermore, we show that both regulatory and coding substitutions contribute to fitness (resistance to ethanol toxicity). Although qualitative changes in protein specificity necessarily derive from coding mutations, these results suggest that regulatory mutations may be the primary source of quantitative changes in protein activity, a possibility overlooked in most analyses of protein evolution.

Genomic signatures accompanying the dietary shift to phytophagy in polyphagan beetles

BACKGROUND

The diversity and evolutionary success of beetles (Coleoptera) are proposed to be related to the diversity of plants on which they feed. Indeed, the largest beetle suborder, Polyphaga, mostly includes plant eaters among its approximately 315,000 species. In particular, plants defend themselves with a diversity of specialized toxic chemicals. These may impose selective pressures that drive genomic diversification and speciation in phytophagous beetles. However, evidence of changes in beetle gene repertoires driven by such interactions remains largely anecdotal and without explicit hypothesis testing.

RESULTS

We explore the genomic consequences of beetle-plant trophic interactions by performing comparative gene family analyses across 18 species representative of the two most species-rich beetle suborders. We contrast the gene contents of species from the mostly plant-eating suborder Polyphaga with those of the mainly predatory Adephaga. We find gene repertoire evolution to be more dynamic, with significantly more adaptive lineage-specific expansions, in the more speciose Polyphaga. Testing the specific hypothesis of adaptation to plant feeding, we identify families of enzymes putatively involved in beetle-plant interactions that underwent adaptive expansions in Polyphaga. There is notable support for the selection hypothesis on large gene families for glutathione S-transferase and carboxylesterase detoxification enzymes.

CONCLUSIONS

Our explicit modeling of the evolution of gene repertoires across 18 species identifies putative adaptive lineage-specific gene family expansions that accompany the dietary shift towards plants in beetles. These genomic signatures support the popular hypothesis of a key role for interactions with plant chemical defenses, and for plant feeding in general, in driving beetle diversification.

Fine-Tuning the Expression of Duplicate Genes by Translational Regulation in Arabidopsis and Maize

Plant genomes are extensively shaped by various types of gene duplication. However, in this active area of investigation, the vast majority of studies focus on the sequence and transcription of duplicate genes, leaving open the question of how translational regulation impacts the expression and evolution of duplicate genes. We explored this issue by analyzing the ribo- and mRNA-seq data sets across six tissue types and stress conditions in Arabidopsis thaliana and maize (Zea mays). We dissected the relative contributions of transcriptional and translational regulation to the divergence in the abundance of ribosome footprint (RF) for different types of duplicate genes. We found that the divergence in RF abundance was largely programmed at the transcription level and that translational regulation plays more of a modulatory role. Intriguingly, translational regulation is characterized by its strong directionality, with the divergence in translational efficiency (TE) globally counteracting the divergence in mRNA abundance, indicating partial buffering of the transcriptional divergence between paralogs by translational regulation. Divergence in TE was associated with several sequence features. The faster-evolving copy in a duplicate pair was more likely to show lower RF abundance, which possibly results from relaxed purifying selection compared with its paralog. A considerable proportion of duplicates displayed differential TE across tissue types and stress conditions, most of which were enriched in photosynthesis, energy production, and translation-related processes. Additionally, we constructed a database TDPDG-DB (http://www.plantdupribo.tk), providing an online platform for data exploration. Overall, our study illustrates the roles of translational regulation in fine-tuning duplicate gene expression in plants.

Variation in transcriptional responses to copper exposure across Daphnia pulex lineages

Copper pollution is pervasive in aquatic habitats and is particularly harmful to invertebrates sensitive to environmental changes such as Daphnia pulex. Mechanisms of toxicity and tolerance to copper are not well understood. We used RNA-sequencing to investigate these mechanisms in three genetically distinct D. pulex clonal lineages with different histories of copper exposure. Upregulated genes after copper exposure were enriched with Gene Ontology (GO) categories involved in digestion, molting and growth, whereas downregulated genes after copper exposure were enriched in the metal-regulatory system, immune response and epigenetic modifications. The three D. pulex clones in our study show largely similar transcriptional patterns in response to copper, with only a total of twenty genes differentially expressed in a single clonal lineages. We also detected lower relative expression of some genes known to be important for copper tolerance, metallothionein and glutathione-S-transferase, in a sensitive lineage sampled from an uncontaminated habitat. Daphnia-specific genes (without orthologs outside the genus) and Daphnia-specific duplications (genes duplicated in the Daphnia lineage) were overrepresented in differentially expressed genes, highlighting an important role for newly emerged genes in tolerating environmental stressors. The results indicate that the D. pulex lineages tested in this study generally respond to copper stress using the same major pathways, but that the more resistant clone with previous copper exposure might be better able to regulate key genes. This finding highlights the important nuances in gene expression among clones, shaped by historical exposure and influencing copper tolerance.

Gene duplications in the E. coli genome: common themes among pathotypes

Background

Gene duplication underlies a significant proportion of gene functional diversity and genome complexity in both eukaryotes and prokaryotes. Although several reports in the literature described the duplication of specific genes in E. coli, a detailed analysis of the extent of gene duplications in this microorganism is needed.

Results

The genomes of the E. coli enteroaggregative strain 042 and other pathogenic strains contain duplications of the gene that codes for the global regulator Hha. To determine whether the presence of additional copies of the hha gene correlates with the presence of other genes, we performed a comparative genomic analysis between E. coli strains with and without hha duplications. The results showed that strains harboring additional copies of the hha gene also encode the yeeR irmA (aec69) gene cluster, which, in turn, is also duplicated in strain 042 and several other strains. The identification of these duplications prompted us to obtain a global map of gene duplications, first in strain 042 and later in other E. coli genomes.

Duplications in the genomes of the enteroaggregative strain 042, the uropathogenic strain CFT073 and the enterohemorrhagic strain O145:H28 have been identified by a BLASTp protein similarity search. This algorithm was also used to evaluate the distribution of the identified duplicates among the genomes of a set of 28 representative E. coli strains. Despite the high genomic diversity of E. coli strains, we identified several duplicates in the genomes of almost all studied pathogenic strains. Most duplicated genes have no known function. Transcriptomic analysis also showed that most of these duplications are regulated by the H-NS/Hha proteins.

Conclusions

Several duplicated genes are widely distributed among pathogenic E. coli strains. In addition, some duplicated genes are present only in specific pathotypes, and others are strain specific. This gene duplication analysis shows novel relationships between E. coli pathotypes and suggests that newly identified genes that are duplicated in a high percentage of pathogenic E. coli isolates may play a role in virulence. Our study also shows a relationship between the duplication of genes encoding regulators and genes encoding their targets.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5683-4) contains supplementary material, which is available to authorized users.

Genome-wide identification and expression profiling of trihelix gene family under abiotic stresses in wheat

BACKGROUND

The trihelix gene family is a plant-specific transcription factor family that plays important roles in plant growth, development, and responses to abiotic stresses. However, to date, no systemic characterization of the trihelix genes has yet been conducted in wheat and its close relatives.

RESULTS

We identified a total of 94 trihelix genes in wheat, as well as 22 trihelix genes in Triticum urartu, 29 in Aegilops tauschii, and 31 in Brachypodium distachyon. We analyzed the chromosomal locations and orthology relations of the identified trihelix genes, and no trihelix gene was found to be located on chromosome 7A, 7B, or 7D of wheat, thereby reflecting the uneven distributions of wheat trihelix genes. Phylogenetic analysis indicated that the 186 identified trihelix proteins in wheat, rice, B. distachyon, and Arabidopsis were clustered into five major clades. The trihelix genes belonging to the same clades usually shared similar motif compositions and exon/intron structural patterns. Five pairs of tandem duplication genes and three pairs of segmental duplication genes were identified in the wheat trihelix gene family, thereby validating the supposition that more intrachromosomal gene duplication events occur in the genome of wheat than in that of other grass species. The tissue-specific expression and differential expression profiling of the identified genes under cold and drought stresses were analyzed by using RNA-seq data. qRT-PCR was also used to confirm the expression profiles of ten selected wheat trihelix genes under multiple abiotic stresses, and we found that these genes mainly responded to salt and cold stresses.

CONCLUSIONS

In this study, we identified trihelix genes in wheat and its close relatives and found that gene duplication events are the main driving force for trihelix gene evolution in wheat. Our expression profiling analysis demonstrated that wheat trihelix genes responded to multiple abiotic stresses, especially salt and cold stresses. The results of our study built a basis for further investigation of the functions of wheat trihelix genes and provided candidate genes for stress-resistant wheat breeding programs.

Landscape of gene transposition-duplication within the Brassicaceae family

We developed the CLfinder-OrthNet pipeline that detects co-linearity among multiple closely related genomes, finds orthologous gene groups, and encodes the evolutionary history of each orthologue group into a representative network (OrthNet). Using a search based on network topology, we identified 1,394 OrthNets that included gene transposition-duplication (tr-d) events, out of 17,432 identified in six Brassicaceae genomes. Occurrences of tr-d shared by subsets of Brassicaceae genomes mirrored the divergence times between the genomes and their repeat contents. The majority of tr-d events resulted in truncated open reading frames (ORFs) in the duplicated loci. However, the duplicates with complete ORFs were significantly more frequent than expected from random events. These were derived from older tr-d events and had a higher chance of being expressed. We also found an enrichment of tr-d events with complete loss of intergenic sequence conservation between the original and duplicated loci. Finally, we identified tr-d events uniquely found in two extremophytes among the six Brassicaceae genomes, including tr-d of SALT TOLERANCE 32 and ZINC TRANSPORTER 3 that relate to their adaptive evolution. CLfinder-OrthNet provides a flexible toolkit to compare gene order, visualize evolutionary paths among orthologues as networks, and identify gene loci that share an evolutionary history.

Copy Number Variation in Domestication

Domesticated plants have long served as excellent models for studying evolution. Many genes and mutations underlying important domestication traits have been identified, and most causal mutations appear to be SNPs. Copy number variation (CNV) is an important source of genetic variation that has been largely neglected in studies of domestication. Ongoing work demonstrates the importance of CNVs as a source of genetic variation during domestication, and during the diversification of domesticated taxa. Here, we review how CNVs contribute to evolutionary processes underlying domestication, and review examples of domestication traits caused by CNVs. We draw from examples in plant species, but also highlight cases in animal systems that could illuminate the roles of CNVs in the domestication process.

Fluorescence in situ hybridization in plants: recent developments and future applications

Fluorescence in situ hybridization (FISH) was developed more than 30 years ago and has been the most paradigm-changing technique in cytogenetic research. FISH has been used to answer questions related to structure, mutation, and evolution of not only individual chromosomes but also entire genomes. FISH has served as an important tool for chromosome identification in many plant species. This review intends to summarize and discuss key technical development and applications of FISH in plants since 2006. The most significant recent advance of FISH is the development and application of probes based on synthetic oligonucleotides (oligos). Oligos specific to a repetitive DNA sequence, to a specific chromosomal region, or to an entire chromosome can be computationally identified, synthesized in parallel, and fluorescently labeled. Oligo probes designed from conserved DNA sequences from one species can be used among genetically related species, allowing comparative cytogenetic mapping of these species. The advances with synthetic oligo probes will significantly expand the applications of FISH especially in non-model plant species. Recent achievements and future applications of FISH and oligo-FISH are discussed.

Linking Duplication of a Calcium Sensor to Salt Tolerance in Eutrema salsugineum

The SALT-OVERLY-SENSITIVE (SOS) pathway in Arabidopsis (Arabidopsis thaliana) functions to prevent the toxic accumulation of sodium in the cytosol when plants are grown in salt-affected soils. In this pathway, the CALCINEURIN B-LIKE10 (AtCBL10) calcium sensor interacts with the AtSOS2 kinase to activate the AtSOS1 plasma membrane sodium/proton exchanger. CBL10 has been duplicated in Eutrema (Eutrema salsugineum), a salt-tolerant relative of Arabidopsis. Because Eutrema maintains growth in salt-affected soils that kill most crop plants, the duplication of CBL10 provides a unique opportunity to functionally test the outcome of gene duplication and its link to plant salt tolerance. In Eutrema, individual down-regulation of the duplicated CBL10 genes (EsCBL10a and EsCBL10b) decreased growth in the presence of salt and, in combination, led to an even greater decrease, suggesting that both genes function in response to salt and have distinct functions. Cross-species complementation assays demonstrated that EsCBL10b has an enhanced ability to activate the SOS pathway while EsCBL10a has a function not performed by AtCBL10 or EsCBL10b Chimeric EsCBL10a/EsCBL10b proteins revealed that the specific functions of the EsCBL10 proteins resulted from changes in the amino terminus. The duplication of CBL10 increased calcium-mediated signaling capacity in Eutrema and conferred increased salt tolerance to salt-sensitive Arabidopsis.

The molecular basis for the neofunctionalization of the juvenile hormone esterase duplication in Drosophila

The Drosophila melanogaster enzymes juvenile hormone esterase (DmJHE) and its duplicate, DmJHEdup, present ideal examples for studying the structural changes involved in the neofunctionalization of enzyme duplicates. DmJHE is a hormone esterase with precise regulation and highly specific activity for its substrate, juvenile hormone. DmJHEdup is an odorant degrading esterase (ODE) responsible for processing various kairomones in antennae. Our phylogenetic analysis shows that the JHE lineage predates the hemi/holometabolan split and that several duplications of JHEs have been templates for the evolution of secreted β-esterases such as ODEs through the course of insect evolution. Our biochemical comparisons further show that DmJHE has sufficient substrate promiscuity and activity against odorant esters for a duplicate to evolve a general ODE function against a range of mid-long chain food esters, as is shown in DmJHEdup. This substrate range complements that of the only other general ODE known in this species, Esterase 6. Homology models of DmJHE and DmJHEdup enabled comparisons between each enzyme and the known structures of a lepidopteran JHE and Esterase 6. Both JHEs showed very similar active sites despite low sequence identity (30%). Both ODEs differed drastically from the JHEs and each other, explaining their complementary substrate ranges. A small number of amino acid changes are identified that may have been involved in the early stages of the neofunctionalization of DmJHEdup. Our results provide key insights into the process of neofunctionalization and the structural changes that can be involved.

Paralogues From the Expanded Tlr11 Gene Family in Mudskipper (Boleophthalmus pectinirostris) Are Under Positive Selection and Respond Differently to LPS/Poly(I:C) Challenge

Toll-like receptors (TLRs) are major molecular pattern recognition receptors, which are essential for triggering a series of innate immune responses against invading pathogens by recognizing their evolutionary conserved molecular patterns. The mudskipper, Boleophthalmus pectinirostris is exceptional among fishes due to its amphibious lifestyle and adaptation to living on mudflats. The whole-genome sequencing of B. pectinirostris has revealed that this species possesses an expansion of Tlr11 family [12 Tlr11 family genes (one tlr21, 4 tlr22, and 7 tlr23)] that we focused on in the present study. The full-length cDNA sequences of the 12 tlrs in B. pectinirostris were cloned and their deduced amino acid sequences possessed a typical TLR domain arrangement. Likelihood tests of selection revealed that these 12 Tlr11 family genes are under diversifying selection. A total of 13 sites were found to be positively selected by more than one evolution model, of which 11 were located in the ligand-binding ectodomain. The observed non-synonymous substitutions may have functional implications in antigen and pathogen recognition specificity. These 12 tlrs were highly expressed in immune-related tissues, i.e. spleen and kidney. Tlr21 and tlr22b transcripts were significantly up-regulated by LPS, whereas tlr22a, tlr22d, tlr23b, tlr23e, tlr23g were significantly up-regulated by poly(I:C) in the spleen or/and kidney, which implies that the expanded Tlr11 family genes may play roles in protecting the fish from the invasion of gram-negative bacteria and double-stranded RNA viruses. The results from the present study suggested that the expansion of Tlr11 family genes in B. pectinirostris may recognize ligands from various pathogens found in the intertidal zone.

Invasive invertebrates associated with highly duplicated gene content

Invasion of alien species has led to serious problems, including the destruction of native ecosystems. In general, invasive species adapt to new environments rapidly, suggesting that they have high genetic diversity that can directly influence environmental adaptability. However, it is not known how genomic architecture causes genetic diversity that leads to invasiveness. Recent studies have showed that the proportion of duplicated genes (P_D ) in whole animal genomes correlate with environmental variability within a habitat. Here, we show that P_D and propagule size significantly explain the differences in species categories (invasive species, noninvasive species, and parasites). P_D correlated negatively with the propagule size. The residual values of regression of P_D on propagule size revealed that the invasive species had higher P_D values and larger propagule size than those of the noninvasive species, whereas the parasites had lower P_D values and smaller propagule size than those of others. There were no correlations between the invasive species and other genomic factors including the genome size, number of genes, and certain gene families. Our results suggest that the P_D values of a genome might be a potential genomic source causing genetic variations for adaptation to diverse environments. The results also showed that the invasiveness status of a species would be predicted by the residual values of regression of P_D on propagule size. Our innovative approach provides a measure to estimate the environmental adaptability of organisms based on genomic data.

Characterization and evolutionary dynamics of complex regions in eukaryotic genomes

Complex regions in eukaryotic genomes are typically characterized by duplications of chromosomal stretches that often include one or more genes repeated in a tandem array or in relatively close proximity. Nevertheless, the repetitive nature of these regions, together with the often high sequence identity among repeats, have made complex regions particularly recalcitrant to proper molecular characterization, often being misassembled or completely absent in genome assemblies. This limitation has prevented accurate functional and evolutionary analyses of these regions. This is becoming increasingly relevant as evidence continues to support a central role for complex genomic regions in explaining human disease, developmental innovations, and ecological adaptations across phyla. With the advent of long-read sequencing technologies and suitable assemblers, the development of algorithms that can accommodate sample heterozygosity, and the adoption of a pangenomic-like view of these regions, accurate reconstructions of complex regions are now within reach. These reconstructions will finally allow for accurate functional and evolutionary studies of complex genomic regions, underlying the generation of genotype-phenotype maps of unprecedented resolution.

The genome of the soybean cyst nematode (Heterodera glycines) reveals complex patterns of duplications involved in the evolution of parasitism genes

BACKGROUND

Heterodera glycines, commonly referred to as the soybean cyst nematode (SCN), is an obligatory and sedentary plant parasite that causes over a billion-dollar yield loss to soybean production annually. Although there are genetic determinants that render soybean plants resistant to certain nematode genotypes, resistant soybean cultivars are increasingly ineffective because their multi-year usage has selected for virulent H. glycines populations. The parasitic success of H. glycines relies on the comprehensive re-engineering of an infection site into a syncytium, as well as the long-term suppression of host defense to ensure syncytial viability. At the forefront of these complex molecular interactions are effectors, the proteins secreted by H. glycines into host root tissues. The mechanisms of effector acquisition, diversification, and selection need to be understood before effective control strategies can be developed, but the lack of an annotated genome has been a major roadblock.

RESULTS

Here, we use PacBio long-read technology to assemble a H. glycines genome of 738 contigs into 123 Mb with annotations for 29,769 genes. The genome contains significant numbers of repeats (34%), tandem duplicates (18.7 Mb), and horizontal gene transfer events (151 genes). A large number of putative effectors (431 genes) were identified in the genome, many of which were found in transposons.

CONCLUSIONS

This advance provides a glimpse into the host and parasite interplay by revealing a diversity of mechanisms that give rise to virulence genes in the soybean cyst nematode, including: tandem duplications containing over a fifth of the total gene count, virulence genes hitchhiking in transposons, and 107 horizontal gene transfers not reported in other plant parasitic nematodes thus far. Through extensive characterization of the H. glycines genome, we provide new insights into H. glycines biology and shed light onto the mystery underlying complex host-parasite interactions. This genome sequence is an important prerequisite to enable work towards generating new resistance or control measures against H. glycines.

A Likely Ancient Genome Duplication in the Speciose Reef-Building Coral Genus, Acropora

Whole-genome duplication (WGD) has been recognized as a significant evolutionary force in the origin and diversification of multiple organisms. Acropora, a speciose reef-building coral genus, is suspected to have originated by polyploidy. Yet, there is no genetic evidence to support this hypothesis. Using comprehensive phylogenomic and comparative genomic approaches, we analyzed six Acroporid genomes and found that a WGD event likely occurred ∼31 million years ago in the most recent common ancestor of Acropora, concurrent with a worldwide coral extinction. We found that duplicated genes were highly enriched in gene regulation functions, including those of stress responses. The functional clusters of duplicated genes are related to the divergence of gene expression patterns during development. Some proteinaceous toxins were generated by WGD in Acropora compared with other cnidarian species. Collectively, this study provides evidence for an ancient WGD event in corals, which helps explain the origin and diversification of Acropora.

•

An ancient genome duplication occurred in the most recent common ancestor of Acropora

•

This WGD event likely occurred between 28 and 36 mya in Acropora

•

The WGD event potentially contributes to the origin and diversification of Acropora

•

Duplications of toxic proteins were found in Acropora following the WGD

Dissecting the Genomic Diversification of Late Embryogenesis Abundant (LEA) Protein Gene Families in Plants

Late embryogenesis abundant (LEA) proteins include eight multigene families that are expressed in response to water loss during seed maturation and in vegetative tissues of desiccation tolerant species. To elucidate LEA proteins evolution and diversification, we performed a comprehensive synteny and phylogenetic analyses of the eight gene families across 60 complete plant genomes. Our integrated comparative genomic approach revealed that synteny conservation and diversification contributed to LEA family expansion and functional diversification in plants. We provide examples that: 1) the genomic diversification of the Dehydrin family contributed to differential evolution of amino acid sequences, protein biochemical properties, and gene expression patterns, and led to the appearance of a novel functional motif in angiosperms; 2) ancient genomic diversification contributed to the evolution of distinct intrinsically disordered regions of LEA_1 proteins; 3) recurrent tandem-duplications contributed to the large expansion of LEA_2; and 4) dynamic synteny diversification played a role on the evolution of LEA_4 and its function on plant desiccation tolerance. Taken together, these results show that multiple evolutionary mechanisms have not only led to genomic diversification but also to structural and functional plasticity among LEA proteins which have jointly contributed to the adaptation of plants to water-limiting environments.

Expansion and Diversification of MFS Transporters in Kluyveromyces marxianus

In yeasts, proteins of the Major Superfamily Transporter selectively bind and allow the uptake of sugars to permit growth on varied substrates. The genome of brewer’s yeast, Saccharomyces cerevisiae, encodes multiple hexose transporters (Hxt) to transport glucose and other MFS proteins for maltose, galactose, and other monomers. For sugar uptake, the dairy yeast, Kluyveromyces lactis, uses Rag1p for glucose, Hgt1 for glucose and galactose, and Lac12 for lactose. In the related industrial species Kluyveromyces marxianus, there are four genes encoding Lac12-like proteins but only one of them, Lac12, can transport lactose. In this study, which initiated with efforts to investigate possible functions encoded by the additional LAC12 genes in K. marxianus, a genome-wide survey of putative MFS sugar transporters was performed. Unexpectedly, it was found that the KHT and the HGT genes are present as tandem arrays of five to six copies, with the precise number varying between isolates. Heterologous expression of individual genes in S. cerevisiae and mutagenesis of single and multiple genes in K. marxianus was performed to establish possible substrates for these transporters. The focus was on the sugar galactose since it was already reported in K. lactis that this hexose was a substrate for both Lac12 and Hgt1. It emerged that three of the four copies of Lac12, four Hgt-like proteins and one Kht-like protein have some capacity to transport galactose when expressed in S. cerevisiae and inactivation of all eight genes was required to completely abolish galactose uptake in K. marxianus. Analysis of the amino acid sequence of all known yeast galactose transporters failed to identify common residues that explain the selectivity for galactose. Instead, the capacity to transport galactose has arisen three different times in K. marxianus via polymorphisms in proteins that are probably ancestral glucose transporters. Although, this is analogous to S. cerevisiae, in which Gal2 is related to glucose transporters, there are not conserved amino acid changes, either with Gal2, or among the K. marxianus galactose transporters. The data highlight how gene duplication and functional diversification has provided K. marxianus with versatile capacity to utilise sugars for growth.

The genomic basis of adaptation to calcareous and siliceous soils in Arabidopsis lyrata

Edaphic conditions are important determinants of plant fitness. While much has been learnt in recent years about plant adaptation to heavy metal contaminated soils, the genomic basis underlying adaptation to calcareous and siliceous substrates remains largely unknown. We performed a reciprocal germination experiment and whole-genome resequencing in natural calcareous and siliceous populations of diploid Arabidopsis lyrata to test for edaphic adaptation and detect signatures of selection at loci associated with soil-mediated divergence. In parallel, genome scans on respective diploid ecotypes from the Arabidopsis arenosa species complex were undertaken, to search for shared patterns of adaptive genetic divergence. Soil ecotypes of A. lyrata display significant genotype-by-treatment responses for seed germination. Sequence (SNPs) and copy-number variants (CNVs) point towards loci involved in ion transport as the main targets of adaptive genetic divergence. Two genes exhibiting high differentiation among soil types in A. lyrata further share trans-specific single nucleotide polymorphisms with A. arenosa. This work applies experimental and genomic approaches to study edaphic adaptation in A. lyrata and suggests that physiological response to elemental toxicity and deficiency underlies the evolution of calcareous and siliceous ecotypes. The discovery of shared adaptive variation between sister species indicates that ancient polymorphisms contribute to adaptive evolution.

Evolution and Expression Characteristics of Receptor-Like Cytoplasmic Protein Kinases in Maize, Rice and Arabidopsis

Receptor-like cytoplasmic protein kinases (RLCKs) are involved in various activities in plant growth and development. We have totally identified 162, 160, and 402 RLCK genes in maize, rice, and Arabidopsis genomes, respectively. Phylogenetic analyses divided 724 RLCK genes into 15 subfamilies and similar structural patterns of kinase activity sites and functional sites were observed within the subfamilies. Furthermore, the structural patterns of intron/exon in the same subfamilies were similar, implicating their close evolutionary relationship. Chromosome distribution indicated that segmental duplication of RLCK genes might be a major mechanism contributing to the expansion of the RLCK superfamilies in maize, rice, and Arabidopsis, respectively. The analysis of the synteny relationship and gene structure indicated that the evolution of most RLCKs in maize were prior to rice and Arabidopsis. Most of the ratio of Ka/Ks is inferior to one, suggesting that RLCK genes have experienced the negative selection in maize, rice and Arabidopsis. Duplication time revealed that the maize was the earliest emergence among these three species. The expression profiles showed that there are some specifically expressed RLCK genes in maize root, leaf, ear, and tassel. These specific expression genes may participate in the developmental regulation of these maize tissues. Our results will be useful in providing new insights into evolution of RLCKs and revealing the regulatory network of maize, rice, and Arabidopsis development.

Comparative analysis of selection mode reveals different evolutionary rate and expression pattern in Arachis duranensis and Arachis ipaënsis duplicated genes

Our results reveal that K_s is a determining factor affecting selective pressure and different evolution and expression patterns are detected between PSGs and NSGs in wild Arachis duplicates. Selective pressure, including purifying (negative) and positive selection, can be detected in organisms. However, studies on comparative evolutionary rates, gene expression patterns and gene features between negatively selected genes (NSGs) and positively selected genes (PSGs) are lagging in paralogs of plants. Arachis duranensis and Arachis ipaënsis are ancestors of the cultivated peanut, an important oil and protein crop. Here, we carried out a series of systematic analyses, comparing NSG and PSG in paralogs, using genome sequences and transcriptome datasets in A. duranensis and A. ipaënsis. We found that synonymous substitution rate (K_s) is a determining factor affecting selective pressure in A. duranensis and A. ipaënsis duplicated genes. Lower expression level, lower gene expression breadth, higher codon bias and shorter polypeptide length were found in PSGs and not in NSGs. The correlation analyses showed that gene expression breadth was positively correlated with polypeptide length and GC content at the first codon site (GC1) in PSGs and NSGs, respectively. There was a negative correlation between expression level and polypeptide length in PSGs. In NSGs, the K_s was positively correlated with expression level, gene expression breadth, GC1, and GC content at the third codon site (GC3), but selective pressure was negatively correlated with expression level, gene expression breadth, polypeptide length, GC1, and GC3 content. The function of most duplicated gene pairs was divergent under drought and nematode stress. Taken together, our results show that different evolution and expression patterns occur between PSGs and NSGs in paralogs of two wild Arachis species.

Comparative analysis of the Pocillopora damicornis genome highlights role of immune system in coral evolution

Comparative analysis of the expanding genomic resources for scleractinian corals may provide insights into the evolution of these organisms, with implications for their continued persistence under global climate change. Here, we sequenced and annotated the genome of Pocillopora damicornis, one of the most abundant and widespread corals in the world. We compared this genome, based on protein-coding gene orthology, with other publicly available coral genomes (Cnidaria, Anthozoa, Scleractinia), as well as genomes from other anthozoan groups (Actiniaria, Corallimorpharia), and two basal metazoan outgroup phlya (Porifera, Ctenophora). We found that 46.6% of P. damicornis genes had orthologs in all other scleractinians, defining a coral ‘core’ genome enriched in basic housekeeping functions. Of these core genes, 3.7% were unique to scleractinians and were enriched in immune functionality, suggesting an important role of the immune system in coral evolution. Genes occurring only in P. damicornis were enriched in cellular signaling and stress response pathways, and we found similar immune-related gene family expansions in each coral species, indicating that immune system diversification may be a prominent feature of scleractinian coral evolution at multiple taxonomic levels. Diversification of the immune gene repertoire may underlie scleractinian adaptations to symbiosis, pathogen interactions, and environmental stress.

Expression of Novel Gene Content Drives Adaptation to Low Iron in the Cyanobacterium Acaryochloris

Variation in genome content is a potent mechanism of microbial adaptation. The genomes of members of the cyanobacterial genus Acaryochloris vary greatly in gene content as a consequence of the idiosyncratic retention of both recent gene duplicates and plasmid-encoded genes acquired by horizontal transfer. For example, the genome of Acaryochloris strain MBIC11017, which was isolated from an iron-limited environment, is enriched in duplicated and novel genes involved in iron assimilation. Here, we took an integrative approach to characterize the adaptation of Acaryochloris MBIC11017 to low environmental iron availability and the relative contributions of the expression of duplicated versus novel genes. We observed that Acaryochloris MBIC11017 grew faster and to a higher yield in the presence of nanomolar concentrations of iron than did a closely related strain. These differences were associated with both a higher rate of iron assimilation and a greater abundance of iron assimilation transcripts. However, recently duplicated genes contributed little to increased transcript dosage; rather, the maintenance of these duplicates in the MBIC11017 genome is likely due to the sharing of ancestral dosage by expression reduction. Instead, novel, horizontally transferred genes are responsible for the differences in transcript abundance. The study provides insights on the mechanisms of adaptive genome evolution and gene expression in Acaryochloris.

Fast Evolution and Lineage-Specific Gene Family Expansions of Aphid Salivary Effectors Driven by Interactions with Host-Plants

Effector proteins play crucial roles in plant-parasite interactions by suppressing plant defenses and hijacking plant physiological responses to facilitate parasite invasion and propagation. Although effector proteins have been characterized in many microbial plant pathogens, their nature and role in adaptation to host plants are largely unknown in insect herbivores. Aphids rely on salivary effector proteins injected into the host plants to promote phloem sap uptake. Therefore, gaining insight into the repertoire and evolution of aphid effectors is key to unveiling the mechanisms responsible for aphid virulence and host plant specialization. With this aim in mind, we assembled catalogues of putative effectors in the legume specialist aphid, Acyrthosiphon pisum, using transcriptomics and proteomics approaches. We identified 3,603 candidate effector genes predicted to be expressed in A. pisum salivary glands (SGs), and 740 of which displayed up-regulated expression in SGs in comparison to the alimentary tract. A search for orthologs in 17 arthropod genomes revealed that SG-up-regulated effector candidates of A. pisum are enriched in aphid-specific genes and tend to evolve faster compared with the whole gene set. We also found that a large fraction of proteins detected in the A. pisum saliva belonged to three gene families, of which certain members show evidence consistent with positive selection. Overall, this comprehensive analysis suggests that the large repertoire of effector candidates in A. pisum constitutes a source of novelties promoting plant adaptation to legumes.

Transcriptome and miRNAs analyses enhance our understanding of the evolutionary advantages of polyploidy

Polyploid organisms have more than two sets of chromosomes, including autopolyploid via intraspecific genome doubling, and allopolyploid via merging genomes of distinct species by hybridization. Polyploid organisms are widespread in plants, indicating that polyploidy has some evolutionary advantages over its diploid ancestor. Actually, polyploidy is always tightly associated with hybrid vigor and adaptation to adverse environmental conditions. However, why polyploidy can develop such advantages is poorly known. MicroRNAs (miRNAs) are endogenous ∼21 nt small RNAs which can play important regulatory roles in animals and plants by targeting mRNAs for cleavage or translational repression. MicroRNAs are essential for cell development, differentiation, signal transduction, and show an adaptive response to biotic and abiotic stresses. Environmental stresses cause plants to over- or under-express certain miRNAs or synthesize new miRNAs to cope with stress. We have here reviewed our current knowledge on the molecular mechanisms, which can account for the evolutionary advantages of polyploidy over its diploid ancestor from genome-wide gene expression and microRNAs expression perspectives.

Comparative genomics and transcriptomics of Chrysolophus provide insights into the evolution of complex plumage coloration

Background

As one of the most recognizable characteristics in birds, plumage color has a high impact on understanding the evolution and mechanisms of coloration. Feather and skin are ideal tissues to explore the genomics and complexity of color patterns in vertebrates. Two species of the genus Chrysolophus, golden pheasant (Chrysolophus pictus) and Lady Amherst's pheasant (Chrysolophus amherstiae), exhibit brilliant colors in their plumage, but with extreme phenotypic differences, making these two species great models to investigate plumage coloration mechanisms in birds.

Results

We sequenced and assembled a genome of golden pheasant with high coverage and annotated 15,552 protein-coding genes. The genome of Lady Amherst's pheasant is sequenced with low coverage. Based on the feather pigment identification, a series of genomic and transcriptomic comparisons were conducted to investigate the complex features of plumage coloration. By identifying the lineage-specific sequence variations in Chrysolophus and golden pheasant against different backgrounds, we found that four melanogenesis biosynthesis genes and some lipid-related genes might be candidate genomic factors for the evolution of melanin and carotenoid pigmentation, respectively. In addition, a study among 47 birds showed some candidate genes related to carotenoid coloration in a broad range of birds. The transcriptome data further reveal important regulators of the two colorations, particularly one splicing transcript of the microphthalmia-associated transcription factor gene for pheomelanin synthesis.

Conclusions

Analysis of the golden pheasant and its sister pheasant genomes, as well as comparison with other avian genomes, are helpful to reveal the underlying regulation of their plumage coloration. The present study provides important genomic information and insights for further studies of avian plumage evolution and diversity.

Adaptive deletion in resistance gene duplications in the malaria vector Anopheles gambiae

While gene copy-number variations play major roles in long-term evolution, their early dynamics remains largely unknown. However, examples of their role in short-term adaptation are accumulating: identical repetitions of a locus (homogeneous duplications) can provide a quantitative advantage, while the association of differing alleles (heterogeneous duplications) allows carrying two functions simultaneously. Such duplications often result from rearrangements of sometimes relatively large chromosome fragments, and even when adaptive, they can be associated with deleterious side effects that should, however, be reduced by subsequent evolution. Here, we took advantage of the unique model provided by the malaria mosquito Anopheles gambiae s.l. to investigate the early evolution of several duplications, heterogeneous and homogeneous, segregating in natural populations from West Africa. These duplications encompass ~200 kb and 11 genes, including the adaptive insecticide resistance ace-1 locus. Through the survey of several populations from three countries over 3-4 years, we showed that an internal deletion of all coamplified genes except ace-1 is currently spreading in West Africa and introgressing from An. gambiae s.s. to An. coluzzii. Both observations provide evidences of its selection, most likely due to reducing the gene-dosage disturbances caused by the excessive copies of the nonadaptive genes. Our study thus provides a unique example of the early adaptive trajectory of duplications and underlines the role of the environmental conditions (insecticide treatment practices and species ecology). It also emphasizes the striking diversity of adaptive responses in these mosquitoes and reveals a worrisome process of resistance/cost trade-off evolution that could impact the control of malaria vectors in Africa.

The transcriptional landscape of polyploid wheat

The coordinated expression of highly related homoeologous genes in polyploid species underlies the phenotypes of many of the world's major crops. Here we combine extensive gene expression datasets to produce a comprehensive, genome-wide analysis of homoeolog expression patterns in hexaploid bread wheat. Bias in homoeolog expression varies between tissues, with ~30% of wheat homoeologs showing nonbalanced expression. We found expression asymmetries along wheat chromosomes, with homoeologs showing the largest inter-tissue, inter-cultivar, and coding sequence variation, most often located in high-recombination distal ends of chromosomes. These transcriptionally dynamic genes potentially represent the first steps toward neo- or subfunctionalization of wheat homoeologs. Coexpression networks reveal extensive coordination of homoeologs throughout development and, alongside a detailed expression atlas, provide a framework to target candidate genes underpinning agronomic traits in wheat.

Shifting the limits in wheat research and breeding using a fully annotated reference genome

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.

Comparative Transcriptome Analysis of Pseudomonas putida KT2440 Revealed Its Response Mechanisms to Elevated Levels of Zinc Stress

The whole-genome transcriptional response of Pseudomonas putida KT2440 to stress-inducing concentrations of zinc was analyzed in this study by RNA sequencing to thoroughly investigate the bacterial cell response to zinc toxicity. The data revealed that different levels of zinc stress strongly affected the transcription of genes from the following categories: metal transport genes, genes involved in membrane homeostasis, oxidative-stress-responding genes, and genes associated with basic cellular metabolism. At the lowest zinc dose, only several genes associated with metal transport and membrane homeostasis were strongly influenced. At the intermediate zinc dose, transcriptional changes of genes belonging to these two categories were highly pronounced. In addition, the intermediate zinc stress produced high levels of oxidative stress, and influenced amino acid metabolism and respiratory chains of P. putida. At the highest zinc dose, the induction of genes responsible for Fe–S cluster biogenesis was the most remarkable feature. Moreover, upregulation of glyoxylate cycle was observed. In summary, the adaptation of the cell envelope, the maintenance of metal homeostasis and intracellular redox status, and the transcriptional control of metabolism are the main elements of stress response, which facilitates the survival of P. putida KT2440 in zinc-polluted environments.