Genome-wide DNA methylation landscape of four Chinese populations and epigenetic variation linked to Tibetan high-altitude adaptation.

DNA methylation (DNAm) is one of the major epigenetic mechanisms in humans and is important in diverse cellular processes. The variation of DNAm in the human population is related to both genetic and environmental factors. However, the DNAm profiles have not been investigated in the Chinese population of diverse ethnicities. Here, we performed double-strand bisulfite sequencing (DSBS) for 32 Chinese individuals representing four major ethnic groups including Han Chinese, Tibetan, Zhuang, and Mongolian. We identified a total of 604,649 SNPs and quantified DNAm at more than 14 million CpGs in the population. We found global DNAm-based epigenetic structure is different from the genetic structure of the population, and ethnic difference only partially explains the variation of DNAm. Surprisingly, non-ethnic-specific DNAm variations showed stronger correlation with the global genetic divergence than these ethnic-specific DNAm. Differentially methylated regions (DMRs) among these ethnic groups were found around genes in diverse biological processes. Especially, these DMR-genes between Tibetan and non-Tibetans were enriched around high-altitude genes including EPAS1 and EGLN1, suggesting DNAm alteration plays an important role in high-altitude adaptation. Our results provide the first batch of epigenetic maps for Chinese populations and the first evidence of the association of epigenetic changes with Tibetans' high-altitude adaptation.

Ruminant-specific genes identified using high-quality genome data and their roles in rumen evolution.

Ruminants comprise a highly successful group of mammals with striking morphological innovations, including the presence of a rumen. Many studies have shown that species-specific or lineage-specific genes (referred to as new genes) play important roles in phenotypic evolution. In this study, we identified 1064 ruminant-specific genes based on the newly assembled high-quality genomes of representative members of two ruminant families and other publically available high-quality genomes. Ruminant-specific genes shared similar evolutionary and expression patterns with new genes found in other mammals, such as primates and rodents. Most new genes were derived from gene duplication and tended to be expressed in the testes or immune-related tissues, but were depleted in the adult brain. We also found that most genes expressed in the rumen were genes predating sheep-sperm whale split (referred to as old genes), but some new genes were also involved in the evolution of the rumen, and contributed more during rumen development than in the adult rumen. Notably, expression levels of members of the ruminant-specific PRD-SPRRII gene family, which are subject to positive selection, varied throughout rumen development and may thus play important roles in the development of the keratin-rich surface of the rumen. Overall, this study generated two novel ruminant genomes and also provided novel insights into the evolution of new mammalian organs.

Production and Characterization of Biotinylated Anti-fenitrothion Nanobodies and Development of Sensitive Fluoroimmunoassay.

A simple and sensitive fluoroimmunoassay (FIA) based on a heavy-chain antibody (VHH) for rapid detection of fenitrothion was developed. A VHH library was constructed from an immunized alpaca, and one clone recognizing fenitrothion (namely, VHHjd8) was achieved after careful biopanning. It was biotinylated by fusing with the Avi tag and biotin ligase to obtain a fusion protein (VHHjd8-BT), showing both binding capacity to fenitrothion and the streptavidin poly-horseradish peroxidase conjugate (SA-polyHRP). Based on a competitive assay format, the absorbance spectrum of oxidized 3,3',5,5'-tetramethylbenzidine generated by SA-polyHRP overlapped the emission spectrum of carbon dots, which resulted in quenching of signals due to the inner-filter effect. The developed FIA showed an IC50 value of 1.4 ng/mL and a limit of detection of 0.03 ng/mL, which exhibited 15-fold improvement compared with conventional enzyme-linked immunosorbent assay. The recovery test of FIA was validated by standard GC-MS/MS, and the results showed good consistency, indicating that the assay is an ideal tool for rapid screening of fenitrothion in bulk food samples.

Molecular mechanisms and topological consequences of drastic chromosomal rearrangements of muntjac deer.

Muntjac deer have experienced drastic karyotype changes during their speciation, making it an ideal model for studying mechanisms and functional consequences of mammalian chromosome evolution. Here we generated chromosome-level genomes for Hydropotes inermis (2n = 70), Muntiacus reevesi (2n = 46), female and male M. crinifrons (2n = 8/9) and a contig-level genome for M. gongshanensis (2n = 8/9). These high-quality genomes combined with Hi-C data allowed us to reveal the evolution of 3D chromatin architectures during mammalian chromosome evolution. We find that the chromosome fusion events of muntjac species did not alter the A/B compartment structure and topologically associated domains near the fusion sites, but new chromatin interactions were gradually established across the fusion sites. The recently borne neo-Y chromosome of M. crinifrons, which underwent male-specific inversions, has dramatically restructured chromatin compartments, recapitulating the early evolution of canonical mammalian Y chromosomes. We also reveal that a complex structure containing unique centromeric satellite, truncated telomeric and palindrome repeats might have mediated muntjacs' recurrent chromosome fusions. These results provide insights into the recurrent chromosome tandem fusion in muntjacs, early evolution of mammalian sex chromosomes, and reveal how chromosome rearrangements can reshape the 3D chromatin regulatory conformations during species evolution.

A towering genome: Experimentally validated adaptations to high blood pressure and extreme stature in the giraffe.

The suite of adaptations associated with the extreme stature of the giraffe has long interested biologists and physiologists. By generating a high-quality chromosome-level giraffe genome and a comprehensive comparison with other ruminant genomes, we identified a robust catalog of giraffe-specific mutations. These are primarily related to cardiovascular, bone growth, vision, hearing, and circadian functions. Among them, the giraffe FGFRL1 gene is an outlier with seven unique amino acid substitutions not found in any other ruminant. Gene-edited mice with the giraffe-type FGFRL1 show exceptional hypertension resistance and higher bone mineral density, both of which are tightly connected with giraffe adaptations to high stature. Our results facilitate a deeper understanding of the molecular mechanism underpinning distinct giraffe traits, and may provide insights into the study of hypertension in humans.

Genomic and experimental data provide new insights into luciferin biosynthesis and bioluminescence evolution in fireflies.

Fireflies are among the most charismatic insects for their spectacular bioluminescence, but the origin and evolution of bioluminescence remain elusive. Especially, the genic basis of luciferin (D-luciferin) biosynthesis and light patterns is largely unknown. Here, we present the high-quality reference genomes of two fireflies Lamprigera yunnana (1053 Mb) and Abscondita terminalis (501 Mb) with great differences in both morphology and luminous behavior. We sequenced the transcriptomes and proteomes of luminous organs of two species. We created the CRISPR/Cas9-induced mutants of Abdominal B gene without luminous organs in the larvae of A. terminalis and sequenced the transcriptomes of mutants and wild-types. Combining gene expression analyses with comparative genomics, we propose a more complete luciferin synthesis pathway, and confirm the convergent evolution of bioluminescence in insects. Using experiments, the function of the firefly acyl-CoA thioesterase (ACOT1) to convert L-luciferin to D-luciferin was validated for the first time. Comparisons of three-dimension reconstruction of luminous organs and their differentially expressed genes among two species suggest that two positive genes in the calcium signaling pathway and structural difference of luminous organs may play an important role in the evolution of flash pattern. Altogether, our results provide important resources for further exploring bioluminescence in insects.

Large-scale ruminant genome sequencing provides insights into their evolution and distinct traits.

The ruminants are one of the most successful mammalian lineages, exhibiting morphological and habitat diversity and containing several key livestock species. To better understand their evolution, we generated and analyzed de novo assembled genomes of 44 ruminant species, representing all six Ruminantia families. We used these genomes to create a time-calibrated phylogeny to resolve topological controversies, overcoming the challenges of incomplete lineage sorting. Population dynamic analyses show that population declines commenced between 100,000 and 50,000 years ago, which is concomitant with expansion in human populations. We also reveal genes and regulatory elements that possibly contribute to the evolution of the digestive system, cranial appendages, immune system, metabolism, body size, cursorial locomotion, and dentition of the ruminants.

Biological adaptations in the Arctic cervid, the reindeer (Rangifer tarandus).

The reindeer is an Arctic species that exhibits distinctive biological characteristics, for which the underlying genetic basis remains largely unknown. We compared the genomes of reindeer against those of other ruminants and nonruminant mammals to reveal the genetic basis of light arrhythmicity, high vitamin D metabolic efficiency, the antler growth trait of females, and docility. We validate that two reindeer vitamin D metabolic genes (CYP27B1 and POR) show signs of positive selection and exhibit higher catalytic activity than those of other ruminants. A mutation upstream of the reindeer CCND1 gene endows an extra functional binding motif of the androgen receptor and thereby may result in female antlers. Furthermore, a mutation (proline-1172→threonine) in reindeer PER2 results in loss of binding ability with CRY1, which may explain circadian arrhythmicity in reindeer.

Genetic basis of ruminant headgear and rapid antler regeneration.

Ruminants are the only extant mammalian group possessing bony (osseous) headgear. We obtained 221 transcriptomes from bovids and cervids and sequenced three genomes representing the only two pecoran lineages that convergently lack headgear. Comparative analyses reveal that bovid horns and cervid antlers share similar gene expression profiles and a common cellular basis developed from neural crest stem cells. The rapid regenerative properties of antler tissue involve exploitation of oncogenetic pathways, and at the same time some tumor suppressor genes are under strong selection in deer. These results provide insights into the evolutionary origin of ruminant headgear as well as mammalian organ regeneration and oncogenesis.

Draft genome of the reindeer (Rangifer tarandus).

Background: The reindeer (Rangifer tarandus) is the only fully domesticated species in the Cervidae family, and it is the only cervid with a circumpolar distribution. Unlike all other cervids, female reindeer, as well as males, regularly grow cranial appendages (antlers, the defining characteristics of cervids). Moreover, reindeer milk contains more protein and less lactose than bovids' milk. A high-quality reference genome of this species will assist efforts to elucidate these and other important features in the reindeer. Findings: We obtained 615 Gb (Gigabase) of usable sequences by filtering the low-quality reads of the raw data generated from the Illumina Hiseq 4000 platform, and a 2.64-Gb final assembly, representing 95.7% of the estimated genome (2.76 Gb according to k-mer analysis), including 92.6% of expected genes according to BUSCO analysis. The contig N50 and scaffold N50 sizes were 89.7 kilo base (kb) and 0.94 mega base (Mb), respectively. We annotated 21 555 protein-coding genes and 1.07 Gb of repetitive sequences by de novo and homology-based prediction. Homology-based searches detected 159 rRNA, 547 miRNA, 1339 snRNA, and 863 tRNA sequences in the genome of R. tarandus. The divergence time between R. tarandus and ancestors of Bos taurus and Capra hircus is estimated to be about 29.5 million years ago. Conclusions: Our results provide the first high-quality reference genome for the reindeer and a valuable resource for studying the evolution, domestication, and other unusual characteristics of the reindeer.