Complexity of avian evolution revealed by family-level genomes.

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

Arabinosylation of cell wall extensin is required for the directional response to salinity in roots.

Soil salinity is a major contributor to crop yield losses. To improve our understanding of root responses to salinity, we developed and exploited a real-time salt-induced tilting assay. This assay follows root growth upon both gravitropic and salt challenges, revealing that root bending upon tilting is modulated by Na+ ions, but not by osmotic stress. Next, we measured this salt-specific response in 345 natural Arabidopsis (Arabidopsis thaliana) accessions and discovered a genetic locus, encoding the cell wall-modifying enzyme EXTENSIN ARABINOSE DEFICIENT TRANSFERASE (ExAD) that is associated with root bending in the presence of NaCl (hereafter salt). Extensins are a class of structural cell wall glycoproteins known as hydroxyproline (Hyp)-rich glycoproteins, which are posttranslationally modified by O-glycosylation, mostly involving Hyp-arabinosylation. We show that salt-induced ExAD-dependent Hyp-arabinosylation influences root bending responses and cell wall thickness. Roots of exad1 mutant seedlings, which lack Hyp-arabinosylation of extensin, displayed increased thickness of root epidermal cell walls and greater cell wall porosity. They also showed altered gravitropic root bending in salt conditions and a reduced salt-avoidance response. Our results suggest that extensin modification via Hyp-arabinosylation is a unique salt-specific cellular process required for the directional response of roots exposed to salinity.

The tuatara genome reveals ancient features of amniote evolution.

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

Dense sampling of bird diversity increases power of comparative genomics.

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity1-4. Sparse taxon sampling has previously been proposed to confound phylogenetic inference5, and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

Publisher Correction: The tuatara genome reveals ancient features of amniote evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Historic Sampling of a Vanishing Beast: Population Structure and Diversity in the Black Rhinoceros.

The black rhinoceros (Diceros bicornis L.) is a critically endangered species historically distributed across sub-Saharan Africa. Hunting and habitat disturbance have diminished both its numbers and distribution since the 19th century, but a poaching crisis in the late 20th century drove them to the brink of extinction. Genetic and genomic assessments can greatly increase our knowledge of the species and inform management strategies. However, when a species has been severely reduced, with the extirpation and artificial admixture of several populations, it is extremely challenging to obtain an accurate understanding of historic population structure and evolutionary history from extant samples. Therefore, we generated and analyzed whole genomes from 63 black rhinoceros museum specimens collected between 1775 and 1981. Results showed that the black rhinoceros could be genetically structured into six major historic populations (Central Africa, East Africa, Northwestern Africa, Northeastern Africa, Ruvuma, and Southern Africa) within which were nested four further subpopulations (Maasailand, southwestern, eastern rift, and northern rift), largely mirroring geography, with a punctuated north-south cline. However, we detected varying degrees of admixture among groups and found that several geographical barriers, most prominently the Zambezi River, drove population discontinuities. Genomic diversity was high in the middle of the range and decayed toward the periphery. This comprehensive historic portrait also allowed us to ascertain the ancestry of 20 resequenced genomes from extant populations. Lastly, using insights gained from this unique temporal data set, we suggest management strategies, some of which require urgent implementation, for the conservation of the remaining black rhinoceros diversity.

A chromosome phased diploid genome assembly of African hunting dog (Lycaon pictus).

The African hunting dog (Lycaon pictus, 2n=78) once ranged over most sub-Saharan ecosystems except its deserts and rainforests. However as a result of (still ongoing) population declines, today they remain only as small fragmented populations. Furthermore, the future of the species remains unclear, due to both anthropogenic pressure as well as interactions with domestic dogs, thus their preservation is a conservation priority. On the tree of life, the hunting dog is basal to Canis and Cuon and forms a crown group with them, making it a useful species for comparative genomic studies. Here, we present a diploid chromosome level assembly of an African hunting dog. Assembled according to VGP guidelines from a combination of PacBio HiFi reads and HiC data, it is phased at the level of individual chromosomes. The maternal (pseudo)haplotype (mat) of our assembly has a length of 2.38 Gbp, and 99.36 % of the sequence is encompassed by 39 chromosomal scaffolds. The rest is included in only 36 unplaced short scaffolds. At the contig level, mat consists of only 166 contigs with an N50 of 39 Mbp. BUSCO analysis showed 95.4 % completeness based on Сarnivora conservative genes (carnivora_odb10). When compared to other available genomes from subtribe Canina, the quality of the assembly is excellent, typically between the 1st and 3rd depending on the parameter used, and a significant improvement on previously published genomes for the species. We hope this assembly will play an important role in future conservation efforts and comparative studies of canid genomes.

NetSurfP-3.0: accurate and fast prediction of protein structural features by protein language models and deep learning.

Recent advances in machine learning and natural language processing have made it possible to profoundly advance our ability to accurately predict protein structures and their functions. While such improvements are significantly impacting the fields of biology and biotechnology at large, such methods have the downside of high demands in terms of computing power and runtime, hampering their applicability to large datasets. Here, we present NetSurfP-3.0, a tool for predicting solvent accessibility, secondary structure, structural disorder and backbone dihedral angles for each residue of an amino acid sequence. This NetSurfP update exploits recent advances in pre-trained protein language models to drastically improve the runtime of its predecessor by two orders of magnitude, while displaying similar prediction performance. We assessed the accuracy of NetSurfP-3.0 on several independent test datasets and found it to consistently produce state-of-the-art predictions for each of its output features, with a runtime that is up to to 600 times faster than the most commonly available methods performing the same tasks. The tool is freely available as a web server with a user-friendly interface to navigate the results, as well as a standalone downloadable package.

Sequencing and de novo assembly of 150 genomes from Denmark as a population reference.

Hundreds of thousands of human genomes are now being sequenced to characterize genetic variation and use this information to augment association mapping studies of complex disorders and other phenotypic traits. Genetic variation is identified mainly by mapping short reads to the reference genome or by performing local assembly. However, these approaches are biased against discovery of structural variants and variation in the more complex parts of the genome. Hence, large-scale de novo assembly is needed. Here we show that it is possible to construct excellent de novo assemblies from high-coverage sequencing with mate-pair libraries extending up to 20 kilobases. We report de novo assemblies of 150 individuals (50 trios) from the GenomeDenmark project. The quality of these assemblies is similar to those obtained using the more expensive long-read technology. We use the assemblies to identify a rich set of structural variants including many novel insertions and demonstrate how this variant catalogue enables further deciphering of known association mapping signals. We leverage the assemblies to provide 100 completely resolved major histocompatibility complex haplotypes and to resolve major parts of the Y chromosome. Our study provides a regional reference genome that we expect will improve the power of future association mapping studies and hence pave the way for precision medicine initiatives, which now are being launched in many countries including Denmark.

Phenotypic Characterization and Comparative Genomic Analysis of Novel Salmonella Bacteriophages Isolated from a Tropical Rainforest.

Salmonella infections across the globe are becoming more challenging to control due to the emergence of multidrug-resistant (MDR) strains. Lytic phages may be suitable alternatives for treating these multidrug-resistant Salmonella infections. Most Salmonella phages to date were collected from human-impacted environments. To further explore the Salmonella phage space, and to potentially identify phages with novel characteristics, we characterized Salmonella-specific phages isolated from the Penang National Park, a conserved rainforest. Four phages with a broad lytic spectrum (kills >5 Salmonella serovars) were further characterized; they have isometric heads and cone-shaped tails, and genomes of ~39,900 bp, encoding 49 CDSs. As the genomes share a <95% sequence similarity to known genomes, the phages were classified as a new species within the genus Kayfunavirus. Interestingly, the phages displayed obvious differences in their lytic spectrum and pH stability, despite having a high sequence similarity (~99% ANI). Subsequent analysis revealed that the phages differed in the nucleotide sequence in the tail spike proteins, tail tubular proteins, and portal proteins, suggesting that the SNPs were responsible for their differing phenotypes. Our findings highlight the diversity of novel Salmonella bacteriophages from rainforest regions, which can be explored as an antimicrobial agent against MDR-Salmonella strains.

INDEL detection, the 'Achilles heel' of precise genome editing: a survey of methods for accurate profiling of gene editing induced indels.

Advances in genome editing technologies have enabled manipulation of genomes at the single base level. These technologies are based on programmable nucleases (PNs) that include meganucleases, zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs) and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated 9 (Cas9) nucleases and have given researchers the ability to delete, insert or replace genomic DNA in cells, tissues and whole organisms. The great flexibility in re-designing the genomic target specificity of PNs has vastly expanded the scope of gene editing applications in life science, and shows great promise for development of the next generation gene therapies. PN technologies share the principle of inducing a DNA double-strand break (DSB) at a user-specified site in the genome, followed by cellular repair of the induced DSB. PN-elicited DSBs are mainly repaired by the non-homologous end joining (NHEJ) and the microhomology-mediated end joining (MMEJ) pathways, which can elicit a variety of small insertion or deletion (indel) mutations. If indels are elicited in a protein coding sequence and shift the reading frame, targeted gene knock out (KO) can readily be achieved using either of the available PNs. Despite the ease by which gene inactivation in principle can be achieved, in practice, successful KO is not only determined by the efficiency of NHEJ and MMEJ repair; it also depends on the design and properties of the PN utilized, delivery format chosen, the preferred indel repair outcomes at the targeted site, the chromatin state of the target site and the relative activities of the repair pathways in the edited cells. These variables preclude accurate prediction of the nature and frequency of PN induced indels. A key step of any gene KO experiment therefore becomes the detection, characterization and quantification of the indel(s) induced at the targeted genomic site in cells, tissues or whole organisms. In this survey, we briefly review naturally occurring indels and their detection. Next, we review the methods that have been developed for detection of PN-induced indels. We briefly outline the experimental steps and describe the pros and cons of the various methods to help users decide a suitable method for their editing application. We highlight recent advances that enable accurate and sensitive quantification of indel events in cells regardless of their genome complexity, turning a complex pool of different indel events into informative indel profiles. Finally, we review what has been learned about PN-elicited indel formation through the use of the new methods and how this insight is helping to further advance the genome editing field.

Normal Neurodevelopment and Fertility in Juvenile Male Rats Exposed to Polyethylene Glycol Following Dosing With PEGylated rFIX (Nonacog Beta Pegol, N9-GP): Evidence from a 10-Week Repeat-Dose Toxicity Study.

N9-GP/Rebinyn®/Refixia® is an approved PEGylated (polyethylene glycol-conjugated) recombinant human factor IX intended for prophylactic and/or on-demand treatment in adults and children with haemophilia B. A juvenile neurotoxicity study was conducted in male rats to evaluate effects on neurodevelopment, sexual maturation, and fertility following repeat-dosing of N9-GP. Male rats were dosed twice weekly from Day 21 of age with N9-GP or vehicle for 10 weeks, followed by a dosing-free recovery period for 13 weeks and terminated throughout the dosing and recovery periods. Overall, dosing N9-GP to juvenile rats did not result in any functional or pathological effects, as measured by neurobehavioural/neurocognitive tests, including motor activity, sensory function, learning and memory as well as growth, sexual maturation, and fertility. This was further supported by the extensive histopathologic evaluation of brain tissue. Exposure and distribution of polyethylene glycol was investigated in plasma, choroid plexus, cerebrospinal fluid, and brain sections. PEG did not cross the blood brain barrier and PEG exposure did not result in any effects on neurodevelopment. In conclusion, dosing of N9-GP to juvenile rats did not identify any effects on growth, sexual maturation and fertility, clinical and histological pathology, or neurodevelopment related to PEG exposure and supports the prophylactic use of N9-GP in children.

Golgi-localized exo-ß1,3-galactosidases involved in cell expansion and root growth in Arabidopsis.

Plant arabinogalactan proteins (AGPs) are a diverse group of cell surface- and wall-associated glycoproteins. Functionally important AGP glycans are synthesized in the Golgi apparatus, but the relationships among their glycosylation levels, processing, and functionalities are poorly understood. Here, we report the identification and functional characterization of two Golgi-localized exo-ß-1,3-galactosidases from the glycosyl hydrolase 43 (GH43) family in Arabidopsis thaliana GH43 loss-of-function mutants exhibited root cell expansion defects in sugar-containing growth media. This root phenotype was associated with an increase in the extent of AGP cell wall association, as demonstrated by Yariv phenylglycoside dye quantification and comprehensive microarray polymer profiling of sequentially extracted cell walls. Characterization of recombinant GH43 variants revealed that the exo-ß-1,3-galactosidase activity of GH43 enzymes is hindered by ß-1,6 branches on ß-1,3-galactans. In line with this steric hindrance, the recombinant GH43 variants did not release galactose from cell wall-extracted glycoproteins or AGP-rich gum arabic. These results indicate that the lack of exo-ß-1,3-galactosidase activity alters cell wall extensibility in roots, a phenotype that could be explained by the involvement of galactosidases in AGP glycan biosynthesis.

Exocyst mutants suppress pollen tube growth and cell wall structural defects of hydroxyproline O-arabinosyltransferase mutants.

HYDROXYPROLINE O-ARABINOSYLTRANSFERASEs (HPATs) initiate a post-translational protein modification (Hyp-Ara) found abundantly on cell wall structural proteins. In Arabidopsis thaliana, HPAT1 and HPAT3 are redundantly required for full pollen fertility. In addition to the lack of Hyp-Ara in hpat1/3 pollen tubes (PTs), we also found broadly disrupted cell wall polymer distributions, particularly the conversion of the tip cell wall to a more shaft-like state. Mutant PTs were slow growing and prone to rupture and morphological irregularities. In a forward mutagenesis screen for suppressors of the hpat1/3 low seed-set phenotype, we identified a missense mutation in exo70a2, a predicted member of the vesicle-tethering exocyst complex. The suppressed pollen had increased fertility, fewer morphological defects and partially rescued cell wall organization. A transcriptional null allele of exo70a2 also suppressed the hpat1/3 fertility phenotype, as did mutants of core exocyst complex member sec15a, indicating that reduced exocyst function bypassed the PT requirement for Hyp-Ara. In a wild-type background, exo70a2 reduced male transmission efficiency, lowered pollen germination frequency and slowed PT elongation. EXO70A2 also localized to the PT tip plasma membrane, consistent with a role in exocyst-mediated secretion. To monitor the trafficking of Hyp-Ara modified proteins, we generated an HPAT-targeted fluorescent secretion reporter. Reporter secretion was partially dependent on EXO70A2 and was significantly increased in hpat1/3 PTs compared with the wild type, but was reduced in the suppressed exo70a2 hpat1/3 tubes.

An improved direct metamobilome approach increases the detection of larger-sized circular elements across kingdoms.

Mobile genetic elements (MGEs) are instrumental in natural prokaryotic genome editing, permitting genome plasticity and allowing microbes to accumulate genetic diversity. MGEs serve as a vast communal gene pool and include DNA elements such as plasmids and bacteriophages (phages) among others. These mobile DNA elements represent a human health risk as they can introduce new traits, such as antibiotic resistance or virulence, to a bacterial strain. Sequencing libraries targeting environmental circular MGEs, referred to as metamobilomes, may broaden our current understanding of the mechanisms behind the mobility, prevalence and content of these elements. However, metamobilomics is affected by a severe bias towards small circular elements, introduced by multiple displacement amplification (MDA). MDA is typically used to overcome limiting DNA quantities after the removal of non-circular DNA during library preparations. By examining the relationship between sequencing coverage and the size of circular MGEs in paired metamobilome datasets with and without MDA, we show that larger circular elements are lost when using MDA. This study is the first to systematically demonstrate that MDA is detrimental to detecting larger-sized plasmids if small plasmids are present. It is also the first to show that MDA can be omitted when using enzyme-based DNA fragmentation and PCR in library preparation kits such as Nextera XT® from Illumina.

Population genomics of grey wolves and wolf-like canids in North America.

North America is currently home to a number of grey wolf (Canis lupus) and wolf-like canid populations, including the coyote (Canis latrans) and the taxonomically controversial red, Eastern timber and Great Lakes wolves. We explored their population structure and regional gene flow using a dataset of 40 full genome sequences that represent the extant diversity of North American wolves and wolf-like canid populations. This included 15 new genomes (13 North American grey wolves, 1 red wolf and 1 Eastern timber/Great Lakes wolf), ranging from 0.4 to 15x coverage. In addition to providing full genome support for the previously proposed coyote-wolf admixture origin for the taxonomically controversial red, Eastern timber and Great Lakes wolves, the discriminatory power offered by our dataset suggests all North American grey wolves, including the Mexican form, are monophyletic, and thus share a common ancestor to the exclusion of all other wolves. Furthermore, we identify three distinct populations in the high arctic, one being a previously unidentified "Polar wolf" population endemic to Ellesmere Island and Greenland. Genetic diversity analyses reveal particularly high inbreeding and low heterozygosity in these Polar wolves, consistent with long-term isolation from the other North American wolves.

Engineering mammalian cells to produce plant-specific N-glycosylation on proteins.

Protein N-glycosylation is an essential and highly conserved posttranslational modification found in all eukaryotic cells. Yeast, plants and mammalian cells, however, produce N-glycans with distinct structural features. These species-specific features not only pose challenges in selecting host cells for production of recombinant therapeutics for human medical use but also provide opportunities to explore and utilize species-specific glycosylation in design of vaccines. Here, we used reverse cross-species engineering to stably introduce plant core α3fucose (α3Fuc) and ß2xylose (ß2Xyl) N-glycosylation epitopes in the mammalian Chinese hamster ovary (CHO) cell line. We used directed knockin of plant core fucosylation and xylosylation genes (AtFucTA/AtFucTB and AtXylT) and targeted knockout of endogenous genes for core fucosylation (fut8) and elongation (B4galt1), for establishing CHO cells with plant N-glycosylation capacities. The engineering was evaluated through coexpression of two human therapeutic N-glycoproteins, erythropoietin (EPO) and an immunoglobulin G (IgG) antibody. Full conversion to the plant-type α3Fuc/ß2Xyl bi-antennary agalactosylated N-glycosylation (G0FX) was demonstrated for the IgG1 produced in CHO cells. These results demonstrate that N-glycosylation in mammalian cells is amenable for extensive cross-kingdom engineering and that engineered CHO cells may be used to produce glycoproteins with plant glycosylation.

Phenolic cross-links: building and de-constructing the plant cell wall.

Covering: Up to 2019Phenolic cross-links and phenolic inter-unit linkages result from the oxidative coupling of two hydroxycinnamates or two molecules of tyrosine. Free dimers of hydroxycinnamates, lignans, play important roles in plant defence. Cross-linking of bound phenolics in the plant cell wall affects cell expansion, wall strength, digestibility, degradability, and pathogen resistance. Cross-links mediated by phenolic substituents are particularly important as they confer strength to the wall via the formation of new covalent bonds, and by excluding water from it. Four biopolymer classes are known to be involved in the formation of phenolic cross-links: lignins, extensins, glucuronoarabinoxylans, and side-chains of rhamnogalacturonan-I. Lignins and extensins are ubiquitous in streptophytes whereas aromatic substituents on xylan and pectic side-chains are commonly assumed to be particular features of Poales sensu lato and core Caryophyllales, respectively. Cross-linking of phenolic moieties proceeds via radical formation, is catalyzed by peroxidases and laccases, and involves monolignols, tyrosine in extensins, and ferulate esters on xylan and pectin. Ferulate substituents, on xylan in particular, are thought to be nucleation points for lignin polymerization and are, therefore, of paramount importance to wall architecture in grasses and for the development of technology for wall disassembly, e.g. for the use of grass biomass for production of 2nd generation biofuels. This review summarizes current knowledge on the intra- and extracellular acylation of polysaccharides, and inter- and intra-molecular cross-linking of different constituents. Enzyme mediated lignan in vitro synthesis for pharmaceutical uses are covered as are industrial exploitation of mutant and transgenic approaches to control cell wall cross-linking.

NetSurfP-2.0: Improved prediction of protein structural features by integrated deep learning.

The ability to predict local structural features of a protein from the primary sequence is of paramount importance for unraveling its function in absence of experimental structural information. Two main factors affect the utility of potential prediction tools: their accuracy must enable extraction of reliable structural information on the proteins of interest, and their runtime must be low to keep pace with sequencing data being generated at a constantly increasing speed. Here, we present NetSurfP-2.0, a novel tool that can predict the most important local structural features with unprecedented accuracy and runtime. NetSurfP-2.0 is sequence-based and uses an architecture composed of convolutional and long short-term memory neural networks trained on solved protein structures. Using a single integrated model, NetSurfP-2.0 predicts solvent accessibility, secondary structure, structural disorder, and backbone dihedral angles for each residue of the input sequences. We assessed the accuracy of NetSurfP-2.0 on several independent test datasets and found it to consistently produce state-of-the-art predictions for each of its output features. We observe a correlation of 80% between predictions and experimental data for solvent accessibility, and a precision of 85% on secondary structure 3-class predictions. In addition to improved accuracy, the processing time has been optimized to allow predicting more than 1000 proteins in less than 2 hours, and complete proteomes in less than 1 day.