Functional dissection and assembly of a small, newly evolved, W chromosome-specific genomic region of the African clawed frog Xenopus laevis.

Genetic triggers for sex determination are frequently co-inherited with other linked genes that may also influence one or more sex-specific phenotypes. To better understand how sex-limited regions evolve and function, we studied a small W chromosome-specific region of the frog Xenopus laevis that contains only three genes (dm-w, scan-w, ccdc69-w) and that drives female differentiation. Using gene editing, we found that the sex-determining function of this region requires dm-w but that scan-w and ccdc69-w are not essential for viability, female development, or fertility. Analysis of mesonephros+gonad transcriptomes during sexual differentiation illustrates masculinization of the dm-w knockout transcriptome, and identifies mostly non-overlapping sets of differentially expressed genes in separate knockout lines for each of these three W-specific gene compared to wildtype sisters. Capture sequencing of almost all Xenopus species and PCR surveys indicate that the female-determining function of dm-w is present in only a subset of species that carry this gene. These findings map out a dynamic evolutionary history of a newly evolved W chromosome-specific genomic region, whose components have distinctive functions that frequently degraded during Xenopus diversification, and evidence the evolutionary consequences of recombination suppression.

Consequences of polyploidy and divergence as revealed by cytogenetic mapping of tandem repeats in African clawed frogs (Xenopus, Pipidae).

Repetitive elements have been identified in several amphibian genomes using whole genome sequencing, but few studies have used cytogenetic mapping to visualize these elements in this vertebrate group. Here, we used fluorescence in situ hybridization and genomic data to map the U1 and U2 small nuclear RNAs and histone H3 in six species of African clawed frog (genus Xenopus), including, from subgenus Silurana, the diploid Xenopus tropicalis and its close allotetraploid relative X. calcaratus and, from subgenus Xenopus, the allotetraploid species X. pygmaeus, X. allofraseri, X. laevis, and X. muelleri. Results allowed us to qualitatively evaluate the relative roles of polyploidization and divergence in the evolution of repetitive elements because our focal species include allotetraploid species derived from two independent polyploidization events - one that is relatively young that gave rise to X. calcaratus and another that is older that gave rise to the other (older) allotetraploids. Our results demonstrated conserved loci number and position of signals in the species from subgenus Silurana; allotetraploid X. calcaratus has twice as many signals as diploid X. tropicalis. However, the content of repeats varied among the other allotetraploid species. We detected almost same number of signals in X. muelleri as in X. calcaratus and same number of signals in X. pygmaeus, X. allofraseri, X. laevis as in the diploid X. tropicalis. Overall, these results are consistent with the proposal that allopolyploidization duplicated these tandem repeats and that variation in their copy number was accumulated over time through reduction and expansion in a subset of the older allopolyploids.

C. albicans UME7 deletion does not have major impacts on white opaque switching, filamentation, or virulence.

C. albicans is the most prevalent human fungal pathogen, and can be especially dangerous to immunocompromised individuals. One key aspect of C. albicans virulence is morphological plasticity. C. albicans can undergo a number of distinct morphological changes and these changes are controlled by complex transcriptional networks. The transcription factor Ume6 is an important member of these networks, playing an essential role mediating filamentation. C. albicans , however encodes a second UME6 homolog, UME7 . UME7 is highly conserved in the CTG fungal clade, but the role of UME7 in C. albicans biology is unknown. Here we truncate and delete C. albicans UME7 . We find Ume7 is dispensable for growth and filamentation. We also find that deletion does not have major consequences on virulence or white opaque switching. Our results suggest that under standard laboratory conditions deletion of UME7 does not have large effects on C. albicans phenotype leaving its role in C. albicans biology undefined.

Characterization of the role of Ume6 C-terminal tail in C. albicans morphological plasticity.

C. albicans is an important human fungal pathogen and filamentation is essential for its virulence. Ume6 is a transcription factor critical for filamentation. Ume6 is composed of three domains, a long N terminal domain, Zn-finger domain, and a C-terminal domain. Previously, it was shown that the Zn-finger domain is essential for filamentation, as removal of this domain led to a lack of filamentation. However, the role for the C-terminal domain has not been defined. We find deletion of the C-terminal domain leads to a filamentation defect and the defect is not as severe as removal of the Zn-finger or ume6 deletion. We mutated a number of residues in the C-terminal domain to try to identify specific residues important for filamentation, but all of our mutants displayed wild type filamentation. Alpha fold predictions suggest the C-terminal domain forms a single alpha helix that is predicted to interact with the Zn-finger domain via hydrogen bond. Our data suggests the C-terminal domain binds the Zn-finger domain and through this interaction is important for filamentation.

A data-driven framework for permeability prediction of natural porous rocks via microstructural characterization and pore-scale simulation.

Understanding the microstructure-property relationships of porous media is of great practical significance, based on which macroscopic physical properties can be directly derived from measurable microstructural informatics. However, establishing reliable microstructure-property mappings in an explicit manner is difficult, due to the intricacy, stochasticity, and heterogeneity of porous microstructures. In this paper, a data-driven computational framework is presented to investigate the inherent microstructure-permeability linkage for natural porous rocks, where multiple techniques are integrated together, including microscopy imaging, stochastic reconstruction, microstructural characterization, pore-scale simulation, feature selection, and data-driven modeling. A large number of 3D digital rocks with a wide porosity range are acquired from microscopy imaging and stochastic reconstruction techniques. A broad variety of morphological descriptors are used to quantitatively characterize pore microstructures from different perspectives, and they compose the raw feature pool for feature selection. High-fidelity lattice Boltzmann simulations are conducted to resolve fluid flow passing through porous media, from which reliable permeability references are obtained. The optimal feature set that best represents permeability is identified through a performance-oriented feature selection process, upon which a cost-effective surrogate model is rapidly fitted to approximate the microstructure-permeability mapping via data-driven modeling. This surrogate model exhibits great advantages over empirical/analytical formulas in terms of prediction accuracy and generalization capacity, which can predict reliable permeability values spanning four orders of magnitude. Besides, feature selection also greatly enhances the interpretability of the data-driven prediction model, from which new insights into the mechanism of how microstructural characteristics determine intrinsic permeability are obtained.

Large-scale functional screen identifies genetic variants with splicing effects in modern and archaic humans.

Humans coexisted and interbred with other hominins which later became extinct. These archaic hominins are known to us only through fossil records and for two cases, genome sequences. Here, we engineer Neanderthal and Denisovan sequences into thousands of artificial genes to reconstruct the pre-mRNA processing patterns of these extinct populations. Of the 5,169 alleles tested in this massively parallel splicing reporter assay (MaPSy), we report 962 exonic splicing mutations that correspond to differences in exon recognition between extant and extinct hominins. Using MaPSy splicing variants, predicted splicing variants, and splicing quantitative trait loci, we show that splice-disrupting variants experienced greater purifying selection in anatomically modern humans than that in Neanderthals. Adaptively introgressed variants were enriched for moderate-effect splicing variants, consistent with positive selection for alternative spliced alleles following introgression. As particularly compelling examples, we characterized a unique tissue-specific alternative splicing variant at the adaptively introgressed innate immunity gene TLR1, as well as a unique Neanderthal introgressed alternative splicing variant in the gene HSPG2 that encodes perlecan. We further identified potentially pathogenic splicing variants found only in Neanderthals and Denisovans in genes related to sperm maturation and immunity. Finally, we found splicing variants that may contribute to variation among modern humans in total bilirubin, balding, hemoglobin levels, and lung capacity. Our findings provide unique insights into natural selection acting on splicing in human evolution and demonstrate how functional assays can be used to identify candidate causal variants underlying differences in gene regulation and phenotype.

Species Delimitation, Phylogenomics, and Biogeography of Sulawesi Flying Lizards: A Diversification History Complicated by Ancient Hybridization, Cryptic Species, and Arrested Speciation.

The biota of Sulawesi is noted for its high degree of endemism and for its substantial levels of in situ biological diversification. While the island's long period of isolation and dynamic tectonic history have been implicated as drivers of the regional diversification, this has rarely been tested in the context of an explicit geological framework. Here, we provide a tectonically informed biogeographical framework that we use to explore the diversification history of Sulawesi flying lizards (the Draco lineatus Group), a radiation that is endemic to Sulawesi and its surrounding islands. We employ a framework for inferring cryptic speciation that involves phylogeographic and genetic clustering analyses as a means of identifying potential species followed by population demographic assessment of divergence-timing and rates of bi-directional migration as means of confirming lineage independence (and thus species status). Using this approach, phylogenetic and population genetic analyses of mitochondrial sequence data obtained for 613 samples, a 50-SNP data set for 370 samples, and a 1249-locus exon-capture data set for 106 samples indicate that the current taxonomy substantially understates the true number of Sulawesi Draco species, that both cryptic and arrested speciations have taken place, and that ancient hybridization confounds phylogenetic analyses that do not explicitly account for reticulation. The Draco lineatus Group appears to comprise 15 species-9 on Sulawesi proper and 6 on peripheral islands. The common ancestor of this group colonized Sulawesi ~11 Ma when proto-Sulawesi was likely composed of two ancestral islands, and began to radiate ~6 Ma as new islands formed and were colonized via overwater dispersal. The enlargement and amalgamation of many of these proto-islands into modern Sulawesi, especially during the past 3 Ma, set in motion dynamic species interactions as once-isolated lineages came into secondary contact, some of which resulted in lineage merger, and others surviving to the present. [Genomics; Indonesia; introgression; mitochondria; phylogenetics; phylogeography; population genetics; reptiles.].

Efficiency of lagoon-based municipal wastewater treatment in removing microplastics.

Municipal wastewater treatment plants act as a sink, but also are a source of microplastics in the environment. A conventional wastewater lagoon system and an activated sludge (AS)-lagoon system in Victoria (Australia) were investigated through a two-year sampling program to understand the fate and transport of MP in such treatment processes. The abundance (>25 µm) and characteristics (size, shape, and colour) of the microplastics present in the various wastewater streams were determined. The mean values of MP in the influent of the two plants were 55.3 ± 38.4 and 42.5 ± 20.1 MP/L, respectively. The dominant MP size of influent and final effluent was <500 µm, with 25-200 µm accounting for >65 % of the total MP; synthetic fibres were the dominant MP in all wastewater streams. Influent MP concentration was significantly higher in summer than in other seasons for both systems, which was related to the lower plant inflow due to less stormwater entering the sewer during summer. The promising MP removal capability of the lagoon system (97 %) was attributed to its lengthy wastewater detention time (total HRT >250 days, including the storage lagoons) that would allow effective separation of MP from the water column via various physical and biological pathways. For the AS-lagoon system, the high MP reduction efficiency (98.4 %) was attributed to the post-secondary treatment of the wastewater with the lagoon system, in which MP was further removed during the month-long detention in the lagoons. The results indicated the potential of such low-energy and low-cost wastewater treatment systems for MP control.

Mitonuclear Interactions and the Origin of Macaque Societies.

In most eukaryotes, aerobic respiration requires interactions between autosomally encoded genes (Ninteract genes) and mitochondrial DNA, RNA, and protein. In species where females are philopatric, contrasting distributions of genetic variation in mitochondrial and nuclear genomes create variation in mitonuclear interactions that may be subject to natural selection. To test this expectation, we turned to a group with extreme female philopatry: the macaque monkeys. We examined four genomic data sets from (1) wild caught and (2) captive populations of rhesus macaque, which is the most widely distributed nonhuman primate, and (3) the stump-tailed macaque and (4) a subspecies of longtail macaque, both of whose mitochondrial DNA is introgressed from a highly diverged ancestor. We identified atypically long runs of homozygosity, low polymorphism, high differentiation, and/or rapid protein evolution associated with Ninteract genes compared with non-Ninteract genes. These metrics suggest a subset of Ninteract genes were independently subject to atypically pervasive natural selection in multiple species. These findings suggest that natural selection on mitonuclear interactions could have influenced several aspects of macaque societies including species diversity, ecological breadth, female-biased adult sex ratio and demography, sexual dimorphism, and mitonuclear phylogenomics.

Investigating the fate and transport of microplastics in a lagoon wastewater treatment system using a multimedia model approach.

Effluents of municipal wastewater treatment plants (WWTPs) are a major source of microplastics (MP) in the terrestrial and aquatic environment; there is growing concern over the environmental and health impacts of MP pollution. In this study, the MP removal (MP cut-off size= 25 µm) in a lagoon-based wastewater treatment system was predicted by developing a model based on the multimedia modelling approach and utilising MP-specific properties for improving the understanding of the fate and transport of MP in such treatment processes. The high MP removal efficiency of the lagoon treatment system as predicted by the model (99.3%) and determined with the site wastewater samples (97%) could be attributed to its high HRT (>200 days, including that for the storage lagoons) that would allow effective MP removal with the system. Evaluation of the model predictions of MP concentration demonstrated reasonable alignment with measured concentrations in the facultative, maturation and winter storage lagoons of the system. Further evaluation of model predictions for various MP size classes (25-100, 100-200, 200-500 and >500 µm) obtained reasonable predictions for MP within the size range of 25-500 µm, indicating that the model is better used for predicting MP within that size range. The sensitivity analysis revealed the model predictions to be sensitive towards the operating/water quality parameters in the order of influent wastewater flowrate, MP concentration in influent wastewater, and MP settling rate in the water column of the lagoon. The study showed the potential of the developed model as a quantitative assessment tool for better management of MP in lagoon-based WWTPs.

Population genomics and subgenome evolution of the allotetraploid frog Xenopus laevis in southern Africa.

Allotetraploid genomes have two distinct genomic components called subgenomes that are derived from separate diploid ancestral species. Many genomic characteristics such as gene function, expression, recombination, and transposable element mobility may differ significantly between subgenomes. To explore the possibility that subgenome population structure and gene flow may differ as well, we examined genetic variation in an allotetraploid frog-the African clawed frog (Xenopus laevis)-over the dynamic and varied habitat of its native range in southern Africa. Using reduced representation genome sequences from 91 samples from 12 localities, we found no strong evidence that population structure and gene flow differed substantially by subgenome. We then compared patterns of population structure in the nuclear genome to the mitochondrial genome using Sanger sequences from 455 samples from 183 localities. Our results provide further resolution to the geographic distribution of mitochondrial and nuclear diversity in this species and illustrate that population structure in both genomes corresponds roughly with variation in seasonal rainfall and with the topography of southern Africa.

Divergent subgenome evolution in the allotetraploid frog Xenopus calcaratus.

Allopolyploid genomes are divided into compartments called subgenomes that are derived from lower ploidy ancestors. In African clawed frogs of the subgenus Xenopus (genus Xenopus), allotetraploid species have two subgenomes (L and S) with morphologically distinct homoeologous chromosomes. In allotetraploid species of the sister subgenus Silurana, independently evolved subgenomes also exist, but their cytogenetics has not been investigated in detail. We used a diverse suite of cytogenetic and molecular FISH techniques on an allotetraploid species in Silurana-Xenopus calcaratus-to explore evolutionary dynamics of chromosome morphology and rearrangements. We find that the subgenomes of X. calcaratus have distinctive characteristics, with a more conserved a-subgenome resembling the closely related genome of the diploid species X. tropicalis, and a more rapidly evolving b-subgenome having more pronounced changes in chromosome structure, including diverged heterochromatic blocks, repetitive sequences, and deletion of a nucleolar secondary constriction. Based on these cytogenetic differences, we propose a chromosome nomenclature for X. calcaratus that may apply to other allotetraploids in subgenus Silurana, depending on as yet unresolved details of their evolutionary origins. These findings highlight the potential for large-scale asymmetry in subgenome evolution following allopolyploidization.

New insights into Xenopus sex chromosome genomics from the Marsabit clawed frog X. borealis.

In many groups, sex chromosomes change frequently but the drivers of their rapid evolution are varied and often poorly characterized. With an aim of further understanding sex chromosome turnover, we investigated the polymorphic sex chromosomes of the Marsabit clawed frog, Xenopus borealis, using genomic data and a new chromosome-scale genome assembly. We confirmed previous findings that 54.1 Mb of chromosome 8L is sex-linked in animals from east Kenya and a laboratory strain, but most (or all) of this region is not sex-linked in natural populations from west Kenya. Previous work suggests possible degeneration of the Z chromosomes in the east population because many sex-linked transcripts of this female heterogametic population have female-biased expression, and we therefore expected this chromosome to not be present in the west population. In contrast, our simulations support a model where most or all of the sex-linked portion of the Z chromosome from the east acquired autosomal segregation in the west, and where much genetic variation specific to the large sex-linked portion of the W chromosome from the east is not present in the west. These recent changes are consistent with the hot-potato model, wherein sex chromosome turnover is favoured by natural selection if it purges a (minimally) degenerate sex-specific sex chromosome, but counterintuitively suggest natural selection failed to purge a Z chromosome that has signs of more advanced and possibly more ancient regulatory degeneration. These findings highlight complex evolutionary dynamics of young, rapidly evolving Xenopus sex chromosomes and set the stage for mechanistic work aimed at pinpointing additional sex-determining genes in this group.

CRISPR-Mediated Genome Editing in the Human Fungal Pathogen C. albicans.

Cas9-mediated genome editing is one tool investigators can use to study fungal pathogens. Such methodologies allow the investigator to examine how fungal cells differ from human cells and thus potentially identify novel therapeutic targets. In this chapter, we describe how CRISPR-mediated genome editing can be used to edit the genome of the most prevalent human fungal pathogen C. albicans. A cassette encoding a fungal optimized Cas9 nuclease and guide RNA is integrated into the C. albicans genome. The guide RNA targets Cas9 to the complementary genome sequence, and Cas9 cleaves the DNA. A repair template encoding whatever changes the investigator wished to make to the genome is co-transformed with the cassette and repairs the break via homologous recombination, thus introducing the change to the genome. The method we describe enables the researcher to edit the C. albicans genome and then efficiently remove the editing machinery and antibiotic resistance markers. This allows one to sequentially edit the C. albicans genome when multiple changes are desired. In addition, we provide notes that provide guidance on how the described protocols can be altered to meet the demands of the researcher. In these notes, we also describe the recent development of a more flexible CRISPR system that has a relaxed PAM site specificity. These and other advancements make CRISPR-mediated genome editing a practical approach when one needs to genetically alter C. albicans.

Evolution of genes involved in the unusual genitals of the bear macaque, Macaca arctoides.

Genital divergence is thought to contribute to reproductive barriers by establishing a "lock-and-key" mechanism for reproductive compatibility. One such example, Macaca arctoides, the bear macaque, has compensatory changes in both male and female genital morphology as compared to close relatives. M. arctoides also has a complex evolutionary history, having extensive introgression between the fascicularis and sinica macaque species groups. Here, phylogenetic relationships were analyzed via whole-genome sequences from five species, including M. arctoides, and two species each from the putative parental species groups. This analysis revealed ~3x more genomic regions supported placement in the sinica species group as compared to the fascicularis species group. Additionally, introgression analysis of the M. arctoides genome revealed it is a mosaic of recent polymorphisms shared with both species groups. To examine the evolution of their unique genital morphology further, the prevalence of candidate genes involved in genital morphology was compared against genome-wide outliers in various population genetic metrics of diversity, divergence, introgression, and selection, while accounting for background variation in recombination rate. This analysis identified 67 outlier genes, including several genes that influence baculum morphology in mice, which were of interest since the bear macaque has the longest primate baculum. The mean of four of the seven population genetic metrics was statistically different in the candidate genes as compared to the rest of the genome, suggesting that genes involved in genital morphology have increased divergence and decreased diversity beyond expectations. These results highlight specific genes that may have played a role in shaping the unique genital morphology in the bear macaque.

Neofunctionalization of a Noncoding Portion of a DNA Transposon in the Coding Region of the Chimerical Sex-Determining Gene dm-W in Xenopus Frogs.

Most vertebrate sex-determining genes (SDGs) emerge as neofunctionalized genes through duplication and/or mutation of ancestral genes that are involved with sexual differentiation. We previously demonstrated dm-W to be the SDG in the African clawed frog Xenopus laevis and found that a portion of this gene emerged from the masculinization gene dmrt1 after allotetraploidization by interspecific hybridization between two ancestral species around 17-18 Ma. dm-W has four exons consisting of a noncoding exon 1, dmrt1-derived exons 2 and 3, and an orphan exon 4 (Ex4) of unknown origin that includes coding sequence (CDS). In this study, we searched for the origin of Ex4 and investigated the function of the CDS of this exon. We found that the Ex4-CDS is derived from a noncoding portion of the hAT-10 family of DNA transposon. Evolutionary analysis of transposons and determination of the Ex4 sequences from three other species indicated that Ex4 was generated before the diversification of most or all extant allotetraploid species in subgenus Xenopus, during which time we hypothesize that transposase activity of this hAT superfamily was active. Using DNA-protein binding and transfection assays, we further demonstrate that the Ex4-encoded amino acid sequence increases the DNA-binding ability and transrepression activity of DM-W. These findings suggest that the conversion of the noncoding transposon sequence to the CDS of dm-W contributed to neofunctionalization of a new chimeric SDG in the ancestor of the allotetraploid Xenopus species, offering new insights into de novo origin and functional evolution of chimerical genes.

Stay-on-Task Exercises as a Tool To Maintain Focus during a CRISPR CURE.

Course-based undergraduate research experiences (CURE) offer the chance for students to experience authentic research investigation in a classroom setting. Such hands-on experiences afford unique opportunities work on a semi-independent research project in an efficient, structured environment. We have developed a CRISPR CURE in which undergraduate and graduate students use in silico, in vitro, and in vivo techniques to edit a fungal genome. During the development of this course, we have found that the asynchronous nature of the CRISPR CURE activities can be disruptive and lead to unproductive class time. To overcome this challenge, we have developed stay-on-task exercises (SOTEs). These short low-stakes assessments provide structured activities that are performed during these asynchronous incubation periods. SOTE activities leverage potentially unproductive class time and complement the CURE learning objectives. We have found SOTEs to be one method of maintaining classroom structure during a CURE. Furthermore, SOTE complexity, length, and subject can be easily modified to match course learning objectives.

Mitonuclear interactions and introgression genomics of macaque monkeys (Macaca) highlight the influence of behaviour on genome evolution.

In most macaques, females are philopatric and males migrate from their natal ranges, which results in pronounced divergence of mitochondrial genomes within and among species. We therefore predicted that some nuclear genes would have to acquire compensatory mutations to preserve compatibility with diverged interaction partners from the mitochondria. We additionally expected that these sex-differences would have distinctive effects on gene flow in the X and autosomes. Using new genomic data from 29 individuals from eight species of Southeast Asian macaque, we identified evidence of natural selection associated with mitonuclear interactions, including extreme outliers of interspecies differentiation and metrics of positive selection, low intraspecies polymorphism and atypically long runs of homozygosity associated with nuclear-encoded genes that interact with mitochondria-encoded genes. In one individual with introgressed mitochondria, we detected a small but significant enrichment of autosomal introgression blocks from the source species of her mitochondria that contained genes which interact with mitochondria-encoded loci. Our analyses also demonstrate that sex-specific demography sculpts genetic exchange across multiple species boundaries. These findings show that behaviour can have profound but indirect effects on genome evolution by influencing how interacting components of different genomic compartments (mitochondria, the autosomes and the sex chromosomes) move through time and space.