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

A Chromosome-level Genome Assembly of the Highly Heterozygous Sea Urchin Echinometra sp. EZ Reveals Adaptation in the Regulatory Regions of Stress Response Genes

Echinometra is the most widespread genus of sea urchin and has been the focus of a wide range of studies in ecology, speciation, and reproduction. However, available genetic data for this genus are generally limited to a few select loci. Here, we present a chromosome-level genome assembly based on 10x Genomics, PacBio, and Hi-C sequencing for Echinometra sp. EZ from the Persian/Arabian Gulf. The genome is assembled into 210 scaffolds totaling 817.8 Mb with an N50 of 39.5 Mb. From this assembly, we determined that the E. sp. EZ genome consists of 2n = 42 chromosomes. BUSCO analysis showed that 95.3% of BUSCO genes were complete. Ab initio and transcript-informed gene modeling and annotation identified 29,405 genes, including a conserved Hox cluster. E. sp. EZ can be found in high-temperature and high-salinity environments, and we therefore compared E. sp. EZ gene families and transcription factors associated with environmental stress response ("defensome") with other echinoid species with similar high-quality genomic resources. While the number of defensome genes was broadly similar for all species, we identified strong signatures of positive selection in E. sp. EZ noncoding elements near genes involved in environmental response pathways as well as losses of transcription factors important for environmental response. These data provide key insights into the biology of E. sp. EZ as well as the diversification of Echinometra more widely and will serve as a useful tool for the community to explore questions in this taxonomic group and beyond.

Genome-wide characterization and function analysis uncovered roles of wheat LIMs in responding to adverse stresses and TaLIM8-4D function as a susceptible gene

The Lin-11, Isl-1, and Mec-3 domains (LIM) transcription factors play essential roles in regulating plant biological processes. Despite that, there is a lack of a full understanding of LIMs in wheat (Triticum aestivum L.). In this study, 28 wheat LIM s (TaLIMs) were identified and designated as TaLIM1-1A to TaLIM12-7D. The cis-regulatory element analysis showed that TaLIMs were rich in elements related to biological and abiotic stresses. Expression profiling analysis showed that certain members of TaLIMs were responsive to biotic and abiotic stresses, such as TaLIM1-1A, TaLIM3-2B, TaLIM8-4D, and TaLIM10-5D, were significantly induced by heat, drought, sodium chloride (NaCl), abscisic acid (ABA) and Fusarium graminearum stresses. Furthermore, the biological function of TaLIM8-4D was analyzed and results showed that it was subcellular localization in the nucleus and could induce weak cell death in Nicotiana benthamiana leaves. Additionally, overexpression of TaLIM8-4D could upregulate plant pathogenesis-related (PR) genes, promoting the infection of hemibiotrophic pathogen, implying that TaLIM8-4D could function as susceptible gene in the nucleus by upregulating PR genes and inducing cell death to promote the colonization of hemibiotrophic agent F. graminearum. Overall, the systematic identification, characterization, expression profiling, evolutionary, and function analyses provided the ability to understand TaLIMs and laid a foundation for the further function study of LIM family members in wheat.

A newly evolved chimeric lysin motif receptor-like kinase in Medicago truncatula spp. tricycla R108 extends its Rhizobia symbiotic partnership

Rhizobial lipochitooligosaccharidic Nod factors (NFs), specified by nod genes, are the primary determinants of host specificity in the legume-Rhizobia symbiosis. We examined the nodulation ability of Medicago truncatula cv Jemalong A17 and M. truncatula ssp. tricycla R108 with the Sinorhizobium meliloti nodF/nodL mutant, which produces modified NFs. We then applied genetic and functional approaches to study the genetic basis and mechanism of nodulation of R108 by this mutant. We show that the nodF/nodL mutant can nodulate R108 but not A17. Using genomics and reverse genetics, we identified a newly evolved, chimeric LysM receptor-like kinase gene in R108, LYK2bis, which is responsible for the phenotype and can allow A17 to gain nodulation with the nodF/nodL mutant. We found that LYK2bis is involved in nodulation by mutants producing nonO-acetylated NFs and interacts with the key receptor protein NFP. Many, but not all, natural S. meliloti and S. medicae strains tested require LYK2bis for efficient nodulation of R108. Our findings reveal that a newly evolved gene in R108, LYK2bis, extends nodulation specificity to mutants producing nonO-acetylated NFs and is important for nodulation by many natural Sinorhizobia. Evolution of this gene may present an adaptive advantage to allow nodulation by a greater variety of strains.

Duplication of a manganese/cadmium transporter gene reduces cadmium accumulation in rice grain

Global contamination of soils with toxic cadmium (Cd) is a serious health threat. Here we found that a tandem duplication of a gene encoding a manganese/Cd transporter, OsNramp5, was responsible for low-Cd accumulation in Pokkali, an old rice cultivar. This duplication doubled the expression of OsNramp5 gene but did not alter its spatial expression pattern and cellular localization. Higher expression of OsNramp5 increased uptake of Cd and Mn into the root cells but decreased Cd release to the xylem. Introgression of this allele into Koshihikari, an elite rice cultivar, through backcrossing significantly reduced Cd accumulation in the grain when cultivated in soil heavily contaminated with Cd but did not affect both grain yield and eating quality. This study not only reveals the molecular mechanism underlying low-Cd accumulation but also provides a useful target for breeding rice cultivars with low-Cd accumulation.

Genome-wide identification and characterization of the bHLH gene family and analysis of their potential relevance to chlorophyll metabolism in Raphanus sativus L

BACKGROUND

Green-fleshed radish (Raphanus sativus L.) is an economically important root vegetable of the Brassicaceae family, and chlorophyll accumulates in its root tissues. It was reported that the basic helix-loop-helix (bHLH) transcription factors play vital roles in the process of chlorophyll metabolism. Nevertheless, a comprehensive study on the bHLH gene family has not been performed in Raphanus sativus L.

RESULTS

In this study, a total of 213 Raphanus sativus L. bHLH (RsbHLH) genes were screened in the radish genome, which were grouped into 22 subfamilies. 204 RsbHLH genes were unevenly distributed on nine chromosomes, and nine RsbHLH genes were located on the scaffolds. Gene structure analysis showed that 25 RsbHLH genes were intron-less. Collineation analysis revealed the syntenic orthologous bHLH gene pairs between radish and Arabidopsis thaliana/Brassica rapa/Brassica oleracea. 162 RsbHLH genes were duplicated and retained from the whole genome duplication event, indicating that the whole genome duplication contributed to the expansion of the RsbHLH gene family. RNA-seq results revealed that RsbHLH genes had a variety of expression patterns at five development stages of green-fleshed radish and white-fleshed radish. In addition, the weighted gene co-expression network analysis confirmed four RsbHLH genes closely related to chlorophyll content.

CONCLUSIONS

A total of 213 RsbHLH genes were identified, and we systematically analyzed their gene structure, evolutionary and collineation relationships, conserved motifs, gene duplication, cis-regulatory elements and expression patterns. Finally, four bHLH genes closely involved in chlorophyll content were identified, which may be associated with the photosynthesis of the green-fleshed radish. The current study would provide valuable information for further functional exploration of RsbHLH genes, and facilitate clarifying the molecular mechanism underlying photosynthesis process in green-fleshed radish.

Convergent copy number increase of genes associated with freshwater colonization in fishes

Copy number variation (CNV) can cause phenotypic changes. However, in contrast to amino acid substitutions and cis-regulatory changes, little is known about the functional categories of genes in which CNV is important for adaptation to novel environments. It is also unclear whether the same genes repeatedly change the copy numbers for adapting to similar environments. Here, we investigate CNV associated with freshwater colonization in fishes, which was observed multiple times across different lineages. Using 48 ray-finned fishes across diverse orders, we identified 23 genes whose copy number increases were associated with freshwater colonization. These genes showed enrichment for peptide receptor activity, hexosyltransferase activity and unsaturated fatty acid metabolism. We further revealed that three of the genes showed copy number increases in freshwater populations compared to marine ancestral populations of the stickleback genus Gasterosteus. These results indicate that copy number increases of genes involved in fatty acid metabolism (FADS2), immune function (PSMB8a) and thyroid hormone metabolism (UGT2) may be important for freshwater colonization at both the inter-order macroevolutionary scale and at the intra-genus microevolutionary scale. Further analysis across diverse taxa will help to understand the role of CNV in the adaptation to novel environments. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Does genome size increase with water depth in marine fishes?

A growing body of research suggests that genome size in animals can be affected by ecological factors. Half a century ago, Ebeling et al. proposed that genome size increases with depth in some teleost fish groups and discussed a number of biological mechanisms that may explain this pattern (e.g., passive accumulation, adaptive acclimation). Using phylogenetic comparative approaches, we revisit this hypothesis based on genome size and ecological data from up to 708 marine fish species in combination with a set of large-scale phylogenies, including a newly inferred tree. We also conduct modeling approaches of trait evolution and implement a variety of regression analyses to assess the relationship between genome size and depth. Our reanalysis of Ebeling et al.'s dataset shows a weak association between these variables, but the overall pattern in their data is driven by a single clade. Although new analyses based on our "all-species" dataset resulted in positive correlations, providing some evidence that genome size evolves as a function of depth, only one subclade consistently yielded statistically significant correlations. By contrast, negative correlations are rare and nonsignificant. All in all, we find modest evidence for an increase in genome size along the depth axis in marine fishes. We discuss some mechanistic explanations for the observed trends.

Differential Expression of Two Copies of the irmA Gene in the Enteroaggregative E. coli Strain 042

Gene duplications significantly impact the gene repertoires of both eukaryotic and prokaryotic microorganisms. The genomes of pathogenic Escherichia coli strains share a group of duplicated genes whose function is mostly unknown. The irmA gene is one of the duplicates encoded in several pathogenic E. coli strains. The function of its gene product was investigated in the uropathogenic E. coli strain CFT073, which contains a single functional copy. The IrmA protein structure mimics that of human interleukin receptors and likely plays a role during infection. The enteroaggregative E. coli strain 042 contains two functional copies of the irmA gene. In the present work, we investigated their biological roles. The irmA_4509 allele is expressed under several growth conditions. Its expression is modulated by the global regulators OxyR and Hha, with optimal expression at 37°C and under nutritional stress conditions. Expression of the irmA_2244 allele can only be detected when the irmA_4509 allele is knocked out. Differences in the promoter regions of both alleles account for their differential expression. Our results show that under several environmental conditions, the expression of the IrmA protein in strain 042 is dictated by the irmA_4509 allele. The irmA_2244 allele appears to play a backup role to ensure IrmA expression when the irmA_4509 allele loses its function. IMPORTANCE Gene duplications occur in prokaryotic genomes at a detectable frequency. In many instances, the biological function of the duplicates is unknown, and hence, the significance of the presence of multiple copies of these genes remains unclear. In pathogenic E. coli isolates, the irmA gene can be present either as a single copy or in two or more copies. We focused our work on studying why a different pathogenic E. coli strain encodes two functional copies of the irmA gene. We show that under several environmental conditions, one of the alleles dictates IrmA expression, and the second remains silent. The latter allele is only expressed when the former is silenced. The presence of more than one functional copy of the irmA gene in some pathogenic E. coli strains can result in sufficient expression of this virulence factor during the infection process.

Adaptative Mechanisms of Halophytic Eutrema salsugineum Encountering Saline Environment

Salt cress (Eutrema salsugineum), an Arabidopsis-related halophyte, can naturally adapt to various harsh climates and soil conditions; thus, it is considered a desirable model plant for deciphering mechanisms of salt and other abiotic stresses. Accumulating evidence has revealed that compared with Arabidopsis, salt cress possesses stomata that close more tightly and more succulent leaves during extreme salt stress, a noticeably higher level of proline, inositols, sugars, and organic acids, as well as stress-associated transcripts in unstressed plants, and they are induced rapidly under stress. In this review, we systematically summarize the research on the morphology, physiology, genome, gene expression and regulation, and protein and metabolite profile of salt cress under salt stress. We emphasize the latest advances in research on the genome adaptive evolution encountering saline environments, and epigenetic regulation, and discuss the mechanisms underlying salt tolerance in salt cress. Finally, we discuss the existing questions and opportunities for future research in halophytic Eutrema. Together, the review fosters a better understanding of the mechanism of plant salt tolerance and provides a reference for the research and utilization of Eutrema as a model extremophile in the future. Furthermore, the prospects for salt cress applied to explore the mechanism of salt tolerance provide a theoretical basis to develop new strategies for agricultural biotechnology.

Impact of Modified Atmospheres on Growth and Metabolism of Meat-Spoilage Relevant Photobacterium spp. as Predicted by Comparative Proteomics

Modified atmosphere packaging (MAP) is a common strategy to selectively prevent the growth of certain species of meat spoiling bacteria. This study aimed to determine the impact of high oxygen MAP (70% O₂, 30% CO₂, red and white meats) and oxygen-free MAP (70% N₂, 30% CO₂, also white meat and seafood) on preventing the growth of spoiling photobacteria on meat. Growth of Photobacterium carnosum and P. phosphoreum was monitored in a meat simulation media under different gas mixtures of nitrogen, oxygen, and carbon dioxide, and samples were taken during exponential growth for a comparative proteomic analysis. Growth under air atmosphere appears optimal, particularly for P. carnosum. Enhanced protein accumulation affected energy metabolism, respiration, oxygen consuming reactions, and lipid usage. However, all the other atmospheres show some degree of growth reduction. An increase in oxygen concentration leads to an increase in enzymes counteracting oxidative stress for both species and enhancement of heme utilization and iron-sulfur cluster assembly proteins for P. phosphoreum. Absence of oxygen appears to switch the metabolism toward fermentative pathways where either ribose (P. phosphoreum) or glycogen (P. carnosum) appear to be the preferred substrates. Additionally, it promotes the use of alternative electron donors/acceptors, mainly formate and nitrate/nitrite. Stress response is manifested as an enhanced accumulation of enzymes that is able to produce ammonia (e.g., carbonic anhydrase, hydroxylamine reductase) and regulate osmotic stress. Our results suggest that photobacteria do not sense the environmental levels of carbon dioxide, but rather adapt to their own anaerobic metabolism. The regulation in presence of carbon dioxide is limited and strain-specific under anaerobic conditions. However, when oxygen at air-like concentration (21%) is present together with carbon dioxide (30%), the oxidative stress appears enhanced compared to air conditions (very low carbon dioxide), as explained if both gases have a synergistic effect. This is further supported by the increase in oxygen concentration in the presence of carbon dioxide. The atmosphere is able to fully inhibit P. carnosum, heavily reduce P. phosphoreum growth in vitro, and trigger diversification of energy production with higher energetic cost, highlighting the importance of concomitant bacteria for their growth on raw meat under said atmosphere.

Genome-Wide Identification and Analysis of the Class III Peroxidase Gene Family in Tobacco (Nicotiana tabacum)

Class III peroxidases (PODs) are plant-specific enzymes that play significant roles in plant physiological processes and stress responses. However, a comprehensive analysis of the POD gene family in tobacco has not yet been conducted. In this study, 210 non-redundant POD gene members (NtPODs) were identified in tobacco (Nicotiana tabacum) and distributed unevenly throughout 24 tobacco chromosomes. Phylogenetic analysis clustered these genes into six subgroups (I-VI). Gene structure and motif analyses showed the structural and functional diversity among the subgroups. Segmental duplication and purifying selection were the main factors affecting NtPOD gene evolution. Our analyses also suggested that NtPODs might be regulated by miRNAs and cis-acting regulatory elements of transcription factors that are involved in various biological processes. In addition, the expression patterns in different tissues and under various stress treatments were investigated. The results showed that the majority of NtPODs had tissue-specific expression patterns and may be involved in many biotic and abiotic responses. qRT-PCR analyses of different tissues and stress treatments were performed to verify transcriptome patterns. Expression of a green fluorescent protein-NtPOD fusion confirmed the plasma membrane localization of NtPOD121 and NtPOD4. Furthermore, 3D structures provided evidences of membrane-bound peroxidase. These findings provide useful information to better understand the evolution of the NtPOD gene family and lay the foundation for further studies on POD gene function in tobacco.

Exploring the SiCCT Gene Family and Its Role in Heading Date in Foxtail Millet

CCT transcription factors are involved in the regulation of photoperiod and abiotic stress in Arabidopsis and rice. It is not clear that how CCT gene family expand and regulate heading date in foxtail millet. In this study, we conducted a systematic analysis of the CCT gene family in foxtail millet. Thirty-nine CCT genes were identified and divided into four subfamilies based on functional motifs. Analysis showed that dispersed duplication played a predominant role in the expansion of CCT genes during evolution. Nucleotide diversity analysis suggested that genes in CONSTANS (COL)-like, CCT MOTIF FAMILY (CMF)-like, and pseudoresponse response regulator (PRR)-like subfamilies were subjected to selection. Fifteen CCT genes were colocalized with previous heading date quantitative trait loci (QTL) and genome-wide association analysis (GWAS) signals. Transgenic plants were then employed to confirm that overexpression of the CCT gene SiPRR37 delayed the heading date and increased plant height. Our study first investigated the characterization and expansion of the CCT family in foxtail millet and demonstrated the role of SiPRR37. These results lay a significant foundation for further research on the function of CCT genes and provide a cue for the regulation of heading date.

The fate of drought-related genes after polyploidization in Arachis hypogaea cv. Tifrunner

Drought stress affects plant growth and development. Cultivated peanut (Arachis hypogaea) was formed by a cross between A. duranensis and A. ipaensis. The drought tolerance of A. duranensis and A. ipaensis is reportedly stronger than that of cultivated peanut. However, there has been little study of drought tolerance genes in Arachis. In this study, we compared drought tolerance genes between A. hypogaea cv. Tifrunner and its diploid donors. We have observed that polyploidization does not generate more drought tolerance genes in A. hypogaea cv. Tifrunner but promotes the loss of many ancient drought tolerance genes. Although putative drought tolerance genes occurred on gene duplication events in A. hypogaea cv. Tifrunner, most copies lacked drought tolerance. These findings suggest that the loss of drought tolerance genes in A. hypogaea cv. Tifrunner could possibly result in weaker drought tolerance. In addition, we have observed that the three Arachis species stochastically lost putative drought tolerance genes. The evolution of drought tolerance genes could possibly have correlated with environmental changes. Our results enhance the current understanding of drought tolerance and polyploidy evolution in Arachis species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01198-0.

Advantages and Challenges of Differential Immune Cell Count Determination in Blood and Milk for Monitoring the Health and Well-Being of Dairy Cows

A key challenge of the 21st century will be to provide the growing world population with a sustainable and secure supply of food. Consequently, the dairy farming’s primary task is to lower milk losses and other inefficiencies associated with diseased cows. Moreover, a shift from curative to preventive health management would be desirable for mastitis and a wide variety of other infectious and non-infectious cattle diseases, some of which are known to have profound negative effects on the performance and well-being of cows. Differential cell counting (DCC), a procedure that aims to determine the proportions of different somatic cell types in raw milk samples, has not only the potential to optimize mastitis diagnostics, but it could furthermore serve as a diagnostic tool for monitoring the general and overall health status of dairy cows. Based on a broad search of the literature, the practical utility of various types of DCC is summarized and discussed in this review. Since it might be of advantage to interpret DCC with the aid of data from studies in humans, differences between the immune systems of humans and dairy cattle, with a special focus on surface marker expression profiles and γδ (gamma delta) T-cell characteristics, are also described.

Long-read and chromosome-scale assembly of the hexaploid wheat genome achieves high resolution for research and breeding

BACKGROUND

The sequencing of the wheat (Triticum aestivum) genome has been a methodological challenge for many years owing to its large size (15.5 Gb), repeat content, and hexaploidy. Many initiatives aiming at obtaining a reference genome of cultivar Chinese Spring have been launched in the past years and it was achieved in 2018 as the result of a huge effort to combine short-read sequencing with many other resources. Reference-quality genome assemblies were then produced for other accessions, but the rapid evolution of sequencing technologies offers opportunities to reach high-quality standards at lower cost.

RESULTS

Here, we report on an optimized procedure based on long reads produced on the Oxford Nanopore Technology PromethION device to assemble the genome of the French bread wheat cultivar Renan.

CONCLUSIONS

We provide the most contiguous chromosome-scale assembly of a bread wheat genome to date. Coupled with an annotation based on RNA-sequencing data, this resource will be valuable for the crop community and will facilitate the rapid selection of agronomically important traits. We also provide a framework to generate high-quality assemblies of complex genomes using ONT.

Venom gene sequence diversity and expression jointly shape diet adaptation in pitvipers

Understanding the joint roles of protein sequence variation and differential expression during adaptive evolution is a fundamental, yet largely unrealized goal of evolutionary biology. Here, we use phylogenetic path analysis to analyze a comprehensive venom-gland transcriptome dataset spanning three genera of pitvipers to identify the functional genetic basis of a key adaptation (venom complexity) linked to diet breadth (DB). The analysis of gene-family-specific patterns reveals that, for genes encoding two of the most important venom proteins (snake venom metalloproteases and snake venom serine proteases), there are direct, positive relationships between sequence diversity (SD), expression diversity (ED), and increased DB. Further analysis of gene-family diversification for these proteins showed no constraint on how individual lineages achieved toxin gene SD in terms of the patterns of paralog diversification. In contrast, another major venom protein family (PLA₂s) showed no relationship between venom molecular diversity and DB. Additional analyses suggest that other molecular mechanisms—such as higher absolute levels of expression—are responsible for diet adaptation involving these venom proteins. Broadly, our findings argue that functional diversity generated through sequence and expression variations jointly determine adaptation in the key components of pitviper venoms, which mediate complex molecular interactions between the snakes and their prey.

Adaptive Variation in Homolog Number Within Transcript Families Promotes Expression Divergence in Reef-Building Coral

Gene expression, especially in multi-species experiments, is used to gain insight into the genetic basis of how organisms adapt and respond to changing environments. However, evolutionary processes which can influence gene expression patterns between species such as the presence of paralogs which arise from gene duplication events are rarely accounted for. Paralogous transcripts can alter the transcriptional output of a gene and thus exclusion of these transcripts can obscure important biological differences between species. To address this issue, we investigated how differences in transcript family size is associated with divergent gene expression patterns in five species of Caribbean reef-building corals. We demonstrate that transcript families that are rapidly evolving in terms of size have increased levels of expression divergence. Additionally, these rapidly evolving transcript families are enriched for multiple biological processes, with genes involved in the coral innate immune system demonstrating pronounced variation in homolog number between species. Overall, this investigation demonstrates the importance of incorporating paralogous transcripts when comparing gene expression across species by influencing both transcriptional output and the number of transcripts within biological processes. As this investigation was based on transcriptome assemblies, additional insights into the relationship between gene duplications and expression patterns will likely emergence once more genome assemblies are available for study.

Genome-Wide Evolution and Comparative Analysis of Superoxide Dismutase Gene Family in Cucurbitaceae and Expression Analysis of Lagenaria siceraria Under Multiple Abiotic Stresses

Superoxide dismutase (SOD) is an important enzyme that serves as the first line of defense in the plant antioxidant system and removes reactive oxygen species (ROS) under adverse conditions. The SOD protein family is widely distributed in the plant kingdom and plays a significant role in plant growth and development. However, the comprehensive analysis of the SOD gene family has not been conducted in Cucurbitaceae. Subsequently, 43 SOD genes were identified from Cucurbitaceae species [Citrullus lanatus (watermelon), Cucurbita pepo (zucchini), Cucumis sativus (cucumber), Lagenaria siceraria (bottle gourd), Cucumis melo (melon)]. According to evolutionary analysis, SOD genes were divided into eight subfamilies (I, II, III, IV, V, VI, VII, VIII). The gene structure analysis exhibited that the SOD gene family had comparatively preserved exon/intron assembly and motif as well. Phylogenetic and structural analysis revealed the functional divergence of Cucurbitaceae SOD gene family. Furthermore, microRNAs 6 miRNAs were predicted targeting 3 LsiSOD genes. Gene ontology annotation outcomes confirm the role of LsiSODs under different stress stimuli, cellular oxidant detoxification processes, metal ion binding activities, SOD activity, and different cellular components. Promoter regions of the SOD family revealed that most cis-elements were involved in plant development, stress response, and plant hormones. Evaluation of the gene expression showed that most SOD genes were expressed in different tissues (root, flower, fruit, stem, and leaf). Finally, the expression profiles of eight LsiSOD genes analyzed by qRT-PCR suggested that these genetic reserves responded to drought, saline, heat, and cold stress. These findings laid the foundation for further study of the role of the SOD gene family in Cucurbitaceae. Also, they provided the potential for its use in the genetic improvement of Cucurbitaceae.

Characterization of SOD and GPX Gene Families in the Soybeans in Response to Drought and Salinity Stresses

Plant stresses causing accumulation of reactive oxidative species (ROS) are scavenged by effective antioxidant defense systems. Therefore, the present study performed genome-wide identification of superoxide dismutase (SOD) and glutathione peroxidase (GPX) gene families in cultivated and wild soybeans, and 11 other legume species. We identified a total of 101 and 95 genes of SOD and GPX, respectively, across thirteen legume species. The highest numbers of SODs and GPXs were identified in cultivated (Glycine max) and wild (Glycine soja). A comparative phylogenetic study revealed highest homology among the SODs and GPXs of cultivated and wild soybeans relative to other legumes. The exon/intron structure, motif and synteny blocks were conserved in both soybean species. According to Ka/Ks, purifying the selection played the major evolutionary role in these gene families, and segmental duplication are major driving force for SODs and GPXs expansion. In addition, the qRT-PCR analysis of the G. max and G. soja SOD and GPX genes revealed significant differential expression of these genes in response to oxidative, drought and salinity stresses in root tissue. In conclusion, our study provides new insights for the evolution of SOD and GPX gene families in legumes, and provides resources for further functional characterization of these genes for multiple stresses.

Genomic Features and Pervasive Negative Selection in Rhodanobacter Strains Isolated from Nitrate and Heavy Metal Contaminated Aquifer

Rhodanobacter species dominate in the Oak Ridge Reservation (ORR) subsurface environments contaminated with acids, nitrate, metal radionuclides, and other heavy metals. To uncover the genomic features underlying adaptations to these mixed-waste environments and to guide genetic tool development, we sequenced the whole genomes of eight Rhodanobacter strains isolated from the ORR site. The genome sizes ranged from 3.9 to 4.2 Mb harboring 3,695 to 4,035 protein-coding genes and GC contents approximately 67%. Seven strains were classified as R. denitrificans and one strain, FW510-R12, as R. thiooxydans based on full length 16S rRNA sequences. According to gene annotation, the top two Cluster of Orthologous Groups (COGs) with high pan-genome expansion rates (Pan/Core gene ratio) were “replication, recombination and repair” and “defense mechanisms.” The denitrifying genes had high DNA homologies except the predicted protein structure variances in NosZ. In contrast, heavy metal resistance genes were diverse with between 7 to 34% of them were located in genomic islands, and these results suggested origins from horizontal gene transfer. Analysis of the methylation patterns in four strains revealed the unique 5mC methylation motifs. Most orthologs (78%) had ratios of nonsynonymous to synonymous substitutions (dN/dS) less than one when compared to the type strain 2APBS1, suggesting the prevalence of negative selection. Overall, the results provide evidence for the important roles of horizontal gene transfer and negative selection in genomic adaptation at the contaminated field site. The complex restriction-modification system genes and the unique methylation motifs in Rhodanobacter strains suggest the potential recalcitrance to genetic manipulation.

IMPORTANCE Despite the dominance of Rhodanobacter species in the subsurface of the contaminated Oak Ridge Reservation (ORR) site, very little is known about the mechanisms underlying their adaptions to the various stressors present at ORR. Recently, multiple Rhodanobacter strains have been isolated from the ORR groundwater samples from several wells with varying geochemical properties. Using Illumina, PacBio, and Oxford Nanopore sequencing platforms, we obtained the whole genome sequences of eight Rhodanobacter strains. Comparison of the whole genomes demonstrated the genetic diversity, and analysis of the long nanopore reads revealed the heterogeneity of methylation patterns in strains isolated from the same well. Although all strains contained a complete set of denitrifying genes, the predicted tertiary structures of NosZ differed. The sequence comparison results demonstrate the important roles of horizontal gene transfer and negative selection in adaptation. In addition, these strains may be recalcitrant to genetic manipulation due to the complex restriction-modification systems and methylations.

Genome-Wide Survey Indicates Diverse Physiological Roles of Dendrobium officinale Calcium-Dependent Protein Kinase Genes

Calcium-dependent protein kinases (CDPKs) are crucial calcium ions (Ca²⁺) sensors in plants with important roles in signal transduction, plant growth, development, and stress responses. Here, we identified 24 genes encoding CDPKs in Dendrobium officinale using genome-wide analysis. The phylogenetic analysis revealed that these genes formed four groups, with similar structures in the same group. The gene expression patterns following hormone treatments and yeast two-hybrid of homologous CDPK gene pairs with Rbohs showed differences, indicating functional divergence between homologous genes. In addition, the rapid accumulation of hydrogen peroxide (H₂O₂) and stomatal closure was observed in response to salicylic acid (SA)/jasmonic acid (JA) stress. Our data showed that CDPK9-2 and CDPK20-4 interacted with Rboh D and Rboh H, respectively, and were implicated in the generation of H₂O₂ and regulation of the stomatal aperture in response to salicylic acid/jasmonic acid treatment. We believe these results can provide a foundation for the functional divergence of homologous genes in D. officinale.

A constitutive stress response is a result of low temperature growth in the Antarctic green alga Chlamydomonas sp. UWO241

The Antarctic green alga Chlamydomonas sp. UWO241 is an obligate psychrophile that thrives in the cold (4-6°C) but is unable to survive at temperatures ≥18°C. Little is known how exposure to heat affects its physiology or whether it mounts a heat stress response in a manner comparable to mesophiles. Here, we dissect the responses of UWO241 to temperature stress by examining its growth, primary metabolome and transcriptome under steady-state low temperature and heat stress conditions. In comparison with Chlamydomonas reinhardtii, UWO241 constitutively accumulates metabolites and proteins commonly considered as stress markers, including soluble sugars, antioxidants, polyamines, and heat shock proteins to ensure efficient protein folding at low temperatures. We propose that this results from life at extreme conditions. A shift from 4°C to a non-permissive temperature of 24°C alters the UWO241 primary metabolome and transcriptome, but growth of UWO241 at higher permissive temperatures (10 and 15°C) does not provide enhanced heat protection. UWO241 also fails to induce the accumulation of HSPs when exposed to heat, suggesting that it has lost the ability to fine-tune its heat stress response. Our work adds to the growing body of research on temperature stress in psychrophiles, many of which are threatened by climate change.

Identification and Comprehensive Genome-Wide Analysis of Glutathione S-Transferase Gene Family in Sweet Cherry (Prunus avium) and Their Expression Profiling Reveals a Likely Role in Anthocyanin Accumulation

Glutathione S-transferases (GSTs) in plants are multipurpose enzymes that are involved in growth and development and anthocyanins transportation. However, members of the GST gene family were not identified in sweet cherry (Prunus avium). To identify the GST genes in sweet cherry, a genome-wide analysis was conducted. In this study, we identified 67 GST genes in P. avium genome and nomenclature according to chromosomal distribution. Phylogenetic tree analysis revealed that PavGST genes were classified into seven chief subfamily: TCHQD, Theta, Phi, Zeta, Lambda, DHAR, and Tau. The majority of the PavGST genes had a relatively well-maintained exon-intron and motif arrangement within the same group, according to gene structure and motif analyses. Gene structure (introns-exons) and conserved motif analysis revealed that the majority of the PavGST genes showed a relatively well-maintained motif and exons-introns configuration within the same group. The chromosomal localization, GO enrichment annotation, subcellular localization, syntenic relationship, Ka/Ks analysis, and molecular characteristics were accomplished using various bioinformatics tools. Mode of gene duplication showed that dispersed duplication might play a key role in the expansion of PavGST gene family. Promoter regions of PavGST genes contain numerous cis-regulatory components, which are involved in multiple stress responses, such as abiotic stress and phytohormones responsive factors. Furthermore, the expression profile of sweet cherry PavGSTs showed significant results under LED treatment. Our findings provide the groundwork for future research into induced LED anthocyanin and antioxidants deposition in sweet cherries.

Adaptation dynamics between copy-number and point mutations

Together, copy-number and point mutations form the basis for most evolutionary novelty, through the process of gene duplication and divergence. While a plethora of genomic data reveals the long-term fate of diverging coding sequences and their cis-regulatory elements, little is known about the early dynamics around the duplication event itself. In microorganisms, selection for increased gene expression often drives the expansion of gene copy-number mutations, which serves as a crude adaptation, prior to divergence through refining point mutations. Using a simple synthetic genetic reporter system that can distinguish between copy-number and point mutations, we study their early and transient adaptive dynamics in real time in Escherichia coli. We find two qualitatively different routes of adaptation, depending on the level of functional improvement needed. In conditions of high gene expression demand, the two mutation types occur as a combination. However, under low gene expression demand, copy-number and point mutations are mutually exclusive; here, owing to their higher frequency, adaptation is dominated by copy-number mutations, in a process we term amplification hindrance. Ultimately, due to high reversal rates and pleiotropic cost, copy-number mutations may not only serve as a crude and transient adaptation, but also constrain sequence divergence over evolutionary time scales.

KOMB: K-core based de novo characterization of copy number variation in microbiomes

Characterizing metagenomes via kmer-based, database-dependent taxonomic classification has yielded key insights into underlying microbiome dynamics. However, novel approaches are needed to track community dynamics and genomic flux within metagenomes, particularly in response to perturbations. We describe KOMB, a novel method for tracking genome level dynamics within microbiomes. KOMB utilizes K-core decomposition to identify Structural variations (SVs), specifically, population-level Copy Number Variation (CNV) within microbiomes. K-core decomposition partitions the graph into shells containing nodes of induced degree at least K, yielding reduced computational complexity compared to prior approaches. Through validation on a synthetic community, we show that KOMB recovers and profiles repetitive genomic regions in the sample. KOMB is shown to identify functionally-important regions in Human Microbiome Project datasets, and was used to analyze longitudinal data and identify keystone taxa in Fecal Microbiota Transplantation (FMT) samples. In summary, KOMB represents a novel graph-based, taxonomy-oblivious, and reference-free approach for tracking CNV within microbiomes. KOMB is open source and available for download at https://gitlab.com/treangenlab/komb.

Dosage sensitivity and exon shuffling shape the landscape of polymorphic duplicates in Drosophila and humans

Despite polymorphic duplicate genes' importance for the early stages of duplicate gene evolution, they are less studied than old gene duplicates. Two essential questions thus remain poorly addressed: how does dosage sensitivity, imposed by stoichiometry in protein complexes or by X chromosome dosage compensation, affect the emergence of complete duplicate genes? Do introns facilitate intergenic and intragenic chimaerism as predicted by the theory of exon shuffling? Here, we analysed new data for Drosophila and public data for humans, to characterize polymorphic duplicate genes with respect to dosage, exon-intron structures and allele frequencies. We found that complete duplicate genes are under dosage constraint induced by protein stoichiometry but potentially tolerated by X chromosome dosage compensation. We also found that in the intron-rich human genome, gene fusions and intragenic duplications extensively use intronic breakpoints generating in-frame proteins, in accordance with the theory of exon shuffling. Finally, we found that only a small proportion of complete or partial duplicates are at high frequencies, indicating the deleterious nature of dosage or gene structural changes. Altogether, we demonstrate how mechanistic factors including dosage sensitivity and exon-intron structure shape the short-term functional consequences of gene duplication.

Genome sequencing and transcriptomic analysis of the Andean killifish Orestias ascotanensis reveals adaptation to high-altitude aquatic life

Orestias ascotanensis (Cyprinodontidae) is a teleost pupfish endemic to springs feeding into the Ascotan saltpan in the Chilean Altiplano (3,700 m.a.s.l.) and represents an opportunity to study adaptations to high-altitude aquatic environments. We have de novo assembled the genome of O. ascotanensis at high coverage. Comparative analysis of the O. ascotanensis genome showed an overall process of contraction, including loss of genes related to G-protein signaling, chemotaxis and signal transduction, while there was expansion of gene families associated with microtubule-based movement and protein ubiquitination. We identified 818 genes under positive selection, many of which are involved in DNA repair. Additionally, we identified novel and conserved microRNAs expressed in O. ascotanensis and its closely-related species, Orestias gloriae. Our analysis suggests that positive selection and expansion of genes that preserve genome stability are a potential adaptive mechanism to cope with the increased solar UV radiation to which high-altitude animals are exposed to.

HSDFinder: A BLAST-Based Strategy for Identifying Highly Similar Duplicated Genes in Eukaryotic Genomes

Gene duplication is an important evolutionary mechanism capable of providing new genetic material for adaptive and nonadaptive evolution. However, bioinformatics tools for identifying duplicate genes are often limited to the detection of paralogs in multiple species or to specific types of gene duplicates, such as retrocopies. Here, we present a user-friendly, BLAST-based web tool, called HSDFinder, which can identify, annotate, categorize, and visualize highly similar duplicate genes (HSDs) in eukaryotic nuclear genomes. HSDFinder includes an online heatmap plotting option, allowing users to compare HSDs among different species and visualize the results in different Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway functional categories. The external software requirements are BLAST, InterProScan, and KEGG. The utility of HSDFinder was tested on various model eukaryotic species, including Chlamydomonas reinhardtii, Arabidopsis thaliana, Oryza sativa, and Zea mays as well as the psychrophilic green alga Chlamydomonas sp. UWO241, and was proven to be a practical and accurate tool for gene duplication analyses. The web tool is free to use at http://hsdfinder.com. Documentation and tutorials can be found via the GitHub: https://github.com/zx0223winner/HSDFinder.

Genome-Wide Identification, Characterization and Function Analysis of Lineage-Specific Genes in the Tea Plant Camellia sinensis

Genes that have no homologous sequences with other species are called lineage-specific genes (LSGs), are common in living organisms, and have an important role in the generation of new functions, adaptive evolution and phenotypic alteration of species. Camellia sinensis var. sinensis (CSS) is one of the most widely distributed cultivars for quality green tea production. The rich catechins in tea have antioxidant, free radical elimination, fat loss and cancer prevention potential. To further understand the evolution and utilize the function of LSGs in tea, we performed a comparative genomics approach to identify Camellia-specific genes (CSGs). Our result reveals that 1701 CSGs were identified specific to CSS, accounting for 3.37% of all protein-coding genes. The majority of CSGs (57.08%) were generated by gene duplication, and the time of duplication occurrence coincide with the time of two genome-wide replication (WGD) events that happened in CSS genome. Gene structure analysis revealed that CSGs have shorter gene lengths, fewer exons, higher GC content and higher isoelectric point. Gene expression analysis showed that CSG had more tissue-specific expression compared to evolutionary conserved genes (ECs). Weighted gene co-expression network analysis (WGCNA) showed that 18 CSGs are mainly associated with catechin synthesis-related pathways, including phenylalanine biosynthesis, biosynthesis of amino acids, pentose phosphate pathway, photosynthesis and carbon metabolism. Besides, we found that the expression of three CSGs (CSS0030246, CSS0002298, and CSS0030939) was significantly down-regulated in response to both types of stresses (salt and drought). Our study first systematically identified LSGs in CSS, and comprehensively analyzed the features and potential functions of CSGs. We also identified key candidate genes, which will provide valuable assistance for further studies on catechin synthesis and provide a molecular basis for the excavation of excellent germplasm resources.

Intraspecific Variation in the Rates of Mutations Causing Structural Variation in Daphnia magna

Mutations that cause structural variation are important sources of genetic variation upon which other evolutionary forces can act, however, they are difficult to observe and therefore few direct estimates of their rate and spectrum are available. Understanding mutation rate evolution, however, requires adding to the limited number of species for which direct estimates are available, quantifying levels of intraspecific variation in mutation rates, and assessing whether rate estimates co-vary across types of mutation. Here, we report structural variation-causing mutation rates (svcMRs) for six categories of mutations (short insertions and deletions, long deletions and duplications, and deletions and duplications at copy number variable sites) from nine genotypes of Daphnia magna collected from three populations in Finland, Germany, and Israel using a mutation accumulation approach. Based on whole-genome sequence data and validated using simulations, we find svcMRs are high (two orders of magnitude higher than base substitution mutation rates measured in the same lineages), highly variable among populations, and uncorrelated across categories of mutation. Furthermore, to assess the impact of scvMRs on the genome, we calculated rates while adjusting for the lengths of events and ran simulations to determine if the mutations occur in genic regions more or less frequently than expected by chance. Our results pose a challenge to most prevailing theories aimed at explaining the evolution of the mutation rate, underscoring the importance of obtaining additional mutation rate estimates in more genotypes, for more types of mutation, in more species, in order to improve our future understanding of mutation rates, their variation, and their evolution.

Cyclic nucleotide gated channel genes (CNGCs) in Rosaceae: genome-wide annotation, evolution and the roles on Valsa canker resistance

In Rosaceae, tandem duplication caused the drastic expansion of CNGC gene family Group I. The members MdCN11 and MdCN19 negatively regulate Valsa canker resistance. Apple (Malus domestica) and pear (Pyrus bretschneideri and P. communis) are important fruit crops in Rosaceae family but are suffering from threats of Valsa canker. Cyclic nucleotide-gated ion channels (CNGCs) take crucial roles in plant immune responses. In the present study, a total of 355 CNGCs was identified from 8 Rosaceae plants. Based on phylogenetic analysis, 540 CNGCs from 18 plants (8 in Rosaceae and 10 others) could be divided into four groups. Group I was greatly expanded in Rosaceae resulted from tandem duplications. A large number of cis-acting regulatory elements (cis-elements) responsive to signals from multiple stresses and hormones were identified in the promoter regions of CNGCs in Malus spp. and Pyrus spp. Expressions of most Group I members were obviously up-regulated in Valsa canker susceptible varieties but not in the resistant ones. Furthermore, overexpression of the MdCN11 and MdCN19 in both apple fruits and 'Duli' (P. betulifolia) suspension cells compromised Valsa canker resistance. Overexpression of MdCN11 induced expression of hypersensitive response (HR)-related genes. In conclusion, tandem duplication resulted in a drastic expansion of CNGC Group I members in Rosaceae. Among these, MdCN11 and MdCN19 negatively regulate the Valsa canker resistance via inducting HR.

Wheat in vivo RNA structure landscape reveals a prevalent role of RNA structure in modulating translational subgenome expression asymmetry

BACKGROUND

Polyploidy, especially allopolyploidy, which entails merging divergent genomes via hybridization and whole-genome duplication (WGD), is a major route to speciation in plants. The duplication among the parental genomes (subgenomes) often leads to one subgenome becoming dominant over the other(s), resulting in subgenome asymmetry in gene content and expression. Polyploid wheats are allopolyploids with most genes present in two (tetraploid) or three (hexaploid) functional copies, which commonly show subgenome expression asymmetry. It is unknown whether a similar subgenome asymmetry exists during translation. We aim to address this key biological question and explore the major contributing factors to subgenome translation asymmetry.

RESULTS

Here, we obtain the first tetraploid wheat translatome and reveal that subgenome expression asymmetry exists at the translational level. We further perform in vivo RNA structure profiling to obtain the wheat RNA structure landscape and find that mRNA structure has a strong impact on translation, independent of GC content. We discover a previously uncharacterized contribution of RNA structure in subgenome translation asymmetry. We identify 3564 single-nucleotide variations (SNVs) across the transcriptomes between the two tetraploid wheat subgenomes, which induce large RNA structure disparities. These SNVs are highly conserved within durum wheat cultivars but are divergent in both domesticated and wild emmer wheat.

CONCLUSIONS

We successfully determine both the translatome and in vivo RNA structurome in tetraploid wheat. We reveal that RNA structure serves as an important modulator of translational subgenome expression asymmetry in polyploids. Our work provides a new perspective for molecular breeding of major polyploid crops.

Genome-wide investigation and comparative analysis of MATE gene family in Rosaceae species and their regulatory role in abiotic stress responses in Chinese pear (Pyrus bretschneideri)

The Multidrug and Toxic Compound Extrusion (MATE) protein belongs to a secondary transporter gene family, which plays a primary role in transporting many kinds of substrates such as organic compounds, secondary metabolites, and phytohormones. MATE protein members exist in both prokaryotes and eukaryotes. However, evolution and comprehensive analysis of the MATE genes has not been performed in Rosaceae species. In the present study, a total of 404 MATEs genes were identified from six Rosaceae genomes (Prunus avium, Pyrus bretschneideri, Prunus persica, Fragaria vesca, Prunus mume, and Malus domestica) and classified into eight main subfamilies (I-VII) based on structural and phylogenetic analysis. Microcollinearity analysis showed that whole-genome duplication events might play a vital role in the expansion of the MATE genes family. The Ka/Ks analysis, chromosomal localization, subcellular localization, and molecular characteristics (length, weight, and pI) were performed using various bioinformatics tools. Furthermore, different subfamilies have different introns-exons structures, cis-acting elements, and conserved motifs analysis, indicating functional divergence in the MATE family. Subsequently, RNA-seq analysis and real-time qRT-PCR were conducted during Chinese pear fruit development. Moreover, PbMATE genes were significantly expressed under hormonal treatments of MeJA (methyl jasmonate), SA (salicylic acid), and ABA (abscisic acid). Overall, our results provide helpful insights into the functions, expansion complexity, and evolutions of the MATE genes in Chinese pear and five Rosaceae species.

Protocol for HSDFinder: Identifying, annotating, categorizing, and visualizing duplicated genes in eukaryotic genomes

Although gene duplications have been documented in many species, the precise numbers of highly similar duplicated genes (HSDs) in eukaryotic nuclear genomes remain largely unknown and can be time-consuming to explore. We developed HSDFinder to identify, categorize, and visualize HSDs in eukaryotic nuclear genomes using protein family domains and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. In contrast to existing tools, HSDFinder allows users to compare HSDs among different species and visualize results in different KEGG pathway functional categories via heatmap plotting. For complete details on the use and execution of this protocol, please refer to Zhang et al. (2021).

Bridging Tumorigenesis and Therapy Resistance With a Non-Darwinian and Non-Lamarckian Mechanism of Adaptive Evolution

Although neo-Darwinian (and less often Lamarckian) dynamics are regularly invoked to interpret cancer's multifarious molecular profiles, they shine little light on how tumorigenesis unfolds and often fail to fully capture the frequency and breadth of resistance mechanisms. This uncertainty frames one of the most problematic gaps between science and practice in modern times. Here, we offer a theory of adaptive cancer evolution, which builds on a molecular mechanism that lies outside neo-Darwinian and Lamarckian schemes. This mechanism coherently integrates non-genetic and genetic changes, ecological and evolutionary time scales, and shifts the spotlight away from positive selection towards purifying selection, genetic drift, and the creative-disruptive power of environmental change. The surprisingly simple use-it or lose-it rationale of the proposed theory can help predict molecular dynamics during tumorigenesis. It also provides simple rules of thumb that should help improve therapeutic approaches in cancer.

Chromosome-scale assembly and evolution of the tetraploid Salvia splendens (Lamiaceae) genome

Polyploidization plays a key role in plant evolution, but the forces driving the fate of homoeologs in polyploid genomes, i.e., paralogs resulting from a whole-genome duplication (WGD) event, remain to be elucidated. Here, we present a chromosome-scale genome assembly of tetraploid scarlet sage (Salvia splendens), one of the most diverse ornamental plants. We found evidence for three WGD events following an older WGD event shared by most eudicots (the γ event). A comprehensive, spatiotemporal, genome-wide analysis of homoeologs from the most recent WGD unveiled expression asymmetries, which could be associated with genomic rearrangements, transposable element proximity discrepancies, coding sequence variation, selection pressure, and transcription factor binding site differences. The observed differences between homoeologs may reflect the first step toward sub- and/or neofunctionalization. This assembly provides a powerful tool for understanding WGD and gene and genome evolution and is useful in developing functional genomics and genetic engineering strategies for scarlet sage and other Lamiaceae species.

Telomere Roles in Fungal Genome Evolution and Adaptation

Telomeres form the ends of linear chromosomes and usually comprise protein complexes that bind to simple repeated sequence motifs that are added to the 3′ ends of DNA by the telomerase reverse transcriptase (TERT). One of the primary functions attributed to telomeres is to solve the “end-replication problem” which, if left unaddressed, would cause gradual, inexorable attrition of sequences from the chromosome ends and, eventually, loss of viability. Telomere-binding proteins also protect the chromosome from 5′ to 3′ exonuclease action, and disguise the chromosome ends from the double-strand break repair machinery whose illegitimate action potentially generates catastrophic chromosome aberrations. Telomeres are of special interest in the blast fungus, Pyricularia, because the adjacent regions are enriched in genes controlling interactions with host plants, and the chromosome ends show enhanced polymorphism and genetic instability. Previously, we showed that telomere instability in some P. oryzae strains is caused by novel retrotransposons (MoTeRs) that insert in telomere repeats, generating interstitial telomere sequences that drive frequent, break-induced rearrangements. Here, we sought to gain further insight on telomeric involvement in shaping Pyricularia genome architecture by characterizing sequence polymorphisms at chromosome ends, and surrounding internalized MoTeR loci (relics) and interstitial telomere repeats. This provided evidence that telomere dynamics have played historical, and likely ongoing, roles in shaping the Pyricularia genome. We further demonstrate that even telomeres lacking MoTeR insertions are poorly preserved, such that the telomere-adjacent sequences exhibit frequent presence/absence polymorphism, as well as exchanges with the genome interior. Using TERT knockout experiments, we characterized chromosomal responses to failed telomere maintenance which suggested that much of the MoTeR relic-/interstitial telomere-associated polymorphism could be driven by compromised telomere function. Finally, we describe three possible examples of a phenomenon known as “Adaptive Telomere Failure,” where spontaneous losses of telomere maintenance drive rapid accumulation of sequence polymorphism with possible adaptive advantages. Together, our data suggest that telomere maintenance is frequently compromised in Pyricularia but the chromosome alterations resulting from telomere failure are not as catastrophic as prior research would predict, and may, in fact, be potent drivers of adaptive polymorphism.

Natural variation in the consequences of gene overexpression and its implications for evolutionary trajectories

Copy number variation through gene or chromosome amplification provides a route for rapid phenotypic variation and supports the long-term evolution of gene functions. Although the evolutionary importance of copy-number variation is known, little is understood about how genetic background influences its tolerance. Here, we measured fitness costs of over 4000 overexpressed genes in 15 Saccharomyces cerevisiae strains representing different lineages, to explore natural variation in tolerating gene overexpression (OE). Strain-specific effects dominated the fitness costs of gene OE. We report global differences in the consequences of gene OE, independent of the amplified gene, as well as gene-specific effects that were dependent on the genetic background. Natural variation in the response to gene OE could be explained by several models, including strain-specific physiological differences, resource limitations, and regulatory sensitivities. This work provides new insight on how genetic background influences tolerance to gene amplification and the evolutionary trajectories accessible to different backgrounds.

The Evolutionary Dynamics of Genetic Incompatibilities Introduced by Duplicated Genes in Arabidopsis thaliana

Although gene duplications provide genetic backup and allow genomic changes under relaxed selection, they may potentially limit gene flow. When different copies of a duplicated gene are pseudofunctionalized in different genotypes, genetic incompatibilities can arise in their hybrid offspring. Although such cases have been reported after manual crosses, it remains unclear whether they occur in nature and how they affect natural populations. Here, we identified four duplicated-gene based incompatibilities including one previously not reported within an artificial Arabidopsis intercross population. Unexpectedly, however, for each of the genetic incompatibilities we also identified the incompatible alleles in natural populations based on the genomes of 1,135 Arabidopsis accessions published by the 1001 Genomes Project. Using the presence of incompatible allele combinations as phenotypes for GWAS, we mapped genomic regions that included additional gene copies which likely rescue the genetic incompatibility. Reconstructing the geographic origins and evolutionary trajectories of the individual alleles suggested that incompatible alleles frequently coexist, even in geographically closed regions, and that their effects can be overcome by additional gene copies collectively shaping the evolutionary dynamics of duplicated genes during population history.

The Role of Gene Duplication in the Divergence of Enzyme Function: A Comparative Approach

Gene duplication is a crucial process involved in the appearance of new genes and functions. It is thought to have played a major role in the growth of enzyme families and the expansion of metabolism at the biosphere's dawn and in recent times. Here, we analyzed paralogous enzyme content within each of the seven enzymatic classes for a representative sample of prokaryotes by a comparative approach. We found a high ratio of paralogs for three enzymatic classes: oxidoreductases, isomerases, and translocases, and within each of them, most of the paralogs belong to only a few subclasses. Our results suggest an intricate scenario for the evolution of prokaryotic enzymes, involving different fates for duplicated enzymes fixed in the genome, where around 20-40% of prokaryotic enzymes have paralogs. Intracellular organisms have a lesser ratio of duplicated enzymes, whereas free-living enzymes show the highest ratios. We also found that phylogenetically close phyla and some unrelated but with the same lifestyle share similar genomic and biochemical traits, which ultimately support the idea that gene duplication is associated with environmental adaptation.

High Complexity and Degree of Genetic Variation in Brettanomyces bruxellensis Population

Genome-wide characterization of genetic variants of a large population of individuals within the same species is essential to have a deeper insight into its evolutionary history as well as the genotype–phenotype relationship. Population genomic surveys have been performed in multiple yeast species, including the two model organisms, Saccharomyces cerevisiae and Schizosaccharomyces pombe. In this context, we sought to characterize at the population level the Brettanomyces bruxellensis yeast species, which is a major cause of wine spoilage and can contribute to the specific flavor profile of some Belgium beers. We have completely sequenced the genome of 53 B. bruxellensis strains isolated worldwide. The annotation of the reference genome allowed us to define the gene content of this species. As previously suggested, our genomic data clearly highlighted that genetic diversity variation is related to ploidy level, which is variable in the B. bruxellensis species. Genomes are punctuated by multiple loss-of-heterozygosity regions, whereas aneuploidies as well as segmental duplications are uncommon. Interestingly, triploid genomes are more prone to gene copy number variation than diploids. Finally, the pangenome of the species was reconstructed and was found to be small with few accessory genes compared with S. cerevisiae. The pangenome is composed of 5,409 ORFs (open reading frames) among which 5,106 core ORFs and 303 ORFs that are variable within the population. All these results highlight the different trajectories of species evolution and consequently the interest of establishing population genomic surveys in more populations.

Exact Probability Landscapes of Stochastic Phenotype Switching in Feed-Forward Loops: Phase Diagrams of Multimodality

Feed-forward loops (FFLs) are among the most ubiquitously found motifs of reaction networks in nature. However, little is known about their stochastic behavior and the variety of network phenotypes they can exhibit. In this study, we provide full characterizations of the properties of stochastic multimodality of FFLs, and how switching between different network phenotypes are controlled. We have computed the exact steady-state probability landscapes of all eight types of coherent and incoherent FFLs using the finite-butter Accurate Chemical Master Equation (ACME) algorithm, and quantified the exact topological features of their high-dimensional probability landscapes using persistent homology. Through analysis of the degree of multimodality for each of a set of 10,812 probability landscapes, where each landscape resides over 10⁵–10⁶ microstates, we have constructed comprehensive phase diagrams of all relevant behavior of FFL multimodality over broad ranges of input and regulation intensities, as well as different regimes of promoter binding dynamics. In addition, we have quantified the topological sensitivity of the multimodality of the landscapes to regulation intensities. Our results show that with slow binding and unbinding dynamics of transcription factor to promoter, FFLs exhibit strong stochastic behavior that is very different from what would be inferred from deterministic models. In addition, input intensity play major roles in the phenotypes of FFLs: At weak input intensity, FFL exhibit monomodality, but strong input intensity may result in up to 6 stable phenotypes. Furthermore, we found that gene duplication can enlarge stable regions of specific multimodalities and enrich the phenotypic diversity of FFL networks, providing means for cells toward better adaptation to changing environment. Our results are directly applicable to analysis of behavior of FFLs in biological processes such as stem cell differentiation and for design of synthetic networks when certain phenotypic behavior is desired.

BLAST from the Past: Impacts of Evolving Approaches on Studies of Evolution by Gene Duplication

In 1970, Susumu Ohno hypothesized that gene duplication was a major reservoir of adaptive innovation. However, it was not until over two decades later that DNA sequencing studies uncovered the ubiquity of gene duplication across all domains of life, highlighting its global importance in the evolution of phenotypic complexity and species diversification. Today, it seems that there are no limits to the study of evolution by gene duplication, as it has rapidly coevolved with numerous experimental and computational advances in genomics. In this perspective, we examine word stem usage in PubMed abstracts to infer how evolving discoveries and technologies have shaped the landscape of studying evolution by gene duplication, leading to a more refined understanding of its role in the emergence of novel phenotypes.

Genomic data support management of anadromous Arctic Char fisheries in Nunavik by highlighting neutral and putatively adaptive genetic variation

Distinguishing neutral and adaptive genetic variation is one of the main challenges in investigating processes shaping population structure in the wild, and landscape genomics can help identify signatures of adaptation to contrasting environments. Arctic Char (Salvelinus alpinus) is an anadromous salmonid and the most harvested fish species by Inuit people, including in Nunavik (Québec, Canada), one of the most recently deglaciated regions in the world. Unlike many other anadromous salmonids, Arctic Char occupy coastal habitats near their natal rivers during their short marine phase restricted to the summer ice-free period. Our main objective was to document putatively neutral and adaptive genomic variation in anadromous Arctic Char populations from Nunavik and bordering regions to inform local fisheries management. We used genotyping by sequencing (GBS) to genotype 18,112 filtered single nucleotide polymorphisms (SNP) in 650 individuals from 23 sampling locations along >2000 km of coastline. Our results reveal a hierarchical genetic structure, whereby neighboring hydrographic systems harbor distinct populations grouped by major oceanographic basins: Hudson Bay, Hudson Strait, Ungava Bay, and Labrador Sea. We found genetic diversity and differentiation to be consistent both with the expected postglacial recolonization history and with patterns of isolation-by-distance reflecting contemporary gene flow. Results from three gene-environment association methods supported the hypothesis of local adaptation to both freshwater and marine environments (strongest associations with sea surface and air temperatures during summer and salinity). Our results support a fisheries management strategy at a regional scale, and other implications for hatchery projects and adaptation to climate change are discussed.

Expanded Diversity and Phylogeny of mer Genes Broadens Mercury Resistance Paradigms and Reveals an Origin for MerA Among Thermophilic Archaea

Mercury (Hg) is a highly toxic element due to its high affinity for protein sulfhydryl groups, which upon binding, can destabilize protein structure and decrease enzyme activity. Prokaryotes have evolved enzymatic mechanisms to detoxify inorganic Hg and organic Hg (e.g., MeHg) through the activities of mercuric reductase (MerA) and organomercury lyase (MerB), respectively. Here, the taxonomic distribution and evolution of MerAB was examined in 84,032 archaeal and bacterial genomes, metagenome assembled genomes, and single-cell genomes. Homologs of MerA and MerB were identified in 7.8 and 2.1% percent of genomes, respectively. MerA was identified in the genomes of 10 archaeal and 28 bacterial phyla previously unknown to code for this functionality. Likewise, MerB was identified in 2 archaeal and 11 bacterial phyla previously unknown to encode this functionality. Surprisingly, homologs of MerB were identified in a number of genomes (∼50% of all MerB-encoding genomes) that did not encode MerA, suggesting alternative mechanisms to detoxify Hg(II) once it is generated in the cytoplasm. Phylogenetic reconstruction of MerA place its origin in thermophilic Thermoprotei (Crenarchaeota), consistent with high levels of Hg(II) in geothermal environments, the natural habitat of this archaeal class. MerB appears to have been recruited to the mer operon relatively recently and likely among a mesophilic ancestor of Euryarchaeota and Thaumarchaeota. This is consistent with the functional dependence of MerB on MerA and the widespread distribution of mesophilic microorganisms that methylate Hg(II) at lower temperature. Collectively, these results expand the taxonomic and ecological distribution of mer-encoded functionalities, and suggest that selection for Hg(II) and MeHg detoxification is dependent not only on the availability and type of mercury compounds in the environment but also the physiological potential of the microbes who inhabit these environments. The expanded diversity and environmental distribution of MerAB identify new targets to prioritize for future research.

Adaptive evolution of sulfite tolerance in Brettanomyces bruxellensis

Brettanomyces bruxellensis is considered one of the most problematic microbes associated with wine production. Sulfur dioxide is commonly used to inhibit the growth of B. bruxellensis and limit the potential wine spoilage. Brettanomyces bruxellensis wine isolates can grow at higher concentrations of this preservative than isolates from other sources. Thus, it has been suggested that the use of sulfite may have selected for B. bruxellensis strains better adapted to survive in the winemaking environment. We utilized laboratory adaptive evolution to determine the potential for this to occur. Three B. bruxellensis strains, representative of known genetic variation within the species, were subjected to increasing sublethal sulfur dioxide concentrations. Individual clones isolated from evolved populations displayed enhanced sulfite tolerance, ranging from 1.6 to 2.5 times higher than the corresponding parental strains. Whole-genome sequencing of sulfite-tolerant clones derived from two of the parental strains revealed structural variations affecting 270 genes. The region containing the sulfite efflux pump encoding gene, SSU1, showed clear copy number variants in all sequenced clones. Regardless of parental strain genetic background, SSU1 copy number changes were reproducibly associated with one SSU1 haplotype. This work clearly demonstrates adaptive evolution of B. bruxellensis when exposed to sublethal sulfites and suggests that, similar to Saccharomyces cerevisiae wine yeast, the mechanism responsible involves the gene SSU1.

Whole-Genome Doubling Affects Pre-miRNA Expression in Plants

Whole-genome doubling (polyploidy) is common in angiosperms. Several studies have indicated that it is often associated with molecular, physiological, and phenotypic changes. Mounting evidence has pointed out that micro-RNAs (miRNAs) may have an important role in whole-genome doubling. However, an integrative approach that compares miRNA expression in polyploids is still lacking. Here, a re-analysis of already published RNAseq datasets was performed to identify microRNAs’ precursors (pre-miRNAs) in diploids (2x) and tetraploids (4x) of five species (Arabidopsis thaliana L., Morus alba L., Brassica rapa L., Isatis indigotica Fort., and Solanum commersonii Dun). We found 3568 pre-miRNAs, three of which (pre-miR414, pre-miR5538, and pre-miR5141) were abundant in all 2x, and were absent/low in their 4x counterparts. They are predicted to target more than one mRNA transcript, many belonging to transcription factors (TFs), DNA repair mechanisms, and related to stress. Sixteen pre-miRNAs were found in common in all 2x and 4x. Among them, pre-miRNA482, pre-miRNA2916, and pre-miRNA167 changed their expression after polyploidization, being induced or repressed in 4x plants. Based on our results, a common ploidy-dependent response was triggered in all species under investigation, which involves DNA repair, ATP-synthesis, terpenoid biosynthesis, and several stress-responsive transcripts. In addition, an ad hoc pre-miRNA expression analysis carried out solely on 2x vs. 4x samples of S. commersonii indicated that ploidy-dependent pre-miRNAs seem to actively regulate the nucleotide metabolism, probably to cope with the increased requirement for DNA building blocks caused by the augmented DNA content. Overall, the results outline the critical role of microRNA-mediated responses following autopolyploidization in plants.

The evolution of opsin genes in five species of mirid bugs: duplication of long-wavelength opsins and loss of blue-sensitive opsins

BACKGROUND

Color vision and phototactic behavior based on opsins are important for the fitness of insects because of their roles in foraging and mate choice. Related topics, including the duplication and loss of opsin genes, have been well investigated in insect orders such as Coleoptera, Lepidoptera, Hymenoptera, Odonata and Orthoptera, and the findings have been used to develop pest management strategies involving light trapping. Mirid bugs of Hemiptera, which are pests that cause heavy economic losses, show capacity for color discrimination and phototaxis. However, the opsins in mirid bugs remain uncharacterized. Herein, we examined five species to investigate the evolution of opsins in the family Miridae.

RESULTS

Using RNA-seq, we identified several contigs showing high identity with opsins, including four contigs in Apolygus lucorum and three contigs each in Adelphocoris suturalis, Adelphocoris fasciaticollis, Adelphocoris lineolatus and Nesidiocoris tenuis. Phylogenetic analyses indicated that one of these genes clustered with ultraviolet-sensitive (UV) opsins and that the others clustered with long-wavelength (LW) opsins, suggesting that duplication of LW opsins and loss of blue light-sensitive (B) opsins occurred in mirid bugs. The existence of introns in the LW opsins of mirid bugs suggested that the duplication events were DNA based. Both LW1 and LW2 opsins of mirid bugs were found to be under strong purifying selection. The LW1 opsins were significantly more highly expressed than the LW2 and UV opsins.

CONCLUSIONS

We identified the opsins of mirid bugs using five selected mirid species as a representative sample. Phylogenetic analyses clustered one of the genes with UV opsins and the others with LW opsins, suggesting the occurrence of LW opsin duplication and B opsin loss during the evolution of mirid bugs. Intron detection suggested that the identified duplication event was DNA based. The evidence of strong purifying selection and the relatively high expression levels suggested that these opsins exhibit fundamental functions in mirid bugs.

Ion Transporters and Osmoregulation in the Kidney of Teleost Fishes as a Function of Salinity

Euryhaline teleosts exhibit major changes in renal function as they move between freshwater (FW) and seawater (SW) environments, thus tolerating large fluctuations in salinity. In FW, the kidney excretes large volumes of water through high glomerular filtration rates (GFR) and low tubular reabsorption rates, while actively reabsorbing most ions at high rates. The excreted product has a high urine flow rate (UFR) with a dilute composition. In SW, GFR is greatly reduced, and the tubules reabsorb as much water as possible, while actively secreting divalent ions. The excreted product has a low UFR, and is almost isosmotic to the blood plasma, with Mg²⁺, SO₄^2–, and Cl^– as the major ionic components. Early studies at the organismal level have described these basic patterns, while in the last two decades, studies of regulation at the cell and molecular level have been implemented, though only in a few euryhaline groups (salmonids, eels, tilapias, and fugus). There have been few studies combining the two approaches. The aim of the review is to integrate known aspects of renal physiology (reabsorption and secretion) with more recent advances in molecular water and solute physiology (gene and protein function of transporters). The renal transporters addressed include the subunits of the Na⁺, K⁺- ATPase (NKA) enzyme, monovalent ion transporters for Na⁺, Cl^–, and K⁺ (NKCC1, NKCC2, CLC-K, NCC, ROMK2), water transport pathways [aquaporins (AQP), claudins (CLDN)], and divalent ion transporters for SO₄^2–, Mg²⁺, and Ca²⁺ (SLC26A6, SLC26A1, SLC13A1, SLC41A1, CNNM2, CNNM3, NCX1, NCX2, PMCA). For each transport category, we address the current understanding at the molecular level, try to synthesize it with classical knowledge of overall renal function, and highlight knowledge gaps. Future research on the kidney of euryhaline fishes should focus on integrating changes in kidney reabsorption and secretion of ions with changes in transporter function at the cellular and molecular level (gene and protein verification) in different regions of the nephrons. An increased focus on the kidney individually and its functional integration with the other osmoregulatory organs (gills, skin and intestine) in maintaining overall homeostasis will have applied relevance for aquaculture.

Homoeolog Inference Methods Requiring Bidirectional Best Hits or Synteny Miss Many Pairs

Homoeologs are pairs of genes or chromosomes in the same species that originated by speciation and were brought back together in the same genome by allopolyploidization. Bioinformatic methods for accurate homoeology inference are crucial for studying the evolutionary consequences of polyploidization, and homoeology is typically inferred on the basis of bidirectional best hit (BBH) and/or positional conservation (synteny). However, these methods neglect the fact that genes can duplicate and move, both prior to and after the allopolyploidization event. These duplications and movements can result in many-to-many and/or nonsyntenic homoeologs-which thus remain undetected and unstudied. Here, using the allotetraploid upland cotton (Gossypium hirsutum) as a case study, we show that conventional approaches indeed miss a substantial proportion of homoeologs. Additionally, we found that many of the missed pairs of homoeologs are broadly and highly expressed. A gene ontology analysis revealed a high proportion of the nonsyntenic and non-BBH homoeologs to be involved in protein translation and are likely to contribute to the functional repertoire of cotton. Thus, from an evolutionary and functional genomics standpoint, choosing a homoeolog inference method which does not solely rely on 1:1 relationship cardinality or synteny is crucial for not missing these potentially important homoeolog pairs.