Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nucleotide”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

A two-component protein condensate of the EGFR cytoplasmic tail and Grb2 regulates Ras activation by SOS at the membrane

We reconstitute a phosphotyrosine-mediated protein condensation phase transition of the ~200 residue cytoplasmic tail of the epidermal growth factor receptor (EGFR) and the adaptor protein, Grb2, on a membrane surface. The phase transition depends on phosphorylation of the EGFR tail, which recruits Grb2, and crosslinking through a Grb2-Grb2 binding interface. The Grb2 Y160 residue plays a structurally critical role in the Grb2-Grb2 interaction, and phosphorylation or mutation of Y160 prevents EGFR:Grb2 condensation. By extending the reconstitution experiment to include the guanine nucleotide exchange factor, SOS, and its substrate Ras, we further find that the condensation state of the EGFR tail controls the ability of SOS, recruited via Grb2, to activate Ras. These results identify an EGFR:Grb2 protein condensation phase transition as a regulator of signal propagation from EGFR to the MAPK pathway.

59 BASIC BIOLOGICAL SCIENCES↗

Transcription factor bHLH121 regulates root cortical aerenchyma formation in maize

Root anatomical phenotypes present a promising yet underexploited avenue to deliver major improvements in yield and climate resilience of crops by improving water and nutrient uptake. For instance, the formation of root cortical aerenchyma (RCA) significantly increases soil exploration and resource capture by reducing the metabolic costs of root tissue. A key bottleneck in studying such phenotypes has been the lack of robust high-throughput anatomical phenotyping platforms. We exploited a phenotyping approach based on laser ablation tomography, termed Anatomics , to quantify variation in RCA formation of 436 diverse maize lines in the field. Results revealed a significant and heritable variation for RCA formation. Genome-wide association studies identified a single-nucleotide polymorphism mapping to a root cortex-expressed gene-encoding transcription factor bHLH121. Functional studies identified that the bHLH121 Mu transposon mutant line and CRISPR/Cas9 loss-of-function mutant line showed reduced RCA formation, whereas an overexpression line exhibited significantly greater RCA formation when compared to the wild-type line. Characterization of these lines under suboptimal water and nitrogen availability in multiple soil environments revealed that bHLH121 is required for RCA formation developmentally as well as under studied abiotic stress. Overall functional validation of the bHLH121 gene’s importance in RCA formation provides a functional marker to select varieties with improved soil exploration and thus yield under suboptimal conditions.

59 BASIC BIOLOGICAL SCIENCES↗

ADP enhances the allosteric activation of eukaryotic elongation factor 2 kinase by calmodulin

Protein translation, one of the most energy-consumptive processes in a eukaryotic cell, requires robust regulation, especially under energy-deprived conditions. A critical component of this regulation is the suppression of translational elongation through reduced ribosome association of the GTPase eukaryotic elongation factor 2 (eEF-2) resulting from its specific phosphorylation by the calmodulin (CaM)-activated α–kinase eEF-2 kinase (eEF-2K). It has been suggested that the eEF-2K response to reduced cellular energy levels is indirect and mediated by the universal energy sensor AMP-activated protein kinase (AMPK) through direct stimulatory phosphorylation and/or downregulation of the eEF-2K-inhibitory nutrient-sensing mTOR pathway. Here, we provide structural, biochemical, and cell-biological evidence of a direct energy-sensing role of eEF-2K through its stimulation by ADP. A crystal structure of the nucleotide-bound complex between CaM and the functional core of eEF-2K phosphorylated at its primary stimulatory site (T348) reveals ADP bound at a unique pocket located on the face opposite that housing the kinase active site. Within this basic pocket (BP), created at the CaM/eEF-2K interface upon complex formation, ADP is stabilized through numerous interactions with both interacting partners. Biochemical analyses using wild-type eEF-2K and specific BP mutants indicate that ADP stabilizes CaM within the active complex, increasing the sensitivity of the kinase to CaM. Induction of energy stress through glycolysis inhibition results in significantly reduced enhancement of phosphorylated eEF-2 levels in cells expressing ADP-binding compromised BP mutants compared to cells expressing wild-type eEF-2K. These results suggest a direct energy-sensing role for eEF-2K through its cooperative interaction with CaM and ADP.

59 BASIC BIOLOGICAL SCIENCES↗

30S subunit recognition and G1405 modification by the aminoglycoside-resistance 16S ribosomal RNA methyltransferase RmtC

Acquired ribosomal RNA (rRNA) methylation has emerged as a significant mechanism of aminoglycoside resistance in pathogenic bacterial infections. Modification of a single nucleotide in the ribosome decoding center by the aminoglycoside-resistance 16S rRNA (m 7 G1405) methyltransferases effectively blocks the action of all 4,6-deoxystreptamine ring-containing aminoglycosides, including the latest generation of drugs. To define the molecular basis of 30S subunit recognition and G1405 modification by these enzymes, we used a S-adenosyl-L-methionine analog to trap the complex in a postcatalytic state to enable determination of a global 3.0 Å cryo-electron microscopy structure of the m 7 G1405 methyltransferase RmtC bound to the mature Escherichia coli 30S ribosomal subunit. This structure, together with functional analyses of RmtC variants, identifies the RmtC N-terminal domain as critical for recognition and docking of the enzyme on a conserved 16S rRNA tertiary surface adjacent to G1405 in 16S rRNA helix 44 (h44). To access the G1405 N7 position for modification, a collection of residues across one surface of RmtC, including a loop that undergoes a disorder-to order transition upon 30S subunit binding, induces significant distortion of h44. This distortion flips G1405 into the enzyme active site where it is positioned for modification by two almost universally conserved RmtC residues. Further, these studies expand our understanding of ribosome recognition by rRNA modification enzymes and present a more complete structural basis for future development of strategies to inhibit m 7 G1405 modification to resensitize bacterial pathogens to aminoglycosides.

59 BASIC BIOLOGICAL SCIENCES↗

Structural and physical features that distinguish tumor-controlling from inactive cancer neoepitopes

Neoepitopes arising from amino acid substitutions due to single nucleotide polymorphisms are targets of T cell immune responses to cancer and are of significant interest in the development of cancer vaccines. However, understanding the characteristics of rare protective neoepitopes that truly control tumor growth has been a challenge, due to their scarcity as well as the challenge of verifying true, neoepitope-dependent tumor control in humans. Taking advantage of recent work in mouse models that circumvented these challenges, here, we compared the structural and physical properties of neoepitopes that range from fully protective to immunologically inactive. As neoepitopes are derived from self-peptides that can induce immune tolerance, we studied not only how the various neoepitopes differ from each other but also from their wild-type counterparts. We identified multiple features associated with protection, including features that describe how neoepitopes differ from self as well as features associated with recognition by diverse T cell receptor repertoires. In conclusion, we demonstrate both the promise and limitations of neoepitope structural analysis and predictive modeling and illustrate important aspects that can be incorporated into neoepitope prediction pipelines.

60 APPLIED LIFE SCIENCES↗

Extraordinary preservation of gene collinearity over three hundred million years revealed in homosporous lycophytes

Homosporous lycophytes (Lycopodiaceae) are a deeply diverged lineage in the plant tree of life, having split from heterosporous lycophytes ( Selaginella and Isoetes ) ~400 Mya. Compared to the heterosporous lineage, Lycopodiaceae has markedly larger genome sizes and remains the last major plant clade for which no chromosome-level assembly has been available. Here, we present chromosomal genome assemblies for two homosporous lycophyte species, the allotetraploid Huperzia asiatica and the diploid Diphasiastrum complanatum . Remarkably, despite that the two species diverged ~350 Mya, around 30% of the genes are still in syntenic blocks. Furthermore, both genomes had undergone independent whole genome duplications, and the resulting intragenomic syntenies have likewise been preserved relatively well. Such slow genome evolution over deep time is in stark contrast to heterosporous lycophytes and is correlated with a decelerated rate of nucleotide substitution. Together, the genomes of H. asiatica and D. complanatum not only fill a crucial gap in the plant genomic landscape but also highlight a potentially meaningful genomic contrast between homosporous and heterosporous species.

59 BASIC BIOLOGICAL SCIENCES↗

Structure of saguaro cactus virus 3′ translational enhancer mimics 5′ cap for eIF4E binding

The genomes of several plant viruses contain RNA structures at their 3' ends called cap-independent translation enhancers (CITEs) that bind the host protein factors such as mRNA 5' cap-binding protein eIF4E for promoting cap-independent genome translation. However, the structural basis of such 5' cap-binding protein recognition by the uncapped RNA remains largely unknown. Here, we have determined the crystal structure of a 3' CITE, panicum mosaic virus-like translation enhancer (PTE) from the saguaro cactus virus (SCV), using a Fab crystallization chaperone. The PTE RNA folds into a three-way junction architecture with a pseudoknot between the purine-rich R domain and pyrimidine-rich Y domain, which organizes the overall structure to protrude out a specific guanine nucleotide, G18, from the R domain that comprises a major interaction site for the eIF4E binding. The superimposable crystal structures of the wild-type, G18A, G18C, and G18U mutants suggest that the PTE scaffold is preorganized with the flipped-out G18 ready to dock into the eIF4E 5' cap-binding pocket. The binding studies with wheat and human eIF4Es using gel electrophoresis and isothermal titration calorimetry, and molecular docking computation for the PTE–eIF4E complex demonstrated that the PTE structure essentially mimics the mRNA 5' cap for eIF4E binding. Such 5' cap mimicry by the uncapped and structured viral RNA highlights how viruses can exploit RNA structures to mimic the host protein-binding partners and bypass the canonical mechanisms for their genome translation, providing opportunities for a better understanding of virus-host interactions and non-canonical translation mechanisms found in many pathogenic RNA viruses.

3' cap-independent translation enhancers↗

Nanoscale dynamics of Dynamin 1 helices reveals squeeze-twist deformation mode critical for membrane fission

Dynamin 1 (Dyn1) GTPase, a principal driver of membrane fission during synaptic endocytosis, self-assembles into short mechanoactive helices cleaving the necks of endocytic vesicles. While structural information about Dyn1 helix is abundant, little is known about the nanoscale dynamics of the helical scaffolding at the moment of fission, complicating mechanistic understanding of Dyn1 action. To address the role of the helix dynamics in fission, we used High-Speed Atomic Force Microscopy (HS-AFM) and fluorescence microscopy to track and compare the spatiotemporal characteristics of the helices formed by wild-type Dyn1 and its K44A mutant impaired in GTP hydrolysis on minimal lipid membrane templates. In the absence of nucleotide, membrane-bound WT Dyn1 and K44A Dyn1 self-assembled into tubular protein scaffolding of similar diameter encaging the lipid bilayer. In both cases, the GTP addition caused scaffold constriction coupled with formation of 20 to 30 nm nanogaps in the protein coverage. While both proteins reached scaffold diameters characteristic for membrane superconstriction causing fission, the fission was detected only with WT Dyn1. We associated the fission activity with the dynamic evolution of the nanogaps: K44A Dyn1 gaps were static, while WT Dyn1 gaps actively evolved via repetitive nonaxisymmetric constriction-bending deformations caused by localized GTP hydrolysis. Modeling of the deformations implicated filament twist as an additional deformation mode which combines with superconstriction to facilitate membrane fission. Our results thus show that the dynamics of the Dyn1 helical scaffold goes beyond radial constriction and involves nonaxisymmetric deformations, where filament twist emerges as a critical driver of membrane fission.

59 BASIC BIOLOGICAL SCIENCES↗

Population structure limits the use of genomic data for predicting phenotypes and managing genetic resources in forest trees

There is overwhelming evidence that forest trees are locally adapted to climate. Thus, genecological models based on population phenotypes have been used to measure local adaptation, infer genetic maladaptation to climate, and guide assisted migration. However, instead of phenotypes, there is increasing interest in using genomic data for gene resource management. We used whole-genome resequencing and common-garden experiments to understand the genetic architecture of adaptive traits in black cottonwood. We studied the potential of using genome-wide association studies (GWAS) and genomic prediction to detect causal loci, identify climate-adapted phenotypes, and inform gene resource management. We analyzed population structure by partitioning phenotypic and genomic (single-nucleotide polymorphism) variation among 840 genotypes collected from 91 stands along 16 rivers. Most phenotypic variation (60 to 81%) occurred among populations and was strongly associated with climate. Population phenotypes were predicted well using genomic data (e.g., predictive abilityr> 0.9) but almost as well using climate or geography (r> 0.8). In contrast, genomic prediction within populations was poor (r< 0.2). We identified many GWAS associations among populations, but most appeared to be spurious based on pooled within-population analyses. Hierarchical partitioning of linkage disequilibrium and haplotype sharing suggested that within-population genomic prediction and GWAS were poor because allele frequencies of causal loci and linked markers differed among populations. Given the urgent need to conserve natural populations and ecosystems, our results suggest that climate variables alone can be used to predict population phenotypes, delineate seed zones and deployment zones, and guide assisted migration.

Science & Technology - Other Topics↗

Comparison of human poly-N-acetyl-lactosamine synthase structure with GT-A fold glycosyltransferases supports a modular assembly of catalytic subsites

Poly-N-acetyl-lactosamine (poly-LacNAc) structures are composed of repeating [-Galβ(1,4)-GlcNAcβ(1,3)-]n glycan extensions. They are found on both N- and O-glycoproteins and glycolipids and play an important role in development, immune function, and human disease. The majority of mammalian poly-LacNAc is synthesized by the alternating iterative action of β1,3-N-acetylglucosaminyltransferase 2 (B3GNT2) and β1,4-galactosyltransferases. B3GNT2 is in the largest mammalian glycosyltransferase family, GT31, but little is known about the structure, substrate recognition, or catalysis by family members. Here we report the structures of human B3GNT2 in complex with UDP:Mg 2+ and in complex with both UDP:Mg 2+ and a glycan acceptor, lacto-N-neotetraose. The B3GNT2 structure conserves the GT-A fold and the DxD motif that coordinates a Mg 2+ ion for binding the UDP-GlcNAc sugar donor. The acceptor complex shows interactions with only the terminal Galβ(1,4)-GlcNAcβ(1,3)- disaccharide unit, which likely explains the specificity for both N- and O-glycan acceptors. Modeling of the UDP-GlcNAc donor supports a direct displacement inverting catalytic mechanism. Comparative structural analysis indicates that nucleotide sugar donors for GT-A fold glycosyltransferases bind in similar positions and conformations without conserving interacting residues, even for enzymes that use the same donor substrate. In contrast, the B3GNT2 acceptor binding site is consistent with prior models suggesting that the evolution of acceptor specificity involves loops inserted into the stable GT-A fold. These observations support the hypothesis that GT-A fold glycosyltransferases employ coevolving donor, acceptor, and catalytic subsite modules as templates to achieve the complex diversity of glycan linkages in biological systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Relaxed tarantula skeletal muscle has two ATP energy-saving mechanisms

Myosin molecules in the relaxed thick filaments of striated muscle have a helical arrangement in which the heads of each molecule interact with each other, forming the interacting-heads motif (IHM). In relaxed mammalian skeletal muscle, this helical ordering occurs only at temperatures >20°C and is disrupted when temperature is decreased. Recent x-ray diffraction studies of live tarantula skeletal muscle have suggested that the two myosin heads of the IHM (blocked heads [BHs] and free heads [FHs]) have very different roles and dynamics during contraction. Here, we explore temperature-induced changes in the BHs and FHs in relaxed tarantula skeletal muscle. We find a change with decreasing temperature that is similar to that in mammals, while increasing temperature induces a different behavior in the heads. At 22.5°C, the BHs and FHs containing ADP.Pi are fully helically organized, but they become progressively disordered as temperature is lowered or raised. Our interpretation suggests that at low temperature, while the BHs remain ordered the FHs become disordered due to transition of the heads to a straight conformation containing Mg.ATP. Above 27.5°C, the nucleotide remains as ADP.Pi, but while BHs remain ordered, half of the FHs become progressively disordered, released semipermanently at a midway distance to the thin filaments while the remaining FHs are docked as swaying heads. We propose a thermosensing mechanism for tarantula skeletal muscle to explain these changes. Our results suggest that tarantula skeletal muscle thick filaments, in addition to having a superrelaxation–based ATP energy-saving mechanism in the range of 8.5–40°C, also exhibit energy saving at lower temperatures (<22.5°C), similar to the proposed refractory state in mammals.

59 BASIC BIOLOGICAL SCIENCES↗

Pan-Coronavirus CRISPR-CasRx Effector System Significantly Reduces Viable Titer in HCoV-OC43, HCoV-229E, and SARS-CoV-2

CRISPR-based technology has become widely used as an antiviral strategy, including as a broad-spectrum human coronavirus (HCoV) therapeutic. In this work, we have designed a CRISPR-CasRx effector system with guide RNAs (gRNAs) that are cross-reactive among several HCoV species. We tested the efficacy of this pan-coronavirus effector system by evaluating the reduction in viral viability associated with different CRISPR targets in HCoV-OC43, HCoV-229E, and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We determined that several CRISPR targets significantly reduce viral titer, despite the presence of single nucleotide polymorphisms in the gRNA when compared with a non-targeting, negative control gRNA. CRISPR targets reduced viral titer between 85% and >99% in HCoV-OC43, between 78% and >99% in HCoV-229E, and between 70% and 94% in SARS-CoV-2 when compared with an untreated virus control. These data establish a proof-of-concept for a pan-coronavirus CRISPR effector system that is capable of reducing viable virus in both Risk Group 2 and Risk Group 3 HCoV pathogens.

60 APPLIED LIFE SCIENCES↗

Divergent selection and climate adaptation fuel genomic differentiation between sister species of Sphagnum (peat moss)

Abstract Background and Aims New plant species can evolve through the reinforcement of reproductive isolation via local adaptation along habitat gradients. Peat mosses (Sphagnaceae) are an emerging model system for the study of evolutionary genomics and have well-documented niche differentiation among species. Recent molecular studies have demonstrated that the globally distributed species Sphagnum magellanicum is a complex of morphologically cryptic lineages that are phylogenetically and ecologically distinct. Here, we describe the architecture of genomic differentiation between two sister species in this complex known from eastern North America: the northern S. diabolicum and the largely southern S. magniae. Methods We sampled plant populations from across a latitudinal gradient in eastern North America and performed whole genome and restriction-site associated DNA sequencing. These sequencing data were then analyzed computationally. Key Results Using sliding-window population genetic analyses we find that differentiation is concentrated within ‘islands’ of the genome spanning up to 400 kb that are characterized by elevated genetic divergence, suppressed recombination, reduced nucleotide diversity and increased rates of non-synonymous substitution. Sequence variants that are significantly associated with genetic structure and bioclimatic variables occur within genes that have functional enrichment for biological processes including abiotic stress response, photoperiodism and hormone-mediated signalling. Demographic modelling demonstrates that these two species diverged no more than 225 000 generations ago with secondary contact occurring where their ranges overlap. Conclusions We suggest that this heterogeneity of genomic differentiation is a result of linked selection and reflects the role of local adaptation to contrasting climatic zones in driving speciation. This research provides insight into the process of speciation in a group of ecologically important plants and strengthens our predictive understanding of how plant populations will respond as Earth’s climate rapidly changes.

58 GEOSCIENCES↗

A variant selection framework for genome graphs

Abstract Motivation Variation graph representations are projected to either replace or supplement conventional single genome references due to their ability to capture population genetic diversity and reduce reference bias. Vast catalogues of genetic variants for many species now exist, and it is natural to ask which among these are crucial to circumvent reference bias during read mapping. Results In this work, we propose a novel mathematical framework for variant selection, by casting it in terms of minimizing variation graph size subject to preserving paths of length α with at most δ differences. This framework leads to a rich set of problems based on the types of variants [e.g. single nucleotide polymorphisms (SNPs), indels or structural variants (SVs)], and whether the goal is to minimize the number of positions at which variants are listed or to minimize the total number of variants listed. We classify the computational complexity of these problems and provide efficient algorithms along with their software implementation when feasible. We empirically evaluate the magnitude of graph reduction achieved in human chromosome variation graphs using multiple α and δ parameter values corresponding to short and long-read resequencing characteristics. When our algorithm is run with parameter settings amenable to long-read mapping (α = 10 kbp, δ = 1000), 99.99% SNPs and 73% SVs can be safely excluded from human chromosome 1 variation graph. The graph size reduction can benefit downstream pan-genome analysis. Availability and implementation https://github.com/AT-CG/VF. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗

MIDAS2: Metagenomic Intra-species Diversity Analysis System

The Metagenomic Intra-Species Diversity Analysis System (MIDAS) is a scalable metagenomic pipeline that identifies single nucleotide variants (SNVs) and gene copy number variants in microbial populations. Here, we present MIDAS2, which addresses the computational challenges presented by increasingly large reference genome databases, while adding functionality for building custom databases and leveraging paired-end reads to improve SNV accuracy. This fast and scalable reengineering of the MIDAS pipeline enables thousands of metagenomic samples to be efficiently genotyped.

59 BASIC BIOLOGICAL SCIENCES↗

Meta2DB: Curated Shotgun Metagenomic Feature Sets and Metadata for Health State Prediction

Meta2DB is a curated metagenomic and metadata database that provides structurally consistent microbiome taxonomy feature count tables for 13 897 samples across 84 studies, 23 disease states, and 34 geographical locations. All samples were uniformly processed using a streamlined metagenomic classification pipeline that employs a unique and comprehensive reference database indexed to contain all sequences across all kingdoms of life that were present in the NCBI Nucleotide (nt) database retrieved on 4 January 2023. This pipeline leverages high-performance computing (HPC) resources at Lawrence Livermore National Laboratory and was used to process 50TB of publicly available raw metagenomic sequence data. Extensive metadata curation was carried out through a combination of manual curation and automated parsing, producing a consistent inter-study metadata table specifically structured to facilitate training of ML models for prediction of human health.

Kok, C [Lawrence Livermore National Laboratory (LL↗

Whole-genome sequencing distinguishes the two most common giant kelp ecomorphs

Abstract Giant kelp, Macrocystis pyrifera, exists as distinct morphological variants—or “ecomorphs”—in different populations, yet the mechanism for this variation is uncertain, and environmental drivers for either adaptive or plastic phenotypes have not been identified. The ecomorphs Macrocystis “pyrifera” and M. “integrifolia” are distributed throughout temperate waters of North and South America with almost no geographic overlap and exhibit an incongruous, non-mirrored, distribution across the equator. This study evaluates the degree of genetic divergence between M. “pyrifera” and M. “integrifolia” across 18 populations in Chile and California using whole-genome sequencing and single-nucleotide polymorphism markers. Our results based on a principal component analysis, admixture clustering by genetic similarity, and phylogenetic inference demonstrate that M. “pyrifera” and M. “integrifolia” are genetically distinguishable. Analyses reveal separation by Northern and Southern Hemispheres and between morphs within hemispheres, suggesting that the convergent “integrifolia” morphology arose separately in each hemisphere. This is the first study to use whole-genome sequencing to understand genetic divergence in giant kelp ecomorphs, identifying 83 potential genes under selection and providing novel insights about Macrocystis evolution that were not evident with previous genetic techniques. Future studies are needed to uncover the environmental forces driving local adaptation and presumed convergent evolution of these morphs.

Environmental Sciences & Ecology↗

Biogeographic patterns in populations of marine Pseudoalteromonas atlantica isolates

Intra-specific genomic diversity is well documented in microbes. The question, however, remains whether natural selection or neutral evolution is the major contributor to this diversity. We undertook this study to estimate genomic diversity in Pseudoalteromonas atlantica populations and whether the diversity, if present, could be attributed to environmental factors or distance effects. Here, we isolated and sequenced twenty-three strains of P. atlantica from three geographically distant deep marine basins and performed comparative genomic analyses to study the genomic diversity of populations among these basins. Average nucleotide identity followed a strictly geographical pattern. In two out of three locations, the strains within the location exhibited >99.5% identity, whereas, among locations, the strains showed <98.11% identity. Phylogenetic and pan-genome analysis also reflected the biogeographical separation of the strains. Strains from the same location shared many accessory genes and clustered closely on the phylogenetic tree. Phenotypic diversity between populations was studied in ten out of twenty-three strains testing carbon and nitrogen source utilization and osmotolerance. A genetic basis for phenotypic diversity could be established in most cases but was apparently not influenced by local environmental conditions. Our study suggests that neutral evolution may have a substantial role in the biodiversity of P. atlantica.

59 BASIC BIOLOGICAL SCIENCES↗