A Chinese indicine pangenome reveals a wealth of novel structural variants introgressed from other Bos species.

Chinese indicine cattle harbor a much higher genetic diversity compared with other domestic cattle, but their genome architecture remains uninvestigated. Using PacBio HiFi sequencing data from 10 Chinese indicine cattle across southern China, we assembled 20 high-quality partially phased genomes and integrated them into a multiassembly graph containing 148.5 Mb (5.6%) of novel sequence. We identified 156,009 high-confidence nonredundant structural variants (SVs) and 206 SV hotspots spanning ∼195 Mb of gene-rich sequence. We detected 34,249 archaic introgressed fragments in Chinese indicine cattle covering 1.93 Gb (73.3%) of the genome. We inferred an average of 3.8%, 3.2%, 1.4%, and 0.5% of introgressed sequence originating, respectively, from banteng-like, kouprey-like, gayal-like, and gaur-like Bos species, as well as 0.6% of unknown origin. Introgression from multiple donors might have contributed to the genetic diversity of Chinese indicine cattle. Altogether, this study highlights the contribution of interspecies introgression to the genomic architecture of an important livestock population and shows how exotic genomic elements can contribute to the genetic variation available for selection.

A sheep pangenome reveals the spectrum of structural variations and their effects on tail phenotypes.

Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here we generated high-quality genome assemblies for 15 individuals from genetically diverse sheep breeds using Pacific Biosciences (PacBio) high-fidelity sequencing, discovering 130.3 Mb nonreference sequences, from which 588 genes were annotated. A total of 149,158 biallelic insertions/deletions, 6531 divergent alleles, and 14,707 multiallelic variations with precise breakpoints were discovered. The SV spectrum is characterized by an excess of derived insertions compared to deletions (94,422 vs. 33,571), suggesting recent active LINE expansions in sheep. Nearly half of the SVs display low to moderate linkage disequilibrium with surrounding single-nucleotide polymorphisms (SNPs) and most SVs cannot be tagged by SNP probes from the widely used ovine 50K SNP chip. We identified 865 population-stratified SVs including 122 SVs possibly derived in the domestication process among 690 individuals from sheep breeds worldwide. A novel 168-bp insertion in the 5' untranslated region (5' UTR) of HOXB13 is found at high frequency in long-tailed sheep. Further genome-wide association study and gene expression analyses suggest that this mutation is causative for the long-tail trait. In summary, we have developed a panel of high-quality de novo assemblies and present a catalog of structural variations in sheep. Our data capture abundant candidate functional variations that were previously unexplored and provide a fundamental resource for understanding trait biology in sheep.

A missense mutation in RRM1 contributes to animal tameness.

The increased tameness to reduce avoidance of human in wild animals has been long proposed as the key step of animal domestication. The tameness is a complex behavior trait and largely determined by genetic factors. However, the underlying genetic mutations remain vague and how they influence the animal behaviors is yet to be explored. Behavior tests of a wild-domestic hybrid goat population indicate the locus under strongest artificial selection during domestication may exert a huge effect on the flight distance. Within this locus, only one missense mutation RRM1I241V which was present in the early domestic goat ~6500 years ago. Genome editing of RRM1I241V in mice showed increased tameness and sociability and reduced anxiety. These behavioral changes induced by RRM1I241V were modulated by the alternation of activity of glutamatergic synapse and some other synapse-related pathways. This study established a link between RRM1I241V and tameness, demonstrating that the complex behavioral change can be achieved by mutations under strong selection during animal domestication.

Chromosome-level genome assembly for takin (Budorcas taxicolor) provides insights into its taxonomic status and genetic diversity.

The takin (Budorcas taxicolor) is one of the largest bovid herbivores in the subfamily Caprinae. The takin is at high risk of extinction, but its taxonomic status and genetic diversity remain unclear. In this study, we constructed the first reference genome of Bu. taxicolor using PacBio long High-Fidelity reads and Hi-C technology. The assembled genome is ~2.95 Gb with a contig N50 of 68.05 Mb, which were anchored onto 25+XY chromosomes. We found that the takin was more closely related to muskox than to other Caprinae species. Compared to the common ancestral karyotype of bovidae (2n = 60), we found the takin (2n = 52) experienced four chromosome fusions and one large translocation. Furthermore, we resequenced nine golden takins from the main distribution area, the Qinling Mountains, and identified 3.3 million single nucleotide polymorphisms. The genetic diversity of takin was very low (θπ = 0.00028 and heterozygosity =0.00038), among the lowest detected in domestic and wild mammals. Takin genomes showed a high inbreeding coefficient (FROH =0.217), suggesting severe inbreeding depression. The demographic history showed that the effective population size of takins declined significantly from ~100,000 years ago. Our results provide valuable information for protection of takins and insights into their evolution.

Chromosome-level Genomes Reveal the Genetic Basis of Descending Dysploidy and Sex Determination in Morus Plants.

Multiple plant lineages have independently evolved sex chromosomes and variable karyotypes to maintain their sessile lifestyles through constant biological innovation. Morus notabilis, a dioecious mulberry species, has the fewest chromosomes among Morus spp., but the genetic basis of sex determination and karyotype evolution in this species has not been identified. In this study, three high-quality genome assemblies were generated for Morus spp. [including dioecious M. notabilis (male and female) and Morus yunnanensis (female)] with genome sizes of 301-329 Mb and were grouped into six pseudochromosomes. Using a combination of genomic approaches, we found that the putative ancestral karyotype of Morus species was close to 14 protochromosomes, and that several chromosome fusion events resulted in descending dysploidy (2n = 2x = 12). We also characterized a ∼ 6.2-Mb sex-determining region on chromosome 3. Four potential male-specific genes, a partially duplicatedDNA helicase gene (named MSDH) and three Ty3_Gypsy long terminal repeat retrotransposons (named MSTG1/2/3), were identified in the Y-linked area and considered to be strong candidate genes for sex determination or differentiation. Population genomic analysis showed that Guangdong accessions in China were genetically similar to Japanese accessions of mulberry. In addition, genomic areas containing selective sweeps that distinguish domesticated mulberry from wild populations in terms of flowering and disease resistance were identified. Our study provides an important genetic resource for sex identification research and molecular breeding in mulberry.

Genomic insights into the phylogeny and biomass-degrading enzymes of rumen ciliates.

Understanding the biodiversity and genetics of gut microbiomes has important implications for host physiology and industrial enzymes, whereas most studies have been focused on bacteria and archaea, and to a lesser extent on fungi and viruses. One group, still underexplored and elusive, is ciliated protozoa, despite its importance in shaping microbiota populations. Integrating single-cell sequencing and an assembly-and-identification pipeline, we acquired 52 high-quality ciliate genomes of 22 rumen morphospecies from 11 abundant morphogenera. With these genomes, we resolved the taxonomic and phylogenetic framework that revised the 22 morphospecies into 19 species spanning 13 genera and reassigned the genus Dasytricha from Isotrichidae to a new family Dasytrichidae. Comparative genomic analyses revealed that extensive horizontal gene transfers and gene family expansion provided rumen ciliate species with a broad array of carbohydrate-active enzymes (CAZymes) to degrade all major kinds of plant and microbial carbohydrates. In particular, the genomes of Diplodiniinae and Ophryoscolecinae species encode as many CAZymes as gut fungi, and ~80% of their degradative CAZymes act on plant cell-wall. The activities of horizontally transferred cellulase and xylanase of ciliates were experimentally verified and were 2-9 folds higher than those of the inferred corresponding bacterial donors. Additionally, the new ciliate dataset greatly facilitated rumen metagenomic analyses by allowing ~12% of the metagenomic sequencing reads to be classified as ciliate sequences.

Genome-Wide DNA Methylation Patterns of Muscle and Tail-Fat in DairyMeade Sheep and Mongolian Sheep.

This study aimed to explore the genome-wide DNA methylation differences between muscle and tail-fat tissues of DairyMeade sheep (thin-tailed, lean carcass) and Mongolian sheep (fat-tailed, fat-deposited carcass). Whole-genome bisulfite sequencing (WGBS) was conducted and the global DNA methylation dynamics were mapped. Generally, CGs had a higher DNA methylation level than CHHs and CHGs, and tail-fat tissues had higher CG methylation levels than muscle tissues. For DNA repeat elements, SINE had the highest methylation level, while Simple had the lowest. When dividing the gene promoter region into small bins (200 bp per bin), the bins near the transcription start site (±200 bp) had the highest CG count per bin but the lowest DNA methylation levels. A series of DMRs were identified in muscle and tail-fat tissues between the two breeds. Among them, the introns of gene CAMK2D (calcium/calmodulin-dependent protein kinase II δ) demonstrated significant DNA methylation level differences between the two breeds in both muscle and tail-fat tissues, and it may play a crucial role in fat metabolism and meat quality traits. This study may provide basic datasets and references for further epigenetic modification studies during sheep genetic improvement.

De Novo Assembly of 20 Chicken Genomes Reveals the Undetectable Phenomenon for Thousands of Core Genes on Microchromosomes and Subtelomeric Regions.

The gene numbers and evolutionary rates of birds were assumed to be much lower than those of mammals, which is in sharp contrast to the huge species number and morphological diversity of birds. It is, therefore, necessary to construct a complete avian genome and analyze its evolution. We constructed a chicken pan-genome from 20 de novo assembled genomes with high sequencing depth, and identified 1,335 protein-coding genes and 3,011 long noncoding RNAs not found in GRCg6a. The majority of these novel genes were detected across most individuals of the examined transcriptomes but were seldomly measured in each of the DNA sequencing data regardless of Illumina or PacBio technology. Furthermore, different from previous pan-genome models, most of these novel genes were overrepresented on chromosomal subtelomeric regions and microchromosomes, surrounded by extremely high proportions of tandem repeats, which strongly blocks DNA sequencing. These hidden genes were proved to be shared by all chicken genomes, included many housekeeping genes, and enriched in immune pathways. Comparative genomics revealed the novel genes had 3-fold elevated substitution rates than known ones, updating the knowledge about evolutionary rates in birds. Our study provides a framework for constructing a better chicken genome, which will contribute toward the understanding of avian evolution and the improvement of poultry breeding.

A near complete genome for goat genetic and genomic research.

BACKGROUND: Goat, one of the first domesticated livestock, is a worldwide important species both culturally and economically. The current goat reference genome, known as ARS1, is reported as the first nonhuman genome assembly using 69× PacBio sequencing. However, ARS1 suffers from incomplete X chromosome and highly fragmented Y chromosome scaffolds. RESULTS: Here, we present a very high-quality de novo genome assembly, Saanen_v1, from a male Saanen dairy goat, with the first goat Y chromosome scaffold based on 117× PacBio long-read sequencing and 118× Hi-C data. Saanen_v1 displays a high level of completeness thanks to the presence of centromeric and telomeric repeats at the proximal and distal ends of two-thirds of the autosomes, and a much reduced number of gaps (169 vs. 773). The completeness and accuracy of the Saanen_v1 genome assembly are also evidenced by more assembled sequences on the chromosomes (2.63 Gb for Saanen_v1 vs. 2.58 Gb for ARS1), a slightly increased mapping ratio for transcriptomic data, and more genes anchored to chromosomes. The eight putative large assembly errors (1 to ~ 7 Mb each) found in ARS1 were amended, and for the first time, the substitution rate of this ruminant Y chromosome was estimated. Furthermore, sequence improvement in Saanen_v1, compared with ARS1, enables us to assign the likely correct positions for 4.4% of the single nucleotide polymorphism (SNP) probes in the widely used GoatSNP50 chip. CONCLUSIONS: The updated goat genome assembly including both sex chromosomes (X and Y) and the autosomes with high-resolution quality will serve as a valuable resource for goat genetic research and applications.

Genomic analyses reveal distinct genetic architectures and selective pressures in Chinese donkeys.

Donkey (Equus asinus) is an important livestock animal in China because of its draft and medicinal value. After a long period of natural and artificial selection, the variety and phenotype of donkeys have become abundant. We clarified the genetic and demographic characteristics of Chinese domestic donkeys and the selection pressures by analyzing 78 whole genomes from 12 breeds. According to population structure, most Chinese domestic donkeys showed a dominant ancestral type. However, the Chinese donkeys still represented a significant geographical distribution trend. In the selective sweep, gene annotation, functional enrichment, and differential expression analyses between large and small donkey groups, we identified selective signals, including NCAPG and LCORL, which are related to rapid growth and large body size. Our findings elucidate the evolutionary history and formation of different donkey breeds and provide theoretical insights into the genetic mechanism underlying breed characteristics and molecular breeding programs of donkey clades.

A new duck genome reveals conserved and convergently evolved chromosome architectures of birds and mammals.

BACKGROUND: Ducks have a typical avian karyotype that consists of macro- and microchromosomes, but a pair of much less differentiated ZW sex chromosomes compared to chickens. To elucidate the evolution of chromosome architectures between ducks and chickens, and between birds and mammals, we produced a nearly complete chromosomal assembly of a female Pekin duck by combining long-read sequencing and multiplatform scaffolding techniques. RESULTS: A major improvement of genome assembly and annotation quality resulted from the successful resolution of lineage-specific propagated repeats that fragmented the previous Illumina-based assembly. We found that the duck topologically associated domains (TAD) are demarcated by putative binding sites of the insulator protein CTCF, housekeeping genes, or transitions of active/inactive chromatin compartments, indicating conserved mechanisms of spatial chromosome folding with mammals. There are extensive overlaps of TAD boundaries between duck and chicken, and also between the TAD boundaries and chromosome inversion breakpoints. This suggests strong natural selection pressure on maintaining regulatory domain integrity, or vulnerability of TAD boundaries to DNA double-strand breaks. The duck W chromosome retains 2.5-fold more genes relative to chicken. Similar to the independently evolved human Y chromosome, the duck W evolved massive dispersed palindromic structures, and a pattern of sequence divergence with the Z chromosome that reflects stepwise suppression of homologous recombination. CONCLUSIONS: Our results provide novel insights into the conserved and convergently evolved chromosome features of birds and mammals, and also importantly add to the genomic resources for poultry studies.

A Hu sheep genome with the first ovine Y chromosome reveal introgression history after sheep domestication.

The Y chromosome plays key roles in male fertility and reflects the evolutionary history of paternal lineages. Here, we present a de novo genome assembly of the Hu sheep with the first draft assembly of ovine Y chromosome (oMSY), using nanopore sequencing and Hi-C technologies. The oMSY that we generated spans 10.6 Mb from which 775 Y-SNPs were identified by applying a large panel of whole genome sequences from worldwide sheep and wild Iranian mouflons. Three major paternal lineages (HY1a, HY1b and HY2) were defined across domestic sheep, of which HY2 was newly detected. Surprisingly, HY2 forms a monophyletic clade with the Iranian mouflons and is highly divergent from both HY1a and HY1b. Demographic analysis of Y chromosomes, mitochondrial and nuclear genomes confirmed that HY2 and the maternal counterpart of lineage C represented a distinct wild mouflon population in Iran that diverge from the direct ancestor of domestic sheep, the wild mouflons in Southeastern Anatolia. Our results suggest that wild Iranian mouflons had introgressed into domestic sheep and thereby introduced this Iranian mouflon specific lineage carrying HY2 to both East Asian and Africa sheep populations.

Comparative Transcriptome Profiling Analysis Uncovers Novel Heterosis-Related Candidate Genes Associated with Muscular Endurance in Mules.

Heterosis has been widely exploited in animal and plant breeding programs to enhance the productive traits of hybrid progeny from two breeds or species. However, its underlying genetic mechanisms remain enigmatic. Transcriptome profiling analysis can be used as a method for exploring the mechanism of heterosis. Here, we performed genome-wide gene expression and alternative splicing (AS) analyses in different tissues (muscle, brain, and skin) from crosses between donkeys and horses. Our results indicated that 86.1% of the differentially expressed genes (DEGs) and 87.2% of the differential alternative splicing (DAS) genes showed over-dominance and dominance in muscle. Further analysis showed that the "muscle contraction" pathway was significantly enriched for both the DEGs and DAS genes in mule muscle tissue. Taken together, these DEGs and DAS genes could provide an index for future studies of the genetic and molecular mechanism of heterosis in the hybrids of donkey and horse.

Chromosome-level reference genome assembly and gene editing of the dead-leaf butterfly Kallima inachus.

The leaf resemblance of Kallima (Nymphalidae) butterflies is an important ecological adaptive mechanism that increases their survival. However, the genetic mechanism underlying ecological adaptation remains unclear owing to a dearth of genomic information. Here, we determined the karyotype (n = 31) of the dead-leaf butterfly Kallima inachus, and generated a high-quality, chromosome-level assembly (568.92 Mb; contig N50: 19.20 Mb). We also identified candidate Z and W chromosomes. To our knowledge, this is the first study to report on these aspects of this species. In the assembled genome, 15,309 protein-coding genes and 49.86% repeat elements were annotated. Phylogenetic analysis showed that K. inachus diverged from Melitaea cinxia (no leaf resemblance), both of which are in Nymphalinae, around 40 million years ago. Demographic analysis indicated that the effective population size of K. inachus decreased during the last interglacial period in the Pleistocene. The wings of adults with the pigmentary gene ebony knocked out using CRISPR/Cas9 showed phenotypes in which the orange dorsal region and entire ventral surface darkened, suggesting its vital role in the ecological adaption of dead-leaf butterflies. Our results provide important genome resources for investigating the genetic mechanism underlying protective resemblance in dead-leaf butterflies and insights into the molecular basis of protective coloration.

Chromosome-Level Reference Genome and Population Genomic Analysis Provide Insights into the Evolution and Improvement of Domesticated Mulberry (Morus alba).

Mulberry (Morus spp.) is the sole plant consumed by the domesticated silkworm. However, the genome of domesticated mulberry has not yet been sequenced, and the ploidy level of this species remains unclear. Here, we report a high-quality, chromosome-level domesticated mulberry (Morus alba) genome. Analysis of genomic data and karyotype analyses confirmed that M. alba is a diploid with 28 chromosomes (2n = 2x = 28). Population genomic analysis based on resequencing of 134 mulberry accessions classified domesticated mulberry into three geographical groups, namely, Taihu Basin of southeastern China (Hu mulberry), northern and southwestern China, and Japan. Hu mulberry had the lowest nucleotide diversity among these accessions and demonstrated obvious signatures of selection associated with environmental adaptation. Further phylogenetic analysis supports a previous proposal that multiple domesticated mulberry accessions previously classified as different species actually belong to one species. This study expands our understanding of genome evolution of the genus Morus and population structure of domesticated mulberry, which would facilitate mulberry breeding and improvement.

Chromosomal-level reference genome of Chinese peacock butterfly (Papilio bianor) based on third-generation DNA sequencing and Hi-C analysis.

BACKGROUND: Papilio bianor Cramer, 1777 (commonly known as the Chinese peacock butterfly) (Insecta, Lepidoptera, Papilionidae) is a widely distributed swallowtail butterfly with a wide number of geographic populations ranging from the southeast of Russia to China, Japan, India, Vietnam, Myanmar, and Thailand. Its wing color consists of both pigmentary colored scales (black, reddish) and structural colored scales (iridescent blue or green dust). A high-quality reference genome of P. bianor is an important foundation for investigating iridescent color evolution, phylogeography, and the evolution of swallowtail butterflies. FINDINGS: We obtained a chromosome-level de novo genome assembly of the highly heterozygous P. bianor using long Pacific Biosciences sequencing reads and high-throughput chromosome conformation capture technology. The final assembly is 421.52 Mb on 30 chromosomes (29 autosomes and 1 Z sex chromosome) with 13.12 Mb scaffold N50. In total, 15,375 protein-coding genes and 233.09 Mb of repetitive sequences were identified. Phylogenetic analyses indicated that P. bianor separated from a common ancestor of swallowtails ∼23.69-36.04 million years ago. Demographic history suggested that the population expansion of this species from the last interglacial period to the last glacial maximum possibly resulted from its decreased natural enemies and its adaptation to climate change during the glacial period. CONCLUSIONS: We present a high-quality chromosome-level reference genome of P. bianor using long-read single-molecule sequencing and Hi-C-based chromatin interaction maps. Our results lay the foundation for exploring the genetic basis of special biological features of P. bianor and also provide a useful data source for comparative genomics and phylogenomics among butterflies and moths.