Inter-individual variation in mean leukocyte telomere length (LTL) is associated with cancer and several age-associated diseases. Here, in a genome-wide meta-analysis of 37,684 individuals with replication of selected variants in a further 10,739 individuals, we identified seven loci, including five novel loci, associated with mean LTL (P<5x10−8). Five of the loci contain genes (TERC, TERT, NAF1, OBFC1, RTEL1) that are known to be involved in telomere biology. Lead SNPs at two loci (TERC and TERT) associate with several cancers and other diseases, including idiopathic pulmonary fibrosis. Moreover, a genetic risk score analysis combining lead variants at all seven loci in 22,233 coronary artery disease cases and 64,762 controls showed an association of the alleles associated with shorter LTL with increased risk of CAD (21% (95% CI: 5–35%) per standard deviation in LTL, p=0.014). Our findings support a causal role of telomere length variation in some age-related diseases.

Using next-generation sequencing technology alone, we have successfully generated and assembled a draft sequence of the giant panda genome. The assembled contigs (2.25 gigabases (Gb)) cover approximately 94% of the whole genome, and the remaining gaps (0.05 Gb) seem to contain carnivore-specific repeats and tandem repeats. Comparisons with the dog and human showed that the panda genome has a lower divergence rate. The assessment of panda genes potentially underlying some of its unique traits indicated that its bamboo diet might be more dependent on its gut microbiome than its own genetic composition. We also identified more than 2.7 million heterozygous single nucleotide polymorphisms in the diploid genome. Our data and analyses provide a foundation for promoting mammalian genetic research, and demonstrate the feasibility for using next-generation sequencing technologies for accurate, cost-effective and rapid de novo assembly of large eukaryotic genomes.

We report improved whole-genome shotgun sequences for the genomes of indica and japonica rice, both with multimegabase contiguity, or almost 1,000-fold improvement over the drafts of 2002. Tested against a nonredundant collection of 19,079 full-length cDNAs, 97.7% of the genes are aligned, without fragmentation, to the mapped super-scaffolds of one or the other genome. We introduce a gene identification procedure for plants that does not rely on similarity to known genes to remove erroneous predictions resulting from transposable elements. Using the available EST data to adjust for residual errors in the predictions, the estimated gene count is at least 38,000–40,000. Only 2%–3% of the genes are unique to any one subspecies, comparable to the amount of sequence that might still be missing. Despite this lack of variation in gene content, there is enormous variation in the intergenic regions. At least a quarter of the two sequences could not be aligned, and where they could be aligned, single nucleotide polymorphism (SNP) rates varied from as little as 3.0 SNP/kb in the coding regions to 27.6 SNP/kb in the transposable elements. A more inclusive new approach for analyzing duplication history is introduced here. It reveals an ancient whole-genome duplication, a recent segmental duplication on Chromosomes 11 and 12, and massive ongoing individual gene duplications. We find 18 distinct pairs of duplicated segments that cover 65.7% of the genome; 17 of these pairs date back to a common time before the divergence of the grasses. More important, ongoing individual gene duplications provide a never-ending source of raw material for gene genesis and are major contributors to the differences between members of the grass family.

Reduced glomerular filtration rate defines chronic kidney disease and is associated with cardiovascular and all-cause mortality. We conducted a meta-analysis of genome-wide association studies for estimated glomerular filtration rate (eGFR), combining data across 133,413 individuals with replication in up to 42,166 individuals. We identify 24 new and confirm 29 previously identified loci. Of these 53 loci, nineteen associate with eGFR among individuals with diabetes. Using bioinformatics, we show that identified genes at eGFR loci are enriched for expression in kidney tissues and in pathways relevant for kidney development and transmembrane transporter activity, kidney structure, and regulation of glucose metabolism. Chromatin state mapping and DNase I hypersensitivity analyses across adult tissues demonstrate preferential mapping of associated variants to regulatory regions in kidney but not extra-renal tissues. These findings suggest that genetic determinants of eGFR are mediated largely through direct effects within the kidney and highlight important cell types and biologic pathways.

Antidiabetic medication may modulate the gut microbiota and thereby alter plasma and faecal bile acid (BA) composition, which may improve metabolic health. Here we show that treatment with Acarbose, but not Glipizide, increases the ratio between primary BAs and secondary BAs and plasma levels of unconjugated BAs in treatment-naive type 2 diabetes (T2D) patients, which may beneficially affect metabolism. Acarbose increases the relative abundances of Lactobacillus and Bifidobacterium in the gut microbiota and depletes Bacteroides, thereby changing the relative abundance of microbial genes involved in BA metabolism. Treatment outcomes of Acarbose are dependent on gut microbiota compositions prior to treatment. Compared to patients with a gut microbiota dominated by Prevotella, those with a high abundance of Bacteroides exhibit more changes in plasma BAs and greater improvement in metabolic parameters after Acarbose treatment. Our work highlights the potential for stratification of T2D patients based on their gut microbiota prior to treatment.

Study of mutant phenotypes is a fundamental method for understanding gene function. The construction of a near-complete collection of yeast knockouts (YKO) and the unique molecular barcodes (or TAGs) that identify each strain has enabled quantitative functional profiling of Saccharomyces cerevisiae. By using these TAGs and the SGA reporter, MFA1pr-HIS3, which facilitates conversion of heterozygous diploid YKO strains into haploid mutants, we have developed a set of highly efficient microarray-based techniques, collectively referred as dSLAM (diploid-based synthetic lethality analysis on microarrays), to probe genome-wide gene-chemical and gene-gene interactions. Direct comparison revealed that these techniques are more robust than existing methods in functional profiling of the yeast genome. Widespread application of these tools will elucidate a comprehensive yeast genetic network.

The utility of CRISPR-Cas9 and TALENs for genome editing may be compromised by their off-target activity. We show that integrase-defective lentiviral vectors (IDLVs) can detect such off-target cleavage with a frequency as low as 1%. In the case of Cas9, we find frequent off-target sites with a one-base bulge or up to 13 mismatches between the single guide RNA (sgRNA) and its genomic target, which refines sgRNA design.

We report a draft sequence for the genome of the domesticated silkworm (Bombyx mori), covering 90.9% of all known silkworm genes. Our estimated gene count is 18,510, which exceeds the 13,379 genes reported for Drosophila melanogaster. Comparative analyses to fruitfly, mosquito, spider, and butterfly reveal both similarities and differences in gene content.