Search

Studying the diversity of energy production pathways is important for understanding the evolutionary relationships between metabolic pathways and their biochemical precursors. The lactate/malate dehydrogenase (LDH/MDH) superfamily has been a model system for structural and functional evolution for a long time. Recently, the type-2 family of LDH/MDH (or LDH2/MDH2 oxidoreductase) has been identified. The LDH2/MDH2 oxidoreductase family is now known to have functionally more diverse enzymes than the LDH/MDH superfamily. In channel catfish, the gene encoding the LDH2/MDH2 oxidoreductase has been found (and was provisionally termed AqE). Homologs of this enzyme are predominantly present in organisms living in an aquatic environment. In this work, we studied the AqE gene distribution among non-tetrapod vertebrates. It was found that the AqE gene is present in the genomes of bony and cartilaginous fish and in the genomes of hagfishes and lampreys. In addition, it has been confirmed that in representatives of Cypriniformes, the AqE gene has been lost. AqE in representatives of Salmoniformes underwent significant deletions, which most likely led to its pseudogenization. In most orders of non-Tetrapoda vertebrates, the AqE gene remains highly conserved, suggesting that the AqE gene in aquatic vertebrates is an essential gene and undergoes rigorous selection. The AqE gene has the highest sequence similarity with the archaeal ComC gene that encodes sulfolactate dehydrogenase (SLDH). Based on the similarity of substrates, the enzyme encoded by the AqE gene is likely involved in the malate–aspartate shuttle mechanism or the biosynthesis of the energy coenzyme M equivalent. [ABSTRACT FROM AUTHOR]

Gene identification and evolution and also their impacts on agronomical traits are important variations between related species which are often because of differences in the cumulative actions of multiple gene products. Understanding genetic variability for agronomic traits is important parameter of breeding programs for broadening the gene pool of different crops. The Short-Chain Dehydrogenase/Reductase (SDR) super-family is one of the largest known protein families and comprises thousands of members found in species ranging from bacteria to humans. SDR involvement has been demonstrated in a variety of primary and secondary metabolisms. Quantitative trait locus (QTL) mapping provides a starting point for dissecting complex traits into its component alleles which may help to quantify relative impacts of alleles on the traits and locates genomic regions responsible for marker-trait association, and provide a foundation of marker-assisted selection which expedites the breeding process given the proper estimation of position and the effects of QTLs. In this manuscript we have reviewed, molecular breeding and gene evolution and domestication in crops, and also survey some important reported genes which have notable impact on agronomical traits in various crops. [ABSTRACT FROM AUTHOR]

Transposable elements (TEs) are major components of eukaryotic genomes and represent approximately 45% of the human genome. TEs can be important sources of novelty in genomes and there is increasing evidence that TEs contribute to the evolution of gene regulation in mammals. Gene duplication is an evolutionary mechanism that also provides new genetic material and opportunities to acquire new functions. To investigate how duplicated genes are maintained in genomes, here, we explored the TE environment of duplicated and singleton genes. We found that singleton genes have more short-interspersed nuclear elements and DNA transposons in their vicinity than duplicated genes, whereas long-interspersed nuclear elements and long-terminal repeat retrotransposons have accumulated more near duplicated genes. We also discovered that this result is highly associated with the degree of essentiality of the genes with an unexpected accumulation of short-interspersed nuclear elements and DNA transposons around the more-essential genes. Our results underline the importance of taking into account the TE environment of genes to better understand how duplicated genes are maintained in genomes. [ABSTRACT FROM AUTHOR]

Background: The SARS-CoV-2 pandemic is one of the greatest global medical and social challenges that have emerged in recent history. Human coronavirus strains discovered during previous SARS outbreaks have been hypothesized to pass from bats to humans using intermediate hosts, e.g. civets for SARS-CoV and camels for MERS-CoV. The discovery of an intermediate host of SARS-CoV-2 and the identification of specific mechanism of its emergence in humans are topics of primary evolutionary importance. In this study we investigate the evolutionary patterns of 11 main genes of SARS-CoV-2. Previous studies suggested that the genome of SARS-CoV-2 is highly similar to the horseshoe bat coronavirus RaTG13 for most of the genes and to some Malayan pangolin coronavirus (CoV) strains for the receptor binding (RB) domain of the spike protein. Results: We provide a detailed list of statistically significant horizontal gene transfer and recombination events (both intergenic and intragenic) inferred for each of 11 main genes of the SARS-CoV-2 genome. Our analysis reveals that two continuous regions of genes S and N of SARS-CoV-2 may result from intragenic recombination between RaTG13 and Guangdong (GD) Pangolin CoVs. Statistically significant gene transfer-recombination events between RaTG13 and GD Pangolin CoV have been identified in region [1215-1425] of gene S and region [534-727] of gene N. Moreover, some statistically significant recombination events between the ancestors of SARS-CoV-2, RaTG13, GD Pangolin CoV and bat CoV ZC45-ZXC21 coronaviruses have been identified in genes ORF1ab, S, ORF3a, ORF7a, ORF8 and N. Furthermore, topology-based clustering of gene trees inferred for 25 CoV organisms revealed a three-way evolution of coronavirus genes, with gene phylogenies of ORF1ab, S and N forming the first cluster, gene phylogenies of ORF3a, E, M, ORF6, ORF7a, ORF7b and ORF8 forming the second cluster, and phylogeny of gene ORF10 forming the third cluster. Conclusions: The results of our horizontal gene transfer and recombination analysis suggest that SARS-CoV-2 could not only be a chimera virus resulting from recombination of the bat RaTG13 and Guangdong pangolin coronaviruses but also a close relative of the bat CoV ZC45 and ZXC21 strains. They also indicate that a GD pangolin may be an intermediate host of this dangerous virus. [ABSTRACT FROM AUTHOR]

Comparative genomic analyses have enormous potential for identifying key genes central to human health phenotypes, including those that promote cancers. In particular, the successful development of novel therapeutics using model species requires phylogenetic analyses to determine molecular homology. Accordingly, we investigate the evolutionary histories of anaplastic lymphoma kinase (ALK)—which can underlie tumorigenesis in neuroblastoma, nonsmall cell lung cancer, and anaplastic large-cell lymphoma—its close relative leukocyte tyrosine kinase (LTK) and their candidate ligands. Homology of ligands identified in model organisms to those functioning in humans remains unclear. Therefore, we searched for homologs of the human genes across metazoan genomes, finding that the candidate ligands Jeb and Hen-1 were restricted to nonvertebrate species. In contrast, the ligand augmentor (AUG) was only identified in vertebrates. We found two ALK -like and four AUG -like protein-coding genes in lamprey. Of these six genes, only one ALK -like and two AUG -like genes exhibited early embryonic expression that parallels model mammal systems. Two copies of AUG are present in nearly all jawed vertebrates. Our phylogenetic analysis strongly supports the presence of previously unrecognized functional convergences of ALK and LTK between actinopterygians and sarcopterygians—despite contemporaneous, highly conserved synteny of ALK and LTK. These findings provide critical guidance regarding the propriety of fish and mammal models with regard to model organism-based investigation of these medically important genes. In sum, our results provide the phylogenetic context necessary for effective investigations of the functional roles and biology of these critically important receptors. [ABSTRACT FROM AUTHOR]

Background: hes genes are chordate homologs of Drosophila genes, hairy and enhancer of split, which encode a basic helix-loop-helix (bHLH) transcriptional repressor with a WRPW motif. Various developmental functions of hes genes, including early embryogenesis and neurogenesis, have been elucidated in vertebrates. However, their orthologous relationships remain unclear partly because of less conservation of relatively short amino acid sequences, the fact that the genome was not analyzed as it is today, and species-specific genome duplication. This results in complicated gene names in vertebrates, which are not consistent in orthologs. We previously revealed that Xenopus frogs have two clusters of hes5, named “the hes5.1 cluster” and “the hes5.3 cluster”, but the origin and the conservation have not yet been revealed. Results: Here, we elucidated the orthologous and paralogous relationships of all hes genes of human, mouse, chicken, gecko, zebrafish, medaka, coelacanth, spotted gar, elephant shark and three species of frogs, Xenopus tropicalis (X. tropicalis), X. laevis, Nanorana parkeri, by phylogenetic and synteny analyses. Any duplicated hes5 were not found in mammals, whereas hes5 clusters in teleost were conserved although not as many genes as the three frog species. In addition, hes5 cluster-like structure was found in the elephant shark genome, but not found in cyclostomata. Conclusion: These data suggest that the hes5 cluster existed in the gnathostome ancestor but became a single gene in mammals. The number of hes5 cluster genes were specifically large in frogs. [ABSTRACT FROM AUTHOR]

Flowering is a vulnerable, but crucial phase in building crop yield. Proper timing of this period is therefore decisive in obtaining optimal yields. However, genetic regulation of flowering integrates many different environmental signals and is therefore extremely complex. This complexity increases in polyploid crops which carry two or more chromosome sets, like wheat, potato or rapeseed. Here, I summarize the current state of knowledge about flowering time gene copies in rapeseed (Brassica napus), an important oil crop with a complex polyploid history and a close relationship to Arabidopsis thaliana. The current data show a high demand for more targeted studies on flowering time genes in crops rather than in models, allowing better breeding designs and a deeper understanding of evolutionary principles. Over evolutionary time, some copies of rapeseed flowering time genes changed or lost their original role, resulting in subfunctionalization of the respective homologs. For useful applications in breeding, such patterns of subfunctionalization need to be identified and better understood.

Flowering is a rigorously timed and morphologically complex shift in plant development. This change depends on endogenous as well as environmental factors. FLOWERING LOCUS T (FT) integrates several cues from different pathways acting as a flowering promoter. Contrary to the role of FT, its paralog TERMINAL FLOWER 1 (TFL1) delays floral transition. Although FT/TFL1 homologs have been studied in model eudicots and monocots, scarce studies are available in non-model monocots like the Orchidaceae. Orchids are very diverse and their floral complexity is translated into a unique aesthetic display, which appeals the ornamental plant market. Nonetheless, orchid trade faces huge limitations due to their long vegetative phase and intractable indoor flowering seasons. Little is known about the genetic basis that control reproductive transition in orchids and, consequently, manipulating their flowering time remains a challenge. In order to contribute to the understanding of the genetic bases that control flowering in orchids we present here the first broad-scale analysis of FT/TFL1-like genes in monocots with an expanded sampling in Orchidaceae. We also compare expression patterns in three selected species and propose hypotheses on the putative role of these genes in their reproductive transition. Our findings show that FT-like genes are by far more diversified than TFL1-like genes in monocots with six subclades in the former and only one in the latter. Within MonFT1, the comparative protein sequences of MonFT1A and MonFT1B suggest that they could have recruited functional roles in delaying flowering, a role typically assigned to TFL1-like proteins. On the other hand, MonFT2 proteins have retained their canonical motifs and roles in promoting flowering transition. This is also shown by their increased expression levels from the shoot apical meristem (SAM) and leaves to inflorescence meristems (IM) and floral buds (FBs). Finally, TFL1-like genes are retained as single copy and often times are lost. Their loss could be linked to the parallel recruitment of MonFT1A and MonFT1B homologs in delaying flowering and maintaining indeterminacy of the inflorescence meristem. These hypotheses lay the foundation for future functional validation in emerging model orchid species and comparative analyses in orchids with high horticultural potential in the market.

Abstract Background The SWEET (Sugars Will Eventually be Exported Transporters) gene family is a recently identified group of sugar transporters that play an indispensable role in sugar efflux, phloem loading, plant-pathogen interaction, nectar secretion, and reproductive tissue development. However, little information on Saccharum SWEET is available for this crop with a complex genetic background. Results In this study, 22 SWEET genes were identified from Saccharum spontaneum Bacterial Artificial Chromosome libraries sequences. Phylogenetic analyses of SWEETs from 11 representative plant species showed that gene expansions of the SWEET family were mainly caused by the recent gene duplication in dicot plants, while these gene expansions were attributed to the ancient whole genome duplication (WGD) in monocot plant species. Gene expression profiles were obtained from RNA-seq analysis. SWEET1a and SWEET2s had higher expression levels in the transitional zone and maturing zone than in the other analyzed zones. SWEET1b was mainly expressed in the leaf tissues and the mature zone of the leaf of both S. spontaneum and S. officinarum, and displayed a peak in the morning and was undetectable in both sclerenchyma and parenchyma cells from the mature stalks of S. officinarum. SsSWEET4a\4b had higher expression levels than SWEET4c and were mainly expressed in the stems of seedlings and mature plants. SWEET13s are recently duplicated genes, and the expression of SWEET13s dramatically increased from the maturing to mature zones. SWEET16b’s expression was not detected in S. officinarum, but displayed a rhythmic diurnal expression pattern. Conclusions Our study revealed the gene evolutionary history of SWEETs in Saccharum and SWEET1b was found to be a sucrose starvation-induced gene involved in the sugar transportation in the high photosynthetic zones. SWEET13c was identified as the key player in the efflux of sugar transportation in mature photosynthetic tissues. SWEET4a\4b were found to be mainly involved in sugar transportation in the stalk. SWEET1a\2a\4a\4b\13a\16b were suggested to be the genes contributing to the differences in sugar contents between S. spontaneum and S. officinarum. Our results are valuable for further functional analysis of SWEET genes and utilization of the SWEET genes for genetic improvement of Saccharum for biofuel production.

Abstract Background Spanins are phage lysis proteins required to disrupt the outer membrane. Phages employ either two-component spanins or unimolecular spanins in this final step of Gram-negative host lysis. Two-component spanins like Rz-Rz1 from phage lambda consist of an integral inner membrane protein: i-spanin, and an outer membrane lipoprotein: o-spanin, that form a complex spanning the periplasm. Two-component spanins exist in three different genetic architectures; embedded, overlapped and separated. In contrast, the unimolecular spanins, like gp11 from phage T1, have an N-terminal lipoylation signal sequence and a C-terminal transmembrane domain to account for the topology requirements. Our proposed model for spanin function, for both spanin types, follows a common theme of the outer membrane getting fused with the inner membrane, effecting the release of progeny virions. Results Here we present a SpaninDataBase which consists of 528 two-component spanins and 58 unimolecular spanins identified in this analysis. Primary analysis revealed significant differences in the secondary structure predictions for the periplasmic domains of the two-component and unimolecular spanin types, as well as within the three different genetic architectures of the two-component spanins. Using a threshold of 40% sequence identity over 40% sequence length, we were able to group the spanins into 143 i-spanin, 125 o-spanin and 13 u-spanin families. More than 40% of these families from each type were singletons, underlining the extreme diversity of this class of lysis proteins. Multiple sequence alignments of periplasmic domains demonstrated conserved secondary structure patterns and domain organization within family members. Furthermore, analysis of families with members from different architecture allowed us to interpret the evolutionary dynamics of spanin gene arrangement. Also, the potential universal role of intermolecular disulfide bonds in two-component spanin function was substantiated through bioinformatic and genetic approaches. Additionally, a novel lipobox motif, AWAC, was identified and experimentally verified. Conclusions The findings from this bioinformatic approach gave us instructive insights into spanin function, evolution, domain organization and provide a platform for future spanin annotation, as well as biochemical and genetic experiments. They also establish that spanins, like viral membrane fusion proteins, adopt different strategies to achieve fusion of the inner and outer membranes.

Abstract Background Hibernation allows animals to survive periods of resource scarcity by reducing their energy expenditure through decreased metabolism. However, hibernators become susceptible to psychrophilic pathogens if they cannot mount an efficient immune response to infection. While Nearctic bats infected with white-nose syndrome (WNS) suffer high mortality, related Palearctic taxa are better able to survive the disease than their Nearctic counterparts. We hypothesised that WNS exerted historical selective pressure in Palearctic bats, resulting in genomic changes that promote infection tolerance. Results We investigated partial sequences of 23 genes related to water metabolism and skin structure function in nine Palearctic and Nearctic hibernating bat species and one non-hibernating species for phylogenetic signals of natural selection. Using maximum likelihood analysis, we found that eight genes were under positive selection and we successfully identified amino acid sites under selection in five encoded proteins. Branch site models revealed positive selection in three genes. Hibernating bats exhibit signals for positive selection in genes ensuring tissue regeneration, wound healing and modulation of the immune response. Conclusion Our results highlight the importance of skin barrier integrity and healing capacity in hibernating bats. The protective role of skin integrity against both pathophysiology and WNS progression, in synergy with down-regulation of the immune reaction in response to the Pseudogymnoascus destructans infection, improves host survival. Our data also suggest that hibernating bat species have evolved into tolerant hosts by reducing the negative impact of skin infection through a set of adaptations, including those at the genomic level.

Abstract Background Most amino acids are encoded by multiple synonymous codons. However, synonymous codons are not used equally, and this biased codon use varies between different organisms. It has previously been shown that both selection acting to increase codon translational efficiency and selection acting to decrease codon biosynthetic cost contribute to differences in codon bias. However, it is unknown how these two factors interact or how they affect molecular sequence evolution. Results Through analysis of 1320 bacterial genomes, we show that bacterial genes are subject to multi-objective selection-driven optimization of codon use. Here, selection acts to simultaneously decrease transcript biosynthetic cost and increase transcript translational efficiency, with highly expressed genes under the greatest selection. This optimization is not simply a consequence of the more translationally efficient codons being less expensive to synthesize. Instead, we show that transfer RNA gene copy number alters the cost-efficiency trade-off of synonymous codons such that, for many species, selection acting on transcript biosynthetic cost and translational efficiency act in opposition. Finally, we show that genes highly optimized to reduce cost and increase efficiency show reduced rates of synonymous and non-synonymous mutation. Conclusions This analysis provides a simple mechanistic explanation for variation in evolutionary rate between genes that depends on selection-driven cost-efficiency optimization of the transcript. These findings reveal how optimization of resource allocation to messenger RNA synthesis is a critical factor that determines both the evolution and composition of genes.

Abstract Background Sall (Spalt-like) proteins are zinc-finger transcription factors involved in a number of biological processes. They have only been studied in a few model organisms, such as Drosophila melanogaster, Caenorhabditis elegans, Schmidtea mediterranea and some vertebrates. Further taxon sampling is critical to understand the evolution and diversification of this protein and its functional roles in animals. Results Using genome and transcriptome mining, we confirmed the presence of sall genes in a range of additional animal taxa, for which their presence had not yet been described. We show that sall genes are broadly conserved across the Bilateria, and likely appeared in the bilaterian stem lineage. Our analysis of the protein domains shows that the characteristic arrangement of the multiple zinc-finger domains is conserved in bilaterians and may represent the ancient arrangement of this family of transcription factors. We also show the existence of a previously unknown zinc-finger domain. In situ hybridization was used to describe the gene expression patterns in embryonic and larval stages in two species of snails: Crepidula fornicata and Lottia gigantea. In L. gigantea, sall presents maternal expression, although later on the expression is restricted to the A and B quadrants during gastrulation and larval stage. In C. fornicata, sall has no maternal expression and it is expressed mainly in the A, C and D quadrants during blastula stages and in an asymmetric fashion during the larval stage. Discussion Our results suggest that the bilaterian common ancestor had a Sall protein with at least six zinc-finger domains. The evolution of Sall proteins in bilaterians might have occurred mostly as a result of the loss of protein domains and gene duplications leading to diversification. The new evidence complements previous studies in highlighting an important role of Sall proteins in bilaterian development. Our results show maternal expression of sall in the snail L. gigantea, but not C. fornicata. The asymmetric expression shown in the ectoderm of the trochophore larva of snails is probably related to shell/mantle development. The observed sall expression in cephalic tissue in snails and some other bilaterians suggests a possible ancestral role of sall in neural development in bilaterians.

Complement C4 is a central protein by acting as pivotal molecule in the activation of the complement system. More than a decade ago, C4 gene duplication had been found in several species including fish, revealing the evolutionary origin of C4 gene. However, the evolutionary pattern and systematic function of C4 are still limited. In this study, C4 D and H types in different species groups were completely diverged. The codon usage of C4 H type in higher vertebrates were much closer to their own genome environment, in contrast to lower vertebrates, suggesting that the evolution may provide the dynamic for homogeneous codon usage between specific gene and genome. Multiple C4 sequence alignment showed that the sequences were conserved among different species. However, sequence similarity was obviously different between species C4 D and H type. Negative selection pressure was found on C4 gene evolution and it may be one of the possible reasons for the sequence broad similarity and conservation among interspecies. Proteins from C4 protein–protein interaction (PPI) network were enriched in more hematopoiesis, infections, diseases and immune–related pathways in human than zebrafish. The result suggested that the functional complexities of C4 isotypes are distinct in species from different evolutionary positions. The simulated C4 protein structures between human and grass carp shared structural similarity and the stereo structures of grass carp C4–MASP–2 protein complexes were further simulated according to a study of human. These results suggested that the interaction between C4 and MASP–2 proteins may also exist in grass carp. Our results can provide an insight for the evolutionary process of C4 and better understanding to the potential mechanism of interaction between C4 and MASP–2 in fish species. • Codon usage of C4 gene is much closer to genomic environment in higher vertebrates. • C4 genes from each species groups underwent negative selection pressures. • Functional complexities of two C4 isotypes tend to distinct in species from different evolutionary positions. • Protein structural analysis showed that C4 may have theoretical interaction with MASP–2 protein in grass carp. [ABSTRACT FROM AUTHOR]

Flowering is a vulnerable, but crucial phase in building crop yield. Proper timing of this period is therefore decisive in obtaining optimal yields. However, genetic regulation of flowering integrates many different environmental signals and is therefore extremely complex. This complexity increases in polyploid crops which carry two or more chromosome sets, like wheat, potato or rapeseed. Here, I summarize the current state of knowledge about flowering time gene copies in rapeseed (Brassica napus), an important oil crop with a complex polyploid history and a close relationship to Arabidopsis thaliana. The current data show a high demand for more targeted studies on flowering time genes in crops rather than in models, allowing better breeding designs and a deeper understanding of evolutionary principles. Over evolutionary time, some copies of rapeseed flowering time genes changed or lost their original role, resulting in subfunctionalization of the respective homologs. For useful applications in breeding, such patterns of subfunctionalization need to be identified and better understood. [ABSTRACT FROM AUTHOR]

The subfamily Pothoideae belongs to the ecologically important plant family Araceae. Here, we report the chloroplast genomes of two species of the subfamily Pothoideae: Anthurium huixtlense (size: 163,116 bp) and Pothos scandens (size: 164,719 bp). The chloroplast genome of P. scandens showed unique contraction and expansion of inverted repeats (IRs), thereby increasing the size of the large single-copy region (LSC: 102,956 bp) and decreasing the size of the small single-copy region (SSC: 6779 bp). This led to duplication of many single-copy genes due to transfer to IR regions from the small single-copy (SSC) region, whereas some duplicate genes became single copy due to transfer to large single-copy regions. The rate of evolution of protein-coding genes was affected by the contraction and expansion of IRs; we found higher mutation rates for genes that exist in single-copy regions as compared to those in IRs. We found a 2.3-fold increase of oligonucleotide repeats in P. scandens when compared with A. huixtlense, whereas amino acid frequency and codon usage revealed similarities. The ratio of transition to transversion mutations was 2.26 in P. scandens and 2.12 in A. huixtlense. Transversion mutations mostly translated in non-synonymous substitutions. The phylogenetic inference of the limited species showed the monophyly of the Araceae subfamilies. Our study provides insight into the molecular evolution of chloroplast genomes in the subfamily Pothoideae and family Araceae. [ABSTRACT FROM AUTHOR]

Recent advances in genetics present unique opportunities for enhancing our understanding of human physiology and disease predisposition through detailed analysis of gene structure, expression, and population variation via examination of data in publicly accessible genome and gene expression repositories. Yet, the vast majority of human genes remain understudied. Here, we show the scope of these genomic and genetic resources by evaluating ZMAT2, a member of a 5-gene family that through May 2020 had been the focus of only 4 peer-reviewed scientific publications. Using analysis of information extracted from public databases, we show that human ZMAT2 is a 6-exon gene and find that it exhibits minimal genetic variation in human populations and in disease states, including cancer. We further demonstrate that the gene and its encoded protein are highly conserved among nonhuman primates and define a cohort of ZMAT2 pseudogenes in the marmoset genome. Collectively, our investigations illustrate how complementary use of genomic, gene expression, and population genetic resources can lead to new insights about human and mammalian biology and evolution, and when coupled with data supporting key roles for ZMAT2 in keratinocyte differentiation and pre-RNA splicing argue that this gene is worthy of further study. [ABSTRACT FROM AUTHOR]

A small phylogenetically conserved sequence of 11,231 bp, termed FAM247, is repeated in human chromosome 22 by segmental duplications. This sequence forms part of diverse genes that span evolutionary time, the protein genes being the earliest as they are present in zebrafish and/or mice genomes, and the long noncoding RNA genes and pseudogenes the most recent as they appear to be present only in the human genome. We propose that the conserved sequence provides a nucleation site for new gene development at evolutionarily conserved chromosomal loci where the FAM247 sequences reside. The FAM247 sequence also carries information in its open reading frames that provides protein exon amino acid sequences; one exon plays an integral role in immune system regulation, specifically, the function of ubiquitin-specific protease (USP18) in the regulation of interferon. An analysis of this multifaceted sequence and the genesis of genes that contain it is presented. [ABSTRACT FROM AUTHOR]

Trichothecenes are sesquiterpene toxins produced by diverse but relatively few fungal species in at least three classes of Ascomycetes: Dothideomycetes, Eurotiomycetes, and Sordariomycetes. Approximately 200 structurally distinct trichothecene analogs have been described, but a given fungal species typically produces only a small subset of analogs. All trichothecenes share a core structure consisting of a four-ring nucleus known as 12,13-epoxytrichothec-9-ene. This structure can be substituted at various positions with hydroxyl, acyl, or keto groups to give rise to the diversity of trichothecene structures that has been described. Over the last 30 years, the genetic and biochemical pathways required for trichothecene biosynthesis in several species of the fungi Fusarium and Trichoderma have been elucidated. In addition, phylogenetic and functional analyses of trichothecene biosynthetic (TRI) genes from fungi in multiple genera have provided insights into how acquisition, loss, and changes in functions of TRI genes have given rise to the diversity of trichothecene structures. These analyses also suggest both divergence and convergence of TRI gene function during the evolutionary history of trichothecene biosynthesis. What has driven trichothecene structural diversification remains an unanswered question. However, insight into the role of trichothecenes in plant pathogenesis of Fusarium species and into plant glucosyltransferases that detoxify the toxins by glycosylating them point to a possible driver. Because the glucosyltransferases can have substrate specificity, changes in trichothecene structures produced by a fungus could allow it to evade detoxification by the plant enzymes. Thus, it is possible that advantages conferred by evading detoxification have contributed to trichothecene structural diversification. Key Points: • TRI genes have evolved by diverse processes: loss, acquisition and changes in function. • Some TRI genes have acquired the same function by convergent evolution. • Some other TRI genes have evolved divergently to have different functions. • Some TRI genes were acquired or resulted from diversification in function of other genes. • Substrate specificity of plant glucosyltransferases could drive trichothecene diversity. [ABSTRACT FROM AUTHOR]

Flowering is a rigorously timed and morphologically complex shift in plant development. This change depends on endogenous as well as environmental factors. FLOWERING LOCUS T (FT) integrates several cues from different pathways acting as a flowering promoter. Contrary to the role of FT , its paralog TERMINAL FLOWER 1 (TFL1) delays floral transition. Although FT / TFL1 homologs have been studied in model eudicots and monocots, scarce studies are available in non-model monocots like the Orchidaceae. Orchids are very diverse and their floral complexity is translated into a unique aesthetic display, which appeals the ornamental plant market. Nonetheless, orchid trade faces huge limitations due to their long vegetative phase and intractable indoor flowering seasons. Little is known about the genetic basis that control reproductive transition in orchids and, consequently, manipulating their flowering time remains a challenge. In order to contribute to the understanding of the genetic bases that control flowering in orchids we present here the first broad-scale analysis of FT / TFL1 - like genes in monocots with an expanded sampling in Orchidaceae. We also compare expression patterns in three selected species and propose hypotheses on the putative role of these genes in their reproductive transition. Our findings show that FT-like genes are by far more diversified than TFL1-like genes in monocots with six subclades in the former and only one in the latter. Within MonFT1 , the comparative protein sequences of MonFT1A and MonFT1B suggest that they could have recruited functional roles in delaying flowering, a role typically assigned to TFL1-like proteins. On the other hand, MonFT2 proteins have retained their canonical motifs and roles in promoting flowering transition. This is also shown by their increased expression levels from the shoot apical meristem (SAM) and leaves to inflorescence meristems (IM) and floral buds (FBs). Finally, TFL1-like genes are retained as single copy and often times are lost. Their loss could be linked to the parallel recruitment of MonFT1A and MonFT1B homologs in delaying flowering and maintaining indeterminacy of the inflorescence meristem. These hypotheses lay the foundation for future functional validation in emerging model orchid species and comparative analyses in orchids with high horticultural potential in the market. [ABSTRACT FROM AUTHOR]

Recently developed CRISPR-mediated base editors, which enable the generation of numerous nucleotide changes in target genomic regions, have been widely adopted for gene correction and generation of crop germplasms containing important gain-of-function genetic variations. However, to engineer target genes with unknown functional SNPs remains challenging. To address this issue, we present here a base-editing-mediated gene evolution (BEMGE) method, employing both Cas9n-based cytosine and adenine base editors as well as a single-guide RNA (sgRNA) library tiling the full-length coding region, for developing novel rice germplasms with mutations in any endogenous gene. To this end, OsALS1 was artificially evolved in rice cells using BEMGE through both Agrobacterium -mediated and particle-bombardment-mediated transformation. Four different types of amino acid substitutions in the evolved OsALS1, derived from two sites that have never been targeted by natural or human selection during rice domestication, were identified, conferring varying levels of tolerance to the herbicide bispyribac-sodium. Furthermore, the P171F substitution identified in a strong OsALS1 allele was quickly introduced into the commercial rice cultivar Nangeng 46 through precise base editing with the corresponding base editor and sgRNA. Collectively, these data indicate great potential of BEMGE in creating important genetic variants of target genes for crop improvement. Base-editing-mediated gene evolution (BEMGE) efficiently drives artificial evolution of target genes in planta , generating a large number of novel alleles in rice in a short time. The P171F substitution in OsALS1 renders rice plants resistant to the herbicide bispyribac-sodium. Important verified single-nucleotide polymorphism can be rapidly introduced into elite rice cultivars through precise base editing. BEMGE is of great value in the generation of important genetic variants of target genes for crop improvement. [ABSTRACT FROM AUTHOR]

Main conclusion: This study provides broad insight into the chloroplast genomes of the subfamily Monsteroideae. The identified polymorphic regions may be suitable for designing unique and robust molecular markers for phylogenetic inference. Monsteroideae is the third largest subfamily (comprises 369 species) and one of the early diverging lineages of the monocot plant family Araceae. The phylogeny of this important subfamily is not well resolved at the species level due to scarcity of genomic resources and suitable molecular markers. Here, we report annotated chloroplast genome sequences of four Monsteroideae species: Spathiphyllum patulinervum, Stenospermation multiovulatum, Monstera adansonii, and Rhaphidophora amplissima. The quadripartite chloroplast genomes (size range 163,335–164,751 bp) consist of a pair of inverted repeats (25,270–25,931 bp), separating a small single copy region (21,448–22,346 bp) from a large single copy region (89,714–91,841 bp). The genomes contain 114 unique genes, including four rRNA genes, 80 protein-coding genes, and 30 tRNA genes. Gene features, amino acid frequencies, codon usage, GC contents, oligonucleotide repeats, and inverted repeats dynamics exhibit similarities among the four genomes. Higher rate of synonymous substitutions was observed as compared to non-synonymous substitutions in 76 protein-coding genes. Positive selection was observed in seven protein-coding genes, including psbK, ndhK, ndhD, rbcL, accD, rps8, and ycf2. Our included species of Araceae showed the monophyly in Monsteroideae and other subfamilies. We report 30 suitable polymorphic regions. The polymorphic regions identified here might be suitable for designing unique and robust markers for inferring the phylogeny and phylogeography among closely related species within the genus Spathiphyllum and among distantly related species within the subfamily Monsteroideae. The chloroplast genomes presented here are a valuable contribution towards understanding the molecular evolutionary dynamics in the family Araceae. [ABSTRACT FROM AUTHOR]

Background: Recent advances in genetics and genomics present unique opportunities for enhancing our understanding of mammalian biology and evolution through detailed multi-species comparative analysis of gene organization and expression. Yet, of the more than 20,000 protein coding genes found in mammalian genomes, fewer than 10% have been examined in any detail. Here we elucidate the power of data available in publicly-accessible genomic and genetic resources by querying them to evaluate Zmat2, a minimally studied gene whose human ortholog has been implicated in spliceosome function and in keratinocyte differentiation. Results: We find extensive conservation in coding regions and overall structure of Zmat2 in 18 mammals representing 13 orders and spanning ~ 165 million years of evolutionary development, and in their encoded proteins. We identify a tandem duplication in the Zmat2 gene and locus in opossum, but not in other monotremes, marsupials, or other mammals, indicating that this event occurred subsequent to the divergence of these species from one another. We also define a collection of Zmat2 pseudogenes in half of the mammals studied, and suggest based on phylogenetic analysis that they each arose independently in the recent evolutionary past. Conclusions: Mammalian Zmat2 genes and ZMAT2 proteins illustrate conservation of structure and sequence, along with the development and diversification of pseudogenes in a large fraction of species. Collectively, these observations also illustrate how the focused identification and interpretation of data found in public genomic and gene expression resources can be leveraged to reveal new insights of potentially high biological significance. [ABSTRACT FROM AUTHOR]

Background: Wheat grain avenin-like proteins (ALPs) belong to a recently discovered class of wheat grain storage protein. ALPs in wheat grains not only have beneficial effects on dough quality but also display antifungal activities, which is a novel observation for wheat storage proteins. Previous studies have shown that ALPs are likely present in the albumin/globulin fractions of total protein extract from wheat flour. However, the accumulation characteristics of these ALPs in the mature wheat grain remains unknown. Results: In the present study, a total of 13 ALPs homologs were isolated and characterized in the albumin/globulin fractions of the wheat protein extract. A combination of multiple techniques including RP-HPLC, SDS-PAGE, MALDI-TOF and peptide sequencing were used for accurate separation and identification of individual ALP homolog. The C-terminal TaALP-by-4AL/7DS, TaALP-by-4AL/7AS/7DS, TaALP-bx/4AL/7AS/7DS, TaALP-ay-7DS, TaALP-ay-4AL, TaALP-ax-4AL, TaALP-ax-7AS, and TaALP-ax-7DS, were separated as individual protein bands from wheat flour for the first time. These unique ALPs peptides were mapped to the latest wheat genome assembly in the IWGSC database. The characteristic defence related proteins present in albumin and globulin fractions, such as protein disulfide-isomerase (PDI), grain softness protein (GSP), alpha-amylase inhibitors (AAIs) and endogenous alpha-amylase/subtilisin inhibitor were also found to co-segregate with these identified ALPs, avenin-3 and α-gliadins. The molecular weight range and the electrophoresis segregation properties of ALPs were characterised in comparison with the proteins containing the tryp_alpha_amyl domain (PF00234) and the gliadin domain (PF13016), which play a role in plant immunity and grain quality. We examined the phylogenetic relationships of the AAIs, GSP, avenin-3, α-gliadins and ALPs, based on the alignment of their functional domains. MALDI-TOF profiling indicated the occurrence of certain post-translations modifications (PTMs) in some ALP subunits. Conclusions: We reported for the first time the complete profiling of ALPs present in the albumin/globulin fractions of wheat grain protein extracts. We concluded that majority of the ALPs homologs are expressed in wheat grains. We found clear evidence of PTMs in several ALPs peptides. The identification of both gliadin domain (PF13016) and Tryp_alpha_amyl domain (PF00234) in the mature forms of ALPs highlighted the multiple functional properties of ALPs in grain quality and disease resistance. [ABSTRACT FROM AUTHOR]

Once believed to be unique features of neoplasia, chimeric RNAs are now being discovered in normal physiology. We speculated that some chimeric RNAs may be functional precursors of genes, and that forming chimeric RNA at the transcriptional level may be a 'trial' mechanism before the functional element is fixed into the genome. Supporting this idea, we identified a chimeric RNA, HNRNPA1L2-SUGT1 (H-S), whose sequence is highly similar to that of a 'pseudogene' MRPS31P5. Sequence analysis revealed that MRPS31P5 transcript is more similar to H-S chimeric RNA than its 'parent' gene, MRPS31. Evolutionarily, H-S precedes MRPS31P5, as it can be detected bioinformatically and experimentally in marmosets, which do not yet possess MRPS31P5 in their genome. Conversely, H-S is minimally expressed in humans, while instead, MRPS31P5 is abundantly expressed. Silencing H-S in marmoset cells resulted in similar phenotype as silencing MRPS31P5 in human cells. In addition, whole transcriptome analysis and candidate downstream target validation revealed common signalling pathways shared by the two transcripts. Interestingly, H-S failed to rescue the phenotype caused by silencing MPRS31P5 in human and rhesus cells, whereas MRPS31P5 can at least partially rescue the phenotype caused by silencing H-S in marmoset cells, suggesting that MRPS31P5 may have further evolved into a distinct entity. Thus, multiple lines of evidence support that MRPS31P5 is not truly a pseudogene of MRPS31, but a likely functional descendent of H-S chimera. Instead being a gene fusion product, H-S is a product of cis-splicing between adjacent genes, while MRPS31P5 is likely produced by genome rearrangement. [ABSTRACT FROM AUTHOR]

BackgroundThe last third of the 20th Century featured an accumulation of research findings that severely challenged the assumptions of the "Modern Synthesis" which provided the foundations for most biological research during that century. The foundations of that "Modernist" biology had thus largely crumbled by the start of the 21st Century. This in turn raises the question of foundations for biology in the 21st Century.ConclusionLike the physical sciences in the first half of the 20th Century, biology at the start of the 21st Century is achieving a substantive maturity of theory, experimental tools, and fundamental findings thanks to relatively secure foundations in genomics. Genomics has also forced biologists to connect evolutionary and molecular biology, because these formerly Balkanized disciplines have been brought together as actors on the genomic stage. Biologists are now addressing the evolution of genetic systems using more than the concepts of population biology alone, and the problems of cell biology using more than the tools of biochemistry and molecular biology alone. It is becoming increasingly clear that solutions to such basic problems as aging, sex, development, and genome size potentially involve elements of biological science at every level of organization, from molecule to population. The new biology knits together genomics, bioinformatics, evolutionary genetics, and other such general-purpose tools to supply novel explanations for the paradoxes that undermined Modernist biology.Open Peer ReviewersThis article was reviewed by W.F. Doolittle, E.V. Koonin, and J.M. Logsdon. For the full reviews, please go to the Reviewers' Comments section.

Abstract Background Protein disulfide isomerase (PDI) and PDI-like proteins contain thioredoxin domains that catalyze protein disulfide bond, inhibit aggregation of misfolded proteins, and function in isomerization during protein folding in endoplasmic reticulum and responses during abiotic stresses.Chinese cabbage is widely recognized as an economically important, nutritious vegetable, but its yield is severely hampered by various biotic and abiotic stresses. Because of, it is prime need to identify those genes whose are responsible for biotic and abiotic stress tolerance. PDI family genes are among of them. Results We have identified 32 PDI genes from the Br135K microarray dataset, NCBI and BRAD database, and in silico characterized their sequences. Expression profiling of those genes was performed using cDNA of plant samples imposed to abiotic stresses; cold, salt, drought and ABA (Abscisic Acid) and biotic stress; Fusarium oxysporum f. sp. conglutinans infection. The Chinese cabbage PDI genes were clustered in eleven groups in phylogeny. Among them, 15 PDI genes were ubiquitously expressed in various organs, while 24 PDI genes were up-regulated under salt and drought stress. By contrast, cold and ABA stress responsive gene number were ten and nine, respectively. In case of F. oxysporum f. sp. conglutinans infection 14 BrPDI genes were highly up-regulated. Interestingly, BrPDI1–1 gene was identified as putative candidate against abiotic (salt and drought) and biotic stresses, BrPDI5–2 gene for ABA stress, and BrPDI1–4, 6–1 and 9–2 were putative candidate genes for both cold and chilling injury stresses. Conclusions Our findings help to elucidate the involvement of PDI genes in stress responses, and they lay the foundation for functional genomics in future studies and molecular breeding of Brassica rapa crops. The stress-responsive PDI genes could be potential resources for molecular breeding of Brassica crops resistant to biotic and abiotic stresses.

Abstract Background Identifying orthologous genes is an initial step required for phylogenetics, and it is also a common strategy employed in functional genetics to find candidates for functionally equivalent genes across multiple species. At the same time, in silico orthology prediction tools often require large computational resources only available on computing clusters. Here we present OrthoReD, an open-source orthology prediction tool with accuracy comparable to published tools that requires only a desktop computer. The low computational resource requirement of OrthoReD is achieved by repeating orthology searches on one gene of interest at a time, thereby generating a reduced dataset to limit the scope of orthology search for each gene of interest. Results The output of OrthoReD was highly similar to the outputs of two other published orthology prediction tools, OrthologID and/or OrthoDB, for the three dataset tested, which represented three phyla with different ranges of species diversity and different number of genomes included. Median CPU time for ortholog prediction per gene by OrthoReD executed on a desktop computer was

The de novo emergence of new genes has been well documented through genomic analyses. However, a functional analysis, especially of very young protein-coding genes, is still largely lacking. Here, we identify a set of house mouse-specific protein-coding genes and assess their translation by ribosome profiling and mass spectrometry data. We functionally analyze one of them, Gm13030, which is specifically expressed in females in the oviduct. The interruption of the reading frame affects the transcriptional network in the oviducts at a specific stage of the estrous cycle. This includes the upregulation of Dcpp genes, which are known to stimulate the growth of preimplantation embryos. As a consequence, knockout females have their second litters after shorter times and have a higher infanticide rate. Given that Gm13030 shows no signs of positive selection, our findings support the hypothesis that a de novo evolved gene can directly adopt a function without much sequence adaptation.

ABSTRACT The origin of novel genes and beneficial functions is of fundamental interest in evolutionary biology. New genes can originate from different mechanisms, including horizontal gene transfer, duplication-divergence, and de novo from noncoding DNA sequences. Comparative genomics has generated strong evidence for de novo emergence of genes in various organisms, but experimental demonstration of this process has been limited to localized randomization in preexisting structural scaffolds. This bypasses the basic requirement of de novo gene emergence, i.e., lack of an ancestral gene. We constructed highly diverse plasmid libraries encoding randomly generated open reading frames and expressed them in Escherichia coli to identify short peptides that could confer a beneficial and selectable phenotype in vivo (in a living cell). Selections on antibiotic-containing agar plates resulted in the identification of three peptides that increased aminoglycoside resistance up to 48-fold. Combining genetic and functional analyses, we show that the peptides are highly hydrophobic, and by inserting into the membrane, they reduce membrane potential, decrease aminoglycoside uptake, and thereby confer high-level resistance. This study demonstrates that randomized DNA sequences can encode peptides that confer selective benefits and illustrates how expression of random sequences could spark the origination of new genes. In addition, our results also show that this question can be addressed experimentally by expression of highly diverse sequence libraries and subsequent selection for specific functions, such as resistance to toxic compounds, the ability to rescue auxotrophic/temperature-sensitive mutants, and growth on normally nonused carbon sources, allowing the exploration of many different phenotypes. IMPORTANCE De novo gene origination from nonfunctional DNA sequences was long assumed to be implausible. However, recent studies have shown that large fractions of genomic noncoding DNA are transcribed and translated, potentially generating new genes. Experimental validation of this process so far has been limited to comparative genomics, in vitro selections, or partial randomizations. Here, we describe selection of novel peptides in vivo using fully random synthetic expression libraries. The peptides confer aminoglycoside resistance by inserting into the bacterial membrane and thereby partly reducing membrane potential and decreasing drug uptake. Our results show that beneficial peptides can be selected from random sequence pools in vivo and support the idea that expression of noncoding sequences could spark the origination of new genes.

Background Anthocyanin compounds playing multiple biological functions can be synthesized in different parts of barley (Hordeum vulgare L.) plant. The diversity of anthocyanin molecules is related with branching the pathway to alternative ways in which dihydroflavonols may be modified either with the help of flavonoid 3′-hydroxylase (F3′H) or flavonoid 3′,5′-hydroxylase (F3′5′H)—the cytochrome P450-dependent monooxygenases. The F3′H and F3′5′H gene families are among the least studied anthocyanin biosynthesis structural genes in barley. The aim of this study was to identify and characterise duplicated copies of the F3′H and F3′5′H genes in the barley genome. Results Four copies of the F3′5′H gene (on chromosomes 4HL, 6HL, 6HS and 7HS) and two copies of the F3′H gene (on chromosomes 1HL and 6HS) were identified in barley genome. These copies have either one or two introns. Amino acid sequences analysis demonstrated the presence of the flavonoid hydroxylase-featured conserved motifs in all copies of the F3′H and F3′5′H genes with the exception of F3′5′H-3 carrying a loss-of-function mutation in a conservative cytochrome P450 domain. It was shown that the divergence between F3′H and F3′5′H genes occurred 129 million years ago (MYA) before the emergence of monocot and dicot plant species. The F3′H copy approximately occurred 80 MYA; the appearance of F3′5′H copies occurred 8, 36 and 91 MYA. qRT-PCR analysis revealed the tissue-specific activity for some copies of the studied genes. The F3′H-1 gene was transcribed in aleurone layer, lemma and pericarp (with an increased level in the coloured pericarp), whereas the F3′H-2 gene was expressed in stems only. The F3′5′H-1 gene was expressed only in the aleurone layer, and in a coloured aleurone its expression was 30-fold higher. The transcriptional activity of F3′5′H-2 was detected in different tissues with significantly higher level in uncoloured genotype in contrast to coloured ones. The F3′5′H-3 gene expressed neither in stems nor in aleurone layer, lemma and pericarp. The F3′5′H-4 gene copy was weakly expressed in all tissues analysed. Conclusion F3′H and F3′5′H-coding genes involved in anthocyanin synthesis in H. vulgare were identified and characterised, from which the copies designated F3′H-1, F3′H-2, F3′5′H-1 and F3′5′H-2 demonstrated tissue-specific expression patterns. Information on these modulators of the anthocyanin biosynthesis pathway can be used in future for manipulation with synthesis of diverse anthocyanin compounds in different parts of barley plant. Finding both the copies with tissue-specific expression and a copy undergoing pseudogenization demonstrated rapid evolutionary events tightly related with functional specialization of the duplicated members of the cytochrome P450-dependent monooxygenases gene families.

For two decades, Solanum bulbocastanum Dun. has been challenging potato geneticists with its durable resistance to Phytophthora infestans Mont. de Bary, the causal agent of late blight, the most devastating potato disease. In this study, highly specific SCAR (sequence characterized amplified region) markers for three S. bulbocastanum genes of broad-spectrum resistance to P. infestans,RB/Rpi-blb1,Rpi-blb2 and Rpi-blb3, were employed to screen a vast collection of wild tuber-bearing species from several series of Solanum L. section Petota Dumort. and also potato hybrids and varieties. In S. bulbocastanum and S. stoloniferum Schlechtd. et Bché, the sequences of the markers Rpi-blb1-820 and Rpi-sto1-890, which recognize the LRR and CC regions respectively differed inconsiderably in the orthologous genes RB/Rpi-blb1 and Rpi-sto1. These two markers were tightly linked in diverse wild Solanum genotypes. These markers as well as RB-226, another marker recognizing the LRR region of the Rpi-blb1 gene, were obviously related to high resistance to P. infestans. One more marker recognizing the CC region of the gene, RB-629, was found in much wider range of Solanum genotypes; it was not related to resistance to P. infestans and is therefore considered non-specific. In addition to Rpi-blb1 = Rpi-sto1, the Rpi-blb2 and Rpi-blb3 genes were registered in both S. bulbocastanum and S. stoloniferum accessions. In most cases, three bulbocastanum genes were not linked together. The markers specific for Rpi-blb1 = Rpi-sto1 genes were found in the Solanum series Bulbocastana (Rydb.) Hawkes (genome designation BB) and Longipedicellata Buk. (genome designation AABB) and were absent from five Pinnatisecta Rydb. species also comprising genome B. The marker for Rpi-blb2 was found in the series Bulbocastana and Longipedicellata and also in three accessions of Pinnatisecta, while the Rpi-blb3 gene, the orthologue of widely spread R2 gene for resistance to P. infestans, was present in all three Solanum series. In contrast with Rpi-blb1 = Rpi-sto1, two other functional RB homologues, Rpi-bt1 in S. bulbocastanum and RBver in S. verrucosum Schlechtd. (genome AA), apparently diverged earlier in the evolution of tuber-bearing Solanum species. Our evidence on functional Rpi-blb1 = Rpi-sto1 genes in three Petota series substantiates the assumption that the B genome of tetraploid Longipedicellata evolved from diploid Bulbocastana rather than from diploid Pinnatisecta progenitors. Similar to wild Solanum species, the markers Rpi-blb1-820 and Rpi-sto1-890 were tightly linked in advanced potato hybrids and established varieties bred from S. stoloniferum as the conventional source of bulbocastanum genes for late blight resistance, but did not always match the markers for the Rpi-blb2 and Rpi-blb3 genes and the marker for the Rysto gene as an independent indicator of stoloniferum germplasm. [ABSTRACT FROM AUTHOR]

Phylogenomics commonly aims to construct evolutionary trees from genomic sequence information. One way to approach this problem is to first estimate event-labeled gene trees (i.e., rooted trees whose non-leaf vertices are labeled by speciation or gene duplication events), and to then look for a species tree which can be reconciled with this tree through a reconciliation map between the trees. In practice, however, it can happen that there is no such map from a given event-labeled tree to any species tree. An important situation where this might arise is where the species evolution is better represented by a network instead of a tree. In this paper, we therefore consider the problem of reconciling event-labeled trees with species networks. In particular, we prove that any event-labeled gene tree can be reconciled with some network and that, under certain mild assumptions on the gene tree, the network can even be assumed to be multi-arc free. To prove this result, we show that we can always reconcile the gene tree with some multi-labeled (MUL-)tree, which can then be "folded up" to produce the desired reconciliation and network. In addition, we study the interplay between reconciliation maps from event-labeled gene trees to MUL-trees and networks. Our results could be useful for understanding how genomes have evolved after undergoing complex evolutionary events such as polyploidy. [ABSTRACT FROM AUTHOR]

Insulin-like growth factor-binding protein-2 (IGFBP-2) plays a key role in regulating growth and development by its affinity with insulin-like growth factors (IGFs). In this study, we cloned the coding sequence (CDS) of IGFBP-2a from the black porgy (Acanthopagrus schlegelii) muscle and identified that the full-length CDS of IGFBP-2a was 882 bp. Real-time quantitative PCR revealed that IGFBP-2a was most abundant in the liver of the black porgy and backcross breed (F1♀×black porgy♂) but remained lower in each tested tissue in self-cross breed (F1♀×F1♂). In addition, the IGFBP-2a expression in the liver of three breeds showed a negative correlation with their growth rates, indicating that the IGFBP-2a played a growth-inhibiting role in the three breeds. We further identified 810 bp IGFBP-2b gene from the draft genome of black porgy. Finally, we examined the IGFBP-2a and IGFBP-2b genes by scanning the genomes of the species of Perciformes and found the IGFBP-2 gene duplication took place earlier than the divergence of perciform species. Interestingly, six positively selected sites were detected in both Perciformes IGFBP-2 genes, although both genes were identified to be under purifying selection. Specially, these positively selected sites were located in the functional domains, suggesting these sites played key roles in the growth of Perciformes. Our study partially explains the molecular basis for the prepotency in black porgy hybrids, which will provide guidance for their cultivation in the future. [ABSTRACT FROM AUTHOR]

Starch is one of two crucial products of photosynthetic carbon-assimilation and mainly functions as the unit of energy storage in most crops such as rice, maize and sorghum, whereas interestingly in sugarcane that unit of energy storage is sucrose. Mature sugarcane stalk tissue has a very large apoplastic volume and contains nearly 700 mM sucrose—which is among the highest recorded sucrose concentrations in plant tissue. We identified 9 genes of starch synthases (SSs) related to the starch synthesis pathway in the genome of S. spontaneum. Based on gene structure and phylogenetic analysis, SSs genes were clustered into five clades and were relatively conserved. In S. spontaneum, the SS is a very ancient gene family, in which, SSIIIa and SSIIIb originated from the ρ-whole genome duplications (WGDs), SSIIb and SsIIc originated from gene duplication after the split of monocots and dicots; GBSSI and GBSSII in Clade V and SSIIa in Clade II were retained from the ε-WGD, and the remaining two SSs (SSI and SSIV) were retained from the very ancient gene duplication event about 355–389 million year ago (Mya). In addition, we found all SS genes were under the influence of strong purification with a Ka/Ks ratio of less than 0.5 in S. spontaneum. In the 5 families, SSIIIa, SSIIb and GBSSII had relatively predominant expression levels in all the examined tissues from the two Saccharum species, indicating the three genes were the fundamental members in the non-storage tissues, leaf or stem, which is in agreement with previous studies. Interestingly, the expression levels of SSs in stems showed significantly higher values in S. spontenum than in S. officinarum at pre-mature and mature stages. These results were negatively correlated with the sucrose levels between the two Saccharum species. At the pre-mature and mature stages, the sucrose contents in stems from S. officinarum were much higher than in stems from S. spontenum, suggesting that SSs involved in the differential of carbohydrate metabolism between the two Saccharum species. Besides, the expression of SSs displayed a clearly consistent trend in line with normal distribution under the diurnal rhythms of S. spontaneum. Moreover, the expression pattern of SSIIIa, SSIIb and GBSSII displayed a clearly consistent trend in both Saccharum species and in maize, rice, which was in accordance with photosynthetic intensity across leaf gradients. This result suggested the functional constraints for the SSs gene family in Gramineae. Our results are valuable for further functional analysis of SSs genes and provided the foundation for carbohydrate metabolism in sugarcane. [ABSTRACT FROM AUTHOR]

Horizontal gene transfer (HGT) is an important factor for the evolution of prokaryotes as well as eukaryotes. According to Walter M. Fitch, two genes are xenologs if they are separated by at least one HGT. This concept is formalized through Fitch relations, which are defined as binary relations that comprise all pairs (x, y) of genes x and y for which y has been horizontally transferred at least once since it diverged from the last common ancestor of x and y. This definition, in particular, preserves the directional character of the transfer. Fitch relations are characterized by a small set of forbidden induced subgraphs on three vertices and can be recognized in linear time. The mathematical characterization of Fitch relations is crucial to understand whether putative xenology relations are at least to some extent "biologically feasible". In this contribution, we provide two novel characterizations of Fitch relations. In particular, these results allow us directly to reconstruct gene trees (together with the location of the horizontal transfer events) that explain the underlying Fitch relation. As a biological side result, we can conclude that the phylogenetic signal to infer these gene trees is entirely contained in those pairs of genes x and y for which no directional transfer has been taken place in the common history of y and the last common ancestor of x and y. In other words, non-HGT events provide the essential information about the gene trees. In addition, we utilize the new characterizations to present an alternative, short and elegant proof of the characterization theorem established by Geiß et al. (J Math Bio 77(5), 2018). [ABSTRACT FROM AUTHOR]

MIRNA (MIR) gene origin and early evolutionary processes, such as hairpin precursor sequence origination, promoter activity acquirement and the sequence of these two processes, are fundamental and fascinating subjects. Three models, including inverted gene duplication, spontaneous evolution and transposon transposition, have been proposed for de novo origination of hairpin precursor sequence. However, these models still open to discussion. In addition, de novo origination of MIR gene promoters has not been well investigated. Here, I systematically investigated the origin of evolutionarily young polyphenol oxidase gene (PPO)-targeting MIRs, including MIR1444, MIR058 and MIR12112, and a genomic region termed AasPPO-as-hp, which contained a hairpin-forming sequence. I found that MIR058 precursors and the hairpin-forming sequence of AasPPO-as-hp originated in an ancient PPO gene through forming short inverted repeats. Palindromic-like sequences and imperfect inverted repeats in the ancient PPO gene contributed to initiate the generation of short inverted repeats probably by causing errors during DNA duplication. Analysis of MIR058 and AasPPO-as-hp promoters showed that they originated in the 3ʹ-flanking region of the ancient PPO gene. Promoter activities were gained by insertion of a CAAT-box and multiple-copper-response element (CuRE)-containing miniature inverted-repeat transposable element (MITE) in the upstream of AT-rich TATA-box-like sequence. Gain of promoter activities occurred before hairpin-forming sequence origination. Sequence comparison of MIR1444, MIR058 and MIR12112 promoters showed frequent birth and death of CuREs, indicating copper could be vital for the origination and evolution of PPO-targeting MIRs. Based on the evidence obtained, a novel model for plant MIR origination and evolution is proposed. [ABSTRACT FROM AUTHOR]

Abstract Insulin-like growth factor 1 (IGF1) actions are essential for somatic growth and tissue repair. IGF1 gene regulation is controlled by many inputs, with growth hormone playing a major role. In most mammals, the 6-exon IGF1 / Igf1 gene produces multiple transcripts via independent activity of its promoters plus alternative RNA splicing and differential polyadenylation. Here, by analyzing public genomic and RNA-sequencing repositories, I have characterized three Australian marsupial IGF1 genes. Koala, Tasmanian devil, and wallaby IGF1 are more complicated than other mammals, as they contain up to 11 exons, and encode multiple mRNAs and predicted protein precursors, including potentially novel isoforms. Moreover, just two of multiple growth hormone-stimulated transcriptional enhancers found in other IGF1 / Igf1 loci are detected in these species. These observations define Australian marsupial IGF1 genes and demonstrate that comprehensive interrogation of genomic and RNA-sequencing resources is an effective strategy for characterizing genes and gene expression in otherwise experimentally intractable organisms. Highlights • Australian marsupial IGF1 genes are more complicated than IGF1 / Igf1 genes from other mammals. • Australian marsupial IGF1 genes encode unique exons that predict novel transcripts and IGF1 protein precursors. • Australian marsupial IGF1 genes retain the two-promoter structure of other mammalian IGF1 / Igf1 genes. • Australian marsupial IGF1 loci contain only two of the GH-activated transcriptional enhancers found in other mammals. • Public genomic resources contain the raw data for understanding genes and their evolution in different organisms. [ABSTRACT FROM AUTHOR]

Many G protein‐coupled receptors have splice variants, with potentially different pharmaceutical properties, expression patterns and roles. The human brain expresses three functional splice variants of the type 2 corticotropin‐releasing hormone: CRHR2α, –β and –γ. CRHR2γ has only been reported in humans, but its phylogenetic distribution, and how and when during mammalian evolution it arose, is unknown. Based on genomic sequence analyses, we predict that a functional CRHR2γ is present in all Old World monkeys and apes, and is unique to these species. CRHR2γ arose by exaptation of an intronic sequence—already present in the common ancestor of primates and rodents—after retrotransposition of a short interspersed nuclear element (SINE) and mutations that created a 5′ donor splice site and in‐frame start codon, 32–43 million years ago. The SINE is not part of the coding sequence, only of the 5′ untranslated region and may therefore play a role in translational regulation. Putative regulatory elements and an alternative transcriptional start site were added earlier to this genomic locus by a DNA transposon. The evolutionary history of CRHR2γ confirms some of the earlier reported principles behind the "birth" of alternative exons. The functional significance of CRHR2γ, particularly in the brain, remains to be showed. [ABSTRACT FROM AUTHOR]

Epigenetic gene control is maintained by chromatin‐associated Polycomb group (PcG) and Trithorax group (TrxG) genes, which act antagonistically via the interplay between PcG and TrxG regulation to generate silenced or permissive transcriptional states. In this study, we searched for PcG/TrxG genes in 180 arthropod genomes, covering all the sequenced arthropod genomes at the time of conducting this study, to perform a global investigation of PcG/TrxG genes in a phylogenetic frame. Results of ancestral state reconstruction analysis revealed that the ancestor of arthropod species has an almost complete repertoire of PcG/TrxG genes, and most of these genes were seldom lost above order level. The domain diversity analysis indicated that the PcG/TrxG genes show variable extent of domain structure changes; some of these changes could be associated with lineage‐specific events. The likelihood ratio tests for selection pressure detected a number of PcG/TrxG genes which underwent episodic positive selection on the branch leading to the insects with holometabolous development. These results suggest that, despite their high conservation across arthropod species, different members of PcG/TrxG genes showed considerable differences in domain structure and sequence divergence in arthropod evolution. Our cross species comparisons using large‐scale genomic data provide insights into divergent evolutionary pattern on highly conserved genes in arthropods. [ABSTRACT FROM AUTHOR]

Background: Sugarcane served as the model plant for discovery of the C4 photosynthetic pathway. Magnesium is the central atom of chlorophyll, and thus is considered as a critical nutrient for plant development and photosynthesis. In plants, the magnesium transporter (MGT) family is composed of a number of membrane proteins, which play crucial roles in maintaining Mg homeostasis. However, to date there is no information available on the genomics of MGTs in sugarcane due to the complexity of the Saccharum genome. Results: Here, we identified 10 MGTs from the Saccharum spontaneum genome. Phylogenetic analysis of MGTs suggested that the MGTs contained at least 5 last common ancestors before the origin of angiosperms. Gene structure analysis suggested that MGTs family of dicotyledon may be accompanied by intron loss and pseudoexon phenomena during evolution. The pairwise synonymous substitution rates corresponding to a divergence time ranged from 142.3 to 236.6 Mya, demonstrating that the MGTs are an ancient gene family in plants. Both the phylogeny and Ks analyses indicated that SsMGT1/SsMGT2 originated from the recent ρWGD, and SsMGT7/SsMGT8 originated from the recent σ WGD. These 4 recently duplicated genes were shown low expression levels and assumed to be functionally redundant. MGT6, MGT9 and MGT10 weredominant genes in the MGT family and werepredicted to be located inthe chloroplast. Of the 3 dominant MGTs, SsMGT6 expression levels were found to be induced in the light period, while SsMGT9 and SsMTG10 displayed high expression levels in the dark period. These results suggested that SsMGT6 may have a function complementary to SsMGT9 and SsMTG10 that follows thecircadian clock for MGT in the leaf tissues of S. spontaneum. MGT3, MGT7 and MGT10 had higher expression levels Insaccharum officinarum than in S. spontaneum, suggesting their functional divergence after the split of S. spontaneum and S. officinarum. Conclusions: This study of gene evolution and expression of MGTs in S. spontaneum provided basis for the comprehensive genomic study of the entire MGT genes family in Saccharum. The results are valuable for further functional analyses of MGT genes and utilization of the MGTs for Saccharum genetic improvement. [ABSTRACT FROM AUTHOR]

Particularly PIN1, PIN protein-mediated rate-limiting auxin distribution plays a critical role in plant differentiation. Although well-characterized in Arabidopsis, little is known about the structural and functional relationship of the PIN1 gene among other plants. Here, we report that the gene structure remained conserved among bryophytes and angiosperms while the gene size varied by ~ 17%. Although the positions were conserved, highly variable intron phase suggests preference for specific regions in the gene sequence for independent events of intron insertion. Significant variation was observed across gene length for insertions and deletions that were mainly localized to the exonic regions flanking intron 1, possibly demarcating the sequences prone to deletions/duplications. The N and C-terminals showed a higher protein sequence similarity (~ 80%) compared to the central hydrophilic loop (~ 26%). In addition to the signature domains and motifs, we identified four novel uncharacterized motifs in the central divergent loop of PIN1 protein. Three different homo-loci, one each on chromosome groups 4, 6, and 7, were identified in wheat each showing dramatically different expression patterns during different plant developmental stages. Virus-induced gene silencing of the TaPIN1 gene resulted up to 26% reduction in plant height. Because of its direct role in controlling plant height along with a higher expression during stem elongation, the TaPIN1 gene can be manipulated to regulate plant height. [ABSTRACT FROM AUTHOR]

Background. Anthocyanin compounds playing multiple biological functions can be synthesized in different parts of barley (Hordeum vulgare L.) plant. The diversity of anthocyanin molecules is related with branching the pathway to alternative ways in which dihydroflavonols may be modified either with the help of flavonoid 30- hydroxylase (F30H) or flavonoid 30,50-hydroxylase (F3050H)-the cytochrome P450- dependent monooxygenases. The F30H and F3050H gene families are among the least studied anthocyanin biosynthesis structural genes in barley. The aim of this study was to identify and characterise duplicated copies of the F30H and F3050H genes in the barley genome. Results. Four copies of the F3050H gene (on chromosomes 4HL, 6HL, 6HS and 7HS) and two copies of the F30H gene (on chromosomes 1HL and 6HS) were identified in barley genome. These copies have either one or two introns. Amino acid sequences analysis demonstrated the presence of the flavonoid hydroxylase-featured conserved motifs in all copies of the F30H and F3050H genes with the exception of F3050H-3 carrying a loss-of-function mutation in a conservative cytochrome P450 domain. It was shown that the divergence between F30H and F3050H genes occurred 129 million years ago (MYA) before the emergence of monocot and dicot plant species. The F30H copy approximately occurred 80 MYA; the appearance of F3050H copies occurred 8, 36 and 91 MYA. qRT-PCR analysis revealed the tissue-specific activity for some copies of the studied genes. The F30H-1 gene was transcribed in aleurone layer, lemma and pericarp (with an increased level in the coloured pericarp), whereas the F30H-2 gene was expressed in stems only. The F3050H-1 gene was expressed only in the aleurone layer, and in a coloured aleurone its expression was 30-fold higher. The transcriptional activity of F3050H-2 was detected in different tissues with significantly higher level in uncoloured genotype in contrast to coloured ones. The F3050H-3 gene expressed neither in stems nor in aleurone layer, lemma and pericarp. The F3050H-4 gene copy was weakly expressed in all tissues analysed. Conclusion. F30H and F3050H-coding genes involved in anthocyanin synthesis in H. vulgare were identified and characterised, from which the copies designated F30H- 1, F30H-2, F3050H-1 and F3050H-2 demonstrated tissue-specific expression patterns. Information on these modulators of the anthocyanin biosynthesis pathway can be used in future for manipulation with synthesis of diverse anthocyanin compounds in different parts of barley plant. Finding both the copies with tissue-specific expression and a copy undergoing pseudogenization demonstrated rapid evolutionary events tightly related with functional specialization of the duplicated members of the cytochrome P450-dependent monooxygenases gene families. [ABSTRACT FROM AUTHOR]