Search

Supergenes offer spectacular examples of long-term balancing selection in nature, but their origin and maintenance remain a mystery. Reduced recombination between arrangements, a critical aspect of many supergenes, protects adaptive multi-trait phenotypes but can lead to mutation accumulation. Mutation accumulation can stabilize the system through the emergence of associative overdominance (AOD), destabilize the system, or lead to new evolutionary outcomes. One outcome is the formation of maladaptive balanced lethal systems, where only heterozygotes remain viable and reproduce. We investigated the conditions under which these different outcomes occur, assuming a scenario of introgression after divergence. We found that AOD aided the invasion of a new supergene arrangement and the establishment of a polymorphism. However, this polymorphism was easily destabilized by further mutation accumulation, which was often asymmetric, disrupting the quasi-equilibrium state. Mechanisms that accelerated degeneration tended to amplify asymmetric mutation accumulation between the supergene arrangements and vice-versa. As the evolution of balanced lethal systems requires symmetric degeneration of both arrangements, this leaves only restricted conditions for their evolution, namely small population sizes and low rates of gene conversion. The dichotomy between the persistence of polymorphism and degeneration of supergene arrangements likely underlies the rarity of balanced lethal systems in nature. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'. [ABSTRACT FROM AUTHOR]

Simple Summary: With the increasing popularity of omics technologies, many scientists are dealing with large sets of genomic data for the first time. While there are many amazing bioinformatics tools available to help analyze this data, one common difficulty these scientists face is that not all gene data or analysis tools use the same genomic nomenclature. Another issue is that many publications still use obsolete or colloquial gene aliases instead of official gene identifiers. Common gene ID conversion tools and other downstream analysis software struggle with these gene aliases. Therefore, we developed a free and publicly available web application, GeneToList, to assist in gene ID disambiguation and gene ID conversion, with a specific focus on a user-friendly interface for the non-bioinformatics-savvy scientist. The increasing incorporation of omics technologies into biomedical research and translational medicine presents challenges to end users of the large and complex datasets that are generated by these methods. A particular challenge in genomics is that the nomenclature for genes is not uniform between large genomic databases or between commonly used genetic analysis tools. Furthermore, outdated genomic nomenclature can still be found amongst scientific communications, including peer-reviewed manuscripts. Therefore, a web application (GeneToList) was developed to assist in gene ID conversion and alias matching, with a specific focus on achieving a user-friendly interface for the non-bioinformatics-savvy scientist. It currently includes gene information for over 38,000 different taxa retrieved from the National Center for Biotechnology and Information (NCBI) Gene resource. Supported databases of gene IDs include NCBI Gene Symbols, NCBI Gene IDs (Entrez IDs), OMIM IDs, HGNC IDs, Ensembl IDs, and 28 other taxa-specific identifiers. GeneToList is available at genetolist.com. The tool is a web application that is compatible with many standard browsers. The gene ID conversion feature of this application was found to outcompete the common gene ID conversion tools. Specifically, it was able to successfully convert all tested IDs, whereas the others were not able to recognize the gene aliases. Therefore, the gene ID disambiguation provided by this application should be beneficial for many scientists dealing with gene data when the uniformity of gene nomenclature is important for downstream analysis. [ABSTRACT FROM AUTHOR]

Protein coding genes terminate with one of three stop codons (TAA, TGA, or TAG) that, like synonymous codons, are not employed equally. With TGA and TAG having identical nucleotide content, analysis of their differential usage provides an unusual window into the forces operating on what are ostensibly functionally identical residues. Across genomes and between isochores within the human genome, TGA usage increases with G + C content but, with a common G + C → A + T mutation bias, this cannot be explained by mutation bias-drift equilibrium. Increased usage of TGA in G + C-rich genomes or genomic regions is also unlikely to reflect selection for the optimal stop codon, as TAA appears to be universally optimal, probably because it has the lowest read-through rate. Despite TAA being favored by selection and mutation bias, as with codon usage bias G + C pressure is the prime determinant of between-species TGA usage trends. In species with strong G + C-biased gene conversion (gBGC), such as mammals and birds, the high usage and conservation of TGA is best explained by an A + T → G + C repair bias. How to explain TGA enrichment in other G + C-rich genomes is less clear. Enigmatically, across bacterial and archaeal species and between human isochores TAG usage is mostly unresponsive to G + C pressure. This unresponsiveness we dub the TAG paradox as currently no mutational, selective, or gBGC model provides a well-supported explanation. That TAG does increase with G + C usage across eukaryotes makes the usage elsewhere yet more enigmatic. We suggest resolution of the TAG paradox may provide insights into either an unknown but common selective preference (probably at the DNA/RNA level) or an unrecognized complexity to the action of gBGC. [ABSTRACT FROM AUTHOR]

Apomictic genera have a complex evolutionary history, including reticulate hybridization, polyploidization and variation in reproduction modes. Restrictions of functional meiosis in polyploid apomictic taxa considerably hamper gene conversion, leaving footprints of past hybridization events in their genomes, which are masked by the homogenization of tandemly organized arrays of nuclear ribosomal DNA in sexual species. Dandelions (Taraxacum) have a highly complex reticulate evolutionary history. Detailed knowledge remains hazy; earlier investigations of rDNA sequences have uncovered the complex pattern derived from evolution but without direct evidence for any particular processes. We investigated the position and number of 45S and 5S rDNA loci in 38 Taraxacum taxa (covering different reproduction modes, geographical regions and putative phylogenetic groups) using fluorescent in situ hybridization (FISH) and measured genome size and GC content. The ITS1-5.8S-ITS2 regions of four sexual and five apomictic taxa were sequenced to investigate inter- and intra-individual variation. Most species considerably differ in the chromosome positions of loci and karyotype patterns, but mostly share the same number of studied loci (45S:5S) in a 1:2 ratio per haploid genome (x = 8), with six exceptions (up to a 4:2 ratio). Genome size (2C) varies sixfold and with GC content partly distinguishes major evolutionary groups. Sexual taxa show limited variation in sequenced regions (with two to eight variants), but apomictic taxa vary significantly (with 20–36 variants). Extensive reticulate evolution in Taraxacum and subsequent phenomena such as genome repatterning and non-effective concerted evolution are probably the cause of the dynamic nature of Taraxacum karyotypes and the large variation in genome size and rDNA sequences. [ABSTRACT FROM AUTHOR]

Progression through meiosis is associated with significant reorganization of chromosome structure, regulated in part by changes in histones and chromatin. Prior studies observed defects in meiotic progression in yeast strains lacking the linker histone H1 or variant histone H2A.Z. To further define the contributions of these chromatin factors, we have conducted genetic and cytological analysis of cells undergoing meiosis in the absence of H1 and H2A.Z. We find that a spore viability defect observed in strains lacking H2A.Z can be partially suppressed if cells also lack histone H1, while the combined loss of both H1 and H2A.Z is associated with elevated gene conversion events. Cytological analysis of Red1 and Rec8 staining patterns indicates that a subset of cells lacking H2A.Z fail to assemble a proper chromosome axis, and the staining pattern of the synaptonemal complex protein Zip1 in htz1D/htz1D cells mimics that of cells deficient for Rec8-dependent meiotic cohesion. Our results suggest a role for H2A.Z in the establishment or maintenance of the meiotic chromosome axis, possibly by promoting the efficient chromosome cohesion. [ABSTRACT FROM AUTHOR]

Finally, residual expression of LAMB3 protein in perilesional and uninvolved patient skin, as detected by western blot, further supports the presence of RM in these clinically restored skin areas (Fig. Common molecular mechanisms of phenotypic alleviation in genetic diseases include loss of heterozygosity (LOH), gene deletion, nondisjunction mutations, as well as nonsense codon readthrough and RNA I trans i -splicing. All authors declare to have no conflicts of interest or affiliations with organizations or entities with a financial interest in or financial conflict with the subject matter or materials discussed in the manuscript and that this information is accurate, complete and up-to-date. [Extracted from the article]

Introgressions of chromosomal segments from related species into wheat are important sources of resistance against fungal diseases. The durability and effectiveness of introgressed resistance genes upon agricultural deployment is highly variable—a phenomenon that remains poorly understood, as the corresponding fungal avirulence genes are largely unknown. Until its breakdown, the Pm17 resistance gene introgressed from rye to wheat provided broad resistance against powdery mildew (Blumeria graminis). Here, we used quantitative trait locus (QTL) mapping to identify the corresponding wheat mildew avirulence effector AvrPm17. It is encoded by two paralogous genes that exhibit signatures of reoccurring gene conversion events and are members of a mildew sublineage specific effector cluster. Extensive haplovariant mining in wheat mildew and related sublineages identified several ancient virulent AvrPm17 variants that were present as standing genetic variation in wheat powdery mildew prior to the Pm17 introgression, thereby paving the way for the rapid breakdown of the Pm17 resistance. QTL mapping in mildew identified a second genetic component likely corresponding to an additional resistance gene present on the 1AL.1RS translocation carrying Pm17. This gene remained previously undetected due to suppressed recombination within the introgressed rye chromosomal segment. We conclude that the initial effectiveness of 1AL.1RS was based on simultaneous introgression of two genetically linked resistance genes. Our results demonstrate the relevance of pathogen-based genetic approaches to disentangling complex resistance loci in wheat. We propose that identification and monitoring of avirulence gene diversity in pathogen populations become an integral part of introgression breeding to ensure effective and durable resistance in wheat. [ABSTRACT FROM AUTHOR]

Analysis of human cancer genome sequences has revealed specific mutational signatures associated with BRCA1-deficient tumors, but the underlying mechanisms remain poorly understood. Here, we show that one-ended DNA double strand breaks (DSBs) converted from CRISPR/Cas9-induced nicks by DNA replication, not two-ended DSBs, cause more characteristic chromosomal aberrations and micronuclei in Brca1-deficient cells than in wild-type cells. BRCA1 is required for efficient homologous recombination of these nick-converted DSBs and suppresses bias towards long tract gene conversion and tandem duplication (TD) mediated by two-round strand invasion in a replication strand asymmetry. However, aberrant repair of these nick-converted one-ended DSBs, not that of two-ended DSBs in Brca1-deficient cells, generates mutational signatures such as small indels with microhomology (MH) at the junctions, translocations and small MH-mediated TDs, resembling those in BRCA1-deficient tumors. These results suggest a major contribution of DNA nicks to mutational signatures associated with BRCA1 deficiency in cancer and the underlying mechanisms. It remains unclear how BRCA1-associated mutational sigantures are generated in cancer. Here, the authors develop nCas9-based approaches to uncover that aberrant repair of one-ended DSBs converted from DNA nicks by replication, not that of two-ended DSBs, induces such mutational signatures. [ABSTRACT FROM AUTHOR]

Blue cone monochromacy (BCM) is an X-linked retinal disorder characterized by low vision, photoaversion, and poor color discrimination. BCM is due to the lack of long-wavelength-sensitive and middle-wavelength-sensitive cone photoreceptor function and caused by mutations in the OPN1LW/OPN1MW gene cluster on Xq28. Here, we investigated the prevalence and the landscape of submicroscopic structural variants (SVs) at single-base resolution in BCM patients. We found that about one-third (n = 73) of the 213 molecularly confirmed BCM families carry an SV, most commonly deletions restricted to the OPN1LW/OPN1MW gene cluster. The structure and precise breakpoints of the SVs were resolved in all but one of the 73 families. Twenty-two families--all from the United States--showed the same SV, and we confirmed a common ancestry of this mutation. In total, 42 distinct SVs were identified, including 40 previously unreported SVs, thereby quadrupling the number of precisely mapped SVs underlying BCM. Notably, there was no "region of overlap" among these SVs. However, 90% of SVs encompass the upstream locus control region, an essential enhancer element. Its minimal functional extent based on deletion mapping in patients was refined to 358 bp. Breakpoint analyses suggest diverse mechanisms underlying SV formation as well as in one case the gene conversion-based exchange of a 142-bp deletion between opsin genes. Using parsimonious assumptions, we reconstructed the composition and copy number of the OPN1LW/OPN1MW gene cluster prior to the mutation event and found evidence that large gene arrays may be predisposed to the occurrence of SVs at this locus. [ABSTRACT FROM AUTHOR]

Arbuscular mycorrhizal fungi (AMF) are part of the most widespread fungal-plant symbiosis. They colonize at least 80% of plant species, promote plant growth and plant diversity. These fungi are multinucleated and contain either one or two haploid nuclear genotypes (monokaryon and dikaryon) identified by the alleles at a putative mating-type locus. This taxon has been considered as an ancient asexual scandal because of the lack of observable sexual structures. Despite identification of a putative mating-type locus and functional activation of genes related to mating when two isolates co-exist, it remains unknown if the AMF life cycle involves a sexual or parasexual stage. We used publicly available genome sequences to test if Rhizophagus irregularis dikaryon genomes display signatures of sexual reproduction in the form of reciprocal recombination patterns, or if they display exclusively signatures of parasexual reproduction involving gene conversion. We used short-read and long-read sequence data to identify nucleus-specific alleles within dikaryons and then compared them to orthologous gene sequences from related monokaryon isolates displaying the same putative MAT-types as the dikaryon. We observed that the two nucleus-specific alleles of the dikaryon A5 are more related to the homolog sequences of monokaryon isolates displaying the same putative MAT-type than between each other. We also observed that these nucleus-specific alleles displayed reciprocal recombination signatures. These results confirm that dikaryon and monokaryon isolates displaying the same putative MAT-type are related in their life-cycle. These results suggest that a genetic exchange mechanism, involving reciprocal recombination in dikaryon genomes, allows AMF to generate genetic diversity. [ABSTRACT FROM AUTHOR]

Chinese chives is a popular herb vegetable and medicine in Asian countries. Southwest China is one of the centers of origin, and the mountainous areas in this region are rich in wild germplasm. In this study, we collected four samples of germplasm from different altitudes: a land race of cultivated Chinese chives (Allium tuberosum), wide-leaf chives and extra-wide-leaf chives (Allium hookeri), and ovoid-leaf chives (Allium funckiaefolium). Leaf metabolites were detected and compared between A. tuberosum and A. hookeri. A total of 158 differentially accumulated metabolites (DAM) were identified by Gas Chromatography—Mass Spectrometry (GC-MS) and Liquid Chromatography—Mass Spectrometry (LC-MS), among which there was a wide range of garlic odor compounds, free amino acids, and sugars. A. hookeri contains a higher content of fructose, garlic odor compounds, and amino acids than A. tuberosum, which is supported by the higher expression level of biosynthetic genes revealed by transcriptome analysis. A. hookeri accumulates the same garlic odor compound precursors that A. tuberosum does (mainly methiin and alliin). We isolated full-length gene sequences of phytochelatin synthase (PCS), γ-glutamyltranspeptidases (GGT), flavin-containing monooxygenase (FMO), and alliinase (ALN). These sequences showed closer relations in phylogenetic analysis between A. hookeri and A. tuberosum (with sequence identities ranging from 86% to 90%) than with Allium cepa or Allium sativum (which had a lower sequence identity ranging from 76% to 88%). Among these assayed genes, ALN, the critical gene controlling the conversion of odorless precursors into odor compounds, was undetected in leaves, bulbs, and roots of A. tuberosum, which could account for its weaker garlic smell. Moreover, we identified a distinct FMO1 gene in extra-wide-leaf A. hookeri that is due to a CDS-deletion and frameshift mutation. These results above reveal the molecular and metabolomic basis of impressive strong odor in wild Chinese chives. [ABSTRACT FROM AUTHOR]

Apart from I GBA i and I PRKRA i , a potential role of pseudogenes has not (yet) caught our attention in the field of movement disorders, although several well-established disease-linked genes also have pseudogenes. The definition of a pseudogene has now broadened to include any DNA sequence that is similar to a known gene and has lost some of its original functionality.25 Therefore, pseudogenes can be translated into proteins, as recently proven for 140 human pseudogenes, and also act at the protein level.26 For instance, pseudogene-derived proteins may have the same activity as the parental proteins, but function in different tissues, cellular compartments, or pathophysiological conditions. Pseudogenes can be divided into two major groups, that is, unprocessed pseudogenes and processed pseudogenes, with the latter representing the majority of pseudogenes in humans (~70%, https://www.gencodegenes.org/human/stats.html). Therefore, the sequence identity for these pseudogenes should be too less to seriously trouble sequence analyses.27 Overall, pseudogenes have emerged as a hot topic because they can challenge genetic testing and trigger mutational events even in biologically functional parental genes. [Extracted from the article]

Patterns of non-uniform usage of synonymous codons vary across genes in an organism and between species across all domains of life. This codon usage bias (CUB) is due to a combination of non-adaptive (e.g. mutation biases) and adaptive (e.g. natural selection for translation efficiency/accuracy) evolutionary forces. Most models quantify the effects of mutation bias and selection on CUB assuming uniform mutational and other non-adaptive forces across the genome. However, non-adaptive nucleotide biases can vary within a genome due to processes such as biased gene conversion (BGC), potentially obfuscating signals of selection on codon usage. Moreover, genome-wide estimates of non-adaptive nucleotide biases are lacking for non-model organisms. We combine an unsupervised learning method with a population genetics model of synonymous coding sequence evolution to assess the impact of intragenomic variation in non-adaptive nucleotide bias on quantification of natural selection on synonymous codon usage across 49 Saccharomycotina yeasts. We find that in the absence of a priori information, unsupervised learning can be used to identify genes evolving under different non-adaptive nucleotide biases. We find that the impact of intragenomic variation in non-adaptive nucleotide bias varies widely, even among closely-related species. We show that the overall strength and direction of translational selection can be underestimated by failing to account for intragenomic variation in non-adaptive nucleotide biases. Interestingly, genes falling into clusters identified by machine learning are also physically clustered across chromosomes. Our results indicate the need for more nuanced models of sequence evolution that systematically incorporate the effects of variable non-adaptive nucleotide biases on codon frequencies. Author summary: Codon usage bias (CUB), or the unequal usage of codons of the same amino acid (i.e. synonymous codons), has been observed in species across all domains of life. CUB is known to be shaped by both non-adaptive (e.g. mutation biases) and adaptive (e.g. natural selection for translation efficiency/accuracy) evolution. A key challenge for researchers is disentangling the role of various processes shaping codon usage, often for the purpose of identifying codons favored by natural selection, sometimes referred to as "optimal" or "preferred" codons. Despite large variation in non-adaptive nucleotide biases within a genome, most methods to quantify natural selection typically ignore this variation for the sake of simplicity. Here, we combine a population genetics model with unsupervised machine learning to identify genes evolving under different non-adaptive nucleotide biases across 49 budding yeasts species. We find that ignoring for variation in non-adaptive nucleotide biases can obfuscate signals of selection on codon usage. Our results indicate the need for more nuanced models of coding sequence evolution. [ABSTRACT FROM AUTHOR]

Long chain acyl-coA synthase (acsl) family genes activate the conversion of long chain fatty acids into acyl-coA to regulate fatty acid metabolism. However, the evolutionary characteristics, tissue expression and nutritional regulation of the acsl gene family are poorly understood in fish. The present study investigated the molecular characterization, tissue expression and nutritional regulation of the acsl gene family in golden pompano (Trachinotus ovatus). The results showed that the coding regions of acsl1, acsl3, acsl4, acsl5 and acsl6 cDNA were 2091 bp, 2142 bp, 2136 bp, 1977 bp and 2007 bp, encoding 697, 714, 712, 659 and 669 amino acids, respectively. Five acsl isoforms divided into two branches, namely, acsl1, acsl5 and acsl6, as well as acsl3 and acsl4. The tissue expression distribution of acsl genes showed that acsl1 and acsl3 are widely expressed in the detected tissues, while acsl4, acsl5 and acsl6 are mainly expressed in the brain. Compared to the fish fed with lard oil diets, the fish fed with soybean oil exhibited high muscular C18 PUFA contents and acsl1 and acsl3 mRNA levels, as well as low muscular SFA contents and acsl4 mRNA levels. High muscular n-3 LC-PUFA contents, and acsl3, acsl4 and acsl6 mRNA levels were observed in the fish fed with fish oil diets compared with those of fish fed with lard oil or soybean oil diets. High n-3 LC-PUFA levels and DHA contents, as well as the acsl3, acsl4 and acsl6 mRNA levels were exhibited in the muscle of fish fed diets with high dietary n-3 LC-PUFA levels. Additionally, the muscular acsl3, acsl4 and acsl6 mRNA expression levels, n-3 LC-PUFA and DHA levels were significantly up-regulated by the increase of dietary DHA proportions. Collectively, the positive relationship among dietary fatty acids, muscular fatty acids and acsl mRNA, indicated that T. ovatus Acsl1 and Acsl3 are beneficial for the C18 PUFA enrichment, and Acsl3, Acsl4 and Acsl6 are for n-3 LC-PUFA and DHA enrichment. The acquisition of fish Acsl potential function in the present study will play the foundation for ameliorating the fatty acids nutrition in farmed fish products. [ABSTRACT FROM AUTHOR]

Background: Gene conversion has an important effect on duplicate genes produced by polyploidization. Poplar (Populus trichocarpa) and willow (Salix brachista) are leading models and excellent green plants in the Salicaceae. Although much attention has been paid to the evolution of duplicated genes in poplar and willow, the role of conversion between duplicates generated from polyploidization remains poorly understood. Results: Here, through genomic synteny analyses, we identified duplicate genes generated by the Salicaceae common tetraploidization (SCT) in the poplar and willow genomes. We estimated that at least 0.58% and 0.25% of poplar and willow duplicates were affected by whole-gene conversion after the poplar-willow divergence, with more (5.73% and 2.66%) affected by partial-gene conversion. Moreover, we found that the converted duplicated genes were unevenly distributed on each chromosome in the two genomes, and the well-preserved homoeologous chromosome regions may facilitate the conversion of duplicates. Notably, we found that conversion maintained the similarity of duplicates, likely contributing to the conservation of certain sequences, but is essentially accelerated the rate of evolution and increased species divergence. In addition, we found that converted duplicates tended to have more similar expression patterns than nonconverted duplicates. We found that genes associated with multigene families were preferentially converted. We also found that the genes encoding conserved structural domains associated with specific traits exhibited a high frequency of conversion. Conclusions: Extensive conversion between duplicate genes generated from the SCT contributes to the diversification of the family Salicaceae and has had long-lasting effects on those genes with important biological functions. [ABSTRACT FROM AUTHOR]

Simple Summary: In this work we present MAGE, an open-source Python package developed for the meta-analysis of gene expression data. It contains functions to convert probes to gene identifiers, and to perform standard meta-analysis, meta-analysis with bootstrap standard errors, and meta-analysis of multiple outcomes, as well as functional enrichment analysis. Additionally, visualizations for every function of this software package are provided. MAGE is available both in a standalone version and as a webserver. MAGE (Meta-Analysis of Gene Expression) is a Python open-source software package designed to perform meta-analysis and functional enrichment analysis of gene expression data. We incorporate standard methods for the meta-analysis of gene expression studies, bootstrap standard errors, corrections for multiple testing, and meta-analysis of multiple outcomes. Importantly, the MAGE toolkit includes additional features for the conversion of probes to gene identifiers, and for conducting functional enrichment analysis, with annotated results, of statistically significant enriched terms in several formats. Along with the tool itself, a web-based infrastructure was also developed to support the features of this package. [ABSTRACT FROM AUTHOR]

Falcons are diverse birds of cultural and economic importance. They have undergone major lineage-specific chromosomal rearrangements, resulting in greatly reduced chromosome counts relative to other birds. Here, we use 10X Genomics linked reads to provide new high-contiguity genomes for two gyrfalcons, a saker falcon, a lanner falcon, three subspecies of peregrine falcons, and the common kestrel. Assisted by a transcriptome sequenced from 22 gyrfalcon tissues, we annotate these genomes for a variety of genomic features, estimate historical demography, and then investigate genomic equilibrium in the context of falcon-specific chromosomal rearrangements. We find that falcon genomes are not in AT–GC equilibrium with a bias in substitutions toward higher AT content; this bias is predominantly but not exclusively driven by hypermutability of CpG sites. Small indels and large structural variants were also biased toward insertions rather than deletions. Patterns of disequilibrium were linked to chromosomal rearrangements: falcons have lost GC content in regions that have fused to larger chromosomes from microchromosomes and gained GC content in regions of macrochromosomes that have translocated to microchromosomes. Inserted bases have accumulated on regions ancestrally belonging to microchromosomes, consistent with insertion-biased gene conversion. We also find an excess of interspersed repeats on regions of microchromosomes that have fused to macrochromosomes. Our results reveal that falcon genomes are in a state of flux. They further suggest that many of the key differences between microchromosomes and macrochromosomes are driven by differences in chromosome size, and indicate a clear role for recombination and biased gene conversion in determining genomic equilibrium. [ABSTRACT FROM AUTHOR]

Aurantiochytrium mangrovei BT3, a heterotrophic marine microalga, has the ability to produce high amounts of polyunsaturated fatty acids and bioactive compounds, which provide many benefits to humans and animals. In this study, (24R)-4α-methyl-5α-stigmasta-7,22-dien-3β-ol (denoted AME1) was isolated from A. mangrovei strain BT3, and its antihyperlipidaemic effects were investigated. From 2 kg of the dried biomass of A. mangrovei BT3, 30 mg of AME1 with 94% purity was achieved. Stimulation of the HepG2 cell line with AME1 decreased intracellular cholesterol and triglycerides in a dose-dependent manner. Hypolipidaemic metabolism analysis revealed that the lipid-lowering effect of AME1 was due to increased expression of the CYP7A1 gene and decreased expression of the gene involved in fatty acid biosynthesis, SCD-1. In addition, AME1 directly stimulated the activity of PPARα receptors and resulted in increased expression of genes involved in fatty acid transport (FATP4) and fatty acid oxidation (CTP-1, MCAD, ACOX) and decreased triglyceride accumulation. The obtained results suggest that AME1 isolated from A. mangrovei BT3 inhibits hepatic lipid accumulation by promoting the expression of genes involved in the conversion of excess cholesterol to bile acids and lipogenesis and inducing the transcriptional activity of PPARα. [ABSTRACT FROM AUTHOR]

We described the population structure of Bordetella pertussis (B. pertussis) in Norway from 1996 to 2019 and determined if there were evolutionary shifts and whether these correlated with changes in the childhood immunization program. We selected 180 B. pertussis isolates, 22 from the whole cell vaccine (WCV) era (1996–1997) and 158 from the acellular vaccine (ACV) era (1998–2019). We conducted whole genome sequencing and determined the distribution and frequency of allelic variants and temporal changes of ACV genes. Norwegian B. pertussis isolates were evenly distributed across a phylogenetic tree that included global strains. We identified seven different allelic profiles of ACV genes (A–F), in which profiles A1, A2, and B dominated (89%), all having pertussis toxin (ptxA) allele 1, pertussis toxin promoter (ptxP) allele 3, and pertactin (prn) allele 2 present. Isolates with ptxP1 and prn1 were not detected after 2007, whereas the prn2 allele likely emerged prior to 1972, and ptxP3 before the early 1980s. Allele conversions of ACV genes all occurred prior to the introduction of ACV. Sixteen percent of our isolates showed mutations within the prn gene. ACV and its booster doses (implemented for children in 2007 and adolescents in 2013) might have contributed to evolvement of a more uniform B. pertussis population, with recent circulating strains having ptxA1, ptxP3, and prn2 present, and an increasing number of prn mutations. These strains clearly deviate from ACV strains (ptxA1, ptxP1, prn1), and this could have implications for vaccine efficiency and, therefore, prevention and control of pertussis. [ABSTRACT FROM AUTHOR]

Leucine-rich repeat proteins and antimicrobial peptides are the key components of the innate immune response to Plasmodium and other microbial pathogens in Anopheles mosquitoes. The APL1 gene of the malaria vector Anopheles funestus has exceptional levels of non-synonymous polymorphism across the range of An. funestus, with an average πn of 0.027 versus a genome-wide average of 0.002, and πn is consistently high in populations across Africa. Elevated APL1 diversity was consistent between the independent pooled-template and target-enrichment datasets, however no link between APL1 diversity and insecticide resistance was observed. Although lacking the diversity of APL1, two further mosquito innate-immunity genes of the gambicin anti-microbial peptide family had πn/πs ratios greater than one, possibly driven by either positive or balancing selection. The cecropin antimicrobial peptides were expressed much more highly than other anti-microbial peptide genes, a result discordant with current models of anti-microbial peptide activity. The observed APL1 diversity likely results from gene conversion between paralogues, as evidenced by shared polymorphisms, overlapping read mappings, and recombination events among paralogues. In conclusion, we hypothesize that higher gene expression of APL1 than its paralogues is correlated with a more open chromatin formation, which enhances gene conversion and elevated diversity at this locus. [ABSTRACT FROM AUTHOR]

Mating-type (P or M) of fission yeast Schizosaccharomyces pombe is determined by the transcriptionally active mat1 cassette and is switched by gene conversion using a donor, either mat2 or mat3, located in an adjacent heterochromatin region (mating-type switching; MTS). In the switching process, heterochromatic donors of genetic information are selected based on the P or M cell type and on the action of two recombination enhancers, SRE2 promoting the use of mat2-P and SRE3 promoting the use of mat3-M, leading to replacement of the content of the expressed mat1 cassette. Recently, we found that the histone H3K4 methyltransferase complex Set1C participates in donor selection, raising the question of how a complex best known for its effects in euchromatin controls recombination in heterochromatin. Here, we report that the histone H2BK119 ubiquitin ligase complex HULC functions with Set1C in MTS, as mutants in the shf1, brl1, brl2 and rad6 genes showed defects similar to Set1C mutants and belonged to the same epistasis group as set1?. Moreover, using H3K4R and H2BK119R histone mutants and a Set1-Y897A catalytic mutant, we found that ubiquitylation of histone H2BK119 by HULC and methylation of histone H3K4 by Set1C are functionally coupled in MTS. Cell-type biases in MTS in these mutants suggested that HULC and Set1C inhibit the use of the SRE3 recombination enhancer in M cells, thus favoring SRE2 and mat2-P. Consistent with this, imbalanced switching in the mutants was traced to compromised association of the directionality factor Swi6 with the recombination enhancers in M cells. Based on their known effects at other chromosomal locations, we speculate that HULC and Set1C control nucleosome mobility and strand invasion near the SRE elements. In addition, we uncovered distinct effects of HULC and Set1C on histone H3K9 methylation and gene silencing, consistent with additional functions in the heterochromatic domain. [ABSTRACT FROM AUTHOR]

BackgroundSiglecs-11 and -16 are members of the sialic acid recognizing Ig-like lectin family, and expressed in same cells. Siglec-11 functions as an inhibitory receptor, whereas Siglec-16 exhibits activating properties. In humans, SIGLEC11 and SIGLEC16 gene sequences are extremely similar in the region encoding the extracellular domain due to gene conversions. Human SIGLEC11 was converted by the nonfunctional SIGLEC16P allele, and the converted SIGLEC11 allele became fixed in humans, possibly because it provides novel neuroprotective functions in brain microglia. However, the detailed evolutionary history of SIGLEC11 and SIGLEC16 in other primates remains unclear.ResultsWe analyzed SIGLEC11 and SIGLEC16 gene sequences of multiple primate species, and examined glycan binding profiles of these Siglecs. The phylogenetic tree demonstrated that gene conversions between SIGLEC11 and SIGLEC16 occurred in the region including the exon encoding the sialic acid binding domain in every primate examined. Functional assays showed that glycan binding preference is similar between Siglec-11 and Siglec-16 in all analyzed hominid species. Taken together with the fact that Siglec-11 and Siglec-16 are expressed in the same cells, Siglec-11 and Siglec-16 are regarded as paired receptors that have maintained similar ligand binding preferences via gene conversions. Relaxed functional constraints were detected on the SIGLEC11 and SIGLEC16 exons that underwent gene conversions, possibly contributing to the evolutionary acceptance of repeated gene conversions. The frequency of nonfunctional SIGLEC16P alleles is much higher than that of SIGLEC16 alleles in every human population.ConclusionsOur findings indicate that Siglec-11 and Siglec-16 have been maintained as paired receptors by repeated gene conversions under relaxed functional constraints in the primate lineage. The high prevalence of the nonfunctional SIGLEC16P allele and the fixation of the converted SIGLEC11 imply that the loss of Siglec-16 and the gain of Siglec-11 in microglia might have been favored during the evolution of human lineage.

In yeast, DNA breaks are usually repaired by homologous recombination (HR). An early step for HR pathways is formation of a heteroduplex, in which a single-strand from the broken DNA molecule pairs with a strand derived from an intact DNA molecule. If the two strands of DNA are not identical, there will be mismatches within the heteroduplex DNA (hetDNA). In wild-type strains, these mismatches are repaired by the mismatch repair (MMR) system, producing a gene conversion event. In strains lacking MMR, the mismatches persist. Most previous studies involving hetDNA formed during mitotic recombination were restricted to one locus. Below, we present a global mapping of hetDNA formed in the MMR-defective mlh1 strain. We find that many recombination events are associated with repair of double-stranded DNA gaps and/or involve Mlh1-independent mismatch repair. Many of our events are not explicable by the simplest form of the double-strand break repair model of recombination.

Gene clusters of recently duplicated genes are hotbeds for evolutionary change. However, our understanding of how mutational mechanisms and evolutionary forces shape the structural and functional evolution of these clusters is hindered by the high sequence identity among the copies, which typically results in their inaccurate representation in genome assemblies. The presumed testis-specific, chimeric gene Sdic originated, and tandemly expanded in Drosophila melanogaster, contributing to increased male-male competition. Using various types of massively parallel sequencing data, we studied the organization, sequence evolution, and functional attributes of the different Sdic copies. By leveraging long-read sequencing data, we uncovered both copy number and order differences from the currently accepted annotation for the Sdic region. Despite evidence for pervasive gene conversion affecting the Sdic copies, we also detected signatures of two episodes of diversifying selection, which have contributed to the evolution of a variety of C-termini and miRNA binding site compositions. Expression analyses involving RNA-seq datasets from 59 different biological conditions revealed distinctive expression breadths among the copies, with three copies being transcribed in females, opening the possibility to a sexually antagonistic effect. Phenotypic assays using Sdic knock-out strains indicated that should this antagonistic effect exist, it does not compromise female fertility. Our results strongly suggest that the genome consolidation of the Sdic gene cluster is more the result of a quick exploration of different paths of molecular tinkering by different copies than a mere dosage increase, which could be a recurrent evolutionary outcome in the presence of persistent sexual selection.

The transport of polarons above the mobility threshold in organic and inorganic polymers is theoretically investigated in the framework of a one-dimensional tight-binding model that includes lattice relaxation. The computational approach is based on parameters for which the model Hamiltonian suitably describes different polymer lattices in the presence of external electric fields. Our findings show that, above critical field strengths, a dissociated polaron moves through the polymer lattice as a free electron performing Bloch oscillations. These critical electric fields are considerably smaller for inorganic lattices in comparison to organic polymers. Interestingly, for inorganic lattices, the free electron propagates preserving charge and spin densities' localization which is a characteristic of a static polaron. Moreover, in the turning points of the spatial Bloch oscillations, transient polaron levels are formed inside the band gap, thus generating a fully characterized polaron structure. For the organic case, on the other hand, no polaron signature is observed: neither in the shape of the distortion-those polaron profile signatures are absent-nor in the energy levels-as no such polaron levels are formed during the simulation. These results solve controversial aspects concerning Bloch oscillations recently reported in the literature and may enlighten the understanding about the charge transport mechanism in polymers above their mobility edge. [ABSTRACT FROM AUTHOR]

[Display omitted] • SP-EP enhanced the production of ethanol, acetate and methane. • SP-EP enriched methanogens such as Methanothrix and Methanoculleus. • EP promoted the conversion of ethanol to acetaldehyde. • SP promoted the production of CO 2 and H 2 from pyruvate. • SP-EP enhanced acetyl-CoA production and hydrogenotrophic methanogenic pathway. In anaerobic digestion (AD) of food waste (FW), the acetogenesis stage is thermodynamically limited and often leads to acid inhibition. This study investigated the impact of optimizing ethanol pre-fermentation (EP) of FW through saccharification pretreatment (SP) on AD by regulating the acidification pathway of substrates. The results showed that the ethanol concentration of FW after SP-EP reached 50.81 g/L, which was 2.11-fold of that in EP. SP-EP enriched hydrolytic and acidogenic bacteria such as Alkaliphilus , Caldicoprobacter , Sphaerochaeta , and Fermentimonas. The analysis showed that the abundance of enzyme-encoding genes catalyzing the conversion of ethanol to acetaldehyde in SP-EP was 1.17-fold higher than control, which further indicated that SP-EP facilitated the conversion of substrate to glucose → ethanol → acetaldehyde → acetate. This process provided abundant acetate for AD, and enriched the acetoclastic methanogens Methanothrix. Consequently, the cumulative methane yield of SP-EP (557 mL/g vs) was 7.49 %–17.79 % higher than the other three treatments. Hydrogenotrophic methanogens such as Methanoculleus were also enriched in SP-EP due to the fact that acetogenesis was accompanied by hydrogen production, which enhanced the hydrogenotrophic methanogenic pathway. This study provides a potential strategy for balancing the processes of acetogenesis and methanogenesis, providing a practical basis for the industrial application of methane production from FW. [ABSTRACT FROM AUTHOR]

Phytoremediation has emerged as a common technique for remediating Cd pollution in farmland soil. Moreover, phosphorus, an essential element for plants, can alter the pectin content of plant cell walls and facilitate the accumulation of Cd in plant tissues, thereby enhancing phytoremediation efficiency. Therefore, pot experiments were conducted in order to investigate the effect of phosphorus levels on Cd extraction, phosphorus transformation and phosphorus-related genes during phytoremediation. The results revealed that an optimal application of suitable phosphate fertilizers elevated the soil's pH and electrical conductivity (EC), facilitated the conversion of soil from insoluble phosphorus into available forms, augmented the release of pertinent enzyme activity, and induced the expression of phosphorus cycling-related genes. These enhancements in soil conditions significantly promoted the growth of ryegrass. When applying phosphorus at a rate of 600 mg/kg, ryegrass exhibited plant height, dry weight, and chlorophyll relative content that were 1.27, 1.26, and 1.18 times higher than those in the control group (P 0), while the Cd content was 1.12 times greater than that of P 0. The potentially toxic elements decline ratio and bioconcentration factor were 42.86% and 1.17 times higher than those of P 0 , respectively. Consequently, ryegrass demonstrated the highest Cd removal efficiency under these conditions. Results from redundancy analysis (RDA) revealed a significant correlation among pH, total phosphorus, heavy metal content, phosphorus forms, soil enzyme activity, and phosphorus-related genes. In conclusion, this study suggests applying an optimal amount of suitable phosphate fertilizers can enhance restoration efficiency, leading to a reduction in soil Cd content and ultimately improving the safety of crop production in farmlands. [Display omitted] • The combination of P and EC-OG improved the growth of ryegrass. • P and EC-OG association drives changes in P morphology and enzymes. • 600 mg/kg P had the best effect on the phytoremediation efficiency. • P conversion and gene expression were optimal when P level at 600 mg/kg. • RDA showed that pH and TP had the greatest influence on the whole process. [ABSTRACT FROM AUTHOR]

Meiotic drive supergenes are complexes of alleles at linked loci that together subvert Mendelian segregation resulting in preferential transmission. In males, the most common mechanism of drive involves the disruption of sperm bearing one of a pair of alternative alleles. While at least two loci are important for male drive--the driver and the target--linked modifiers can enhance drive, creating selection pressure to suppress recombination. In this work, we investigate the evolution and genomic consequences of an autosomal, multilocus, male meiotic drive system, Segregation Distorter (SD) in the fruit fly, Drosophila melanogaster. In African populations, the predominant SD chromosome variant, SD-Mal, is characterized by two overlapping, paracentric inversions on chromosome arm 2R and nearly perfect (~100%) transmission. We study the SD-Mal system in detail, exploring its components, chromosomal structure, and evolutionary history. Our findings reveal a recent chromosome-scale selective sweep mediated by strong epistatic selection for haplotypes carrying Sd, the main driving allele, and one or more factors within the double inversion. While most SD-Mal chromosomes are homozygous lethal, SD-Mal haplotypes can recombine with other, complementing haplotypes via crossing over, and with wildtype chromosomes via gene conversion. SD-Mal chromosomes have nevertheless accumulated lethal mutations, excess non-synonymous mutations, and excess transposable element insertions. Therefore, SD-Mal haplotypes evolve as a small, semi-isolated subpopulation with a history of strong selection. These results may explain the evolutionary turnover of SD haplotypes in different populations around the world and have implications for supergene evolution broadly. [ABSTRACT FROM AUTHOR]

The Lyme disease spirochete Borrelia burgdorferi , encodes an elaborate antigenic variation system that promotes the ongoing variation of a major surface lipoprotein, VlsE. Changes in VlsE are continual and always one step ahead of the host acquired immune system, which requires 1–2 weeks to generate specific antibodies. By the time this happens, new VlsE variants have arisen that escape immunosurveillance, providing an avenue for persistent infection. This antigenic variation system is driven by segmental gene conversion events that transfer information from a series of silent cassettes (vls2-16) to the expression locus, vlsE. The molecular details of this process remain elusive. Recombinational switching at vlsE is RecA-independent and the only required factor identified to date is the RuvAB branch migrase. In this work we have used next generation long-read sequencing to analyze the effect of several DNA replication/recombination/repair gene disruptions on the frequency of gene conversions at vlsE and report a requirement for the mismatch repair protein MutL. Site directed mutagenesis of mutL suggests that the putative MutL endonuclease activity is required for recombinational switching at vlsE. This is the first report of an unexpected essential role for MutL in a bacterial recombination system and expands the known function of this protein as well as our knowledge of the details of the novel recombinational switching mechanism for vlsE variation. [ABSTRACT FROM AUTHOR]

The assumption that conservation of sequence implies the action of purifying selection is central to diverse methodologies to infer functional importance. GC-biased gene conversion (gBGC), a meiotic mismatch repair bias strongly favouring GC over AT, can in principle mimic the action of selection, this being thought to be especially important in mammals. As mutation is GC→AT biased, to demonstrate that gBGC does indeed cause false signals requires evidence that an AT-rich residue is selectively optimal compared to its more GC-rich allele, while showing also that the GC-rich alternative is conserved. We propose that mammalian stop codon evolution provides a robust test case. Although in most taxa TAA is the optimal stop codon, TGA is both abundant and conserved in mammalian genomes. We show that this mammalian exceptionalism is well explained by gBGC mimicking purifying selection and that TAA is the selectively optimal codon. Supportive of gBGC, we observe (i) TGA usage trends are consistent at the focal stop codon and elsewhere (in UTR sequences); (ii) that higher TGA usage and higher TAA→TGA substitution rates are predicted by a high recombination rate; and (iii) across species the difference in TAA <-> TGA substitution rates between GC-rich and GC-poor genes is largest in genomes that possess higher between-gene GC variation. TAA optimality is supported both by enrichment in highly expressed genes and trends associated with effective population size. High TGA usage and high TAA→TGA rates in mammals are thus consistent with gBGC's predicted ability to "drive" deleterious mutations and supports the hypothesis that sequence conservation need not be indicative of purifying selection. A general trend for GC-rich trinucleotides to reside at frequencies far above their mutational equilibrium in high recombining domains supports the generality of these results. Is sequence conservation a sign of purifying selection and hence functional importance? This analysis of why mammals use and conserve the most error-prone stop codon suggests not, consistent with GC-biased gene conversion's predicted ability to "drive" deleterious mutations and supporting the hypothesis that sequence conservation need not be indicative of purifying selection. [ABSTRACT FROM AUTHOR]

Paucilactobacillus wasatchensis can use gluconate (GLCN) as well as galactose as an energy source and because sodium GLCN can be added during salting of Cheddar cheese to reduce calcium lactate crystal formation, our primary objective was to determine if the presence of GLCN in cheese is another risk factor for unwanted gas production leading to slits in cheese. A secondary objective was to calculate the amount of CO 2 produced during storage and to relate this to the amount of gas-forming substrate that was utilized. Ribose was added to promote growth of Pa. wasatchensis WDC04 (P.wa WDC04) to high numbers during storage. Cheddar cheese was made with lactococcal starter culture with addition of P.wa WDC04 on 3 separate occasions. After milling, the curd was divided into six 10-kg portions. To the curd was added (A) salt, or salt plus (B) 0.5% galactose + 0.5% ribose (similar to previous studies), (C) 1% sodium GLCN, (D) 1% sodium GLCN + 0.5% ribose, (E) 2% sodium GLCN, (F) 2% sodium GLCN + 0.5% ribose. A vat of cheese without added P.wa WDC04 was made using the same milk and a block of cheese used as an additional control. Cheeses were cut into 900-g pieces, vacuum packaged and stored at 12°C for 16 wk. Each month the bags were examined for gas production and cheese sampled and tested for lactose, galactose and GLCN content, and microbial numbers. In the control cheese, P.wa WDC04 remained undetected (i.e., <104 cfu/g), whereas in cheeses A, C, and E it increased to 107 cfu/g, and when ribose was included with salting (cheeses B, D, and F) increased to 108 cfu/g. The amount of gas (measured as headspace height or calculated as mmoles of CO 2) during 16 wk storage was increased by adding P.wa WDC04 into the milk, and by adding galactose or GLCN to the curd. Galactose levels in cheese B were depleted by 12 wk while no other cheeses had residual galactose. Except for cheese D, the other cheeses with GLCN added (C, E and F) showed little decline in GLCN levels until wk 12, even though gas was being produced starting at wk 4. Based on calculations of CO 2 in headspace plus CO 2 dissolved in cheese, galactose and GLCN added to cheese curd only accounted for about half of total gas production. It is proposed that CO 2 was also produced by decarboxylation of amino acids. Although P.wa WDC04 does not have all the genes for complete conversion and decarboxylation of the amino acids in cheese, this can be achieved in conjunction with starter culture lactococcal. Adding GLCN to curd can now be considered another confirmed risk factor for unwanted gas production during storage of Cheddar cheese that can lead to slits and cracks in cheese. Putative risk factors now include having a community of bacteria in cheese leading to decarboxylation of amino acids and release of CO 2 as well autolysis of the starter culture that would provide a supply of ribose that can promote growth of Pa. wasatchensis. [ABSTRACT FROM AUTHOR]

Key message: A major gene controls flowering pattern in peanut, possibly encoding a TFL1-like. It was subjected to gain/loss events of a deletion and changes in mRNA expression levels, partly explaining the evolution of flowering pattern in Arachis. Flowering pattern (FP) is a major characteristic differentiating the two subspecies of cultivated peanut (Arachis hypogaea L.). Subsp. fastigiata possessing flowers on the mainstem (MSF) and a sequential FP, whereas subsp. hypogaea lacks MSF and exhibits an alternate FP. FP is considered the main contributor to plant adaptability, and evidence indicates that its diversification occurred during the several thousand years of domestication. However, the genetic mechanism that controls FP in peanut is unknown. We investigated the genetics of FP in a recombinant inbred population, derivatives of an A. hypogaea by A. fastigiata cross. Lines segregated 1:1 for FP, indicating a single gene effect. Using Axiom_Arachis2 SNP-array, FP was mapped to a small segment in chromosome B02, wherein a Terminal Flowering 1-like (AhTFL1) gene with a 1492 bp deletion was found in the fastigiata line, leading to a truncated protein. Remapping FP in the RIL population with the AhTFL1 indel as a marker increased the LOD score from 53.3 to 158.8 with no recombination in the RIL population. The same indel was found co-segregating with the phenotype in two independent EMS-mutagenized M2 families, suggesting a hotspot for gene conversion. Also, AhTFL1 was significantly less expressed in the fastigiata line compared to hypogaea and in flowering than non-flowering branches. Sequence analysis of the AhTFL1 in peanut world collections indicated significant conservation, supporting the putative role of AhTFL1 in peanut speciation during domestication and modern cultivation. [ABSTRACT FROM AUTHOR]

In bony vertebrates, skeletal mineralization relies on the secretory calcium-binding phosphoproteins (Scpp) family whose members are acidic extracellular proteins posttranslationally regulated by the Fam20°C kinase. As scpp genes are absent from the elephant shark genome, they are currently thought to be specific to bony fishes (osteichthyans). Here, we report a scpp gene present in elasmobranchs (sharks and rays) that evolved from local tandem duplication of sparc-L 5′ exons and show that both genes experienced recent gene conversion in sharks. The elasmobranch scpp is remarkably similar to the osteichthyan scpp members as they share syntenic and gene structure features, code for a conserved signal peptide, tyrosine-rich and aspartate/glutamate-rich regions, and harbor putative Fam20°C phosphorylation sites. In addition, the catshark scpp is coexpressed with sparc-L and fam20°C in tooth and scale ameloblasts, similarly to some osteichthyan scpp genes. Despite these strong similarities, molecular clock and phylogenetic data demonstrate that the elasmobranch scpp gene originated independently from the osteichthyan scpp gene family. Our study reveals convergent events at the sparc-L locus in the two sister clades of jawed vertebrates, leading to parallel diversification of the skeletal biomineralization toolkit. The molecular evolution of sparc-L and its coexpression with fam20°C in catshark ameloblasts provides a unifying genetic basis that suggests that all convergent scpp duplicates inherited similar features from their sparc-L precursor. This conclusion supports a single origin for the hypermineralized outer odontode layer as produced by an ancestral developmental process performed by Sparc-L, implying the homology of the enamel and enameloid tissues in all vertebrates. [ABSTRACT FROM AUTHOR]

Genes of the Major Histocompatibility Complex (MHC) form a key component of vertebrate adaptive immunity, as they code for molecules which bind antigens of intra- and extracellular pathogens (MHC class I and II, respectively) and present them to T cell receptors. In general, MHC genes are hyper-polymorphic and high MHC diversity is often maintained within natural populations (via balancing selection) and within individuals (via gene duplications). Because of its complex architecture with tandems of duplicated genes, characterization of MHC region in non-model vertebrate species still poses a major challenge. Here, we combined de novo genome assembly and high-throughput sequencing to characterize MHC polymorphism in a rallid bird species, the Eurasian coot Fulica atra. An analysis of genome assembly indicated high duplication rate at MHC-I, which was also supported by targeted sequencing of peptide-binding exons (at least five MHC-I loci genotyped). We found high allelic richness at both MHC-I and MHC-II, although signature of diversifying selection and recombination (gene conversion) was much stronger at MHC-II. Our results indicate that Eurasian coot retains extraordinary polymorphism at both MHC classes (when compared to other non-passerine bird species), although they may be subject to different evolutionary mechanism. [ABSTRACT FROM AUTHOR]

Microbial consortia-based consolidated bioprocessing (CBP) is a promising trend in biomass biorefinery, but still faces challenges in terms of complicated microbial structure, low conversion efficiency and operation instability. This study constructed an artificial consortium with improved biomass conversion capability and good operation stability. Diversity of the microbial community structure and gene functions of the domesticated consortium was then analyzed, finding that it had a simplified microbial structure and aggregation of functional genes related to conversion of cellulosic materials for biofuel production. Finally, CBP of wheat straw was performed using the domesticated consortium, ethanol and solvent production with the highest yield ever reported at 0.37 g/g and 0.60 g/g respectively were achieved. Our results further highlight the potential of this domesticated consortium in lignocellulosic biomass biorefinery in comparison with previously reported microbial consortia. Overall, this study provides a guidance on the artificial construction of simplified functional consortia for producing valuable chemicals in a sustainable way. [Display omitted] [ABSTRACT FROM AUTHOR]

Lytic and lysogenic infections are the main strategies used by viruses to interact with microbial hosts. The genetic information of prophages provides insights into the nature of phages and their potential influences on hosts. Here, the siphovirus vB_MoxS-R1 was induced from a Microbacterium strain isolated from an estuarine Synechococcus culture. vB_MoxS-R1 has a high replication capability, with an estimated burst size of 2000 virions per cell. vB_MoxS-R1 represents a novel phage genus-based genomic analysis. Six transcriptional regulator (TR) genes were predicted in the vB_MoxS-R1 genome. Four of these TR genes are involved in stress responses, virulence and amino acid transportation in bacteria, suggesting that they may play roles in regulating the host cell metabolism in response to external environmental changes. A glycerophosphodiester phosphodiesterase gene related to phosphorus acquisition was also identified in the vB_MoxS-R1 genome. The presence of six TR genes and the phosphorus-acquisition gene suggests that prophage vB_MoxS-R1 has the potential to influence survival and adaptation of its host during lysogeny. Possession of four endonuclease genes in the prophage genome suggests that vB_MoxS-R1 is likely involved in DNA recombination or gene conversion and further influences host evolution. [ABSTRACT FROM AUTHOR]

The strength of selection varies among populations and across the genome, but the determinants of efficacy of selection remain unclear. In this study, we used whole-genome sequencing data from 467 Boechera stricta accessions to quantify the strength of selection and characterize the pattern of local adaptation. We found low genetic diversity on 0-fold degenerate sites and conserved non-coding sites, indicating functional constraints on these regions. The estimated distribution of fitness effects and the proportion of fixed substitutions suggest relaxed negative and positive selection in B. stricta. Among the four population groups, the NOR and WES groups have smaller effective population size (Ne), higher proportions of effectively neutral sites, and lower rates of adaptive evolution compared with UTA and COL groups, reflecting the effect of Ne on the efficacy of natural selection. We also found weaker selection on GC-biased sites compared with GC-conservative (unbiased) sites, suggested that GC-biased gene conversion has affected the strength of selection in B. stricta. We found mixed evidence for the role of the recombination rate on the efficacy of selection. The positive and negative selection was stronger in high-recombination regions compared with low-recombination regions in COL but not in other groups. By scanning the genome, we found different subsets of selected genes suggesting differential adaptation among B. stricta groups. These results show that differences in effective population size, nucleotide composition, and recombination rate are important determinants of the efficacy of selection. This study enriches our understanding of the roles of natural selection and local adaptation in shaping genomic variation. [ABSTRACT FROM AUTHOR]

The dog leukocyte antigen (DLA) class I genomic region is located on chromosome 12, and the class I genomic region is composed of at least two distinct haplotypic gene structures, DLA-88–DLA-12 and DLA-88–DLA-88L. However, detailed information of the genomic differences among DLA-88, DLA-12, and DLA-88L are still lacking at the full-length gene level, and therefore, DLA allelic sequences classified for each of these loci are limited in number so far. In this study, we determined the DNA sequence of a 95-kb DLA class I genomic region including DLA-88, DLA-12/88L, and DLA-64 with three DLA homozygous dogs and of 37 full-length allelic gene sequences for DLA-88 and DLA-12/88L loci in 26 DLA class I homozygous dogs. Nucleotide diversity profiles of the 95-kb regions and sequence identity scores of the allelic sequences suggested that DLA-88L is a hybrid gene generated by interlocus and/or intralocus gene conversion between DLA-88 and DLA-12. The putative minimum conversion tract was estimated to be at least an 850-bp segment in length located from the 5´flanking untranslated region to the end of intron 2. In addition, at least one DLA-12 allele (DLA-12*004:01) was newly generated by interlocus gene conversion. In conclusion, the analysis for the occurrence of gene conversion within the dog DLA class I region revealed intralocus gene conversion tracts in 17 of 27 DLA-88 alleles and two of 10 DLA-12 alleles, suggesting that intralocus gene conversion has played an important role in expanding DLA allelic variations. [ABSTRACT FROM AUTHOR]

The Long non-coding RNA (lncRNA) expression profile data of ten samples including human Mesenchymal Stem Cell (MSC) adipogenic differentiation 0, 3, and 6 days from the GEO database, and then perform gene ID conversion, BLAST comparison, and annotation marking. Finally, group A (treatment group on day 3 of differentiation and control group on day 0 of differentiation) obtained a total of 1180 mRNA and 185 lncRNA; group B (treatment group on day 6 of differentiation and control group on day 0 of differentiation). A total of 1376 mRNA and 206 lncRNA were obtained. Finally, we processed the differential lncRNAs and mRNAs obtained in the two groups, and obtained 113 shared differential lncRNAs to further predict the targeted miRNA, a total of 815 lncRNA-miRNA pairs. The targeted mRNA was further predicted, and the grouped differential mRNAs were combined to obtain 64 differential mRNAs. In the end, we obtained 216 ceRNAs containing 26 lncRNAs, 27 miRNAs and 64 mRNAs. We found that the mRNAs in the ceRNA network were mainly enriched with 45 Gene Ontology (GO) terms, mainly including glucose homeostasis mechanism and insulin stimulation response. 69 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were mainly enriched. It mainly includes many pathways related to lipid metabolism such as Adenosine 5′-monophosphate (AMP)-activated protein kinase (AMPK), Rap1, cAMP, mitogen-activated protein kinase (MAPK), Ras, hypoxia inducible factor-1 (HIF-1), PI3K-Akt, insulin signaling and so on. In the end, we identified 216 ceRNA regulatory relationships related to obesity research. Our research provides a clearer direction for understanding the molecular mechanism of obesity, the screening and determination of drug targets biomarkers in the future. [ABSTRACT FROM AUTHOR]

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at dN / dS  = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination. [ABSTRACT FROM AUTHOR]

The soil organic carbon (C) cycle is primarily mediated by soil microorganisms and their genes that function in the C cycle (C‐cycle genes), both of which are strongly affected by land cover disturbance. However, the mechanism underlying microbially mediated soil C loss after conversion of primary natural broadleaf forests (BF) to plantation forests (PF) and secondary forests (SF) remains unknown. Here, we measured soil physicochemical properties and soil microbial community properties, and examined their linkages with microbial C‐cycle genes. Forest conversion dramatically decreased the richness of the soil fungal community but not of the bacterial community, and altered the composition of both communities. Analysis of C‐cycle genes revealed that the abundance of genes associated with C fixation, methane metabolism, and C degradation decreased by 51.3%, 57.9%, and 67.0%, respectively with the conversion of BF to PF; and by 6.3%, 4.1%, and 15.6%, respectively, with the conversion of BF to SF. The reductions in the abundance of C‐cycle genes, especially the reduction of hemicellulose‐ and lignin‐degradation genes, were primarily associated with the declines in the abundance of forest conversion‐sensitive microbes indexed by operational taxonomic units (fsOTUs, β = 0.41). fsOTUs were taxonomically diverse and included members frequently co‐occurring with numerous other microbes in the microbial communities, indicating that the manipulation of fsOTUs by forest management could improve soil fertility and soil C sequestration. Forest conversion‐induced shifts in fsOTUs abundance were associated with changes in soil potassium permanganate oxidizable organic carbon (PXC) concentration, dissolved organic carbon (DOC) concentration, and soil pH. Our results indicate that alterations in soil substrate supply (e.g., DOC and PXC) and soil pH induced by forest conversion may strongly shape fsOTUs structure and decrease the abundance of hemicellulose and lignin degradation genes, and consequently increase C loss. [ABSTRACT FROM AUTHOR]

Myeloproliferative syndromes (MPS) are hematologic malignancies due to the expansion of an abnormal hematopoietic stem cell. They include chronic myeloid leukemia (CML) and non‐CML MPS such as polycythemia vera, essential thrombocythemia and primary myelofibrosis. The latter are distinguished by somatic pathogenic variants affecting JAK2, CALR, or MPL genes. Apparent germline pathogenic variants have been reported in the general population. Here, we found that two gnomAD data‐sets report more homozygotes than expected for the JAK2 c.1849G > T(Val617Phe) variant. We propose that somatic gene conversion can explain the presence of those unexpected homozygotes in normal populations. Consistently, homozygous individuals are older than 65 years. We also found a lower‐than‐expected frequency of the JAK2 variant in younger individuals suggesting that somatic mutation can underlie its presence in (at least some) heterozygotes. Regarding pathogenic variants in MPL and CALR, they are also present in the gnomAD data‐sets explored. However, we cannot conclude that such seemingly germline variants are in fact somatic alterations. These results suggest that apparently normal individuals bearing MPS‐related variants can be subclinical/undiagnosed MPS cases of somatic origin. It would be interesting to assess the hematologic phenotype of such individuals and the presence of the relevant variants in other tissues. [ABSTRACT FROM AUTHOR]

Simple Summary: Harnessing the natural and induced diversity existing in plant genetic resources is fundamental for building future crops more sober in fertilizers, water, and pesticides that can cope with climate instability while yielding healthier and more nutritious products. The essence of plant breeding is to combine favorable traits in crossing parental varieties to select novel performing associations amongst the progenies. These associations are the product of recombination between the parental chromosomes occurring during meiosis, mainly by a reciprocal DNA exchange called Cross Over (CO). However, recombination does not occur randomly along the chromosomes, and COs are limited in number often hampering the desired associations of favorable traits. This review surveys the recent advances in methods for achieving a stimulation and/or a redistribution of meiotic COs along the parental chromosomes and targeting COs specifically at desired chromosomal sites. Reciprocal (cross-overs = COs) and non-reciprocal (gene conversion) DNA exchanges between the parental chromosomes (the homologs) during meiotic recombination are, together with mutation, the drivers for the evolution and adaptation of species. In plant breeding, recombination combines alleles from genetically diverse accessions to generate new haplotypes on which selection can act. In recent years, a spectacular progress has been accomplished in the understanding of the mechanisms underlying meiotic recombination in both model and crop plants as well as in the modulation of meiotic recombination using different strategies. The latter includes the stimulation and redistribution of COs by either modifying environmental conditions (e.g., T°), harnessing particular genomic situations (e.g., triploidy in Brassicaceae), or inactivating/over-expressing meiotic genes, notably some involved in the DNA double-strand break (DSB) repair pathways. These tools could be particularly useful for shuffling diversity in pre-breeding generations. Furthermore, thanks to the site-specific properties of genome editing technologies the targeting of meiotic recombination at specific chromosomal regions nowadays appears an attainable goal. Directing COs at desired chromosomal positions would allow breaking linkage situations existing between favorable and unfavorable alleles, the so-called linkage drag, and accelerate genetic gain. This review surveys the recent achievements in the manipulation of meiotic recombination in plants that could be integrated into breeding schemes to meet the challenges of deploying crops that are more resilient to climate instability, resistant to pathogens and pests, and sparing in their input requirements. [ABSTRACT FROM AUTHOR]

Porcine 17-hydroxysteroid dehydrogenase type 14 (HSD17B14) and FSH reporter (FSHR) genes play important roles in the metabolism of steroid hormones and the apoptosis of ovarian granulosa cells (GCs). Our bioinformatics analyses and the dual luciferase reporter assays indicated that porcine miR-20b and miR-31 target the 3′-UTR region of HSD17B14 gene, and miR-31 also targets the 3′-UTR region of FSHR gene. Overexpression of porcine HSD17B14 gene promoted the conversion from estradiol (E2) to estrone (E1) and increased the apoptosis of porcine GCs. Overexpression of miR-20b down-regulated the mRNA and protein expression level of HSD17B14 gene, decreased the concentration of progesterone (P4) and E1, increased E2, as well as reduced apoptosis of GCs. Moreover, overexpression of miR-31 also down-regulated the protein expression level of HSD17B14 gene, decreased the concentration of P4 and E1, and increased E2. However, miR-31 promoted apoptosis of GCs by targeting to the 3′-UTR of porcine FSHR gene. Taken together, we found that both porcine miR-20b and miR-31 target HSD17B14 gene, but miR-31 also targets FSHR gene to regulate the metabolism of steroid hormones and the apoptosis of porcine ovarian GCs. These findings expand the epigenetic regulatory mechanism of porcine miR-31 and miR-20b in ovarian GCs. • Both HSD17B14 and FSHR are the target genes of miR-31 in porcine GCs. • HSD17B14 is a target gene of miR-20b in porcine GCs. • Both miR-31 and miR-20b affected steroid hormone metabolism. • miR-31 promoted GC apoptosis, while miR-20b decreased GC apoptosis. [ABSTRACT FROM AUTHOR]

Background: Congenital adrenal hyperplasia (CAH) is a monogenic disorder caused by genetic diversity in the CYP21A2 gene, with 21-hydroxylase deficiency (21-OHD) as the most common type. Early sex assignment and early diagnosis of different genetic variations with a proper technique are important to reduce mortality and morbidity. Proper early sex identification reduces emotional, social, and psychological stress. Aim: Detection of a spectrum of aberrations in the CYP21A2 gene, including copy number variations, gene conversion, chimeric genes, and point variations. Methods: The CYP21A2 gene was screened using MLPA assay in 112 unrelated Egyptian children with 21-OHD CAH (33 males and 79 females). Results: In the studied group, 79.5% were diagnosed within the first month of life. 46.8% of the genetic females were misdiagnosed as males. Among the copy number variation results, large deletions in 15.4% and three types of chimeric genes in 9% (CH-1, CH-7, and CAH-X CH-1) were detected. Regarding gene dosage, one copy of CYP21A2 was found in 5 cases (4.5%), three copies were detected in 7 cases (6.3%), and one case (0.9%) showed four copies. Eight common genetic variants were identified, I2G, large deletions, large gene conversion (LGC), I172N, F306 + T, -113 SNP, 8bp Del, and exon 6 cluster (V237E and M239K) with an allelic frequency of 32.62%, 15.45%, 7.30%, 3.00%, 2.58%, 2.15%, 0.86%, and 0.86%, respectively. Conclusion: High prevalence of copy number variations highlights the added value of using MLPA in routine laboratory diagnosis of CAH patients. [ABSTRACT FROM AUTHOR]

Abstract Background As rice (Oryza sativa) is the staple food of more than half the world’s population, rice production contributes greatly to global food security. Rice blast caused by the fungus Magnaporthe oryzae (M. oryzae) is a devastating disease that affects rice yields and grain quality, resulting in substantial economic losses annually. Because the fungus evolves rapidly, the resistance conferred by most the single blast-resistance genes is broken after a few years of intensive agricultural use. Therefore, effective resistance breeding in rice requires continual enrichment of the reservoir of resistance genes, alleles, or QTLs. Seed banks represent a rich source of genetic diversity; however, they have not been extensively used to identify novel genes and alleles. Results We carried out a large-scale screen for novel blast-resistance alleles in 1883 rice varieties from major rice-producing areas across China. Of these, 361 varieties showed at least moderate resistance to natural infection by rice blast at rice blast nurseries in Enshi and Yichang, Hubei Province. We used sequence-based allele mining to amplify and sequence the allelic variants of the major rice blast-resistance genes at the Pi2/Pi9 locus of chromosome 6 from the 361 blast-resistant varieties, and the full-length coding region of this gene could be amplified from 107 varieties. Thirteen novel Pi9 alleles (named Pi9-Type1 to Pi9-Type13) were identified in these 107 varieties based on comparison to the Pi9 referenced sequence. Based on the sequencing results, the Pi2/Pi9 locus of the 107 varieties was divided into 15 genotypes (including three different genotypes of Pi9-Type5). Fifteen varieties, each representing one genotype, were evaluated for resistance to 34 M. oryzae isolates. The alleles from seven varieties with the highest resistance and widest resistance spectra were selected for transformation into the susceptible variety J23B to construct near-isogenic lines (NILs). These NILs showed resistance in a field test in Enshi and Yichang, indicating that the seven novel rice blast-resistance tandem-repeat regions at the Pi2/Pi9 locus of chromosome 6 could potentially serve as a genetic resource for molecular breeding of resistance to rice blast. Conclusions The thirteen novel Pi9 alleles identified in this study expand the list of available of blast-resistance alleles. Seven tandem-repeat regions of the Pi2/Pi9 locus from different donors were characterized as broad-spectrum rice blast-resistance fragments; these donors enrich the genetic resources available for rice blast-resistance breeding programs.

The Lyme disease spirochete Borrelia burgdorferi, encodes an elaborate antigenic variation system that promotes the ongoing variation of a major surface lipoprotein, VlsE. Changes in VlsE are continual and always one step ahead of the host acquired immune system, which requires 1–2 weeks to generate specific antibodies. By the time this happens, new VlsE variants have arisen that escape immunosurveillance, providing an avenue for persistent infection. This antigenic variation system is driven by segmental gene conversion events that transfer information from a series of silent cassettes (vls2-16) to the expression locus, vlsE. The molecular details of this process remain elusive. Recombinational switching at vlsE is RecA-independent and the only required factor identified to date is the RuvAB branch migrase. In this work we have used next generation long-read sequencing to analyze the effect of several DNA replication/recombination/repair gene disruptions on the frequency of gene conversions at vlsE and report a requirement for the mismatch repair protein MutL. Site directed mutagenesis of mutL suggests that the putative MutL endonuclease activity is required for recombinational switching at vlsE. This is the first report of an unexpected essential role for MutL in a bacterial recombination system and expands the known function of this protein as well as our knowledge of the details of the novel recombinational switching mechanism for vlsE variation.

Our understanding of the chronology of human evolution relies on the "molecular clock" provided by the steady accumulation of substitutions on an evolutionary lineage. Recent analyses of human pedigrees have called this understanding into question by revealing unexpectedly low germline mutation rates, which imply that substitutions accrue more slowly than previously believed. Translating mutation rates estimated from pedigrees into substitution rates is not as straightforward as it may seem, however. We dissect the steps involved, emphasizing that dating evolutionary events requires not "a mutation rate" but a precise characterization of how mutations accumulate in development in males and females-knowledge that remains elusive.

Events in primate evolution are often dated by assuming a constant rate of substitution per unit time, but the validity of this assumption remains unclear. Among mammals, it is well known that there exists substantial variation in yearly substitution rates. Such variation is to be expected from differences in life history traits, suggesting it should also be found among primates. Motivated by these considerations, we analyze whole genomes from 10 primate species, including Old World Monkeys (OWMs), New World Monkeys (NWMs), and apes, focusing on putatively neutral autosomal sites and controlling for possible effects of biased gene conversion and methylation at CpG sites. We find that substitution rates are up to 64% higher in lineages leading from the hominoid-NWM ancestor to NWMs than to apes. Within apes, rates are ∼2% higher in chimpanzees and ∼7% higher in the gorilla than in humans. Substitution types subject to biased gene conversion show no more variation among species than those not subject to it. Not all mutation types behave similarly, however; in particular, transitions at CpG sites exhibit a more clocklike behavior than do other types, presumably because of their nonreplicative origin. Thus, not only the total rate, but also the mutational spectrum, varies among primates. This finding suggests that events in primate evolution are most reliably dated using CpG transitions. Taking this approach, we estimate the human and chimpanzee divergence time is 12.1 million years,​ and the human and gorilla divergence time is 15.1 million years​.