Is genome size influenced by colonization of new environments in dipteran species?

The chromosome-level assembly of the wild diploid alfalfa genome provides insights into the full landscape of genomic variations between cultivated and wild alfalfa

Alfalfa (Medicago sativa L.) is one of the most important forage legumes in the world, including autotetraploid (M. sativa ssp. sativa) and diploid alfalfa (M. sativa ssp. caerulea, progenitor of autotetraploid alfalfa). Here, we reported a high-quality genome of ZW0012 (diploid alfalfa, 769 Mb, contig N50 = 5.5 Mb), which was grouped into the Northern group in population structure analysis, suggesting that our genome assembly filled a major gap among the members of M. sativa complex. During polyploidization, large phenotypic differences occurred between diploids and tetraploids, and the genetic information underlying its massive phenotypic variations remains largely unexplored. Extensive structural variations (SVs) were identified between ZW0012 and XinJiangDaYe (an autotetraploid alfalfa with released genome). We identified 71 ZW0012-specific PAV genes and 1296 XinJiangDaYe-specific PAV genes, mainly involved in defence response, cell growth, and photosynthesis. We have verified the positive roles of MsNCR1 (a XinJiangDaYe-specific PAV gene) in nodulation using an Agrobacterium rhizobia-mediated transgenic method. We also demonstrated that MsSKIP23_1 and MsFBL23_1 (two XinJiangDaYe-specific PAV genes) regulated leaf size by transient overexpression and virus-induced gene silencing analysis. Our study provides a high-quality reference genome of an important diploid alfalfa germplasm and a valuable resource of variation landscape between diploid and autotetraploid, which will facilitate the functional gene discovery and molecular-based breeding for the cultivars in the future.

Genomic analysis of Nypa fruticans elucidates its intertidal adaptations and early palm evolution

Nypa fruticans (Wurmb), a mangrove palm species with origins dating back to the Late Cretaceous period, is a unique species for investigating long-term adaptation strategies to intertidal environments and the early evolution of palms. Here, we present a chromosome-level genome sequence and assembly for N. fruticans. We integrated the genomes of N. fruticans and other palm family members for a comparative genomic analysis, which confirmed that the common ancestor of all palms experienced a whole-genome duplication event around 89 million years ago, shaping the distinctive characteristics observed in this clade. We also inferred a low mutation rate for the N. fruticans genome, which underwent strong purifying selection and evolved slowly, thus contributing to its stability over a long evolutionary period. Moreover, ancient duplicates were preferentially retained, with critical genes having experienced positive selection, enhancing waterlogging tolerance in N. fruticans. Furthermore, we discovered that the pseudogenization of Early Methionine-labelled 1 (EM1) and EM6 in N. fruticans underly its crypto-vivipary characteristics, reflecting its intertidal adaptation. Our study provides valuable genomic insights into the evolutionary history, genome stability, and adaptive evolution of the mangrove palm. Our results also shed light on the long-term adaptation of this species and contribute to our understanding of the evolutionary dynamics in the palm family.

Comparative genomics reveals the dynamic evolutionary history of cement protein genes of barnacles from intertidal to deep-sea hydrothermal vents

Thoracican barnacles are a diverse group of marine organisms for which the availability of genome assemblies is currently limited. In this study, we sequenced the genomes of two neolepadoid species (Ashinkailepas kermadecensis, Imbricaverruca yamaguchii) from hydrothermal vents, in addition to two intertidal species. Genome sizes ranged from 481 to 1054 Mb, with repetitive sequence contents of 21.2% to 50.7%. Concordance rates of orthologs and heterozygosity rates were between 82.4% and 91.7% and between 1.0% and 2.1%, respectively, indicating high genetic diversity and heterozygosity. Based on phylogenomic analyses, we revised the nomenclature of cement genes encoding cement proteins that are not homologous to any known proteins. The major cement gene, CP100A, was found in all thoracican species, including vent-associated neolepadoids, and was hypothesised to be essential for thoracican settlement. Duplicated genes, CP100B and CP100C, were found only in balanids, suggesting potential functional redundancy or acquisition of new functions associated with the calcareous base. An ancestor of CP52 genes was duplicated dynamically among lepadids, pollicipedids with multiple copies on a single scaffold, and balanids with multiple sequential repeats of the conserved regions, but no CP52 genes were found in neolepadoids, providing insights into cement gene evolution among thoracican lineages. This study enhances our understanding of the adhesion mechanisms of thoracicans in underwater environments. The newly sequenced genomes provide opportunities for studying their evolution and ecology, shedding light on their adaptation to diverse marine environments, and contributing to our knowledge of barnacle biology with valuable genomic resources for further studies in this field.

Compositae-ParaLoss-1272: A complementary sunflower-specific probe set reduces paralogs in phylogenomic analyses of complex systems

Premise

A family-specific probe set for sunflowers, Compositae-1061, enables family-wide phylogenomic studies and investigations at lower taxonomic levels, but may lack resolution at genus to species levels, especially in groups complicated by polyploidy and hybridization.

Methods

We developed a Hyb-Seq probe set, Compositae-ParaLoss-1272, that targets orthologous loci in Asteraceae. We tested its efficiency across the family by simulating target enrichment sequencing in silico. Additionally, we tested its effectiveness at lower taxonomic levels in the historically complex genus Packera. We performed Hyb-Seq with Compositae-ParaLoss-1272 for 19 Packera taxa that were previously studied using Compositae-1061. The resulting sequences from each probe set, plus a combination of both, were used to generate phylogenies, compare topologies, and assess node support.

Results

We report that Compositae-ParaLoss-1272 captured loci across all tested Asteraceae members, had less gene tree discordance, and retained longer loci than Compositae-1061. Most notably, Compositae-ParaLoss-1272 recovered substantially fewer paralogous sequences than Compositae-1061, with only ~5% of the recovered loci reporting as paralogous, compared to ~59% with Compositae-1061.

Discussion

Given the complexity of plant evolutionary histories, assigning orthology for phylogenomic analyses will continue to be challenging. However, we anticipate Compositae-ParaLoss-1272 will provide improved resolution and utility for studies of complex groups and lower taxonomic levels in the sunflower family.

Epigenomic divergence correlates with sequence polymorphism in Arabidopsis paralogs

Processes affecting rates of sequence polymorphism are fundamental to the evolution of gene duplicates. The relationship between gene activity and sequence polymorphism can influence the likelihood that functionally redundant gene copies are co-maintained in stable evolutionary equilibria vs other outcomes such as neofunctionalization. Here, we investigate genic variation in epigenome-associated polymorphism rates in Arabidopsis thaliana and consider whether these affect the evolution of gene duplicates. We compared the frequency of sequence polymorphism and patterns of genetic differentiation between genes classified by exon methylation patterns: unmethylated (unM), gene-body methylated (gbM), and transposon-like methylated (teM) states, which reflect divergence in gene expression. We found that the frequency of polymorphism was higher in teM (transcriptionally repressed, tissue-specific) genes and lower in gbM (active, constitutively expressed) genes. Comparisons of gene duplicates were largely consistent with genome-wide patterns - gene copies that exhibit teM accumulate more variation, evolve faster, and are in chromatin states associated with reduced DNA repair. This relationship between expression, the epigenome, and polymorphism may lead to the breakdown of equilibrium states that would otherwise maintain genetic redundancies. Epigenome-mediated polymorphism rate variation may facilitate the evolution of novel gene functions in duplicate paralogs maintained over evolutionary time.

Ace and ace-like genes of invasive redlegged earth mite: copy number variation, target-site mutations, and their associations with organophosphate insensitivity

BACKGROUND

Invasive Australian populations of redlegged earth mite, Halotydeus destructor (Tucker), are evolving increasing organophosphate resistance. In addition to the canonical ace gene, the target gene of organophosphates, the H. destructor genome contains many radiated ace-like genes that vary in copy number and amino acid sequence. In this work, we characterise copy number and target-site mutation variation at the canonical ace and ace-like genes and test for potential associations with organophosphate insensitivity. This was achieved through comparisons of whole-genome pool-seq data from alive and dead mites following organophosphate exposure.

RESULTS

A combination of increased copy number and target-site mutations at the canonical ace was associated with organophosphate insensitivity in H. destructor. Resistant populations were segregating for G119S, A201S, F331Y at the canonical ace. A subset of populations also had copy numbers of canonical ace > 2, which potentially helps overexpress proteins carrying these target-site mutations. Haplotypes possessing different copy numbers and target-site mutations of the canonical ace gene may be under selection across H. destructor populations. We also detected some evidence that increases in copy number of radiated ace-like genes are associated with organophosphate insensitivity, which might suggest potential roles in sequestration or breakdown of organophosphates.

CONCLUSION

Different combinations of target-site mutations and (or) copy number variation in the canonical ace and ace-like genes may provide non-convergent ways for H. destructor to respond to organophosphate selection. However, these changes may only play a partial role in organophosphate insensitivity, which appears to have a polygenic architecture. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

RANBP2 evolution and human disease

Ran-binding protein 2 (RANBP2)/Nup358 is a nucleoporin and a key component of the nuclear pore complex. Through its multiple functions (e.g., SUMOylation, regulation of nucleocytoplasmic transport) and subcellular localizations (e.g., at the nuclear envelope, kinetochores, annulate lamellae), it is involved in many cellular processes. RANBP2 dysregulation or mutation leads to the development of human pathologies, such as acute necrotizing encephalopathy 1, cancer, neurodegenerative diseases, and it is also involved in viral infections. The chromosomal region containing the RANBP2 gene is highly dynamic, with high structural variation and recombination events that led to the appearance of a gene family called RANBP2 and GCC2 Protein Domains (RGPD), with multiple gene loss/duplication events during ape evolution. Although RGPD homoplasy and maintenance during evolution suggest they might confer an advantage to their hosts, their functions are still unknown and understudied. In this review, we discuss the appearance and importance of RANBP2 in metazoans and its function-related pathologies, caused by an alteration of its expression levels (through promotor activity, post-transcriptional, or post-translational modifications), its localization, or genetic mutations.

Origins and evolution of biological novelty

Understanding the origins and impacts of novel traits has been a perennial interest in many realms of ecology and evolutionary biology. Here, we build on previous evolutionary and philosophical treatments of this subject to encompass novelties across biological scales and eco-evolutionary perspectives. By defining novelties as new features at one biological scale that have emergent effects at other biological scales, we incorporate many forms of novelty that have previously been treated in isolation (such as novelty from genetic mutations, new developmental pathways, new morphological features, and new species). Our perspective is based on the fundamental idea that the emergence of a novelty, at any biological scale, depends on its environmental and genetic context. Through this lens, we outline a broad array of generative mechanisms underlying novelty and highlight how genomic tools are transforming our understanding of the origins of novelty. Lastly, we present several case studies to illustrate how novelties across biological scales and systems can be understood based on common mechanisms of change and their environmental and genetic contexts. Specifically, we highlight how gene duplication contributes to the evolution of new complex structures in visual systems; how genetic exchange in symbiosis alters functions of both host and symbiont, resulting in a novel organism; and how hybridisation between species can generate new species with new niches.

Biogeography and the evolution of acoustic communication in the polyploid North American grey treefrog complex

After polyploid species are formed, interactions between diploid and polyploid lineages may generate additional diversity in novel cytotypes and phenotypes. In anurans, mate choice by acoustic communication is the primary method by which individuals identify their own species and assess suitable mates. As such, the evolution of acoustic signals is an important mechanism for contributing to reproductive isolation and diversification in this group. Here, we estimate the biogeographical history of the North American grey treefrog complex, consisting of the diploid Hyla chrysoscelis and the tetraploid Hyla versicolor, focusing specifically on the geographical origin of whole genome duplication and the expansion of lineages out of glacial refugia. We then test for lineage-specific differences in mating signals by applying comparative methods to a large acoustic data set collected over 52 years that includes >1500 individual frogs. Along with describing the overall biogeographical history and call diversity, we found evidence that the geographical origin of H. versicolor and the formation of the midwestern polyploid lineage are both associated with glacial limits, and that the southwestern polyploid lineage is associated with a shift in acoustic phenotype relative to the diploid lineage with which they share a mitochondrial lineage. In H. chrysoscelis, we see that acoustic signals are largely split by Eastern and Western lineages, but that northward expansion along either side of the Appalachian Mountains is associated with further acoustic diversification. Overall, results of this study provide substantial clarity on the evolution of grey treefrogs as it relates to their biogeography and acoustic communication.

Tracing 100 million years of grass genome evolutionary plasticity

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

Identification and characterization of sperm motility-initiating substance-2 gene in internally fertilizing Cynops species

Sperm motility-initiating substance (SMIS) is an oviductal protein critical for internal fertilization in urodeles. It contributes to the establishment of various reproductive modes in amphibians and is thus a unique research model for the gene evolution of gamete-recognizing ligands that have diversified among animal species. In this study, a paralogous SMIS gene, smis2, was identified via the RNA sequencing of the oviduct of the newt, Cynops pyrrhogaster. The base sequence of the smis2 gene was homologous (˃90%) to that of the original smis gene (smis1), and deduced amino acid sequences of both genes conserved six cysteine residues essential for the cysteine knot motif. Furthermore, smis2 complementary DNA was identified in the oviduct of Cynops ensicauda, and the base substitution patterns also suggested that the smis gene was duplicated in the Salamandridae. Nonsynonymous/synonymous substitution ratios of smis1 and smis2 genes were 0.79 and 2.6, respectively, suggesting that smis2 gene evolution was independently driven by positive selection. Amino acid substitutions were concentrated in the cysteine knot motif of SMIS2. The smis2 gene was expressed in some organs in addition to the oviduct; in contrast, SMIS1 was only expressed in the oviduct. The SMIS2 protein was suggested to be produced and secreted at least in the oviduct and redundantly act in sperm. These results suggest that smis1 plays the original role in the oviduct, whereas smis2 may undergo neofunctionalization, which rarely occurs in gene evolution.

The blackgrass genome reveals patterns of non-parallel evolution of polygenic herbicide resistance

Globally, weedy plants are a major constraint to sustainable crop production. Much of the success of weeds rests with their ability to rapidly adapt in the face of human-mediated management of agroecosystems. Alopecurus myosuroides (blackgrass) is a widespread and impactful weed affecting agriculture in Europe. Here we report a chromosome-scale genome assembly of blackgrass and use this reference genome to explore the genomic/genetic basis of non-target site herbicide resistance (NTSR). Based on our analysis of F2 seed families derived from two distinct blackgrass populations with the same NTSR phenotype, we demonstrate that the trait is polygenic and evolves from standing genetic variation. We present evidence that selection for NTSR has signatures of both parallel and non-parallel evolution. There are parallel and non-parallel changes at the transcriptional level of several stress- and defence-responsive gene families. At the genomic level, however, the genetic loci underpinning NTSR are different (non-parallel) between seed families. We speculate that variation in the number, regulation and function of stress- and defence-related gene families enable weedy species to rapidly evolve NTSR via exaptation of genes within large multi-functional gene families. These results provide novel insights into the potential for, and nature of plant adaptation in rapidly changing environments.

Where whole-genome duplication is most beneficial: Adaptation of mangroves to a wide salinity range between land and sea

Whole-genome duplication (WGD) is believed to increase the chance of adaptation to a new environment. This conjecture may apply particularly well to new environments that are not only different but also more variable than ancestral habitats. One such prominent environment is the interface between land and sea, which has been invaded by woody plants, collectively referred as mangroves, multiple times. Here, we use two distantly related mangrove species (Avicennia marina and Rhizophora apiculata) to explore the effects of WGD on the adaptive process. We found that a high proportion of duplicated genes retained after WGD have acquired derived differential expression in response to salt gradient treatment. The WGD duplicates differentially expressed in at least one copy usually (>90%) diverge from their paralogues' expression profiles. Furthermore, both species evolved in parallel to have one paralogue expressed at a high level in both fresh water and hypersaline conditions but at a lower level at medium salinity. The pattern contrasts with the conventional view of monotone increase/decrease as salinity increases. Differentially expressed copies have thus probably acquired a new role in salinity tolerance. Our results indicate that the WGD duplicates may have evolved to function collaboratively in coping with different salinity levels, rather than specializing in the intermediate salinity optimal for mangrove plants. In conclusion, WGD and the retained duplicates appear to be an effective solution for adaptation to new and unstable environments.

Identification of candidate genes in cotton associated with specific seed traits and their initial functional characterization in Arabidopsis

Oilseed crops are used to produce vegetable oil to satisfy the requirements of humans and livestock. Cotton (Gossypium spp.) is of great economic value because it is used as both an important textile commodity and a nutrient-rich resource. Cottonseed oil is rich in polyunsaturated fatty acids and does not contain trans fatty acids; hence, it is considered a healthy vegetable oil. However, research on the genetic basis for cottonseed protein content, oil production, and fatty acid composition is lacking. Here, we investigated the protein content, oil content, and fatty acid composition in terms of oleic acid (C18:1) and linoleic acid (C18:2) in mature cottonseeds from 318 Gossypium hirsutum accessions. Moreover, we examined the dynamic change of protein content and lipid composition including palmitic acid (C16:0), stearic acid (C18:0), oleic acid (C18:1), linoleic acid (C18:2), and linolenic acid (C18:3) in developing seeds from 258 accessions at 10 and 20 days post-anthesis. Then, we conducted a genome-wide association study and identified 152 trait-associated loci and 64 candidate genes responsible for protein and oil-related contents in mature cottonseeds and ovules. Finally, six candidate genes were experimentally validated to be involved in the regulation of fatty acid biosynthesis through heterologous expression in Arabidopsis. These results comprise a solid foundation for expanding our understanding of lipid biosynthesis in cotton, which will help breeders manipulate protein and oil contents to make it a fully developed 'fiber, food, and oil crop'.

Retrieving a disrupted gene encoding phospholipase A for fibre enhancement in allotetraploid cultivated cotton

After polyploidization originated from one interspecific hybridization event in Gossypium, Gossypium barbadense evolved to produce extra-long staple fibres than Gossypium hirsutum (Upland cotton), which produces a higher fibre yield. The genomic diversity between G. barbadense and G. hirsutum thus provides a genetic basis for fibre trait variation. Recently, rapid accumulation of gene disruption or deleterious mutation was reported in allotetraploid cotton genomes, with unknown impacts on fibre traits. Here, we identified gene disruptions in allotetraploid G. hirsutum (18.14%) and G. barbadense (17.38%) through comparison with their presumed diploid progenitors. Relative to conserved genes, these disrupted genes exhibited faster evolution rate, lower expression level and altered gene co-expression networks. Within a module regulating fibre elongation, a hub gene experienced gene disruption in G. hirsutum after polyploidization, with a 2-bp deletion in the coding region of GhNPLA1D introducing early termination of translation. This deletion was observed in all of the 34 G. hirsutum landraces and 36 G. hirsutum cultivars, but not in 96% of 57 G. barbadense accessions. Retrieving the disrupted gene GhNPLA1D using its homoeolog GhNPLA1A achieved longer fibre length in G. hirsutum. Further enzyme activity and lipids analysis confirmed that GhNPLA1A encodes a typical phospholipase A and promotes cotton fibre elongation via elevating intracellular levels of linolenic acid and 34:3 phosphatidylinositol. Our work opens a strategy for identifying disrupted genes and retrieving their functions in ways that can provide valuable resources for accelerating fibre trait enhancement in cotton breeding.

The chromosome-level genome provides insight into the molecular mechanism underlying the tortuous-branch phenotype of Prunus mume

Plant with naturally twisted branches is referred to as a tortuous-branch plant, which have extremely high ornamental value due to their zigzag shape and the natural twisting of their branches. Prunus mume is an important woody ornamental plant. However, the molecular mechanism underlying this unique trait in Prunus genus is unknown. Here, we present a chromosome-level genome assembly of the cultivated P. mume var. tortuosa created using Oxford Nanopore combined with Hi-C scaffolding, which resulted in a 237.8 Mb genome assembly being anchored onto eight pseudochromosomes. Molecular dating indicated that P. mume is the most recently differentiated species in Prunus. Genes associated with cell division, development and plant hormones play essential roles in the formation of tortuous branch trait. A putative regulatory pathway for the tortuous branch trait was constructed based on gene expression levels. Furthermore, after transferring candidate PmCYCD genes into Arabidopsis thaliana, we found that seedlings overexpressing these genes exhibited curled rosette leaves. Our results provide insights into the evolutionary history of recently differentiated species in Prunus genus, the molecular basis of stem morphology, and the molecular mechanism underlying the tortuous branch trait and highlight the utility of multi-omics in deciphering the properties of P. mume plant architecture.

A lineage-specific arginine in POS1 is required for fruit size control in Physaleae (Solanaceae) via gene co-option

Solanaceae have important economic value mainly due to their edible fruits. Physalis organ size 1/cytokinin response factor 3 (POS1/CRF3), a unique gene in Solanaceae, is involved in fruit size variation in Physalis but not in Solanum. However, the underlying mechanisms remain elusive. Here, we found that POS1/CRF3 was likely created via the fusion of CRF7 and CRF8 duplicates. Multiple genetic manipulations revealed that only POS1 and Capsicum POS1 (CaPOS1) functioned in fruit size control via the positive regulation of cell expansion. Comparative studies in a phylogenetic framework showed the directional enhancement of POS1-like expression in the flowers and fruits of Physaleae and the specific gain of certain interacting proteins associated with cell expansion by POS1 and CaPOS1. A lineage-specific single nucleotide polymorphism (SNP) caused the 68th amino acid histidine in the POS1 orthologs of non-Physaleae (Nicotiana and Solanum) to change to arginine in Physaleae (Physalis and Capsicum). Substituting the arginine in Physaleae POS1-like by histidine completely abolished their function in the fruits and the protein-protein interaction (PPI) with calreticulin-3. Transcriptomic comparison revealed the potential downstream pathways of POS1, including the brassinosteroid biosynthesis pathway. However, POS1-like may have functioned ancestrally in abiotic stress within Solanaceae. Our work demonstrated that heterometric expression and a SNP caused a single amino acid change to establish new PPIs, which contributed to the co-option of POS1 in multiple regulatory pathways to regulate cell expansion and thus fruit size in Physaleae. These results provide new insights into fruit morphological evolution and fruit yield control.

Catalogue of the Diptera (Insecta) of Morocco— an annotated checklist, with distributions and a bibliography

The faunistic knowledge of the Diptera of Morocco recorded from 1787 to 2021 is summarized and updated in this first catalogue of Moroccan Diptera species. A total of 3057 species, classified into 948 genera and 93 families (21 Nematocera and 72 Brachycera), are listed. Taxa (superfamily, family, genus and species) have been updated according to current interpretations, based on reviews in the literature, the expertise of authors and contributors, and recently conducted fieldwork. Data to compile this catalogue were primarily gathered from the literature. In total, 1225 references were consulted and some information was also obtained from online databases. Each family was reviewed and the checklist updated by the respective taxon expert(s), including the number of species that can be expected for that family in Morocco. For each valid species, synonyms known to have been used for published records from Morocco are listed under the currently accepted name. Where available, distribution within Morocco is also included. One new combination is proposed: Assuaniamelanoleuca (Séguy, 1941), comb. nov. (Chloropidae).

Silencing of the BtTPS genes by transgenic plant-mediated RNAi to control Bemisia tabaci MED

BACKGROUND

Whitefly (Bemisia tabaci) is a typical pest that causes severe damage to hundreds of agricultural crops. The trehalose-6-phosphate synthase (TPS) genes, as the key genes in the insect trehalose synthesis pathway, are important for insect growth and development. The whitefly TPS genes may be a main reason for the severe damage and may represent potential targets for the control of whiteflies.

RESULTS

In this study, we identified and cloned three TPS genes from B. tabaci MED and found that the BtTPS1 and BtTPS2 genes showed higher expression levels than the BtTPS3 gene. Then, RNA interference (RNAi) of BtTPS1 and BtTPS2 resulted in significant mortality and influenced the expression of related genes involved in energy metabolism and chitin biosynthesis in whitefly adults. Finally, the transgenic tobacco plants showed a significant effect on B. tabaci, and knockdown of BtTPS1 or BtTPS2 led to retarded growth and low hatchability in whitefly nymphs, and caused 90% mortality and decreased the fecundity in whitefly adults. Additionally, the transgenic tobacco with combinatorial RNAi of BtTPS1 and BtTPS2 showed a better efficacy against whiteflies than individual silencing.

CONCLUSION

Our results suggest that silencing of the BtTPS genes can compromise the growth and development of whiteflies, offering not only a new option for whitefly control but also a secure and environmentally friendly management strategy.

Genome-wide identification and transcriptome-based expression profiling of the Sox gene family in the spinyhead croaker (Collichthys lucidus)

Sox genes encode transcription factors with a high-mobility group (HMG) box, playing critical roles in the initiation and maintenance of a variety of developmental processes, such as sex determination and differentiation. In the present study, we identified 26 Sox genes in the genome of spinyhead croaker Collichthys lucidus (Richardson, 1844) with homology-based analysis of the HMG box. The transcriptome-based expression profiles revealed that the expression of the Sox gene in gonads began to differ between sexes when the body length was 2.74 ± 0.24 cm. At that time, three Sox genes (Sox11b, Sox8a and Sox19) were significantly upregulated, accompanied by the downregulation of 12 Sox genes in the ovary, and six Sox genes were temporarily significantly upregulated in the testis. Afterwards, the expression profile of Sox genes changed only with a small amplitude in both the ovary and testis. For adult tissues, huge differences were observed in the expression profiles of Sox genes between ovaries and testes, as well as small differences in somatic tissues between sexes. These results provide clues to further decipher the role of Sox genes in the processes of sex determination and differentiation in spinyhead croaker and other teleosts.

The genome sequence provides insights into salt tolerance of Achnatherum splendens (Gramineae), a constructive species of alkaline grassland

Achnatherum splendens Trin. (Gramineae) is a constructive species of the arid grassland ecosystem in Northwest China and is a major forage grass. It has good tolerance of salt and drought stress in alkaline habitats. Here, we report its chromosome-level genome, determined through a combination of Illumina HiSeq sequencing, PacBio sequencing and Hi-C technology. The final assembly of the ~1.17 Gb genome sequence had a super-scaffold N50 of 40.3 Mb. A total of 57 374 protein-coding genes were annotated, of which 54 426 (94.5%) genes have functional protein annotations. Approximately 735 Mb (62.37%) of the assembly were identified as repetitive elements, and among these, LTRs (40.53%) constitute the highest proportion, having made a major contribution to the expansion of genome size in A. splendens. Phylogenetic analysis revealed that A. splendens diverged from the Brachypodium distachyon-Hordeum vulgare-Aegilops tauschii subclade around 37 million years ago (Ma) and that a clade comprising these four species diverged from the Phyllostachys edulis clade ~47 Ma. Genomic synteny indicates that A. splendens underwent an additional species-specific whole-genome duplication (WGD) 18-20 Ma, which further promoted an increase in copies of numerous saline-alkali-related gene families in the A. splendens genome. By transcriptomic analysis, we further found that many of these duplicated genes from this extra WGD exhibited distinct functional divergence in response to salt stress. This WGD, therefore, contributed to the strong resistance to salt stress and widespread arid adaptation of A. splendens.

Evolutionary and functional analysis of two-component system in chickpea reveals CaRR13, a TypeB RR, as positive regulator of symbiosis

The critical role of cytokinin in early nodulation in legumes is well known. In our study, exogenous cytokinin application to roots of the important crop legume, chickpea (Cicer arietinum L.), led to the formation of pseudo-nodules even in the absence of rhizobia. Hence, a genome-wide analysis of the cytokinin signalling, two-component system (TCS) genes, was conducted in chickpea, Medicago and Cajanus cajan. The integrated phylogenetic, evolutionary and expression analysis of the TCS genes was carried out, which revealed that histidine kinases (HKs) were highly conserved, whereas there was diversification leading to neofunctionalization at the level of response regulators (RRs) especially the TypeB RRs. Further, the functional role of the CaHKs in nodulation was established by complementation of the sln1Δ mutant of yeast and cre1 mutants of (Medicago) which led to restoration of the nodule-deficient phenotype. Additionally, the highest expressing TypeB RR of chickpea, CaRR13, was functionally characterized. Its localization in the nucleus and its Y1H assay-based interaction with the promoter of the early nodulation gene CaNSP2 indicated its role as a transcription factor regulating early nodulation. Overexpression, RNAi lines and complementation of cre1 mutants with CaRR13 revealed its critical involvement as an important signalling molecule regulating early events of nodule organogenesis in chickpea.

Ancient rapid functional differentiation and fixation of the duplicated members in rice Dof genes after whole genome duplication

Whole genome duplication (WGD) in plants is typically followed by genomic downsizing, where large portions of the new genome are lost. Whether this downsizing is accompanied by increased or decreased evolutionary rates of the remaining genes is poorly known, not least because homeolog pairings are often obscured by chromosomal rearrangement. Here, we use the newly published genome from a sedge, namely Kobresia littledalei, and CRISPR/Cas-9 editing to investigate how the Rho WGD event 70 million years ago (MYA) affected transcription factor evolutionary rates, fates, and function in rice (Oryza sativa) and sorghum (Sorghum bicolor). We focus on the 30-member DNA-binding with one zinc finger (Dof) transcription factor family in both crops due to their agronomic importance. Using the known speciation dates of rice from Kobresia (97 MYA) and sorghum (50 MYA), we find that rates of amino acid substitution in the critical Dof domain region were over twofold higher during the 20-million-year period following the WGD than before or afterward. Through comparison of synteny blocks, we report that at least 11% of Dof genes were purged from 70 to 50 MYA, while only 6% have been lost in the most recent 50-million-year interval. CRISPR/Cas9 editing revealed widespread fitness-related defects in flowering and lack of redundancy of paired members, as well as significant differences in expression between gene pairs. Together these findings demonstrate the strength of Dof genes as a model for deep evolutionary study and offer one of the most detailed portraits yet of the Rho WGD impact on a gene lineage.

Evolutionary and expression dynamics of LRR-RLKs and functional establishment of KLAVIER homolog in shoot mediated regulation of AON in chickpea symbiosis

Chickpea shoot exogenously treated with cytokinin showed stunted phenotype of root, shoot and significantly reduced nodule numbers. Genome-wide identification of LRR-RLKs in chickpea and Medicago resulted in 200 and 371 genes respectively. Gene duplication analysis revealed that LRR-RLKs family expanded through segmental duplications in chickpea and tandem duplications in Medicago. Expression profiling of LRR-RLKs revealed their involvement in cytokinin signaling and plant organ development. Overexpression of KLAVIER ortholog of chickpea, Ca_LRR-RLK147, in roots revealed its localization in the membrane but showed no effect on root nodulation despite increased cle peptide levels. Two findings (i) drastic effect on nodule number by exogenous cytokinin treatment to only shoot and restoration to normal nodulation by treatment to both root and shoot tissue and (ii) no effect on nodule number by overexpression of Ca_LRR-RLK147 establishes the fact that despite presence of cle peptides in root, the function of Ca_LRR-RLK147 was shoot mediated during AON.

Nonadaptive molecular evolution of seminal fluid proteins in Drosophila

Seminal fluid proteins (SFPs) are a group of reproductive proteins that are among the most evolutionarily divergent known. As SFPs can impact male and female fitness, these proteins have been proposed to evolve under postcopulatory sexual selection (PCSS). However, the fast change of the SFPs can also result from nonadaptive evolution, and the extent to which selective constraints prevent SFPs rapid evolution remains unknown. Using intra‐ and interspecific sequence information, along with genomics and functional data, we examine the molecular evolution of approximately 300 SFPs in Drosophila. We found that 50–57% of the SFP genes, depending on the population examined, are evolving under relaxed selection. Only 7–12% showed evidence of positive selection, with no evidence supporting other forms of PCSS, and 35–37% of the SFP genes were selectively constrained. Further, despite associations of positive selection with gene location on the X chromosome and protease activity, the analysis of additional genomic and functional features revealed their lack of influence on SFPs evolving under positive selection. Our results highlight a lack of sufficient evidence to claim that most SFPs are driven to evolve rapidly by PCSS while identifying genomic and functional attributes that influence different modes of SFPs evolution.

Comparative analysis of the superoxide dismutase gene family in Cetartiodactyla

Cetacea, whales, dolphins and porpoises form an order of mammals adapted to aquatic life. Their transition to an aquatic habitat resulted in exceptional protection against cellular insults, including oxidative and osmotic stress. Here, we considered the structure and molecular evolution of the superoxide dismutase (SOD) gene family, which encodes essential enzymes in the mammalian antioxidant system, in the superorder Cetartiodactyla. To this end, we juxtaposed cetaceans and their closest extant relatives (order Artiodactyla). We identified 94 genes in 23 species, of which 70 are bona fide intact genes. Although the SOD gene family is conserved in Cetartiodactyla, lineage-specific gene duplications and deletions were observed. Phylogenetic analyses show that the SOD2 subfamily diverged from a clade containing SOD1 and SOD3, suggesting that cytoplasmic, extracellular and mitochondrial SODs have started down independent evolutionary paths. Specific-amino acid changes (e.g. K130N in SOD2) that may enhance ROS elimination were identified in cetaceans. In silico analysis suggests that the core transcription factor repertoire of cetartiodactyl SOD genes may include Sp1, NF-κB, Nrf2 and AHR. Putative transcription factors binding sites responding to hypoxia were (e.g. Suppressor of Hairless; Su(H)) found in the cetacean SOD1 gene. We found significant evidence for positive selection in cetaceans using codon models. Cetaceans with different diving abilities also show divergent evolution of SOD1 and SOD2. Our genome-wide analysis of SOD genes helps clarify their relationship and evolutionary trajectory and identify putative functional changes in cetaceans.

The worldwide invasion of Drosophila suzukii is accompanied by a large increase of transposable element load and a small number of putatively adaptive insertions

Transposable Elements (TEs) are ubiquitous and mobile repeated sequences. They are major determinants of host fitness. Here, we characterized the TE content of the spotted wing fly Drosophila suzukii. Using a recently improved genome assembly, we reconstructed TE sequences de novo, and found that TEs occupy 47% of the genome and are mostly located in gene poor regions. The majority of TE insertions segregate at low frequencies, indicating a recent and probably ongoing TE activity. To explore TE dynamics in the context of biological invasions, we studied variation of TE abundance in genomic data from 16 invasive and six native populations of D. suzukii. We found a large increase of the TE load in invasive populations correlated with a reduced Watterson estimate of genetic diversity θ̂w a proxy of effective population size. We did not find any correlation between TE contents and bioclimatic variables, indicating a minor effect of environmentally induced TE activity. A genome-wide association study revealed that ca. 2,000 genomic regions are associated with TE abundance. We did not find, however, any evidence in such regions of an enrichment for genes known to interact with TE activity (e.g. transcription factor encoding genes or genes of the piRNA pathway). Finally, the study of TE insertion frequencies revealed 15 putatively adaptive TE insertions, six of them being likely associated with the recent invasion history of the species.

Evolution of rapid blue-light response linked to explosive diversification of ferns in angiosperm forests

Ferns appear in the fossil record some 200 Myr before angiosperms. However, as angiosperm-dominated forest canopies emerged in the Cretaceous period there was an explosive diversification of modern (leptosporangiate) ferns, which thrived in low, blue-enhanced light beneath angiosperm canopies. A mechanistic explanation for this transformative event in the diversification of ferns has remained elusive. We used physiological assays, transcriptome analysis and evolutionary bioinformatics to investigate a potential connection between the evolution of enhanced stomatal sensitivity to blue light in modern ferns and the rise of angiosperm-dominated forests in the geological record. We demonstrate that members of the largest subclade of leptosporangiate ferns, Polypodiales, have significantly faster stomatal response to blue light than more ancient fern lineages and a representative angiosperm. We link this higher sensitivity to levels of differentially expressed genes in blue-light signaling, particularly in the cryptochrome (CRY) signaling pathway. Moreover, CRYs of the Polypodiales examined show gene duplication events between 212.9-196.9 and 164.4-151.8 Ma, when angiosperms were emerging, which are lacking in other major clades of extant land plants. These findings suggest that evolution of stomatal blue-light sensitivity helped modern ferns exploit the shady habitat beneath angiosperm forest canopies, fueling their Cretaceous hyperdiversification.

Comprehensive analysis of structural, functional, and evolutionary dynamics of Leucine Rich Repeats-RLKs in Thinopyrum elongatum

Leucine Rich Repeats-receptor-like protein kinases (LRR-RLKs) regulate several critical biological processes ranging from growth and development to stress response. Thinopyrum elongatum harbours many desirable traits such as biotic and abiotic stress resistance and therefore commonly used by wheat breeders. In the present investigation, in-silico analysis of LRR-RLKs yielded 589 genes of which 431 were membrane surface RLKs and 158 were receptor like cytoplasmic kinases. An insight into the gene and protein structure revealed quite a conserved nature of these proteins within subgroups. A large expansion in LRR-RLKs was due to tandem and segmental duplication event. Maximum number of tandem and segmentally duplicated pairs was observed in LRR-VI and LRR-XII subfamily, respectively. Furthermore, syntenic analyses revealed that chromosome 6 harboured more (48) tandem duplicated genes while chromosome 7 possessed more (47) segmentally duplicated genes. A detailed analysis about the gene duplication events coupled with expression profiles during Fusarium graminearum infection and water deficiency unravelled the expansion of the gene family with sub functionalization and neofunctionalization. Interaction network analysis showed that LRR-RLKs can heterodimerize upon ligand binding to perform various plant functional attributes.

Genes derived from ancient polyploidy have higher genetic diversity and are associated with domestication in Brassica rapa

Many crops are polyploid or have a polyploid ancestry. Recent phylogenetic analyses have found that polyploidy often preceded the domestication of crop plants. One explanation for this observation is that increased genetic diversity following polyploidy may have been important during the strong artificial selection that occurs during domestication. In order to test the connection between domestication and polyploidy, we identified and examined candidate genes associated with the domestication of the diverse crop varieties of Brassica rapa. Like all 'diploid' flowering plants, B. rapa has a diploidized paleopolyploid genome and experienced many rounds of whole genome duplication (WGD). We analyzed transcriptome data of more than 100 cultivated B. rapa accessions. Using a combination of approaches, we identified > 3000 candidate genes associated with the domestication of four major B. rapa crop varieties. Consistent with our expectation, we found that the candidate genes were significantly enriched with genes derived from the Brassiceae mesohexaploidy. We also observed that paleologs were significantly more diverse than non-paleologs. Our analyses find evidence for that genetic diversity derived from ancient polyploidy played a key role in the domestication of B. rapa and provide support for its importance in the success of modern agriculture.

The soybean plasma membrane-localized cation/H+ exchanger GmCHX20a plays a negative role under salt stress

Cation/H+ -exchanger (CHX) perform diverse functions in plants, including being a part of the protective mechanisms to cope with salt stress. GmCHX1 has been previously identified as the causal gene in a major salt-tolerance quantitative trait locus (QTL) in soybean, but little is known about another close paralog, GmCHX20a, found in the same QTL. In this study, GmCHX20a was characterized along with GmCHX1. The expression patterns of the two genes and the direction of Na+ flux directed by overexpression of these two transporters are different, suggesting that they are functionally distinct. The ectopic expression of GmCHX20a led to an increase in salt sensitivity and osmotic tolerance, which was consistent with its role in increasing Na+ uptake into the root. Although this seems counter-intuitive, it may in fact be part of the mechanism by which soybean could counter act the effects of osmotic stress, which is commonly manifested in the initial stage of salinity stress. On the other hand, GmCHX1 from salt-tolerant soybean was shown to protect plants via Na+ exclusion under salt stress. Taken together these results suggest that GmCHX20a and GmCHX1 might work complementally through a concerted effort to address both osmotic stress and ionic stress as a result of elevated salinity.

Genome-wide identification and expression analysis of the NF-Y transcription factor family in Populus

In the past few years, many studies have reported that the transcription factor Nuclear Factor Y (NF-Y) gene family plays important roles in embryonic development, photosynthesis, flowering time regulation and stress response, in various plants. Although the NF-Y gene family has been systematically studied in many species, little is known about NF-Y genes in Populus. In this study, the NF-Y gene family in the Populus genome was identified and its structural characteristics were described. Fifty-two NF-Y genes were authenticated in the Populus trichocarpa genome and categorized into three subfamilies (NF-YA/B/C) by phylogenetic analysis. Chromosomal localization of these genes revealed that they were distributed randomly across 17 of the 19 chromosomes. Segmental duplication played a vital role in the amplification of Populus NF-Y gene family. Moreover, microsynteny analysis indicated that, among Populus trichocarpa, Arabidopsis thaliana, Vitis vinifera and Carica papaya, NF-Y duplicated regions were more conserved between Populus trichocarpa and Vitis vinifera. Redundant stress-related cis-elements were also found in the promoters of most 13 NF-YA genes and their expression levels varied widely following drought, salt, ABA and cold treatments. Subcellular localization experiments in tobacco showed that PtNF-YA3 was localized in nucleus and cytomembrane, while PtNF-YA4 was only in the nucleus in tobacco. According to the transcriptional activity experiments, neither of them had transcriptional activity in yeast. In summary, a comprehensive analysis of the Populus NF-Y gene family was performed to establish a theoretical basis for further functional studies on this family.

Characterisation of Medicago truncatula CLE34 and CLE35 in nitrate and rhizobia regulation of nodulation

Legumes form a symbiosis with atmospheric nitrogen (N₂ )-fixing soil rhizobia, resulting in new root organs called nodules that enable N₂ -fixation. Nodulation is a costly process that is tightly regulated by the host through autoregulation of nodulation (AON) and nitrate-dependent regulation of nodulation. Both pathways require legume-specific CLAVATA/ESR-related (CLE) peptides. Nitrogen-induced nodulation-suppressing CLE peptides have not previously been investigated in Medicago truncatula, for which only rhizobia-induced MtCLE12 and MtCLE13 have been characterised. Here, we report on novel peptides MtCLE34 and MtCLE35 in nodulation control. The nodulation-suppressing CLE peptides of five legume species were classified into three clades based on sequence homology and phylogeny. This approached identified MtCLE34 and MtCLE35 and four new CLE peptide orthologues of Pisum sativum. Whereas MtCLE12 and MtCLE13 are induced by rhizobia, MtCLE34 and MtCLE35 respond to both rhizobia and nitrate. MtCLE34 was identified as a pseudogene lacking a functional CLE-domain. MtCLE35 was found to inhibit nodulation in a SUNN- and RDN1-dependent manner via overexpression analysis. Together, our findings indicate that MtCLE12 and MtCLE13 have a specific role in AON, while MtCLE35 regulates nodule numbers in response to both rhizobia and nitrate. MtCLE34 likely had a similar role to MtCLE35, but its function was lost due to a premature nonsense mutation.

Copy number variations and young duplicate genes have high methylation levels in sticklebacks

Gene duplication is an important driver of genomic diversity that can promote adaptive evolution. However, like most mutations, a newly duplicated gene is often deleterious and removed from the genome by drift or natural selection. The early molecular changes that occur soon after duplication therefore may influence the long-term survival of gene duplicates, but relatively little empirical data exist on the events near the onset of duplication before mutations have time to accumulate. In this study, we contrast gene expression and DNA methylation levels of duplicate genes in the threespine stickleback, Gasterosteus aculeatus, including recently emerged duplications that segregate as copy number variations (CNVs). We find that younger duplicate genes have higher levels of promoter methylation than older genes, and that gene CNVs have higher promoter methylation than non-CNVs. These results suggest preferential duplication of highly methylated genes or rapid methylation changes soon after duplication. We also find a negative association between methylation and expression, providing a putative role for methylation in suppressing transcription that compensates for increases in gene copy numbers and promoting paralog retention. We propose that methylation contributes to the longevity of young duplicate genes, extending the window of opportunity for functional divergence via mutation.

Genome evolution in Oryza allopolyploids of various ages: Insights into the process of diploidization

The prevalence and recurrence of whole-genome duplication in plants and its major role in evolution have been well recognized. Despite great efforts, many aspects of genome evolution, particularly the temporal progression of genomic responses to allopolyploidy and the underlying mechanisms, remain poorly understood. The rice genus Oryza consists of both recently formed and older allopolyploid species, representing an attractive system for studying the genome evolution after allopolyploidy. In this study, through screening BAC libraries and sequencing and annotating the targeted BAC clones, we generated orthologous genomic sequences surrounding the DEP1 locus, a major grain yield QTL in cultivated rice, from four Oryza polyploids of various ages and their likely diploid genome donors or close relatives. Based on sequenced DEP1 region and published data from three other genomic regions, we investigated the temporal evolutionary dynamics of four polyploid genomes at both genetic and expression levels. In the recently formed BBCC polyploid, Oryza minuta, genome dominance was not observed and its short-term responses to allopolyploidy are mainly manifested as a high proportion of homoeologous gene pairs showing unequal expression. This could partly be explained by parental legacy, rewiring of divergent regulatory networks and epigenetic modulation. Moreover, we detected an ongoing diploidization process in this genus, and suggest that the expression divergence driven by changes of selective constraint probably plays a big role in the long-term diploidization. These findings add novel insights into our understanding of genome evolution after allopolyploidy, and could facilitate crop improvements through hybridization and polyploidization.

Evolution of the Brassicaceae-specific MS5-Like family and neofunctionalization of the novel MALE STERILITY 5 gene essential for male fertility in Brassica napus

New genes (or lineage-specific genes) can facilitate functional innovations. MALE STERILITY 5 (MS5) in Brassica napus is a fertility-related new gene, which has two wild-type alleles (BnMS5a and BnMS5c ) and two mutant alleles (BnMS5b and BnMS5d ) that could induce male sterility. Here, we studied the history and functional evolution of MS5 homologs in plants by phylogenetic analysis and molecular genetic experiments. We identified 727 MS5 homologs and found that they define a Brassicaceae-specific gene family that has expanded partly via multiple tandem gene duplications and also probably transpositions. The MS5 in B. napus is inherited from a basic diploid ancestor of B. rapa. Molecular genetic experiments indicate that BnMS5a and BnMS5c are functionally distinct in B. napus and that BnMS5d can inhibit BnMS5a in B. napus in a dosage-dependent manner. The BnMS5a protein can move in coordination with meiotic telomeres and interact with the nuclear envelope protein SUN1, with a possible crucial role in meiotic chromosome behavior. In summary, BnMS5 belongs to a Brassicaceae-specific new gene family, and has gained a novel function that is essential for male fertility in B. napus through neofunctionalization that has likely occurred since the origin of B. rapa.

Genetics of autoimmunity in plants: an evolutionary genetics perspective

Autoimmunity in plants has been found in numerous hybrids as a form of hybrid necrosis and mutant panels. Uncontrolled cell death is a main cellular outcome of autoimmunity, which negatively impacts growth. Its occurrence highlights the vulnerable nature of the plant immune system. Genetic investigation of autoimmunity in hybrid plants revealed that extreme variation in the immune receptor repertoire is a major contributor, reflecting an evolutionary conundrum that plants face in nature. In this review, we discuss natural variation in the plant immune system and its contribution to fitness. The value of autoimmunity genetics lies in its ability to identify combinations of a natural immune receptor and its partner that are predisposed to triggering autoimmunity. The network of immune components for autoimmunity becomes instrumental in revealing mechanistic details of how immune receptors recognize cellular invasion and activate signaling. The list of autoimmunity-risk variants also allows us to infer evolutionary processes contributing to their maintenance in the natural population. Our approach to autoimmunity, which integrates mechanistic understanding and evolutionary genetics, has the potential to serve as a prognosis tool to optimize immunity in crops.

The effect of pelage, background, and distance on predator detection and the evolution of primate color vision

Primates' predators, such as carnivore mammals, usually rely on camouflage to increase proximity to prey and start a predatory attempt. Camouflage depends not only on the color pattern presented by a predator's pelage but also on the background scene in which the hunting takes place. Another factor that influences camouflage effectiveness is prey's color vision since a given camouflage strategy might not work for all visual phenotypes. Still, little research has been made on the effect of primate visual phenotype on predator detection. Here, we investigate the effects of natural pelages, background scenarios, visual phenotypes, and detection distances on predator detection. We used photographs of taxidermized carnivores (ocelots, cougars, and lesser grisons) as detection stimuli, taken in three different natural scenarios (forest, savanna, and grassland), and at two viewing distances (near and far). On a touchscreen monitor, sets of four images (only one containing a hidden animal) were randomly presented to 39 human males (19 dichromats and 20 trichromats). We found that trichromats, when compared to dichromats, present a lower latency and a higher accuracy of carnivore detection for some conditions tested. We also found that pelage color, background scenario, and detection distance interact to influence the effectiveness of camouflage. Our results suggest that trichromacy might be even more advantageous for carnivore detection than thought before, since it facilitates detection of mammals with diverse pelage colorations, in environments with different phytophysiognomies, and at longer distances. We also propose that the higher rates of dichromacy found in modern human societies could have resulted from a relaxation in predation.

Oxidative stress drives divergent evolution of the glutathione peroxidase (GPX) gene family in mammals

The molecular basis for adaptations to extreme environments can now be understood by interrogating the ever-increasing number of sequenced genomes. Mammals such as cetaceans, bats, and highland species can protect themselves from oxidative stress, a disruption in the balance of reactive oxygen species, which results in oxidative injury and cell damage. Here, we consider the evolution of the glutathione peroxidase (GPX) family of antioxidant enzymes by interrogating publicly available genome data from 70 mammalian species from all major clades. We identified 8 GPX subclasses ubiquitous to all mammalian groups. Mammalian GPX gene families resolved into the GPX4/7/8 and GPX1/2/3/5/6 groups and are characterized by several instances of gene duplication and loss, indicating a dynamic process of gene birth and death in mammals. Seven of the eight GPX subfamilies (all but GPX7) were under positive selection, with the residues under selection located at or close to active sites or at the dimer interface. We also reveal evidence of a correlation between ecological niches (e.g. high oxidative stress) and the divergent selection and gene copy number of GPX subclasses. Notably, a convergent expansion of GPX1 was observed in several independent lineages of mammals under oxidative stress and may be important for avoiding oxidative damage. Collectively, this study suggests that the GPX gene family has shaped the adaption of mammals to stressful environments.

Genome-scale molecular principles of mRNA half-life regulation in yeast

Precise control of protein and messenger RNA (mRNA) degradation is essential for cellular metabolism and homeostasis. Controlled and specific degradation of both molecular species necessitates their engagements with the respective degradation machineries; this engagement involves a disordered/unstructured segment of the substrate traversing the degradation tunnel of the machinery and accessing the catalytic sites. However, while molecular factors influencing protein degradation have been extensively explored on a genome scale, and in multiple organisms, such a comprehensive understanding remains missing for mRNAs. Here, we analyzed multiple genome-scale experimental yeast mRNA half-life data in light of experimentally derived mRNA secondary structures and protein binding data, along with high-resolution X-ray crystallographic structures of the RNase machines. Results unraveled a consistent genome-scale trend that mRNAs comprising longer terminal and/or internal unstructured segments have significantly shorter half-lives; the lengths of the 5'-terminal, 3'-terminal, and internal unstructured segments that affect mRNA half-life are compatible with molecular structures of the 5' exo-, 3' exo-, and endoribonuclease machineries. Sequestration into ribonucleoprotein complexes elongates mRNA half-life, presumably by burying ribonuclease engagement sites under oligomeric interfaces. After gene duplication, differences in terminal unstructured lengths, proportions of internal unstructured segments, and oligomerization modes result in significantly altered half-lives of paralogous mRNAs. Side-by-side comparison of molecular principles underlying controlled protein and mRNA degradation in yeast unravels their remarkable mechanistic similarities and suggests how the intrinsic structural features of the two molecular species, at two different levels of the central dogma, regulate their half-lives on genome scale.

Why PGD2 has different functions from PGE2

Prostaglandin (PG) D₂ and PGE₂ are positional isomers; however, they sometimes exhibit opposite physiological functions, such as in cancer development. Because DP receptors are considered to be a duplicated copy of EP2 receptors, PGD₂ and PGE₂ cross-react with both receptors. These prostanoids may act as biased agonists for each receptor. In reviewing this field, a hypothesis was proposed to explain the opposed effects of these prostanoids from the viewpoints of the evolution of, mutations in, and biased activities of their receptors. Previous findings showing more mutations/variations in DP receptors than EP2 receptors among individuals worldwide indicate that DP receptors are still in a rapid evolutionary stage. The opposing effects of these prostanoids on cancer development may be attributed to the biased activity of PGE₂ for DP receptors, which may incidentally develop during the process of the old ligand, PGE₂ gaining selectivity to newly diverged DP receptors.

Molecular basis for neofunctionalization of duplicated E3 ubiquitin ligases underlying adaptation to drought tolerance in Arabidopsis thaliana

Multigene families in plants expanded from ancestral genes via gene duplication mechanisms constitute a significant fraction of the coding genome. Although most duplicated genes are lost over time, many are retained in the genome. Clusters of tandemly arrayed genes are commonly found in the plant genome where they can promote expansion of gene families. In the present study, promoter fusion to the GUS reporter gene was used to examine the promoter architecture of duplicated E3 ligase genes that are part of group C in the Arabidopsis thaliana ATL family. Acquisition of gene expression by AtATL78, possibly generated from defective AtATL81 expression, is described. AtATL78 expression was purportedly enhanced by insertion of a TATA box within the core promoter region after a short tandem duplication that occurred during evolution of Brassicaceae lineages. This gene is associated with an adaptation to drought tolerance of A. thaliana. These findings also suggest duplicated genes could serve as a reservoir of tacit genetic information, and expression of these duplicated genes is activated upon acquisition of core promoter sequences. Remarkably, drought transcriptome profiling in response to rehydration suggests that ATL78-dependent gene expression predominantly affects genes with root-specific activities.

Out of the testis, into the ovary: biased outcomes of gene duplication and deletion in Drosophila

Gene turnover is a key source of adaptive variation. Yet most evolutionary studies have focused on gene duplication, dismissing gene deletion as a mechanism that simply eradicates redundancy. Here, I use genome-scale sequence and multi-tissue expression data from Drosophila melanogaster and Drosophila pseudoobscura to simultaneously assess the evolutionary outcomes of gene duplication and deletion in Drosophila. I find that gene duplication is more frequent than gene deletion in both species, indicating that it may play a more important role in Drosophila evolution. However, examination of several genic properties reveals that genes likely possess distinct functions after duplication that diverge further before deletion, suggesting that loss of redundancy cannot explain a majority of gene deletion events in Drosophila. Moreover, in addition to providing support for the well-known "out of the testis" origin of young duplicate genes, analyses of gene expression profiles uncover a preferential bias against deletion of old ovary-expressed genes. Therefore, I propose a novel "into the ovary" hypothesis for gene deletion in Drosophila, in which gene deletion may promote adaptation by salvaging genes that contribute to the evolution of female reproductive phenotypes. Under this combined "out of the testis, into the ovary" evolutionary model, gene duplication and deletion work in concert to generate and maintain a balanced repertoire of genes that promote sex-specific adaptation in Drosophila.

Intraspecific variation in the karyotype length and genome size of fungus-farming ants (genus Mycetophylax), with remarks on procedures for the estimation of genome size in the Formicidae by flow cytometry

Ants (Formicidae) present considerable diversity in chromosome numbers, which vary from n = 1 to n = 60, although this variation is not proportional to that in genome size, for which estimates range from 0.18 pg to 0.77 pg. Intraspecific variation in the chromosome number and karyotype structure has been reported among species, although the variation among populations of the same species has received much less attention, and there are few data on genome size. Here, we studied the karyotype length and genome size of different populations of the fungus-farming ants Mycetophylax conformis (Mayr, 1884) and Mycetophylax morschi (Emery, 1888). We also provide remarks on procedure for the estimation of ant genome size by Flow Cytometry (FCM) analysis. Chromosome number and morphology did not vary among the populations of M. conformis or the cytotypes of M. morschi, but karyotype length and genome size were significantly distinct among the populations of these ants. Our results on the variation in karyotype length and genome size among M. morschi and M. conformis populations reveal considerable diversity that would be largely overlooked by more traditional descriptions of karyotypes, which were also supported by the estimates of genome size obtained using flow cytometry. Changes in the amount of DNA reflect variation in the fine structure of the chromosomes, which may represent the first steps of karyotype evolution and may occur previously to any changes in the chromosome number.

The genome of a cave plant, Primulina huaijiensis, provides insights into adaptation to limestone karst habitats

Although whole genome duplication (WGD) has been suggested to facilitate adaptive evolution and diversification, the role of specific WGD events in promoting diversification and adaptation in angiosperms remains poorly understood. Primulina, a species-rich genus with > 180 species associated with limestone karst habitat, constitutes an ideal system for studying the impact of WGD events on speciation and evolutionary adaptation. We sequenced and assembled a chromosome-level genome of the cave-dwelling species P. huaijiensis to study gene family expansion and gene retention following WGDs. We provide evidence that P. huaijiensis has undergone two WGDs since the γ triplication event shared by all eudicots. In addition to a WGD shared by almost all Lamiales (L event), we identified a lineage-specific WGD (D event) that occurred in the early Miocene around 20.6-24.2 Myr ago and that is shared by almost the entire subtribe Didymocarpinae. We found that gene retentions following the D event led to gene family proliferation (e.g. WRKYs) that probably facilitated adaptation to the high salinity and drought stress in limestone karst. Our study highlights the role of lineage-specific WGD in species diversification and adaptation of plants from special habitats.

Polyploidization events shaped the transcription factor repertoires in legumes (Fabaceae)

Transcription factors (TFs) are essential for plant growth and development. Several legumes (e.g. soybean) are rich sources of protein and oil and have great economic importance. Here we report a phylogenomic analysis of TF families in legumes and their potential association with important traits (e.g. nitrogen fixation). We used TF DNA-binding domains to systematically screen the genomes of 15 leguminous and five non-leguminous species. Transcription factor orthologous groups (OGs) were used to estimate OG sizes in ancestral nodes using a gene birth-death model, which allowed the identification of lineage-specific expansions. The OG analysis and rate of synonymous substitutions show that major TF expansions are strongly associated with whole-genome duplication (WGD) events in the legume (approximately 58 million years ago) and Glycine (approximately 13 million years ago) lineages, which account for a large fraction of the Phaseolus vulgaris and Glycine max TF repertoires. Of the 3407 G. max TFs, 1808 and 676 have homeologs within single syntenic regions in Phaseolus vulgaris and Vitis vinifera, respectively. We found a trend for TFs expanded in legumes to be preferentially transcribed in roots and nodules, supporting their recruitment early in the evolution of nodulation in the legume clade. Some families also showed count differences between G. max and the wild soybean Glycine soja, including genes located within important quantitative trait loci. Our findings strongly support the roles of two WGDs in shaping the TF repertoires in the legume and Glycine lineages, and these are probably related to important aspects of legume and soybean biology.

The genetic adaptations of fall armyworm Spodoptera frugiperda facilitated its rapid global dispersal and invasion

The fall armyworm (Spodoptera frugiperda) is a lepidopteran insect pest that causes huge economic losses. This notorious insect pest has rapidly spread over the world in the past few years. However, the mechanisms of rapid dispersal are not well understood. Here, we report a chromosome-level assembled genome of the fall armyworm, named the ZJ-version, using PacBio and Hi-C technology. The sequenced individual was a female collected from the Zhejiang province of China and had high heterozygosity. The assembled genome size of ZJ-version was 486 Mb, containing 361 contigs with an N50 of 1.13 Mb. Hi-C scaffolding further assembled the genome into 31 chromosomes and a portion of W chromosome, representing 97.4% of all contigs and resulted in a chromosome-level genome with scaffold N50 of 16.3 Mb. The sex chromosomes were identified by genome resequencing of a single male pupa and a single female pupa. About 28% of the genome was annotated as repeat sequences, and 22,623 protein-coding genes were identified. Comparative genomics revealed the expansion of the detoxification-associated gene families, chemoreception-associated gene families, nutrition metabolism and transport system gene families in the fall armyworm. Transcriptomic and phylogenetic analyses focused on these gene families revealed the potential roles of the genes in polyphagia and invasion of fall armyworm. The high-quality of the fall armyworm genome provides an important genomic resource for further explorations of the mechanisms of polyphagia and insecticide resistance, as well as for pest management of fall armyworm.

Evolutionary strategies drive a balance of the interacting gene products for the CBL and CIPK gene families

Genes encoding interacting proteins tend to be co-retained after whole-genome duplication (WGD). The preferential retention after WGD has been explained by the gene balance hypothesis (GBH). However, small-scale duplications could independently occur in the connected gene families. Certain evolutionary strategies might keep the dosage balanced. Here, we examined the gene duplication, interaction and expression patterns of calcineurin B-like (CBL) and CBL-interacting protein kinase (CIPK) gene families to understand the underlying principles. The ratio of the CBL and CIPK gene numbers evolved from 5 : 7 in Physcomitrella to 10 : 26 in Arabidopsis, and retrotransposition, tandem duplication, and WGDs contributed to the expansion. Two pairs of CBLs and six pairs of CIPKs were retained after the α WGD in Arabidopsis, in which specific interaction patterns were identified. In some cases, two retained CBLs (CIPKs) might compete to interact with a sole CIPK (CBL). Results of gene expression analyses indicated that the relatively over-retained duplicates tend to show asymmetric expression, thus avoiding competition. In conclusion, our results suggested that the highly specific interaction, together with the differential gene expression pattern, jointly maintained the balanced dosage for the interacting CBL and CIPK proteins.

Nuclear phylogenomic analyses of asterids conflict with plastome trees and support novel relationships among major lineages

PREMISE

Discordance between nuclear and organellar phylogenies (cytonuclear discordance) is a well-documented phenomenon at shallow evolutionary levels but has been poorly investigated at deep levels of plant phylogeny. Determining the extent of cytonuclear discordance across major plant lineages is essential not only for elucidating evolutionary processes, but also for evaluating the currently used framework of plant phylogeny, which is largely based on the plastid genome.

METHODS

We present a phylogenomic examination of a major angiosperm clade (Asteridae) based on sequence data from the nuclear, plastid, and mitochondrial genomes as a means of evaluating currently accepted relationships inferred from the plastome and exploring potential sources of genomic conflict in this group.

RESULTS

We recovered at least five instances of well-supported cytonuclear discordance concerning the placements of major asterid lineages (i.e., Ericales, Oncothecaceae, Aquifoliales, Cassinopsis, and Icacinaceae). We attribute this conflict to a combination of incomplete lineage sorting and hybridization, the latter supported in part by previously inferred whole-genome duplications.

CONCLUSIONS

Our results challenge several long-standing hypotheses of asterid relationships and have implications for morphological character evolution and for the importance of ancient whole-genome duplications in early asterid evolution. These findings also highlight the value of reevaluating broad-scale angiosperm and green-plant phylogeny with nuclear genomic data.

Origin and adaptation of green-sensitive (RH2) pigments in vertebrates

One of the critical times for the survival of animals is twilight where the most abundant visible lights are between 400 and 550 nanometres (nm). Green-sensitive RH2 pigments help nonmammalian vertebrate species to better discriminate wavelengths in this blue-green region. Here, evaluation of the wavelengths of maximal absorption (λmax s) of genetically engineered RH2 pigments representing 13 critical stages of vertebrate evolution revealed that the RH2 pigment of the most recent common ancestor of vertebrates had a λmax of 503 nm, while the 12 ancestral pigments exhibited an expanded range in λmax s between 474 and 524 nm, and present-day RH2 pigments have further expanded the range to ~ 450-530 nm. During vertebrate evolution, eight out of the 16 significant λmax shifts (or |Δλmax | ≥ 10 nm) of RH2 pigments identified were fully explained by the repeated mutations E122Q (twice), Q122E (thrice) and M207L (twice), and A292S (once). Our data indicated that the highly variable λmax s of teleost RH2 pigments arose from gene duplications followed by accelerated amino acid substitution.