The distribution of coastal fish eDNA sequences in the Anthropocene

Aim

Coastal fishes have a fundamental role in marine ecosystem functioning and contributions to people, but face increasing threats due to climate change, habitat degradation and overexploitation. The extent to which human pressures are impacting coastal fish biodiversity in comparison with geographic and environmental factors at large spatial scale is still under scrutiny. Here, we took advantage of environmental DNA (eDNA) metabarcoding to investigate the relationship between fish biodiversity, including taxonomic and genetic components, and environmental but also socio-economic factors.

Location

Tropical, temperate and polar coastal areas.

Time period

Present day.

Major taxa studied

Marine fishes.

Methods

We analysed fish eDNA in 263 stations (samples) in 68 sites distributed across polar, temperate and tropical regions. We modelled the effect of environmental, geographic and socio-economic factors on α- and β-diversity. We then computed the partial effect of each factor on several fish biodiversity components using taxonomic molecular units (MOTU) and genetic sequences. We also investigated the relationship between fish genetic α- and β-diversity measured from our barcodes, and phylogenetic but also functional diversity.

Results

We show that fish eDNA MOTU and sequence α- and β-diversity have the strongest correlation with environmental factors on coastal ecosystems worldwide. However, our models also reveal a negative correlation between biodiversity and human dependence on marine ecosystems. In areas with high dependence, diversity of all fish, cryptobenthic fish and large fish MOTUs declined steeply. Finally, we show that a sequence diversity index, accounting for genetic distance between pairs of MOTUs, within and between communities, is a reliable proxy of phylogenetic and functional diversity.

Main conclusions

Together, our results demonstrate that short eDNA sequences can be used to assess climate and direct human impacts on marine biodiversity at large scale in the Anthropocene and can further be extended to investigate biodiversity in its phylogenetic and functional dimensions.     

Trade-offs and synergies between ecosystem productivity and stability in temperate grasslands

Aim

It is crucial to monitor how the productivity of grasslands varies with its temporal stability for management of these ecosystems. However, identifying the direction of the productivity–stability relationship remains challenging because ecological stability has multiple components that can display neutral, positive or negative covariations. Furthermore, evidence suggests that the direction of the productivity–stability relationship depends on the biotic interactions and abiotic conditions that underlie ecosystem productivity and stability. We decipher the relationships between grassland productivity and two components of its stability in four habitat types with contrasting environments and flora.

Location

France.

Time period

2000–2020.

Major taxa

Grassland plant species.

Methods

We used c. 20,000 vegetation plots spread across French permanent grasslands and remotely sensed vegetation indices to quantify grassland productivity and temporal stability. We decomposed stability into constancy (i.e., temporal invariability) and resistance (i.e., maximum deviation from average) and deciphered the direct and indirect effects of abiotic (namely growing season length and nitrogen input) and biotic (namely plant taxonomic diversity, trait diversity and community-weighted mean traits) factors on productivity–stability relationships using structural equation models.

Results

We found a positive relationship between productivity and constancy and a negative relationship between productivity and resistance in all habitats. Abiotic factors had stronger effects on productivity and stability compared with biotic factors. A longer growing season enhanced grassland productivity and constancy. Nitrogen input had positive and negative effects on grassland productivity and resistance, respectively. Trait values affected the constancy and resistance of grassland more than taxonomic and trait diversity, with effects varying from one habitat to another. Productivity was not related to any biotic factor.

Main conclusions

Our findings reveal how vital it is to consider both the multiple components of stability and the interaction between environment and biodiversity to gain an understanding of the relationships between productivity and stability in real-world ecosystems, which is a crucial step for sustainable grassland management.     

Canopy tree mortality depends on the proportion of crown exposed to sunlight,but this effect varies with species' wood density

Understanding what drives changes in tree mortality as well as the covariates influencing trees' response is a research priority to predict forest responses to global change. Here, we combined drone photogrammetry and ground-based data to assess the influence of crown exposure to light (relative to total crown area), growth deviations (relative to conspecifics), tree size, and species' wood density (as a surrogate for light-demanding and shade-tolerant life-history strategies) on the mortality of 984 canopy trees in an Amazon terra firme forest. Trees with lower wood density were less prone to die when their proportion of crown was more exposed to sunlight, but this relationship with relative crown exposure weakened and slightly reversed as wood density increased. Trees growing less than their species average had higher mortality, especially when the species' wood density decreased. The role of wood density in determining the survival of canopy trees under varying light conditions indicates differential responses of light-demanding versus shade-tolerant species. Our results highlight the importance of accounting for life-history strategies, via plant functional types, in vegetation dynamic models aiming to predict forest demography under a rapidly changing climate. Abstract in Spanish is available with online material.     

Construction of Chemical Libraries of Volatile Compounds by Combinatorial Synthesis of Homologous Mixtures: Alk-4-en-1-ols,Alk-4-enals and Methyl Alk-4-enoates

A compound library of sixty six linear compounds, eleven representatives of six molecular families: (E)- and (Z)-isomers of alk-4-en-1-ols, alk-4-enals, and methyl alk-4-enoates, was prepared by combinatorial syntheses to allow the creation of a mass spectral database directly usable for their identification in GC/MS analyses. We demonstrate here that compound libraries can be prepared by combinatorial syntheses using long linear synthetic sequences, i. e., eight step in the case of 4-enals. The resulting mixtures of homologues are still perfectly exploitable to deliver the requested information such as clean mass spectra and good gas chromatographic retention indices.     

Statistical shape modelling of the thoracic spine for the development of pedicle screw insertion guides

Spinal fixation and fusion are surgical procedures undertaken to restore stability in the spine and restrict painful or degenerative motion. Malpositioning of pedicle screws during these procedures can result in major neurological and vascular damage. Patient-specific surgical guides offer clear benefits, reducing malposition rates by up to 25%. However, they suffer from long lead times and the manufacturing process is dependent on third-party specialists. The development of a standard set of surgical guides may eliminate the issues with the manufacturing process. To evaluate the feasibility of this option, a statistical shape model (SSM) was created and used to analyse the morphological variations of the T4–T6 vertebrae in a population of 90 specimens from the Visible Korean Human dataset (50 females and 40 males). The first three principal components, representing 39.7% of the variance within the population, were analysed. The model showed high variability in the transverse process (~ 4 mm) and spinous process (~ 4 mm) and relatively low variation (< 1 mm) in the vertebral lamina. For a Korean population, a standardised set of surgical guides would likely need to align with the lamina where the variance in the population is lower. It is recommended that this standard set of surgical guides should accommodate pedicle screw diameters of 3.5–6 mm and transverse pedicle screw angles of 3.5°–12.4°.

     

A robust approach to estimate relative phytoplankton cell abundances from metagenomes

Phytoplankton account for >45% of global primary production, and have an enormous impact on aquatic food webs and on the entire Earth System. Their members are found among prokaryotes (cyanobacteria) and multiple eukaryotic lineages containing chloroplasts. Genetic surveys of phytoplankton communities generally consist of PCR amplification of bacterial (16S), nuclear (18S) and/or chloroplastic (16S) rRNA marker genes from DNA extracted from environmental samples. However, our appreciation of phytoplankton abundance or biomass is limited by PCR-amplification biases, rRNA gene copy number variations across taxa, and the fact that rRNA genes do not provide insights into metabolic traits such as photosynthesis. Here, we targeted the photosynthetic gene psbO from metagenomes to circumvent these limitations: the method is PCR-free, and the gene is universally and exclusively present in photosynthetic prokaryotes and eukaryotes, mainly in one copy per genome. We applied and validated this new strategy with the size-fractionated marine samples collected by Tara Oceans, and showed improved correlations with flow cytometry and microscopy than when based on rRNA genes. Furthermore, we revealed unexpected features of the ecology of these ecosystems, such as the high abundance of picocyanobacterial aggregates and symbionts in the ocean, and the decrease in relative abundance of phototrophs towards the larger size classes of marine dinoflagellates. To facilitate the incorporation of psbO in molecular-based surveys, we compiled a curated database of >18,000 unique sequences. Overall, psbO appears to be a promising new gene marker for molecular-based evaluations of entire phytoplankton communities.     

First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.     

INDIVIDUALITY IN THE VOICE OF FUR SEAL FEMALES: AN ANALYSIS STUDY OF THE PUP ATTRACTION CALL IN ARCTOCEPHALUS TROPICALIS

Like most otariids species, the Subantarctic fur seal breeds on land in large, dense colonies. Pups are confronted by the long and repetitive absences of their mother throughout lactation. At each mother's return, pups have to find her among several hundreds of congeners. This recognition process mainly relies on acoustic signals. We performed an acoustic analysis on 125 calls from 20 females recorded during the 1999–2000 breeding season on Amsterdam Island (Indian Ocean). Ten variables were measured in both temporal and frequency domains. To find the acoustic parameters supporting individual signature, we assessed the differences between individuals using Kruskall-Wallis univariate analysis of variance. For each variable, we also calculated the potential of individuality coding (PIC) as the ratio between the between-individual coefficient of variation and the mean value of the within-individual coefficients of variation. We found that the frequency spectrum, the characteristics of the frequency modulation of the initial and middle part of the call and the call duration exhibit an important individual stereotypy (PIC values ranging between 1.5 and 3), whereas features relative to amplitude and the frequency modulation of the final part of the call are weakly individualized (PIC values between 1 and 1.2).     

Characterization of an Escherichia coli mutant, feeA, displaying resistance to the calmodulin inhibitor 48/80 and reduced expression of the rare tRNA3Leu

We previously described a mutation feeB1 conferring a temperature-sensitive filamentation phenotype and resistance to the calmodulin inhibitor 48/80 in Escherichia coli, which constitutes a single base change in the acceptor stem of the rare tRNA₃^Leu recognizing CUA codons. We now describe a second mutant, feeA1, unlinked to feeB, but displaying a similar phenotype, 48/80 resistance and a reduced growth rate at the permissive temperature, 30°C, and temperature-sensitive, forming short filaments at 42°C. In the feeA mutant, tRNA₃^Leu expression (but not that of tRNA₁^Leu) was reduced approximately fivefold relative to the wild type. We previously showed that the synthesis of β-galactosidase, which unusually requires the translation of 6-CUA codons, was substantially reduced, particularly at 42°C, in feeB mutants. The feeA mutant also shows drastically reduced synthesis of β-galactosidase at the non-permissive temperature and reduced levels even at the permissive temperature. We also show that increased copy numbers of the abundant tRNA₁^Leu, which can also read CUA codons at low efficiency, suppressed the effects of feeA1 under some conditions, providing further evidence that the mutant was deficient in CUA translation. This, and the previous study, demonstrates that mutations which either reduce the activity of tRNA₃^Leu or the cellular amount of tRNA₃^Leu confer resistance to the drug 48/80, with concomitant inhibition of cell division at 42°C.