The nuclear fertility restorer gene Rf5 in HA-R9, originating from the wild sunflower species Helianthus annuus, is able to restore the widely used PET1 cytoplasmic male sterility in sunflowers. Previous mapping placed Rf5 at an interval of 5.8 cM on sunflower chromosome 13, distal to a rust resistance gene R11 at a 1.6 cM genetic distance in an SSR map. In the present study, publicly available SNP markers were further mapped around Rf5 and R11 using 192 F2 individuals, reducing the Rf5 interval from 5.8 to 0.8 cM. Additional SNP markers were developed in the target region of the two genes from the whole-genome resequencing of HA-R9, a donor line carrying Rf5 and R11. Fine mapping using 3517 F3 individuals placed Rf5 at a 0.00071 cM interval and the gene co-segregated with SNP marker S13_216392091. Similarly, fine mapping performed using 8795 F3 individuals mapped R11 at an interval of 0.00210 cM, co-segregating with two SNP markers, S13_225290789 and C13_181790141. Sequence analysis identified Rf5 as a pentatricopeptide repeat-encoding gene. The high-density map and diagnostic SNP markers developed in this study will accelerate the use of Rf5 and R11 in sunflower breeding.

Walking is a simple form of exercise, widely promoted for its health benefits. Self-reported walking pace has been associated with a range of cardiorespiratory and cancer outcomes, and is a strong predictor of mortality. Here we perform a genome-wide association study of self-reported walking pace in 450,967 European ancestry UK Biobank participants. We identify 70 independent associated loci (P < 5 × 10-8), 11 of which are novel. We estimate the SNP-based heritability as 13.2% (s.e. = 0.21%), reducing to 8.9% (s.e. = 0.17%) with adjustment for body mass index. Significant genetic correlations are observed with cardiometabolic, respiratory and psychiatric traits, educational attainment and all-cause mortality. Mendelian randomization analyses suggest a potential causal link of increasing walking pace with a lower cardiometabolic risk profile. Given its low heritability and simple measurement, these findings suggest that self-reported walking pace is a pragmatic target for interventions aiming for general benefits on health.

A new method to improve the efficiency of flanking sequence identification by genome walking was developed based on an expanded, sequential list of criteria for selecting candidate enzymes, plus several other optimization steps. These criteria include: step (1) initially choosing the most appropriate restriction enzyme according to the average fragment size produced by each enzyme determined using in silico digestion of genomic DNA, step (2) evaluating the in silico frequency of fragment size distribution between individual chromosomes, step (3) selecting those enzymes that generate fragments with the majority between 100 bp and 3,000 bp, step (4) weighing the advantages and disadvantages of blunt-end sites vs. cohesive-end sites, step (5) elimination of methylation sensitive enzymes with methylation-insensitive isoschizomers, and step (6) elimination of enzymes with recognition sites within the binary vector sequence (T-DNA and plasmid backbone). Step (7) includes the selection of a second restriction enzyme with highest number of recognition sites within regions not covered by the first restriction enzyme. Step (8) considers primer and adapter sequence optimization, selecting the best adapter-primer pairs according to their hairpin/dimers and secondary structure. In step (9), the efficiency of genomic library development was improved by column-filtration of digested DNA to remove restriction enzyme and phosphatase enzyme, and most important, to remove small genomic fragments (<100 bp) lacking the T-DNA insertion, hence improving the chance of ligation between adapters and fragments harbouring a T-DNA. Two enzymes, NsiI and NdeI, fit these criteria for the Arabidopsis thaliana genome. Their efficiency was assessed using 54 T(3) lines from an Arabidopsis SK enhancer population. Over 70% success rate was achieved in amplifying the flanking sequences of these lines. This strategy was also tested with Brachypodium distachyon to demonstrate its applicability to other larger genomes.

Chromosome organization is crucial for genome function. Here, we present a method for visualizing chromosomal DNA at super-resolution and then integrating Hi-C data to produce three-dimensional models of chromosome organization. Using the super-resolution microscopy methods of OligoSTORM and OligoDNA-PAINT, we trace 8 megabases of human chromosome 19, visualizing structures ranging in size from a few kilobases to over a megabase. Focusing on chromosomal regions that contribute to compartments, we discover distinct structures that, in spite of considerable variability, can predict whether such regions correspond to active (A-type) or inactive (B-type) compartments. Imaging through the depths of entire nuclei, we capture pairs of homologous regions in diploid cells, obtaining evidence that maternal and paternal homologous regions can be differentially organized. Finally, using restraint-based modeling to integrate imaging and Hi-C data, we implement a method-integrative modeling of genomic regions (IMGR)-to increase the genomic resolution of our traces to 10 kb.

Walking catfish (Clarias batrachus) is a freshwater fish capable of air-breathing and locomotion on land. It usually inhabits various low-oxygen habitats, burrows inside the mudflat, and sometimes "walks" to search for suitable environments during summer. It has evolved accessory air-breathing organs for respiring air and corresponding mechanisms to survive in such challenging environments. Thereby, it serves as a great model for understanding adaptations to terrestrial life.

The generation of BAC/PAC contigs in targeted genome regions is a powerful method to establish high-resolution physical maps. In domestic animal species the generation of such contigs is typically initiated with the screening of libraries with probes derived from human genes that are expected to be located in the region of interest by comparative mapping. However, in many instances the available gene-derived probes are too far apart to allow the cloning of BAC/PAC contigs larger than a few hundred kb. High resolution physical mapping allows to estimate the sizes of gaps and to control the orientation of the individual sub-contigs, which helps to avoid errors during the assembly of smaller contigs into final Mb-sized contigs. The recently constructed porcine IMNpRH2 panel allowed us to use this approach for the construction of high-resolution physical maps of SSC 6q1.2.

We describe persistent one-way walking of Drosophila melanogaster in a circular arena. Wild-type Canton-S adult flies walked in one direction, counter-clockwise or clockwise, for minutes, whereas white-eyed mutant [Formula: see text] changed directions frequently. Locomotion in the circular arena could be classified into four components: counter-clockwise walking, clockwise walking, nondirectional walking and pausing. Genetic analysis revealed that while wild-type genetic background was associated with reduced directional change and reduced numbers of one-way (including counter-clockwise and clockwise) and nondirectional walks, the white ([Formula: see text]) locus promoted persistent one-way walking by increasing the maximal duration of one-way episodes. The promoting effect of [Formula: see text] was further supported by the observations that (1) [Formula: see text] duplicated to the Y chromosome, (2) four genomic copies of mini-white inserted on the autosomes, and (3) pan-neuronal overexpression of the White protein increased the maximal duration of one-way episodes, and that RNAi knockdown of [Formula: see text] in the neurons decreased the maximal duration of one-way episodes. These results suggested a pleiotropic function of [Formula: see text] in promoting persistent one-way walking in the circular arena.

Large repeat rich genomes present challenges for assembly using short read technologies. The 32 Gb axolotl genome is estimated to contain ~19 Gb of repetitive DNA making an assembly from short reads alone effectively impossible. Indeed, this model species has been sequenced to 20× coverage but the reads could not be conventionally assembled. Using an alternative strategy, we have assembled subsets of these reads into scaffolds describing over 19,000 gene models. We call this method Virtual Genome Walking as it locally assembles whole genome reads based on a reference transcriptome, identifying exons and iteratively extending them into surrounding genomic sequence. These assemblies are then linked and refined to generate gene models including upstream and downstream genomic, and intronic, sequence. Our assemblies are validated by comparison with previously published axolotl bacterial artificial chromosome (BAC) sequences. Our analyses of axolotl intron length, intron-exon structure, repeat content and synteny provide novel insights into the genic structure of this model species. This resource will enable new experimental approaches in axolotl, such as ChIP-Seq and CRISPR and aid in future whole genome sequencing efforts. The assembled sequences and annotations presented here are freely available for download from https://tinyurl.com/y8gydc6n . The software pipeline is available from https://github.com/LooseLab/iterassemble .

Djulis (Chenopodium formosanum Koidz.) is a crop grown since antiquity in Taiwan. It is a BCD-genome hexaploid (2n = 6x = 54) domesticated form of lambsquarters (C. album L.) and a relative of the allotetraploid (AABB) C. quinoa. As with quinoa, djulis seed contains a complete protein profile and many nutritionally important vitamins and minerals. While still sold locally in Taiwanese markets, its traditional culinary uses are being lost as diets of younger generations change. Moreover, indigenous Taiwanese peoples who have long safeguarded djulis are losing their traditional farmlands. We used PacBio sequencing and Hi-C-based scaffolding to produce a chromosome-scale, reference-quality assembly of djulis. The final genome assembly spans 1.63 Gb in 798 scaffolds, with 97.8% of the sequence contained in 27 scaffolds representing the nine haploid chromosomes of each sub-genome of the species. Benchmarking of universal, single-copy orthologs indicated that 98.5% of the conserved orthologous genes for Viridiplantae are complete within the assembled genome, with 92.9% duplicated, as expected for a polyploid. A total of 67.8% of the assembly is repetitive, with the most common repeat being Gypsy long terminal repeat retrotransposons, which had significantly expanded in the B sub-genome. Gene annotation using Iso-Seq data from multiple tissues identified 75,056 putative gene models. Comparisons to quinoa showed strong patterns of synteny which allowed for the identification of homoeologous chromosomes, and sub-genome-specific sequences were used to assign homoeologs to each sub-genome. These results represent the first hexaploid genome assembly and the first assemblies of the C and D genomes of the Chenopodioideae subfamily.

Eukaryotic pathogens of humans often evade the immune system by switching the expression of surface proteins encoded by subtelomeric gene families. To determine if plant pathogenic fungi use a similar mechanism to avoid host defenses, we sequenced the 14 chromosome ends of the rice blast pathogen, Magnaporthe oryzae. One telomere is directly joined to ribosomal RNA-encoding genes, at the end of the approximately 2 Mb rDNA array. Two are attached to chromosome-unique sequences, and the remainder adjoin a distinct subtelomere region, consisting of a telomere-linked RecQ-helicase (TLH) gene flanked by several blocks of tandem repeats. Unlike other microbes, M.oryzae exhibits very little gene amplification in the subtelomere regions-out of 261 predicted genes found within 100 kb of the telomeres, only four were present at more than one chromosome end. Therefore, it seems unlikely that M.oryzae uses switching mechanisms to evade host defenses. Instead, the M.oryzae telomeres have undergone frequent terminal truncation, and there is evidence of extensive ectopic recombination among transposons in these regions. We propose that the M.oryzae chromosome termini play more subtle roles in host adaptation by promoting the loss of terminally-positioned genes that tend to trigger host defenses.

Poa annua L. is a globally distributed grass with economic and horticultural significance as a weed and as a turfgrass. This dual significance, and its phenotypic plasticity and ecological adaptation, have made P. annua an intriguing plant for genetic and evolutionary studies. Because of the lack of genomic resources and its allotetraploid (2n = 4x = 28) nature, a reference genome sequence would be a valuable asset to better understand the significance and polyploid origin of P. annua. Here we report a genome assembly with scaffolds representing the 14 haploid chromosomes that are 1.78 Gb in length with an N50 of 112 Mb and 96.7% of BUSCO orthologs. Seventy percent of the genome was identified as repetitive elements, 91.0% of which were Copia- or Gypsy-like long-terminal repeats. The genome was annotated with 76,420 genes spanning 13.3% of the 14 chromosomes. The two subgenomes originating from Poa infirma (Knuth) and Poa supina (Schrad) were sufficiently divergent to be distinguishable but syntenic in sequence and annotation with repetitive elements contributing to the expansion of the P. infirma subgenome.

Long-range and highly accurate de novo assembly from short-read data is one of the most pressing challenges in genomics. Recently, it has been shown that read pairs generated by proximity ligation of DNA in chromatin of living tissue can address this problem, dramatically increasing the scaffold contiguity of assemblies. Here, we describe a simpler approach ("Chicago") based on in vitro reconstituted chromatin. We generated two Chicago data sets with human DNA and developed a statistical model and a new software pipeline ("HiRise") that can identify poor quality joins and produce accurate, long-range sequence scaffolds. We used these to construct a highly accurate de novo assembly and scaffolding of a human genome with scaffold N50 of 20 Mbp. We also demonstrated the utility of Chicago for improving existing assemblies by reassembling and scaffolding the genome of the American alligator. With a single library and one lane of Illumina HiSeq sequencing, we increased the scaffold N50 of the American alligator from 508 kbp to 10 Mbp.

Holoprosencephaly (HPE) is the most common developmental field defect in patterning of the human prosencephalon and associated craniofacial structures. The genetics is complex, with 12 loci defined on 11 chromosomes. We defined a locus for HPE (HPE8) on human chromosome 14q13 between markers D14S49 and AFM205XG5, by mapping deletion intervals of affected subjects with proximal chromosome 14q interstitial cytogenetic deletions. A 35-BAC contig was built by chromosome walking. By annotation of the 2.82-Mb minimal critical region, we identified 28 possible genes. Seven genes were expressed in human fetal brain: NPAS3, SNX6, C14ORF11, C14ORF10, PAX9, NKX2.1, and C14ORF19, the last an apparent gene fragment. Molecular embryology, animal modeling, and human mutation studies were reported elsewhere for PAX9 and NKX2.1. We focused on three genes, SNX6, NPAS3, and C14ORF11, as potential candidates for HPE. Genomic structure, human expression patterns, protein cellular localization, and embryonic expression patterns of orthologous murine genes were determined, showing that the three genes have properties similar to those of known HPE genes.

Substantial morphological variation in land plants remains inaccessible to genetic analysis because current models lack variation in important ecological and agronomic traits. The genus Gilia was historically a model for biosystematics studies and includes variation in morphological traits that are poorly understood at the genetic level. We assembled a chromosome-scale reference genome of G. yorkii and used it to investigate genome evolution in the Polemoniaceae. We performed QTL (quantitative trait loci) mapping in a G. yorkii×G. capitata interspecific population for traits related to inflorescence architecture and flower color. The genome assembly spans 2.75 Gb of the estimated 2.80-Gb genome, with 96.7% of the sequence contained in the nine largest chromosome-scale scaffolds matching the haploid chromosome number. Gilia yorkii experienced at least one round of whole-genome duplication shared with other Polemoniaceae after the eudicot paleohexaploidization event. We identified QTL linked to variation in inflorescence architecture and petal color, including a candidate for the major flower color QTL-a tandem duplication of flavanol 3',5'-hydroxylase. Our results demonstrate the utility of Gilia as a forward genetic model for dissecting the evolution of development in plants including the causal loci underlying inflorescence architecture transitions.

Methicillin resistance determinant mecA is generally transferred by SCCmec elements. However, the mecA gene might not be carried by a SCCmec in a Staphylococcus haemolyticus clinical isolate, WCH1, as no cassette chromosome recombinase genes were detected. Therefore, the genetic context of mecA in WCH1 was investigated.

The walking catfish Clarias magur (Hamilton, 1822) (magur) is an important catfish species inhabiting the Indian subcontinent. It is considered as a highly nutritious food fish and has the capability to walk to some distance, and survive a considerable period without water. Assembly, scaffolding and several rounds of iterations resulted in 3,484 scaffolds covering ∼94% of estimated genome with 9.88 Mb largest scaffold, and N50 1.31 Mb. The genome possessed 23,748 predicted protein encoding genes with annotation of 19,279 orthologous genes. A total of 166 orthologous groups represented by 222 genes were found to be unique for this species. The Computational Analysis of gene Family Evolution (CAFE) analysis revealed expansion of 207 gene families and 100 gene families have rapidly evolved. Genes specific to important environmental and terrestrial adaptation, viz. urea cycle, vision, locomotion, olfactory and vomeronasal receptors, immune system, anti-microbial properties, mucus, thermoregulation, osmoregulation, air-breathing, detoxification, etc. were identified and critically analysed. The analysis clearly indicated that C. magur genome possessed several unique and duplicate genes similar to that of terrestrial or amphibians' counterparts in comparison to other teleostean species. The genome information will be useful in conservation genetics, not only for this species but will also be very helpful in such studies in other catfishes.

Background: Infantile hypotonia with psychomotor retardation and characteristic facies 2 (IHPRF2) is a rare autosomal recessive neurodevelopmental disorder caused by mutations in the UNC80 gene. It is characterized by severe global developmental delay, poor or absent speech and absent or limited walking abilities. The current study explored a case of a Chinese patient with IHPRF2 caused by a novel splicing variant of UNC80. Case Report: The proband is a 8-year-old Chinese male manifested with global developmental delay, severe truncal hypotonia, absent speech and intellectual disability. SNP array analysis revealed a uniparental isodisomy of the entire chromosome 2 [UPD(2)] in the proband. Whole exome sequencing (WES) subsequently identified a novel mutation c.5609-4G>A in the UNC80 gene, which was inherited from his mother and was confirmed by Sanger sequencing, indicating that UPD(2) was of maternal origin. Conclusion: A novel UNC80 homozygous splicing variant c.5609-4G>A associated with maternal UPD(2) was identified. These findings indicate that UPD poses a high risk of autosomal recessive diseases, and provides information on the variant spectrum for UNC80. Our findings elucidate on understanding of the genotype-phenotype associations that occur in IHPRF2 patients.

A key question in developmental biology addresses the mechanism of asymmetric cell division. Asymmetry is crucial for generating cellular diversity required for development in multicellular organisms. As one of the potential mechanisms, chromosomally borne epigenetic difference between sister cells that changes mating/cell type has been demonstrated only in the Schizosaccharomyces pombe fission yeast. For technical reasons, it is nearly impossible to determine the existence of such a mechanism operating during embryonic development of multicellular organisms. Our work addresses whether such an epigenetic mechanism causes asymmetric cell division in the recently sequenced fission yeast, S. japonicus (with 36% GC content), which is highly diverged from the well-studied S. pombe species (with 44% GC content). We find that the genomic location and DNA sequences of the mating-type loci of S. japonicus differ vastly from those of the S. pombe species. Remarkably however, similar to S. pombe, the S. japonicus cells switch cell/mating type after undergoing two consecutive cycles of asymmetric cell divisions: only one among four "granddaughter" cells switches. The DNA-strand-specific epigenetic imprint at the mating-type locus1 initiates the recombination event, which is required for cellular differentiation. Therefore the S. pombe and S. japonicus mating systems provide the first two examples in which the intrinsic chirality of double helical structure of DNA forms the primary determinant of asymmetric cell division. Our results show that this unique strand-specific imprinting/segregation epigenetic mechanism for asymmetric cell division is evolutionary conserved. Motivated by these findings, we speculate that DNA-strand-specific epigenetic mechanisms might have evolved to dictate asymmetric cell division in diploid, higher eukaryotes as well.

We describe a new cell line, Clarias dussumieri fin (ClDuF), from the caudal fin of C. dussumieri using the explant technique followed by cryopreservation. The cryopreserved CiDuF cells were validated for quality and other characteristics. They showed typical epithelial morphology in vitro and epithelial cells outgrew their fibroblast cells after the fifth passage. ClDuF cells had a characteristic sigmoid curve with population doubling in 24 h. Immunotyping of the ClDuF cells against cytokeratin suggested the epithelial lineage. Chromosome analysis showed normal diploid (2n = 50) numbers and the cells did not contain any contamination, including Mycoplasma and other microbes. Partial sequencing of fragments of mitochondrial 16s rRNA and COI genes of ClDuF confirmed that the cell line was initiated from C. dussumieri. Cells at the 10th and 25th passages had more than 80% and 70% viability in the culture, respectively, after 6 months of storage at LN2 . These ClDuF cells were morphologically identical to the cells before freezing and the genetic resource of C. dussumieri was preserved. The species-specific cells can serve as a valuable source for virus isolation, conservation and cloning of somatic cells.

It has been suggested that the haplotypic relationship between microsatellite markers and single nucleotide polymorphisms (SNPs) is of considerable importance, as microsatellite markers can potentially be incorporated into haplotypes containing SNPs to increase marker density across a region of interest. However, SNPs and microsatellite markers have different mutation rates and durations, and it is conceivable that the linkage disequilibrium (LD) patterns between the genetic markers may considerably differ. We assessed the LD patterns using 1,661 SNPs and 65 microsatellite markers along chromosome 22 and investigated whether common patterns of LD between the two genetic markers are deduced from the results. The results demonstrated that the patterns of LD among microsatellite markers varied considerably and the LD runs of SNPs and microsatellite markers showed distinct patterns. Microsatellite markers have a much higher mutation rate and the evolution of microsatellite markers is a more complex process which has distinct mutation properties from those of SNPs. We consider that these might contribute to the different LD patterns between the two genetic markers. Therefore, it would seem inadvisable to make assumptions about persistence of LD across even a relatively small genetic distance among microsatellite markers and to construct mixed marker haplotypes/LD maps employing microsatellite markers.