Reduced LYNX1 expression in transcriptome of human iPSC-derived neural progenitors modeling fragile X syndrome.

Lack of FMR1 protein results in fragile X syndrome (FXS), which is the most common inherited intellectual disability syndrome and serves as an excellent model disease to study molecular mechanisms resulting in neuropsychiatric comorbidities. We compared the transcriptomes of human neural progenitors (NPCs) generated from patient-derived induced pluripotent stem cells (iPSCs) of three FXS and three control male donors. Altered expression of RAD51C, PPIL3, GUCY1A2, MYD88, TRAPPC4, LYNX1, and GTF2A1L in FXS NPCs suggested changes related to triplet repeat instability, RNA splicing, testes development, and pathways previously shown to be affected in FXS. LYNX1 is a cholinergic brake of tissue plasminogen activator (tPA)-dependent plasticity, and its reduced expression was consistent with augmented tPA-dependent radial glial process growth in NPCs derived from FXS iPSC lines. There was evidence of human iPSC line donor-dependent variation reflecting potentially phenotypic variation. NPCs derived from an FXS male with concomitant epilepsy expressed differently several epilepsy-related genes, including genes shown to cause the auditory epilepsy phenotype in the murine model of FXS. Functional enrichment analysis highlighted regulation of insulin-like growth factor pathway in NPCs modeling FXS with epilepsy. Our results demonstrated potential of human iPSCs in disease modeling for discovery and development of therapeutic interventions by showing early gene expression changes in FXS iPSC-derived NPCs consistent with the known pathophysiological changes in FXS and by revealing disturbed FXS progenitor growth linked to reduced expression of LYNX1, suggesting dysregulated cholinergic system.

MircoRNA in Extracellular Vesicles from Patients with Pulmonary Arterial Hypertension Alters Endothelial Angiogenic Response.

Pulmonary arterial hypertension (PAH) is characterized by a progressive elevation of pulmonary pressure leading to right ventricular dysfunction and is associated with a poor prognosis. Patients with PAH have increased numbers of circulating extracellular vesicles (EVs) and altered expression of circulating microRNAs (miRs). The study aimed to evaluate the miR profile contained within purified EVs derived from the plasma of PAH patients as compared to healthy controls (HC). Circulating EVs, purified from platelet-free plasma were analyzed using flow cytometry, western blot, and electron microscopy. Total RNA isolated from EVs was subjected to Microarray analysis using GeneChip miRNA 4.0 Array and bioinformatics tools. Overexpression and inhibition of miRs were conducted in human pulmonary artery endothelial cells (hPAECs) that had been incubated previously with either PAH- or HC-derived EVs. Cell proliferation (MTT assay) and angiogenesis (tube formation assay) were tested in hPAECs to determine miR functionality. MiR profiling revealed 370 heats while comparing PAH and HC groups, 22 of which were found to be down-regulated and 6 were up-regulated in the PAH EVs. Among the altered miRs, miR-486-5p was overexpressed, while miR-26a-5p was downregulated in PAH EVs compared to HC EVs. Inhibition of mir-486-5p or overexpression of miR-26a-5p in hPAECs post-exposure of PAH EVs abrogated proangiogenic and proliferative effects posed by PAH EVs contrary to HC EVs. The angiogenic and proliferative effects of the miRs from PAH EVs were observed to be mediated through nuclear factor (NF)-κB activation. PAH EVs carry and present an altered miR profile that can be targeted to restrict angiogenesis and reduce pulmonary endothelium activation. Further studies concerning miRs from circulating heterogeneous EVs in PAH patients are warranted to understand their potential as targets for treatment in PAH.

Urokinase plasminogen activator mediates changes in human astrocytes modeling fragile X syndrome.

The function of astrocytes intertwines with the extracellular matrix, whose neuron and glial cell-derived components shape neuronal plasticity. Astrocyte abnormalities have been reported in the brain of the mouse model for fragile X syndrome (FXS), the most common cause of inherited intellectual disability, and a monogenic cause of autism spectrum disorder. We compared human FXS and control astrocytes generated from human induced pluripotent stem cells and we found increased expression of urokinase plasminogen activator (uPA), which modulates degradation of extracellular matrix. Several pathways associated with uPA and its receptor function were activated in FXS astrocytes. Levels of uPA were also increased in conditioned medium collected from FXS hiPSC-derived astrocyte cultures and correlated inversely with intracellular Ca2+ responses to activation of L-type voltage-gated calcium channels in human astrocytes. Increased uPA augmented neuronal phosphorylation of TrkB within the docking site for the phospholipase-Cγ1 (PLCγ1), indicating effects of uPA on neuronal plasticity. Gene expression changes during neuronal differentiation preceding astrogenesis likely contributed to properties of astrocytes with FXS-specific alterations that showed specificity by not affecting differentiation of adenosine triphosphate (ATP)-responsive astrocyte population. To conclude, our studies identified uPA as an important regulator of astrocyte function and demonstrated that increased uPA in human FXS astrocytes modulated astrocytic responses and neuronal plasticity.

Integrated Molecular Analysis of Undifferentiated Uterine Sarcomas Reveals Clinically Relevant Molecular Subtypes.

PURPOSE: Undifferentiated uterine sarcomas (UUS) are rare, extremely deadly, sarcomas with no effective treatment. The goal of this study was to identify novel intrinsic molecular UUS subtypes using integrated clinical, histopathologic, and molecular evaluation of a large, fully annotated, patient cohort. EXPERIMENTAL DESIGN: Fifty cases of UUS with full clinicopathologic annotation were analyzed for gene expression (n = 50), copy-number variation (CNV, n = 40), cell morphometry (n = 39), and protein expression (n = 22). Gene ontology and network enrichment analysis were used to relate over- and underexpressed genes to pathways and further to clinicopathologic and phenotypic findings. RESULTS: Gene expression identified four distinct groups of tumors, which varied in their clinicopathologic parameters. Gene ontology analysis revealed differential activation of pathways related to genital tract development, extracellular matrix (ECM), muscle function, and proliferation. A multivariable, adjusted Cox proportional hazard model demonstrated that RNA group, mitotic index, and hormone receptor expression influence patient overall survival (OS). CNV arrays revealed characteristic chromosomal changes for each group. Morphometry demonstrated that the ECM group, the most aggressive, exhibited a decreased cell density and increased nuclear area. A cell density cutoff of 4,300 tumor cells per mm2 could separate ECM tumors from the remaining cases with a sensitivity of 83% and a specificity of 94%. IHC staining of MMP-14, Collagens 1 and 6, and Fibronectin proteins revealed differential expression of these ECM-related proteins, identifying potential new biomarkers for this aggressive sarcoma subgroup. CONCLUSIONS: Molecular evaluation of UUS provides novel insights into the biology, prognosis, phenotype, and possible treatment of these tumors.

Increased Calcium Influx through L-type Calcium Channels in Human and Mouse Neural Progenitors Lacking Fragile X Mental Retardation Protein.

The absence of FMR1 protein (FMRP) causes fragile X syndrome (FXS) and disturbed FMRP function is implicated in several forms of human psychopathology. We show that intracellular calcium responses to depolarization are augmented in neural progenitors derived from human induced pluripotent stem cells and mouse brain with FXS. Increased calcium influx via nifedipine-sensitive voltage-gated calcium (Cav) channels contributes to the exaggerated responses to depolarization and type 1 metabotropic glutamate receptor activation. The ratio of L-type/T-type Cav channel expression is increased in FXS progenitors and correlates with enhanced progenitor differentiation to glutamate-responsive cells. Genetic reduction of brain-derived neurotrophic factor in FXS mouse progenitors diminishes the expression of Cav channels and activity-dependent responses, which are associated with increased phosphorylation of the phospholipase C-γ1 site within TrkB receptors and changes of differentiating progenitor subpopulations. Our results show developmental effects of increased calcium influx via L-type Cav channels in FXS neural progenitors.

Large-scale suppression of recombination predates genomic rearrangements in Neurospora tetrasperma.

A common feature of eukaryote genomes is large chromosomal regions where recombination is absent or strongly reduced, but the factors that cause this reduction are not well understood. Genomic rearrangements have often been implicated, but they may also be a consequence of recombination suppression rather than a cause. In this study, we generate eight high-quality genomic data sets of the filamentous ascomycete Neurospora tetrasperma, a fungus that lacks recombination over most of its largest chromosome. The genomes surprisingly reveal collinearity of the non-recombining regions and although large inversions are enriched in these regions, we conclude these inversions to be derived and not the cause of the suppression. To our knowledge, this is the first time that non-recombining, genic regions as large as 86% of a full chromosome (or 8 Mbp), are shown to be collinear. These findings are of significant interest for our understanding of the evolution of sex chromosomes and other supergene complexes.

Determinants of the Efficacy of Natural Selection on Coding and Noncoding Variability in Two Passerine Species.

Population genetic theory predicts that selection should be more effective when the effective population size (Ne) is larger, and that the efficacy of selection should correlate positively with recombination rate. Here, we analyzed the genomes of ten great tits and ten zebra finches. Nucleotide diversity at 4-fold degenerate sites indicates that zebra finches have a 2.83-fold larger Ne. We obtained clear evidence that purifying selection is more effective in zebra finches. The proportion of substitutions at 0-fold degenerate sites fixed by positive selection (α) is high in both species (great tit 48%; zebra finch 64%) and is significantly higher in zebra finches. When α was estimated on GC-conservative changes (i.e., between A and T and between G and C), the estimates reduced in both species (great tit 22%; zebra finch 53%). A theoretical model presented herein suggests that failing to control for the effects of GC-biased gene conversion (gBGC) is potentially a contributor to the overestimation of α, and that this effect cannot be alleviated by first fitting a demographic model to neutral variants. We present the first estimates in birds for α in the untranslated regions, and found evidence for substantial adaptive changes. Finally, although purifying selection is stronger in high-recombination regions, we obtained mixed evidence for α increasing with recombination rate, especially after accounting for gBGC. These results highlight that it is important to consider the potential confounding effects of gBGC when quantifying selection and that our understanding of what determines the efficacy of selection is incomplete.

Introgression maintains the genetic integrity of the mating-type determining chromosome of the fungus Neurospora tetrasperma.

Genome evolution is driven by a complex interplay of factors, including selection, recombination, and introgression. The regions determining sexual identity are particularly dynamic parts of eukaryotic genomes that are prone to molecular degeneration associated with suppressed recombination. In the fungus Neurospora tetrasperma, it has been proposed that this molecular degeneration is counteracted by the introgression of nondegenerated DNA from closely related species. In this study, we used comparative and population genomic analyses of 92 genomes from eight phylogenetically and reproductively isolated lineages of N. tetrasperma, and its three closest relatives, to investigate the factors shaping the evolutionary history of the genomes.We found that suppressed recombination extends across at least 6 Mbp (∼ 63%) of the mating-type (mat) chromosome in N. tetrasperma and is associated with decreased genetic diversity, which is likely the result primarily of selection at linked sites. Furthermore, analyses of molecular evolution revealed an increased mutational load in this region, relative to recombining regions. However, comparative genomic and phylogenetic analyses indicate that the mat chromosomes are temporarily regenerated via introgression from sister species; six of eight lineages show introgression into one of their mat chromosomes, with multiple Neurospora species acting as donors. The introgressed tracts have been fixed within lineages, suggesting that they confer an adaptive advantage in natural populations, and our analyses support the presence of selective sweeps in at least one lineage. Thus, these data strongly support the previously hypothesized role of introgression as a mechanism for the maintenance of mating-type determining chromosomal regions.

Evolutionary signals of selection on cognition from the great tit genome and methylome.

For over 50 years, the great tit (Parus major) has been a model species for research in evolutionary, ecological and behavioural research; in particular, learning and cognition have been intensively studied. Here, to provide further insight into the molecular mechanisms behind these important traits, we de novo assemble a great tit reference genome and whole-genome re-sequence another 29 individuals from across Europe. We show an overrepresentation of genes related to neuronal functions, learning and cognition in regions under positive selection, as well as increased CpG methylation in these regions. In addition, great tit neuronal non-CpG methylation patterns are very similar to those observed in mammals, suggesting a universal role in neuronal epigenetic regulation which can affect learning-, memory- and experience-induced plasticity. The high-quality great tit genome assembly will play an instrumental role in furthering the integration of ecological, evolutionary, behavioural and genomic approaches in this model species.

The Effects of Background and Interference Selection on Patterns of Genetic Variation in Subdivided Populations.

It is well known that most new mutations that affect fitness exert deleterious effects and that natural populations are often composed of subpopulations (demes) connected by gene flow. To gain a better understanding of the joint effects of purifying selection and population structure, we focus on a scenario where an ancestral population splits into multiple demes and study neutral diversity patterns in regions linked to selected sites. In the background selection regime of strong selection, we first derive analytic equations for pairwise coalescent times and FST as a function of time after the ancestral population splits into two demes and then construct a flexible coalescent simulator that can generate samples under complex models such as those involving multiple demes or nonconservative migration. We have carried out extensive forward simulations to show that the new methods can accurately predict diversity patterns both in the nonequilibrium phase following the split of the ancestral population and in the equilibrium between mutation, migration, drift, and selection. In the interference selection regime of many tightly linked selected sites, forward simulations provide evidence that neutral diversity patterns obtained from both the nonequilibrium and equilibrium phases may be virtually indistinguishable for models that have identical variance in fitness, but are nonetheless different with respect to the number of selected sites and the strength of purifying selection. This equivalence in neutral diversity patterns suggests that data collected from subdivided populations may have limited power for differentiating among the selective pressures to which closely linked selected sites are subject.

Plant-mediated gene silencing restricts growth of the potato late blight pathogen Phytophthora infestans.

Phytophthora infestans is an oomycete that causes severe damage to potato, and is well known for its ability to evolve rapidly in order to overcome resistant potato varieties. An RNA silencing strategy was evaluated here to clarify if small interfering RNA homologous to selected genes in P. infestans could be targeted from the plant host to reduce the magnitude of the infection. As a proof-of-concept, a hairpin RNA (hp-RNA) construct using the GFP marker gene was designed and introduced in potato. At 72 hpi, a 55-fold reduction of the signal intensity of a corresponding GFP expressing P. infestans strain on leaf samples of transgenic plants, compared with wild-type potato, was detected. This suggests that an RNA interference construct in the potato host could be processed and target a transcript of the pathogen. Three genes important in the infection process of P. infestans, PiGPB1, PiCESA2, and PiPEC, together with PiGAPDH taking part in basic cell maintenance were subsequently tested using an analogous transgenic strategy. Out of these gene candidates, the hp-PiGPB1 targeting the G protein ß-subunit (PiGPB1) important for pathogenicity resulted in most restricted disease progress. Further, Illumina sequencing of inoculated transgenic potato leaves revealed sRNAs of 24/25 nt size homologous to the PiGPB1 gene in the transgenic plants indicating post-transcriptional silencing of the target gene. The work demonstrates that a host-induced gene-silencing approach is functional against P. infestans but is highly dependent on target gene for a successful outcome. This finding broadens the arsenal of control strategies to this important plant disease.

Intron evolution in Neurospora: the role of mutational bias and selection.

We used comparative and population genomics to study intron evolutionary dynamics in the fungal model genus Neurospora. For our investigation, we used well-annotated genomes of N. crassa, N. discreta, and N. tetrasperma, and 92 resequenced genomes of N. tetrasperma from natural populations. By analyzing the four well-annotated genomes, we identified 9495 intron sites in 7619 orthologous genes. Our data supports nonhomologous end joining (NHEJ) and tandem duplication as mechanisms for intron gains in the genus and the RT-mRNA process as a mechanism for intron loss. We found a moderate intron gain rate (5.78-6.89 × 10(-13) intron gains per nucleotide site per year) and a high intron loss rate (7.53-13.76 × 10(-10) intron losses per intron sites per year) as compared to other eukaryotes. The derived intron gains and losses are skewed to high frequencies, relative to neutral SNPs, in natural populations of N. tetrasperma, suggesting that selection is involved in maintaining a high intron turnover. Furthermore, our analyses of the association between intron population-level frequency and genomic features suggest that selection is involved in shaping a 5' intron position bias and a low intron GC content. However, intron sequence analyses suggest that the gained introns were not exposed to recent selective sweeps. Taken together, this work contributes to our understanding of the importance of mutational bias and selection in shaping the intron distribution in eukaryotic genomes.

Fragmentation of tRNA in Phytophthora infestans asexual life cycle stages and during host plant infection.

BACKGROUND: The oomycete Phytophthora infestans possesses active RNA silencing pathways, which presumably enable this plant pathogen to control the large numbers of transposable elements present in its 240 Mb genome. Small RNAs (sRNAs), central molecules in RNA silencing, are known to also play key roles in this organism, notably in regulation of critical effector genes needed for infection of its potato host. RESULTS: To identify additional classes of sRNAs in oomycetes, we mapped deep sequencing reads to transfer RNAs (tRNAs) thereby revealing the presence of 19-40 nt tRNA-derived RNA fragments (tRFs). Northern blot analysis identified abundant tRFs corresponding to half tRNA molecules. Some tRFs accumulated differentially during infection, as seen by examining sRNAs sequenced from P. infestans-potato interaction libraries. The putative connection between tRF biogenesis and the canonical RNA silencing pathways was investigated by employing hairpin RNA-mediated RNAi to silence the genes encoding P. infestans Argonaute (PiAgo) and Dicer (PiDcl) endoribonucleases. By sRNA sequencing we show that tRF accumulation is PiDcl1-independent, while Northern hybridizations detected reduced levels of specific tRNA-derived species in the PiAgo1 knockdown line. CONCLUSIONS: Our findings extend the sRNA diversity in oomycetes to include fragments derived from non-protein-coding RNA transcripts and identify tRFs with elevated levels during infection of potato by P. infestans.

A global multilocus analysis of the model fungus Neurospora reveals a single recent origin of a novel genetic system.

The large diversity of mating systems observed in the fungal kingdom underlines the importance of mating system change in fungal evolution. The selfing species Neurospora tetrasperma has evolved a novel method of achieving self-fertility by a mating system referred to as pseudohomothallism. However, little is known about the origin of N. tetrasperma and its relationship to the self-sterile, heterothallic, Neurospora species. In this study, we used a combination of phylogenetic and population genetic analyses to reconstruct the evolutionary history of N. tetrasperma and its heterothallic relatives. We sequenced 9 unlinked nuclear loci from 106 strains of N. tetrasperma sampled from across the globe, and a sample of 28 heterothallic strains of Neurospora. Our analyses provide strong support for monophyly of N. tetrasperma, but reject the monophyly of N. crassa. We estimate that N. tetrasperma is of a recent origin and that it diverged from the heterothallic species ∼1 million years ago. We also extend previous findings on the diversification within the N. tetrasperma clade, with 10 lineages identified. Taken together, these findings indicate that N. tetrasperma is younger than has been previously reported and that a rapid diversification of lineages has occurred within the N. tetrasperma clade.

Quantifying functional heterothallism in the pseudohomothallic ascomycete Neurospora tetrasperma.

Neurospora tetrasperma is a pseudohomothallic filamentous ascomycete that has evolved from heterothallic ancestors. Throughout its life cycle, it is predominantly heterokaryotic for mating type, and thereby self-fertile. However, studies of N. tetrasperma have revealed the occasional production of self-sterile asexual and sexual spores of a single-mating type, indicating that it can be functionally heterothallic. Here, we report the extensive sampling and isolation of natural, heterokaryotic, strains of N. tetrasperma from the United Kingdom (UK): 99 strains were collected from Surrey, England, and four from Edinburgh, Scotland. We verified by phylogenetic analyses that these strains belong to N. tetrasperma. We isolated cultures from single germinated asexual spores (conidia) from 17 of these newly sampled UK strains from Surrey, and 16 previously sampled strains of N. tetrasperma from New Zealand (NZ). Our results show that the N. tetrasperma strains from the UK population produced a significantly greater proportion of self-sterile, homokaryotic conidia than the NZ population: the proportion of homokaryotic conidia was 42.6 % (133/312 spores) and 15.3 % (59/386) from the UK and the NZ populations, respectively. Although homokaryons recovered from several strains show a bias for one of the mating types, the total ratio of mat A to mat a mating type in homokaryons (UK: 72/61, NZ 28/31) did not deviate significantly from the expected 1:1 ratio for either of these populations. These results indicate that different populations exhibit differences in their life cycle characteristics, and that a higher degree of outcrossing might be expected from the UK population. This study points to the importance of studying multiple strains and populations when investigating life history traits of an organism with a complex life cycle, as previously undetected differences between populations may be revealed.

Large-scale introgression shapes the evolution of the mating-type chromosomes of the filamentous ascomycete Neurospora tetrasperma.

The significance of introgression as an evolutionary force shaping natural populations is well established, especially in animal and plant systems. However, the abundance and size of introgression tracts, and to what degree interspecific gene flow is the result of adaptive processes, are largely unknown. In this study, we present medium coverage genomic data from species of the filamentous ascomycete Neurospora, and we use comparative genomics to investigate the introgression landscape at the genomic level in this model genus. We revealed one large introgression tract in each of the three investigated phylogenetic lineages of Neurospora tetrasperma (sizes of 5.6 Mbp, 5.2 Mbp, and 4.1 Mbp, respectively). The tract is located on the chromosome containing the locus conferring sexual identity, the mating-type (mat) chromosome. The region of introgression is confined to the region of suppressed recombination and is found on one of the two mat chromosomes (mat a). We used Bayesian concordance analyses to exclude incomplete lineage sorting as the cause for the observed pattern, and multilocus genealogies from additional species of Neurospora show that the introgression likely originates from two closely related, freely recombining, heterothallic species (N. hispaniola and N. crassa/N. perkinsii). Finally, we investigated patterns of molecular evolution of the mat chromosome in Neurospora, and we show that introgression is correlated with reduced level of molecular degeneration, consistent with a shorter time of recombination suppression. The chromosome specific (mat) and allele specific (mat a) introgression reported herein comprise the largest introgression tracts reported to date from natural populations. Furthermore, our data contradicts theoretical predictions that introgression should be less likely on sex-determining chromosomes. Taken together, the data presented herein advance our general understanding of introgression as a force shaping eukaryotic genomes.

Coalescent-based analysis distinguishes between allo- and autopolyploid origin in Shepherd's Purse (Capsella bursa-pastoris).

Polyploidization plays an important role in plant speciation. The most recent estimates report that up to 15% of angiosperm speciation events and 31% in ferns are accompanied by changes in ploidy level. Polyploids can arise either through autopolyploidy, when the sets of chromosomes originate from a single species, or through allopolyploidy, when they originate from different species. In this study, we used two different coalescent-based methods to determine the date and mode of the polyploidization event that led to the tetraploid cosmopolitan weed, Capsella bursa-pastoris. We sampled 78 C. bursa-pastoris accessions, and 53 and 43 accessions from the only two other members of this genus, C. grandiflora and C. rubella, respectively, and sequenced these accessions at 14 unlinked nuclear loci with locus-specific primers in order to be able to distinguish the two homeologues in the tetraploid. A large fraction of fixed differences between homeologous genes in C. bursa-pastoris are segregating as polymorphisms in C. grandiflora, consistent with an autopolyploid origin followed by disomic inheritance. To test this, we first estimated the demographic parameters of an isolation-with-migration model in a pairwise fashion between C. grandiflora and both genomes of C. bursa-pastoris and used these parameters in coalescent simulations to test the mode of origin of C. bursa-pastoris. Second, we used Approximate Bayesian Computation to compare an allopolyploid and an autopolyploid model. Both analyses led to the conclusion that C. bursa-pastoris originated less than 1 Ma by doubling of the C. grandiflora genome.