Biological ammonium transporters: evolution and diversification.

Although ammonium is the preferred nitrogen source for microbes and plants, in animal cells it is a toxic product of nitrogen metabolism that needs to be excreted. Thus, ammonium movement across biological membranes, whether for uptake or excretion, is a fundamental and ubiquitous biological process catalysed by the superfamily of the Amt/Mep/Rh transporters. A remarkable feature of the Amt/Mep/Rh family is that they are ubiquitous and, despite sharing low amino acid sequence identity, are highly structurally conserved. Despite sharing a common structure, these proteins have become involved in a diverse range of physiological process spanning all domains of life, with reports describing their involvement in diverse biological processes being published regularly. In this context, we exhaustively present their range of biological roles across the domains of life and after explore current hypotheses concerning their evolution to help to understand how and why the conserved structure fulfils diverse physiological functions.

Insights into sequence characteristics and evolutionary history of DGATs in arthropods.

Triacylglycerol (TAG) is crucial in animal energy storage and membrane biogenesis. The conversion of diacylglycerol (DAG) to triacylglycerol (TAG) is catalyzed by diacylglycerol acyltransferase enzymes (DGATs), which are encoded by genes belonging to two distinct gene families. Although arthropods are known to possess DGATs activities and utilize the glycerol-3-phosphate pathway and MAG pathway for TAG biosynthesis, the sequence characterization and evolutionary history of DGATs in arthropods remains unclear. This study aimed to comparatively evaluate genomic analyses of DGATs in 13 arthropod species and 14 outgroup species. We found that arthropods lack SOAT2 genes within the DGAT1 family, while DGAT2, MOGAT3, AWAT1, and AWAT2 were absent from in DGAT2 family. Gene structure and phylogenetic analyses revealed that DGAT1 and DGAT2 genes come from different gene families. The expression patterns of these genes were further analyzed in crustaceans, demonstrating the importance of DGAT1 in TAG biosynthesis. Additionally, we identified the DGAT1 gene in Swimming crab (P. trituberculatus) undergoes a mutually exclusive alternative splicing event in the molt stages. Our newly determined DGAT inventory data provide a more complete scenario and insights into the evolutionary dynamics and functional diversification of DGATs in arthropods.

MpTGA, together with MpNPR, regulates sexual reproduction and independently affects oil body formation in Marchantia polymorpha.

In angiosperms, basic leucine-zipper (bZIP) TGACG-motif-binding (TGA) transcription factors (TFs) regulate developmental and stress-related processes, the latter often involving NON EXPRESSOR OF PATHOGENESIS-RELATED GENES (NPR) coregulator interactions. To gain insight into their functions in an early diverging land-plant lineage, the single MpTGA and sole MpNPR genes were investigated in the liverwort Marchantia polymorpha. We generated Marchantia MpTGA and MpNPR knockout and overexpression mutants and conducted morphological, transcriptomic and expression studies. Furthermore, we investigated MpTGA interactions with wild-type and mutagenized MpNPR and expanded our analyses including TGA TFs from two streptophyte algae. Mptga mutants fail to induce the switch from vegetative to reproductive development and lack gametangiophore formation. MpTGA and MpNPR proteins interact and Mpnpr mutant analysis reveals a novel coregulatory NPR role in sexual reproduction. Additionally, MpTGA acts independently of MpNPR as a repressor of oil body (OB) formation and can thereby affect herbivory. The single MpTGA TF exerts a dual role in sexual reproduction and OB formation in Marchantia. Common activities of MpTGA/MpNPR in sexual development suggest that coregulatory interactions were established after emergence of land-plant-specific NPR genes and contributed to the diversification of TGA TF functions during land-plant evolution.

Comparative Studies on Duplicated foxl2 Paralogs in Spotted Knifejaw Oplegnathus punctatus Show Functional Diversification.

As a member of the forkhead box L gene family, foxl2 plays a significant role in gonadal development and the regulation of reproduction. During the evolution of deuterostome, whole genome duplication (WGD)-enriched lineage diversifications and regulation mechanisms occurs. However, only limited research exists on foxl2 duplication in teleost or other vertebrate species. In this study, two foxl2 paralogs, foxl2 and foxl2l, were identified in the transcriptome of spotted knifejaw (Oplegnathus punctatus), which had varying expressions in the gonads. The foxl2 was expressed higher in the ovary, while foxl2l was expressed higher in the testis. Phylogenetic reconstruction, synteny analysis, and the molecular evolution test confirmed that foxl2 and foxl2l likely originated from the first two WGD. The expression patterns test using qRT-PCR and ISH as well as motif scan analysis revealed evidence of potentially functional divergence between the foxl2 and foxl2l paralogs in spotted knifejaw. Our results indicate that foxl2 and foxl2l may originate from the first two WGD, be active in transcription, and have undergone functional divergence. These results shed new light on the evolutionary trajectories of foxl2 and foxl2l and highlights the need for further detailed functional analysis of these two duplicated paralogs.

Bivalves Present the Largest and Most Diversified Repertoire of Toll-Like Receptors in the Animal Kingdom, Suggesting Broad-Spectrum Pathogen Recognition in Marine Waters.

Toll-like receptors (TLRs) are the most widespread class of membrane-bound innate immune receptors, responsible of specific pathogen recognition and production of immune effectors through the activation of intracellular signaling cascades. The repertoire of TLRs was analyzed in 85 metazoans, enriched on molluscan species, an underrepresented phylum in previous studies. Following an ancient evolutionary origin, suggested by the presence of TLR genes in Anthozoa (Cnidaria), these receptors underwent multiple independent gene family expansions, the most significant of which occurred in bivalve molluscs. Marine mussels (Mytilus spp.) had the largest TLR repertoire in the animal kingdom, with evidence of several lineage-specific expanded TLR subfamilies with different degrees of orthology conservation within bivalves. Phylogenetic analyses revealed that bivalve TLR repertoires were more diversified than their counterparts in deuterostomes or ecdysozoans. The complex evolutionary history of TLRs, characterized by lineage-specific expansions and losses, along with episodic positive selection acting on the extracellular recognition domains, suggests that functional diversification might be a leading evolutionary force. We analyzed a comprehensive transcriptomic data set from Mytilus galloprovincialis and built transcriptomic correlation clusters with the TLRs expressed in gills and in hemocytes. The implication of specific TLRs in different immune pathways was evidenced, as well as their specific modulation in response to different biotic and abiotic stimuli. We propose that, in a similar fashion to the remarkable functional specialization of vertebrate TLRs, the expansion of the TLR gene family in bivalves attends to a functional specification motivated by the biological particularities of these organisms and their living environment.

Phytochrome higher order mutants reveal a complex set of light responses in the moss Physcomitrium patens.

Phytochromes are photoreceptors enabling plants to respond to various light conditions. Independent gene duplications resulted in small phytochrome families in mosses, ferns and seed plants. This phytochrome diversity is hypothesised to be critical for sensing and adapting to different light conditions, but experimental evidence for this idea is lacking for mosses and ferns. The moss model species Physcomitrium patens contains seven phytochromes grouped into three clades, PHY1/3, PHY2/4 and PHY5. Here, we used CRISPR/Cas9-generated single and higher order mutants to investigate their role in light regulation of protonema and gametophore growth, protonema branching and induction of gametophores. We found both specific and partially overlapping roles for the three phytochrome clades in regulating these responses in different light conditions. PHY1/3 clade phytochromes act as primary far-red light receptors, while PHY5 clade phytochromes are the primary red light receptors. PHY2/4 clade phytochromes have functions in both red and far-red light. We also observed that PHY1/3 and PHY2/4 clade phytochromes promote gametophore growth in simulated canopy shade and also play a role in blue light. Similar to seed plants, gene duplications in the phytochrome lineage in mosses were followed by functional diversification into red and far-red light-sensing phytochromes.

Mutation-driven parallel evolution in emergence of ACE2-utilizing sarbecoviruses.

Mutation and recombination are two major genetic mechanisms that drive the evolution of viruses. They both exert an interplay during virus evolution, in which mutations provide a first ancestral source of genetic diversity for subsequent recombination. Sarbecoviruses are a group of evolutionarily related ß-coronaviruses including human severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2 and a trove of related animal viruses called SARS-like CoVs (SL-CoVs). This group of members either use or not use angiotensin-converting enzyme 2 (ACE2) as their entry receptor, which has been linked to the properties of their spike protein receptor binding domains (RBDs). This raises an outstanding question regarding how ACE2 binding originated within sarbecoviruses. Using a combination of analyses of phylogenies, ancestral sequences, structures, functions and molecular dynamics, we provide evidence in favor of an evolutionary scenario, in which three distinct ancestral RBDs independently developed the ACE2 binding trait via parallel amino acid mutations. In this process, evolutionary intermediate RBDs might be firstly formed through loop extensions to offer key functional residues accompanying point mutations to remove energetically unfavorable interactions and to change the dynamics of the functional loops, all required for ACE2 binding. Subsequent optimization in the context of evolutionary intermediates led to the independent emergence of ACE2-binding RBDs in the SARS-CoV and SARS-CoV-2 clades of Asian origin and the clade comprising SL-CoVs of European and African descent. These findings will help enhance our understanding of mutation-driven evolution of sarbecoviruses in their early history.

Flavan-3-ol Galloylation-Related Functional Gene Cluster and the Functional Diversification of SCPL Paralogs in Camellia sp.

The high accumulation of galloylated flavan-3-ols in Camellia sp. is a noteworthy phenomenon. We identified a flavan-3-ol galloylation-related functional gene cluster in tannin-rich plant Camellia sp., which included UGT84A22 and SCPL-AT gene clusters. We investigated the possible correlation between the accumulation of metabolites and the expression of SCPL-ATs and UGT84A22. The results revealed that C. sinensis, C. ptilophylla, and C. oleifera accumulated galloylated cis-flavan-3-ols (EGCG), galloylated trans-flavan-3-ols (GCG), and hydrolyzed tannins, respectively; however, C. nitidissima did not accumulate any galloylated compounds. C. nitidissima exhibited no expression of SCPL-AT or UGT84A22, whereas the other three species of Camellia exhibited various expression patterns. This indicated that the functions of the paralogs of SCPL-AT vary. Enzymatic analysis revealed that SCPL5 was neofunctionalized as a noncatalytic chaperone paralog, a type of chaerone-like protein, associating with flavan-3-ol galloylation; moreover, CsSCPL4 was subfunctionalized in association with the galloylation of cis- and trans-flavan-3-ols. In C. nitidissima, an SCPL4 homolog was noted with mutations in two cysteine residues forming a disulfide bond, which suggested that this homolog was defunctionalized. The findings of this study improve our understanding of the functional diversification of SCPL paralogs in Camellia sp.

Mutation-driven evolution of antibacterial function in an ancestral antifungal scaffold: Significance for peptide engineering.

Mutation-driven evolution of novel function on an old gene has been documented in many development- and adaptive immunity-related genes but is poorly understood in immune effector molecules. Drosomycin-type antifungal peptides (DTAFPs) are a family of defensin-type effectors found in plants and ecdysozoans. Their primitive function was to control fungal infection and then co-opted for fighting against bacterial infection in plants, insects, and nematodes. This provides a model to study the structural and evolutionary mechanisms behind such functional diversification. In the present study, we determined the solution structure of mehamycin, a DTAFP from the Northern root-knot nematode Meloidogyne hapla with antibacterial activity and an 18-mer insert, and studied the mutational effect through using a mutant with the insert deleted. Mehamycin adopts an expected cysteine-stabilized α-helix and ß-sheet fold in its core scaffold and the inserted region, called single Disulfide Bridge-linked Domain (abbreviated as sDBD), forms an extended loop protruding from the scaffold. The latter folds into an amphipathic architecture stabilized by one disulfide bridge, which likely confers mehamycin a bacterial membrane permeability. Deletion of the sDBD remarkably decreased the ability but accompanying an increase in thermostability, indicative of a structure-function trade-off in the mehamycin evolution. Allosteric analysis revealed an interior interaction between the two domains, which might promote point mutations at some key sites of the core domain and ultimately give rise to the emergence of antibacterial function. Our work may be valuable in guiding protein engineering of mehamycin to improve its activity and stability.

Functional diversification of heat shock factors.

Heat shock transcription factors (HSFs) are widely known as master regulators of the heat shock response. In invertebrates, a single heat shock factor, HSF1, is responsible for the maintenance of protein homeostasis. In vertebrates, seven members of the HSF family have been identified, namely HSF1, HSF2, HSF3, HSF4, HSF5, HSFX, and HSFY, of which HSF1 and HSF2 are clearly associated with heat shock response, while HSF4 is involved in development. Other members of the family have not yet been studied as extensively. Besides their role in cellular proteostasis, HSFs influence a plethora of biological processes such as aging, development, cell proliferation, and cell differentiation, and they are implicated in several pathologies such as neurodegeneration and cancer. This is achieved by regulating the expression of a great variety of genes including chaperones. Here, we review our current knowledge on the function of HSF family members and important aspects that made possible the functional diversification of HSFs.

Genome-Wide Classification and Evolutionary Analysis Reveal Diverged Patterns of Chalcone Isomerase in Plants.

Flavonoids as a class of important secondary metabolites are widely present in land plants, and chalcone isomerase (CHI) is the key rate-limiting enzyme that participates in catalyzing the stereospecific isomerization of chalcones to yield their corresponding flavanones. However, the phylogenetic dynamics and functional divergence of CHI family genes during the evolutionary path of green plants remains poorly understood. Here, a total of 122 CHI genes were identified by performing a genome-wide survey of 15 representative green plants from the most ancestral basal plant chlorophyte algae to higher angiosperm plants. Phylogenetic, orthologous groups (OG) classification, and genome structure analysis showed that the CHI family genes have evolved into four distinct types (types I-IV) containing eight OGs after gene duplication, and further studies indicated type III CHIs consist of three subfamilies (FAP1, FAP2, and FAP3). The phylogeny showed FAP3 CHIs as an ancestral out-group positioned on the outer layers of the main branch, followed by type IV CHIs, which are placed in an evolutionary intermediate between FAP3 CHIs and bona fide CHIs (including type I and type II). The results imply a potential intrinsic evolutionary connection between CHIs existing in the green plants. The amino acid substitutions occurring in several residues have potentially affected the functional divergence between CHI proteins. This is supported by the analysis of transcriptional divergence and cis-acting element analysis. Evolutionary dynamics analyses revealed that the differences in the total number of CHI family genes in each plant are primarily attributed to the lineage-specific expansion by natural selective forces. The current studies provide a deeper understanding of the phylogenetic relationships and functional diversification of CHI family genes in green plants, which will guide further investigation on molecular characteristics and biological functions of CHIs.

Evolution and functional diversification of catalase genes in the green lineage.

BACKGROUND: Catalases (CATs) break down hydrogen peroxide into water and oxygen to prevent cellular oxidative damage, and play key roles in the development, biotic and abiotic stresses of plants. However, the evolutionary relationships of the plant CAT gene family have not been systematically reported. RESULTS: Here, we conducted genome-wide comparative, phylogenetic, and structural analyses of CAT orthologs from 29 out of 31 representative green lineage species to characterize the evolution and functional diversity of CATs. We found that CAT genes in land plants were derived from core chlorophytes and detected a lineage-specific loss of CAT genes in Fabaceae, suggesting that the CAT genes in this group possess divergent functions. All CAT genes were split into three major groups (group α, ß1, and ß2) based on the phylogeny. CAT genes were transferred from bacteria to core chlorophytes and charophytes by lateral gene transfer, and this led to the independent evolution of two types of CAT genes: α and ß types. Ten common motifs were detected in both α and ß groups, and ß CAT genes had five unique motifs, respectively. The findings of our study are inconsistent with two previous hypotheses proposing that (i) new CAT genes are acquired through intron loss and that (ii) the Cys-343 residue is highly conserved in plants. We found that new CAT genes in most higher plants were produced through intron acquisition and that the Cys-343 residue was only present in monocots, Brassicaceae and Pp_CatX7 in P. patens, which indicates the functional specificity of the CATs in these three lineages. Finally, our finding that CAT genes show high overall sequence identity but that individual CAT genes showed developmental stage and organ-specific expression patterns suggests that CAT genes have functionally diverged independently. CONCLUSIONS: Overall, our analyses of the CAT gene family provide new insights into their evolution and functional diversification in green lineage species.

Conservation and Diversity of miR166 Family Members From Highbush Blueberry (Vaccinium corymbosum) and Their Potential Functions in Abiotic Stress.

MicroRNA166 (miR166) is highly conserved and has diverse functions across plant species. The highbush blueberry (Vaccinium corymbosum) genome is thought to harbor 10 miRNA166 loci (Vco-miR166), but the extent of their evolutionary conservation or functional diversification remains unknown. In this study, we identified six additional Vco-miR166 loci based on conserved features of the miR166 family. Phylogenetic analyses showed that mature Vco-miR166s and their precursor cluster in several clades are evolutionary conserved with diverse species. The cis-regulatory elements in the Vco-miR166 promoters indicated functions related to different phytohormones and defense responses. We also identified putative targets of vco-miR166s, which targeted the same gene families, suggesting the functional conservation and diversification of Vco-miR166 family members. Furthermore, we examined the accumulation patterns of six mature Vco-miR166s in response to abiotic stresses by stem-loop reverse RT-qPCR, which revealed their upregulation under freezing, cold, and heat stress, while they were downregulated by drought compared to control growth conditions. However, Vco-miR166 members showed different expression patterns when exposed to salt stress. These results showed that conserved Vco-miR166 family members display functional diversification but also coordinately influence plant responses to abiotic stress.

The glutathione peroxidase family: Discoveries and mechanism.

The discoveries leading to our present understanding of the glutathione peroxidases (GPxs) are recalled. The cytosolic GPx, now GPx1, was first described by Mills in 1957 and claimed to depend on selenium by Rotruck et al., in 1972. With the determination of a stoichiometry of one selenium per subunit, GPx1 was established as the first selenoenzyme of vertebrates. In the meantime, the GPxs have grown up to a huge family of enzymes that prevent free radical formation from hydroperoxides and, thus, are antioxidant enzymes, but they are also involved in regulatory processes or synthetic functions. The kinetic mechanism of the selenium-containing GPxs is unusual in neither showing a defined KM nor any substrate saturation. More recently, the reaction mechanism has been investigated by the density functional theory and nuclear magnetic resonance of model compounds mimicking the reaction cycle. The resulting concept sees a selenolate oxidized to a selenenic acid. This very fast reaction results from a concerted dual attack on the hydroperoxide bond, a nucleophilic one by the selenolate and an electrophilic one by a proton that is unstably bound in the reaction center. Postulated intermediates have been identified either in the native enzymes or in model compounds.

Functional Differentiation of BnVTE4 Gene Homologous Copies in α-Tocopherol Biosynthesis Revealed by CRISPR/Cas9 Editing.

Tocopherols are essential nutrients for human health known as vitamin E. Vitamin E deficiency can have a profound effect on human health, including the central nervous system and cardiovascular and immune protection. Multiple enzymatic steps are involved in the conversion between different forms of tocopherols. Among them, γ-tocopherol methyltransferase encoded by gene VTE4 catalyzes the conversion of γ- to α-tocopherol or δ- to ß-tocopherol isoforms. However, the gene copies and their functional contribution of VTE4 homologs in Brassica napus were not elucidated. To this end, different mutation combinations of four putative BnVTE4 homologous copies were generated by using CRISPR/Cas9 genome editing technology. Editing of those BnVTE4 homologs led to a significant change of the α-tocopherol content and the ratio between α- and γ-tocopherol compared with wide-type control. Analysis of the different combinations of BnVTE4-edited homologs revealed that the contribution of the BnVTE4 individual gene displayed obvious functional differentiation in α-tocopherol biosynthesis. Their contribution could be in order of VTE4.C02-2 (BnaC02G0331100ZS) > VTE4.A02-1 (BnaA02G0247300ZS) > VTE4.A02-2 (BnaA02G0154300ZS). Moreover, the VTE4.A02-1 and VTE4.A02-2 copies might have severe functional redundancies in α-tocopherol biosynthesis. Overall, this study systemically studied the different effects of BnVTE4 homologs, which provided a theoretical basis for breeding high α-tocopherol content oilseed rape.

Evolution, family expansion, and functional diversification of plant aldehyde dehydrogenases.

Aldehyde dehydrogenases (ALDHs) act as "aldehyde scavengers" in plants, eliminating reactive aldehydes and hence performing a crucial part in response to stress. ALDH has been specified multiple activities since its identification in the mammalian system 72 years ago. But the most widely researched role in plants is their engagement in stress tolerance. Multiple ALDH families are found in both animals and plants, and many genes are substantially conserved within these two evolutionary diverse taxa, yet both have their unique members/families. A total of twenty-four ALDH protein family has been reported across organisms, where plants contain fourteen families. Surprisingly, the number of genes in the ALDH superfamily has risen in the higher plants because of genome duplication and expansion, indicating the functional versatilely. Observed expansion in the ALDH isoforms might provide high plasticity in their actions to achieve diversified roles in the plant. The physiological importance and functional diversity of ALDHs including plant development and environmental stress adaptability, and their evolution in plants has been studied extensively. Future investigations need to focus on evaluating the individual and interconnecting function of multiple ALDH isoforms across organisms in providing plants with proper development, maturation, and adaptability against harsh environmental conditions.

Alternative splicing diversifies the transcriptome and proteome of the rice blast fungus during host infection.

Alternative splicing (AS) contributes to diversifying and regulating cellular responses to environmental conditions and developmental cues by differentially producing multiple mRNA and protein isoforms from a single gene. Previous studies on AS in pathogenic fungi focused on profiling AS isoforms under a limited number of conditions. We analysed AS profiles in the rice blast fungus Magnaporthe oryzae, a global threat to rice production, using high-quality transcriptome data representing its vegetative growth (mycelia) and multiple host infection stages. We identified 4,270 AS isoforms derived from 2,413 genes, including 499 genes presumably regulated by infection-specific AS. AS appears to increase during infection, with 32.7% of the AS isoforms being produced during infection but absent in mycelia. Analysis of the isoforms observed at each infection stage showed that 636 AS isoforms were more abundant than corresponding annotated mRNAs, especially after initial hyphal penetration into host cell. Many such dominant isoforms were predicted to encode regulatory proteins such as transcription factors and phospho-transferases. We also identified the genes encoding distinct proteins via AS and confirmed the translation of some isoforms via a proteomic analysis, suggesting potential AS-mediated neo-functionalization of some genes during infection. Comprehensive profiling of the pattern of genome-wide AS during multiple stages of rice-M. oryzae interaction established a foundational resource that will help investigate the role and regulation of AS during rice infection.

Two alpha isopropylmalate synthase isozymes with similar kinetic properties are extant in the yeast Lachancea kluyveri.

The first committed step in the leucine biosynthetic pathway is catalyzed by α-isopropylmalate synthase (α-IPMS, EC 2.3.3.13), which in the Saccaromycotina subphylum of Ascomycete yeasts is frequently encoded by duplicated genes. Following a gene duplication event, the two copies may be preserved presumably because the encoded proteins diverge in either functional properties and/or cellular localization. The genome of the petite-negative budding yeast Lachancea kluyveri includes two SAKL0E10472 (LkLEU4) and SAKL0F05170 g (LkLEU4BIS) paralogous genes, which are homologous to other yeast α-IPMS sequences. Here, we investigate whether these paralogous genes encode functional α-IPMS isozymes and whether their functions have diverged. Molecular phylogeny suggested that the LkLeu4 isozyme is located in the mitochondria and LkLeu4BIS in the cytosol. Comparison of growth rates, leucine intracellular pools and mRNA levels, indicate that the LkLeu4 isozyme is the predominant α-IPMS enzyme during growth on glucose as carbon source. Determination of the kinetic parameters indicates that the isozymes have similar affinities for the substrates and for the feedback inhibitor leucine. Thus, the diversification of the physiological roles of the genes LkLEU4 and LkLEU4BIS involves preferential transcription of the LkLEU4 gene during growth on glucose and different subcellular localization, although ligand interactions have not diverged.

The GA 20-Oxidase Encoding Gene MSD1 Controls the Main Stem Elongation in Medicago truncatula.

Plant height is an important agronomic trait that is closely related to biomass yield and crop production. Despite legumes comprise one of the largest monophyletic families that are second only to grasses in terms of economic and nutritional values, due to an ancient genome duplication event, most legume plants have complex genomes, thus the molecular mechanisms that determine plant height are less known in legumes. Here, we report the identification and characterization of MAIN STEM DWARF1 (MSD1), which is required for the plant height in the model legume Medicago truncatula. Loss of function of MSD1 leads to severely reduced main stem height but normal lateral branch elongation in M. truncatula. Histological analysis revealed that the msd1-1 main stem has shorter internodes with reduced cell size and number compared with the wild type, indicating that MSD1 affects cell elongation and cell proliferation. MSD1 encodes a putative GA 20-oxidase that is expressed at significantly higher levels in the main shoot apex than in the lateral shoot apices, suggesting that MSD1 expression is associated with its effect on the main stem elongation. UPLC-MS/MS analysis showed that GA9 and GA4, two identified products of the GA 20-oxidase, were severely reduced in msd1-1, and the dwarf phenotype of msd1-1 could be rescued by supplementation with gibberellic acid GA3, confirming that MSD1 functions as a biologically active GA 20-oxidase. Moreover, we found that disruption of either MtGA20ox7 or MtGA20ox8, homologs of MSD1, has little effects on the elongation of the main stem, while the msd1-1 mtga20ox7-1 mtga20ox8 triple mutants exhibits a severe short main shoot and lateral branches, as well as reduced leaf size, suggesting that MSD1 and its homologs MtGA20ox7 and MtGA20ox8, redundantly regulate M. truncatula shoot elongation and leaf development. Taken together, our findings demonstrate the molecular mechanism of MSD1-mediated regulation of main stem elongation in M. truncatula and provide insights into understanding the functional diversity of GA 20-oxidases in optimizing plant architecture in legumes.

Intraspecific Transcriptome Variation and Sex-Biased Expression in Anopheles arabiensis.

The magnitude and functional patterns of intraspecific transcriptional variation in the anophelines, including those of sex-biased genes underlying sex-specific traits relevant for malaria transmission, remain understudied. As a result, how changes in expression levels drive adaptation in these species is poorly understood. We sequenced the female, male, and larval transcriptomes of three populations of Anopheles arabiensis from Burkina Faso. One-third of the genes were differentially expressed between populations, often involving insecticide resistance-related genes in a sample type-specific manner, and with the females showing the largest number of differentially expressed genes. At the genomic level, the X chromosome appears depleted of differentially expressed genes compared with the autosomes, chromosomes harboring inversions do not exhibit evidence for enrichment of such genes, and genes that are top contributors to functional enrichment patterns of population differentiation tend to be clustered in the genome. Further, the magnitude of variation for the sex expression ratio across populations did not substantially differ between male- and female-biased genes, except for some populations in which male-limited expressed genes showed more variation than their female counterparts. In fact, female-biased genes exhibited a larger level of interpopulation variation than male-biased genes, both when assayed in males and females. Beyond uncovering the extensive adaptive potential of transcriptional variation in An. Arabiensis, our findings suggest that the evolutionary rate of changes in expression levels on the X chromosome exceeds that on the autosomes, while pointing to female-biased genes as the most variable component of the An. Arabiensis transcriptome.