Stefan A. Haas

Bats are the only mammals capable of self-powered flight, an evolutionary innovation based on the transformation of forelimbs into wings. The bat wing is characterized by an extreme elongation of the second to fifth digits with a wing membrane called the chiropatagium connecting them. Here we investigated the developmental and cellular origin of this structure by comparing bat and mouse limbs using omics tools and single-cell analyses. Despite the substantial morphological differences between the species, we observed an overall conservation of cell populations and gene expression patterns including interdigital apoptosis. Single-cell analyses of micro-dissected embryonic chiropatagium identified a specific fibroblast population, independent of apoptosis-associated interdigital cells, as the origin of this tissue. These distal cells express a conserved gene programme including the transcription factors MEIS2 and TBX3, which are commonly known to specify and pattern the early proximal limb. Transgenic ectopic expression of MEIS2 and TBX3 in mouse distal limb cells resulted in the activation of genes expressed during wing development and phenotypic changes related to wing morphology, such as the fusion of digits. Our results elucidate fundamental molecular mechanisms of bat wing development and illustrate how drastic morphological changes can be achieved through repurposing of existing developmental programmes during evolution.

Intersexuality in female moles: Female moles are intersexual and develop masculinizing ovotestes, a distinctive trait among mammals. Real et al. investigated the origin of this trait by sequencing the Iberian mole genome and applying comparative strategies that integrate transcriptomic, epigenetic, and chromatin interaction data. They identified mole-specific genomic rearrangements that alter the three-dimensional regulatory landscape of the androgen-converting gene CYP17A1 and the pro-testicular factor gene FGF9 , both of which show distinct expression patterns in mole gonads. The use of transgenic mice confirms the capability of these factors to increase circulating testosterone levels and to induce gonadal masculinization. This study highlights how integrative approaches can reveal the phenotypic impact of genomic variation. Science , this issue p. 208

Virus replication displays a large cell-to-cell heterogeneity; yet, not all sources of this variability are known. Here, we study the effect of defective interfering (DI) particle (DIP) co-infection on cell-to-cell variability in influenza A virus (IAV) replication. DIPs contain a large internal deletion in one of their eight viral RNAs (vRNA) and are, thus, defective in virus replication. Moreover, they interfere with virus replication. Using single-cell isolation and reverse transcription polymerase chain reaction, we uncovered a large between-cell heterogeneity in the DI vRNA content of infected cells, which was confirmed for DI mRNAs by single-cell RNA sequencing. A high load of intracellular DI vRNAs and DI mRNAs was found in low-productive cells, indicating their contribution to the large cell-to-cell variability in virus release. Furthermore, we show that the magnitude of host cell mRNA expression (some factors may inhibit virus replication), but not the ribosome content, may further affect the strength of single-cell virus replication. Finally, we show that the load of viral mRNAs (facilitating viral protein production) and the DI mRNA content are, independently from one another, connected with single-cell virus production. Together, these insights advance single-cell virology research toward the elucidation of the complex multi-parametric origin of the large cell-to-cell heterogeneity in virus infections.

Long-range gene regulation involves physical proximity between enhancers and promoters to generate precise patterns of gene expression in space and time. However, in some cases proximity coincides with gene activation, whereas in others preformed topologies already exist before activation. In this study, we investigate the preformed configuration underlying the regulation of theShhgene by its unique limb enhancer, theZRS, in vivoduring mouse development. Abrogating the constitutive transcription covering theZRSregion led to a shift within theShh-ZRScontacts and a moderate reduction inShhtranscription. Deletion of the CTCF binding sites around theZRSresulted in a loss of theShh-ZRSpreformed interaction and a 50% decrease inShhexpression but no phenotype, suggesting an additional, CTCF-independent mechanism of promoter-enhancer communication. This residual activity, however, was diminished by combining the loss of CTCF binding with a hypomorphic ZRS allele resulting in severeShhloss-of-function and digit agenesis. Our results indicate that the preformed chromatin structure of theShhlocus is sustained by multiple components and acts to reinforce enhancer-promoter communication for robust transcription.

RLIM, also known as RNF12, is an X-linked E3 ubiquitin ligase acting as a negative regulator of LIM-domain containing transcription factors and participates in X-chromosome inactivation (XCI) in mice. We report the genetic and clinical findings of 84 individuals from nine unrelated families, eight of whom who have pathogenic variants in RLIM (RING finger LIM domain-interacting protein). A total of 40 affected males have X-linked intellectual disability (XLID) and variable behavioral anomalies with or without congenital malformations. In contrast, 44 heterozygous female carriers have normal cognition and behavior, but eight showed mild physical features. All RLIM variants identified are missense changes co-segregating with the phenotype and predicted to affect protein function. Eight of the nine altered amino acids are conserved and lie either within a domain essential for binding interacting proteins or in the C-terminal RING finger catalytic domain. In vitro experiments revealed that these amino acid changes in the RLIM RING finger impaired RLIM ubiquitin ligase activity. In vivo experiments in rlim mutant zebrafish showed that wild type RLIM rescued the zebrafish rlim phenotype, whereas the patient-specific missense RLIM variants failed to rescue the phenotype and thus represent likely severe loss-of-function mutations. In summary, we identified a spectrum of RLIM missense variants causing syndromic XLID and affecting the ubiquitin ligase activity of RLIM, suggesting that enzymatic activity of RLIM is required for normal development, cognition and behavior.

Small cell lung cancer (SCLC) patient-derived xenografts (PDX) can be generated from biopsies or circulating tumor cells (CTC), though scarcity of tissue and low efficiency of tumor growth have previously limited these approaches. Applying an established clinical–translational pipeline for tissue collection and an automated microfluidic platform for CTC enrichment, we generated 17 biopsy-derived PDXs and 17 CTC-derived PDXs in a 2-year timeframe, at 89% and 38% efficiency, respectively. Whole-exome sequencing showed that somatic alterations are stably maintained between patient tumors and PDXs. Early-passage PDXs maintain the genomic and transcriptional profiles of the founder PDX. In vivo treatment with etoposide and platinum (EP) in 30 PDX models demonstrated greater sensitivity in PDXs from EP-naïve patients, and resistance to EP corresponded to increased expression of a MYC gene signature. Finally, serial CTC-derived PDXs generated from an individual patient at multiple time points accurately recapitulated the evolving drug sensitivities of that patient's disease. Collectively, this work highlights the translational potential of this strategy. Significance: Effective translational research utilizing SCLC PDX models requires both efficient generation of models from patients and fidelity of those models in representing patient tumor characteristics. We present approaches for efficient generation of PDXs from both biopsies and CTCs, and demonstrate that these models capture the mutational landscape and functional features of the donor tumors. Cancer Discov; 8(5); 600–15. ©2018 AACR. This article is highlighted in the In This Issue feature, p. 517

Pulmonary large-cell neuroendocrine carcinomas (LCNECs) have similarities with other lung cancers, but their precise relationship has remained unclear. Here we perform a comprehensive genomic (n = 60) and transcriptomic (n = 69) analysis of 75 LCNECs and identify two molecular subgroups: “type I LCNECs” with bi-allelic TP53 and STK11/KEAP1 alterations (37%), and “type II LCNECs” enriched for bi-allelic inactivation of TP53 and RB1 (42%). Despite sharing genomic alterations with adenocarcinomas and squamous cell carcinomas, no transcriptional relationship was found; instead LCNECs form distinct transcriptional subgroups with closest similarity to SCLC. While type I LCNECs and SCLCs exhibit a neuroendocrine profile with ASCL1^high/DLL3^high/NOTCH^low, type II LCNECs bear TP53 and RB1 alterations and differ from most SCLC tumors with reduced neuroendocrine markers, a pattern of ASCL1^low/DLL3^low/NOTCH^high, and an upregulation of immune-related pathways. In conclusion, LCNECs comprise two molecularly defined subgroups, and distinguishing them from SCLC may allow stratified targeted treatment of high-grade neuroendocrine lung tumors.

Genetic epilepsies are caused by mutations in a range of different genes, many of them encoding ion channels, receptors or transporters. While the number of detected variants and genes increased dramatically in the recent years, pleiotropic effects have also been recognized, revealing that clinical syndromes with various degrees of severity arise from a single gene, a single mutation, or from different mutations showing similar functional defects. Accordingly, several genes coding for GABAA receptor subunits have been linked to a spectrum of benign to severe epileptic disorders and it was shown that a loss of function presents the major correlated pathomechanism. Here, we identified six variants in GABRA3 encoding the α3-subunit of the GABAA receptor. This gene is located on chromosome Xq28 and has not been previously associated with human disease. Five missense variants and one microduplication were detected in four families and two sporadic cases presenting with a range of epileptic seizure types, a varying degree of intellectual disability and developmental delay, sometimes with dysmorphic features or nystagmus. The variants co-segregated mostly but not completely with the phenotype in the families, indicating in some cases incomplete penetrance, involvement of other genes, or presence of phenocopies. Overall, males were more severely affected and there were three asymptomatic female mutation carriers compared to only one male without a clinical phenotype. X-chromosome inactivation studies could not explain the phenotypic variability in females. Three detected missense variants are localized in the extracellular GABA-binding NH2-terminus, one in the M2-M3 linker and one in the M4 transmembrane segment of the α3-subunit. Functional studies in Xenopus laevis oocytes revealed a variable but significant reduction of GABA-evoked anion currents for all mutants compared to wild-type receptors. The degree of current reduction correlated partially with the phenotype. The microduplication disrupted GABRA3 expression in fibroblasts of the affected patient. In summary, our results reveal that rare loss-of-function variants in GABRA3 increase the risk for a varying combination of epilepsy, intellectual disability/developmental delay and dysmorphic features, presenting in some pedigrees with an X-linked inheritance pattern. 10.1093/brain/awx236_video1 awx236media1 5636589232001

Here we present the 15 pseudochromosomes of sweet potato, Ipomoea batatas, the seventh most important crop in the world and the fourth most significant in China. By using a novel haplotyping method based on genome assembly, we have produced a half haplotype-resolved genome from ~296 Gb of paired-end sequence reads amounting to roughly 67-fold coverage. By phylogenetic tree analysis of homologous chromosomes, it was possible to estimate the time of two recent whole-genome duplication events as occurring about 0.8 and 0.5 million years ago. This half haplotype-resolved hexaploid genome represents the first successful attempt to investigate the complexity of chromosome sequence composition directly in a polyploid genome, using sequencing of the polyploid organism itself rather than any of its simplified proxy relatives. Adaptation and application of our approach should provide higher resolution in future genomic structure investigations, especially for similarly complex genomes.

Homeotic genes code for key transcription factors (HOX-TFs) that pattern the animal body plan. During embryonic development, Hox genes are expressed in overlapping patterns and function in a partially redundant manner. In vitro biochemical screens probing the HOX-TF sequence specificity revealed largely overlapping sequence preferences, indicating that co-factors might modulate the biological function of HOX-TFs. However, due to their overlapping expression pattern, high protein homology, and insufficiently specific antibodies, little is known about their genome-wide binding preferences. In order to overcome this problem, we virally expressed tagged versions of limb-expressed posterior HOX genes (HOXA9-13, and HOXD9-13) in primary chicken mesenchymal limb progenitor cells (micromass). We determined the effect of each HOX-TF on cellular differentiation (chondrogenesis) and gene expression and found that groups of HOX-TFs induce distinct regulatory programs. We used ChIP-seq to determine their individual genome-wide binding profiles and identified between 12,721 and 28,572 binding sites for each of the nine HOX-TFs. Principal Component Analysis (PCA) of binding profiles revealed that the HOX-TFs are clustered in two subgroups (Group 1: HOXA/D9, HOXA/D10, HOXD12, and HOXA13 and Group 2: HOXA/D11 and HOXD13), which are characterized by differences in their sequence specificity and by the presence of cofactor motifs. Specifically, we identified CTCF binding sites in Group 1, indicating that this subgroup of HOX-proteins cooperates with CTCF. We confirmed this interaction by an independent biological assay (Proximity Ligation Assay) and demonstrated that CTCF is a novel HOX cofactor that specifically associates with Group 1 HOX-TFs, pointing towards a possible interplay between HOX-TFs and chromatin architecture. Hox genes encode transcription factors that determine the vertebrate body plan and pattern structures and organs in the developing embryo. Despite decades of effort and research on Hox genes, little is known about the HOX-DNA binding properties in vivo. This lack of knowledge is mainly due to the absence of appropriate antibodies to distinguish between different HOX transcription factors. Here, we adapt a cell culture system that allows us to investigate HOX-DNA binding on a genome-wide scale. With this approach, we define and directly compare the genome-wide binding sites of nine posterior HOXA and HOXD transcription factors. We report that the in vivo HOX binding specificity differs from the in vitro specificity and find that HOX-TFs largely rely on co-factor binding and not only on direct HOX-DNA binding. Finally, we identify a novel HOX co-factor, a genome architecture protein, CTCF suggesting a possible interplay between HOX-TF function and chromatin architecture.