Deep learning of the retina enables phenome- and genome-wide analyses of the microvasculature.

Background: The microvasculature, the smallest blood vessels in the body, has key roles in maintenance of organ health as well as tumorigenesis. The retinal fundus is a window for human in vivo non-invasive assessment of the microvasculature. Large-scale complementary machine learning-based assessment of the retinal vasculature with phenome-wide and genome-wide analyses may yield new insights into human health and disease. Methods: We utilized 97,895 retinal fundus images from 54,813 UK Biobank participants. Using convolutional neural networks to segment the retinal microvasculature, we calculated fractal dimension (FD) as a measure of vascular branching complexity, and vascular density. We associated these indices with 1,866 incident ICD-based conditions (median 10y follow-up) and 88 quantitative traits, adjusting for age, sex, smoking status, and ethnicity. Results: Low retinal vascular FD and density were significantly associated with higher risks for incident mortality, hypertension, congestive heart failure, renal failure, type 2 diabetes, sleep apnea, anemia, and multiple ocular conditions, as well as corresponding quantitative traits. Genome-wide association of vascular FD and density identified 7 and 13 novel loci respectively, which were enriched for pathways linked to angiogenesis (e.g., VEGF, PDGFR, angiopoietin, and WNT signaling pathways) and inflammation (e.g., interleukin, cytokine signaling). Conclusions: Our results indicate that the retinal vasculature may serve as a biomarker for future cardiometabolic and ocular disease and provide insights on genes and biological pathways influencing microvascular indices. Moreover, such a framework highlights how deep learning of images can quantify an interpretable phenotype for integration with electronic health records, biomarker, and genetic data to inform risk prediction and risk modification

Molecular Phylogeny and Evolution of Parabasalia with Improved Taxon Sampling and New Protein Markers of Actin and Elongation Factor-1α

BACKGROUND: Inferring the evolutionary history of phylogenetically isolated, deep-branching groups of taxa-in particular determining the root-is often extraordinarily difficult because their close relatives are unavailable as suitable outgroups. One of these taxonomic groups is the phylum Parabasalia, which comprises morphologically diverse species of flagellated protists of ecological, medical, and evolutionary significance. Indeed, previous molecular phylogenetic analyses of members of this phylum have yielded conflicting and possibly erroneous inferences. Furthermore, many species of Parabasalia are symbionts in the gut of termites and cockroaches or parasites and therefore formidably difficult to cultivate, rendering available data insufficient. Increasing the numbers of examined taxa and informative characters (e.g., genes) is likely to produce more reliable inferences. PRINCIPAL FINDINGS: Actin and elongation factor-1α genes were identified newly from 22 species of termite-gut symbionts through careful manipulations and seven cultured species, which covered major lineages of Parabasalia. Their protein sequences were concatenated and analyzed with sequences of previously and newly identified glyceraldehyde-3-phosphate dehydrogenase and the small-subunit rRNA gene. This concatenated dataset provided more robust phylogenetic relationships among major groups of Parabasalia and a more plausible new root position than those previously reported. CONCLUSIONS/SIGNIFICANCE: We conclude that increasing the number of sampled taxa as well as the addition of new sequences greatly improves the accuracy and robustness of the phylogenetic inference. A morphologically simple cell is likely the ancient form in Parabasalia as opposed to a cell with elaborate flagellar and cytoskeletal structures, which was defined as most basal in previous inferences. Nevertheless, the evolution of Parabasalia is complex owing to several independent multiplication and simplification events in these structures. Therefore, systematics based solely on morphology does not reflect the evolutionary history of parabasalids

Meiotic Recombination Intermediates Are Resolved with Minimal Crossover Formation during Return-to-Growth, an Analogue of the Mitotic Cell Cycle

Accurate segregation of homologous chromosomes of different parental origin (homologs) during the first division of meiosis (meiosis I) requires inter-homolog crossovers (COs). These are produced at the end of meiosis I prophase, when recombination intermediates that contain Holliday junctions (joint molecules, JMs) are resolved, predominantly as COs. JM resolution during the mitotic cell cycle is less well understood, mainly due to low levels of inter-homolog JMs. To compare JM resolution during meiosis and the mitotic cell cycle, we used a unique feature of Saccharomyces cerevisiae, return to growth (RTG), where cells undergoing meiosis can be returned to the mitotic cell cycle by a nutritional shift. By performing RTG with ndt80 mutants, which arrest in meiosis I prophase with high levels of interhomolog JMs, we could readily monitor JM resolution during the first cell division of RTG genetically and, for the first time, at the molecular level. In contrast to meiosis, where most JMs resolve as COs, most JMs were resolved during the first 1.5–2 hr after RTG without producing COs. Subsequent resolution of the remaining JMs produced COs, and this CO production required the Mus81/Mms4 structure-selective endonuclease. RTG in sgs1-ΔC795 mutants, which lack the helicase and Holliday junction-binding domains of this BLM homolog, led to a substantial delay in JM resolution; and subsequent JM resolution produced both COs and NCOs. Based on these findings, we suggest that most JMs are resolved during the mitotic cell cycle by dissolution, an Sgs1 helicase-dependent process that produces only NCOs. JMs that escape dissolution are mostly resolved by Mus81/Mms4-dependent cleavage that produces both COs and NCOs in a relatively unbiased manner. Thus, in contrast to meiosis, where JM resolution is heavily biased towards COs, JM resolution during RTG minimizes CO formation, thus maintaining genome integrity and minimizing loss of heterozygosity

Gene–Environment Interactions at Nucleotide Resolution

Interactions among genes and the environment are a common source of phenotypic variation. To characterize the interplay between genetics and the environment at single nucleotide resolution, we quantified the genetic and environmental interactions of four quantitative trait nucleotides (QTN) that govern yeast sporulation efficiency. We first constructed a panel of strains that together carry all 32 possible combinations of the 4 QTN genotypes in 2 distinct genetic backgrounds. We then measured the sporulation efficiencies of these 32 strains across 8 controlled environments. This dataset shows that variation in sporulation efficiency is shaped largely by genetic and environmental interactions. We find clear examples of QTN:environment, QTN: background, and environment:background interactions. However, we find no QTN:QTN interactions that occur consistently across the entire dataset. Instead, interactions between QTN only occur under specific combinations of environment and genetic background. Thus, what might appear to be a QTN:QTN interaction in one background and environment becomes a more complex QTN:QTN:environment:background interaction when we consider the entire dataset as a whole. As a result, the phenotypic impact of a set of QTN alleles cannot be predicted from genotype alone. Our results instead demonstrate that the effects of QTN and their interactions are inextricably linked both to genetic background and to environmental variation

RNA Methylation by the MIS Complex Regulates a Cell Fate Decision in Yeast

For the yeast Saccharomyces cerevisiae, nutrient limitation is a key developmental signal causing diploid cells to switch from yeast-form budding to either foraging pseudohyphal (PH) growth or meiosis and sporulation. Prolonged starvation leads to lineage restriction, such that cells exiting meiotic prophase are committed to complete sporulation even if nutrients are restored. Here, we have identified an earlier commitment point in the starvation program. After this point, cells, returned to nutrient-rich medium, entered a form of synchronous PH development that was morphologically and genetically indistinguishable from starvation-induced PH growth. We show that lineage restriction during this time was, in part, dependent on the mRNA methyltransferase activity of Ime4, which played separable roles in meiotic induction and suppression of the PH program. Normal levels of meiotic mRNA methylation required the catalytic domain of Ime4, as well as two meiotic proteins, Mum2 and Slz1, which interacted and co-immunoprecipitated with Ime4. This MIS complex (Mum2, Ime4, and Slz1) functioned in both starvation pathways. Together, our results support the notion that the yeast starvation response is an extended process that progressively restricts cell fate and reveal a broad role of post-transcriptional RNA methylation in these decisions

Phylogeny of Parasitic Parabasalia and Free-Living Relatives Inferred from Conventional Markers vs. Rpb1, a Single-Copy Gene

Parabasalia are single-celled eukaryotes (protists) that are mainly comprised of endosymbionts of termites and wood roaches, intestinal commensals, human or veterinary parasites, and free-living species. Phylogenetic comparisons of parabasalids are typically based upon morphological characters and 18S ribosomal RNA gene sequence data (rDNA), while biochemical or molecular studies of parabasalids are limited to a few axenically cultivable parasites. These previous analyses and other studies based on PCR amplification of duplicated protein-coding genes are unable to fully resolve the evolutionary relationships of parabasalids. As a result, genetic studies of Parabasalia lag behind other organisms.Comparing parabasalid EF1α, α-tubulin, enolase and MDH protein-coding genes with information from the Trichomonas vaginalis genome reveals difficulty in resolving the history of species or isolates apart from duplicated genes. A conserved single-copy gene encodes the largest subunit of RNA polymerase II (Rpb1) in T. vaginalis and other eukaryotes. Here we directly sequenced Rpb1 degenerate PCR products from 10 parabasalid genera, including several T. vaginalis isolates and avian isolates, and compared these data by phylogenetic analyses. Rpb1 genes from parabasalids, diplomonads, Parabodo, Diplonema and Percolomonas were all intronless, unlike intron-rich homologs in Naegleria, Jakoba and Malawimonas.The phylogeny of Rpb1 from parasitic and free-living parabasalids, and conserved Rpb1 insertions, support Trichomonadea, Tritrichomonadea, and Hypotrichomonadea as monophyletic groups. These results are consistent with prior analyses of rDNA and GAPDH sequences and ultrastructural data. The Rpb1 phylogenetic tree also resolves species- and isolate-level relationships. These findings, together with the relative ease of Rpb1 isolation, make it an attractive tool for evaluating more extensive relationships within Parabasalia

Optimizing the Design of Oligonucleotides for Homology Directed Gene Targeting

BACKGROUND: Gene targeting depends on the ability of cells to use homologous recombination to integrate exogenous DNA into their own genome. A robust mechanistic model of homologous recombination is necessary to fully exploit gene targeting for therapeutic benefit. METHODOLOGY/PRINCIPAL FINDINGS: In this work, our recently developed numerical simulation model for homology search is employed to develop rules for the design of oligonucleotides used in gene targeting. A Metropolis Monte-Carlo algorithm is used to predict the pairing dynamics of an oligonucleotide with the target double-stranded DNA. The model calculates the base-alignment between a long, target double-stranded DNA and a probe nucleoprotein filament comprised of homologous recombination proteins (Rad51 or RecA) polymerized on a single strand DNA. In this study, we considered different sizes of oligonucleotides containing 1 or 3 base heterologies with the target; different positions on the probe were tested to investigate the effect of the mismatch position on the pairing dynamics and stability. We show that the optimal design is a compromise between the mean time to reach a perfect alignment between the two molecules and the stability of the complex. CONCLUSION AND SIGNIFICANCE: A single heterology can be placed anywhere without significantly affecting the stability of the triplex. In the case of three consecutive heterologies, our modeling recommends using long oligonucleotides (at least 35 bases) in which the heterologous sequences are positioned at an intermediate position. Oligonucleotides should not contain more than 10% consecutive heterologies to guarantee a stable pairing with the target dsDNA. Theoretical modeling cannot replace experiments, but we believe that our model can considerably accelerate optimization of oligonucleotides for gene therapy by predicting their pairing dynamics with the target dsDNA

Recombinant α-actinin subunit antigens of Trichomonas vaginalis as potential vaccine candidates in protecting against trichomoniasis

BACKGROUND: Human trichomoniasis caused by Trichomonas vaginalis is one of the most common sexually transmitted diseases with more than 200 million cases worldwide. It has caused a series of health problems to patients. For prevention and control of infectious diseases, vaccines are usually considered as one of the most cost-efficient tools. However, until now, work on the development of T. vaginalis vaccines is still mainly focused on the screening of potential immunogens. Alpha-actinin characterized by high immunogenicity in T. vaginalis was suggested as a promising candidate. Therefore, the purpose of this study was to evaluate the protective potency of recombinant α-actinin against T. vaginalis infection in a mouse intraperitoneal model. METHODS: Two selected coding regions of α-actinin (ACT-F, 14-469 aa and ACT-T, 462-844 aa) amplified from cDNA were cloned into pET-32a (+) expression vector and transfected into BL21 cells. After induction with IPTG and purification with electroelution, the two recombinant fusion proteins were emulsified in Freund's adjuvant (FA) and used to immunize BALB/C mice. Following intraperitoneal inoculation with T. vaginalis, the survival rate of mice was monitored for the assessment of protective potency. After immunization, the antibody level in mouse serum was assessed by ELISA, splenocyte proliferation response was detected with CCK8 and cytokines in the supernatant of splenocytes were quantified with a cytometric bead-based assay. RESULTS: We successfully obtained purified ACT-F (70.33 kDa) and ACT-T (61.7kDa). Both recombinant proteins could provide significant protection against T. vaginalis challenge, especially ACT-T (with 100% protection within one month). Meanwhile, high levels of specific total IgG and subtypes (IgG1 > IgG2a) were detected in sera from the immunized mice. Our results also revealed a statistically significant increase in splenocyte proliferation and related cytokine (IFN-γ, IL-6, IL-17A and IL-10) production after repeated stimulation with the corresponding antigens in vitro. CONCLUSIONS: Immunization with both ACT-F and ACT-T could confer partial to complete protection and trigger strong Th1/Th2 mixed humoral and cellular immune responses in the mouse host. This suggested that recombinant α-actinin subunit antigens may be promising vaccine candidates against trichomoniasis

Srf1 Is a Novel Regulator of Phospholipase D Activity and Is Essential to Buffer the Toxic Effects of C16:0 Platelet Activating Factor

During Alzheimer's Disease, sustained exposure to amyloid-β42 oligomers perturbs metabolism of ether-linked glycerophospholipids defined by a saturated 16 carbon chain at the sn-1 position. The intraneuronal accumulation of 1-O-hexadecyl-2-acetyl-sn-glycerophosphocholine (C16:0 PAF), but not its immediate precursor 1-O-hexadecyl-sn-glycerophosphocholine (C16:0 lyso-PAF), participates in signaling tau hyperphosphorylation and compromises neuronal viability. As C16:0 PAF is a naturally occurring lipid involved in cellular signaling, it is likely that mechanisms exist to protect cells against its toxic effects. Here, we utilized a chemical genomic approach to identify key processes specific for regulating the sensitivity of Saccharomyces cerevisiae to alkyacylglycerophosphocholines elevated in Alzheimer's Disease. We identified ten deletion mutants that were hypersensitive to C16:0 PAF and five deletion mutants that were hypersensitive to C16:0 lyso-PAF. Deletion of YDL133w, a previously uncharacterized gene which we have renamed SRF1 (Spo14 Regulatory Factor 1), resulted in the greatest differential sensitivity to C16:0 PAF over C16:0 lyso-PAF. We demonstrate that Srf1 physically interacts with Spo14, yeast phospholipase D (PLD), and is essential for PLD catalytic activity in mitotic cells. Though C16:0 PAF treatment does not impact hydrolysis of phosphatidylcholine in yeast, C16:0 PAF does promote delocalization of GFP-Spo14 and phosphatidic acid from the cell periphery. Furthermore, we demonstrate that, similar to yeast cells, PLD activity is required to protect mammalian neural cells from C16:0 PAF. Together, these findings implicate PLD as a potential neuroprotective target capable of ameliorating disruptions in lipid metabolism in response to accumulating oligomeric amyloid-β42