We used two-dimensional quantitative trait locus analysis to identify interacting genetic loci that contribute to the native airway constrictor hyperresponsiveness to methacholine that characterizes A/J mice, relative to C57BL/6J mice. We quantified airway responsiveness to intravenous methacholine boluses in eighty-eight (C57BL/6J X A/J) F₂ and twenty-seven (A/J X C57BL/6J) F₂ mice as well as ten A/J mice and six C57BL/6J mice; all studies were performed in male mice. Mice were genotyped at 384 SNP markers, and from these data two-QTL analyses disclosed one pair of interacting loci on chromosomes 11 and 18; the homozygous A/J genotype at each locus constituted the genetic interaction linked to the hyperresponsive A/J phenotype. Bioinformatic network analysis of potential interactions among proteins encoded by genes in the linked regions disclosed two high priority subnetworks--Myl7, Rock1, Limk2; and Npc1, Npc1l1. Evidence in the literature supports the possibility that either or both networks could contribute to the regulation of airway constrictor responsiveness. Together, these results should stimulate evaluation of the genetic contribution of these networks in the regulation of airway responsiveness in humans.

The human Y-chromosome does not recombine across its male-specific part and is therefore an excellent marker of human migrations. It also plays an important role in male fertility. However, its evolution is difficult to fully understand because of repetitive sequences, inverted repeats and the potentially large role of gene conversion. Here we perform an evolutionary analysis of 62 Y-chromosomes of Danish descent sequenced using a wide range of library insert sizes and high coverage, thus allowing large regions of these chromosomes to be well assembled. These include 17 father-son pairs, which we use to validate variation calling. Using a recent method that can integrate variants based on both mapping and de novo assembly, we genotype 10898 SNVs and 2903 indels (max length of 27241 bp) in our sample and show by father-son concordance and experimental validation that the non-recurrent SNP and indel variation on the Y chromosome tree is called very accurately. This includes variation called in a 0.9 Mb centromeric heterochromatic region, which is by far the most variable in the Y chromosome. Among the variation is also longer sequence-stretches not present in the reference genome but shared with the chimpanzee Y chromosome. We analyzed 2.7 Mb of large inverted repeats (palindromes) for variation patterns among the two palindrome arms and identified 603 mutation and 416 gene conversions events. We find clear evidence for GC-biased gene conversion in the palindromes (and a balancing AT mutation bias), but irrespective of this, also a strong bias towards gene conversion towards the ancestral state, suggesting that palindromic gene conversion may alleviate Muller's ratchet. Finally, we also find a large number of large-scale gene duplications and deletions in the palindromic regions (at least 24) and find that such events can consist of complex combinations of simultaneous insertions and deletions of long stretches of the Y chromosome.

Aneuploidy is a hallmark of cancer and underlies genetic disorders characterized by severe developmental defects, yet the molecular mechanisms explaining its effects on cellular physiology remain elusive. Here we show, using a series of human cells with defined aneuploid karyotypes, that gain of a single chromosome increases genomic instability. Next-generation sequencing and SNP-array analysis reveal accumulation of chromosomal rearrangements in aneuploids, with break point junction patterns suggestive of replication defects. Trisomic and tetrasomic cells also show increased DNA damage and sensitivity to replication stress. Strikingly, we find that aneuploidy-induced genomic instability can be explained by the reduced expression of the replicative helicase MCM2-7. Accordingly, restoring near-wild-type levels of chromatin-bound MCM helicase partly rescues the genomic instability phenotypes. Thus, gain of chromosomes triggers replication stress, thereby promoting genomic instability and possibly contributing to tumorigenesis.

The genus Liriodendron belongs to the family Magnoliaceae, which resides within the magnoliids, an early diverging lineage of the Mesangiospermae. However, the phylogenetic relationship of magnoliids with eudicots and monocots has not been conclusively resolved and thus remains to be determined1-6. Liriodendron is a relict lineage from the Tertiary with two distinct species-one East Asian (L. chinense (Hemsley) Sargent) and one eastern North American (L. tulipifera Linn)-identified as a vicariad species pair. However, the genetic divergence and evolutionary trajectories of these species remain to be elucidated at the whole-genome level7. Here, we report the first de novo genome assembly of a plant in the Magnoliaceae, L. chinense. Phylogenetic analyses suggest that magnoliids are sister to the clade consisting of eudicots and monocots, with rapid diversification occurring in the common ancestor of these three lineages. Analyses of population genetic structure indicate that L. chinense has diverged into two lineages-the eastern and western groups-in China. While L. tulipifera in North America is genetically positioned between the two L. chinense groups, it is closer to the eastern group. This result is consistent with phenotypic observations that suggest that the eastern and western groups of China may have diverged long ago, possibly before the intercontinental differentiation between L. chinense and L. tulipifera. Genetic diversity analyses show that L. chinense has tenfold higher genetic diversity than L. tulipifera, suggesting that the complicated regions comprising east-west-orientated mountains and the Yangtze river basin (especially near 30° N latitude) in East Asia offered more successful refugia than the south-north-orientated mountain valleys in eastern North America during the Quaternary glacial period.

The auxin-inducible degron (AID) system enables rapid depletion of target proteins within the cell by applying the natural auxin IAA. The AID system is useful for investigating the physiological functions of essential proteins; however, this system generally requires high dose of auxin to achieve effective depletion in vertebrate cells. Here, we describe a super-sensitive AID system that incorporates the synthetic auxin derivative 5-Ad-IAA and its high-affinity-binding partner OsTIR1F74A. The super-sensitive AID system enabled more than a 1000-fold reduction of the AID inducer concentrations in chicken DT40 cells. To apply this system to various mammalian cell lines including cancer cells containing multiple sets of chromosomes, we utilized a single-step method where CRISPR/Cas9-based gene knockout is combined with insertion of a pAID plasmid. The single-step method coupled with the super-sensitive AID system enables us to easily and rapidly generate AID-based conditional knockout cells in a wide range of vertebrate cell lines. Our improved method that incorporates the super-sensitive AID system and the single-step method provides a powerful tool for elucidating the roles of essential genes.

It is widely known that castration has a significant effect on the accumulation of adipose tissue. microRNAs (miRNAs) are known to be involved in fat deposition and to be regulated by the androgen-induced androgen receptor (AR). However, there is little understanding of the relationship between miRNAs and fat deposition after castration. In this study, the high-throughput SOLiD sequencing approach was used to identify and characterize miRNA expression in backfat from intact and castrated full-sib male 23-week-old pigs. The patterns of adipogenesis and fat deposition were compared between castrated and intact male pigs.

LC16m8 (m8), a highly attenuated vaccinia virus (VAC) strain, was developed as a smallpox vaccine, and its safety and immunogenicity have been confirmed. Here, we aimed to develop a system that recovers infectious m8 from a bacterial artificial chromosome (BAC) that retains the full-length viral genomic DNA (m8-BAC system). The infectious virus was successfully recovered from a VAC-BAC plasmid, named pLC16m8-BAC. Furthermore, the bacterial replicon-free virus was generated by intramolecular homologous recombination and was successfully recovered from a modified VAC-BAC plasmid, named pLC16m8.8S-BAC. Also, the growth of the recovered virus was indistinguishable from that of authentic m8. The full genome sequence of the plasmid, which harbors identical inverted terminal repeats (ITR) to that of authentic m8, was determined by long-read next-generation sequencing (NGS). The ITR contains x 18 to 32 of the 70 and x 30 to 45 of 54 base pair tandem repeats, and the number of tandem repeats was different between the ITR left and right. Since the virus recovered from pLC16m8.8S-BAC was expected to retain the identical viral genome to that of m8, including the ITR, a reference-based alignment following a short-read NGS was performed to validate the sequence of the recovered virus. Based on the pattern of coverage depth in the ITR, no remarkable differences were observed between the virus and m8, and the other region was confirmed to be identical as well. In summary, this new system can recover the virus, which is geno- and phenotypically indistinguishable from authentic m8.

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.

Chromosome-level assemblies are indispensable for accurate gene prediction, synteny assessment, and understanding higher-order genome architecture. Reference and draft genomes of key helminth species have been published, but little is yet known about the biology of their chromosomes. Here, we present the complete genome of the tapeworm Hymenolepis microstoma, providing a reference quality, end-to-end assembly that represents the first fully assembled genome of a spiralian/lophotrochozoan, revealing new insights into chromosome evolution.

To investigate how exogenous DNA concatemerizes to form episomal artificial chromosomes (ACs), acquire equal segregation ability and maintain stable holocentromeres, we injected DNA sequences with different features, including sequences that are repetitive or complex, and sequences with different AT-contents, into the gonad of Caenorhabditis elegans to form ACs in embryos, and monitored AC mitotic segregation. We demonstrated that AT-poor sequences (26% AT-content) delayed the acquisition of segregation competency of newly formed ACs. We also co-injected fragmented Saccharomyces cerevisiae genomic DNA, differentially expressed fluorescent markers and ubiquitously expressed selectable marker to construct a less repetitive, more complex AC. We sequenced the whole genome of a strain which propagates this AC through multiple generations, and de novo assembled the AC sequences. We discovered CENP-AHCP-3 domains/peaks are distributed along the AC, as in endogenous chromosomes, suggesting a holocentric architecture. We found that CENP-AHCP-3 binds to the unexpressed marker genes and many fragmented yeast sequences, but is excluded in the yeast extremely high-AT-content centromeric and mitochondrial DNA (> 83% AT-content) on the AC. We identified A-rich motifs in CENP-AHCP-3 domains/peaks on the AC and on endogenous chromosomes, which have some similarity with each other and similarity to some non-germline transcription factor binding sites.

The human oral microbiome is formed early in development. Its composition is influenced by environmental factors including diet, substance use, oral health, and overall health and disease. The influence of human genes on the composition and stability of the oral microbiome is still poorly understood. We studied both environmental and genetic characteristics on the oral microbiome in a large twin sample as well as in a large cohort of unrelated individuals. We identify several significantly heritable features of the oral microbiome. The heritability persists in twins even when their cohabitation changes. The heritability of these traits correlates with the cumulative genetic contributions of over half a million single nucleotide sequence variants measured in a different population of unrelated individuals. Comparison of same-sex and opposite sex cotwins showed no significant differences. We show that two new loci on chromosomes 7 and 12 are associated with the most heritable traits.

Fatty acid-binding proteins (FABPs) are members of the intracellular lipid-binding protein (iLBP) family and are involved in reversibly binding intracellular hydrophobic ligands and trafficking them throughout cellular compartments, including the peroxisomes, mitochondria, endoplasmic reticulum and nucleus. FABPs are small, structurally conserved cytosolic proteins consisting of a water-filled, interior-binding pocket surrounded by ten anti-parallel beta sheets, forming a beta barrel. At the superior surface, two alpha-helices cap the pocket and are thought to regulate binding. FABPs have broad specificity, including the ability to bind long-chain (C16-C20) fatty acids, eicosanoids, bile salts and peroxisome proliferators. FABPs demonstrate strong evolutionary conservation and are present in a spectrum of species including Drosophila melanogaster, Caenorhabditis elegans, mouse and human. The human genome consists of nine putatively functional protein-coding FABP genes. The most recently identified family member, FABP12, has been less studied.

During the life cycle of retroviruses, establishment of a productive infection requires stable joining of a DNA copy of the viral RNA genome into host cell chromosomes. Retroviruses are thus promising vectors for the efficient and stable delivery of genes in therapeutic protocols. Integration of retroviral DNA is catalyzed by the viral enzyme integrase (IN), and one salient feature of retroviral DNA integration is its lack of specificity, as many chromosomal sites can serve as targets for integration. Despite the promise for success in the clinic, one major drawback of the retrovirus-based vector is that any unintended insertion events from the therapy can potentially lead to deleterious effects in patients, as demonstrated by the development of malignancies in both animal and human studies. One approach to directing integration into predetermined DNA sites is fusing IN to a sequence-specific DNA-binding protein, which results in a bias of integration near the recognition site of the fusion partner. Encouraging results have been generated in vitro and in vivo using fusion protein constructs of human immunodeficiency virus type 1 IN and E2C, a designed polydactyl zinc-finger protein that specifically recognizes an 18-base pair DNA sequence. This review focuses on the method for preparing infectious virions containing the IN fusion proteins and on the quantitative PCR assays for determining integration site specificity. Efforts to engineer IN to recognize specific target DNA sequences within the genome may lead to development of effective retroviral vectors that can safely deliver gene-based therapeutics in a clinical setting.

Human papillomavirus (HPV) causes 5% of all cancers and frequently integrates into host chromosomes. The HPV oncoproteins E6 and E7 are necessary but insufficient for cancer formation, indicating that additional secondary genetic events are required. Here, we investigate potential oncogenic impacts of virus integration. Analysis of 105 HPV-positive oropharyngeal cancers by whole-genome sequencing detects virus integration in 77%, revealing five statistically significant sites of recurrent integration near genes that regulate epithelial stem cell maintenance (i.e., SOX2, TP63, FGFR, MYC) and immune evasion (i.e., CD274). Genomic copy number hyperamplification is enriched 16-fold near HPV integrants, and the extent of focal host genomic instability increases with their local density. The frequency of genes expressed at extreme outlier levels is increased 86-fold within ±150 kb of integrants. Across 95% of tumors with integration, host gene transcription is disrupted via intragenic integrants, chimeric transcription, outlier expression, gene breaking, and/or de novo expression of noncoding or imprinted genes. We conclude that virus integration can contribute to carcinogenesis in a large majority of HPV-positive oropharyngeal cancers by inducing extensive disruption of host genome structure and gene expression.

Tenosynovial giant cell tumors (TSGCTs) are rare tumors arising in tendons or the synoviae of joints and bursae. The localized type is benign while the diffuse type shows expansive growth leading to greater morbidity and is therefore considered locally aggressive. Typical recurrent chromosomal aberrations are found in the majority of TSCGT and the CSF1 gene is frequently involved. In this article, we describe a newly identified gene fusion mediated by an inversion in a case of diffuse TSGCT. Multicolor-fluorescence in situ hybridization (FISH) molecular karyotyping identified a pericentric inversion of chromosome 1 in 7 out of 17 analyzed cells 46,XX,inv(1)(p13.3q24.3) [7]/46,XX [10], and with interphase FISH the involvement the CSF1 locus was detected. After performing transcriptome sequencing analysis for fusion detection, only one out of five fusion gene algorithms detected a fusion involving the CSF1 gene product. The resulting chimera fuses a sequence from a human endogenous retrovirus (HERV) gene to CSF1 Exon 6 on chromosome 1, abrogating the regulatory element of the 3' untranslated region of the CSF1 gene. This new translocation involving Exon 6 of the CSF1 gene fused to 1q24.1, supports the hypothesis that a mutated CSF1 protein is likely to play a vital role in the pathogenesis of TSGCT. The role of the HERV partner identified as a translocation partner, however, remains unclear. Our data add to the complexity of involved translocation partners in TSGCT and point to the potential difficulty of identifying fusion partners in tumor diagnostics using transcriptome sequencing when HERV or other repeat elements are involved.

Viruses need to hijack the translational machinery of the host cell for a productive infection to happen. However, given the dynamic landscape of tRNA pools among tissues, it is unclear whether different viruses infecting different tissues have adapted their codon usage toward their tropism. Here, we collect the coding sequences of 502 human-infecting viruses and determine that tropism explains changes in codon usage. Using the tRNA abundances across 23 human tissues from The Cancer Genome Atlas (TCGA), we build an in silico model of translational efficiency that validates the correspondence of the viral codon usage with the translational machinery of their tropism. For instance, we detect that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is specifically adapted to the upper respiratory tract and alveoli. Furthermore, this correspondence is specifically defined in early viral proteins. The observed tissue-specific translational efficiency could be useful for the development of antiviral therapies and vaccines.

Studies of the Y chromosome in primates, rodents and carnivores provide compelling evidence that the male specific region of Y (MSY) contains functional genes, many of which have specialized roles in spermatogenesis and male-fertility. Little similarity, however, has been found between the gene content and sequence of MSY in different species. This hinders the discovery of species-specific male fertility genes and limits our understanding about MSY evolution in mammals. Here, a detailed MSY gene catalogue was developed for the horse--an odd-toed ungulate. Using direct cDNA selection from horse testis, and sequence analysis of Y-specific BAC clones, 37 horse MSY genes/transcripts were identified. The genes were mapped to the MSY BAC contig map, characterized for copy number, analyzed for transcriptional profiles by RT-PCR, examined for the presence of ORFs, and compared to other mammalian orthologs. We demonstrate that the horse MSY harbors 20 X-degenerate genes with known orthologs in other eutherian species. The remaining 17 genes are acquired or novel and have so far been identified only in the horse or donkey Y chromosomes. Notably, 3 transcripts were found in the heterochromatic part of the Y. We show that despite substantial differences between the sequence, gene content and organization of horse and other mammalian Y chromosomes, the functions of MSY genes are predominantly related to testis and spermatogenesis. Altogether, 10 multicopy genes with testis-specific expression were identified in the horse MSY, and considered likely candidate genes for stallion fertility. The findings establish an important foundation for the study of Y-linked genetic factors governing fertility in stallions, and improve our knowledge about the evolutionary processes that have shaped Y chromosomes in different mammalian lineages.

We report the analysis of 335 microsatellite markers genotyped in 110 multiplex families with autism. All families include at least two "affected" siblings, at least one of whom has autism; the remaining affected sibs carry diagnoses of either Asperger syndrome or pervasive developmental disorder. Affected sib-pair analysis yielded multipoint maximum LOD scores (MLS) that reach the accepted threshold for suggestive linkage on chromosomes 5, X, and 19. Nominal evidence for linkage (point-wise P<.05) was obtained on chromosomes 2, 3, 4, 8, 10, 11, 12, 15, 16, 18, and 20, and secondary loci were found on chromosomes 5 and 19. Analysis of families sharing alleles at the putative X chromosomal linked locus and one or more other putative linked loci produced an MLS of 3.56 for the DXS470-D19S174 marker combination. In an effort to increase power to detect linkage, scan statistics were used to evaluate the significance of peak LOD scores based on statistical evidence at adjacent marker loci. This analysis yielded impressive evidence for linkage to autism and autism-spectrum disorders with significant genomewide P values <.05 for markers on chromosomes 5 and 8 and with suggestive linkage evidence for a marker on chromosome 19.

Sex chromosomes evolve once recombination is halted between a homologous pair of chromosomes. The dominant model of sex chromosome evolution posits that recombination is suppressed between emerging X and Y chromosomes in order to resolve sexual conflict. Here we test this model using whole genome and transcriptome resequencing data in the guppy, a model for sexual selection with many Y-linked colour traits. We show that although the nascent Y chromosome encompasses nearly half of the linkage group, there has been no perceptible degradation of Y chromosome gene content or activity. Using replicate wild populations with differing levels of sexually antagonistic selection for colour, we also show that sexual selection leads to greater expansion of the non-recombining region and increased Y chromosome divergence. These results provide empirical support for longstanding models of sex chromosome catalysis, and suggest an important role for sexual selection and sexual conflict in genome evolution.

Androgenetic complete hydatidiform moles are human pregnancies with no embryos and affect 1 in every 1,400 pregnancies. They have mostly androgenetic monospermic genomes with all the chromosomes originating from a haploid sperm and no maternal chromosomes. Androgenetic complete hydatidiform moles were described in 1977, but how they occur has remained an open question. We identified bi-allelic deleterious mutations in MEI1, TOP6BL/C11orf80, and REC114, with roles in meiotic double-strand breaks formation in women with recurrent androgenetic complete hydatidiform moles. We investigated the occurrence of androgenesis in Mei1-deficient female mice and discovered that 8% of their oocytes lose all their chromosomes by extruding them with the spindles into the first polar body. We demonstrate that Mei1-/- oocytes are capable of fertilization and 5% produce androgenetic zygotes. Thus, we uncover a meiotic abnormality in mammals and a mechanism for the genesis of androgenetic zygotes that is the extrusion of all maternal chromosomes and their spindles into the first polar body.