Enrichment and physiological characterization of a novel comammox Nitrospira indicates ammonium inhibition of complete nitrification

The recent discovery of bacteria within the genus Nitrospira capable of complete ammonia oxidation (comammox) demonstrated that the sequential oxidation of ammonia to nitrate via nitrite can also be performed within a single bacterial cell. Although comammox Nitrospira exhibit a wide distribution in natural and engineered ecosystems, information on their physiological properties is scarce due to the limited number of cultured representatives. Additionally, most available genomic information is derived from metagenomic sequencing and high-quality genomes of Nitrospira in general are limited. In this study, we obtained a high (90%) enrichment of a novel comammox species, tentatively named “Candidatus Nitrospira kreftii”, and performed a detailed genomic and physiological characterization. The complete genome of “Ca. N. kreftii” allowed reconstruction of its basic metabolic traits. Similar to Nitrospira inopinata, the enrichment culture exhibited a very high ammonia affinity (K_{m(app)_NH3} ≈ 0.040 ± 0.01 µM), but a higher nitrite affinity (K_{m(app)_NO2}- = 12.5 ± 4.0 µM), indicating an adaptation to highly oligotrophic environments. Furthermore, we observed partial inhibition of ammonia oxidation at ammonium concentrations as low as 25 µM. This inhibition of “Ca. N. kreftii” indicates that differences in ammonium tolerance rather than affinity could potentially be a niche determining factor for different comammox Nitrospira.Subject terms: Bacterial genomics, Environmental microbiology, Bacterial physiology     

Selective Whole Genome Amplification for Resequencing Target Microbial Species from Complex Natural Samples

Population genomic analyses have demonstrated power to address major questions in evolutionary and molecular microbiology. Collecting populations of genomes is hindered in many microbial species by the absence of a cost effective and practical method to collect ample quantities of sufficiently pure genomic DNA for next-generation sequencing. Here we present a simple method to amplify genomes of a target microbial species present in a complex, natural sample. The selective whole genome amplification (SWGA) technique amplifies target genomes using nucleotide sequence motifs that are common in the target microbe genome, but rare in the background genomes, to prime the highly processive phi29 polymerase. SWGA thus selectively amplifies the target genome from samples in which it originally represented a minor fraction of the total DNA. The post-SWGA samples are enriched in target genomic DNA, which are ideal for population resequencing. We demonstrate the efficacy of SWGA using both laboratory-prepared mixtures of cultured microbes as well as a natural host–microbe association. Targeted amplification of Borrelia burgdorferi mixed with Escherichia coli at genome ratios of 1:2000 resulted in >10⁵-fold amplification of the target genomes with <6.7-fold amplification of the background. SWGA-treated genomic extracts from Wolbachia pipientis-infected Drosophila melanogaster resulted in up to 70% of high-throughput resequencing reads mapping to the W. pipientis genome. By contrast, 2–9% of sequencing reads were derived from W. pipientis without prior amplification. The SWGA technique results in high sequencing coverage at a fraction of the sequencing effort, thus allowing population genomic studies at affordable costs.     

Entire nucleotide sequences of Gossypium raimondii and G. arboreum mitochondrial genomes revealed A‐genome species as cytoplasmic donor of the allotetraploid species

• Cotton (Gossypium spp.) is commonly grouped into eight diploid genomic groups, designated A–G and K, and an allotetraploid genomic group, AD. Gossypium raimondii (D₅) and G. arboreum (A₂) are the putative contributors to the progenitor of G. hirsutum (AD₁), the economically important fibre‐producing cotton species.
• Mitochondrial DNA from week‐old etiolated seedlings was extracted from isolated organelles using discontinuous sucrose density gradient method. Mitochondrial genomes were sequenced, assembled, annotated and analysed in orderly.
• Gossypium raimondii (D₅) and G. arboreum (A₂) mitochondrial genomes were provided in this study. The mitochondrial genomes of two diploid species harboured circular genome of 643,914 bp (D₅) and 687,482 bp (A₂), respectively. They differ in size and number of repeat sequences, both contain illuminating triplicate sequences with 7317 and 10,246 bp, respectively, demonstrating dynamic difference and rearranged genome organisations. Comparing the D₅ and A₂ mitogenomes with mitogenomes of tetraploid Gossypium species (AD₁, G. hirsutum; AD₂, G. barbadense), a shared 11 kbp fragment loss was detected in allotetraploid species, three regions shared by G. arboreum (A₂), G. hirsutum (AD₁) and G. barbadense (AD₂), while eight regions were specific to G. raimondii (D₅). The presence/absence variations and gene‐based phylogeny supported that A‐genome is a cytoplasmic donor to the progenitor of allotetraploid species G. hirsutum and G. barbadense.
• The results present structure variations and phylogeny of Gossypium mitochondrial genome evolution.

     

Incipient Genome Differentiation in Gossypium. I. Chromosomes 14, 15, 16, 19 and 20 Assessed in G. HIRSUTUM, G. RAIMONDII and G. LOBATUM by Means of Seven a-D Translocations

The genus Gossypium is favorable for study of genome divergence at several levels. Early stages of divergence have been studied among four D genomes by comparing chiasma frequencies (reciprocal exchanges) between pairs of genomes and between individual counterpart chromosomes marked by heterozygous translocations. D₅ (G. raimondii) shows barely detectable differentiation from from D_h (G. hirsutum), whereas D₇ (G. lobatum) is considerably less closely related to D_h than is D₅. Fragmentary data suggest that D_2–2 (G. harknessii) falls between D₅ and D₇ in its relationship to D_h. Since chiasma frequencies in individual chromosomes and marked regions exhibit the same order of relationships as their corresponding whole genomes, it is concluded that the genome differentiation is generalized (i.e., nucleus-wide) rather than localized in specific chromosomes or chromosome regions. Estimates of relationships based on reciprocal exchange frequencies agree with those based upon preferential synapsis in allohexaploids reported previously. Since preferential synapsis and reciprocal exchange frequencies reveal the same order of relationships, it is concluded that to some extent they reflect common underlying changes in chromosome properties, despite recent evidence that synapsis and crossing over are under independent genetic control.     

Activation of the β-carbonic anhydrase from the protozoan pathogen Trichomonas vaginalis with amines and amino acids

We report the first activation study of the β-class carbonic anhydrase (CA, EC 4.2.1.1) encoded in the genome of the protozoan pathogen Trichomonas vaginalis, TvaCA1. Among 24 amino acid and amine activators investigated, derivatives incorporating a second carboxylic moiety, such as L-Asp, L- and D-Glu, were devoid of activating effects up to concentrations of 50 µM within the assay system, whereas the corresponding compounds with a CONH₂ moiety, i.e. L-Gln and L-Asn showed modest activating effects, with activation constants in the range of 26.9 − 32.5 µM. Moderate activation was observed with L- and D-DOPA, histamine, dopamine, serotonin, (2-Aminoethyl)pyridine/piperazine and morpholine (K_A‘s ranging between 8.3 and 14.5 µM), while the best activators were L-and D-Trp, L-and D-Tyr and 4-amino-Phe, which showed K_A‘s ranging between 3.0 and 5.1 µM. Understanding in detail the activation mechanism of β-CAs may be relevant for the design of enzyme activity modulators with potential clinical significance.     

The Oxytricha trifallax Macronuclear Genome: A Complex Eukaryotic Genome with 16,000 Tiny Chromosomes

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor “silent” germline micronuclear genome by a process of “unscrambling” and fragmentation. The tiny macronuclear “nanochromosomes” typically encode single, protein-coding genes (a small portion, 10%, encode 2–8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.     

Genome sequencing of a single cell of the widely distributed marine subsurface Dehalococcoidia,phylum Chloroflexi

Bacteria of the class Dehalococcoidia (DEH), phylum Chloroflexi, are widely distributed in the marine subsurface, yet metabolic properties of the many uncultivated lineages are completely unknown. This study therefore analysed genomic content from a single DEH cell designated ‘DEH-J10'' obtained from the sediments of Aarhus Bay, Denmark. Real-time PCR showed the DEH-J10 phylotype was abundant in upper sediments but was absent below 160 cm below sea floor. A 1.44 Mbp assembly was obtained and was estimated to represent up to 60.8% of the full genome. The predicted genome is much larger than genomes of cultivated DEH and appears to confer metabolic versatility. Numerous genes encoding enzymes of core and auxiliary beta-oxidation pathways were identified, suggesting that this organism is capable of oxidising various fatty acids and/or structurally related substrates. Additional substrate versatility was indicated by genes, which may enable the bacterium to oxidise aromatic compounds. Genes encoding enzymes of the reductive acetyl-CoA pathway were identified, which may also enable the fixation of CO₂ or oxidation of organics completely to CO₂. Genes encoding a putative dimethylsulphoxide reductase were the only evidence for a respiratory terminal reductase. No evidence for reductive dehalogenase genes was found. Genetic evidence also suggests that the organism could synthesise ATP by converting acetyl-CoA to acetate by substrate-level phosphorylation. Other encoded enzymes putatively conferring marine adaptations such as salt tolerance and organo-sulphate sulfohydrolysis were identified. Together, these analyses provide the first insights into the potential metabolic traits that may enable members of the DEH to occupy an ecological niche in marine sediments.     

Chromosome-scale haplotype-phased genome assemblies of the male and female lines of wild asparagus (Asparagus kiusianus), a dioecious plant species

Asparagus kiusianus is a disease-resistant dioecious plant species and a wild relative of garden asparagus (Asparagus officinalis). To enhance A. kiusianus genomic resources, advance plant science, and facilitate asparagus breeding, we determined the genome sequences of the male and female lines of A. kiusianus. Genome sequence reads obtained with a linked-read technology were assembled into four haplotype-phased contig sequences (∼1.6 Gb each) for the male and female lines. The contig sequences were aligned onto the chromosome sequences of garden asparagus to construct pseudomolecule sequences. Approximately 55,000 potential protein-encoding genes were predicted in each genome assembly, and ∼70% of the genome sequence was annotated as repetitive. Comparative analysis of the genomes of the two species revealed structural and sequence variants between the two species as well as between the male and female lines of each species. Genes with high sequence similarity with the male-specific sex determinant gene in A. officinalis, MSE1/AoMYB35/AspTDF1, were presented in the genomes of the male line but absent from the female genome assemblies. Overall, the genome sequence assemblies, gene sequences, and structural and sequence variants determined in this study will reveal the genetic mechanisms underlying sexual differentiation in plants, and will accelerate disease-resistance breeding in garden asparagus.     

Analysis of Transposable Elements in the Genome of Asparagus officinalis from High Coverage Sequence Data

Asparagus officinalis is an economically and nutritionally important vegetable crop that is widely cultivated and is used as a model dioecious species to study plant sex determination and sex chromosome evolution. To improve our understanding of its genome composition, especially with respect to transposable elements (TEs), which make up the majority of the genome, we performed Illumina HiSeq2000 sequencing of both male and female asparagus genomes followed by bioinformatics analysis. We generated 17 Gb of sequence (12×coverage) and assembled them into 163,406 scaffolds with a total cumulated length of 400 Mbp, which represent about 30% of asparagus genome. Overall, TEs masked about 53% of the A. officinalis assembly. Majority of the identified TEs belonged to LTR retrotransposons, which constitute about 28% of genomic DNA, with Ty1/copia elements being more diverse and accumulated to higher copy numbers than Ty3/gypsy. Compared with LTR retrotransposons, non-LTR retrotransposons and DNA transposons were relatively rare. In addition, comparison of the abundance of the TE groups between male and female genomes showed that the overall TE composition was highly similar, with only slight differences in the abundance of several TE groups, which is consistent with the relatively recent origin of asparagus sex chromosomes. This study greatly improves our knowledge of the repetitive sequence construction of asparagus, which facilitates the identification of TEs responsible for the early evolution of plant sex chromosomes and is helpful for further studies on this dioecious plant.     

Incipient Genome Differentiation in Gossypium. III. Comparison of Chromosomes of G. HIRSUTUM and Asiatic Diploids Using Heterozygous Translocations

Hybrids between upland cotton (G. hirsutum, genome constitution 2A_hD_h) and either A-genome or D-genome diploid species exhibit 26 paired and 13 unpaired chromosomes at metaphase I. The A_h and D_h genomes are therefore considered homoeologous with those of the respective diploids. Previous studies, nevertheless, revealed a low level of ("incipient") differentiation between D_h and various diploid D genomes. The diploid A genomes have been regarded as more closely homologous to A_h on the basis of low preferential pairing and autotetraploid segregation ratios in allohexaploids.—The present study addressed the following questions: Are the diploid A genomes differentiated from A_h in meiotic homology? If so, is the differentiation manifested equally by all 13 chromosomes or is it localized in certain chromosomes?—Three diploid A-genome lines representing G. herbaceum and G. arboreum were hybridized by in ovulo culture of embryos (1) with a standard line of G. hirsutum, which differs from G. herbaceum by two and from G. arboreum by three naturally occurring reciprocal translocations involving chromosomes 1–5, and (2) with six lines homozygous for experimental translocations involving chromosomes 6, 7, 10, 11, 12 and 13. Chiasma frequencies in hybrids were compared with those in appropriate G. hirsutum controls. In every comparison overall chiasma frequencies were slightly lower in the hybrids. Therefore A_h appears to be differentiated from the diploid A genomes. No localized differentiation was detected in chromosomes marked by experimental translocations. The differentiation may be localized mainly in chromosomes 4 and 5.     

Chromosomal-level genome assembly of the orchid tree Bauhinia variegata (Leguminosae; Cercidoideae) supports the allotetraploid origin hypothesis of Bauhinia

Cercidoideae, one of the six subfamilies of Leguminosae, contains one genus Cercis with its chromosome number 2n = 14 and all other genera with 2n = 28. An allotetraploid origin hypothesis for the common ancestor of non-Cercis genera in this subfamily has been proposed; however, no chromosome-level genomes from Cercidoideae have been available to test this hypothesis. Here, we conducted a chromosome-level genome assembly of Bauhinia variegata to test this hypothesis. The assembled genome is 326.4 Mb with the scaffold N50 of 22.1 Mb and contains 37,996 protein-coding genes. The Ks distribution between gene pairs in the syntenic regions indicates two whole-genome duplications (WGDs): one is B. variegata-specific, and the other is shared among core eudicots. Although Ks between gene pairs generated by the recent WGD in Bauhinia is greater than that between Bauhinia and Cercis, the WGD was not detected in Cercis, which can be explained by an accelerated evolutionary rate in Bauhinia after divergence from Cercis. Ks distribution and phylogenetic analysis for gene pairs generated by the recent WGD in Bauhinia and their corresponding orthologs in Cercis support the allopolyploidy origin hypothesis of Bauhinia. The genome of B. variegata also provides a genomic resource for dissecting genetic basis of its ornamental traits.     

Chromosome-scale genome assembly of the transformation-amenable common wheat cultivar ‘Fielder’

We have established a high-quality, chromosome-level genome assembly for the hexaploid common wheat cultivar ‘Fielder’, an American, soft, white, pastry-type wheat released in 1974 and known for its amenability to Agrobacterium tumefaciens-mediated transformation and genome editing. Accurate, long-read sequences were obtained using PacBio circular consensus sequencing with the HiFi approach. Sequence reads from 16 SMRT cells assembled using the hifiasm assembler produced assemblies with N50 greater than 20 Mb. We used the Omni-C chromosome conformation capture technique to order contigs into chromosome-level assemblies, resulting in 21 pseudomolecules with a cumulative size of 14.7 and 0.3 Gb of unanchored contigs. Mapping of published short reads from a transgenic wheat plant with an edited seed-dormancy gene, TaQsd1, identified four positions of transgene insertion into wheat chromosomes. Detection of guide RNA sequences in pseudomolecules provided candidates for off-target mutation induction. These results demonstrate the efficiency of chromosome-scale assembly using PacBio HiFi reads and their application in wheat genome-editing studies.     

Calcium Regulates Molecular Interactions of Otoferlin with Soluble NSF Attachment Protein Receptor (SNARE) Proteins Required for Hair Cell Exocytosis

Mutations in otoferlin, a C2 domain-containing ferlin family protein, cause non-syndromic hearing loss in humans (DFNB9 deafness). Furthermore, transmitter secretion of cochlear inner hair cells is compromised in mice lacking otoferlin. In the present study, we show that the C2F domain of otoferlin directly binds calcium (K_D = 267 μm) with diminished binding in a pachanga (D1767G) C2F mouse mutation. Calcium was found to differentially regulate binding of otoferlin C2 domains to target SNARE (t-SNARE) proteins and phospholipids. C2D–F domains interact with the syntaxin-1 t-SNARE motif with maximum binding within the range of 20–50 μm Ca²⁺. At 20 μm Ca²⁺, the dissociation rate was substantially lower, indicating increased binding (K_D = ∼10⁻⁹) compared with 0 μm Ca²⁺ (K_D = ∼10⁻⁸), suggesting a calcium-mediated stabilization of the C2 domain·t-SNARE complex. C2A and C2B interactions with t-SNAREs were insensitive to calcium. The C2F domain directly binds the t-SNARE SNAP-25 maximally at 100 μm and with reduction at 0 μm Ca²⁺, a pattern repeated for C2F domain interactions with phosphatidylinositol 4,5-bisphosphate. In contrast, C2F did not bind the vesicle SNARE protein synaptobrevin-1 (VAMP-1). Moreover, an antibody targeting otoferlin immunoprecipitated syntaxin-1 and SNAP-25 but not synaptobrevin-1. As opposed to an increase in binding with increased calcium, interactions between otoferlin C2F domain and intramolecular C2 domains occurred in the absence of calcium, consistent with intra-C2 domain interactions forming a “closed” tertiary structure at low calcium that “opens” as calcium increases. These results suggest a direct role for otoferlin in exocytosis and modulation of calcium-dependent membrane fusion.     

Utilization of Exogenous Inorganic Carbon Species in Photosynthesis by Chlorella pyrenoidosa

The nature of the inorganic carbon utilized during photosynthesis by Chlorella pyrenoidosa was investigated using three experimental techniques (open gas analysis system with “artificial leaf” or “aqueous” chambers and O₂ electrode system) to measure carbon assimilation. Photosynthesis was studied as a function of pH and CO₂ concentration. The CO₂ concentration was inadequate to meet the requirements of photosynthesis only when HCO₃⁻ was added at high pH. Under all other conditions, the low and constant K_m (CO₂), in contrast to the highly variable K_m (HCO₃⁻), suggested that CO₂ was the major species utilized.     

AtnMat2, a nuclear-encoded maturase required for splicing of group-II introns in Arabidopsis mitochondria

Mitochondria (mt) in plants house about 20 group-II introns, which lie within protein-coding genes required in both organellar genome expression and respiration activities. While in nonplant systems the splicing of group-II introns is mediated by proteins encoded within the introns themselves (known as “maturases”), only a single maturase ORF (matR) has retained in the mitochondrial genomes in plants; however, its putative role(s) in the splicing of organellar introns is yet to be established. Clues to other proteins are scarce, but these are likely encoded within the nucleus as there are no obvious candidates among the remaining ORFs within the mtDNA. Intriguingly, higher plants genomes contain four maturase-related genes, which exist in the nucleus as self-standing ORFs, out of the context of their evolutionary-related group-II introns “hosts.” These are all predicted to reside within mitochondria and may therefore act “in-trans” in the splicing of organellar-encoded introns. Here, we analyzed the intracellular locations of the four nuclear-encoded maturases in Arabidopsis and established the roles of one of these genes, At5g46920 (AtnMat2), in the splicing of several mitochondrial introns, including the single intron within cox2, nad1 intron2, and nad7 intron2.     

Alu recombination-mediated structural deletions in the chimpanzee genome

With more than 1.2 million copies, Alu elements are one of the most important sources of structural variation in primate genomes. Here, we compare the chimpanzee and human genomes to determine the extent of Alu recombination-mediated deletion (ARMD) in the chimpanzee genome since the divergence of the chimpanzee and human lineages (~6 million y ago). Combining computational data analysis and experimental verification, we have identified 663 chimpanzee lineage-specific deletions (involving a total of ~771 kb of genomic sequence) attributable to this process. The ARMD events essentially counteract the genomic expansion caused by chimpanzee-specific Alu inserts. The RefSeq databases indicate that 13 exons in six genes, annotated as either demonstrably or putatively functional in the human genome, and 299 intronic regions have been deleted through ARMDs in the chimpanzee lineage. Therefore, our data suggest that this process may contribute to the genomic and phenotypic diversity between chimpanzees and humans. In addition, we found four independent ARMD events at orthologous loci in the gorilla or orangutan genomes. This suggests that human orthologs of loci at which ARMD events have already occurred in other nonhuman primate genomes may be “at-risk” motifs for future deletions, which may subsequently contribute to human lineage-specific genetic rearrangements and disorders.