The zebrafish grime mutant uncovers an evolutionarily conserved role for Tmem161b in the control of cardiac rhythm.

The establishment of cardiac function in the developing embryo is essential to ensure blood flow and, therefore, growth and survival of the animal. The molecular mechanisms controlling normal cardiac rhythm remain to be fully elucidated. From a forward genetic screen, we identified a unique mutant, grime, that displayed a specific cardiac arrhythmia phenotype. We show that loss-of-function mutations in tmem161b are responsible for the phenotype, identifying Tmem161b as a regulator of cardiac rhythm in zebrafish. To examine the evolutionary conservation of this function, we generated knockout mice for Tmem161b. Tmem161b knockout mice are neonatal lethal and cardiomyocytes exhibit arrhythmic calcium oscillations. Mechanistically, we find that Tmem161b is expressed at the cell membrane of excitable cells and live imaging shows it is required for action potential repolarization in the developing heart. Electrophysiology on isolated cardiomyocytes demonstrates that Tmem161b is essential to inhibit Ca2+ and K+ currents in cardiomyocytes. Importantly, Tmem161b haploinsufficiency leads to cardiac rhythm phenotypes, implicating it as a candidate gene in heritable cardiac arrhythmia. Overall, these data describe Tmem161b as a highly conserved regulator of cardiac rhythm that functions to modulate ion channel activity in zebrafish and mice.

Localised Collagen2a1 secretion supports lymphatic endothelial cell migration in the zebrafish embryo.

The lymphatic vasculature develops primarily from pre-existing veins. A pool of lymphatic endothelial cells (LECs) first sprouts from cardinal veins followed by migration and proliferation to colonise embryonic tissues. Although much is known about the molecular regulation of LEC fate and sprouting during early lymphangiogenesis, we know far less about the instructive and permissive signals that support LEC migration through the embryo. Using a forward genetic screen, we identified mbtps1 and sec23a, components of the COP-II protein secretory pathway, as essential for developmental lymphangiogenesis. In both mutants, LECs initially depart the cardinal vein but then fail in their ongoing migration. A key cargo that failed to be secreted in both mutants was a type II collagen (Col2a1). Col2a1 is normally secreted by notochord sheath cells, alongside which LECs migrate. col2a1a mutants displayed defects in the migratory behaviour of LECs and failed lymphangiogenesis. These studies thus identify Col2a1 as a key cargo secreted by notochord sheath cells and required for the migration of LECs. These findings combine with our current understanding to suggest that successive cell-to-cell and cell-matrix interactions regulate the migration of LECs through the embryonic environment during development.

mafba is a downstream transcriptional effector of Vegfc signaling essential for embryonic lymphangiogenesis in zebrafish.

The lymphatic vasculature plays roles in tissue fluid balance, immune cell trafficking, fatty acid absorption, cancer metastasis, and cardiovascular disease. Lymphatic vessels form by lymphangiogenesis, the sprouting of new lymphatics from pre-existing vessels, in both development and disease contexts. The apical signaling pathway in lymphangiogenesis is the VEGFC/VEGFR3 pathway, yet how signaling controls cellular transcriptional output remains unknown. We used a forward genetic screen in zebrafish to identify the transcription factor mafba as essential for lymphatic vessel development. We found that mafba is required for the migration of lymphatic precursors after their initial sprouting from the posterior cardinal vein. mafba expression is enriched in sprouts emerging from veins, and we show that mafba functions cell-autonomously during lymphatic vessel development. Mechanistically, Vegfc signaling increases mafba expression to control downstream transcription, and this regulatory relationship is dependent on the activity of SoxF transcription factors, which are essential for mafba expression in venous endothelium. Here we identify an indispensable Vegfc-SoxF-Mafba pathway in lymphatic development.

Kidney organoids from human iPS cells contain multiple lineages and model human nephrogenesis.

The human kidney contains up to 2 million epithelial nephrons responsible for blood filtration. Regenerating the kidney requires the induction of the more than 20 distinct cell types required for excretion and the regulation of pH, and electrolyte and fluid balance. We have previously described the simultaneous induction of progenitors for both collecting duct and nephrons via the directed differentiation of human pluripotent stem cells. Paradoxically, although both are of intermediate mesoderm in origin, collecting duct and nephrons have distinct temporospatial origins. Here we identify the developmental mechanism regulating the preferential induction of collecting duct versus kidney mesenchyme progenitors. Using this knowledge, we have generated kidney organoids that contain nephrons associated with a collecting duct network surrounded by renal interstitium and endothelial cells. Within these organoids, individual nephrons segment into distal and proximal tubules, early loops of Henle, and glomeruli containing podocytes elaborating foot processes and undergoing vascularization. When transcription profiles of kidney organoids were compared to human fetal tissues, they showed highest congruence with first trimester human kidney. Furthermore, the proximal tubules endocytose dextran and differentially apoptose in response to cisplatin, a nephrotoxicant. Such kidney organoids represent powerful models of the human organ for future applications, including nephrotoxicity screening, disease modelling and as a source of cells for therapy.

Identification of miR-29b targets using 3-cyanovinylcarbazole containing mimics.

MicroRNAs (miRNAs) are highly conserved ∼22 nt small noncoding RNAs that bind partially complementary sequences in target transcripts. MicroRNAs regulate both translation and transcript stability, and play important roles in development, cellular homeostasis, and disease. There are limited approaches available to agnostically identify microRNA targets transcriptome-wide, and methods using miRNA mimics, which in principle identify direct miRNA:transcript pairs, have low sensitivity and specificity. Here, we describe a novel method to identify microRNA targets using miR-29b mimics containing 3-cyanovinylcarbazole (CNVK), a photolabile nucleoside analog. We demonstrate that biotin-tagged, CNVK-containing miR-29b (CNVK-miR-29b) mimics are nontoxic in cell culture, associate with endogenous mammalian Argonaute2, are sensitive for known targets and recapitulate endogenous transcript destabilization. Partnering CNVK-miR-29b with ultra-low-input RNA sequencing, we recover ∼40% of known miR-29b targets and find conservation of the focal adhesion and apoptotic target pathways in mouse and human. We also identify hundreds of novel targets, including NRAS, HOXA10, and KLF11, with a validation rate of 71% for a subset of 73 novel target transcripts interrogated using a high-throughput luciferase assay. Consistent with previous reports, we show that both endogenous miR-29b and CNVK-miR-29b are trafficked to the nucleus, but find no evidence of nuclear-specific miR-29b transcript binding. This may indicate that miR-29b nuclear sequestration is a regulatory mechanism in itself. We suggest that CNVK-containing small RNA mimics may find applicability in other experimental models.

Myosin Vb is required for correct trafficking of N-cadherin and cardiac chamber ballooning.

BACKGROUND: During heart morphogenesis, the cardiac chambers undergo ballooning: a process involving regionalized elongation of cardiomyocytes. Cardiomyocyte shape changes require reorganization of the actin cytoskeleton; however, the genetic regulation of this process is not well understood. RESULTS: From a forward genetic screen, we identified the zebrafish uq 23ks mutant which manifests chamber ballooning defects. Whole-genome sequencing-mapping identified a truncating mutation in the gene, myo5b. myo5b encodes an atypical myosin required for endosome recycling and, consistent with this, increased vesicles were observed in myo5b mutant cardiomyocytes. Expression of RFP-Rab11a (a recycling endosome marker) confirmed increased recycling endosomes in cardiomyocytes of myo5b mutants. To investigate potential cargo of MyoVb-associated vesicles, we examined the adherens junction protein, N-cadherin. N-cadherin appeared mispatterned at cell junctions, and an increase in the number of intracellular particles was also apparent. Co-localization with RFP-Rab11a confirmed increased N-cadherin-positive recycling endosomes, demonstrating N-cadherin trafficking is perturbed in myo5b mutants. Finally, phalloidin staining showed disorganized F-actin in myo5b cardiomyocytes, suggesting the cytoskeleton fails to remodel, obstructing chamber ballooning. CONCLUSIONS: MyoVb is required for cardiomyocyte endosomal recycling and appropriate N-cadherin localization during the onset of chamber ballooning. Cardiomyocytes lacking MyoVb are unable to reorganize their actin cytoskeleton, resulting in failed chamber ballooning. Developmental Dynamics 248:284-295, 2019. © 2019 Wiley Periodicals, Inc.

Whole exome sequencing in patients with white matter abnormalities.

Here we report whole exome sequencing (WES) on a cohort of 71 patients with persistently unresolved white matter abnormalities with a suspected diagnosis of leukodystrophy or genetic leukoencephalopathy. WES analyses were performed on trio, or greater, family groups. Diagnostic pathogenic variants were identified in 35% (25 of 71) of patients. Potentially pathogenic variants were identified in clinically relevant genes in a further 7% (5 of 71) of cases, giving a total yield of clinical diagnoses in 42% of individuals. These findings provide evidence that WES can substantially decrease the number of unresolved white matter cases. Ann Neurol 2016;79:1031-1037.

Neuronal cell lines as model dorsal root ganglion neurons: A transcriptomic comparison.

BACKGROUND: Dorsal root ganglion neuron-derived immortal cell lines including ND7/23 and F-11 cells have been used extensively as in vitro model systems of native peripheral sensory neurons. However, while it is clear that some sensory neuron-specific receptors and ion channels are present in these cell lines, a systematic comparison of the molecular targets expressed by these cell lines with those expressed in intact peripheral neurons is lacking. RESULTS: In this study, we examined the expression of RNA transcripts in the human neuroblastoma-derived cell line, SH-SY5Y, and two dorsal root ganglion hybridoma cell lines, F-11 and ND7/23, using Illumina next-generation sequencing, and compared the results with native whole murine dorsal root ganglions. The gene expression profiles of these three cell lines did not resemble any specific defined dorsal root ganglion subclass. The cell lines lacked many markers for nociceptive sensory neurons, such as the Transient receptor potential V1 gene, but expressed markers for both myelinated and unmyelinated neurons. Global gene ontology analysis on whole dorsal root ganglions and cell lines showed similar enrichment of biological process terms across all samples. CONCLUSIONS: This paper provides insights into the receptor repertoire expressed in common dorsal root ganglion neuron-derived cell lines compared with whole murine dorsal root ganglions, and illustrates the limits and potentials of these cell lines as tools for neuropharmacological exploration.

Evolution of equine influenza virus in vaccinated horses.

Influenza A viruses are characterized by their ability to evade host immunity, even in vaccinated individuals. To determine how prior immunity shapes viral diversity in vivo, we studied the intra- and interhost evolution of equine influenza virus in vaccinated horses. Although the level and structure of genetic diversity were similar to those in naïve horses, intrahost bottlenecks may be more stringent in vaccinated animals, and mutations shared among horses often fall close to putative antigenic sites.

Transmission of equine influenza virus during an outbreak is characterized by frequent mixed infections and loose transmission bottlenecks.

The ability of influenza A viruses (IAVs) to cross species barriers and evade host immunity is a major public health concern. Studies on the phylodynamics of IAVs across different scales - from the individual to the population - are essential for devising effective measures to predict, prevent or contain influenza emergence. Understanding how IAVs spread and evolve during outbreaks is critical for the management of epidemics. Reconstructing the transmission network during a single outbreak by sampling viral genetic data in time and space can generate insights about these processes. Here, we obtained intra-host viral sequence data from horses infected with equine influenza virus (EIV) to reconstruct the spread of EIV during a large outbreak. To this end, we analyzed within-host viral populations from sequences covering 90% of the infected yards. By combining gene sequence analyses with epidemiological data, we inferred a plausible transmission network, in turn enabling the comparison of transmission patterns during the course of the outbreak and revealing important epidemiological features that were not apparent using either approach alone. The EIV populations displayed high levels of genetic diversity, and in many cases we observed distinct viral populations containing a dominant variant and a number of related minor variants that were transmitted between infectious horses. In addition, we found evidence of frequent mixed infections and loose transmission bottlenecks in these naturally occurring populations. These frequent mixed infections likely influence the size of epidemics.

Evolution of an Eurasian avian-like influenza virus in naïve and vaccinated pigs.

Influenza viruses are characterized by an ability to cross species boundaries and evade host immunity, sometimes with devastating consequences. The 2009 pandemic of H1N1 influenza A virus highlights the importance of pigs in influenza emergence, particularly as intermediate hosts by which avian viruses adapt to mammals before emerging in humans. Although segment reassortment has commonly been associated with influenza emergence, an expanded host-range is also likely to be associated with the accumulation of specific beneficial point mutations. To better understand the mechanisms that shape the genetic diversity of avian-like viruses in pigs, we studied the evolutionary dynamics of an Eurasian Avian-like swine influenza virus (EA-SIV) in naïve and vaccinated pigs linked by natural transmission. We analyzed multiple clones of the hemagglutinin 1 (HA1) gene derived from consecutive daily viral populations. Strikingly, we observed both transient and fixed changes in the consensus sequence along the transmission chain. Hence, the mutational spectrum of intra-host EA-SIV populations is highly dynamic and allele fixation can occur with extreme rapidity. In addition, mutations that could potentially alter host-range and antigenicity were transmitted between animals and mixed infections were commonplace, even in vaccinated pigs. Finally, we repeatedly detected distinct stop codons in virus samples from co-housed pigs, suggesting that they persisted within hosts and were transmitted among them. This implies that mutations that reduce viral fitness in one host, but which could lead to fitness benefits in a novel host, can circulate at low frequencies.

Restriction of V3 region sequence divergence in the HIV-1 envelope gene during antiretroviral treatment in a cohort of recent seroconverters.

BACKGROUND: Dynamic changes in Human Immunodeficiency Virus 1 (HIV-1) sequence diversity and divergence are associated with immune control during primary infection and progression to AIDS. Consensus sequencing or single genome amplification sequencing of the HIV-1 envelope (env) gene, in particular the variable (V) regions, is used as a marker for HIV-1 genome diversity, but population diversity is only minimally, or semi-quantitatively sampled using these methods. RESULTS: Here we use second generation deep sequencing to determine inter-and intra-patient sequence heterogeneity and to quantify minor variants in a cohort of individuals either receiving or not receiving antiretroviral treatment following seroconversion; the SPARTAC trial. We show, through a cross-sectional study of sequence diversity of the env V3 in 30 antiretroviral-naive patients during primary infection that considerable population structure diversity exists, with some individuals exhibiting highly constrained plasma virus diversity. Diversity was independent of clinical markers (viral load, time from seroconversion, CD4 cell count) of infection. Serial sampling over 60 weeks of non-treated individuals that define three initially different diversity profiles showed that complex patterns of continuing HIV-1 sequence diversification and divergence could be readily detected. Evidence for minor sequence turnover, emergence of new variants and re-emergence of archived variants could be inferred from this analysis. Analysis of viral divergence over the same time period in patients who received short (12 weeks, ART12) or long course antiretroviral therapy (48 weeks, ART48) and a non-treated control group revealed that ART48 successfully suppressed viral divergence while ART12 did not have a significant effect. CONCLUSIONS: Deep sequencing is a sensitive and reliable method for investigating the diversity of the env V3 as an important component of HIV-1 genome diversity. Detailed insights into the complex early intra-patient dynamics of env V3 diversity and divergence were explored in antiretroviral-naïve recent seroconverters. Long course antiretroviral therapy, initiated soon after seroconversion and administered for 48 weeks, restricts HIV-1 divergence significantly. The effect of ART12 and ART48 on clinical markers of HIV infection and progression is currently investigated in the SPARTAC trial.

Evolutionary dynamics of local pandemic H1N1/2009 influenza virus lineages revealed by whole-genome analysis.

Virus gene sequencing and phylogenetics can be used to study the epidemiological dynamics of rapidly evolving viruses. With complete genome data, it becomes possible to identify and trace individual transmission chains of viruses such as influenza virus during the course of an epidemic. Here we sequenced 153 pandemic influenza H1N1/09 virus genomes from United Kingdom isolates from the first (127 isolates) and second (26 isolates) waves of the 2009 pandemic and used their sequences, dates of isolation, and geographical locations to infer the genetic epidemiology of the epidemic in the United Kingdom. We demonstrate that the epidemic in the United Kingdom was composed of many cocirculating lineages, among which at least 13 were exclusively or predominantly United Kingdom clusters. The estimated divergence times of two of the clusters predate the detection of pandemic H1N1/09 virus in the United Kingdom, suggesting that the pandemic H1N1/09 virus was already circulating in the United Kingdom before the first clinical case. Crucially, three clusters contain isolates from the second wave of infections in the United Kingdom, two of which represent chains of transmission that appear to have persisted within the United Kingdom between the first and second waves. This demonstrates that whole-genome analysis can track in fine detail the behavior of individual influenza virus lineages during the course of a single epidemic or pandemic.

ahctf1 and kras mutations combine to amplify oncogenic stress and restrict liver overgrowth in a zebrafish model of hepatocellular carcinoma.

The nucleoporin (NUP) ELYS, encoded by AHCTF1, is a large multifunctional protein with essential roles in nuclear pore assembly and mitosis. Using both larval and adult zebrafish models of hepatocellular carcinoma (HCC), in which the expression of an inducible mutant kras transgene (krasG12V) drives hepatocyte-specific hyperplasia and liver enlargement, we show that reducing ahctf1 gene dosage by 50% markedly decreases liver volume, while non-hyperplastic tissues are unaffected. We demonstrate that in the context of cancer, ahctf1 heterozygosity impairs nuclear pore formation, mitotic spindle assembly, and chromosome segregation, leading to DNA damage and activation of a Tp53-dependent transcriptional programme that induces cell death and cell cycle arrest. Heterozygous expression of both ahctf1 and ranbp2 (encoding a second nucleoporin), or treatment of heterozygous ahctf1 larvae with the nucleocytoplasmic transport inhibitor, Selinexor, completely blocks krasG12V-driven hepatocyte hyperplasia. Gene expression analysis of patient samples in the liver hepatocellular carcinoma (LIHC) dataset in The Cancer Genome Atlas shows that high expression of one or more of the transcripts encoding the 10 components of the NUP107-160 subcomplex, which includes AHCTF1, is positively correlated with worse overall survival. These results provide a strong and feasible rationale for the development of novel cancer therapeutics that target ELYS function and suggest potential avenues for effective combinatorial treatments.

Intrahost evolutionary dynamics of canine influenza virus in naive and partially immune dogs.

The patterns and dynamics of evolution in acutely infecting viruses within individual hosts are largely unknown. To this end, we investigated the intrahost variation of canine influenza virus (CIV) during the course of experimental infections in naïve and partially immune dogs and in naturally infected dogs. Tracing sequence diversity in the gene encoding domain 1 of the hemagglutinin (HA1) protein over the time course of infection provided information on the patterns and processes of intrahost viral evolution and revealed some of the effects of partial host immunity. Viral populations sampled on any given day were generally characterized by mean pairwise genetic diversities between 0.1 and 0.2% and by mutational spectra that changed considerably on different days. Some observed mutations may have affected antigenicity or host range, including reversions of CIV host-associated mutations. Patterns of sequence diversity differed between naïve and vaccinated dogs, with some presumably antigenic mutations transiently reaching high frequency in the latter. CIV populations are therefore characterized by the rapid generation and clearance of genetic diversity. Potentially advantageous mutations arise readily during the course of single infections and may give rise to antigenic escape or host range variants.

Intra- and interhost evolutionary dynamics of equine influenza virus.

Determining the evolutionary basis of cross-species transmission and immune evasion is key to understanding the mechanisms that control the emergence of either new viruses or novel antigenic variants with pandemic potential. The hemagglutinin glycoprotein of influenza A viruses is a critical host range determinant and a major target of neutralizing antibodies. Equine influenza virus (EIV) is a significant pathogen of the horse that causes periodical outbreaks of disease even in populations with high vaccination coverage. EIV has also jumped the species barrier and emerged as a novel respiratory pathogen in dogs, canine influenza virus. We studied the dynamics of equine influenza virus evolution in horses at the intrahost level and how this evolutionary process is affected by interhost transmission in a natural setting. To this end, we performed clonal sequencing of the hemagglutinin 1 gene derived from individual animals at different times postinfection. Our results show that despite the population consensus sequence remaining invariant, genetically distinct subpopulations persist during the course of infection and are also transmitted, with some variants likely to change antigenicity. We also detected a natural case of mixed infection in an animal infected during an outbreak of equine influenza, raising the possibility of reassortment between different strains of virus. In sum, our data suggest that transmission bottlenecks may not be as narrow as originally perceived and that the genetic diversity required to adapt to new host species may be partially present in the donor host and potentially transmitted to the recipient host.

The RNA helicase Ddx21 controls Vegfc-driven developmental lymphangiogenesis by balancing endothelial cell ribosome biogenesis and p53 function.

The development of a functional vasculature requires the coordinated control of cell fate, lineage differentiation and network growth. Cellular proliferation is spatiotemporally regulated in developing vessels, but how this is orchestrated in different lineages is unknown. Here, using a zebrafish genetic screen for lymphatic-deficient mutants, we uncover a mutant for the RNA helicase Ddx21. Ddx21 cell-autonomously regulates lymphatic vessel development. An established regulator of ribosomal RNA synthesis and ribosome biogenesis, Ddx21 is enriched in sprouting venous endothelial cells in response to Vegfc-Flt4 signalling. Ddx21 function is essential for Vegfc-Flt4-driven endothelial cell proliferation. In the absence of Ddx21, endothelial cells show reduced ribosome biogenesis, p53 and p21 upregulation and cell cycle arrest that blocks lymphangiogenesis. Thus, Ddx21 coordinates the lymphatic endothelial cell response to Vegfc-Flt4 signalling by balancing ribosome biogenesis and p53 function. This mechanism may be targetable in diseases of excessive lymphangiogenesis such as cancer metastasis or lymphatic malformation.

ORAI1 and ORAI3 in Breast Cancer Molecular Subtypes and the Identification of ORAI3 as a Hypoxia Sensitive Gene and a Regulator of Hypoxia Responses.

The remodeling of specific calcium-permeable ion channels is a feature of some breast cancer subtypes. ORAI1 is a protein that forms a calcium-permeable ion channel responsible for store-operated calcium entry (SOCE) in a variety of cell types. ORAI3, a related isoform, is not a regulator of SOCE in most cell types. However, ORAI3 does control SOCE in many estrogen receptor-positive breast cancer cell lines, where it also controls proliferation. ORAI1 is a well-characterized regulator of the proliferation and migration of many basal breast cancer cells; however, the role of ORAI3 in these types of breast cancer cells remains unclear. Here, we sought to define ORAI1 and ORAI3 expression in breast cancer cell lines of different molecular subtypes and assess the potential role and regulation of ORAI3 in basal breast cancer cells. Our study demonstrates that elevated ORAI1 is a feature of basal-like breast cancers, while elevated ORAI3 is a feature of luminal breast cancers. Intriguingly, we found that ORAI3 is over-expressed in the mesenchymal subtype of triple-negative breast cancer. Given this, we assessed ORAI3 levels in the presence of two inducers of the mesenchymal phenotype, hypoxia and epidermal growth factor (EGF). Hypoxia induced ORAI3 levels in basal breast cancer cell lines through a pathway involving hypoxia-inducible factor-1 alpha (HIF1α. The silencing of ORAI3 attenuated hypoxia-associated phosphorylation of the EGF receptor (EGFR) and the expression of genes associated with cell migration and inflammatory/immune responses in the MDA-MB-468 model of basal breast cancer. Although elevated ORAI3 levels were not associated with survival; basal, estrogen receptor-negative and triple-negative breast cancers with high ORAI3 and low ORAI1 levels were associated with poorer clinical outcomes. This study defines ORAI3 as a potential fine-tuner for processes relevant to the progression of basal breast cancers.

The Alternative Splicing Regulator Nova2 Constrains Vascular Erk Signaling to Limit Specification of the Lymphatic Lineage.

The correct assignment of cell fate within fields of multipotent progenitors is essential for accurate tissue diversification. The first lymphatic vessels arise from pre-existing veins after venous endothelial cells become specified as lymphatic progenitors. Prox1 specifies lymphatic fate and labels these progenitors; however, the mechanisms restricting Prox1 expression and limiting the progenitor pool remain unknown. We identified a zebrafish mutant that displayed premature, expanded, and prolonged lymphatic specification. The gene responsible encodes the regulator of alternative splicing, Nova2. In zebrafish and human endothelial cells, Nova2 selectively regulates pre-mRNA splicing for components of signaling pathways and phosphoproteins. Nova2-deficient endothelial cells display increased Mapk/Erk signaling, and Prox1 expression is dynamically controlled by Erk signaling. We identify a mechanism whereby Nova2-regulated splicing constrains Erk signaling, thus limiting lymphatic progenitor cell specification. This identifies the capacity of a factor that tunes mRNA splicing to control assignment of cell fate during vascular differentiation.