Stacking up RADSeq assembly programs: From complete hit to completely abysmal

Decreasing sequencing costs have driven a rapid expansion of novel genotyping methods. One of these methods is the exploitation of restriction enzyme cut sites to generate genome‐wide but reduced representation sequencing libraries (RRLs), alternatively termed genotyping by sequencing or restriction‐site associated DNA sequencing. Without a reference genome, the resulting short sequence reads must be assembled de novo. There are many possible assembly programs, most not explicitly developed for RRL data, and we know little of their effectiveness. In this issue of Molecular Ecology Resources, LaCava et al. (2020) systematically evaluate six commonly used programs and two commonly varied parameters for complete and accurate assembly of RRLs, using simulated double digests of Homo sapiens and Arabidopsis thaliana genomes with varied mutation rates and types. The authors find substantial variation in performance across assembly programs. The most consistently high‐performing assembler is infrequently used in their literature survey (CD‐HIT; Li and Godzik, 2006), while several others fail to produce complete, accurate assemblies under many conditions. LaCava et al. additionally recommend best practices in parameter choice and evaluation of future assembly programs—advice that molecular ecologists working to assemble sequences of all kinds should take to heart.     

Low invasive estrous synchronization protocol for wild animals: an example with melengestrol acetate in brown brocket deer (Mazama gouazoubira)

Deer are sensitive to stressful stimuli by handling and their reproductive physiology could be altered by these procedures, making it necessary to develop less invasive protocols for ART. Melengestrol acetate (MGA), a synthetic progestin administered orally, appears as an alternative for estrous synchronization protocols (ESP), such as reported in cattle. Firstly, we compared two MGA doses (0.5 and 1.0 mg/day/animal), which would have suppression effect in estrous behavior (EB). Eight females were randomly and equally distributed in Group 1 (G1) and Group 2 (G2), which received 0.5 and 1.0 mg/day/animal respectively for 15 days (D1 to D15). Two cloprostenol (CP) applications were performed on D0 and D11. Estrus detection (ED) was performed every day. All females from G1 displayed estrus during treatment period, whereas all females from G2 displayed estrus after treatment, suggesting a suppressive effect of 1.0 mg in the EB. Once the suppressive MGA dose (1.0 mg) was defined, we used this dose for assessing ESP. The same eight females received 1.0 mg/animal for eight days (D-8 to D-1), followed by 0.25 mg of estradiol benzoate on D-8 and 265 μg of CP on D0. Feces for fecal progesterone metabolites (FPM) measurement were collected from D0 until seven days after the last day of estrus. Seven females displayed estrus between 12 and 72 h after CP application, which was followed by a significant increase in FPM levels (except female MG6), suggesting the formation of corpus luteum. After ED, females were placed with a fertile male to assess the fertility of the protocol. Pregnancy was confirmed by ultrasound 30 days after mating in 3/6 individuals. Although the low effectiveness of MGA protocol, it should be considered as a promising alternative in deer ESP since this protocol has less stressful effect on the animal during reproductive management when compared to other ESP.     

Coupling between tolerance and resistance for two related Eimeria parasite species

Resistance (host capacity to reduce parasite burden) and tolerance (host capacity to reduce impact on its health for a given parasite burden) manifest two different lines of defense. Tolerance can be independent from resistance, traded off against it, or the two can be positively correlated because of redundancy in underlying (immune) processes. We here tested whether this coupling between tolerance and resistance could differ upon infection with closely related parasite species. We tested this in experimental infections with two parasite species of the genus Eimeria. We measured proxies for resistance (the (inverse of) number of parasite transmission stages (oocysts) per gram of feces at the day of maximal shedding) and tolerance (the slope of maximum relative weight loss compared to day of infection on number of oocysts per gram of feces at the day of maximal shedding for each host strain) in four inbred mouse strains and four groups of F1 hybrids belonging to two mouse subspecies, Mus musculus domesticus and Mus musculus musculus. We found a negative correlation between resistance and tolerance against Eimeria falciformis, while the two are uncoupled against Eimeria ferrisi. We conclude that resistance and tolerance against the first parasite species might be traded off, but evolve more independently in different mouse genotypes against the latter. We argue that evolution of the host immune defenses can be studied largely irrespective of parasite isolates if resistance–tolerance coupling is absent or weak (E. ferrisi) but host–parasite coevolution is more likely observable and best studied in a system with negatively correlated tolerance and resistance (E. falciformis).     

Genetic insights into the dark matter of the mammalian gut microbiota through targeted genome reconstruction

Whole metagenomic shotgun (WMS) sequencing has dramatically enhanced our ability to study microbial genomics. The possibility to unveil the genetic makeup of bacteria that cannot be easily isolated has significantly expanded our microbiological horizon. Here, we report an approach aimed at uncovering novel bacterial species by the use of targeted WMS sequencing. Employing in silico data retrieved from metabolic modelling to formulate a chemically defined medium (CDM), we were able to isolate and subsequently sequence the genomes of six putative novel species of bacteria from the gut of non-human primates.     

Recent loss of the Dim2 DNA methyltransferase decreases mutation rate in repeats and changes evolutionary trajectory in a fungal pathogen

DNA methylation is found throughout all domains of life, yet the extent and function of DNA methylation differ among eukaryotes. Strains of the plant pathogenic fungus Zymoseptoria tritici appeared to lack cytosine DNA methylation (5mC) because gene amplification followed by Repeat-Induced Point mutation (RIP) resulted in the inactivation of the dim2 DNA methyltransferase gene. 5mC is, however, present in closely related sister species. We demonstrate that inactivation of dim2 occurred recently as some Z. tritici isolates carry a functional dim2 gene. Moreover, we show that dim2 inactivation occurred by a different path than previously hypothesized. We mapped the genome-wide distribution of 5mC in strains with or without functional dim2 alleles. Presence of functional dim2 correlates with high levels of 5mC in transposable elements (TEs), suggesting a role in genome defense. We identified low levels of 5mC in strains carrying non-functional dim2 alleles, suggesting that 5mC is maintained over time, presumably by an active Dnmt5 DNA methyltransferase. Integration of a functional dim2 allele in strains with mutated dim2 restored normal 5mC levels, demonstrating de novo cytosine methylation activity of Dim2. To assess the importance of 5mC for genome evolution, we performed an evolution experiment, comparing genomes of strains with high levels of 5mC to genomes of strains lacking functional dim2. We found that presence of a functional dim2 allele alters nucleotide composition by promoting C to T transitions (C→T) specifically at CpA (CA) sites during mitosis, likely contributing to TE inactivation. Our results show that 5mC density at TEs is a polymorphic trait in Z. tritici populations that can impact genome evolution.     

SYMBIODINIUM (PYRRHOPHYTA) GENOME SIZES (DNA CONTENT) ARE SMALLEST AMONG DINOFLAGELLATES1

Using flow cytometric analysis of fluorescence, we measured the genome sizes of 18 cultured “free‐living” species and 29 Symbiodinium spp. isolates cultured from stony corals, gorgonians, anemones, jellyfish, and giant clams. Genome size directly correlated with cell size, as documented previously for most eukaryotic cell lines. Among the smallest of dinoflagellates, Symbiodinium spp. (6–15 μm) possessed the lowest DNA content that we measured (1.5–4.8 pg·cell^?1). Bloom‐forming or potentially harmful species in the genera Alexandrium, Karenia, Pfiesteria, and Prorocentrum possessed genomes approximately 2 to 50 times larger in size. A phylogenetic analysis indicated that genome/cell size has apparently increased and decreased repeatedly during the evolution of dinoflagellates. In contrast, genome sizes were relatively consistent across distantly and closely related Symbiodinium spp. This may be the product of intracellular host habitats imposing strong selective pressures that have restricted symbiont size.     

The binding and release of oxygen and hydrogen peroxide are directed by a hydrophobic tunnel in cholesterol oxidase

The usage by enzymes of specific binding pathways for gaseous substrates or products is debated. The crystal structure of the redox enzyme cholesterol oxidase, determined at sub-angstrom resolution, revealed a hydrophobic tunnel that may serve as a binding pathway for oxygen and hydrogen peroxide. This tunnel is formed by a cascade of conformational rearrangements and connects the active site with the exterior surface of the protein. To elucidate the relationship between this tunnel and gas binding and release, three mutant enzymes were constructed to block the tunnel or its putative gate. Mutation of the proposed gating residue Asn485 to Asp or tunnel residue Phe359 or Gly347 to Trp or Asn reduces the catalytic efficiency of oxidation. The K mO 2 increases from 300 +/- 35 microM for the wild-type enzyme to 617 +/- 15 microM for the F359W mutant. The k cat for the F359W mutant-catalyzed reaction decreases 13-fold relative to that of the wild-type-catalyzed reaction. The N485D and G347N mutants could not be saturated with oxygen. Transfer of hydride from the sterol to the flavin prosthetic group is no longer rate-limiting for these tunnel mutants. The steady-state kinetics of both wild-type and tunnel mutant enzymes are consistent with formation of a ternary complex of steroid and oxygen during catalysis. Furthermore, kinetic cooperativity with respect to molecular oxygen is observed with the tunnel mutants, but not with the wild-type enzyme. A rate-limiting conformational change for binding and release of oxygen and hydrogen peroxide, respectively, is consistent with the cooperative kinetics. In the atomic-resolution structure of F359W, the indole ring of the tryptophan completely fills the tunnel and is observed in only a single conformation. The size of the indole is proposed to limit conformational rearrangement of residue 359 that leads to tunnel opening in the wild-type enzyme. Overall, these results substantiate the functional importance of the tunnel for substrate binding and product release.     

Meta‐analysis of global livestock urine‐derived nitrous oxide emissions from agricultural soils

Nitrous oxide (N₂O) is an air pollutant of major environmental concern, with agriculture representing 60% of anthropogenic global N₂O emissions. Much of the N₂O emissions from livestock production systems result from transformation of N deposited to soil within animal excreta. There exists a substantial body of literature on urine patch N₂O dynamics, we aimed to identify key controlling factors influencing N₂O emissions and to aid understanding of knowledge gaps to improve GHG reporting and prioritize future research. We conducted an extensive literature review and random effect meta‐analysis (using REML) of results to identify key relationships between multiple potential independent factors and global N₂O emissions factors (EFs) from urine patches. Mean air temperature, soil pH and ruminant animal species (sheep or cow) were significant factors influencing the EFs reviewed. However, several factors that are known to influence N₂O emissions, such as animal diet and urine composition, could not be considered due to the lack of reported data. The review highlighted a widespread tendency for inadequate metadata and uncertainty reporting in the published studies, as well as the limited geographical extent of investigations, which are more often conducted in temperate regions thus far. Therefore, here we give recommendations for factors that are likely to affect the EFs and should be included in all future studies, these include the following: soil pH and texture; experimental set‐up; direct measurement of soil moisture and temperature during the study period; amount and composition of urine applied; animal type and diet; N₂O emissions with a measure of uncertainty; data from a control with zero‐N application and meteorological data.