Weak Effect of Gypsy Retrotransposon Bursts on Sonneratia alba Salt Stress Gene Expression

Transposable elements (TEs) are an important source of genetic diversity and can be co-opted for the regulation of host genes. However, to what extent the pervasive TE colonization of plant genomes has contributed to stress adaptation remains controversial. Plants inhabiting harsh environments in nature provide a unique opportunity to answer this question. We compared TE compositions and their evolutionary dynamics in the genomes of two mangrove species: the pioneer Sonneratia alba and its less salt-tolerant relative S. caseolaris. Age distribution, strength of purifying selection and the removal rate of LTR (long terminal repeat) retrotransposons were estimated. Phylogenetic analysis of LTR retrotransposons and their distribution in the genome of S. alba were surveyed. Small RNA sequencing and whole-genome bisulfite sequencing was conducted using leaves of S. alba. Expression pattern of LTR retrotransposons and their nearby genes were examined using RNA-seq data of S. alba under different salt treatments. S. alba possesses more TEs than S. caseolaris. Particularly, many more young Gypsy LTR retrotransposons have accumulated in S. alba than in S. caseolaris despite an increase in purifying selection against TE insertions. The top two most abundant Gypsy families in S. alba preferentially insert in gene-poor regions. They are under relaxed epigenetic repression, probably due to the presence of CHROMO domains in their 3′-ends. Although a considerable number of TEs in S. alba showed differential expression under salt stress, only four copies were significantly correlated with their nearby genes in expression levels. One such TE-gene pair involves Abscisic acid 8'-hydroxylase 3 functioning in abscisic acid catabolism. This study sheds light on the evolutionary dynamics and potential function of TEs in an extremophile. Our results suggest that the conclusion on co-option of TEs should be cautious even though activation of TEs by stress might be prevalent.

Identification, Mapping, and Genetic Diversity of Novel Conserved Cross-Species Epitopes of RhopH2 in Plasmodium knowlesi With Plasmodium vivax

Malaria is a major public health concern, and any tangible intervention during the pre-elimination phase can result in a significant reduction in infection rates. Recent studies have reported that antigens producing cross-protective immunity can play an important role as vaccines and halt malaria transmission in different endemic regions. In this study, we studied the genetic diversity, natural selection, and discovered novel conserved epitopes of a high molecular weight rhoptry protein 2 (RhopH2) in clinical samples of Plasmodium knowlesi and Plasmodium vivax cross-protective domains, which has been proven to produce cross-protective immunity in both species. We found low levels of nucleotide diversity (P. knowlesi; π ~ 0.0093, SNPs = 49 and P. vivax π ~ 0.0014, SNPs = 23) in P. knowlesi (n = 40) and P. vivax (n = 65) samples in the PkRhopH2 cross-protective domain. Strong purifying selection was observed for both species (P. knowlesi; dS - dN = 2.41, p < 0.009, P. vivax; dS - dN = 1.58, p < 0.050). In silico epitope prediction in P. knowlesi identified 10 potential epitopes, of which 7 epitopes were 100% conserved within clinical samples. Of these epitopes, an epitope with 10 amino acids (QNSKHFKKEK) was found to be fully conserved within all P. knowlesi and P. vivax clinical samples and 80%–90% conservation within simian malaria ortholog species, i.e., P. coatneyi and P. cynomolgi. Phylogenetic analysis of the PkRhopH2 cross-protective domain showed geographical clustering, and three subpopulations of P. knowlesi were identified of which two subpopulations originated from Sarawak, Malaysian Borneo, and one comprised only the laboratory lines from Peninsular Malaysia. This study suggests that RhopH2 could be an excellent target for cross-protective vaccine development with potential for outwitting strain as well as species-specific immunity. However, more detailed studies on genetic diversity using more clinical samples from both species as well as the functional role of antibodies specific to the novel conserved epitope identified in this study can be explored for protection against infection.

Proteome evaluation of homolog abundance patterns in Arachis hypogaea cv. Tifrunner

Background Cultivated peanut (Arachis hypogaea, AABB genome), an allotetraploid from a cross between A. duranensis (AA genome) and A. ipaensis (BB genome), is an important oil and protein crop with released genome and RNA-seq sequence datasets. These datasets provide the molecular foundation for studying gene expression and evolutionary patterns. However, there are no reports on the proteomic data of A. hypogaea cv. Tifrunner, which limits understanding of its gene function and protein level evolution. Results This study sequenced the A. hypogaea cv. Tifrunner leaf and root proteome using the tandem mass tag technology. A total of 4803 abundant proteins were identified. The 364 differentially abundant proteins were estimated by comparing protein abundances between leaf and root proteomes. The differentially abundant proteins enriched the photosystem process. The number of biased abundant homeologs between the two sub-genomes A (87 homeologs in leaf and root) and B (69 and 68 homeologs in leaf and root, respectively) was not significantly different. However, homeologous proteins with biased abundances in different sub-genomes enriched different biological processes. In the leaf, homeologs biased to sub-genome A enriched biosynthetic and metabolic process, while homeologs biased to sub-genome B enriched iron ion homeostasis process. In the root, homeologs with biased abundance in sub-genome A enriched inorganic biosynthesis and metabolism process, while homeologs with biased abundance in sub-genome B enriched organic biosynthesis and metabolism process. Purifying selection mainly acted on paralogs and homeologs. The selective pressure values were negatively correlated with paralogous protein abundance. About 77.42% (24/31) homeologous and 80% (48/60) paralogous protein pairs had asymmetric abundance, and several protein pairs had conserved abundances in the leaf and root tissues. Conclusions This study sequenced the proteome of A. hypogaea cv. Tifrunner using the leaf and root tissues. Differentially abundant proteins were identified, and revealed functions. Paralog abundance divergence and homeolog bias abundance was elucidated. These results indicate that divergent abundance caused retention of homologs in A. hypogaea cv. Tifrunner.

A (GCC) Repeat in the Untranslated Region of Human SBF1 Departs from Hardy-Weinberg Equilibrium in Human and Links to Late-onset Neurocognitive Disorder.

The human SBF1 (SET binding factor 1) gene, alternatively known as MTMR5, is predominantly expressed in the brain, and its epigenetic dysregulation is linked to late-onset neurocognitive disorders (NCDs), such as Alzheimer’s disease. This gene contains a (GCC)-repeat at the interval between +1 and +60 of the transcription start site (SBF1-202 ENST00000380817.8). Sequencing of the SBF1 (GCC)-repeat in a sample of 542 Iranian individuals, consisting of late-onset NCDs (N=260) and controls (N=282) revealed a predominantly bi-allelic locus for this STR, consisting of 8 and 9 repeats, with allele frequencies ranging from 0.39 to 0.55, and four other alleles with frequencies of <0.03 across the two groups. Overall heterozygosity for the observed alleles was significantly less than expected in the NCD and control groups, at 22.3% and 16.31%, respectively (p=0.000). Specifically, the heterozygous 8/9 genotype was significantly less than expected in both case and control groups (Hardy-Weinberg disequilibrium, p=0.000), and significantly enriched in the NCD group (Yates corrected p=0.001). Skewed heterozygous genotypes were also detected for other allele combinations, such as 6/8 vs 6/9 across groups (p=0.000). Bioinformatics studies revealed that the number of (GCC)-repeats may change the RNA secondary structure and interaction sites across human exon 1. This STR was specifically expanded beyond 2-repeats in primates. In conclusion, we report a novel biological phenomenon in which there is indication of purifying selection against heterozygous genotypes at a STR locus in human, and skewed genotype compartment in late-onset NCD vs. controls. In view of the location of this STR in the 5′ UTR, RNA/RNA or RNA/DNA heterodimer formation of the involved genotypes and possible deleterious downstream events should be considered.

Multiple horizontal acquisitions of plant genes in the whitefly Bemisia tabaci

The extent to which horizontal gene transfer (HGT) has shaped eukaryote evolution remains an open question. Two recent studies reported four plant-like genes acquired through two HGT events by the whitefly Bemisia tabaci, a major agricultural pest (Lapadula et al. 2020; Xia et al. 2021). Here, we performed a systematic search for plant-to-insect HGT in B. tabaci and uncovered a total of 50 plant-like genes deriving from at least 24 independent HGT events. Most of these genes are present in three cryptic B. tabaci species, show high level of amino-acid identity to plant genes (mean = 64%), are phylogenetically nested within plant sequences, and are expressed and evolve under purifying selection. The predicted functions of these genes suggest that most of them are involved in plant-insect interactions. Thus, substantial plant-to-insect HGT may have facilitated the evolution of B. tabaci towards adaptation to a large host spectrum. Our study shows that eukaryote-to-eukaryote HGT may be relatively common in some lineages and it provides new candidate genes that may be targeted to improve current control strategies against whiteflies.

Mutation bias reflects natural selection in Arabidopsis thaliana

Since the first half of the twentieth century, evolutionary theory has been dominated by the idea that mutations occur randomly with respect to their consequences1. Here we test this assumption with large surveys of de novo mutations in the plant Arabidopsis thaliana. In contrast to expectations, we find that mutations occur less often in functionally constrained regions of the genome—mutation frequency is reduced by half inside gene bodies and by two-thirds in essential genes. With independent genomic mutation datasets, including from the largest Arabidopsis mutation accumulation experiment conducted to date, we demonstrate that epigenomic and physical features explain over 90% of variance in the genome-wide pattern of mutation bias surrounding genes. Observed mutation frequencies around genes in turn accurately predict patterns of genetic polymorphisms in natural Arabidopsis accessions (r = 0.96). That mutation bias is the primary force behind patterns of sequence evolution around genes in natural accessions is supported by analyses of allele frequencies. Finally, we find that genes subject to stronger purifying selection have a lower mutation rate. We conclude that epigenome-associated mutation bias2 reduces the occurrence of deleterious mutations in Arabidopsis, challenging the prevailing paradigm that mutation is a directionless force in evolution.

Molecular Characterization Revealed the Role of Thaumatin-Like Proteins of Bread Wheat in Stress Response

Thaumatin-like proteins (TLPs) are related to pathogenesis-related-5 (PR-5) family and involved in stress response. Herein, a total of 93 TLP genes were identified in the genome of Triticum aestivum. Further, we identified 26, 27, 39, and 37 TLP genes in the Brachypodium distachyon, Oryza sativa, Sorghum bicolor, and Zea mays genomes for comparative characterization, respectively. They could be grouped into small and long TLPs with conserved thaumatin signature motif. Tightly clustered genes exhibited conserved gene and protein structure. The physicochemical analyses suggested significant differences between small and long TLPs. Evolutionary analyses suggested the role of duplication events and purifying selection in the expansion of the TLP gene family. Expression analyses revealed the possible roles of TLPs in plant development and abiotic and fungal stress response. Recombinant expression of TaTLP2-B in Saccharomyces cerevisiae provided significant tolerance against cold, heat, osmotic, and salt stresses. The results depicted the importance of TLPs in cereal crops that would be highly useful in future crop improvement programs.

RAD gene family analysis in cotton provides some key genes for flowering and stress tolerance in upland cotton G. hirsutum

Background RADIALIS (RAD), belongs to the MYB gene family and regulates a variety of functions including floral dorsoventral asymmetry in Antirrhinum majus and development of fruit proteins in Solanum lycopersicum. RAD genes contain an SNF2_N superfamily domain. Here, we comprehensively identified 68 RAD genes from six different species including Arabidopsis and five species of cotton. Results Phylogenetic analysis classified RAD genes into five groups. Gene structure, protein motifs and conserved amino acid residues indicated that GhRAD genes were highly conserved during the evolutionary process. Chromosomal location information showed that GhRAD genes were distributed unevenly on different chromosomes. Collinearity and selection pressure analysis indicated RAD gene family expansion in G. hirsutum and G. barbadense with purifying selection pressure. Further, various growth and stress related promotor cis-acting elements were observed. Tissue specific expression level indicated that most GhRAD genes were highly expressed in roots and flowers (GhRAD2, GhRAD3, GhRAD4 and GhRAD11). Next, GhRAD genes were regulated by phytohormonal stresses (JA, BL and IAA). Moreover, Ghi-miRN1496, Ghi-miR1440, Ghi-miR2111b, Ghi-miR2950a, Ghi-miR390a, Ghi-miR390b and Ghi-miR7495 were the miRNAs targeting most of GhRAD genes. Conclusions Our study revealed that RAD genes are evolutionary conserved and might be involved in different developmental processes and hormonal stress response. Data presented in our study could be used as the basis for future studies of RAD genes in cotton.

Complete Chloroplast Genome Sequence of Fagus longipetiolata Seemen (Fagaceae): Genome Structure, Adaptive Evolution, and Phylogenetic Relationships

Fagus longipetiolata Seemen is a deciduous tree of the Fagus genus in Fagaceae, which is endemic to China. In this study, we successfully sequenced the cp genome of F. longipetiolata, compared the cp genomes of the Fagus genus, and reconstructed the phylogeny of Fagaceae. The results showed that the cp genome of F. longipetiolata was 158,350 bp, including a pair of inverted repeat (IRA and IRB) regions with a length of 25,894 bp each, a large single-copy (LSC) region of 87,671 bp, and a small single-copy (SSC) region of 18,891 bp. The genome encoded 131 unique genes, including 81 protein-coding genes, 37 transfer RNA genes (tRNAs), 8 ribosomal RNA genes (rRNAs), and 5 pseudogenes. In addition, 33 codons and 258 simple sequence repeats (SSRs) were identified. The cp genomes of Fagus were relatively conserved, especially the IR regions, which showed the best conservation, and no inversions or rearrangements were found. The five regions with the largest variations were the rps12, rpl32, ccsA, trnW-CCA, and rps3 genes, which spread over in LSC and SSC. The comparison of gene selection pressure indicated that purifying selection was the main selective pattern maintaining important biological functions in Fagus cp genomes. However, the ndhD, rpoA, and ndhF genes of F. longipetiolata were affected by positive selection. Phylogenetic analysis revealed that F. longipetiolata and F. engleriana formed a close relationship, which partially overlapped in their distribution in China. Our analysis of the cp genome of F. longipetiolata would provide important genetic information for further research into the classification, phylogeny and evolution of Fagus.

Fragmentation in mitochondrial genomes in relation to elevated sequence divergence and extreme rearrangements

Background A single circular mitochondrial (mt) genome is a common feature across most metazoans. The mt-genome includes protein-coding genes involved in oxidative phosphorylation, as well as RNAs necessary for translation of mt-RNAs, whose order and number are highly conserved across animal clades, with few known exceptions of alternative mt-gene order or mt-genome architectures. One such exception consists of the fragmented mitochondrial genome, a type of genome architecture where mt-genes are split across two or more mt-chromosomes. However, the origins of mt-genome fragmentation and its effects on mt-genome evolution are unknown. Here, we investigate these origin and potential mechanisms underlying mt-genome fragmentation, focusing on a genus of booklice, Liposcelis, which exhibits elevated sequence divergence, frequent rearrangement of mt-gene order, and fragmentation of the mt genome, and compare them to other Metazoan clades. Results We found this genus Liposcelis exhibits very low conservation of mt-gene order across species, relative to other metazoans. Levels of gene order rearrangement were, however, unrelated to whether or not mt-genomes were fragmented or intact, suggesting mitochondrial genome fragmentation is not affecting mt-gene order directly. We further investigated possible mechanisms underpinning these patterns and revealed very high conservation of non-coding sequences at the edges of multiple recombination regions across populations of one particular Liposcelis species, supportive of a hypothesis that mt-fragmentation arises from recombination errors between mt-genome copies. We propose these errors may arise as a consequence of a heightened mutation rate in clades exhibiting mt-fragmentation. Consistent with this, we observed a striking pattern across three Metazoan phyla (Arthropoda, Nematoda, Cnidaria) characterised by members exhibiting high levels of mt-gene order rearrangement and cases of mt-fragmentation, whereby the mt-genomes of species more closely related to species with fragmented mt-genomes diverge more rapidly despite experiencing strong purifying selection. Conclusions We showed that contrary to expectations, mt-genome fragmentation is not correlated with the increase in mt-genome rearrangements. Furthermore, we present evidence that fragmentation of the mt-genome may be part of a general relaxation of a natural selection on the mt-genome, thus providing new insights into the origins of mt-genome fragmentation and evolution.