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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 mitogenomics of the Decapoda reveals evolutionary heterogeneity in architecture and composition.

The emergence of cost-effective and rapid sequencing approaches has resulted in an exponential rise in the number of mitogenomes on public databases in recent years, providing greater opportunity for undertaking large-scale comparative genomic and systematic research. Nonetheless, current datasets predominately come from small and disconnected studies on a limited number of related species, introducing sampling biases and impeding research of broad taxonomic relevance. This study contributes 21 crustacean mitogenomes from several under-represented decapod infraorders including Polychelida and Stenopodidea, which are used in combination with 225 mitogenomes available on NCBI to investigate decapod mitogenome diversity and phylogeny. An overview of mitochondrial gene orders (MGOs) reveals a high level of genomic variability within the Decapoda, with a large number of MGOs deviating from the ancestral arthropod ground pattern and unevenly distributed among infraorders. Despite the substantial morphological and ecological variation among decapods, there was limited evidence for correlations between gene rearrangement events and species ecology or lineage specific nucleotide substitution rates. Within a phylogenetic context, predicted scenarios of rearrangements show some MGOs to be informative synapomorphies for some taxonomic groups providing strong independent support for phylogenetic relationships. Additional comparisons for a range of mitogenomic features including nucleotide composition, strand asymmetry, unassigned regions and codon usage indicate several clade-specific trends that are of evolutionary and ecological interest.

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.

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.