Longitudinal Analysis of CCR5 and CXCR4 Usage in a Cohort of Antiretroviral Therapy-Na?ve Subjects with Progressive HIV-1 Subtype C Infection

HIV-1 subtype C (C-HIV) is responsible for most HIV-1 cases worldwide. Although the pathogenesis of C-HIV is thought to predominantly involve CCR5-restricted (R5) strains, we do not have a firm understanding of how frequently CXCR4-using (X4 and R5X4) variants emerge in subjects with progressive C-HIV infection. Nor do we completely understand the molecular determinants of coreceptor switching by C-HIV variants. Here, we characterized a panel of HIV-1 envelope glycoproteins (Envs) (n = 300) cloned sequentially from plasma of 21 antiretroviral therapy (ART)-naïve subjects who experienced progression from chronic to advanced stages of C-HIV infection, and show that CXCR4-using C-HIV variants emerged in only one individual. Mutagenesis studies and structural models suggest that the evolution of R5 to X4 variants in this subject principally involved acquisition of an “Ile-Gly” insertion in the gp120 V3 loop and replacement of the V3 “Gly-Pro-Gly” crown with a “Gly-Arg-Gly” motif, but that the accumulation of additional gp120 “scaffold” mutations was required for these V3 loop changes to confer functional effects. In this context, either of the V3 loop changes could confer possible transitional R5X4 phenotypes, but when present together they completely abolished CCR5 usage and conferred the X4 phenotype. Our results show that the emergence of CXCR4-using strains is rare in this cohort of untreated individuals with advanced C-HIV infection. In the subject where X4 variants did emerge, alterations in the gp120 V3 loop were necessary but not sufficient to confer CXCR4 usage.     

The Effect of Genomic Inversions on Estimation of Population Genetic Parameters from SNP Data

In recent years it has emerged that structural variants have a substantial impact on genomic variation. Inversion polymorphisms represent a significant class of structural variant, and despite the challenges in their detection, data on inversions in the human genome are increasing rapidly. Statistical methods for inferring parameters such as the recombination rate and the selection coefficient have generally been developed without accounting for the presence of inversions. Here we exploit new software for simulating inversions in population genetic data, invertFREGENE, to assess the potential impact of inversions on such methods. Using data simulated by invertFREGENE, as well as real data from several sources, we test whether large inversions have a disruptive effect on widely applied population genetics methods for inferring recombination rates, for detecting selection, and for controlling for population structure in genome-wide association studies (GWAS). We find that recombination rates estimated by LDhat are biased downward at inversion loci relative to the true contemporary recombination rates at the loci but that recombination hotspots are not falsely inferred at inversion breakpoints as may have been expected. We find that the integrated haplotype score (iHS) method for detecting selection appears robust to the presence of inversions. Finally, we observe a strong bias in the genome-wide results of principal components analysis (PCA), used to control for population structure in GWAS, in the presence of even a single large inversion, confirming the necessity to thin SNPs by linkage disequilibrium at large physical distances to obtain unbiased results.     

A computational framework to investigate retinal haemodynamics and tissue stress

Biomechanics and Modeling in Mechanobiology - The process of vision begins in the retina, yet the role of biomechanical forces in the retina is relatively unknown and only recently being explored....     

Phylodynamic Inference of Bacterial Outbreak Parameters Using Nanopore Sequencing

Nanopore sequencing and phylodynamic modeling have been used to reconstruct the transmission dynamics of viral epidemics, but their application to bacterial pathogens has remained challenging. Cost-effective bacterial genome sequencing and variant calling on nanopore platforms would greatly enhance surveillance and outbreak response in communities without access to sequencing infrastructure. Here, we adapt random forest models for single nucleotide polymorphism (SNP) polishing developed by Sanderson and colleagues (2020. High precision Neisseria gonorrhoeae variant and antimicrobial resistance calling from metagenomic nanopore sequencing. Genome Res. 30(9):1354–1363) to estimate divergence and effective reproduction numbers (R_e) of two methicillin-resistant Staphylococcus aureus (MRSA) outbreaks from remote communities in Far North Queensland and Papua New Guinea (PNG; n = 159). Successive barcoded panels of S. aureus isolates (2 × 12 per MinION) sequenced at low coverage (>5× to 10×) provided sufficient data to accurately infer genotypes with high recall when compared with Illumina references. Random forest models achieved high resolution on ST93 outbreak sequence types (>90% accuracy and precision) and enabled phylodynamic inference of epidemiological parameters using birth–death skyline models. Our method reproduced phylogenetic topology, origin of the outbreaks, and indications of epidemic growth (R_e > 1). Nextflow pipelines implement SNP polisher training, evaluation, and outbreak alignments, enabling reconstruction of within-lineage transmission dynamics for infection control of bacterial disease outbreaks on portable nanopore platforms. Our study shows that nanopore technology can be used for bacterial outbreak reconstruction at competitive costs, providing opportunities for infection control in hospitals and communities without access to sequencing infrastructure, such as in remote northern Australia and PNG.     

Enhanced protein domain discovery using taxonomy

Background  

It is well known that different species have different protein domain repertoires, and indeed that some protein domains are kingdom specific. This information has not yet been incorporated into statistical methods for finding domains in sequences of amino acids.     

In addition to the SH3 binding region, multiple regions within the N-terminal noncatalytic portion of the cAMP-specific phosphodiesterase, PDE4A5, contribute to its intracellular targeting

The long cyclic AMP (cAMP)-specific phosphodiesterase isoform, PDE4A5 (PDE4A subfamily isoform variant 5), when transiently expressed in COS-7 cells, was shown in subcellular fractionation studies to be associated with both membrane and cytosol fractions, with immunofluorescence analyses identifying PDE4A5 as associated both with ruffles at the cell margin and also at a distinct perinuclear localisation. Deletion of the first nine amino acids of PDE4A5 (1) ablated its ability to interact with the SH3 domain of the tyrosyl kinase, LYN; (2) reduced, but did not ablate, membrane association; and (3) disrupted the focus of PDE4A5 localisation within ruffles at the cell margin. This deleted region contained a Class I SH3 binding motif of similar sequence to those identified by screening a phage display library with the LYN-SH3 domain. Truncation to remove the PDE4A5 isoform-specific N-terminal region caused a further reduction in membrane association and ablated localisation at the cell margin. Progressive truncation to delete the PDE4A long isoform common region and then the long isoform-specific UCR1 did not cause any further change in membrane association or intracellular distribution. However, deletion up to the super-short form splice junction generated an entirely soluble 'core' PDE4A species. We propose that multiple sites in the N-terminal noncatalytic portion of PDE4A5 have the potential to associate with intracellular structures and thus define its intracellular localisation. At least two such sites lie within the PDE4A5 isoform-specific N-terminal region and these appear to be primarily responsible for targeting PDE4A5 to, and organising it within, the cell margin; one is an SH3 binding motif able to interact with LYN kinase and the other lies within the C-terminal portion of the PDE4A5 unique region. A third membrane association region is located within the N-terminal portion of UCR2 and appears to be primarily responsible for targeting to the perinuclear region. Progressive N-terminal truncation, to delete defined regions of PDE4A5, identified activity changes occurring upon deletion of the SH3 binding site region and then upon deletion of the membrane association site region located within UCR2. This suggests that certain of these anchor sites may not only determine intracellular targeting but may also transduce regulatory effects on PDE4A5 activity.     

Regeneration of halved embryonic lower first mouse molars: correlation with the distribution pattern of non dividing IDE cells, the putative organizers of morphogenetic units, the cusps

Recently we demonstrated that non-cycling, cap-stage, mouse molar inner dental epithelial (IDE) cells corresponding to the primary enamel knot (EK) area underwent a coordinated temporo-spatial patterning leading to their patchy irregular segregation at the tips of the forming cusps. These non-cycling cells were suggested to perhaps represent the organizers of the morphogenetic units (OMU), the cusps. The present study has analyzed the regenerative capacity of halved cap-stage first lower mouse molars through three dimensional (3D) reconstructions. Partial regeneration of the anterior half and possible complete regeneration of the posterior half were documented. Using BrdU (5-bromo-2'-deoxyuridine) labeling and 3D reconstructions of the IDE, we have correlated the patterns of cusp regeneration with the distribution of BrdU negative IDE cells. These data support a morphogenetic role for the non-cycling IDE cells.