Neutral vs. non-neutral genetic footprints of Plasmodium falciparum multiclonal infections.

At a time when effective tools for monitoring malaria control and eradication efforts are crucial, the increasing availability of molecular data motivates their application to epidemiology. The multiplicity of infection (MOI), defined as the number of genetically distinct parasite strains co-infecting a host, is one key epidemiological parameter for evaluating malaria interventions. Estimating MOI remains a challenge for high-transmission settings where individuals typically carry multiple co-occurring infections. Several quantitative approaches have been developed to estimate MOI, including two cost-effective ones relying on molecular data: i) THE REAL McCOIL method is based on putatively neutral single nucleotide polymorphism loci, and ii) the varcoding method is a fingerprinting approach that relies on the diversity and limited repertoire overlap of the var multigene family encoding the major Plasmodium falciparum blood-stage antigen PfEMP1 and is therefore under selection. In this study, we assess the robustness of the MOI estimates generated with these two approaches by simulating P. falciparum malaria dynamics under three transmission conditions using an extension of a previously developed stochastic agent-based model. We demonstrate that these approaches are complementary and best considered across distinct transmission intensities. While varcoding can underestimate MOI, it allows robust estimation, especially under high transmission where repertoire overlap is extremely limited from frequency-dependent selection. In contrast, THE REAL McCOIL often considerably overestimates MOI, but still provides reasonable estimates for low and moderate transmission. Regardless of transmission intensity, results for THE REAL McCOIL indicate that an inaccurate tail at high MOI values is generated, and that at high transmission, an apparently reasonable estimated MOI distribution can arise from some degree of compensation between overestimation and underestimation. As many countries pursue malaria elimination targets, defining the most suitable approach to estimate MOI based on sample size and local transmission intensity is highly recommended for monitoring the impact of intervention programs.

Ocean warming threatens key trophic interactions supporting a commercial fishery in a climate change hotspot.

Worldwide, rising ocean temperatures are causing declines and range shifts in marine species. The direct effects of climate change on the biology of marine organisms are often well documented; yet, knowledge on the indirect effects, particularly through trophic interactions, is largely lacking. We provide evidence of ocean warming decoupling critical trophic interactions supporting a commercially important mollusc in a climate change hotspot. Dietary assessments of the Australian blacklip abalone (Haliotis rubra) indicate primary dependency on a widespread macroalgal species (Phyllospora comosa) which we show to be in state of decline due to ocean warming, resulting in abalone biomass reductions. Niche models suggest further declines in P. comosa over the coming decades and ongoing risks to H. rubra. This study highlights the importance of studies from climate change hotspots and understanding the interplay between climate and trophic interactions when determining the likely response of marine species to environmental changes.

First Genome of Labyrinthula sp., an Opportunistic Seagrass Pathogen, Reveals Novel Insight into Marine Protist Phylogeny, Ecology and CAZyme Cell-Wall Degradation.

Labyrinthula spp. are saprobic, marine protists that also act as opportunistic pathogens and are the causative agents of seagrass wasting disease (SWD). Despite the threat of local- and large-scale SWD outbreaks, there are currently gaps in our understanding of the drivers of SWD, particularly surrounding Labyrinthula spp. virulence and ecology. Given these uncertainties, we investigated the Labyrinthula genus from a novel genomic perspective by presenting the first draft genome and predicted proteome of a pathogenic isolate Labyrinthula SR_Ha_C, generated from a hybrid assembly of Nanopore and Illumina sequences. Phylogenetic and cross-phyla comparisons revealed insights into the evolutionary history of Stramenopiles. Genome annotation showed evidence of glideosome-type machinery and an apicoplast protein typically found in protist pathogens and parasites. Proteins involved in Labyrinthula SR_Ha_C's actin-myosin mode of transport, as well as carbohydrate degradation were also prevalent. Further, CAZyme functional predictions revealed a repertoire of enzymes involved in breakdown of cell-wall and carbohydrate storage compounds common to seagrasses. The relatively low number of CAZymes annotated from the genome of Labyrinthula SR_Ha_C compared to other Labyrinthulea species may reflect the conservative annotation parameters, a specialized substrate affinity and the scarcity of characterized protist enzymes. Inherently, there is high probability for finding both unique and novel enzymes from Labyrinthula spp. This study provides resources for further exploration of Labyrinthula spp. ecology and evolution, and will hopefully be the catalyst for new hypothesis-driven SWD research revealing more details of molecular interactions between the Labyrinthula genus and its host substrate.

Cytosine methylation of rice mitochondrial DNA from grain and leaf tissues.

MAIN CONCLUSION: The rice leaf mitochondrial DNA is  more methylated compared to the rice grain mitochondrial DNA. The old rice leaf mitochondrial DNA has also a higher methylation level than the young rice leaf mitochondrial DNA. The presence of DNA methylation in rice organelles has not been well characterized. We have previously shown that cytosine methylation of chloroplast DNA is different between leaf and grain, and varies between young and old leaves in rice. However, the variation in cytosine methylation of mitochondrial DNA is still poorly characterized. In this study, we have investigated cytosine methylation of mitochondrial DNA in the rice grain and leaf. Based on CpG, CHG, and CHH methylation analyses, the leaf mitochondrial DNA was found to be  more methylated compared to the grain mitochondrial DNA. The methylation of the leaf mitochondrial DNA was also higher in old compared to young leaves. Differences in methylation were observed at different cytosine positions of the mitochondrial DNA between grain and leaf, although there were also positions with a similar level of high methylation in all the tissues examined. The differentially methylated cytosine positions in rice mitochondrial DNA were observed mostly in the intergenic region and in some mitochondrial-specific genes involved in ATP production, transcription, and translation. The functional importance of cytosine methylation in the life cycle of rice mitochondria is still to be determined.

Comparative sequence and methylation analysis of chloroplast and amyloplast genomes from rice.

KEY MESSAGE: Grain amyloplast and leaf chloroplast DNA sequences are identical in rice plants but are differentially methylated. The leaf chloroplast DNA becomes more methylated as the rice plant ages. Rice is an important crop worldwide. Chloroplasts and amyloplasts are critical organelles but the amyloplast genome is poorly studied. We have characterised the sequence and methylation of grain amyloplast DNA and leaf chloroplast DNA in rice. We have also analysed the changes in methylation patterns in the chloroplast DNA as the rice plant ages. Total genomic DNA from grain, old leaf and young leaf tissues were extracted from the Oryza sativa ssp. indica cv. MR219 and sequenced using Illumina Miseq. Sequence variant analysis revealed that the amyloplast and chloroplast DNA of MR219 were identical to each other. However, comparison of CpG and CHG methylation between the identical amyloplast and chloroplast DNA sequences indicated that the chloroplast DNA from rice leaves collected at early ripening stage was more methylated than the amyloplast DNA from the grains of the same plant. The chloroplast DNA became more methylated as the plant ages so that chloroplast DNA from young leaves was less methylated overall than amyloplast DNA. These differential methylation patterns were primarily observed in organelle-encoded genes related to photosynthesis followed by those involved in transcription and translation.

Local and regional scale habitat heterogeneity contribute to genetic adaptation in a commercially important marine mollusc (Haliotis rubra) from southeastern Australia.

Characterising adaptive genetic divergence among conspecific populations is often achieved by studying genetic variation across defined environmental gradients. In marine systems this is challenging due to a paucity of information on habitat heterogeneity at local and regional scales and a dependency on sampling regimes that are typically limited to broad longitudinal and latitudinal environmental gradients. As a result, the spatial scales at which selection processes operate and the environmental factors that contribute to genetic adaptation in marine systems are likely to be unclear. In this study we explore patterns of adaptive genetic structuring in a commercially- harvested abalone species (Haliotis rubra) from southeastern Australia, using a panel of genome-wide SNP markers (5,239 SNPs), and a sampling regime informed by marine LiDAR bathymetric imagery and 20-year hindcasted oceanographic models. Despite a lack of overall genetic structure across the sampling distribution, significant genotype associations with heterogeneous habitat features were observed at local and regional spatial scales, including associations with wave energy, ocean current, sea surface temperature, and geology. These findings provide insights into the potential resilience of the species to changing marine climates and the role of migration and selection on recruitment processes, with implications for conservation and fisheries management. This study points to the spatial scales at which selection processes operate in marine systems and highlights the benefits of geospatially-informed sampling regimes for overcoming limitations associated with marine population genomic research.

More evolution underground: Accelerated mitochondrial substitution rate in Australian burrowing freshwater crayfishes (Decapoda: Parastacidae).

To further understand the evolutionary history and mitogenomic features of Australia's highly distinctive freshwater crayfish fauna, we utilized a recently described rapid mitogenome sequencing pipeline to generate 24 new crayfish mitogenomes including a diversity of burrowing crayfish species and the first for Astacopsis gouldi, the world's largest freshwater invertebrate. Whole mitogenome-based phylogeny estimates using both Bayesian and Maximum Likelihood methods substantially strengthen existing hypotheses for systematic relationships among Australian freshwater crayfish with evidence of pervasive diversifying selection and accelerated mitochondrial substitution rate among the members of the clade representing strongly burrowing crayfish that may reflect selection pressures for increased energy requirement for adaptation to terrestrial environment and a burrowing lifestyle. Further, gene rearrangements are prevalent in the burrowing crayfish mitogenomes involving both tRNA and protein coding genes. In addition, duplicated control regions were observed in two closely related Engaeus species, together with evidence for concerted evolution. This study significantly adds to the understanding of Australian freshwater crayfish evolutionary relationships and suggests a link between mitogenome evolution and adaptation to terrestrial environments and a burrowing lifestyle in freshwater crayfish.

ORDER within the chaos: Insights into phylogenetic relationships within the Anomura (Crustacea: Decapoda) from mitochondrial sequences and gene order rearrangements.

The infraorder Anomura consists of a morphologically and ecologically heterogeneous group of decapod crustaceans, and has attracted interest from taxonomists for decades attempting to find some order out of the seemingly chaotic diversity within the group. Species-level diversity within the Anomura runs the gamut from the "hairy" spindly-legged yeti crab found in deep-sea hydrothermal vent environments to the largest known terrestrial invertebrate, the robust coconut or robber crab. Owing to a well-developed capacity for parallel evolution, as evidenced by the occurrence of multiple independent carcinization events, Anomura has long tested the patience and skill of both taxonomists attempting to find order, and phylogeneticists trying to establish stable hypotheses of evolutionary inter-relationships. In this study, we performed genome skimming to recover the mitogenome sequences of 12 anomuran species including the world's largest extant invertebrate, the robber crab (Birgus latro), thereby over doubling these resources for this group, together with 8 new brachyuran mitogenomes. Maximum-likelihood (ML) and Bayesian-inferred (BI) phylogenetic reconstructions based on amino acid sequences from mitogenome protein-coding genes provided strong support for the monophyly of the Anomura and Brachyura and their sister relationship, consistent with previous studies. The majority of relationships within families were supported and were largely consistent with current taxonomic classifications, whereas many relationships at higher taxonomic levels were unresolved. Nevertheless, we have strong support for a polyphyletic Paguroidea and recovered a well-supported clade of a subset of paguroids (Diogenidae + Coenobitidae) basal to all other anomurans, though this requires further testing with greater taxonomic sampling. We also introduce a new feature to the MitoPhAST bioinformatics pipeline (https://github.com/mht85/MitoPhAST) that enables the extraction of mitochondrial gene order (MGO) information directly from GenBank files and clusters groups based on common MGOs. Using this tool, we compared MGOs across the Anomura and Brachyura, identifying Anomura as a taxonomic "hot spot" with high variability in MGOs among congeneric species from multiple families while noting the broad association of highly-rearranged MGOs with several anomuran lineages inhabiting extreme niches. We also demonstrate the value of MGOs as a source of novel synapomorphies for independently reinforcing tree-based relationships and for shedding light on relationships among challenging groups such as the Aegloidea and Lomisoidea that were unresolved in phylogenetic reconstructions. Overall, this study contributes a substantial amount of new genetic material for Anomura and attempts to further resolve anomuran evolutionary relationships where possible based on a combination of sequence and MGO information. The new feature in MitoPhAST adds to the relatively limited number of bioinformatics tools available for MGO analyses, which can be utilized widely across animal groups.

Characterisation of 12 microsatellite loci in the Vietnamese commercial clam Lutraria rhynchaena Jonas 1844 (Heterodonta: Bivalvia: Mactridae) through next-generation sequencing.

The marine clam Lutraria rhynchaena is gaining popularity as an aquaculture species in Asia. Lutraria populations are present in the wild throughout Vietnam and several stocks have been established and translocated for breeding and aquaculture grow-out purposes. In this study, we demonstrate the feasibility of utilising Illumina next-generation sequencing technology to streamline the identification and genotyping of microsatellite loci from this clam species. Based on an initial partial genome scan, 48 microsatellite markers with similar melting temperatures were identified and characterised. The 12 most suitable polymorphic loci were then genotyped using 51 individuals from a population in Quang Ninh Province, North Vietnam. Genetic variation was low (mean number of alleles per locus = 2.6; mean expected heterozygosity = 0.41). Two loci showed significant deviation from Hardy-Weinberg equilibrium (HWE) and the presence of null alleles, but there was no evidence of linkage disequilibrium among loci. Three additional populations were screened (n = 7-36) to test the geographic utility of the 12 loci, which revealed 100 % successful genotyping in two populations from central Vietnam (Nha Trang). However, a second population from north Vietnam (Co To) could not be successfully genotyped and morphological evidence and mitochondrial variation suggests that this population represents a cryptic species of Lutraria. Comparisons of the Qang Ninh and Nha Trang populations, excluding the 2 loci out of HWE, revealed statistically significant allelic variation at 4 loci. We reported the first microsatellite loci set for the marine clam Lutraria rhynchaena and demonstrated its potential in differentiating clam populations. Additionally, a cryptic species population of Lutraria rhynchaena was identified during initial loci development, underscoring the overlooked diversity of marine clam species in Vietnam and the need to genetically characterise population representatives prior to microsatellite development. The rapid identification and validation of microsatellite loci using next-generation sequencing technology warrant its integration into future microsatellite loci development for key aquaculture species in Vietnam and more generally, aquaculture countries in the South East Asia region.

MitoPhAST, a new automated mitogenomic phylogeny tool in the post-genomic era with a case study of 89 decapod mitogenomes including eight new freshwater crayfish mitogenomes.

The increased rate at which complete mitogenomes are being sequenced and their increasing use for phylogenetic studies have resulted in a bioinformatic bottleneck in preparing and utilising such data for phylogenetic analysis. Hence, we present MitoPhAST, an automated tool that (1) identifies annotated protein-coding gene features and generates a standardised, concatenated and partitioned amino acid alignment directly from complete/partial GenBank/EMBL-format mitogenome flat files, (2) generates a maximum likelihood phylogenetic tree using optimised protein models and (3) reports various mitochondrial genes and sequence information in a table format. To demonstrate the capacity of MitoPhAST in handling a large dataset, we used 81 publicly available decapod mitogenomes, together with eight new complete mitogenomes of Australian freshwater crayfishes, including the first for the genus Gramastacus, to undertake an updated test of the monophyly of the major groups of the order Decapoda and their phylogenetic relationships. The recovered phylogenetic trees using both Bayesian and ML methods support the results of studies using fragments of mtDNA and nuclear markers and other smaller-scale studies using whole mitogenomes. In comparison to the fragment-based phylogenies, nodal support values are generally higher despite reduced taxon sampling suggesting there is value in utilising more fully mitogenomic data. Additionally, the simple table output from MitoPhAST provides an efficient summary and statistical overview of the mitogenomes under study at the gene level, allowing the identification of missing or duplicated genes and gene rearrangements. The finding of new mtDNA gene rearrangements in several genera of Australian freshwater crayfishes indicates that this group has undergone an unusually high rate of evolutionary change for this organelle compared to other major families of decapod crustaceans. As a result, freshwater crayfishes are likely to be a useful model for studies designed to understand the evolution of mtDNA rearrangements. We anticipate that our bioinformatics pipeline will substantially help mitogenome-based studies increase the speed, accuracy and efficiency of phylogenetic studies utilising mitogenome information. MitoPhAST is available for download at https://github.com/mht85/MitoPhAST.

Performance of SNP barcodes to determine genetic diversity and population structure of Plasmodium falciparum in Africa.

Panels of informative biallelic single nucleotide polymorphisms (SNPs) have been proposed to be an economical method to fast-track the population genetic analysis of Plasmodium falciparum in malaria-endemic areas. Whilst used successfully in low-transmission areas where infections are monoclonal and highly related, we present the first study to evaluate the performance of these 24- and 96-SNP molecular barcodes in African countries, characterised by moderate-to-high transmission, where multiclonal infections are prevalent. For SNP barcodes it is generally recommended that the SNPs chosen i) are biallelic, ii) have a minor allele frequency greater than 0.10, and iii) are independently segregating, to minimise bias in the analysis of genetic diversity and population structure. Further, to be standardised and used in many population genetic studies, these barcodes should maintain characteristics i) to iii) across various iv) geographies and v) time points. Using haplotypes generated from the MalariaGEN P. falciparum Community Project version six database, we investigated the ability of these two barcodes to fulfil these criteria in moderate-to-high transmission African populations in 25 sites across 10 countries. Predominantly clinical infections were analysed, with 52.3% found to be multiclonal, generating high proportions of mixed-allele calls (MACs) per isolate thereby impeding haplotype construction. Of the 24- and 96-SNPs, loci were removed if they were not biallelic and had low minor allele frequencies in all study populations, resulting in 20- and 75-SNP barcodes respectively for downstream population genetics analysis. Both SNP barcodes had low expected heterozygosity estimates in these African settings and consequently biased analyses of similarity. Both minor and major allele frequencies were temporally unstable. These SNP barcodes were also shown to identify weak genetic differentiation across large geographic distances based on Mantel Test and DAPC. These results demonstrate that these SNP barcodes are vulnerable to ascertainment bias and as such cannot be used as a standardised approach for malaria surveillance in moderate-to-high transmission areas in Africa, where the greatest genomic diversity of P. falciparum exists at local, regional and country levels.

A paradoxical population structure of var DBLα types in Africa.

The var multigene family encodes the P. falciparum erythrocyte membrane protein 1 (PfEMP1), which is important in host-parasite interaction as a virulence factor and major surface antigen of the blood stages of the parasite, responsible for maintaining chronic infection. Whilst important in the biology of P. falciparum, these genes (50 to 60 genes per parasite genome) are routinely excluded from whole genome analyses due to their hyper-diversity, achieved primarily through recombination. The PfEMP1 head structure almost always consists of a DBLα-CIDR tandem. Categorised into different groups (upsA, upsB, upsC), different head structures have been associated with different ligand-binding affinities and disease severities. We study how conserved individual DBLα types are at the country, regional, and local scales in Sub-Saharan Africa. Using publicly-available sequence datasets and a novel ups classification algorithm, cUps, we performed an in silico exploration of DBLα conservation through time and space in Africa. In all three ups groups, the population structure of DBLα types in Africa consists of variants occurring at rare, low, moderate, and high frequencies. Non-rare variants were found to be temporally stable in a local area in endemic Ghana. When inspected across different geographical scales, we report different levels of conservation; while some DBLα types were consistently found in high frequencies in multiple African countries, others were conserved only locally, signifying local preservation of specific types. Underlying this population pattern is the composition of DBLα types within each isolate DBLα repertoire, revealed to also consist of a mix of types found at rare, low, moderate, and high frequencies in the population. We further discuss the adaptive forces and balancing selection, including host genetic factors, potentially shaping the evolution and diversity of DBLα types in Africa.

Measuring changes in Plasmodium falciparum census population size in response to sequential malaria control interventions.

Here we introduce a new endpoint "census population size" to evaluate the epidemiology and control of Plasmodium falciparum infections, where the parasite, rather than the infected human host, is the unit of measurement. To calculate census population size, we rely on a definition of parasite variation known as multiplicity of infection (MOIvar), based on the hyper-diversity of the var multigene family. We present a Bayesian approach to estimate MOIvar from sequencing and counting the number of unique DBLα tags (or DBLα types) of var genes, and derive from it census population size by summation of MOIvar in the human population. We track changes in this parasite population size and structure through sequential malaria interventions by indoor residual spraying (IRS) and seasonal malaria chemoprevention (SMC) from 2012 to 2017 in an area of high-seasonal malaria transmission in northern Ghana. Following IRS, which reduced transmission intensity by > 90% and decreased parasite prevalence by ~40-50%, significant reductions in var diversity, MOIvar, and population size were observed in ~2,000 humans across all ages. These changes, consistent with the loss of diverse parasite genomes, were short lived and 32-months after IRS was discontinued and SMC was introduced, var diversity and population size rebounded in all age groups except for the younger children (1-5 years) targeted by SMC. Despite major perturbations from IRS and SMC interventions, the parasite population remained very large and retained the var population genetic characteristics of a high-transmission system (high var diversity; low var repertoire similarity) demonstrating the resilience of P. falciparum to short-term interventions in high-burden countries of sub-Saharan Africa.

Mitogenome of the extinct Desert 'rat-kangaroo' times the adaptation to aridity in macropodoids.

The evolution of Australia's distinctive marsupial fauna has long been linked to the onset of continent-wide aridity. However, how this profound climate change event affected the diversification of extant lineages is still hotly debated. Here, we assemble a DNA sequence dataset of Macropodoidea-the clade comprising kangaroos and their relatives-that incorporates a complete mitogenome for the Desert 'rat-kangaroo', Caloprymnus campestris. This enigmatic species went extinct nearly 90 years ago and is known from a handful of museum specimens. Caloprymnus is significant because it was the only macropodoid restricted to extreme desert environments, and therefore calibrates the group's specialisation for increasingly arid conditions. Our robustly supported phylogenies nest Caloprymnus amongst the bettongs Aepyprymnus and Bettongia. Dated ancestral range estimations further reveal that the Caloprymnus-Bettongia lineage originated in nascent xeric settings during the middle to late Miocene, ~ 12 million years ago (Ma), but subsequently radiated into fragmenting mesic habitats after the Pliocene to mid-Pleistocene. This timeframe parallels the ancestral divergences of kangaroos in woodlands and forests, but predates their adaptive dispersal into proliferating dry shrublands and grasslands from the late Miocene to mid-Pleistocene, after ~ 7 Ma. We thus demonstrate that protracted changes in both climate and vegetation likely staged the emergence of modern arid zone macropodoids.

Unravelling var complexity: Relationship between DBLα types and var genes in Plasmodium falciparum.

The enormous diversity and complexity of var genes that diversify rapidly by recombination has led to the exclusion of assembly of these genes from major genome initiatives (e.g., Pf6). A scalable solution in epidemiological surveillance of var genes is to use a small 'tag' region encoding the immunogenic DBLα domain as a marker to estimate var diversity. As var genes diversify by recombination, it is not clear the extent to which the same tag can appear in multiple var genes. This relationship between marker and gene has not been investigated in natural populations. Analyses of in vitro recombination within and between var genes have suggested that this relationship would not be exclusive. Using a dataset of publicly-available assembled var sequences, we test this hypothesis by studying DBLα-var relationships for four study sites in four countries: Pursat (Cambodia) and Mae Sot (Thailand), representing low malaria transmission, and Navrongo (Ghana) and Chikwawa (Malawi), representing high malaria transmission. In all study sites, DBLα-var relationships were shown to be predominantly 1-to-1, followed by a second largest proportion of 1-to-2 DBLα-var relationships. This finding indicates that DBLα tags can be used to estimate not just DBLα diversity but var gene diversity when applied in a local endemic area. Epidemiological applications of this result are discussed.

Transcriptomic and Histological Analysis of the Greentail Prawn (Metapenaeus bennettae) Following Light Crude Oil Exposure.

Oil spills pose a significant threat to marine biodiversity. Crude oil can partition into sediments where it may be persistent, placing benthic species such as decapods at particular risk of exposure. Transcriptomic and histological tools are often used to investigate the effects of hydrocarbon exposure on marine organisms following oil spill events, allowing for the identification of metabolic pathways impacted by oil exposure. However, there is limited information available for decapod crustaceans, many of which carry significant economic value. In the present study, we assess the sublethal impacts of crude oil exposure in the commercially important Australian greentail prawn (Metapenaeus bennettae) using transcriptomic and histological analyses. Prawns exposed to light, unweathered crude oil "spiked" sediments for 90 h were transferred to clean sediments for a further 72 h to assess recovery. Chemical analyses indicated that polycyclic aromatic hydrocarbons increased by approximately 65% and 91% in prawn muscle following 24 and 90 h of exposure, respectively, and significantly decreased during 24- and 72-h recovery periods. Transcriptomic responses followed an exposure and recovery pattern with innate immunity and nutrient metabolism transcripts significantly lowered in abundance after 24 h of exposure and were higher in abundance after 72 h of recovery. In addition, transcription/translation, cellular responses, and DNA repair pathways were significantly impacted after 24 h of exposure and recovered after 72 h of recovery. However, histological alterations such as tubule atrophy indicated an increase in severity after 24 and 72 h of recovery. The present study provides new insights into the sublethal impacts of crude oil exposure in greentail prawns and identifies molecular pathways altered by exposure. We expect these findings to inform future management associated with oil extraction activity and spills. Environ Toxicol Chem 2022;41:2162-2180. © 2022 John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Discovery Genome-Wide Association Study of Body Composition in 4,386 Adults From the UK Biobank's Pilot Imaging Enhancement Study.

Body composition (fat, skeletal muscle and bone mass) is an important determinant of overall health and risk of endocrine disorders such as type 2 diabetes and osteoporosis. Although diet and physical activity are strongly implicated, body composition is also heritable. We conducted a discovery genome-wide association study on 31 phenotypes from the three-compartment body composition model (fat, lean and bone mass) in a set of 4 386 individuals (n = 2 109 males, n = 2 294 females) from the UK Biobank pilot imaging enhancement program that underwent a dual energy X-ray absorptiometry (DXA) scan for assessment of body composition and genetic screening. From 6 137 607 imputed single nucleotide polymorphisms (SNPs) we identified 17 body composition loci (P<5.0 x 10-8). GWAS from the combined dataset identified four statistically significant SNPs (rs7592270, rs145972737, rs13212044, rs77772562). In sex-stratified GWAS, 10 male specific SNPs across all traits were identified and five female specific SNPs. Of the 17 SNPs, six were in or close to a gene where there was a plausible functional connection. Three SNPs (rs7592270, rs77772562 and rs7552312) were correlated with obesity phenotypes, one SNP (rs2236705) with lean phenotypes and two with bone mass phenotypes (rs112098641 and rs113380185). These results highlight candidate genes and biological pathways related to body composition, including glucose metabolism and estrogen regulation, that are of interest to replicate in future studies.

Dataset for genome sequencing and de novo assembly of the Vietnamese bighead catfish (Clarias macrocephalus Günther, 1864).

Freshwater catfish of the genus Clarias, known as the airbreathing catfish, are widespread and important for food security through small scale inland fisheries and aquaculture. Limited genomic data are available for this important group of fishes. The bighead catfish (Clarias macrocephalus) is a commercial aquaculture species in southeast Asia used for aquaculture and threatened in its natural environment through habitat destruction, over-exploitation and competition from other introduced species of Clarias. Despite its commercial importance and threats to natural populations, public databases do not include any genomic data for C. macrocephalus. We present the first genomic data for the bighead catfish from Illumina sequencing. A total of 128 Gb of sequence data in paired-end 150 bp reads were assembled de novo, generating a final assembly of 883 Mbp contained in 27,833 scaffolds (N50 length: 80.8 kbp) with BUSCO completeness assessments of 96.3% and 87.6% based on metazoan and Actinopterygii ortholog datasets, respectively. Annotation of the genome predicted 21,124 gene sequences, which were assigned putative functions based on homology to existing protein sequences in public databases. Raw fastq reads and the final version of the genome assembly have been deposited in the NCBI (BioProject: PRJNA604477, WGS: JAAGKR000000000, SRA: SRR11188453). The complete C. macrocephalus mitochondrial genome was also recovered from the same sequence read dataset and is available on NCBI (accession: MT109097), representing the first mitogenome for this species. Lastly, we find an expansion of the mb and ora1 genes thought to be associated with adaptations to air-breathing and a semi-terrestrial life style in this genus of catfish.

De Novo assembly and characterisation of the greentail prawn (Metapenaeus bennettae) hepatopancreas transcriptome - identification of stress response and detoxification transcripts.

Crude oil is a key contaminant in aquatic environments entering via natural and anthropogenic sources, causing toxicity in marine organisms. Traditionally, biomarkers have been utilised to determine crude oil exposure and effects in aquatic organisms, however advances in genomic technologies has led to increased adoption of transcriptomic approaches for identifying response and detoxification pathways following contaminant exposure. This study presents the first transcriptome for the greentail prawn (Metapenaeus bennettae), a commercially targeted benthic decapod crustacean from eastern and south-eastern Australia. The Trinity generated de novo assembly, after redundancy clustering, resulted in 86,401 contigs, of these 22,252 displayed strong homology to transcripts in the NCBI's non-redundant protein, Swiss-Prot and TrEMBL databases. Furthermore, Gene Ontology was assigned to 15,079 annotated contigs and KEGG Orthology was identified for 1318 annotated contigs. Transcripts encoding common biomarkers utilised to determine crude oil exposure were identified, including those for detoxification phase I and II enzymes; with 40 transcripts encoding for members of the cytochrome P450 gene family and 8 transcripts encoding glutathione S-Transferases (GSTs). Transcripts encoding oxidative stress enzymes including superoxide dismutase (SOD), catalase (CAT), glutathione peroxidase (GPx) and metallothionein (MT) were identified, as well as stress induced proteins including crustacean hyperglycemic hormone (CHH) and heat shock proteins (Hsps). The annotated transcriptome of the greentail prawn and the identification of detoxification and stress response transcripts, provides a necessary resource for future studies geared toward characterising differential transcriptomic patterns and molecular pathways after exposure to crude oil in this and other crustacean species of environmental and commercial importance.